[MamboJambo]

Hello to all,


In compressibleInterFoam of OF2006, the temperature equation is defined as:

    fvScalarMatrix TEqn
    (
        fvm::ddt(rho, T) + fvm::div(rhoPhi, T) - fvm::Sp(contErr, T)
      - fvm::laplacian(turbulence.alphaEff(), T)
      + (
             divUp()// - contErr/rho*p
          +  (fvc::ddt(rho, K) + fvc::div(rhoPhi, K))() - contErr*K
        )
       *(
           alpha1()/mixture.thermo1().Cv()()
         + alpha2()/mixture.thermo2().Cv()()
        )
     ==
        fvOptions(rho, T)
    );

Afterwards two functions are called:

    mixture.correctThermo();
     mixture.correct();

The first function will create a he field from the values of pressure and temperature for each phase and call the correct() function.
In the correct() function the temperature is being recalculated from the energy from line 80-88 in the heRhoThermo.C with a Newton-Raphson method [https://caefn.com/openfoam/temperature-calculation]

My question is: If we are already solving for temperature, why do we need to recalculate the values from energy?

Edit:
The temperature used for the calculation is from obtained from the mixture object which will get the temperature field by inheriting from the basicThermo class.

This object has two AutoPtr pointers for each phase. Each of which will create they own temperature field since they inherit from basicThermo. Are the lines 80-88 just updating the temperature values of these fields? By why from energy and not directly from the calculated temperature?

Or am I getting the wrong idea?