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&
9 
10  {
11  tmp<fvScalarMatrix> E1Eqn(phase1.heEqn());
12 
13  if (E1Eqn.valid())
14  {
15  E1Eqn =
16  (
17  E1Eqn
18  ==
19  *heatTransfer[phase1.name()]
20  + alpha1*rho1*(U1&g)
21  + fvOptions(alpha1, rho1, phase1.thermo().he())
22  );
23 
24  E1Eqn->relax();
25  fvOptions.constrain(E1Eqn());
26  E1Eqn->solve();
27  }
28  }
29 
30  {
31  tmp<fvScalarMatrix> E2eqn(phase2.heEqn());
32 
33  if (E2eqn.valid())
34  {
35  E2eqn =
36  (
37  E2eqn
38  ==
39  *heatTransfer[phase2.name()]
40  + alpha2*rho2*(U2&g)
41  + fvOptions(alpha2, rho2, phase2.thermo().he())
42  );
43 
44  E2eqn->relax();
45  fvOptions.constrain(E2eqn());
46  E2eqn->solve();
47  }
48  }
49 }
50 
51 fluid.correctThermo();
52 
53 Info<< phase1.name() << " min/max T "
54  << min(phase1.thermo().T()).value()
55  << " - "
56  << max(phase1.thermo().T()).value()
57  << endl;
58 
59 Info<< phase2.name() << " min/max T "
60  << min(phase2.thermo().T()).value()
61  << " - "
62  << max(phase2.thermo().T()).value()
63  << endl;
phaseModel & phase1
Definition: createFields.H:12
phaseModel & phase2
Definition: createFields.H:13
volVectorField & U1
Definition: createFields.H:18
multiphaseSystem & fluid
Definition: createFields.H:10
volVectorField & U2
Definition: createFields.H:23
messageStream Info
fv::IOoptionList & fvOptions
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
rho1
Definition: pEqn.H:124
phaseSystem::heatTransferTable & heatTransfer
Definition: EEqns.H:7
rho2
Definition: pEqn.H:125
volScalarField & alpha1
Definition: createFields.H:15
alpha2
Definition: alphaEqn.H:112
const dimensionedVector & g
E1Eqn
Definition: EEqns.H:33
autoPtr< phaseSystem::heatTransferTable > heatTransferPtr(fluid.heatTransfer())
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)