[sazabi2001]

Hi:

I'm now working on natural convection problem and I'd like to try changing the reference density 'rho0' during running.

I know it will be tricky, since changing the 'rho0' will mean changing the Enthalpy state of the whole problem. But before looking into the mathematics, I'd like to know whether it is possible, or in a manageable difficult that I can read and change the 'rho0' in the solver (maybe by building an extra rho0Eqn.H) during ruining.

For example, i'd like to use the average density calculated during real-time solving to replace the static 'rho0' I put in the thermophysicalProperties

Now the problem is: the rho0 value sits in the thermophysicalProperties / equationOfState, and some other things (thermo.H basicThermo.H rhoThermo.h heThermo.H and heRhoThermo.h) will read it. can I construct it like a uniformDimensionedScalarField in the [0] folder, read it and make it affect those thermo.H programs?

Regards