EEqns.H
Go to the documentation of this file.
1 {
2  fluid.correctEnergyTransport();
3 
4  autoPtr<phaseSystem::heatTransferTable>
5  heatTransferPtr(fluid.heatTransfer());
6 
7  phaseSystem::heatTransferTable& heatTransfer = heatTransferPtr();
8 
10  {
11  phaseModel& phase = phases[phasei];
12 
13  const volScalarField& alpha = phase;
14  const volScalarField& rho = phase.rho();
15  const volVectorField& U = phase.U();
16 
17  tmp<fvScalarMatrix> EEqn(phase.heEqn());
18 
19  if (EEqn.valid())
20  {
21  EEqn =
22  (
23  EEqn
24  ==
25  *heatTransfer[phase.name()]
26  + alpha*rho*(U&g)
27  + fvOptions(alpha, rho, phase.thermo().he())
28  );
29 
30  EEqn->relax();
31  fvOptions.constrain(EEqn());
32  EEqn->solve();
33  }
34  }
35 }
36 
37 fluid.correctThermo();
38 
40 {
41  phaseModel& phase = phases[phasei];
42 
43  Info<< phase.name() << " min/max T "
44  << min(phase.thermo().T()).value()
45  << " - "
46  << max(phase.thermo().T()).value()
47  << endl;
48 }
forAll(phases, phasei)
Definition: EEqns.H:9
multiphaseSystem & fluid
Definition: createFields.H:10
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)==combustion->Sh()+radiation->Sh(thermo)+parcels.Sh(he)+surfaceFilm.Sh()+fvOptions(rho, he))
label phasei
Definition: pEqn.H:37
messageStream Info
fv::IOoptionList & fvOptions
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:52
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
phaseSystem::heatTransferTable & heatTransfer
Definition: EEqns.H:7
GeometricField< vector, fvPatchField, volMesh > volVectorField
Definition: volFieldsFwd.H:55
const dimensionedVector & g
autoPtr< phaseSystem::heatTransferTable > heatTransferPtr(fluid.heatTransfer())
multiphaseSystem::phaseModelList & phases
Definition: createFields.H:11
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
U
Definition: pEqn.H:82