createFields.H
Go to the documentation of this file.
2 
3 Info<< "Reading displacement field D\n" << endl;
5 (
6  IOobject
7  (
8  "D",
9  runTime.timeName(),
10  mesh,
11  IOobject::MUST_READ,
12  IOobject::AUTO_WRITE
13  ),
14  mesh
15 );
16 
17 Info<< "Reading/setting displacement correction field Dcorr\n" << endl;
18 volVectorField Dcorr
19 (
20  IOobject
21  (
22  "Dcorr",
23  runTime.timeName(),
24  mesh
25  ),
26  D
27 );
28 Dcorr = Zero;
29 
30 Info<< "Calculating stress field sigmaD\n" << endl;
31 volSymmTensorField sigmaD
32 (
33  IOobject
34  (
35  "sigmaD",
36  runTime.timeName(),
37  mesh,
38  IOobject::NO_READ,
39  IOobject::NO_WRITE
40  ),
41  mu*twoSymm(fvc::grad(D)) + (lambda*I)*tr(fvc::grad(D))
42 );
43 
44 volTensorField sigmaExp
45 (
46  IOobject
47  (
48  "sigmaExp",
49  runTime.timeName(),
50  mesh,
51  IOobject::NO_READ,
52  IOobject::NO_WRITE
53  ),
54  (lambda - mu)*fvc::grad(Dcorr) + mu*fvc::grad(Dcorr)().T()
55  + (lambda*I)*tr(fvc::grad(Dcorr))
56 );
GeometricField< symmTensor, fvPatchField, volMesh > volSymmTensorField
Definition: volFieldsFwd.H:62
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcGrad.C:52
GeometricField< tensor, fvPatchField, volMesh > volTensorField
Definition: volFieldsFwd.H:63
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
dimensionedScalar lambda(viscosity->lookup("lambda"))
GeometricField< vector, fvPatchField, volMesh > volVectorField
Definition: volFieldsFwd.H:59
fvMesh & mesh
dimensionedSymmTensor twoSymm(const dimensionedSymmTensor &dt)
static const Identity< scalar > I
Definition: Identity.H:93
static const zero Zero
Definition: zero.H:97
dimensionedScalar tr(const dimensionedSphericalTensor &dt)
volScalarField mu(IOobject("mu", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), viscosity->nu() *rhoInfValue)
const volScalarField & T
messageStream Info