BUG: Adding relative flux to the mesh

2021-02-15 12:32:41 -08:00 · 2021-02-15 12:32:41 -08:00 · 0990d30b73
commit 0990d30b73
parent 48fdec46c7
1 changed files with 6 additions and 1 deletions
--- a/applications/solvers/multiphase/interPhaseChangeFoam/alphaEqnSubCycle.H
+++ b/applications/solvers/multiphase/interPhaseChangeFoam/alphaEqnSubCycle.H
@ -9,7 +9,12 @@
        phic += (interface.cAlpha()*icAlpha)*fvc::interpolate(mag(U));
    }

-    volScalarField divU(fvc::div(phi));
+    volScalarField divU
+    (
+        mesh.moving()
+      ? fvc::div(phi + mesh.phi())
+      : fvc::div(phi)
+    );

    if (nAlphaSubCycles > 1)
    {