fix relperm bug
This commit is contained in:
parent
c4d1268278
commit
b53b7984d6
@ -804,14 +804,14 @@ void ScaLBL_ColorModel::Run(){
|
||||
|
||||
double flow_rate_A_connected = Vol_nc*(vAc_x*dir_x + vAc_y*dir_y + vAc_z*dir_z);
|
||||
double flow_rate_B_connected = Vol_wc*(vBc_x*dir_x + vBc_y*dir_y + vBc_z*dir_z);
|
||||
double flow_rate_A_disconnected = Vol_nd*(vAd_x*dir_x + vAd_y*dir_y + vAd_z*dir_z);
|
||||
double flow_rate_B_disconnected = Vol_wd*(vBd_x*dir_x + vBd_y*dir_y + vBd_z*dir_z);
|
||||
double flow_rate_A_total = (Vol_nd+Vol_nc)*(vAd_x*dir_x + vAd_y*dir_y + vAd_z*dir_z + vAc_x*dir_x + vAc_y*dir_y + vAc_z*dir_z);
|
||||
double flow_rate_B_total = (Vol_wd+Vol_wc)*(vBd_x*dir_x + vBd_y*dir_y + vBd_z*dir_z + vBc_x*dir_x + vBc_y*dir_y + vBc_z*dir_z);
|
||||
|
||||
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 = h*h*muA*(flow_rate_A_connected+flow_rate_A_disconnected)/(force_mag);
|
||||
double kBeff = h*h*muB*(flow_rate_B_connected+flow_rate_B_disconnected)/(force_mag);
|
||||
double kAeff = h*h*muA*(flow_rate_A_total)/(force_mag);
|
||||
double kBeff = h*h*muB*(flow_rate_B_total)/(force_mag);
|
||||
double viscous_pressure_drop = (rhoA*volA + rhoB*volB)*force_mag;
|
||||
double Mobility = muA/muB;
|
||||
|
||||
|
Loading…
Reference in New Issue
Block a user