ENH: maxDeltaxyz - simplified the delta calc and reset the default coeff to 2

This commit is contained in:
Andrew Heather
2017-12-15 17:25:22 +00:00
parent 627f312776
commit 9c1bb8d707

View File

@ -64,14 +64,14 @@ void Foam::LESModels::maxDeltaxyz::calcDelta()
const point& fc = faceC[facei]; const point& fc = faceC[facei];
const vector& n = faceN[facei]; const vector& n = faceN[facei];
scalar tmp = magSqr(n*(n & (fc - cc))); scalar tmp = mag(n & (fc - cc));
if (tmp > deltaMaxTmp) if (tmp > deltaMaxTmp)
{ {
deltaMaxTmp = tmp; deltaMaxTmp = tmp;
} }
} }
hmax[celli] = deltaCoeff_*Foam::sqrt(deltaMaxTmp); hmax[celli] = deltaCoeff_*deltaMaxTmp;
} }
if (nD == 3) if (nD == 3)
@ -113,7 +113,7 @@ Foam::LESModels::maxDeltaxyz::maxDeltaxyz
dict.optionalSubDict(type() + "Coeffs").lookupOrDefault<scalar> dict.optionalSubDict(type() + "Coeffs").lookupOrDefault<scalar>
( (
"deltaCoeff", "deltaCoeff",
1 2
) )
) )
{ {