Tubrulence model diffusion term ???

Tobi · June 5, 2012, 06:37

hi all together,

i have a question to the RNG-kEpsilon equation respectively to the diffusion term.

Well in the code there is the following equation:

Code:

// Dissipation equation
    tmp<fvScalarMatrix> epsEqn
    (
        fvm::ddt(rho_, epsilon_)
      + fvm::div(phi_, epsilon_)
      - fvm::laplacian(DepsilonEff(), epsilon_)
      ==
        (C1_ - R)*G*epsilon_/k_
      - fvm::SuSp(((2.0/3.0)*C1_ + C3_)*rho_*divU, epsilon_)
      - fvm::Sp(C2_*rho_*epsilon_/k_, epsilon_)
    );

especially to that code!!!

Code:

- fvm::laplacian(DepsilonEff(), epsilon_)

the diffusivity is calculated in the header file like:

Code:

mut_/sigmaEps_ + mu()

Thats all clear but I have that equation in my literature (pdf)
and (please correct me) if I set the diffusion term to the left and transform it I get:

Code:

mu() - mut_/sigmaEps

Hope I asked clearly.
Thx in advance,
tobi

eysteinn · June 5, 2012, 07:22

Quote:

Originally Posted by Tobi

hi all together,

i have a question to the RNG-kEpsilon equation respectively to the diffusion term.

Well in the code there is the following equation:

Code:

// Dissipation equation
    tmp<fvScalarMatrix> epsEqn
    (
        fvm::ddt(rho_, epsilon_)
      + fvm::div(phi_, epsilon_)
      - fvm::laplacian(DepsilonEff(), epsilon_)
      ==
        (C1_ - R)*G*epsilon_/k_
      - fvm::SuSp(((2.0/3.0)*C1_ + C3_)*rho_*divU, epsilon_)
      - fvm::Sp(C2_*rho_*epsilon_/k_, epsilon_)
    );

especially to that code!!!

Code:

- fvm::laplacian(DepsilonEff(), epsilon_)

the diffusivity is calculated in the header file like:

Code:

mut_/sigmaEps_ + mu()

Thats all clear but I have that equation in my literature (pdf)
and (please correct me) if I set the diffusion term to the left and transform it I get:

Code:

mu() - mut_/sigmaEps

Hope I asked clearly.
Thx in advance,
tobi

Hmm . . . my literature says nu+nu_t/sigma_epsilon not nu - nu_t/sigma_epsilon for the diffusion term in the epsilon eq..

Tobi · June 5, 2012, 08:21

Quote:

Originally Posted by eysteinn

Hmm . . . my literature says nu+nu_t/sigma_epsilon not nu - nu_t/sigma_epsilon for the diffusion term in the epsilon eq..

Well thx for your replay!
I see that in the diffusion term for k its like " nu + nu_t/sigma_k"

I ll have a look at my bible "Ferziger"

"May be its just a mistake while copying things

"

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
fv3 term in Spalart Allmaras model	bfa	OpenFOAM Programming & Development	1	July 16, 2013 14:02
Include Viscous Work term (in fluid model)	mehshan.javaid	Main CFD Forum	0	September 22, 2011 11:43
Buoyancy term in Spalart-Allmaras model	nickvinn	Main CFD Forum	0	November 17, 2010 09:41
Diffusion Flames with Flamelet model ,little question.	AdidaKK	Main CFD Forum	11	November 3, 2009 03:01
How to model diffusion	Aboulard	CFX	3	May 7, 2007 10:01