solveSolid.H
Go to the documentation of this file.
1{
3 (
4 fvm::ddt(betav*rho, h)
5 - (
6 thermo.isotropic()
7 ? fvm::laplacian(betav*thermo.alpha(), h, "laplacian(alpha,h)")
8 : fvm::laplacian(betav*taniAlpha(), h, "laplacian(alpha,h)")
9 )
10 ==
12 );
13
14 hEqn.relax();
15
16 fvOptions.constrain(hEqn);
17
19 {
20 fvMatrixAssemblyPtr->addFvMatrix(hEqn);
21 }
22 else
23 {
24 Info<< "\nSolving for solid region "<< solidRegions[i].name() << endl;
25
26 if (finalIter)
27 {
28 mesh.data::add("finalIteration", true);
29 }
30
31 hEqn.solve(mesh.solver(h.select(finalIter)));
32
33 fvOptions.correct(h);
34
35 thermo.correct();
36
37 Info<< "Min/max T:" << min(thermo.T()).value() << ' '
38 << max(thermo.T()).value() << endl;
39
40 if (finalIter)
41 {
42 mesh.data::remove("finalIteration");
43 }
44 }
45}
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))
dynamicFvMesh & mesh
PtrList< fvMesh > solidRegions(solidNames.size())
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:45
volScalarField & h
tmp< volSymmTensorField > taniAlpha
const volScalarField & betav