CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (http://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   Addition of Rotation correction to SA Turbulence Model (http://www.cfd-online.com/Forums/openfoam-programming-development/101527-addition-rotation-correction-sa-turbulence-model.html)

kalyangoparaju May 7, 2012 02:29

Addition of Rotation correction to SA Turbulence Model
 
Guys,

I am trying to simulate wing tip vortices by modified the existing SA model by adding a rotational correction to it. This was proposed by Dacles and Mariani. It is the SA-R model in the following.

http://turbmodels.larc.nasa.gov/spalart.html

According to the above link, the only major difference is the way in which S_tilda is calculated. Finally, S_tilda will have to look like

S_tilda=Omega+Crot*min(0,S-Omega)+nu_tilda*fv2/(k*d)^2

I tried to code this into openfoam the following way.

const volScalarField omega(::sqrt(2.0)*(mag(skew(fvc::grad(U_)))));
const volScalarField S(::sqrt(2.0)*(mag(symm(fvc::grad(U_)))));
const volScalarField Stilda
(
omega + scalar(2) * min(omega-omega,S-omega)
+ fv2(chi, fv1)*nuTilda_/sqr(kappa_*d_)
);


The only reason I have omega-omega is so that I don't have to define 0 as a dimentionedScalar. can someone tell me if I am doing this right.

I've tested this version and the results seem to be wrong. The results are right with the actual ( unmodified )implementation of SA model .

Kalyan


All times are GMT -4. The time now is 02:05.