From e8d45b1d3c68d04486cb20dfbf83580ef02e512b Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Fri, 3 Jul 2015 18:04:34 +0100 Subject: [PATCH] reactingTwoPhaseEulerFoam: Replace x*pos(x) with posPart(x) and same for neg --- .../multiphase/reactingTwoPhaseEulerFoam/pU/UEqns.H | 4 ++-- .../HeatAndMassTransferPhaseSystem.C | 8 ++++---- 2 files changed, 6 insertions(+), 6 deletions(-) diff --git a/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/pU/UEqns.H b/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/pU/UEqns.H index 49fa1bd63..82b929548 100644 --- a/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/pU/UEqns.H +++ b/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/pU/UEqns.H @@ -11,8 +11,8 @@ volScalarField Kd(fluid.Kd()); { volScalarField Vm(fluid.Vm()); - const volScalarField dmdt12(fluid.dmdt()*pos(fluid.dmdt())); - const volScalarField dmdt21(fluid.dmdt()*neg(fluid.dmdt())); + const volScalarField dmdt12(posPart(fluid.dmdt())); + const volScalarField dmdt21(negPart(fluid.dmdt())); { U1Eqn = diff --git a/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/phaseSystems/PhaseSystems/HeatAndMassTransferPhaseSystem/HeatAndMassTransferPhaseSystem.C b/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/phaseSystems/PhaseSystems/HeatAndMassTransferPhaseSystem/HeatAndMassTransferPhaseSystem.C index b95f81aa0..71c5a0672 100644 --- a/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/phaseSystems/PhaseSystems/HeatAndMassTransferPhaseSystem/HeatAndMassTransferPhaseSystem.C +++ b/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/phaseSystems/PhaseSystems/HeatAndMassTransferPhaseSystem/HeatAndMassTransferPhaseSystem.C @@ -199,8 +199,8 @@ Foam::HeatAndMassTransferPhaseSystem::momentumTransfer() const const volVectorField& U2(pair.phase2().U()); const volScalarField dmdt(this->dmdt(pair)); - const volScalarField dmdt12(dmdt*pos(dmdt)); - const volScalarField dmdt21(dmdt*neg(dmdt)); + const volScalarField dmdt12(posPart(dmdt)); + const volScalarField dmdt21(negPart(dmdt)); *eqns[pair.phase1().name()] += fvm::Sp(dmdt21, U1) - dmdt21*U2; *eqns[pair.phase2().name()] += dmdt12*U1 - fvm::Sp(dmdt12, U2); @@ -315,8 +315,8 @@ Foam::HeatAndMassTransferPhaseSystem::heatTransfer() const const volScalarField& K2(phase2.K()); const volScalarField dmdt(this->dmdt(pair)); - const volScalarField dmdt12(dmdt*pos(dmdt)); - const volScalarField dmdt21(dmdt*neg(dmdt)); + const volScalarField dmdt12(posPart(dmdt)); + const volScalarField dmdt21(negPart(dmdt)); const volScalarField& Tf(*Tf_[pair]); *eqns[phase1.name()] +=