calculateStress.H
Go to the documentation of this file.
1 if (runTime.writeTime())
2 {
4  (
5  IOobject
6  (
7  "sigma",
8  runTime.timeName(),
9  mesh,
10  IOobject::NO_READ,
11  IOobject::AUTO_WRITE
12  ),
13  sigmaD
14  );
15 
16  if (thermo.thermalStress())
17  {
18  sigma = sigma - I*(threeKalpha*thermo.T());
19  }
20 
21  volScalarField sigmaEq
22  (
23  IOobject
24  (
25  "sigmaEq",
26  runTime.timeName(),
27  mesh,
28  IOobject::NO_READ,
29  IOobject::AUTO_WRITE
30  ),
31  sqrt((3.0/2.0)*magSqr(dev(sigma)))
32  );
33 
34  Info<< "Max sigmaEq = " << max(sigmaEq).value()
35  << endl;
36 
37  runTime.write();
38 }
GeometricField< symmTensor, fvPatchField, volMesh > volSymmTensorField
Definition: volFieldsFwd.H:61
fluidReactionThermo & thermo
Definition: createFields.H:28
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
engineTime & runTime
dimensionedScalar sqrt(const dimensionedScalar &ds)
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:57
const dimensionedScalar sigma
Stefan-Boltzmann constant: default SI units: [W/m^2/K^4].
dynamicFvMesh & mesh
dimensionedSymmTensor dev(const dimensionedSymmTensor &dt)
static const Identity< scalar > I
Definition: Identity.H:93
dimensioned< scalar > magSqr(const dimensioned< Type > &)
messageStream Info