motionSmootherAlgoCheck::checkMesh: improved minVol test and removed unused tests

motionSmootherAlgoCheck::checkMesh is used by snappyHexMesh to check the mesh
after snapping and morphing.  The minVol test which checks for collapsed cells
is now relative to the cube of the minimum bounding box length so that it is
less dependent on the size of the geometry and less likely to need changing for
very small geometries.

The default value is set in
etc/caseDicts/mesh/generation/meshQualityDict
etc/caseDicts/mesh/generation/meshQualityDict.cfg

//- Minimum cell pyramid volume relative to min bounding box length^3
//  Set to a fraction of the smallest cell volume expected.
//  Set to very negative number (e.g. -1e30) to disable.
minVol 1e-10;

The unused minArea and minTriangleTwist tests have been removed
This commit is contained in:
Henry Weller
2022-02-24 16:44:34 +00:00
parent a34ddef943
commit 32e7b52c81
20 changed files with 110 additions and 301 deletions

View File

@ -433,8 +433,8 @@ meshQualityControls
// Set to 180 to disable. // Set to 180 to disable.
maxConcave 80; maxConcave 80;
//- Minimum pyramid volume. Is absolute volume of cell pyramid. //- Minimum cell pyramid volume relative to min bounding box length^3
// Set to a sensible fraction of the smallest cell volume expected. // Set to a fraction of the smallest cell volume expected.
// Set to very negative number (e.g. -1e30) to disable. // Set to very negative number (e.g. -1e30) to disable.
minVol -1e30; minVol -1e30;
@ -449,9 +449,6 @@ meshQualityControls
// 1 = regular tet // 1 = regular tet
minTetQuality 1e-30; minTetQuality 1e-30;
//- Minimum absolute face area. Set to <0 to disable.
minArea -1;
//- Minimum face twist. Set to <-1 to disable. dot product of face normal //- Minimum face twist. Set to <-1 to disable. dot product of face normal
// and face centre triangles normal // and face centre triangles normal
minTwist 0.02; minTwist 0.02;
@ -465,9 +462,6 @@ meshQualityControls
//- minVolRatio (0 -> 1) //- minVolRatio (0 -> 1)
minVolRatio 0.01; minVolRatio 0.01;
// must be >0 for Fluent compatibility
minTriangleTwist -1;
} }

View File

@ -25,10 +25,10 @@ maxInternalSkewness 4;
// Set to 180 to disable. // Set to 180 to disable.
maxConcave 80; maxConcave 80;
//- Minimum pyramid volume. Is absolute volume of cell pyramid. //- Minimum cell pyramid volume relative to min bounding box length^3
// Set to a sensible fraction of the smallest cell volume expected. // Set to a fraction of the smallest cell volume expected.
// Set to very negative number (e.g. -1e30) to disable. // Set to very negative number (e.g. -1e30) to disable.
minVol 1e-13; minVol 1e-10;
//- Minimum quality of the tet formed by the face-centre //- Minimum quality of the tet formed by the face-centre
// and variable base point minimum decomposition triangles and // and variable base point minimum decomposition triangles and
@ -39,9 +39,6 @@ minVol 1e-13;
// 1 = regular tet // 1 = regular tet
minTetQuality 1e-15; minTetQuality 1e-15;
//- Minimum face area. Set to <0 to disable.
minArea -1;
//- Minimum face twist. Set to <-1 to disable. dot product of face normal //- Minimum face twist. Set to <-1 to disable. dot product of face normal
// and face centre triangles normal // and face centre triangles normal
minTwist 0.02; minTwist 0.02;
@ -61,20 +58,9 @@ minFaceWeight 0.05;
//- Volume ratio of neighbouring cells (0 -> 1) //- Volume ratio of neighbouring cells (0 -> 1)
minVolRatio 0.01; minVolRatio 0.01;
//- Per triangle normal compared to average normal. Like face twist
// but now per (face-centre decomposition) triangle. Must be >0 for Fluent
// compatibility
minTriangleTwist -1;
nSmoothScale 4; nSmoothScale 4;
errorReduction 0.75; errorReduction 0.75;
//- If >0 : preserve cells with all points on the surface if the
// resulting volume after snapping (by approximation) is larger than
// minVolCollapseRatio times old volume (i.e. not collapsed to flat cell).
// If <0 : delete always.
//minVolCollapseRatio 0.1;
// ************************************************************************* // // ************************************************************************* //

View File

