exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Calculation of cell aspect ratio sensitised to 1- and 2-D cases

Browse source
This commit is contained in:
Hrvoje Jasak 2012-03-27 11:53:21 +01:00
parent 0ee9ddb2b2
commit 190c4f8eae
1 changed files with 20 additions and 7 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

27
src/OpenFOAM/meshes/primitiveMesh/primitiveMeshCheck/primitiveMeshCheck.C
View file

@ -258,13 +258,26 @@ bool Foam::primitiveMesh::checkClosedCells
// Calculate the aspect ration as the maximum of Cartesian component // Calculate the aspect ration as the maximum of Cartesian component
// aspect ratio to the total area hydraulic area aspect ratio // aspect ratio to the total area hydraulic area aspect ratio
scalar aspectRatio = max scalar minCmpt = VGREAT;
( scalar maxCmpt = -VGREAT;
cmptMax(sumMagClosed[cellI]) for (direction dir = 0; dir < vector::nComponents; dir++)
/(cmptMin(sumMagClosed[cellI]) + VSMALL), {
1.0/6.0*cmptSum(sumMagClosed[cellI])/ if (meshD[dir] == 1)
Foam::pow(Foam::max(vols[cellI], SMALL), 2.0/3.0) {
); minCmpt = min(minCmpt, sumMagClosed[cellI][dir]);
maxCmpt = max(maxCmpt, sumMagClosed[cellI][dir]);
}
}
scalar aspectRatio = maxCmpt/(minCmpt + VSMALL);
if (nDims == 3)
{
aspectRatio = max
(
aspectRatio,
1.0/6.0*cmptSum(sumMagClosed[cellI])/pow(vols[cellI], 2.0/3.0)
);
}
maxAspectRatio = max(maxAspectRatio, aspectRatio); maxAspectRatio = max(maxAspectRatio, aspectRatio);