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 +
fvModels.source(alpha, rho, 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
layerAndWeight max(const layerAndWeight &a, const layerAndWeight &b)
volScalarField alpha(IOobject("alpha", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), lambda *max(Ua &U, zeroSensitivity))
Ostream & endl(Ostream &os)
Add newline and flush stream.
GeometricField< vector, fvPatchField, volMesh > volVectorField
fvScalarMatrix EEqn(fvm::ddt(rho, he)+mvConvection->fvmDiv(phi, he)+fvc::ddt(rho, K)+fvc::div(phi, K)+(he.name()=="e" ? mvConvection->fvcDiv(fvc::absolute(phi, rho, U), p/rho) :-dpdt)+thermophysicalTransport->divq(he)==reaction->Qdot()+fvModels.source(rho, he))
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Foam::fvConstraints & fvConstraints
layerAndWeight min(const layerAndWeight &a, const layerAndWeight &b)
bool constrain(fvMatrix< Type > &eqn) const
Apply constraints to an equation.
Foam::fvModels & fvModels
phaseSystem::phaseModelList & phases
const dimensionedVector & g