From 7394e7ff331c36b86ab1c1cd7677b8a05b82e7f2 Mon Sep 17 00:00:00 2001 From: James McClure Date: Tue, 29 Aug 2023 19:34:23 -0400 Subject: [PATCH] fix bug in capillary flux computation --- models/ColorModel.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/models/ColorModel.cpp b/models/ColorModel.cpp index ae3c9fe7..88caa2ea 100644 --- a/models/ColorModel.cpp +++ b/models/ColorModel.cpp @@ -521,7 +521,7 @@ void ScaLBL_ColorModel::Initialize() { double IFT = 6.0 * alpha; double CrossSectionalArea = (double)(nprocx * (Nx - 2) * nprocy * (Ny - 2)); - flux = Mask->Porosity() * CrossSectionalArea * (Ny - 2) * IFT * + flux = Mask->Porosity() * CrossSectionalArea * IFT * capillary_number / MuB; if (rank == 0) printf(" flux=%f \n", flux);