File: kineticEnergyLimiter.H

package info (click to toggle)
openfoam 4.1%2Bdfsg1-1
  • links: PTS, VCS
  • area: main
  • in suites: stretch
  • size: 163,028 kB
  • ctags: 58,990
  • sloc: cpp: 830,760; sh: 10,227; ansic: 8,215; xml: 745; lex: 437; awk: 194; sed: 91; makefile: 77; python: 18
file content (19 lines) | stat: -rw-r--r-- 704 bytes parent folder | download | duplicates (2) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
 
if (!(runTime.timeIndex() % 5))
{
    volScalarField kineticEnergy(magSqr(Dcorr));
    dimensionedScalar intKineticEnergy = fvc::domainIntegrate(kineticEnergy);
    Info<< "kineticEnergy = " << intKineticEnergy.value() << endl;

    volScalarField kineticPower((Dcorr - Dcorr.oldTime()) & Dcorr);
    dimensionedScalar intKineticPower = fvc::domainIntegrate(kineticPower);
    Info<< "kineticPower = " << intKineticPower.value() << endl;

    scalar smi = intKineticPower.value()/(intKineticEnergy.value() + VSMALL);
    Info<< "smi = " << smi << endl;

    if (smi < -SMALL)
    {
        Info<< "Resetting Dcorr to 0" << endl;
        Dcorr == dimensionedVector("0", Dcorr.dimensions(), Zero);
    }
}