From 1e79d0245c45c763b9bc044ce045c7b6c5c3ce5a Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Thu, 27 Jul 2023 10:54:52 +0100 Subject: [PATCH] MomentumTransferPhaseSystem::partialElimination: Removed the phase velocity mean adjustment The phase velocity mean adjustment was introduced for consistency with phase flux mean adjustment which is necessary to ensure the mean flux divergence is preserved. However for systems with very high drag it has proved preferable to not adjust the velocity of the phases to conserve momentum rather than ensure consistency with the fluxes. --- .../MomentumTransferPhaseSystem.C | 21 ++++--------------- 1 file changed, 4 insertions(+), 17 deletions(-) diff --git a/applications/modules/multiphaseEuler/phaseSystems/PhaseSystems/MomentumTransferPhaseSystem/MomentumTransferPhaseSystem.C b/applications/modules/multiphaseEuler/phaseSystems/PhaseSystems/MomentumTransferPhaseSystem/MomentumTransferPhaseSystem.C index 4f71dc00b2..a732ec4be3 100644 --- a/applications/modules/multiphaseEuler/phaseSystems/PhaseSystems/MomentumTransferPhaseSystem/MomentumTransferPhaseSystem.C +++ b/applications/modules/multiphaseEuler/phaseSystems/PhaseSystems/MomentumTransferPhaseSystem/MomentumTransferPhaseSystem.C @@ -1344,22 +1344,6 @@ void Foam::MomentumTransferPhaseSystem::partialElimination phaseSystem::phaseModelList& phases = this->phaseModels_; - // Calculate the mean velocity from the current velocity - // of the moving phases - volVectorField Um(this->movingPhases()[0]*this->movingPhases()[0].U()); - - for - ( - label movingPhasei=1; - movingPhaseimovingPhases().size(); - movingPhasei++ - ) - { - Um += - this->movingPhases()[movingPhasei] - *this->movingPhases()[movingPhasei].U(); - } - // Remove the drag contributions from the velocity and flux of the phases // in preparation for the partial elimination of these terms forAll(phases, i) @@ -1509,7 +1493,8 @@ void Foam::MomentumTransferPhaseSystem::partialElimination } } - this->setMixtureU(Um); + // Adjust the phase-fluxes such that the mean flux + // obtained from the pressure solution is retained this->setMixturePhi(alphafs, this->phi()); } @@ -1643,6 +1628,8 @@ void Foam::MomentumTransferPhaseSystem::partialEliminationf } } + // Adjust the phase-fluxes such that the mean flux + // obtained from the pressure solution is retained this->setMixturePhi(alphafs, this->phi()); }