From f3fac5bbffd61c75431cba372da029ac776685d1 Mon Sep 17 00:00:00 2001 From: Henry Date: Wed, 14 Mar 2012 18:18:48 +0000 Subject: [PATCH] multiphaseEulerFoam: Improved consistency between implicit and explicit parts of the virtual mass --- .../solvers/multiphase/multiphaseEulerFoam/UEqns.H | 13 +++++++++---- 1 file changed, 9 insertions(+), 4 deletions(-) diff --git a/applications/solvers/multiphase/multiphaseEulerFoam/UEqns.H b/applications/solvers/multiphase/multiphaseEulerFoam/UEqns.H index bce95f7d5a..b612e6b5d9 100644 --- a/applications/solvers/multiphase/multiphaseEulerFoam/UEqns.H +++ b/applications/solvers/multiphase/multiphaseEulerFoam/UEqns.H @@ -15,12 +15,17 @@ forAllIter(PtrDictionary, fluid.phases(), iter) phasei, new fvVectorMatrix ( - (scalar(1) + fluid.Cvm(phase)/phase.rho())* + fvm::ddt(alpha, U) + + fvm::div(phase.phiAlpha(), U) + - fvm::Sp(fvc::ddt(alpha) + fvc::div(phase.phiAlpha()), U) + + + (alpha/phase.rho())*fluid.Cvm(phase)* ( - fvm::ddt(alpha, U) - + fvm::div(phase.phiAlpha(), U) - - fvm::Sp(fvc::ddt(alpha) + fvc::div(phase.phiAlpha()), U) + fvm::ddt(U) + + fvm::div(phase.phi(), U) + - fvm::Sp(fvc::div(phase.phi()), U) ) + - fvm::laplacian(alpha*nuEff, U) - fvc::div (