From 913eba05f89431608b520f49445e4dc9373b3611 Mon Sep 17 00:00:00 2001 From: graham Date: Tue, 20 Oct 2009 13:35:57 +0100 Subject: [PATCH] Adding harmonic (commented) and geometric means of cell sizes to reduce the impact of large step changes in cell size. --- .../conformalVoronoiMesh/conformalVoronoiMesh.C | 13 +++++++++++-- 1 file changed, 11 insertions(+), 2 deletions(-) diff --git a/src/conformalVoronoiMesh/conformalVoronoiMesh/conformalVoronoiMesh.C b/src/conformalVoronoiMesh/conformalVoronoiMesh/conformalVoronoiMesh.C index bf2ef6e50b..e3445e68e9 100644 --- a/src/conformalVoronoiMesh/conformalVoronoiMesh/conformalVoronoiMesh.C +++ b/src/conformalVoronoiMesh/conformalVoronoiMesh/conformalVoronoiMesh.C @@ -2523,8 +2523,18 @@ void Foam::conformalVoronoiMesh::move() > cvMeshControls().cosAlignmentAcceptanceAngle() ) { + // Arithmetic mean + // scalar targetCellSize = + // 0.5*(vA->targetCellSize() + vB->targetCellSize()); + + // Geometric mean scalar targetCellSize = - 0.5*(vA->targetCellSize() + vB->targetCellSize()); + sqrt(vA->targetCellSize()*vB->targetCellSize()); + + // Harmonic mean + // scalar targetCellSize = + // 2.0*(vA->targetCellSize()*vB->targetCellSize()) + // /(vA->targetCellSize() + vB->targetCellSize()); scalar targetFaceArea = sqr(targetCellSize); @@ -2712,4 +2722,3 @@ void Foam::conformalVoronoiMesh::move() // ************************************************************************* // -