[Maffei]

Hello everyone!

I am writing a solver (let's name it as firstSolver) where I get the initial variables from another solver (let's name it as primalSolver). The primalSolver is called using inside the firstSolver:

system(myPrimalSolver);

Then, the primalSolver saves in the current runTime and, forcing the last runtime value in runTime object, the firstSolver reads what was calculated in the primalSolver. I know that it is not the best practice, however, using both solvers in one will create a problem that I pretend to solve in the future.

For U, p and phi variables, everything works fine. The problem lies in the FluidThemo object that I can not force to update to the loaded variables. So, when I calculate things using:

thermo.T()
thermo.he()
thermo.rho()

I have different values from what was calculated in the primalSolver. The best that I could do is force the variables by doing:

volScalarField TRead (IOobject("T", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);

thermo.T() = Tread;
thermo.he() = Tread*thermo.Cp();
thermo.rho() = thermo.p()/(Rgas*thermo.T())

But I still have differences from the reader field around 1e-1. This increases the residuals from my equation which is terrible for me because I need this information very precisely.

Does anyone have any idea how can I force the values in the thermo class?