applications/solvers/solidMechanics/elasticIncrAcpSolidFoam/calculateDivDSigmaExp.H

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