applications/solvers/solidMechanics/elasticOrthoSolidFoam/calculateDivSigmaExp.H

   1 if(divSigmaExpMethod == "standard")
   2   {
   3     //- calculating the full gradient has good convergence and no high freq oscillations
   4     divSigmaExp =
   5       fvc::div(
   6                (C && symm(gradU))
   7                - (K & gradU),
   8                "div(sigma)"
   9                );
  10   }
  11  else if(divSigmaExpMethod == "surface")
  12    {
  13      //- this form seems to have the best convergence
  14      divSigmaExp =
  15        fvc::div(mesh.magSf()*
  16                 (
  17                  (n&(Cf && fvc::interpolate(symm(gradU))))
  18                  - (n&(Kf & fvc::interpolate(gradU)))
  19                  )
  20                 );
  21    }
  22  else if(divSigmaExpMethod == "laplacian")
  23    {
  24      //- can cause high freq oscillations and slow convergence
  25      divSigmaExp =
  26        fvc::div(C && symm(epsilon), "div(sigma)")
  27        - fvc::laplacian(K,U, "laplacian(K,U)");
  28    }
  29  else
  30    {
  31      FatalError << "divSigmaExp method " << divSigmaExpMethod << " not found!" << endl;
  32    }