zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-12-28 03:37:59 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

multiphaseEulerFoam: Corrected the handling of drag for more than 2 phases

Browse Source
This commit is contained in:
Henry
2012-03-14 18:19:18 +00:00
parent f3fac5bbff
commit 0d0dfbafe3
1 changed files with 12 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

14
applications/solvers/multiphase/multiphaseEulerFoam/pEqn.H
View File

@ -76,6 +76,7 @@
+ fvc::ddtPhiCorr(rAUs[phasei], alpha, phase.U(), phase.phi())
+ rAlphaAUfs[phasei]*(g & mesh.Sf())
);
mrfZones.relativeFlux(phiHbyAs[phasei]);
multiphaseSystem::dragModelTable::const_iterator dmIter =
@ -89,11 +90,20 @@
++dmIter, ++dcIter
)
{
const phaseModel *phase2Ptr = &dmIter()->phase2();
if (phase2Ptr == &phase)
const phaseModel *phase2Ptr = NULL;
if (&phase == &dmIter()->phase1())
{
phase2Ptr = &dmIter()->phase2();
}
else if (&phase == &dmIter()->phase2())
{
phase2Ptr = &dmIter()->phase1();
}
else
{
continue;
}
phiHbyAs[phasei] +=
fvc::interpolate