solidRegionDiffusionNo.H
Go to the documentation of this file.
1 scalar DiNum = -GREAT;
2 
4 {
5  // Note: do not use setRegionSolidFields.H to avoid double registering Cp
6  //#include "setRegionSolidFields.H"
7  const solidThermo& thermo = thermos[i];
8 
9  tmp<volScalarField> magKappa;
10  if (thermo.isotropic())
11  {
12  magKappa = thermo.kappa();
13  }
14  else
15  {
16  magKappa = mag(thermo.Kappa());
17  }
18 
19  tmp<volScalarField> tcp = thermo.Cp();
20  const volScalarField& cp = tcp();
21 
22  tmp<volScalarField> trho = thermo.rho();
23  const volScalarField& rho = trho();
24 
25  DiNum = max
26  (
28  (
29  solidRegions[i],
30  runTime,
31  rho*cp,
32  magKappa()
33  ),
34  DiNum
35  );
36 
37 }
scalar DiNum
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
forAll(solidRegions, i)
tmp< volScalarField > trho
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:52
psiReactionThermo & thermo
Definition: createFields.H:31
scalar solidRegionDiffNo(const fvMesh &mesh, const Time &runTime, const volScalarField &Cprho, const volScalarField &kappa)
PtrList< solidThermo > thermos(solidRegions.size())
const volScalarField & cp
tmp< volScalarField > tcp
PtrList< fvMesh > solidRegions(solidsNames.size())
dimensioned< scalar > mag(const dimensioned< Type > &)