applications/solvers/solidMechanics/elasticThermalSolidFoam/calculateSigmaExp.H

   1 if(sigmaExpMethod == "standard")
   2 {
   3     sigmaExp = fvc::div
   4       (
   5        mu*gradU.T() + lambda*(I*tr(gradU)) - (mu + lambda)*gradU,
   6        "div(sigma)"
   7        );
   8 }
   9 else if(sigmaExpMethod == "surface")
  10 {
  11      sigmaExp = fvc::div
  12      (
  13          muf*(mesh.Sf() & fvc::interpolate(gradU.T()))
  14        + lambdaf*(mesh.Sf() & I*fvc::interpolate(tr(gradU)))
  15        - (muf + lambdaf)*(mesh.Sf() & fvc::interpolate(gradU))
  16     );
  17 }
  18 else if(sigmaExpMethod == "decompose")
  19 {
  20      surfaceTensorField shearGradU = ((I - n*n)&fvc::interpolate(gradU));
  21
  22      sigmaExp = fvc::div
  23      (
  24          mesh.magSf()
  25         *
  26          (
  27             - (muf + lambdaf)*(fvc::snGrad(U)&(I - n*n))
  28             + lambdaf*tr(shearGradU&(I - n*n))*n
  29             + muf*(shearGradU&n)
  30          )
  31     );
  32 }
  33 else if(sigmaExpMethod == "expLaplacian")
  34 {
  35     sigmaExp =
  36        - fvc::laplacian(mu + lambda, U, "laplacian(DU,U)")
  37        + fvc::div
  38          (
  39              mu*gradU.T()
  40            + lambda*(I*tr(gradU)),
  41              "div(sigma)"
  42          );
  43 }
  44 else
  45 {
  46      FatalError << "sigmaExp method " << sigmaExpMethod << " not found!" << endl;
  47 }