44 int main(
int argc,
char *argv[])
51 #include "createControls.H" 52 #include "createFields.H" 53 #include "createFieldRefs.H" 57 Info<<
"\nCalculating displacement field\n" <<
endl;
66 scalar initialResidual = 0;
70 if (
thermo.thermalStress())
97 if (
thermo.thermalStress())
104 initialResidual = DEqn.solve().max().initialResidual();
108 divSigmaExp =
fvc::div(DEqn.flux());
132 #include "calculateStress.H" 134 Info<<
"ExecutionTime = " <<
runTime.elapsedCpuTime() <<
" s" 135 <<
" ClockTime = " <<
runTime.elapsedClockTime() <<
" s"
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
fvMatrix< scalar > fvScalarMatrix
tmp< GeometricField< Type, fvPatchField, volMesh > > d2dt2(const GeometricField< Type, fvPatchField, volMesh > &vf)
const dimensionedScalar & kappa
Coulomb constant: default SI units: [N.m2/C2].
GeometricField< tensor, fvPatchField, volMesh > volTensorField
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Switch compactNormalStress(stressControl.lookup("compactNormalStress"))
Ostream & endl(Ostream &os)
Add newline and flush stream.
rhoReactionThermo & thermo
tmp< GeometricField< Type, fvPatchField, volMesh > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
dimensionedSymmTensor twoSymm(const dimensionedSymmTensor &dt)
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
GeometricField< scalar, fvPatchField, volMesh > volScalarField
static const Identity< scalar > I
dimensionedScalar tr(const dimensionedSphericalTensor &dt)
const dimensionedScalar & mu
Atomic mass unit.
fvMatrix< vector > fvVectorMatrix
dimensionedScalar lambda(laminarTransport.lookup("lambda"))
scalar Cp(const scalar p, const scalar T) const
Execute application functionObjects to post-process existing results.