Foam::tmp<Foam::volScalarField>
Foam::kineticTheoryModels::granularPressureModels::SyamlalRogersOBrien::
granularPressureCoeff
(
const phaseModel& phase1,
const phaseModel& phase2,
const volScalarField& Theta1,
const volScalarField& Theta2,
const volScalarField& g0,
const dimensionedScalar& e
) const
{
if (&phase1 != &phase2)
{
return tmp<volScalarField>
(
new volScalarField
(
IOobject
(
"granularPressureCoeff",
phase1.time().timeName(),
phase1.mesh(),
IOobject::NO_READ,
IOobject::NO_WRITE,
false
),
phase1.mesh(),
dimensionedScalar("0", dimDensity*sqr(dimVelocity), 0.0)
)
);
}
return 2.0*phase1.rho()*Theta1*(1.0 + e)*sqr(phase1)*g0;
}
Foam::tmp<Foam::volScalarField> Foam::blendingMethods::noBlending::f2
(
const phaseModel& phase1,
const phaseModel& phase2
) const
{
const fvMesh& mesh(phase1.mesh());
return
tmp<volScalarField>
(
new volScalarField
(
IOobject
(
"f",
mesh.time().timeName(),
mesh
),
mesh,
dimensionedScalar
(
"f",
dimless,
phase1.name() == continuousPhase_
)
)
);
}
Foam::phasePair::phasePair
(
const phaseModel& phase1,
const phaseModel& phase2,
const bool ordered
)
:
phasePairKey(phase1.name(), phase2.name(), ordered),
phase1_(phase1),
phase2_(phase2),
g_(phase1.mesh().lookupObject<uniformDimensionedVectorField>("g"))
{}
Foam::tmp<Foam::volScalarField> Foam::blendingMethods::noBlending::f2
(
const phaseModel& phase1,
const phaseModel& phase2
) const
{
const fvMesh& mesh(phase1.mesh());
return volScalarField::New
(
"f",
mesh,
dimensionedScalar(dimless, phase1.name() == continuousPhase_)
);
}