From 26f2ce9c9e1339459dc17c7f2cb382f6108cebd6 Mon Sep 17 00:00:00 2001 From: Will Bainbridge Date: Wed, 14 Mar 2018 17:06:35 +0000 Subject: [PATCH] reactingMultiphaseEulerFoam: Bug fix to divergence handling --- .../multiphaseSystem/multiphaseSystem.C | 15 +++------------ 1 file changed, 3 insertions(+), 12 deletions(-) diff --git a/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/multiphaseSystem/multiphaseSystem.C b/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/multiphaseSystem/multiphaseSystem.C index e794c47031..93d7ec3fb5 100644 --- a/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/multiphaseSystem/multiphaseSystem.C +++ b/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/multiphaseSystem/multiphaseSystem.C @@ -189,8 +189,6 @@ void Foam::multiphaseSystem::solveAlphas() ); - volScalarField divU(fvc::div(fvc::absolute(phi_, phases().first().U()))); - forAll(phases(), phasei) { phaseModel& phase = phases()[phasei]; @@ -209,20 +207,13 @@ void Foam::multiphaseSystem::solveAlphas() mesh_ ), mesh_, - dimensionedScalar("Sp", divU.dimensions(), 0.0) + dimensionedScalar("Sp", dimless/dimTime, 0) ); volScalarField::Internal Su ( - IOobject - ( - "Su", - mesh_.time().timeName(), - mesh_ - ), - // Divergence term is handled explicitly to be - // consistent with the explicit transport solution - divU*min(alpha, scalar(1)) + "Su", + min(alpha, scalar(1))*fvc::div(fvc::absolute(phi_, phase.U())) ); if (phase.divU().valid())