From 8fe96e5dee9a4a700cb5715dc30c4cd7f0b8db73 Mon Sep 17 00:00:00 2001
From: sergio <sergio>
Date: Mon, 19 Aug 2019 13:38:36 -0700
Subject: [PATCH] ENH: Changing tolerance overlap. Taking out filter by small
overlap
---
.../cellVolumeWeight/cellVolumeWeightCellCellStencil.C | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/src/overset/cellCellStencil/cellVolumeWeight/cellVolumeWeightCellCellStencil.C b/src/overset/cellCellStencil/cellVolumeWeight/cellVolumeWeightCellCellStencil.C
index 4f02c9d466c..60464e71eb5 100644
--- a/src/overset/cellCellStencil/cellVolumeWeight/cellVolumeWeightCellCellStencil.C
+++ b/src/overset/cellCellStencil/cellVolumeWeight/cellVolumeWeightCellCellStencil.C
@@ -49,7 +49,7 @@ namespace cellCellStencils
}
Foam::scalar
-Foam::cellCellStencils::cellVolumeWeight::defaultOverlapTolerance_ = 1e-9;
+Foam::cellCellStencils::cellVolumeWeight::defaultOverlapTolerance_ = 1e-6;
// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
@@ -953,7 +953,7 @@ bool Foam::cellCellStencils::cellVolumeWeight::update()
// Interpolate. Check if enough overlap
scalar v = mesh_.V()[cellI];
scalar overlapVol = sum(allWeights[cellI]);
- if (overlapVol > (1.0-overlapTolerance_)*v)
+ if (overlapVol > overlapTolerance_*v)
{
allCellTypes[cellI] = INTERPOLATED;
}
@@ -971,7 +971,7 @@ bool Foam::cellCellStencils::cellVolumeWeight::update()
}
}
-
+/*
// Knock out cell with insufficient interpolation weights
forAll(allCellTypes, cellI)
{
@@ -989,7 +989,7 @@ bool Foam::cellCellStencils::cellVolumeWeight::update()
}
}
}
-
+*/
if (debug)
{
tmp<volScalarField> tfld
--
GitLab