From 721da9c038b191549c3374afc8bcb863ca5f922f Mon Sep 17 00:00:00 2001 From: Henry Date: Wed, 26 Oct 2011 14:35:34 +0100 Subject: [PATCH] face: add special treatment for the sweptVol of a triangular face --- src/OpenFOAM/meshes/meshShapes/face/face.C | 21 +++++++++++++++++++++ 1 file changed, 21 insertions(+) diff --git a/src/OpenFOAM/meshes/meshShapes/face/face.C b/src/OpenFOAM/meshes/meshShapes/face/face.C index 995c75cdc0..3211df1514 100644 --- a/src/OpenFOAM/meshes/meshShapes/face/face.C +++ b/src/OpenFOAM/meshes/meshShapes/face/face.C @@ -655,6 +655,27 @@ Foam::scalar Foam::face::sweptVol const pointField& newPoints ) const { + if (size() == 3) + { + return + ( + triPointRef + ( + oldPoints[operator[](0)], + oldPoints[operator[](1)], + oldPoints[operator[](2)] + ).sweptVol + ( + triPointRef + ( + newPoints[operator[](0)], + newPoints[operator[](1)], + newPoints[operator[](2)] + ) + ) + ); + } + scalar sv = 0; // Calculate the swept volume by breaking the face into triangles and