[chiven]

I am wondering why the temperature equation of buoyantBoussinesqPisoFoam/Simple is considering both the laminar and turbulent Prandtl numbers attached in the follows.
According to the energy conservation equation, it should be enough only considering the turbulent one for a turbulent case.
Can anyone give me an explanation? Thanks in advance.
Chiven

{
volScalarField kappaEff
(
"kappaEff",
turbulence->nu()/Pr + turbulence->nut()/Prt
);
fvScalarMatrix TEqn
(
fvm::ddt(T)
+ fvm::div(phi, T)
- fvm::laplacian(kappaEff, T)
);
TEqn.relax();
TEqn.solve();
rhok = 1.0 - beta*(T - TRef);
}