From 58f6258d7a006ab07489ce8ccec644c9c339714b Mon Sep 17 00:00:00 2001
From: sergio <sergio>
Date: Fri, 27 Sep 2019 15:16:00 -0700
Subject: [PATCH] BUG: Correcting delta in fa consistent with fv.
Solve issue 1431 Gitlab. delta() calculation on the patch is consistent
with not applying non-orthogonality on the patch.
---
src/finiteArea/faMesh/faPatches/faPatch/faPatch.C | 5 +++--
1 file changed, 3 insertions(+), 2 deletions(-)
diff --git a/src/finiteArea/faMesh/faPatches/faPatch/faPatch.C b/src/finiteArea/faMesh/faPatches/faPatch/faPatch.C
index b0f16927219..199c2bfc5ce 100644
--- a/src/finiteArea/faMesh/faPatches/faPatch/faPatch.C
+++ b/src/finiteArea/faMesh/faPatches/faPatch/faPatch.C
@@ -20,7 +20,7 @@ License
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
for more details.
- You should have received a copy of the GNU General Public License
+// You should have received a copy of the GNU General Public License
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
\*---------------------------------------------------------------------------*/
@@ -424,7 +424,8 @@ Foam::tmp<Foam::vectorField> Foam::faPatch::edgeFaceCentres() const
Foam::tmp<Foam::vectorField> Foam::faPatch::delta() const
{
- return edgeCentres() - edgeFaceCentres();
+ return edgeNormals()*(edgeNormals() & (edgeCentres() - edgeFaceCentres()));
+ //return edgeCentres() - edgeFaceCentres();
}
--
GitLab