[Jsembla]

Hello all! I hope your're doing well.

I'm trying to add a turbulence model correction to Menter's SST that works in a very similar manner to Spalar-Allmaras' QCR but uses the specific dissipation rate as a parameter. Thus, I'm coding the model on flow_diffusion.cpp but I can't seem to find a way to bring omega into it. I see that SU2 uses pre-accumulation to store variables and bring them when in need (as far as I understand), for example, in that same script we have:

AD::StartPreacc();
AD::SetPreaccIn(V_i, nDim+9); AD::SetPreaccIn(V_j, nDim+9);
AD::SetPreaccIn(Coord_i, nDim); AD::SetPreaccIn(Coord_j, nDim);
AD::SetPreaccIn(PrimVar_Grad_i, nDim+1, nDim);
AD::SetPreaccIn(PrimVar_Grad_j, nDim+1, nDim);
AD::SetPreaccIn(turb_ke_i); AD::SetPreaccIn(turb_ke_j);
AD::SetPreaccIn(TauWall_i); AD::SetPreaccIn(TauWall_j);
AD::SetPreaccIn(Normal, nDim);

Therefore, think I should achieve something like:

AD::SetPreaccIn(omega_i); AD::SetPreaccIn(omega_j);

However, my attempts to do so have not yet been successful. Is there a way to do this? Maybe not by the means of pre-accumulation but by any other.

Thank you in advance, I hope you have a nice day.

Jsembla.

[bigfootedrockmidget]

For all SU2 related issues, please go to:
https://www.cfd-online.com/Forums/su2/


This thread is continued here:

Using pre-accumulation SU2