applications/solvers/lagrangian/LTSReactingParcelFoam/setPressureWork.H

   1 DpDt == dimensionedScalar("zero", DpDt.dimensions(), 0.0);
   2
   3 if (pressureWork)
   4 {
   5     surfaceScalarField phiU("phiU", phi/fvc::interpolate(rho));
   6
   7     DpDt += fvc::div(phiU*fvc::interpolate(p)) - p*fvc::div(phiU);
   8
   9     if (pressureWorkTimeDerivative)
  10     {
  11         DpDt += fvc::ddt(p);
  12     }
  13 }