CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Problem with LaunderSharma compressible turbulence model wrong formulation (http://www.cfd-online.com/Forums/openfoam-solving/60135-problem-laundersharma-compressible-turbulence-model-wrong-formulation.html)

 luca July 3, 2006 03:39

maybe I'm not right, but in La

maybe I'm not right, but in LaunderSharma compressible model Epsilon equation, we have that :

... + fvm::SuSp((C3 - 2.0/3.0*C1)*rho_*divU, epsilon_)

and in the default turbulenceProperties file we have:

C3 = -0.33

I think that we have to change the default value of C3 as:

C3 = 0.33

or we have to change the formulation of Epsilon eq. as:

... + fvm::SuSp((- C3 - 2.0/3.0*C1)*rho_*divU, epsilon_) or

... - fvm::SuSp(( C3 + 2.0/3.0*C1)*rho_*divU, epsilon_)

because the term C3*rho_*divU has to be positive in RHS of Epsilon equation.

I'm right ?!?

Luca

 mattijs July 3, 2006 05:34

Please post a bug report (in O

Please post a bug report (in OpenFOAM-bugs) if you think it is incorrect.

 hjasak July 3, 2006 05:37

Nope - this is to do with the

Nope - this is to do with the C3 term in engines simulations, where it turns out that a "small positive constant" C3 actually needs to be negative to get the right behaviour. Find an old-ish engines modelling Thesis from Imperial College and you will find the details.

Hrv

 luca July 3, 2006 06:21

so why we have different formu

so why we have different formulation in epsilon Eqn. between LaunderSharmaKE and kEpsilon turbulence model, (except f2() term)

kEpsilon

C1*G*epsilon_/k_
- fvm::SuSp(((2.0/3.0)*C1 + C3)*rho_*divU, epsilon_)
- fvm::Sp(C2*rho_*epsilon_/k_, epsilon_)

LaunderSharmaKE

C1*G*epsilon_/k_
+ fvm::SuSp((C3 - 2.0/3.0*C1)*rho_*divU, epsilon_)
- fvm::Sp(C2*f2()*rho_*epsilon_/k_, epsilon_)

Luca

 hjasak July 3, 2006 06:24

In that case: hmm. What does

In that case: hmm. What does the original paper say?

Hrv

 All times are GMT -4. The time now is 12:49.