EEqns.H
Go to the documentation of this file.
1 for (int Ecorr=0; Ecorr<nEnergyCorrectors; Ecorr++)
2 {
3  fluid.correctEnergyTransport();
4 
5  autoPtr<phaseSystem::heatTransferTable>
6  heatTransferPtr(fluid.heatTransfer());
7 
8  phaseSystem::heatTransferTable&
9  heatTransfer = heatTransferPtr();
10 
11  if (!phase1.isothermal())
12  {
14  (
15  phase1.heEqn()
16  ==
17  *heatTransfer[phase1.name()]
18  + alpha1*rho1*(U1&g)
19  + fvOptions(alpha1, rho1, thermo1.he())
20  );
21 
22  E1Eqn.relax();
23  fvOptions.constrain(E1Eqn);
24  E1Eqn.solve();
25  fvOptions.correct(thermo1.he());
26  }
27 
28  if (!phase2.isothermal())
29  {
31  (
32  phase2.heEqn()
33  ==
34  *heatTransfer[phase2.name()]
35  + alpha2*rho2*(U2&g)
36  + fvOptions(alpha2, rho2, phase2.thermoRef().he())
37  );
38 
39  E2Eqn.relax();
40  fvOptions.constrain(E2Eqn);
41  E2Eqn.solve();
42  fvOptions.correct(thermo2.he());
43  }
44 
45  fluid.correctThermo();
46  fluid.correct();
47 }
rhoThermo & thermo2
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:42
fv::options & fvOptions
multiphaseSystem & fluid
Definition: createFields.H:11
E1Eqn
Definition: EEqns.H:34
E2Eqn
Definition: EEqns.H:67
rhoThermo & thermo1
const volScalarField & alpha1
volVectorField & U1
alpha2
Definition: alphaEqn.H:115
phaseModel & phase1
volVectorField & U2
const dimensionedScalar & rho2
Definition: createFields.H:40
phaseModel & phase2
const dimensionedVector & g
const dimensionedScalar & rho1
Definition: createFields.H:39