[jameswilson620]

All,

In buoyantPimpleFoamSolver a few lines in EEqn.H read:

            fvc::div
            (
                fvc::absolute(phi/fvc::interpolate(rho), U),
                p,
                "div(phiv,p)"
            )

I understand dividing by the interpolated surfaceScalarField rhof to get the volume flux for the divergence operator but why didn't we stop there?

For example:

fvc::div(phi/fvc::interpolate(rho),p) .... (1)

In other words, what is the provided code doing that cannot be accomplished by (1)? Or, what is the purpose of fvc::absolute() (what are it's arguments and it's output)?

Also, I suppose providing "div(phiv,p)" in the argument of fvc::div() in the solver gives a name to reference in fvSchemes? I have a question in another post related to this issue: http://www.cfd-online.com/Forums/mai...ansitions.html

Thanks, James

[alexeym]

Hi,

1. fvc::absolute does nothing for static mesh, see http://foam.sourceforge.net/docs/cpp...9881bf6671d8bf.

2. as surfaceScalarField fvc::absolute(phi/fvc::interpolate(rho), U) has rather fancy name it is necessary explicitly set the name of operation, so fvc::div can decide what scheme to use for discretization.