zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

ENH: Adding deltaCoeff comparison code - requires fvMesh.

Browse Source
This commit is contained in:
graham
2010-12-20 11:45:35 +00:00
parent a8727486a4
commit 53529c1112
2 changed files with 26 additions and 15 deletions
Download Patch File Download Diff File Expand all files Collapse all files

1
src/mesh/conformalVoronoiMesh/conformalVoronoiMesh/conformalVoronoiMesh.H
View File

@ -59,6 +59,7 @@ SourceFiles
#include "unitConversion.H"
#include "transform.H"
#include "volFields.H"
#include "surfaceFields.H"
#include "fvMesh.H"
#include "Histogram.H"
#include "labelPair.H"

40
src/mesh/conformalVoronoiMesh/conformalVoronoiMesh/conformalVoronoiMeshCalcDualMesh.C
View File

@ -1405,11 +1405,11 @@ Foam::label Foam::conformalVoronoiMesh::checkPolyMeshQuality
removeUnusedCells(owner, neighbour, cellCentres);
polyMesh pMesh
fvMesh mesh
(
IOobject
(
"cvMesh_temporary",
Foam::polyMesh::defaultRegion,
runTime_.timeName(),
runTime_,
IOobject::NO_READ,
@ -1431,23 +1431,33 @@ Foam::label Foam::conformalVoronoiMesh::checkPolyMeshQuality
patchSizes[p],
patchStarts[p],
p,
pMesh.boundaryMesh()
mesh.boundaryMesh()
);
}
pMesh.addPatches(patches);
mesh.addPatches(patches);
pMesh.overrideCellCentres(cellCentres);
surfaceScalarField preCellCentreShiftDeltaCoeffs = mesh.deltaCoeffs();
mesh.clearOut();
surfaceScalarField postCellCentreShiftDeltaCoeffs = mesh.deltaCoeffs();
surfaceScalarField deltaCoeffRatio =
postCellCentreShiftDeltaCoeffs/preCellCentreShiftDeltaCoeffs;
Info<< "deltaCoeffRatio min " << min(deltaCoeffRatio) << nl
<< "deltaCoeffRatio max " << max(deltaCoeffRatio) << endl;
timeCheck();
labelHashSet wrongFaces(pMesh.nFaces()/100);
labelHashSet wrongFaces(mesh.nFaces()/100);
Info << endl;
DynamicList<label> checkFaces(pMesh.nFaces());
DynamicList<label> checkFaces(mesh.nFaces());
const vectorField& fAreas = pMesh.faceAreas();
const vectorField& fAreas = mesh.faceAreas();
scalar faceAreaLimit = SMALL;
@ -1468,7 +1478,7 @@ Foam::label Foam::conformalVoronoiMesh::checkPolyMeshQuality
motionSmoother::checkMesh
(
false,
pMesh,
mesh,
cvMeshControls().cvMeshDict().subDict("meshQualityControls"),
checkFaces,
wrongFaces
@ -1476,7 +1486,7 @@ Foam::label Foam::conformalVoronoiMesh::checkPolyMeshQuality
label nInvalidPolyhedra = 0;
const cellList& cells = pMesh.cells();
const cellList& cells = mesh.cells();
forAll(cells, cI)
{
@ -1502,7 +1512,7 @@ Foam::label Foam::conformalVoronoiMesh::checkPolyMeshQuality
forAllConstIter(labelHashSet, wrongFaces, iter)
{
const face f = pMesh.faces()[iter.key()];
const face f = mesh.faces()[iter.key()];
forAll(f, fPtI)
{
@ -1512,13 +1522,13 @@ Foam::label Foam::conformalVoronoiMesh::checkPolyMeshQuality
// // Limit connected cells
// labelHashSet limitCells(pMesh.nCells()/100);
// labelHashSet limitCells(mesh.nCells()/100);
// const labelListList& ptCells = pMesh.pointCells();
// const labelListList& ptCells = mesh.pointCells();
// forAllConstIter(labelHashSet, wrongFaces, iter)
// {
// const face f = pMesh.faces()[iter.key()];
// const face f = mesh.faces()[iter.key()];
// forAll(f, fPtI)
// {
@ -1533,7 +1543,7 @@ Foam::label Foam::conformalVoronoiMesh::checkPolyMeshQuality
// }
// }
// const labelListList& cellPts = pMesh.cellPoints();
// const labelListList& cellPts = mesh.cellPoints();
// forAllConstIter(labelHashSet, limitCells, iter)
// {