Ejemplo n.º 1
0
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());
}
Ejemplo n.º 2
0
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());
}
Ejemplo n.º 3
0
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));
}
Ejemplo n.º 4
0
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())
        );
}