TEqn.H
Go to the documentation of this file.
1{
3 (
4 fvm::ddt(rho, T) + fvm::div(rhoPhi, T)
5 - fvm::Sp(contErr, T)
6 - fvm::laplacian(turbulence.alphaEff(), T)
7 + (
8 (
9 fvc::div(fvc::absolute(phi, U), p)
10 + fvc::ddt(rho, K) + fvc::div(rhoPhi, K)
11 )()() - contErr*K
12 )
13 *(
14 alpha1/mixture.thermo1().Cv()
15 + alpha2/mixture.thermo2().Cv()
16 )()()
17 ==
19 + pos(Srho)
20 *(
21 surfaceFilm.Sh()()/mixture.thermo1().Cp()()
22 + surfaceFilm.Tref*Srho
23 )
24 );
25
26 TEqn.relax();
27
28 fvOptions.constrain(TEqn);
29
30 TEqn.solve();
31
32 fvOptions.correct(T);
33
34 mixture.correctThermo();
35 mixture.correct();
36}
CGAL::Exact_predicates_exact_constructions_kernel K
rhoPhi
Definition: rhoEqn.H:10
fv::options & fvOptions
surfaceScalarField & phi
const volScalarField & alpha1
const volScalarField & alpha2
U
Definition: pEqn.H:72
volScalarField & p
const volScalarField & T
regionModels::surfaceFilmModel & surfaceFilm
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))
dimensionedScalar pos(const dimensionedScalar &ds)
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:45
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\n"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture
Definition: createFields.H:39