tmp<surfaceInterpolationScheme<Type> > scheme
(
const surfaceScalarField& faceFlux,
const word& name
)
{
return surfaceInterpolationScheme<Type>::New
(
faceFlux.mesh(),
faceFlux,
faceFlux.mesh().interpolationScheme(name)
);
}
Foam::combustionModel::combustionModel
(
const dictionary& combustionProperties,
const hsCombustionThermo& thermo,
const compressible::turbulenceModel& turbulence,
const surfaceScalarField& phi,
const volScalarField& rho
)
:
combustionModelCoeffs_
(
combustionProperties.subDict
(
word(combustionProperties.lookup("combustionModel")) + "Coeffs"
)
),
thermo_(thermo),
turbulence_(turbulence),
mesh_(phi.mesh()),
phi_(phi),
rho_(rho),
stoicRatio_(thermo.lookup("stoichiometricAirFuelMassRatio")),
s_(thermo.lookup("stoichiometricOxygenFuelMassRatio")),
qFuel_(thermo_.lookup("qFuel")),
composition_(thermo.composition())
{}
tmp<GeometricField<Type, fvPatchField, volMesh> >
DDt
(
const surfaceScalarField& phi,
const GeometricField<Type, fvPatchField, volMesh>& psi
)
{
tmp<GeometricField<Type, fvPatchField, volMesh> > ddtDivPhiPsi
= fvc::ddt(psi) + fvc::div(phi, psi);
if (phi.mesh().moving())
{
return ddtDivPhiPsi - fvc::div(phi + phi.mesh().phi())*psi;
}
else
{
return ddtDivPhiPsi - fvc::div(phi)*psi;
}
}
void Foam::fvc::makeRelative
(
surfaceScalarField& phi,
const volVectorField& U
)
{
if (phi.mesh().moving())
{
phi -= fvc::meshPhi(U);
}
}
void Foam::fvc::makeAbsolute
(
surfaceScalarField& phi,
const volVectorField& U
)
{
if (phi.mesh().moving())
{
phi += fvc::meshPhi(U);
}
}
void Foam::fvc::makeRelative
(
surfaceScalarField& phi,
const dimensionedScalar& rho,
const volVectorField& U
)
{
if (phi.mesh().moving())
{
phi -= rho*fvc::meshPhi(rho, U);
}
}
//JK interpolate(rho) is declared but not defined
void Foam::fvc::makeAbsolute
(
surfaceScalarField& phi,
const volScalarField& rho,
const volVectorField& U
)
{
if (phi.mesh().moving())
{
phi += fvc::interpolate(rho)*fvc::meshPhi(rho, U);
}
}
tmp<surfaceInterpolationScheme<Type> > scheme
(
const surfaceScalarField& faceFlux,
Istream& streamData
)
{
return surfaceInterpolationScheme<Type>::New
(
faceFlux.mesh(),
faceFlux,
streamData
);
}
Foam::combustionModel::combustionModel
(
const dictionary& combustionProps,
hsCombustionThermo& thermo,
const compressible::turbulenceModel& turbulence,
const surfaceScalarField& phi,
const volScalarField& rho
)
:
coeffs_(dictionary::null),
thermo_(thermo),
turbulence_(turbulence),
mesh_(phi.mesh()),
phi_(phi),
rho_(rho)
{}
