[JD_PM]

I want to understand how to make the CavitatingFoam solver (code here: https://cpp.openfoam.org/dev/cavitat...8C_source.html) non-isothermal.

To do so, I am comparing it to a similar solver that actually implements non-isothermal conditions: CompressibleInterFoam (code here: https://cpp.openfoam.org/dev/compres...8C_source.html).

By comparison, I see that the key has to be in how the pressure-velocity PIMPLE corrector loop is written for CompressibleInterFoam i.e. the key has to be in the following piece of code

// --- Pressure-velocity PIMPLE corrector loop
while (pimple.loop())
{
if (pimple.firstPimpleIter() || moveMeshOuterCorrectors)
{
// Store divU from the previous mesh so that it can be mapped
// and used in correctPhi to ensure the corrected phi has the
// same divergence
tmp<volScalarField> divU;

if (correctPhi)
{
// Construct and register divU for mapping
divU = new volScalarField
(
"divU0",
fvc::div(fvc::absolute(phi, U))
);
}

fvModels.preUpdateMesh();

mesh.update();

if (mesh.changing())
{
gh = (g & mesh.C()) - ghRef;
ghf = (g & mesh.Cf()) - ghRef;

MRF.update();

if (correctPhi)
{
#include "correctPhi.H"
}

mixture.correct();

if (checkMeshCourantNo)
{
#include "meshCourantNo.H"
}
}

divU.clear();
}

fvModels.correct();

#include "alphaControls.H"
#include "compressibleAlphaEqnSubCycle.H"

turbulence.correctPhasePhi();

#include "UEqn.H"
#include "TEqn.H"


Might you please hint me where should I focus on studying (moveMeshOuterCorrectors, ...?) in order to understand how to modify CavitatingFoam solver code in order to make it non-isothermal?


Thank you!