EEqns.H
Go to the documentation of this file.
1for (int Ecorr=0; Ecorr<nEnergyCorrectors; Ecorr++)
2{
3 fluid.correctEnergyTransport();
4
5 autoPtr<phaseSystem::heatTransferTable>
6 heatTransferPtr(fluid.heatTransfer());
7
8 phaseSystem::heatTransferTable& heatTransfer = heatTransferPtr();
9
10 forAll(fluid.anisothermalPhases(), anisothermalPhasei)
11 {
12 phaseModel& phase = fluid.anisothermalPhases()[anisothermalPhasei];
13
14 const volScalarField& alpha = phase;
15 const volScalarField& rho = phase.rho();
16 const volVectorField& U = phase.U();
17
18 fvScalarMatrix EEqn
19 (
20 phase.heEqn()
21 ==
22 *heatTransfer[phase.name()]
23 + alpha*rho*(U&g)
24 + fvOptions(alpha, rho, phase.thermoRef().he())
25 );
26
27 EEqn.relax();
28 fvOptions.constrain(EEqn);
29 EEqn.solve();
30 fvOptions.correct(phase.thermoRef().he());
31 }
32
33 fluid.correctThermo();
34 fluid.correct();
35}
36
37
39{
40 phaseModel& phase = phases[phasei];
41
42 Info<< phase.name() << " min/max T "
43 << min(phase.thermo().T()).value()
44 << " - "
45 << max(phase.thermo().T()).value()
46 << endl;
47}
Y[inertIndex] max(0.0)
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))
const uniformDimensionedVectorField & g
fv::options & fvOptions
int nEnergyCorrectors(pimpleDict.getOrDefault< int >("nEnergyCorrectors", 1))
twoPhaseSystem & fluid
U
Definition: pEqn.H:72
label phasei
Definition: pEqn.H:27
volScalarField & alpha
multiphaseSystem::phaseModelList & phases
Definition: createFields.H:12
#define forAll(list, i)
Loop across all elements in list.
Definition: stdFoam.H:333