@ -14,19 +14,10 @@ maxInternalSkewness 4;
maxConcave 80; maxConcave 80;
// Minimum cell pyramid volume; case dependent minVol 1e-10;
minVol 1e-13;
minTetQuality -1e30; minTetQuality -1e30;
// if >0 : preserve single cells with all points on the surface if the
// resulting volume after snapping (by approximation) is larger than
// minVolCollapseRatio times old volume (i.e. not collapsed to flat cell).
// If <0 : delete always.
//minVolCollapseRatio 0.5;
minArea -1;
minTwist 0.02; minTwist 0.02;
minDeterminant 0.001; minDeterminant 0.001;
@ -35,8 +26,6 @@ minFaceWeight 0.05;
minVolRatio 0.01; minVolRatio 0.01;
minTriangleTwist -1;
nSmoothScale 4; nSmoothScale 4;
errorReduction 0.75; errorReduction 0.75;

View File

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org \\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2011-2020 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2022 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -76,10 +76,6 @@ bool Foam::motionSmootherAlgo::checkMesh
( (
dict.lookup<scalar>("maxConcave", true) dict.lookup<scalar>("maxConcave", true)
); );
const scalar minArea
(
dict.lookup<scalar>("minArea", true)
);
const scalar maxIntSkew const scalar maxIntSkew
( (
dict.lookup<scalar>("maxInternalSkewness", true) dict.lookup<scalar>("maxInternalSkewness", true)
@ -100,10 +96,6 @@ bool Foam::motionSmootherAlgo::checkMesh
( (
dict.lookup<scalar>("minTwist", true) dict.lookup<scalar>("minTwist", true)
); );
const scalar minTriangleTwist
(
dict.lookup<scalar>("minTriangleTwist", true)
);
scalar minFaceFlatness = -1.0; scalar minFaceFlatness = -1.0;
dict.readIfPresent("minFaceFlatness", minFaceFlatness, true); dict.readIfPresent("minFaceFlatness", minFaceFlatness, true);
const scalar minDet const scalar minDet
@ -141,10 +133,12 @@ bool Foam::motionSmootherAlgo::checkMesh
if (minVol > -great) if (minVol > -great)
{ {
const scalar refVol = pow3(mesh.bounds().minDim());
polyMeshCheck::checkFacePyramids polyMeshCheck::checkFacePyramids
( (
report, report,
minVol, minVol*refVol,
mesh, mesh,
mesh.cellCentres(), mesh.cellCentres(),
mesh.points(), mesh.points(),
@ -209,28 +203,6 @@ bool Foam::motionSmootherAlgo::checkMesh
nWrongFaces = nNewWrongFaces; nWrongFaces = nNewWrongFaces;
} }
if (minArea > -small)
{
polyMeshCheck::checkFaceArea
(
report,
minArea,
mesh,
mesh.faceAreas(),
checkFaces,
&wrongFaces
);
label nNewWrongFaces = returnReduce(wrongFaces.size(), sumOp<label>());
Info<< " faces with area < "
<< setw(5) << minArea
<< " m^2 : "
<< nNewWrongFaces-nWrongFaces << endl;
nWrongFaces = nNewWrongFaces;
}
if (maxIntSkew > 0 || maxBounSkew > 0) if (maxIntSkew > 0 || maxBounSkew > 0)
{ {
polyMeshCheck::checkFaceSkewness polyMeshCheck::checkFaceSkewness
@ -333,32 +305,6 @@ bool Foam::motionSmootherAlgo::checkMesh
nWrongFaces = nNewWrongFaces; nWrongFaces = nNewWrongFaces;
} }
if (minTriangleTwist > -1)
{
// Pout<< "Checking triangle twist: dot product of consecutive triangle"
// << " normals resulting from face-centre decomposition" << endl;
polyMeshCheck::checkTriangleTwist
(
report,
minTriangleTwist,
mesh,
mesh.faceAreas(),
mesh.faceCentres(),
mesh.points(),
checkFaces,
&wrongFaces
);
label nNewWrongFaces = returnReduce(wrongFaces.size(), sumOp<label>());
Info<< " faces with triangle twist < "
<< setw(5) << minTriangleTwist
<< " : "
<< nNewWrongFaces-nWrongFaces << endl;
nWrongFaces = nNewWrongFaces;
}
if (minFaceFlatness > -small) if (minFaceFlatness > -small)
{ {
polyMeshCheck::checkFaceFlatness polyMeshCheck::checkFaceFlatness

View File

@ -249,12 +249,10 @@ meshQualityControls
maxConcave 80; maxConcave 80;
minVol 1e-13; minVol 1e-13;
minTetQuality -1; minTetQuality -1;
minArea -1;
minTwist 0.01; minTwist 0.01;
minDeterminant 0.001; minDeterminant 0.001;
minFaceWeight 0.05; minFaceWeight 0.05;
minVolRatio 0.01; minVolRatio 0.01;
minTriangleTwist -1;
nSmoothScale 4; nSmoothScale 4;
errorReduction 0.75; errorReduction 0.75;
relaxed relaxed

View File

@ -388,10 +388,10 @@ meshQualityControls
// Set to 180 to disable. // Set to 180 to disable.
maxConcave 80; maxConcave 80;
//- Minimum pyramid volume. Is absolute volume of cell pyramid. //- Minimum cell pyramid volume relative to min bounding box length^3
// Set to a sensible fraction of the smallest cell volume expected. // Set to a fraction of the smallest cell volume expected.
// Set to very negative number (e.g. -1e30) to disable. // Set to very negative number (e.g. -1e30) to disable.
minVol 1e-13; minVol 1e-10;
//- Minimum quality of the tet formed by the face-centre //- Minimum quality of the tet formed by the face-centre
// and variable base point minimum decomposition triangles and // and variable base point minimum decomposition triangles and
@ -403,9 +403,6 @@ meshQualityControls
// 1 = regular tet // 1 = regular tet
minTetQuality -1; // 1e-30; minTetQuality -1; // 1e-30;
//- Minimum face area. Set to <0 to disable.
minArea -1;
//- Minimum face twist. Set to <-1 to disable. dot product of face normal //- Minimum face twist. Set to <-1 to disable. dot product of face normal
// and face centre triangles normal // and face centre triangles normal
minTwist 0.01; minTwist 0.01;
@ -420,14 +417,11 @@ meshQualityControls
//- minVolRatio (0 -> 1) //- minVolRatio (0 -> 1)
minVolRatio 0.01; minVolRatio 0.01;
// must be >0 for Fluent compatibility
minTriangleTwist -1;
// Advanced // Advanced
//- Number of error distribution iterations //- Number of error distribution iterations
nSmoothScale 4; nSmoothScale 4;
//- Amount to scale back displacement at error points //- Amount to scale back displacement at error points
errorReduction 0.75; errorReduction 0.75;

View File

@ -265,10 +265,10 @@ meshQualityControls
// Set to 180 to disable. // Set to 180 to disable.
maxConcave 80; maxConcave 80;
//- Minimum pyramid volume. Is absolute volume of cell pyramid. //- Minimum cell pyramid volume relative to min bounding box length^3
// Set to a sensible fraction of the smallest cell volume expected. // Set to a fraction of the smallest cell volume expected.
// Set to very negative number (e.g. -1e30) to disable. // Set to very negative number (e.g. -1e30) to disable.
minVol 1e-13; minVol 1e-10;
//- Minimum quality of the tet formed by the face-centre //- Minimum quality of the tet formed by the face-centre
// and variable base point minimum decomposition triangles and // and variable base point minimum decomposition triangles and
@ -279,9 +279,6 @@ meshQualityControls
// 1 = regular tet // 1 = regular tet
minTetQuality 1e-30; minTetQuality 1e-30;
//- Minimum face area. Set to <0 to disable.
minArea -1;
//- Minimum face twist. Set to <-1 to disable. dot product of face normal //- Minimum face twist. Set to <-1 to disable. dot product of face normal
// and face centre triangles normal // and face centre triangles normal
minTwist 0.05; minTwist 0.05;
@ -296,14 +293,11 @@ meshQualityControls
//- minVolRatio (0 -> 1) //- minVolRatio (0 -> 1)
minVolRatio 0.01; minVolRatio 0.01;
// must be >0 for Fluent compatibility
minTriangleTwist -1;
// Advanced // Advanced
//- Number of error distribution iterations //- Number of error distribution iterations
nSmoothScale 4; nSmoothScale 4;
//- Amount to scale back displacement at error points //- Amount to scale back displacement at error points
errorReduction 0.75; errorReduction 0.75;
} }

View File

@ -47,7 +47,6 @@ geometry
// Settings for the castellatedMesh generation. // Settings for the castellatedMesh generation.
castellatedMeshControls castellatedMeshControls
{ {
// Refinement parameters // Refinement parameters
// ~~~~~~~~~~~~~~~~~~~~~ // ~~~~~~~~~~~~~~~~~~~~~
@ -277,14 +276,12 @@ meshQualityControls
//- Minimum projected area v.s. actual area. Set to -1 to disable. //- Minimum projected area v.s. actual area. Set to -1 to disable.
minFlatness 0.5; minFlatness 0.5;
//- Minimum pyramid volume. Is absolute volume of cell pyramid. //- Minimum cell pyramid volume relative to min bounding box length^3
// Set to a sensible fraction of the smallest cell volume expected. // Set to a fraction of the smallest cell volume expected.
// Set to very negative number (e.g. -1e30) to disable. // Set to very negative number (e.g. -1e30) to disable.
minVol 1e-13; minVol 1e-10;
minTetQuality 1e-30;
//- Minimum face area. Set to <0 to disable. minTetQuality 1e-30;
minArea -1;
//- Minimum face twist. Set to <-1 to disable. dot product of face normal //- Minimum face twist. Set to <-1 to disable. dot product of face normal
// and face centre triangles normal // and face centre triangles normal
@ -300,10 +297,6 @@ meshQualityControls
//- minVolRatio (0 -> 1) //- minVolRatio (0 -> 1)
minVolRatio 0.01; minVolRatio 0.01;
// must be >0 for Fluent compatibility
minTriangleTwist -1;
// Advanced // Advanced
//- Number of error distribution iterations //- Number of error distribution iterations

View File

@ -324,10 +324,10 @@ meshQualityControls
// Set to 180 to disable. // Set to 180 to disable.
maxConcave 80; maxConcave 80;
//- Minimum pyramid volume. Is absolute volume of cell pyramid. //- Minimum cell pyramid volume relative to min bounding box length^3
// Set to a sensible fraction of the smallest cell volume expected. // Set to a fraction of the smallest cell volume expected.
// Set to very negative number (e.g. -1e30) to disable. // Set to very negative number (e.g. -1e30) to disable.
minVol 1e-13; minVol 1e-10;
//- Minimum quality of the tet formed by the face-centre //- Minimum quality of the tet formed by the face-centre
// and variable base point minimum decomposition triangles and // and variable base point minimum decomposition triangles and
@ -338,9 +338,6 @@ meshQualityControls
// 1 = regular tet // 1 = regular tet
minTetQuality 1e-30; minTetQuality 1e-30;
//- Minimum face area. Set to <0 to disable.
minArea -1;
//- Minimum face twist. Set to <-1 to disable. dot product of face normal //- Minimum face twist. Set to <-1 to disable. dot product of face normal
// and face centre triangles normal // and face centre triangles normal
minTwist 0.05; minTwist 0.05;
@ -355,25 +352,14 @@ meshQualityControls
//- minVolRatio (0 -> 1) //- minVolRatio (0 -> 1)
minVolRatio 0.01; minVolRatio 0.01;
// must be >0 for Fluent compatibility
minTriangleTwist -1;
//- If >0 : preserve single cells with all points on the surface if the
// resulting volume after snapping (by approximation) is larger than
// minVolCollapseRatio times old volume (i.e. not collapsed to flat cell).
// If <0 : delete always.
// minVolCollapseRatio 0.5;
// Advanced // Advanced
//- Number of error distribution iterations //- Number of error distribution iterations
nSmoothScale 4; nSmoothScale 4;
//- Amount to scale back displacement at error points //- Amount to scale back displacement at error points
errorReduction 0.75; errorReduction 0.75;
// Optional : some meshing phases allow usage of relaxed rules. // Optional : some meshing phases allow usage of relaxed rules.
// See e.g. addLayersControls::nRelaxedIter. // See e.g. addLayersControls::nRelaxedIter.
relaxed relaxed

View File

@ -35,13 +35,10 @@ maxConcave 80;
// 1 = regular tet // 1 = regular tet
minTetQuality 1e-30; minTetQuality 1e-30;
//- Minimum pyramid volume. Is absolute volume of cell pyramid. //- Minimum cell pyramid volume relative to min bounding box length^3
// Set to a sensible fraction of the smallest cell volume expected. // Set to a fraction of the smallest cell volume expected.
// Set to very negative number (e.g. -1e30) to disable. // Set to very negative number (e.g. -1e30) to disable.
minVol 1e-20; minVol 1e-10;
//- Minimum face area. Set to <0 to disable.
minArea -1;
//- Minimum face twist. Set to <-1 to disable. dot product of face normal //- Minimum face twist. Set to <-1 to disable. dot product of face normal
// and face centre triangles normal // and face centre triangles normal
@ -57,8 +54,5 @@ minFaceWeight 0.02;
//- minVolRatio (0 -> 1) //- minVolRatio (0 -> 1)
minVolRatio 0.01; minVolRatio 0.01;
//must be >0 for Fluent compatibility
minTriangleTwist -1;
// ************************************************************************* // // ************************************************************************* //

View File

@ -35,13 +35,10 @@ maxConcave 80;
// 1 = regular tet // 1 = regular tet
minTetQuality 1e-30; minTetQuality 1e-30;
//- Minimum pyramid volume. Is absolute volume of cell pyramid. //- Minimum cell pyramid volume relative to min bounding box length^3
// Set to a sensible fraction of the smallest cell volume expected. // Set to a fraction of the smallest cell volume expected.
// Set to very negative number (e.g. -1e30) to disable. // Set to very negative number (e.g. -1e30) to disable.
minVol 0; minVol 0;
//- Minimum face area. Set to <0 to disable.
minArea -1;
//- Minimum face twist. Set to <-1 to disable. dot product of face normal //- Minimum face twist. Set to <-1 to disable. dot product of face normal
// and face centre triangles normal // and face centre triangles normal
@ -57,8 +54,5 @@ minFaceWeight 0.02;
//- minVolRatio (0 -> 1) //- minVolRatio (0 -> 1)
minVolRatio 0.01; minVolRatio 0.01;
//must be >0 for Fluent compatibility
minTriangleTwist -1;
// ************************************************************************* // // ************************************************************************* //

View File

@ -35,13 +35,10 @@ maxConcave 80;
// 1 = regular tet // 1 = regular tet
minTetQuality 1e-30; minTetQuality 1e-30;
//- Minimum pyramid volume. Is absolute volume of cell pyramid. //- Minimum cell pyramid volume relative to min bounding box length^3
// Set to a sensible fraction of the smallest cell volume expected. // Set to a fraction of the smallest cell volume expected.
// Set to very negative number (e.g. -1e30) to disable. // Set to very negative number (e.g. -1e30) to disable.
minVol 0; minVol 0;
//- Minimum face area. Set to <0 to disable.
minArea -1;
//- Minimum face twist. Set to <-1 to disable. dot product of face normal //- Minimum face twist. Set to <-1 to disable. dot product of face normal
// and face centre triangles normal // and face centre triangles normal
@ -57,8 +54,5 @@ minFaceWeight 0.02;
//- minVolRatio (0 -> 1) //- minVolRatio (0 -> 1)
minVolRatio 0.01; minVolRatio 0.01;
//must be >0 for Fluent compatibility
minTriangleTwist -1;
// ************************************************************************* // // ************************************************************************* //

View File

@ -26,10 +26,10 @@ maxInternalSkewness 4;
// Set to 180 to disable. // Set to 180 to disable.
maxConcave 80; maxConcave 80;
//- Minimum pyramid volume. Is absolute volume of cell pyramid. //- Minimum cell pyramid volume relative to min bounding box length^3
// Set to a sensible fraction of the smallest cell volume expected. // Set to a fraction of the smallest cell volume expected.
// Set to very negative number (e.g. -1e30) to disable. // Set to very negative number (e.g. -1e30) to disable.
minVol 1e-13; minVol 1e-10;
//- Minimum quality of the tet formed by the face-centre //- Minimum quality of the tet formed by the face-centre
// and variable base point minimum decomposition triangles and // and variable base point minimum decomposition triangles and
@ -40,9 +40,6 @@ minVol 1e-13;
// 1 = regular tet // 1 = regular tet
minTetQuality 1e-30; minTetQuality 1e-30;
//- Minimum face area. Set to <0 to disable.
minArea -1;
//- Minimum face twist. Set to <-1 to disable. dot product of face normal //- Minimum face twist. Set to <-1 to disable. dot product of face normal
// and face centre triangles normal // and face centre triangles normal
minTwist 0.05; minTwist 0.05;
@ -57,16 +54,6 @@ minFaceWeight 0.05;
//- minVolRatio (0 -> 1) //- minVolRatio (0 -> 1)
minVolRatio 0.01; minVolRatio 0.01;
//must be >0 for Fluent compatibility
minTriangleTwist -1;
//- If >0 : preserve single cells with all points on the surface if the
// resulting volume after snapping (by approximation) is larger than
// minVolCollapseRatio times old volume (i.e. not collapsed to flat cell).
// If <0 : delete always.
//minVolCollapseRatio 0.5;
// Advanced // Advanced
//- Number of error distribution iterations //- Number of error distribution iterations

View File

@ -26,10 +26,10 @@ maxInternalSkewness 4;
// Set to 180 to disable. // Set to 180 to disable.
maxConcave 80; maxConcave 80;
//- Minimum pyramid volume. Is absolute volume of cell pyramid. //- Minimum cell pyramid volume relative to min bounding box length^3
// Set to a sensible fraction of the smallest cell volume expected. // Set to a fraction of the smallest cell volume expected.
// Set to very negative number (e.g. -1e30) to disable. // Set to very negative number (e.g. -1e30) to disable.
minVol 1e-13; minVol 1e-10;
//- Minimum quality of the tet formed by the face-centre //- Minimum quality of the tet formed by the face-centre
// and variable base point minimum decomposition triangles and // and variable base point minimum decomposition triangles and
@ -40,9 +40,6 @@ minVol 1e-13;
// 1 = regular tet // 1 = regular tet
minTetQuality 1e-30; minTetQuality 1e-30;
//- Minimum face area. Set to <0 to disable.
minArea -1;
//- Minimum face twist. Set to <-1 to disable. dot product of face normal //- Minimum face twist. Set to <-1 to disable. dot product of face normal
// and face centre triangles normal // and face centre triangles normal
minTwist 0.05; minTwist 0.05;
@ -57,16 +54,6 @@ minFaceWeight 0.05;
//- minVolRatio (0 -> 1) //- minVolRatio (0 -> 1)
minVolRatio 0.01; minVolRatio 0.01;
//must be >0 for Fluent compatibility
minTriangleTwist -1;
//- If >0 : preserve single cells with all points on the surface if the
// resulting volume after snapping (by approximation) is larger than
// minVolCollapseRatio times old volume (i.e. not collapsed to flat cell).
// If <0 : delete always.
//minVolCollapseRatio 0.5;
// Advanced // Advanced
//- Number of error distribution iterations //- Number of error distribution iterations
@ -74,8 +61,6 @@ nSmoothScale 4;
//- Amount to scale back displacement at error points //- Amount to scale back displacement at error points
errorReduction 0.75; errorReduction 0.75;
// Optional : some meshing phases allow usage of relaxed rules. // Optional : some meshing phases allow usage of relaxed rules.
// See e.g. addLayersControls::nRelaxedIter. // See e.g. addLayersControls::nRelaxedIter.
relaxed relaxed

View File

@ -26,10 +26,10 @@ maxInternalSkewness 4;
// Set to 180 to disable. // Set to 180 to disable.
maxConcave 80; maxConcave 80;
//- Minimum pyramid volume. Is absolute volume of cell pyramid. //- Minimum cell pyramid volume relative to min bounding box length^3
// Set to a sensible fraction of the smallest cell volume expected. // Set to a fraction of the smallest cell volume expected.
// Set to very negative number (e.g. -1e30) to disable. // Set to very negative number (e.g. -1e30) to disable.
minVol 1e-13; minVol 1e-10;
//- Minimum quality of the tet formed by the face-centre //- Minimum quality of the tet formed by the face-centre
// and variable base point minimum decomposition triangles and // and variable base point minimum decomposition triangles and
@ -40,9 +40,6 @@ minVol 1e-13;
// 1 = regular tet // 1 = regular tet
minTetQuality 1e-30; minTetQuality 1e-30;
//- Minimum face area. Set to <0 to disable.
minArea -1;
//- Minimum face twist. Set to <-1 to disable. dot product of face normal //- Minimum face twist. Set to <-1 to disable. dot product of face normal
// and face centre triangles normal // and face centre triangles normal
minTwist 0.05; minTwist 0.05;
@ -57,16 +54,6 @@ minFaceWeight 0.05;
//- minVolRatio (0 -> 1) //- minVolRatio (0 -> 1)
minVolRatio 0.01; minVolRatio 0.01;
//must be >0 for Fluent compatibility
minTriangleTwist -1;
//- If >0 : preserve single cells with all points on the surface if the
// resulting volume after snapping (by approximation) is larger than
// minVolCollapseRatio times old volume (i.e. not collapsed to flat cell).
// If <0 : delete always.
//minVolCollapseRatio 0.5;
// Advanced // Advanced
//- Number of error distribution iterations //- Number of error distribution iterations

View File

@ -513,10 +513,10 @@ meshQualityControls
// Set to 180 to disable. // Set to 180 to disable.
maxConcave 80; maxConcave 80;
// Minimum pyramid volume. Is absolute volume of cell pyramid. //- Minimum cell pyramid volume relative to min bounding box length^3
// Set to a sensible fraction of the smallest cell volume expected. // Set to a fraction of the smallest cell volume expected.
// Set to very negative number (e.g. -1e30) to disable. // Set to very negative number (e.g. -1e30) to disable.
minVol 1e-13; minVol 1e-10;
// Minimum quality of the tet formed by the face-centre // Minimum quality of the tet formed by the face-centre
// and variable base point minimum decomposition triangles and // and variable base point minimum decomposition triangles and
@ -528,9 +528,6 @@ meshQualityControls
// 1 = regular tet // 1 = regular tet
minTetQuality -1e-30; minTetQuality -1e-30;
// Minimum face area. Set to <0 to disable.
minArea -1;
// Minimum face twist. Set to <-1 to disable. dot product of face normal // Minimum face twist. Set to <-1 to disable. dot product of face normal
// and face centre triangles normal // and face centre triangles normal
minTwist 0.05; minTwist 0.05;
@ -545,19 +542,11 @@ meshQualityControls
// minVolRatio (0 -> 1) // minVolRatio (0 -> 1)
minVolRatio 0.01; minVolRatio 0.01;
// must be >0 for Fluent compatibility
minTriangleTwist -1;
//- If >0 : preserve single cells with all points on the surface if the
// resulting volume after snapping (by approximation) is larger than
// minVolCollapseRatio times old volume (i.e. not collapsed to flat cell).
// If <0 : delete always.
// minVolCollapseRatio 0.5;
// Advanced // Advanced
// Number of error distribution iterations // Number of error distribution iterations
nSmoothScale 4; nSmoothScale 4;
// amount to scale back displacement at error points // amount to scale back displacement at error points
errorReduction 0.75; errorReduction 0.75;

View File

@ -26,10 +26,10 @@ maxInternalSkewness 4;
// Set to 180 to disable. // Set to 180 to disable.
maxConcave 80; maxConcave 80;
//- Minimum pyramid volume. Is absolute volume of cell pyramid. //- Minimum cell pyramid volume relative to min bounding box length^3
// Set to a sensible fraction of the smallest cell volume expected. // Set to a fraction of the smallest cell volume expected.
// Set to very negative number (e.g. -1e30) to disable. // Set to very negative number (e.g. -1e30) to disable.
minVol 1e-13; minVol 1e-10;
//- Minimum quality of the tet formed by the face-centre //- Minimum quality of the tet formed by the face-centre
// and variable base point minimum decomposition triangles and // and variable base point minimum decomposition triangles and
@ -40,9 +40,6 @@ minVol 1e-13;
// 1 = regular tet // 1 = regular tet
minTetQuality 1e-30; minTetQuality 1e-30;
//- Minimum face area. Set to <0 to disable.
minArea -1;
//- Minimum face twist. Set to <-1 to disable. dot product of face normal //- Minimum face twist. Set to <-1 to disable. dot product of face normal
// and face centre triangles normal // and face centre triangles normal
minTwist 0.05; minTwist 0.05;
@ -57,25 +54,14 @@ minFaceWeight 0.05;
//- minVolRatio (0 -> 1) //- minVolRatio (0 -> 1)
minVolRatio 0.01; minVolRatio 0.01;
//must be >0 for Fluent compatibility
minTriangleTwist -1;
//- If >0 : preserve single cells with all points on the surface if the
// resulting volume after snapping (by approximation) is larger than
// minVolCollapseRatio times old volume (i.e. not collapsed to flat cell).
// If <0 : delete always.
//minVolCollapseRatio 0.5;
// Advanced // Advanced
//- Number of error distribution iterations //- Number of error distribution iterations
nSmoothScale 4; nSmoothScale 4;
//- Amount to scale back displacement at error points //- Amount to scale back displacement at error points
errorReduction 0.75; errorReduction 0.75;
// Optional : some meshing phases allow usage of relaxed rules. // Optional : some meshing phases allow usage of relaxed rules.
// See e.g. addLayersControls::nRelaxedIter. // See e.g. addLayersControls::nRelaxedIter.
relaxed relaxed

View File

@ -388,10 +388,10 @@ meshQualityControls
// Set to 180 to disable. // Set to 180 to disable.
maxConcave 80; maxConcave 80;
//- Minimum pyramid volume. Is absolute volume of cell pyramid. //- Minimum cell pyramid volume relative to min bounding box length^3
// Set to a sensible fraction of the smallest cell volume expected. // Set to a fraction of the smallest cell volume expected.
// Set to very negative number (e.g. -1e30) to disable. // Set to very negative number (e.g. -1e30) to disable.
minVol 1e-13; minVol 1e-10;
//- Minimum quality of the tet formed by the face-centre //- Minimum quality of the tet formed by the face-centre
// and variable base point minimum decomposition triangles and // and variable base point minimum decomposition triangles and
@ -403,9 +403,6 @@ meshQualityControls
// 1 = regular tet // 1 = regular tet
minTetQuality -1; // 1e-30; minTetQuality -1; // 1e-30;
//- Minimum face area. Set to <0 to disable.
minArea -1;
//- Minimum face twist. Set to <-1 to disable. dot product of face normal //- Minimum face twist. Set to <-1 to disable. dot product of face normal
// and face centre triangles normal // and face centre triangles normal
minTwist 0.01; minTwist 0.01;
@ -420,14 +417,11 @@ meshQualityControls
//- minVolRatio (0 -> 1) //- minVolRatio (0 -> 1)
minVolRatio 0.01; minVolRatio 0.01;
// must be >0 for Fluent compatibility
minTriangleTwist -1;
// Advanced // Advanced
//- Number of error distribution iterations //- Number of error distribution iterations
nSmoothScale 4; nSmoothScale 4;
//- Amount to scale back displacement at error points //- Amount to scale back displacement at error points
errorReduction 0.75; errorReduction 0.75;

View File

@ -53,32 +53,56 @@ blocks
hex (6 7 11 10 18 19 23 22) (60 24 1) simpleGrading (1 1 1) hex (6 7 11 10 18 19 23 22) (60 24 1) simpleGrading (1 1 1)
); );
patches defaultPatch
{
type empty;
}
boundary
( (
patch inlet inlet
( {
(0 12 16 4) type patch;
(4 16 20 8) faces
) (
patch outlet (0 12 16 4)
( (4 16 20 8)
(7 19 15 3) );
(11 23 19 7) }
)
wall lowerWall outlet
( {
(0 1 13 12) type patch;
(1 5 17 13) faces
(5 6 18 17) (
(2 14 18 6) (7 19 15 3)
(2 3 15 14) (11 23 19 7)
) );
patch atmosphere }
(
(8 20 21 9) lowerWall
(9 21 22 10) {
(10 22 23 11) type wall;
) faces
(
(0 1 13 12)
(1 5 17 13)
(5 6 18 17)
(2 14 18 6)
(2 3 15 14)
);
}
atmosphere
{
type patch;
faces
(
(8 20 21 9)
(9 21 22 10)
(10 22 23 11)
);
}
); );

View File

@ -289,9 +289,10 @@ meshQualityControls
// Set to 180 to disable. // Set to 180 to disable.
maxConcave 80; maxConcave 80;
//- Minimum pyramid volume. Is absolute volume of cell pyramid. //- Minimum cell pyramid volume relative to min bounding box length^3
// Set to a fraction of the smallest cell volume expected.
// Set to very negative number (e.g. -1e30) to disable. // Set to very negative number (e.g. -1e30) to disable.
minVol 1e-20; minVol 1e-10;
//- Minimum quality of the tet formed by the face-centre //- Minimum quality of the tet formed by the face-centre
// and variable base point minimum decomposition triangles and // and variable base point minimum decomposition triangles and
@ -302,9 +303,6 @@ meshQualityControls
// 1 = regular tet // 1 = regular tet
minTetQuality 1e-30; minTetQuality 1e-30;
//- Minimum face area. Set to <0 to disable.
minArea -1;
//- Minimum face twist. Set to <-1 to disable. dot product of face normal //- Minimum face twist. Set to <-1 to disable. dot product of face normal
// and face centre triangles normal // and face centre triangles normal
minTwist 0.02; minTwist 0.02;
@ -319,14 +317,11 @@ meshQualityControls
//- minVolRatio (0 -> 1) //- minVolRatio (0 -> 1)
minVolRatio 0.01; minVolRatio 0.01;
// must be >0 for Fluent compatibility
minTriangleTwist -1;
// Advanced // Advanced
//- Number of error distribution iterations //- Number of error distribution iterations
nSmoothScale 4; nSmoothScale 4;
//- Amount to scale back displacement at error points //- Amount to scale back displacement at error points
errorReduction 0.75; errorReduction 0.75;
} }