solidDiffusionNo.H
Go to the documentation of this file.
1 scalar DiNum = -GREAT;
2 
3 {
4  tmp<volScalarField> magKappa;
5  if (thermo.isotropic())
6  {
7  magKappa = thermo.kappa();
8  }
9  else
10  {
11  magKappa = mag(thermo.Kappa());
12  }
13 
14  tmp<volScalarField> tcp = thermo.Cp();
15  const volScalarField& cp = tcp();
16 
17  tmp<volScalarField> trho = thermo.rho();
18  const volScalarField& rho = trho();
19 
20  DiNum = max
21  (
22  solidRegionDiffNo
23  (
24  mesh,
25  runTime,
26  rho*cp,
27  magKappa()
28  ),
29  DiNum
30  );
31 }
Foam::solidRegionDiffNo
scalar solidRegionDiffNo(const fvMesh &mesh, const Time &runTime, const volScalarField &Cprho, const volScalarField &kappa)
Foam::mag
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
Foam::max
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
Definition: hashSets.C:40
runTime
engineTime & runTime
Definition: createEngineTime.H:13
cp
const volScalarField & cp
Definition: solidDiffusionNo.H:15
Foam::volScalarField
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:72
thermo
psiReactionThermo & thermo
Definition: createFields.H:28
mesh
dynamicFvMesh & mesh
Definition: createDynamicFvMesh.H:6
DiNum
scalar DiNum
Definition: solidDiffusionNo.H:1
rho
const volScalarField & rho
Definition: solidDiffusionNo.H:18
trho
tmp< volScalarField > trho
Definition: solidDiffusionNo.H:17
tcp
tmp< volScalarField > tcp
Definition: solidDiffusionNo.H:14