[Tobi]

Hi all,

I am trying to find the place where the temperature is derived from the enthalpy.

- Solver: rhoSimpleFoam
- thermodynamics: hPsiThermo
- transport: hConstTransport

Therefor I had a look at the solver and the enthalpy equation h:

{
    fvScalarMatrix hEqn
    (
        fvm::div(phi, h)
      - fvm::Sp(fvc::div(phi), h)
      - fvm::laplacian(turbulence->alphaEff(), h) 
     ==
      - fvc::div(phi, 0.5*magSqr(U), "div(phi,K)")
    );

    hEqn.relax();

    hEqn.solve();

    thermo.correct();

Here I think the step I am searching is in the function "correct".
That leads me to the thermomodel -> hPsiThermo.C.

From the correct() function you get into the calculate() function. And there (I think) you get the temperature derived from the enthalpy. Its described in the *.H file that T is calculated there.

I found that function:

        TCells[celli] = mixture_.TH(hCells[celli], TCells[celli]);

And would find the definition of that function. With grep I get into the specieThermo.H file. Here I found that:

        // Energy->temperature  inversion functions

            //- Temperature from Enthalpy given an initial temperature T0
            inline scalar TH(const scalar H, const scalar T0) const;

So I am very sure that I am on the right way. But thats just a function prototype and not a function declaration. I searched everywhere but do not find the definition of that function.

Where can i find it?
Can someone give me an advice?

Thanks in advance
Tobi

[Tobi]

I should have a look at the specieThermoI.H
So solved Thanks for reading

[mturcios777]

If you look in specieThermoI.H you will find the function TH. Looking at it you will see that the template is calling another function T and feeding it the proper coefficients (whether its enthalpy or energy you are inverting).

The function T is also in the file, near the very beginning. Here is where temperature is obtained by inverting the enthalpy/energy based on the the species thermo (its a simple Newton's method if I read the code right).

EDIT: Hehe Tobi, you must have been thinking while I was typing ;-)

[Tobi]

Hihi,

well, nevertheless thanks for reading and anwering
You are right. After I read the specieThermoI.H file I saw that TH is calling a other function T() and thats defined in the beginning of the thermo model

Wow .. cool ... hihi. I am getting familiar and familiar with OF code.