From b5b2fd70b65925233d1dfa514cba867f4cdc73e8 Mon Sep 17 00:00:00 2001 From: Timo Koch Date: Wed, 25 Apr 2018 17:26:25 +0200 Subject: [PATCH] [fickslaw][network] Consider that transmissibilities can be zero --- dumux/discretization/cellcentered/tpfa/fickslaw.hh | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/dumux/discretization/cellcentered/tpfa/fickslaw.hh b/dumux/discretization/cellcentered/tpfa/fickslaw.hh index dd8bd86300..e21792b7a5 100644 --- a/dumux/discretization/cellcentered/tpfa/fickslaw.hh +++ b/dumux/discretization/cellcentered/tpfa/fickslaw.hh @@ -242,7 +242,8 @@ private: sumTi += outsideTi; sumXTi += outsideTi*outsideVolVars.moleFraction(phaseIdx, compIdx); } - return sumXTi/sumTi; + + return sumTi > 0 ? sumXTi/sumTi : 0; } //! compute the density at branching facets for network grids as arithmetic mean -- GitLab