Foam::tmp<Foam::volScalarField>
Foam::kineticTheoryModels::radialModels::Gao::g0
(
const phaseModel& phase1,
const phaseModel& phase2
) const
{
volScalarField g1(1.0/(1.0 - cbrt(min(phase1/phase1.alphaMax(), 0.999))));
volScalarField g2(1.0/(1.0 - cbrt(min(phase2/phase2.alphaMax(), 0.999))));
return (phase1.d()*g1 + phase2.d()*g2)/(phase1.d() + phase2.d());
}
Foam::tmp<Foam::volScalarField>
Foam::kineticTheoryModels::radialModels::Gao::g0prime
(
const phaseModel& phase1,
const phaseModel& phase2
) const
{
volScalarField aByaMax
(
cbrt(min(max(phase1, scalar(1e-3))/phase1.alphaMax(), 0.999))
);
return
(1.0/(3.0*phase1.alphaMax()))/sqr(aByaMax - sqr(aByaMax))
*phase1.d()/(phase1.d() + phase2.d());
}
Foam::tmp<Foam::volScalarField>
Foam::kineticTheoryModels::radialModels::SinclairJackson::g0prime
(
const phaseModel& phase1,
const phaseModel& phase2
) const
{
volScalarField aByaMax
(
cbrt
(
min
(
max(phase1, scalar(1e-3)),
alphaMinFriction_
)/phase1.alphaMax()
)
);
return (1.0/(3*phase1.alphaMax()))/sqr(aByaMax - sqr(aByaMax));
}
Foam::tmp<Foam::volScalarField>
Foam::kineticTheoryModels::radialModels::SinclairJackson::g0
(
const phaseModel& phase1,
const phaseModel& phase2
) const
{
return
1.0
/(
1.0
- cbrt(min(phase1, alphaMinFriction_)/phase1.alphaMax())
);
}