TEqn.H
Go to the documentation of this file.
1{
3 (
4 fvm::ddt(rho, T) + fvm::div(rhoPhi, T) - fvm::Sp(contErr, T)
5 - fvm::laplacian(turbulence.alphaEff(), T)
6 + (
7 divUp()// - contErr/rho*p
8 + (fvc::ddt(rho, K) + fvc::div(rhoPhi, K))() - contErr*K
9 )
10 *(
11 alpha1()/mixture.thermo1().Cv()()
12 + alpha2()/mixture.thermo2().Cv()()
13 )
14 ==
16 );
17
18 TEqn.relax();
19
20 fvOptions.constrain(TEqn);
21
22 TEqn.solve();
23
24 fvOptions.correct(T);
25
26 mixture.correctThermo();
27 mixture.correct();
28}
CGAL::Exact_predicates_exact_constructions_kernel K
rhoPhi
Definition: rhoEqn.H:10
fv::options & fvOptions
const volScalarField & alpha1
const volScalarField & alpha2
const volScalarField & T
volScalarField::Internal contErr((fvc::ddt(rho)+fvc::div(rhoPhi) -(fvOptions(alpha1, mixture.thermo1().rho())&rho1) -(fvOptions(alpha2, mixture.thermo2().rho())&rho2))())
compressible::turbulenceModel & turbulence
fvScalarMatrix TEqn(fvm::ddt(T)+fvm::div(phi, T) - fvm::laplacian(alphaEff, T)==radiation->ST(rhoCpRef, T)+fvOptions(T))
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:45
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\n"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture
Definition: createFields.H:39