Incompressible Viscous Turbulent Heat Transfer Energy Equation

Generally, incompressible viscous turbulent heat transfer energy equation is
[latex]\dfrac{\partial T}{\partial t} + \nabla\cdot(\boldsymbol{U} T) =
\alpha \nabla^2 T + \dfrac{\nu}{c_\mathrm{p}} [\nabla \boldsymbol{U} +(\nabla \boldsymbol{U})^\mathrm{T}]:\nabla \boldsymbol{U}
ddt(T) + div(U T) = alpha laplacian(T) + ((nu/cp (grad(U) + (grad(U)).T)) && grad(U))

where [latex]\alpha = \dfrac{k}{\rho c_\mathrm{p}}[/latex]

When turbulent stress considered, the equation becomes
[latex]\dfrac{\partial T}{\partial t} + \nabla\cdot(\boldsymbol{U} T) =
\alpha_\mathrm{eff}\nabla^2 T + \dfrac{\nu_\mathrm{eff}}{c_\mathrm{p}} [\nabla \boldsymbol{U} +(\nabla \boldsymbol{U})^\mathrm{T}]:\nabla \boldsymbol{U}

where [latex]\alpha_\mathrm{eff}= \nu/Pr + \nu_t/Pr_t[/latex], $latex Pr$ is Prandl number, $latex \nu_t$ is turbulent viscosity and $latex Pr_t$ is turbulent Prandl number.

With OpenFOAM, the incompressible viscous turbulence energy equation is written

    volScalarField alphaEff
        turbulence->nu()/Pr + turbulence->nut()/Prt

    volTensorField gradU = fvc::grad(U);

    volTensorField tauEff = turbulence->nuEff() / cp * (gradU + gradU.T());

        // Add energy equation
        fvScalarMatrix TEqn
          + fvm::div(phi, T)
          - fvm::laplacian(alphaEff, T)
          - (tauEff && gradU)


Leave a Reply

Your email address will not be published. Required fields are marked *

To create code blocks or other preformatted text, indent by four spaces:

    This will be displayed in a monospaced font. The first four 
    spaces will be stripped off, but all other whitespace
    will be preserved.
    Markdown is turned off in code blocks:
     [This is not a link](

To create not a block, but an inline code span, use backticks:

Here is some inline `code`.

For more help see


This site uses Akismet to reduce spam. Learn how your comment data is processed.