From 52d60f0e01bb8633dc4701ecf7bdd3757a60d1a5 Mon Sep 17 00:00:00 2001 From: James E McClure Date: Sun, 28 Oct 2018 10:18:06 -0400 Subject: [PATCH] force tuning for auto relperm --- models/ColorModel.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/models/ColorModel.cpp b/models/ColorModel.cpp index c4011ddd..51fac042 100644 --- a/models/ColorModel.cpp +++ b/models/ColorModel.cpp @@ -552,7 +552,7 @@ void ScaLBL_ColorModel::Run(){ if (lead_timesteps > 5000){ Fz *= capillary_number / Ca; if (Fz > 1e-3) Fz = 1e-3; // impose ceiling for stability - if (Fz < 1e-7) Fz = 1e-6; // impose floor so we don't do something super dumb + if (Fz < 1e-7) Fz = 1e-7; // impose floor so we don't do something super dumb if (Fz*vA_z < 0.f || Fz*vB_z < 0.f) Fz *= 2.f; // bigger forces are needed if flow rate is "too small" tolerance = fabs(capillary_number - Ca) / capillary_number ;