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];
22 *heatTransfer[phase.name()]
24 +
fvOptions(alpha, rho, phase.thermoRef().he())
30 fvOptions.correct(phase.thermoRef().he());
33 fluid.correctThermo();
42 Info<< phase.name() <<
" min/max T " 43 <<
min(phase.thermo().T()).value()
45 <<
max(phase.thermo().T()).value()
fvMatrix< scalar > fvScalarMatrix
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
multiphaseSystem::phaseModelList & phases
Ostream & endl(Ostream &os)
Add newline and flush stream.
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) - fvm::laplacian(turbulence->alphaEff(), he)==Qdot+radiation->Sh(thermo, he)+parcels.Sh(he)+surfaceFilm.Sh()+fvOptions(rho, he))
GeometricField< vector, fvPatchField, volMesh > volVectorField
GeometricField< scalar, fvPatchField, volMesh > volScalarField
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
const dimensionedVector & g
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].