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

reactingEulerFoam: improvements to population balance modeling

Browse Source 
allow renormalization of sizeGroup volume fractions for restarts involving
initial conditions with a slight degree of unboundedness

Patch contributed by Institute of Fluid Dynamics, Helmholtz-Zentrum Dresden - Rossendorf
(HZDR) and VTT Technical Research Centre of Finland Ltd.
This commit is contained in:
Henry Weller
2018-01-04 17:20:28 +00:00
parent 9b58f5226b
commit 79207b6147
1 changed files with 29 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files

41
applications/solvers/multiphase/reactingEulerFoam/phaseSystems/diameterModels/velocityGroup/velocityGroup.C
View File

@ -145,12 +145,6 @@ Foam::diameterModels::velocityGroup::fSum() const
sumSizeGroups += sizeGroups_[i];
}
Info<< phase_.name() << " sizeGroups-sum volume fraction, min, max = "
<< sumSizeGroups.weightedAverage(phase_.mesh().V()).value()
<< ' ' << min(sumSizeGroups).value()
<< ' ' << max(sumSizeGroups).value()
<< endl;
return tsumSizeGroups;
}
@ -225,7 +219,7 @@ Foam::diameterModels::velocityGroup::velocityGroup
),
phase.mesh()
),
formFactor_("formFactor", dimless, diameterProperties.lookup("formFactor")),
formFactor_("formFactor", dimless, diameterProperties),
sizeGroups_
(
diameterProperties.lookup("sizeGroups"),
@ -314,12 +308,21 @@ Foam::diameterModels::velocityGroup::velocityGroup
{
if
(
mag
phase_.mesh().solverDict(popBalName_).lookupOrDefault<Switch>
(
1.0
- fSum_.weightedAverage(fSum_.mesh().V()).value()
"renormalize",
false
)
>= 1e-5
)
{
renormalize();
}
fSum_ = fSum();
if
(
mag(1 - fSum_.weightedAverage(fSum_.mesh().V()).value()) >= 1e-5
|| mag(1 - max(fSum_).value()) >= 1e-5
|| mag(1 - min(fSum_).value()) >= 1e-5
)
@ -327,7 +330,9 @@ Foam::diameterModels::velocityGroup::velocityGroup
FatalErrorInFunction
<< "Initial values of the sizeGroups belonging to velocityGroup "
<< this->phase().name()
<< " must add to unity." << endl
<< " must add to unity. The sizeGroup fractions can be"
<< " renormalized by setting the renormalize switch"
<< " in the fvSolution subdictionary " << popBalName_ << "." << endl
<< exit(FatalError);
}
@ -367,8 +372,20 @@ void Foam::diameterModels::velocityGroup::postSolve()
d_ = dsm();
Info<< this->phase().name() << " Sauter mean diameter, min, max = "
<< d_.weightedAverage(d_.mesh().V()).value()
<< ' ' << min(d_).value()
<< ' ' << max(d_).value()
<< endl;
fSum_ = fSum();
Info<< phase_.name() << " sizeGroups-sum volume fraction, min, max = "
<< fSum_.weightedAverage(phase_.mesh().V()).value()
<< ' ' << min(fSum_).value()
<< ' ' << max(fSum_).value()
<< endl;
if
(
phase_.mesh().solverDict(popBalName_).lookupOrDefault<Switch>