From c8a63f2da3668549a6d076c12698544abb321b43 Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Thu, 14 Sep 2023 20:52:41 +0100 Subject: [PATCH] multiphaseEuler::phaseSystemSolve: Corrected handling of dilatation effects in the presence of stationary phases --- .../phaseSystems/phaseSystem/phaseSystemSolve.C | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/applications/modules/multiphaseEuler/phaseSystems/phaseSystem/phaseSystemSolve.C b/applications/modules/multiphaseEuler/phaseSystems/phaseSystem/phaseSystemSolve.C index 2484c499d9..4bc83f933b 100644 --- a/applications/modules/multiphaseEuler/phaseSystems/phaseSystem/phaseSystemSolve.C +++ b/applications/modules/multiphaseEuler/phaseSystems/phaseSystem/phaseSystemSolve.C @@ -182,19 +182,19 @@ void Foam::phaseSystem::solve(const PtrList& rAs) ) ); - forAll(movingPhases(), movingPhasej) + forAll(phases(), phasej) { - const phaseModel& phase2 = movingPhases()[movingPhasej]; + const phaseModel& phase2 = phases()[phasej]; const volScalarField& alpha2 = phase2; if (&phase2 != &phase) { - if (phase.divU().valid()) + if (!phase.stationary() && phase.divU().valid()) { dgdt += alpha2()*phase.divU()()(); } - if (phase2.divU().valid()) + if (!phase2.stationary() && phase2.divU().valid()) { dgdt -= alpha()*phase2.divU()()(); }