Files
openfoam/applications/solvers/lagrangian/LTSReactingParcelFoam/setPressureWork.H
2011-03-25 12:13:58 +00:00

14 lines
292 B
C

DpDt == dimensionedScalar("zero", DpDt.dimensions(), 0.0);
if (pressureWork)
{
surfaceScalarField phiU("phiU", phi/fvc::interpolate(rho));
DpDt += fvc::div(phiU*fvc::interpolate(p)) - p*fvc::div(phiU);
if (pressureWorkTimeDerivative)
{
DpDt += fvc::ddt(p);
}
}