diff --git a/applications/modules/multiphaseEuler/phaseSystems/phaseInterface/phaseInterface/phaseInterface.C b/applications/modules/multiphaseEuler/phaseSystems/phaseInterface/phaseInterface/phaseInterface.C index 2b8a73021d..06a38a8ff7 100644 --- a/applications/modules/multiphaseEuler/phaseSystems/phaseInterface/phaseInterface/phaseInterface.C +++ b/applications/modules/multiphaseEuler/phaseSystems/phaseInterface/phaseInterface/phaseInterface.C @@ -508,7 +508,18 @@ Foam::tmp Foam::phaseInterface::rho() const Foam::tmp Foam::phaseInterface::magUr() const { - return mag(phase1().U() - phase2().U()); + if (phase1().stationary()) + { + return mag(phase2().U()); + } + else if (phase2().stationary()) + { + return mag(phase1().U()); + } + else + { + return mag(phase1().U() - phase2().U()); + } }