From 3516b97a4dc671707f9f5cc62b8d0eb7ed11e454 Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Mon, 4 Sep 2023 13:43:25 +0100 Subject: [PATCH] multiphaseEuler::setRDeltaT: Loop over the moving phases only to obtain maxPhi --- .../modules/multiphaseEuler/multiphaseEuler/setRDeltaT.C | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/applications/modules/multiphaseEuler/multiphaseEuler/setRDeltaT.C b/applications/modules/multiphaseEuler/multiphaseEuler/setRDeltaT.C index 44edde958c..cfd37ec593 100644 --- 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