From c064448239bd7e79a388ab94e6f79c5a08744614 Mon Sep 17 00:00:00 2001
From: Hrvoje Jasak <h.jasak@wikki.co.uk>
Date: Wed, 25 Sep 2019 16:42:10 +0100
Subject: [PATCH] Stabilise fvcReconstruct in the presence of zero area faces. 
 Experimental

---
 src/finiteVolume/finiteVolume/fvc/fvcReconstruct.C | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/src/finiteVolume/finiteVolume/fvc/fvcReconstruct.C b/src/finiteVolume/finiteVolume/fvc/fvcReconstruct.C
index 06982ba1e..166789409 100644
--- a/src/finiteVolume/finiteVolume/fvc/fvcReconstruct.C
+++ b/src/finiteVolume/finiteVolume/fvc/fvcReconstruct.C
@@ -120,7 +120,9 @@ reconstruct
     // G  -> G
 
     // Calculate sum of the directional fluxes
-    const surfaceScalarField magSfSqr = sqr(mesh.magSf());
+    const surfaceScalarField magSfSqr =
+        sqr(mesh.magSf() + dimensionedScalar("vsmall", dimArea, 1e-100));
+
     const GeometricField<GradType, fvPatchField, volMesh> fluxTimesNormal =
         surfaceSum((mesh.Sf()/magSfSqr)*ssf);