1 for (
int Ecorr=0; Ecorr<nEnergyCorrectors; Ecorr++)
3 fluid.correctEnergyTransport();
5 autoPtr<phaseSystem::heatTransferTable>
6 heatTransferPtr(
fluid.heatTransfer());
8 phaseSystem::heatTransferTable& heatTransfer = heatTransferPtr();
10 forAll(
fluid.anisothermalPhases(), anisothermalPhasei)
12 phaseModel& phase =
fluid.anisothermalPhases()[anisothermalPhasei];
15 tmp<volScalarField> tRho = phase.rho();
17 tmp<volVectorField> tU = phase.U();
24 *heatTransfer[phase.name()]
26 +
fvOptions(alpha, rho, phase.thermoRef().he())
32 fvOptions.correct(phase.thermoRef().he());
35 fluid.correctThermo();
36 fluid.correctContinuityError();
44 Info<< phase.name() <<
" min/max T " 45 <<
min(phase.thermo().T()).value()
47 <<
max(phase.thermo().T()).value()
fvMatrix< scalar > fvScalarMatrix
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
Ostream & endl(Ostream &os)
Add newline and flush stream.
GeometricField< vector, fvPatchField, volMesh > volVectorField
GeometricField< scalar, fvPatchField, volMesh > volScalarField
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
phaseSystem::phaseModelList & phases
volScalarField alpha(IOobject("alpha", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), lambda *max(Ua &U, zeroSensitivity))
fvScalarMatrix EEqn(fvm::ddt(rho, he)+mvConvection->fvmDiv(phi, he)+fvc::ddt(rho, K)+fvc::div(phi, K)+(he.name()=="e" ? fvc::div(fvc::absolute(phi/fvc::interpolate(rho), U), p, "div(phiv,p)") :-dpdt)+thermophysicalTransport->divq(he)==combustion->Qdot()+radiation->Sh(thermo, he)+parcels.Sh(he)+surfaceFilm.Sh()+fvOptions(rho, he))
const dimensionedVector & g