[TommyM]

Hi All,
I am trying to switch on the "transonic" option in reactingFoam. It should be exist in this solver since I have found it in the solver implementation.
Anyway, when I try to run the solver, it completely ignore the "div(phiv, p)" term in the fvSchemes, so I don't know if I am using it correctly.
Does anyone know something about it?
Thanks in advance.


Tommy

[clapointe]

I am only aware of its usage in the energy equation (below) that is used if internalEnergy is chosen. It is not used in the pressure equation (but phid is).

          he.name() == "e"
          ? fvc::div
            (
                fvc::absolute(phi/fvc::interpolate(rho), U),
                p,
                "div(phiv,p)"
            )
          : -dpdt

Caelan

[TommyM]

Hi Caelan,
Thanks for your reply.
I have chosen sensibleEnthalpy so the "phiv" term does not come into my equations, as you said.
I am trying to use reactingFoam to simulate chemical reactions in a supersonic flow (Ma~1.8).
Anyway, I should switch my energy equation to "internalEnergy" in order to do that (as far as I know, "e" equation is preferable to "h" equation in supersonic flows).
Can you give me your opinion about this, please?

Tommy





P.S.: I know that the PIMPLE-based solvers are not the optimal choice for this kind of problems, but at the moment I have not been able to implement reactions in rhoCentralFoam.

[clapointe]

I do not have much experience using reactingFoam for (tran)sonic problems, so I will refrain from commenting on e vs h thermo. Perhaps someone else could be more help (or you could test yourself). For density-based reacting solvers there are various ones floating around -- e.g. hyFoam, various threads about reacting rhoCentralFoam, etc.

Caelan

[TommyM]

Yes, I already read these threads, in particular reacting rhoCentral, but the implementation of chemical reactions in the central-upwind scheme's perspective is not so straightforward (to me).
Thank you anyway for your help.


Tommy