FreeFOAM The Cross-Platform CFD Toolkit
kineticEnergyLimiter.H
Go to the documentation of this file.
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 }
20 
21 // ************************ vim: set sw=4 sts=4 et: ************************ //