sonicLiquidFoam: Updated pressure equation to use Dp

2025-11-28 03:28:01 +00:00 · 2012-05-03 15:14:32 +01:00
parent f09f545dfb
commit 2e09a6628e
1 changed files with 2 additions and 1 deletions
--- a/applications/solvers/compressible/sonicFoam/sonicLiquidFoam/sonicLiquidFoam.C
+++ b/applications/solvers/compressible/sonicFoam/sonicLiquidFoam/sonicLiquidFoam.C
@ -85,13 +85,14 @@ int main(int argc, char *argv[])
            );

            phi = (rhoO/psi)*phid;
+            volScalarField Dp("Dp", rho*rAU);

            fvScalarMatrix pEqn
            (
                fvm::ddt(psi, p)
              + fvc::div(phi)
              + fvm::div(phid, p)
-              - fvm::laplacian(rho*rAU, p)
+              - fvm::laplacian(Dp, p)
            );

            pEqn.solve();