[agustinvo]

Hi everyone,


I was looking at the divq term of the ThermophysicalTransportModels library and I don't understand the current implementation:

// Return heat flux source as an implicit energy correction
// to the temperature gradient flux
return
    -correction(fvm::laplacian(this->alpha()*this->alphaEff(), he))
    -fvc::laplacian(this->alpha()*this->kappaEff(), this->thermo().T());

I have some questions regarding this:

1. If [math]\alpha_{eff}[\math] is usually employed (as far as I have used) as the effective diffusion of temperature, why is it employed in the laplacian term of he?
2. What does the correction function do?
3. I suppose there is a sort of stability reason to have a fvm and a fvc term, however I don't see the point on it here.

[CorbinMG]

Hi Agustín,

In regards to (2), Some mysteries (but not all) of the correction() function are explained here: Foam::correction() function

Best regards,
Corbin