kineticEnergyLimiter.H
Go to the documentation of this file.
1if (!(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(Dcorr.dimensions(), Zero);
18 }
19}
engineTime & runTime