45 int main(
int argc,
char *argv[])
52 #include "createControls.H" 53 #include "createFields.H" 54 #include "createFieldRefs.H" 58 Info<<
"\nCalculating displacement field\n" <<
endl;
67 scalar initialResidual = 0;
73 if (
thermo.thermalStress())
101 if (
thermo.thermalStress())
108 initialResidual = DEqn.solve().max().initialResidual();
112 divSigmaExp =
fvc::div(DEqn.flux());
136 #include "calculateStress.H" 138 Info<<
"ExecutionTime = " <<
runTime.elapsedCpuTime() <<
" s" 139 <<
" 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)
fluidReactionThermo & thermo
dimensionedScalar lambda(laminarTransport.lookup("lambda"))
GeometricField< tensor, fvPatchField, volMesh > volTensorField
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Switch compactNormalStress(stressControl.lookup("compactNormalStress"))
virtual void correct()
Correct the fvModels.
Ostream & endl(Ostream &os)
Add newline and flush stream.
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
Foam::fvConstraints & fvConstraints
tmp< fvMatrix< Type > > d2dt2(const GeometricField< Type, fvPatchField, volMesh > &field) const
Return source for an equation with a second time derivative.
dimensionedScalar tr(const dimensionedSphericalTensor &dt)
const dimensionedScalar mu
Atomic mass unit.
bool constrain(fvMatrix< Type > &eqn) const
Apply constraints to an equation.
Foam::fvModels & fvModels
fvMatrix< vector > fvVectorMatrix
Execute application functionObjects to post-process existing results.