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

"

Thread Tools	Search this Thread
Show Printable Version Email this Page	Search this Thread: Advanced Search
Display Modes
Linear Mode Switch to Hybrid Mode Switch to Threaded Mode

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