0000751: rhoCentralFoam: diffusion coefficients in the two last terms of energy corrector equation

Description

In my opinion, the thermal diffusion coefficients in the two last terms of the diffusive energy corrector equation in rhoCentralFoam should both be effective. If I understood correctly, the diffusion would be in principle fully taken into account by the first diffusion term. The last two terms are there to drive consistency between energy and temperature, so that e --> CvT + 0.5magSqr(U), and the two terms would cancel each other out in a fully converged solution ( e=Cv*T 0.5magSqr(U) ). Am I correct in this?

If so, shouldn't the second diffusion term have a diffusion coefficient alphaEff (turbulent + laminar) instead of alpha (only laminar)? At the moment the thermal conductivity coefficient (k) is turbulent + laminar. Here's the code:

If I replace the alpha with alphaEff (or comment out the two last terms) in my modified rhoCentralFoam with specie transport for a system I'm considering, I get temperatures consistent with a reactingFoam solution. Am I missing something, or is there a bug?

There is certainly an issue with the code, either the k should be laminar only and the turbulent part handled by e or if the k is actuall kEff then the alpha sholud be alphaEff to cancel. So the improntant question is what energy/temperature does the turbulence "diffuse"?