Stabilise fvcReconstruct in the presence of zero area faces. Experimental

2019-09-25 16:42:10 +01:00 · 2019-09-25 16:42:10 +01:00 · c064448239
commit c064448239
parent dd0f2bb14e
1 changed files with 3 additions and 1 deletions
--- 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);