songtianlun/openfoam
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

ENH: cyclicAMI - added message if incorrect angle specified for rotational case

Browse Source
This commit is contained in:
andy 2014-01-21 10:59:24 +00:00
parent e7201c0326
commit fdc785d671
1 changed files with 30 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

30
src/meshTools/AMIInterpolation/patches/cyclicAMI/cyclicAMIPolyPatch/cyclicAMIPolyPatch.C
View File

@ -186,6 +186,36 @@ void Foam::cyclicAMIPolyPatch::calcTransforms
rotationAngle_ *= -1;
}
scalar areaError =
max(errorPos, errorNeg)/(mag(area1) + ROOTVSMALL);
if (areaError > matchTolerance())
{
WarningIn
(
"void Foam::cyclicAMIPolyPatch::calcTransforms"
"("
"const primitivePatch&, "
"const pointField&, "
"const vectorField&, "
"const pointField&, "
"const vectorField&"
")"
)
<< "Patch areas are not consistent within "
<< 100*matchTolerance()
<< " % indicating a possible error in the specified "
<< "angle of rotation" << nl
<< " owner patch : " << name() << nl
<< " neighbour patch : " << neighbPatch().name()
<< nl
<< " angle : "
<< radToDeg(rotationAngle_) << " deg" << nl
<< " area error : " << 100*areaError << " %"
<< " match tolerance : " << matchTolerance()
<< endl;
}
if (debug)
{
scalar theta = radToDeg(rotationAngle_);