zhyang-dev/OpenFOAM-6
1
0
Fork 0
mirror of https://github.com/OpenFOAM/OpenFOAM-6.git synced 2025-12-08 06:57:46 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

snappyHexMesh::layerParameters: Increased the maximum number of layer expansion iterations to 100

Browse Source 
Resolves bug-report https://bugs.openfoam.org/view.php?id=2938
This commit is contained in:
Henry Weller
2018-05-22 20:03:27 +01:00
parent 314e464f7d
commit 5caa5a3e9b
1 changed files with 24 additions and 18 deletions
Download Patch File Download Diff File Expand all files Collapse all files

40
src/mesh/snappyHexMesh/snappyHexMeshDriver/layerParameters/layerParameters.C
View File

@ -45,24 +45,25 @@ Foam::scalar Foam::layerParameters::layerExpansionRatio
{ {
if (n <= 1) if (n <= 1)
{ {
return 1.0; return 1;
} }
const label maxIters = 20;
const scalar tol = 1e-8; const scalar tol = 1e-8;
if (mag(n - totalOverFirst) < tol) if (mag(n - totalOverFirst) < tol)
{ {
return 1.0; return 1;
} }
const label maxIters = 100;
// Calculate the bounds of the solution // Calculate the bounds of the solution
scalar minR; scalar minR;
scalar maxR; scalar maxR;
if (totalOverFirst < n) if (totalOverFirst < n)
{ {
minR = 0.0; minR = 0;
maxR = pow(totalOverFirst/n, 1/(n-1)); maxR = pow(totalOverFirst/n, 1/(n-1));
} }
else else
@ -77,9 +78,9 @@ Foam::scalar Foam::layerParameters::layerExpansionRatio
for (label i = 0; i < maxIters; ++i) for (label i = 0; i < maxIters; ++i)
{ {
const scalar prevr = r; const scalar prevr = r;
const scalar fx = pow(r, n) - totalOverFirst*r - (1 - totalOverFirst); const scalar fx = pow(r, n) - totalOverFirst*r - (1 - totalOverFirst);
const scalar dfx = n*pow(r, n - 1) - totalOverFirst; const scalar dfx = n*pow(r, n - 1) - totalOverFirst;
r -= fx/dfx; r -= fx/dfx;
if (mag(r - prevr) < tol) if (mag(r - prevr) < tol)
@ -87,6 +88,7 @@ Foam::scalar Foam::layerParameters::layerExpansionRatio
break; break;
} }
} }
return r; return r;
} }
@ -274,7 +276,7 @@ Foam::layerParameters::layerParameters
forAllConstIter(labelHashSet, patchIDs, patchiter) forAllConstIter(labelHashSet, patchIDs, patchiter)
{ {
label patchi = patchiter.key(); const label patchi = patchiter.key();
numLayers_[patchi] = numLayers_[patchi] =
readLabel(layerDict.lookup("nSurfaceLayers")); readLabel(layerDict.lookup("nSurfaceLayers"));
@ -396,8 +398,8 @@ Foam::scalar Foam::layerParameters::layerThickness
else else
{ {
return firstLayerThickess return firstLayerThickess
*(1.0 - pow(expansionRatio, nLayers)) *(1 - pow(expansionRatio, nLayers))
/(1.0 - expansionRatio); /(1 - expansionRatio);
} }
} }
break; break;
@ -410,10 +412,11 @@ Foam::scalar Foam::layerParameters::layerThickness
} }
else else
{ {
scalar invExpansion = 1.0 / expansionRatio; const scalar invExpansion = 1.0/expansionRatio;
return finalLayerThickess return finalLayerThickess
*(1.0 - pow(invExpansion, nLayers)) *(1 - pow(invExpansion, nLayers))
/(1.0 - invExpansion); /(1 - invExpansion);
} }
} }
break; break;
@ -504,7 +507,7 @@ Foam::scalar Foam::layerParameters::firstLayerThickness
case FINAL_AND_TOTAL: case FINAL_AND_TOTAL:
{ {
scalar r = layerExpansionRatio const scalar r = layerExpansionRatio
( (
nLayers, nLayers,
firstLayerThickess, firstLayerThickess,
@ -512,18 +515,21 @@ Foam::scalar Foam::layerParameters::firstLayerThickness
totalThickness, totalThickness,
expansionRatio expansionRatio
); );
return finalLayerThickess/pow(r, nLayers-1); return finalLayerThickess/pow(r, nLayers-1);
} }
break; break;
case TOTAL_AND_EXPANSION: case TOTAL_AND_EXPANSION:
{ {
scalar r = finalLayerThicknessRatio const scalar r = finalLayerThicknessRatio
( (
nLayers, nLayers,
expansionRatio expansionRatio
); );
scalar finalThickness = r*totalThickness;
const scalar finalThickness = r*totalThickness;
return finalThickness/pow(expansionRatio, nLayers-1); return finalThickness/pow(expansionRatio, nLayers-1);
} }
break; break;
@ -554,13 +560,13 @@ Foam::scalar Foam::layerParameters::finalLayerThicknessRatio
{ {
return return
pow(expansionRatio, nLayers - 1) pow(expansionRatio, nLayers - 1)
*(1.0 - expansionRatio) *(1 - expansionRatio)
/(1.0 - pow(expansionRatio, nLayers)); /(1 - pow(expansionRatio, nLayers));
} }
} }
else else
{ {
return 0.0; return 0;
} }
} }