DPMDyMFoam, DPMDyMFoam: Corrected support for closed-domain simulations

Also added support for extrapolated pressure boundary conditions.

applications/solvers/lagrangian/DPMFoam/DPMDyMFoam/pEqn.H

@@ -1,6 +1,5 @@
 {
-    volVectorField HbyA("HbyA", Uc);
+    volVectorField HbyA(constrainHbyA(rAUc*UcEqn.H(), Uc, p));
-    HbyA = rAUc*UcEqn.H();
 
     surfaceScalarField phiHbyA
     (
@@ -8,7 +7,6 @@
         (
            fvc::flux(HbyA)
          + alphacf*rAUcf*fvc::ddtCorr(Uc, Ucf)
-         + phicForces
         )
     );
 
@@ -19,6 +17,8 @@
         fvc::makeAbsolute(phiHbyA, Uc);
     }
 
+    phiHbyA += phicForces;
+
     // Update the pressure BCs to ensure flux consistency
     constrainPressure(p, Uc, phiHbyA, rAUcf);
 

applications/solvers/lagrangian/DPMFoam/pEqn.H

@@ -1,6 +1,5 @@
 {
-    volVectorField HbyA("HbyA", Uc);
+    volVectorField HbyA(constrainHbyA(rAUc*UcEqn.H(), Uc, p));
-    HbyA = rAUc*UcEqn.H();
 
     surfaceScalarField phiHbyA
     (
@@ -8,10 +7,16 @@
         (
            fvc::flux(HbyA)
          + alphacf*rAUcf*fvc::ddtCorr(Uc, phic)
-         + phicForces
         )
     );
 
+    if (p.needReference())
+    {
+        adjustPhi(phiHbyA, Uc, p);
+    }
+
+    phiHbyA += phicForces;
+
     // Update the pressure BCs to ensure flux consistency
     constrainPressure(p, Uc, phiHbyA, rAUcf);