[Daedalus]

Hello!


This chunk of code represents viscous correction for energy equation in rhoCentralFoam solver in OpenFOAM-dev.

if (!inviscid)
{ 

solve
(
fvm::ddt(rho, e) - fvc::ddt(rho, e)
  - fvm::laplacian(turbulence->alphaEff(), e)  

);
thermo.correct();
rhoE = rho*(e + 0.5*magSqr(U));
}

Here, turbulence->alphaEff() is defined as kappaEff/Cp and using energy definition as sensibleInternalEnergy in constant/thermophysicalProperties file, e=Cv*T.
So, above equation effectively translates into




But, it should be only

.


Since, it would not make much difference if is constant. But when using Janaf tables at high temperature, these two gives different solutions because varies with temperature. Which one is correct? Can somebody help?


I have same question regarding calculation of wallHeatFlux. Should it be

calcWallHeatFlux(turbModel.alphaEff(), turbModel.transport().he())
or
calcWallHeatFlux(turbModel.kappaEff(), turbModel.transport().T())


in file https://github.com/OpenFOAM/OpenFOAM...wallHeatFlux.C at line 220



Thanks