solveSolid.H
Go to the documentation of this file.
1{
2 for (int nonOrth=0; nonOrth<=nNonOrthCorr; ++nonOrth)
3 {
5 (
6 (
7 thermo.isotropic()
8 ? -fvm::laplacian(betav*thermo.alpha(), h, "laplacian(alpha,h)")
9 : -fvm::laplacian(betav*taniAlpha(), h, "laplacian(alpha,h)")
10 )
11 ==
13 );
14
15 hEqn.relax();
16
17 fvOptions.constrain(hEqn);
18
19 if (coupled)
20 {
21 fvMatrixAssemblyPtr->addFvMatrix(hEqn);
22 }
23 else
24 {
25 Info<< "\nSolving for solid region "
26 << solidRegions[i].name() << endl;
27
28 hEqn.solve();
29 fvOptions.correct(h);
30 thermo.correct();
31
32 Info<< "Min/max T:" << min(thermo.T()).value() << ' '
33 << max(thermo.T()).value() << endl;
34 }
35 }
36}
Y[inertIndex] max(0.0)
fv::options & fvOptions
Basic thermodynamics type based on the use of fitting functions for cp, h, s obtained from the templa...
autoPtr< fvMatrix< scalar > > fvMatrixAssemblyPtr
bool coupled(solutionDict.getOrDefault("coupledEnergyField", false))
PtrList< fvMesh > solidRegions(solidNames.size())
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:45
messageStream Info
Information stream (stdout output on master, null elsewhere)
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:372
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
Definition: hashSets.C:33
const int nNonOrthCorr
volScalarField & h
tmp< volSymmTensorField > taniAlpha
const volScalarField & betav