zhyang-dev/OpenFOAM-6
1
0
Fork 0
mirror of https://github.com/OpenFOAM/OpenFOAM-6.git synced 2025-12-08 06:57:46 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

reactingTwoPhaseEulerFoam: Replace x*pos(x) with posPart(x) and same for neg

Browse Source
This commit is contained in:
Henry Weller
2015-07-03 18:04:34 +01:00
parent 43030b7fbc
commit e8d45b1d3c
2 changed files with 6 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
applications/solvers/multiphase/reactingTwoPhaseEulerFoam/pU/UEqns.H
View File

@ -11,8 +11,8 @@ volScalarField Kd(fluid.Kd());
{ {
volScalarField Vm(fluid.Vm()); volScalarField Vm(fluid.Vm());
const volScalarField dmdt12(fluid.dmdt()*pos(fluid.dmdt())); const volScalarField dmdt12(posPart(fluid.dmdt()));
const volScalarField dmdt21(fluid.dmdt()*neg(fluid.dmdt())); const volScalarField dmdt21(negPart(fluid.dmdt()));
{ {
U1Eqn = U1Eqn =

8
applications/solvers/multiphase/reactingTwoPhaseEulerFoam/phaseSystems/PhaseSystems/HeatAndMassTransferPhaseSystem/HeatAndMassTransferPhaseSystem.C
View File

@ -199,8 +199,8 @@ Foam::HeatAndMassTransferPhaseSystem<BasePhaseSystem>::momentumTransfer() const
const volVectorField& U2(pair.phase2().U()); const volVectorField& U2(pair.phase2().U());
const volScalarField dmdt(this->dmdt(pair)); const volScalarField dmdt(this->dmdt(pair));
const volScalarField dmdt12(dmdt*pos(dmdt)); const volScalarField dmdt12(posPart(dmdt));
const volScalarField dmdt21(dmdt*neg(dmdt)); const volScalarField dmdt21(negPart(dmdt));
*eqns[pair.phase1().name()] += fvm::Sp(dmdt21, U1) - dmdt21*U2; *eqns[pair.phase1().name()] += fvm::Sp(dmdt21, U1) - dmdt21*U2;
*eqns[pair.phase2().name()] += dmdt12*U1 - fvm::Sp(dmdt12, U2); *eqns[pair.phase2().name()] += dmdt12*U1 - fvm::Sp(dmdt12, U2);
@ -315,8 +315,8 @@ Foam::HeatAndMassTransferPhaseSystem<BasePhaseSystem>::heatTransfer() const
const volScalarField& K2(phase2.K()); const volScalarField& K2(phase2.K());
const volScalarField dmdt(this->dmdt(pair)); const volScalarField dmdt(this->dmdt(pair));
const volScalarField dmdt12(dmdt*pos(dmdt)); const volScalarField dmdt12(posPart(dmdt));
const volScalarField dmdt21(dmdt*neg(dmdt)); const volScalarField dmdt21(negPart(dmdt));
const volScalarField& Tf(*Tf_[pair]); const volScalarField& Tf(*Tf_[pair]);
*eqns[phase1.name()] += *eqns[phase1.name()] +=