47 lines
1 KiB
C++
47 lines
1 KiB
C++
{
|
|
fvScalarMatrix eEqn
|
|
(
|
|
fvm::ddt(rho, e)
|
|
+ fvm::div(phi, e)
|
|
- fvm::laplacian(turbulence->alphaEff(), e)
|
|
==
|
|
- fvm::SuSp
|
|
(
|
|
p*fvc::div
|
|
(
|
|
phi/fvc::interpolate(rho)
|
|
+ fvc::meshPhi(rho, U)
|
|
)/e,
|
|
e
|
|
)
|
|
// viscous heating?
|
|
);
|
|
|
|
eEqn.relax
|
|
(
|
|
mesh.solutionDict().equationRelaxationFactor(e.select(pimple.finalIter()))
|
|
);
|
|
|
|
eEqn.solve();
|
|
|
|
// Bound the energy using TMin and TMax
|
|
{
|
|
dimensionedScalar Tstd("Tstd", dimTemperature, specie::Tstd());
|
|
|
|
volScalarField Cv = thermo.Cv();
|
|
volScalarField R = thermo.Cp() - Cv;
|
|
|
|
e = Foam::min(e, TMax*Cv + R*Tstd);
|
|
e = Foam::max(e, TMin*Cv + R*Tstd);
|
|
e.correctBoundaryConditions();
|
|
}
|
|
|
|
thermo.correct();
|
|
|
|
// Recalculate compressibility
|
|
psis = thermo.psi()/thermo.Cp()*thermo.Cv();
|
|
|
|
// Recalculate density
|
|
rho = thermo.rho();
|
|
rho.correctBoundaryConditions();
|
|
}
|