createFields.H
Go to the documentation of this file.
2
3Info<< "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
17Info<< "Reading/setting displacement correction field Dcorr\n" << endl;
19(
20 IOobject
21 (
22 "Dcorr",
23 runTime.timeName(),
24 mesh
25 ),
26 D
27);
28Dcorr *= 0.0;
29
30Info<< "Calculating stress field sigmaD\n" << endl;
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
44volTensorField 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);
const volScalarField & mu
dynamicFvMesh & mesh
engineTime & runTime
GeometricField< vector, fvPatchField, volMesh > volVectorField
Definition: volFieldsFwd.H:83
dimensionedScalar tr(const dimensionedSphericalTensor &dt)
dimensionedSymmTensor twoSymm(const dimensionedSymmTensor &dt)
messageStream Info
Information stream (stdout output on master, null elsewhere)
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:372
GeometricField< tensor, fvPatchField, volMesh > volTensorField
Definition: volFieldsFwd.H:87
GeometricField< symmTensor, fvPatchField, volMesh > volSymmTensorField
Definition: volFieldsFwd.H:86
const dimensionedScalar & D
dimensionedScalar lambda("lambda", dimTime/sqr(dimLength), laminarTransport)