42 void Foam::IDDESDelta::calcDelta()
49 wallDistReflection wallNorm(
mesh());
53 tmp<volScalarField> faceToFacenMax
74 scalar deltaMaxTmp = 0.0;
77 const vector nCell = n[cellI];
80 label faceI = cFaces[cFaceI];
82 scalar tmp = (faceCentre - faceCentreTwo) & nCell;
83 if (tmp > deltaMaxTmp)
88 faceToFacenMax()[cellI] = deltaMaxTmp;
99 max(cw_*wallDist(
mesh()).
y(), cw_*hmax),
108 <<
"Case is 2D, LES is not strictly applicable\n"
115 max(
max(cw_*wallDist(
mesh()).
y(), cw_*hmax), faceToFacenMax()),
122 <<
"Case is not 3D or 2D, LES is not strictly applicable"
130 Foam::IDDESDelta::IDDESDelta
161 if (mesh_.changing())