[Maffei]

Hello everyone!

I am programming a solver where I need to reload the compressible turbulent field. I can successfully realize the reading, however, the values diverge when I compare the turbulent read variables (mut and alphat) and the one calculated with the read turbulence field (let's say k and omega). I am sure about the code for calculating the variable mut because I get it from the turbulence model code (just rewrite it In my solver). I can't say the same to alphat but I am using (mut+mu)/Prandtl. Also I have changed the time write to binary with a precision of 30 cases.

The same happens when I create the turbulence object with:

Info<< "Creating turbulence model\n" << endl;
autoPtr<compressible::turbulenceModel> turbulence
(
    compressible::turbulenceModel::New
    (
        rho,
        U,
        phi,
        thermo
    )
);

and made the comparison using turbulence->aphat() and turbulence->mut() after run the turbulence->correct().

For what I programming, I really need to be able to compute the mut and alphat from the turbulent field and be consistent with the calculated primitive field.

Why am I having these diverging values?