applications/solvers/stressAnalysis/solidEquilibriumDisplacementFoam/kineticEnergyLimiter.H

   1 if (!(runTime.timeIndex() % 5))
   2 {
   3     volScalarField kineticEnergy = magSqr(Dcorr);
   4     dimensionedScalar intKineticEnergy = fvc::domainIntegrate(kineticEnergy);
   5     Info<< "kineticEnergy = " << intKineticEnergy.value() << endl;
   6
   7     volScalarField kineticPower = (Dcorr - Dcorr.oldTime()) & Dcorr;
   8     dimensionedScalar intKineticPower = fvc::domainIntegrate(kineticPower);
   9     Info<< "kineticPower = " << intKineticPower.value() << endl;
  10
  11     scalar smi = intKineticPower.value()/(intKineticEnergy.value() + VSMALL);
  12     Info<< "smi = " << smi << endl;
  13
  14     if (smi < -SMALL)
  15     {
  16         Info<< "Resetting Dcorr to 0" << endl;
  17         Dcorr == dimensionedVector("0", Dcorr.dimensions(), vector::zero);
  18     }
  19 }