zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

reactingMultiphaseEulerFoam: Updated using new pair "loop" structure

Browse Source
This commit is contained in:
Henry Weller
2018-01-01 16:58:42 +00:00
parent a9cb40b55b
commit 87c507d8c8
2 changed files with 18 additions and 40 deletions
Download Patch File Download Diff File Expand all files Collapse all files

22
applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/pUf/UEqns.H
View File

@ -60,22 +60,16 @@ PtrList<fvVectorMatrix> UEqns(phases.size());
const phasePair& pair(fluid.phasePairs()[VmIter.key()]);
const phaseModel* phase1 = &pair.phase1();
const phaseModel* phase2 = &pair.phase2();
forAllConstIter(phasePair, pair, iter)
if (pair.contains(phase))
{
if (phase1 == &phase)
{
UEqns[phasei] +=
Vm
*(
UgradUs[phasei]
- (UgradUs[phase2->index()] & phase2->U())
);
}
const phaseModel& otherPhase = pair.other(phase);
Swap(phase1, phase2);
UEqns[phasei] +=
Vm
*(
UgradUs[phasei]
- (UgradUs[otherPhase.index()] & otherPhase.U())
);
}
}

36
applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/pUf/pEqn.H
View File

@ -194,18 +194,10 @@ while (pimple.correct())
const phasePair& pair(fluid.phasePairs()[KdfIter.key()]);
const phaseModel* phase1 = &pair.phase1();
const phaseModel* phase2 = &pair.phase2();
forAllConstIter(phasePair, pair, iter)
if (pair.contains(phase))
{
if (phase1 == &phase)
{
phiHbyAs[phasei] +=
rAUfs[phasei]*Kf*MRF.absolute(phase2->phi());
}
Swap(phase1, phase2);
phiHbyAs[phasei] +=
rAUfs[phasei]*Kf*MRF.absolute(pair.other(phase).phi());
}
}
@ -221,22 +213,14 @@ while (pimple.correct())
const phasePair& pair(fluid.phasePairs()[VmfIter.key()]);
const phaseModel* phase1 = &pair.phase1();
const phaseModel* phase2 = &pair.phase2();
forAllConstIter(phasePair, pair, iter)
if (pair.contains(phase))
{
if (phase1 == &phase)
{
phiHbyAs[phasei] +=
rAUfs[phasei]
*(
Vmf*byDt(MRF.absolute(phase1->phi()().oldTime()))
+ Vmf*ddtPhis[phase2->index()]
);
}
Swap(phase1, phase2);
phiHbyAs[phasei] +=
rAUfs[phasei]
*(
Vmf*byDt(MRF.absolute(phase.phi().oldTime()))
+ Vmf*ddtPhis[pair.other(phase).index()]
);
}
}