alphaSuSp.H

Go to the documentation of this file.

tmp<volScalarField::Internal> Sp
(
    volScalarField::Internal::New
    (
        "Sp",
        mesh,
        dimensionedScalar(dgdt.dimensions(), 0)
    )
);

tmp<volScalarField::Internal> Su
(
    volScalarField::Internal::New
    (
        "Su",
        mesh,
        dimensionedScalar(dgdt.dimensions(), 0)
    )
);

// Phase change alpha1 source
Pair<tmp<volScalarField::Internal>> phaseChangeS
(
    phaseChange.Salpha(alpha1)
);

if (phaseChangeS[0].valid())
{
    Su = phaseChangeS[0];
    Sp = phaseChangeS[1];
}

volScalarField::Internal& SpRef = Sp.ref();
volScalarField::Internal& SuRef = Su.ref();

forAll(dgdt, celli)
{
    if (dgdt[celli] > 0.0)
    {
        SpRef[celli] -= dgdt[celli]/max(1.0 - alpha1[celli], 1e-4);
        SuRef[celli] += dgdt[celli]/max(1.0 - alpha1[celli], 1e-4);
    }
    else if (dgdt[celli] < 0.0)
    {
        SpRef[celli] += dgdt[celli]/max(alpha1[celli], 1e-4);
    }
}

tmp<volScalarField> divU
(
    mesh.moving()
  ? fvc::div(phiCN() + mesh.phi())
  : fvc::div(phiCN())
);