applications/solvers/solidMechanics/elasticOrthoAcpSolidFoam/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   }