[er99]

Hi,

I want to use reactingFoam in OpenFOAM 4.x while accounting for the effects of gravity. To do this, I have chosen to include the source terms offered by fvOptions buyoancyForce and buoyancyEnergy.

Example usage for the buyoancyForce option is:
fields (U); // Name of velocity field

Example usage for the buyoancyEnergy option is:
fields (h); // Name of energy field

I used U for the velocity field and he for the energy field. The latter was chosen because in reactingFoam's energy equation (EEqn.H) you see an extra fvOptions term with rho and he: + fvOptions(rho, he).

However, when running the code, I got this error message in the logfile:

--> FOAM Warning :
From function virtual void Foam::fv:: option::checkApplied() const
in file cfdTools/general/fvOptions/fvOption.C at line 118
Source gravityEnergyEqtn defined for field he but never used


The code ran to the end but since I got the above error message, I'm assuming this means that the gravitational term in the energy equation has not been solved for and therefore my results are not the ones I need.

I also tried with field "e" but this gave me the same error message. However, when I use the field "h" it works. This is confusing to me as reactingFoam's energy equation has the line:
he.name() == "e"


Can someone help me understand why it works when "h" is selected as the field?

Thanks in advance!

[LuckyTran]

he is not a valid field. You'll find he in solvers because the energy equation can be written in h or e. If you dig deeper into the codes you'll find condition statements that handle cases when the fieldname is h or e.

Go to constant/thermophysialProperties dict and check the entry for energy. If it says sensiblEnthalpy then you are using an energy equation based on h. If it says sensibleInternalEnergy then you can use e.