applications/solvers/solidMechanics/elasticAcpSolidFoam/createFields.H

   1     Info<< "Reading field U\n" << endl;
   2     volVectorField U
   3     (
   4         IOobject
   5         (
   6             "U",
   7             runTime.timeName(),
   8             mesh,
   9             IOobject::MUST_READ,
  10             IOobject::AUTO_WRITE
  11         ),
  12         mesh
  13     );
  14
  15     volTensorField gradU //= fvc::grad(U);
  16     (
  17         IOobject
  18         (
  19             "grad(U)",
  20             runTime.timeName(),
  21             mesh,
  22             IOobject::NO_READ,
  23             IOobject::NO_WRITE
  24         ),
  25         mesh,
  26         dimensionedTensor("zero", dimless, tensor::zero)
  27     );
  28     surfaceVectorField snGradU // = fvc::snGrad(U);
  29     (
  30         IOobject
  31         (
  32             "snGrad(U)",
  33             runTime.timeName(),
  34             mesh,
  35             IOobject::NO_READ,
  36             IOobject::NO_WRITE
  37         ),
  38         mesh,
  39         dimensionedVector("zero", dimless, vector::zero)
  40     );
  41
  42     volVectorField V
  43     (
  44         IOobject
  45         (
  46             "V",
  47             runTime.timeName(),
  48             mesh,
  49             IOobject::READ_IF_PRESENT,
  50             IOobject::AUTO_WRITE
  51         ),
  52         fvc::ddt(U)
  53     );
  54
  55     volSymmTensorField epsilon
  56     (
  57         IOobject
  58         (
  59             "epsilon",
  60             runTime.timeName(),
  61             mesh,
  62             IOobject::READ_IF_PRESENT,
  63             IOobject::AUTO_WRITE
  64         ),
  65         mesh,
  66         dimensionedSymmTensor("zero", dimless, symmTensor::zero)
  67     );
  68
  69     volSymmTensorField sigma
  70     (
  71         IOobject
  72         (
  73             "sigma",
  74             runTime.timeName(),
  75             mesh,
  76             IOobject::READ_IF_PRESENT,
  77             IOobject::AUTO_WRITE
  78         ),
  79         mesh,
  80         dimensionedSymmTensor("zero", dimForce/dimArea, symmTensor::zero)
  81     );
  82
  83       volVectorField divSigmaExp
  84       (
  85         IOobject
  86         (
  87          "divSigmaExp",
  88          runTime.timeName(),
  89          mesh,
  90          IOobject::NO_READ,
  91          IOobject::NO_WRITE
  92          ),
  93         mesh,
  94         dimensionedVector("zero", dimForce/dimVolume, vector::zero)
  95        );
  96
  97     constitutiveModel rheology(sigma, U);
  98
  99     surfaceVectorField traction
 100     (
 101         IOobject
 102         (
 103             "traction",
 104             runTime.timeName(),
 105             mesh,
 106             IOobject::NO_READ,
 107             IOobject::NO_WRITE
 108         ),
 109         mesh,
 110         dimensionedVector("zero", dimForce/dimArea, vector::zero)
 111     );
 112
 113   // for aitken relaxation
 114    volVectorField aitkenDelta
 115     (
 116         IOobject
 117         (
 118             "aitkenDelta",
 119             runTime.timeName(),
 120             mesh,
 121             IOobject::NO_READ,
 122             IOobject::NO_WRITE
 123         ),
 124         mesh,
 125         dimensionedVector("zero", dimLength, vector::zero)
 126     );
 127 // aitken relaxation factor
 128 scalar aitkenInitialRes = 1.0;
 129 scalar aitkenTheta = 0.1;
 130
 131
 132 //    volVectorField resid
 133 //     (
 134 //         IOobject
 135 //         (
 136 //             "resid",
 137 //             runTime.timeName(),
 138 //             mesh,
 139 //             IOobject::NO_READ,
 140 //             IOobject::AUTO_WRITE
 141 //         ),
 142 //         mesh,
 143 //      dimensionedVector("zero", dimless, vector::zero)
 144 //     );