From d53e3bb2649b77a90cc547bc02825fa1da4b5909 Mon Sep 17 00:00:00 2001 From: JamesEMcclure Date: Tue, 24 Sep 2019 14:49:28 -0400 Subject: [PATCH] fix density in relperm --- models/ColorModel.cpp | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/models/ColorModel.cpp b/models/ColorModel.cpp index a31586fe..44fa8fee 100644 --- a/models/ColorModel.cpp +++ b/models/ColorModel.cpp @@ -812,11 +812,11 @@ void ScaLBL_ColorModel::Run(){ double kAeff_connected = h*h*muA*flow_rate_A_connected/(force_mag); double kBeff_connected = h*h*muB*flow_rate_B_connected/(force_mag); - double kAeff_disconnected = h*h*muA*flow_rate_A_disconnected/(rhoA*force_mag); - double kBeff_disconnected = h*h*muB*flow_rate_B_disconnected/(rhoB*force_mag); + double kAeff_disconnected = h*h*muA*flow_rate_A_disconnected/(force_mag); + double kBeff_disconnected = h*h*muB*flow_rate_B_disconnected/(force_mag); - double kAeff = h*h*muA*(flow_rate_A)/(rhoA*force_mag); - double kBeff = h*h*muB*(flow_rate_B)/(rhoB*force_mag); + double kAeff = h*h*muA*(flow_rate_A)/(force_mag); + double kBeff = h*h*muB*(flow_rate_B)/(force_mag); double viscous_pressure_drop = (rhoA*volA + rhoB*volB)*force_mag; double Mobility = muA/muB;