[Ale22pic]

Hi everybody, it is my first thread, even if I'm always on this forum reading other threads in order to learn something new.

Anyway, I'm working on Openfoam V2112 for my PhD project and more precisely on a solver developed here in my group, called rhoEnergyFoam, which is based on rhoCentralFoam but obviously, there are some modifications.

The solver works perfectly on V6 and V2006 versions, but now we are trying to use it on V2112 too, but I've found that there is a problem in rhoCentralFoam in V2112.

If you do whatever rhoCentralFoam tutorial, and asking to the code to give you in the output the specific internal energy "e" (Info << "e \n" << e << endl; in rhoCentralFoam.C) you'll find that the energy is negative. This value seems to not be a problem for rhoCentralFoam, but for my solver, it is since I got a square root of this value.

Reading the createFields.H in rhoCentralFoam the specific internal energy "e" is defined as follows:

Info<< "Reading thermophysical properties\n" << endl;

autoPtr<psiThermo> pThermo
(
    psiThermo::New(mesh)
);
psiThermo& thermo = pThermo();

volScalarField& e = thermo.he();

Does anybody have the same issues? How can be solved? And also does anybody know how it is defined "e"? I've searched in all the OpenFOAM files.

Thanks in advance for your time.

[Alishaha2]

Hello Alessio,


Probably it is from your reference temperature (Tref).
if you get negative energy values it means for your case the reference temperature is bellow the standard temperature (300 K).

you can fix it in thermophysicalProperties file by setting:


mixture
{
thermodynamics
{
Tref 0; // or etc.
}
}


Best,
Ali