[c_underwood]

I am trying to create a new solver which is effectively rhoCentralFoam, but with a way of giving a min energy to the simulation.
The reason I want to do this is because I am trying to simulate a Gas Nozzle and the temperature keeps dropping below 0K, which causes a crash.

I found on here:
Instability solving low-density plume expansion with rhoCentralFoam
where changing line 218 in openfoam6 from:
e = rhoE/rho - 0.5*magSqr(U);
to
e = max(rhoE/rho - 0.5*magSqr(U),minE);

I now want to be able to add minE to the control dict, and then I will be able to enter my value in there. However, I can't work out how the control dict is read in, and whether I can make a modification for just my new solver?

If someone could point me in the right direction for adding this minE term and reading it into my new solver that would be great.

[raumpolizei]

Hello c_underwood,

it is strange that your temperature goes below 0 K, are you getting negative energies or something alike in your domain? I don't think it's a too good idea to bind the energy or temperature in such a manner. Nonetheless, regarding your issue, you could do something like this (assuming thermo is a reference to your thermodynamics object in your simulation):

// In your createFields.H - of1612
dimensionedScalar minEnergy(static_cast<dictionary>(thermo).lookup("eMin"));

the entry eMin will have to be specified in your case/constant/thermodynamicProperties file, which is in my opinion better than the controlDict.

[permianchain]

PermianChain Technologies Inc. is operating the world’s first Blockchain-as-a-Service (BaaS) platform for natural resources, enhancing transparency, efficiency and profitability for our stakeholders. We discover the effective ways to harness the blockchain’s potential, using data science and AI to digitise, tokenize and capitalise on proven but undeveloped natural resources. If you’d like to hear more about how PermianChain is enhancing oil and gas trading and investment, join our mailing list today.