applications/solvers/solidMechanics/elasticIncrAcpSolidFoam/setDeltaT.H

   1 if (dynamicTimeStep)
   2 {
   3     if
   4     (
   5         //(maxEffTraction < 0.999*CzLaw.sigmaMax().value())
   6         (returnReduce(maxEffTractionFraction, maxOp<scalar>()) < 0.99)
   7       //&& (cohesivePatchU.size() == 0)
   8       && (mag(runTime.deltaT().value() - deltaTmax) < SMALL)
   9     )
  10     {
  11         runTime.setDeltaT(deltaTmax);
  12     }
  13     else
  14     {
  15         scalar newDeltaT = deltaTmin;
  16
  17         if (newDeltaT/runTime.deltaT().value() < 0.5)
  18         {
  19             newDeltaT = 0.5*runTime.deltaT().value();
  20             Info << "Reducing time step" << nl;
  21         }
  22
  23         runTime.setDeltaT(newDeltaT);
  24     }
  25
  26     Pout << "Current time step size: "
  27         << runTime.deltaT().value() << " s" << endl;
  28
  29     scalar maxDT = runTime.deltaT().value();
  30
  31     if(mag(returnReduce(maxDT, maxOp<scalar>()) - runTime.deltaT().value()) > SMALL)
  32     {
  33         FatalError << "Processors have different time-steps!"
  34             << exit(FatalError);
  35     }
  36 }