applications/solvers/solidMechanics/icoFsiElasticNonLinULSolidFoam/solveSolidBackward.H

   1     {
   2 #       include "readSolidMechanicsControls.H"
   3
   4         int iCorr = 0;
   5         lduMatrix::solverPerformance solverPerf;
   6         scalar initialResidual = 0;
   7
   8         lduMatrix::debug = 0;
   9
  10 #       include "backwardCoeffs.H"
  11
  12         do
  13         {
  14             DU.storePrevIter();
  15
  16             fvVectorMatrix DUEqn
  17             (
  18                 Cn*rho*fvm::ddt(DU)
  19               - Co*rho*DV.oldTime()
  20               + Coo*rho*DV.oldTime().oldTime()
  21              ==
  22                 fvm::laplacian(2*mu + lambda, DU, "laplacian(DDU,DU)")
  23               - fvc::laplacian(mu + lambda, DU, "laplacian(DDU,DU)")
  24               + fvc::div
  25                 (
  26                     mu*gradDU.T()
  27                   + lambda*(I*tr(gradDU))
  28                   + mu*(gradDU&gradDU.T())
  29                   + 0.5*lambda*(I*tr(gradDU & gradDU.T()))
  30                   + (sigma & DF.T())
  31                   + (DSigma & DF.T()),
  32                     "div(sigma)"
  33                 )
  34             );
  35
  36             solverPerf = DUEqn.solve();
  37
  38             DU.relax();
  39
  40             if(iCorr == 0)
  41             {
  42                 initialResidual = solverPerf.initialResidual();
  43             }
  44
  45             gradDU = fvc::grad(DU);
  46
  47             DF = gradDU.T();
  48
  49 #           include "calculateDSigma.H"
  50         }
  51         while
  52         (
  53             solverPerf.initialResidual() > convergenceTolerance
  54          && ++iCorr < nCorr
  55         );
  56
  57         Info << "Solving for " << DU.name()
  58             << ", Initial residual = " << initialResidual
  59             << ", Final residual = " << solverPerf.initialResidual()
  60             << ", No outer iterations " << iCorr << endl;
  61
  62         DV = fvc::ddt(DU);
  63
  64         lduMatrix::debug = 1;
  65     }
  66