solveSolid.H
Go to the documentation of this file.
1 {
2  for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
3  {
4  fvScalarMatrix hEqn
5  (
6  (
7  thermo.isotropic()
8  ? -fvm::laplacian(betav*thermo.alpha(), h, "laplacian(alpha,h)")
9  : -fvm::laplacian(betav*taniAlpha(), h, "laplacian(alpha,h)")
10  )
11  ==
12  fvOptions(rho, h)
13  );
14 
15  hEqn.relax();
16 
17  fvOptions.constrain(hEqn);
18 
19  hEqn.solve();
20 
21  fvOptions.correct(h);
22  }
23 }
24 
25 thermo.correct();
26 
27 Info<< "Min/max T:" << min(thermo.T()).value() << ' '
28  << max(thermo.T()).value() << endl;
const volScalarField & betav
tmp< GeometricField< Type, fvPatchField, volMesh > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
Definition: fvcLaplacian.C:45
messageStream Info
tmp< volSymmTensorField > taniAlpha
const int nNonOrthCorr
fv::IOoptionList & fvOptions
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:42
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
psiReactionThermo & thermo
Definition: createFields.H:32
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
const dimensionedScalar h
Planck constant.
Definition: createFields.H:6