setRegionSolidFields.H
Go to the documentation of this file.
1fvMesh& mesh = solidRegions[i];
2solidThermo& thermo = thermos[i];
3
4tmp<volScalarField> trho = thermo.rho();
5const volScalarField& rho = trho();
6
7tmp<volScalarField> tcp = thermo.Cp();
8const volScalarField& cp = tcp();
9
10tmp<volSymmTensorField> taniAlpha;
11if (!thermo.isotropic())
12{
13 volSymmTensorField& aniAlpha = aniAlphas[i];
14 tmp<volVectorField> tkappaByCp = thermo.Kappa()/cp;
15 const coordinateSystem& coodSys = coordinates[i];
16
17 aniAlpha.primitiveFieldRef() =
18 coodSys.transformPrincipal
19 (
20 mesh.cellCentres(),
21 tkappaByCp()
22 );
23
24 aniAlpha.correctBoundaryConditions();
25
26 taniAlpha = tmp<volSymmTensorField>
27 (
28 new volSymmTensorField(aniAlpha)
29 );
30}
31
32
33volScalarField& h = thermo.he();
34
35const volScalarField& betav = betavSolid[i];
36
37fv::options& fvOptions = solidHeatSources[i];
Basic thermodynamics type based on the use of fitting functions for cp, h, s obtained from the templa...
PtrList< coordinateSystem > coordinates(solidRegions.size())
PtrList< solidThermo > thermos(solidRegions.size())
PtrList< volSymmTensorField > aniAlphas(solidRegions.size())
PtrList< volScalarField > betavSolid(solidRegions.size())
PtrList< fv::options > solidHeatSources(solidRegions.size())
PtrList< fvMesh > solidRegions(solidNames.size())
tmp< volScalarField > tcp
fv::options & fvOptions
tmp< volScalarField > trho
volScalarField & h
fvMesh & mesh
const volScalarField & rho
tmp< volSymmTensorField > taniAlpha
const volScalarField & betav
const volScalarField & cp