multiphaseEuler::setRDeltaT: Loop over the moving phases only to obtain maxPhi

2023-09-04 13:43:25 +01:00
parent 58f82fccbd
commit 3516b97a4d
1 changed files with 2 additions and 2 deletions
--- a/applications/modules/multiphaseEuler/multiphaseEuler/setRDeltaT.C
+++ b/applications/modules/multiphaseEuler/multiphaseEuler/setRDeltaT.C
@ -47,9 +47,9 @@ void Foam::solvers::multiphaseEuler::setRDeltaT()

    surfaceScalarField maxPhi("maxPhi", phi);

-    forAll(phases, phasei)
+    forAll(movingPhases, movingPhasei)
    {
-        maxPhi = max(maxPhi, mag(phases[phasei].phi()));
+        maxPhi = max(maxPhi, mag(movingPhases[movingPhasei].phi()));
    }

    // Set the reciprocal time-step from the local Courant number