TEqn.H
Go to the documentation of this file.
1{
2 radiation->correct();
3 rhoCp = rho*fluid.Cp();
4
5 const surfaceScalarField rhoCpPhi(fvc::interpolate(fluid.Cp())*rhoPhi);
6
7 const volScalarField kappaEff
8 (
9 "kappaEff",
10 fluid.kappa() + fluid.Cp()*turbulence->mut()/fluid.Prt()
11 );
12
13 fvScalarMatrix TEqn
14 (
15 fvm::ddt(rhoCp, T)
16 + fvm::div(rhoCpPhi, T, "div(phi,T)")
17 - fvm::Sp(fvc::ddt(rhoCp) + fvc::div(rhoCpPhi), T)
18 - fvm::laplacian(kappaEff, T, "laplacian(kappa,T)")
19 ==
20 fluid.heatTransfer(T)
21 + radiation->ST(T)
22 + fvOptions(rhoCp, T)
23 );
24
25 TEqn.relax();
26
27 fvOptions.constrain(TEqn);
28
29 TEqn.solve();
30
31 fvOptions.correct(T);
32
33 fluid.correct();
34
35 Info<< "min/max(T) = "
36 << min(T).value() << ", " << max(T).value() << endl;
37}
Y[inertIndex] max(0.0)
rhoPhi
Definition: rhoEqn.H:10
fv::options & fvOptions
twoPhaseSystem & fluid
const volScalarField & T
autoPtr< radiation::radiationModel > radiation(radiation::radiationModel::New(T))
compressible::turbulenceModel & turbulence
fvScalarMatrix TEqn(fvm::ddt(T)+fvm::div(phi, T) - fvm::laplacian(alphaEff, T)==radiation->ST(rhoCpRef, T)+fvOptions(T))
const surfaceScalarField rhoCpPhi(fvc::interpolate(fluid.Cp()) *rhoPhi)
rhoCp
Definition: TEqn.H:3
kappaEff
Definition: TEqn.H:10