From 5571dbf107cd6860aee7d5484001105e57063e12 Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Thu, 7 Jul 2016 14:36:10 +0100 Subject: [PATCH] reactingTwoPhaseEulerFoam/pUf/UEqns: Changed naming convention for the dmdt's for consistency with HeatAndMassTransferPhaseSystem.C Resolves bug-report http://bugs.openfoam.org/view.php?id=2141 --- .../reactingTwoPhaseEulerFoam/pUf/UEqns.H | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/pUf/UEqns.H b/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/pUf/UEqns.H index 00a0fb8cfe..546c369e5a 100644 --- a/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/pUf/UEqns.H +++ b/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/pUf/UEqns.H @@ -18,14 +18,14 @@ fvVectorMatrix U2Eqn(U2, rho2.dimensions()*U2.dimensions()*dimVol/dimTime); + MRF.DDt(U2) ); - const volScalarField dmdt12(posPart(fluid.dmdt())); - const volScalarField dmdt21(negPart(fluid.dmdt())); + const volScalarField dmdt21(posPart(fluid.dmdt())); + const volScalarField dmdt12(negPart(fluid.dmdt())); { U1Eqn = ( fvm::div(alphaRhoPhi1, U1) - fvm::Sp(fvc::div(alphaRhoPhi1), U1) - + fvm::Sp(dmdt12, U1) - dmdt12*U2 + + fvm::Sp(dmdt21, U1) - dmdt21*U2 + MRF.DDt(alpha1*rho1, U1) + phase1.turbulence().divDevRhoReff(U1) + Vm*(UgradU1 - (UgradU2 & U2)) @@ -40,7 +40,7 @@ fvVectorMatrix U2Eqn(U2, rho2.dimensions()*U2.dimensions()*dimVol/dimTime); U2Eqn = ( fvm::div(alphaRhoPhi2, U2) - fvm::Sp(fvc::div(alphaRhoPhi2), U2) - - fvm::Sp(dmdt21, U2) + dmdt21*U1 + - fvm::Sp(dmdt12, U2) + dmdt12*U1 + MRF.DDt(alpha2*rho2, U2) + phase2.turbulence().divDevRhoReff(U2) + Vm*(UgradU2 - (UgradU1 & U1))