From 0fc4b79dde0b7d7a4595f4a1d975dade0f89f038 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Mon, 6 Jan 2020 13:59:16 +0000
Subject: [PATCH] ENH: checkMesh: write minPyrVol. Fixes #1543.
---
.../mesh/manipulation/checkMesh/checkMesh.C | 1 +
.../mesh/manipulation/checkMesh/writeFields.C | 63 +++++++++++++++++++
2 files changed, 64 insertions(+)
diff --git a/applications/utilities/mesh/manipulation/checkMesh/checkMesh.C b/applications/utilities/mesh/manipulation/checkMesh/checkMesh.C
index 5c0c64747c8..216e2b1ce47 100644
--- a/applications/utilities/mesh/manipulation/checkMesh/checkMesh.C
+++ b/applications/utilities/mesh/manipulation/checkMesh/checkMesh.C
@@ -161,6 +161,7 @@ int main(int argc, char *argv[])
"cellVolume",
"cellVolumeRatio",
"minTetVolume",
+ "minPyrVolume",
"cellRegion",
"wallDistance"
});
diff --git a/applications/utilities/mesh/manipulation/checkMesh/writeFields.C b/applications/utilities/mesh/manipulation/checkMesh/writeFields.C
index 1623278c6b7..7f86c28bd9b 100644
--- a/applications/utilities/mesh/manipulation/checkMesh/writeFields.C
+++ b/applications/utilities/mesh/manipulation/checkMesh/writeFields.C
@@ -472,6 +472,69 @@ void Foam::writeFields
minTetVolume.write();
}
+ // minPyrVolume
+ if (selectedFields.found("minPyrVolume"))
+ {
+ volScalarField minPyrVolume
+ (
+ IOobject
+ (
+ "minPyrVolume",
+ mesh.time().timeName(),
+ mesh,
+ IOobject::NO_READ,
+ IOobject::AUTO_WRITE,
+ false
+ ),
+ mesh,
+ dimensionedScalar("minPyrVolume", dimless, GREAT),
+ zeroGradientFvPatchScalarField::typeName
+ );
+
+ // Get owner and neighbour pyr volumes
+ scalarField ownPyrVol(mesh.nFaces());
+ scalarField neiPyrVol(mesh.nInternalFaces());
+ primitiveMeshTools::facePyramidVolume
+ (
+ mesh,
+ mesh.points(),
+ mesh.cellCentres(),
+
+ ownPyrVol,
+ neiPyrVol
+ );
+
+ // Get min pyr vol per cell
+ scalarField& cellFld = minPyrVolume.ref();
+ cellFld = GREAT;
+
+ const labelUList& own = mesh.owner();
+ const labelUList& nei = mesh.neighbour();
+
+ // Internal faces
+ forAll(own, facei)
+ {
+ cellFld[own[facei]] = min(cellFld[own[facei]], ownPyrVol[facei]);
+ cellFld[nei[facei]] = min(cellFld[nei[facei]], neiPyrVol[facei]);
+ }
+
+ // Patch faces
+ for (const auto& fvp : minPyrVolume.boundaryField())
+ {
+ const labelUList& fc = fvp.patch().faceCells();
+
+ forAll(fc, i)
+ {
+ const label meshFacei = fvp.patch().start();
+ cellFld[fc[i]] = min(cellFld[fc[i]], ownPyrVol[meshFacei]);
+ }
+ }
+
+ minPyrVolume.correctBoundaryConditions();
+ Info<< " Writing minPyrVolume to " << minPyrVolume.name() << endl;
+ minPyrVolume.write();
+ }
+
if (selectedFields.found("cellRegion"))
{
volScalarField cellRegion
--
GitLab