applications/solvers/solidMechanics/elasticOrthoAcpSolidFoam/calculateDivSigmaExp.H

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