[harshad88]

Dear all,

I am trying to develop a new solver for premixed flame simulations using reactingFoam. I have a question about the EEqn.H file in the default code. It reads:

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

What is the meaning of fvc::absolute(phi/fvc::interpolate(rho), U), and the line : -dpdt ???

Your help would be really appreciated. Thanks

[wouter]

Hello harshad88,

I am no expert.
It is part of a C or C++ construct A?b:c meaning
if A is true than b else c in this case c = "-dpdt"

hope this helps
Wouter

[harshad88]

Thanks for your reply Wouter.

Could you also tell me what the meaning of fvc::absolute(phi/fvc::interpolate(rho), U) is? I cant understand what this terms corresponds to in the energy equation.

Thanks

[wouter]

As i wrote I am no expert but I think this means phiv because of the "div(phiv,p)" but I have no clue what this means.

Hope someone else can take over
Wouter

[alundilong]

Originally Posted by harshad88

Thanks for your reply Wouter.

Could you also tell me what the meaning of fvc::absolute(phi/fvc::interpolate(rho), U) is? I cant understand what this terms corresponds to in the energy equation.

Thanks

fvc::absolute(phi/fvc::interpolate(rho), U) return absolute face flux.
if the mesh is moving, it will add up the contribution of moving mesh,
which mean the face flux is
phi/fvc::interpolate(rho) + fvc::meshPhi(U)
;
if static mesh, it will return phi/fvc::interpolate(rho)

[Dan1788]

Thanks Yijin, That was really helpful

[massive_turbulence]

Originally Posted by wouter

Hello harshad88,

I am no expert.
It is part of a C or C++ construct A?b:c meaning
if A is true than b else c in this case c = "-dpdt"

hope this helps
Wouter

Simply called the "ternary operator" as part of C++.