All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
solveSolid.H
Go to the documentation of this file.
1 {
2  while (pimple.correctNonOrthogonal())
3  {
4  fvScalarMatrix eEqn
5  (
6  fvm::ddt(rho, e)
7  + thermo.divq(e)
8  ==
9  fvModels.source(rho, e)
10  );
11 
12  eEqn.relax();
13 
15 
16  eEqn.solve();
17 
19  }
20 }
21 
22 thermo.correct();
23 
24 Info<< "Min/max T:" << min(thermo.T()).value() << ' '
25  << max(thermo.T()).value() << endl;
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:42
fluidReactionThermo & thermo
Definition: createFields.H:28
pimpleNoLoopControl & pimple
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
Definition: fvcDdt.C:45
Foam::fvConstraints & fvConstraints
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
bool constrain(fvMatrix< Type > &eqn) const
Apply constraints to an equation.
Foam::fvModels & fvModels
messageStream Info
const dimensionedScalar e
Elementary charge.
Definition: doubleScalar.H:105