setRegionSolidFields.H
Go to the documentation of this file.
1 fvMesh& mesh = solidRegions[i];
2 solidThermo& thermo = thermos[i];
3 
4 tmp<volScalarField> trho = thermo.rho();
5 const volScalarField& rho = trho();
6 
7 tmp<volScalarField> tcp = thermo.Cp();
8 const volScalarField& cp = tcp();
9 
10 tmp<volSymmTensorField> taniAlpha;
11 if (!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 
34 
36 
37 fv::options& fvOptions = solidHeatSources[i];
fvOptions
fv::options & fvOptions
Definition: setRegionSolidFields.H:37
cp
const volScalarField & cp
Definition: setRegionSolidFields.H:8
tcp
tmp< volScalarField > tcp
Definition: setRegionSolidFields.H:7
thermo
Basic thermodynamics type based on the use of fitting functions for cp, h, s obtained from the templa...
solidHeatSources
PtrList< fv::options > solidHeatSources(solidRegions.size())
h
volScalarField & h
Planck constant.
Definition: setRegionSolidFields.H:33
trho
tmp< volScalarField > trho
Definition: setRegionSolidFields.H:4
Foam::volSymmTensorField
GeometricField< symmTensor, fvPatchField, volMesh > volSymmTensorField
Definition: volFieldsFwd.H:65
coordinates
PtrList< coordinateSystem > coordinates(solidRegions.size())
Foam::volScalarField
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:57
betav
const volScalarField & betav
Definition: setRegionSolidFields.H:35
aniAlphas
PtrList< volSymmTensorField > aniAlphas(solidRegions.size())
solidRegions
PtrList< fvMesh > solidRegions(solidNames.size())
mesh
fvMesh & mesh
Definition: setRegionSolidFields.H:1
rho
const volScalarField & rho
Definition: setRegionSolidFields.H:5
betavSolid
PtrList< volScalarField > betavSolid(solidRegions.size())
taniAlpha
tmp< volSymmTensorField > taniAlpha
Definition: setRegionSolidFields.H:10
thermos
PtrList< solidThermo > thermos(solidRegions.size())