3 volScalarField kineticEnergy(magSqr(Dcorr));
4 dimensionedScalar intKineticEnergy = fvc::domainIntegrate(kineticEnergy);
5 Info<<
"kineticEnergy = " << intKineticEnergy.value() << endl;
7 volScalarField kineticPower((Dcorr - Dcorr.oldTime()) & Dcorr);
8 dimensionedScalar intKineticPower = fvc::domainIntegrate(kineticPower);
9 Info<<
"kineticPower = " << intKineticPower.value() << endl;
11 scalar smi = intKineticPower.value()/(intKineticEnergy.value() + VSMALL);
12 Info<<
"smi = " << smi << endl;
16 Info<<
"Resetting Dcorr to 0" << endl;
17 Dcorr == dimensionedVector(Dcorr.dimensions(), Zero);