From 1956b71d0a6c1e86b96992adabafc10662f51b9b Mon Sep 17 00:00:00 2001 From: Bhargav Akula Date: Wed, 8 Apr 2026 08:35:01 -0400 Subject: [PATCH] reveted some other edits so that we can make the recilinear heterogeneous case working again --- README.md | 2 +- lib/simulation.py | 22 +++++++++++----------- lib/test.py | 2 +- lib/water.py | 6 +++--- 4 files changed, 16 insertions(+), 16 deletions(-) diff --git a/README.md b/README.md index 66c3a31..de16f88 100644 --- a/README.md +++ b/README.md @@ -24,7 +24,7 @@ This research aims to view the effects of surfactant in EOR while accounting for - Analyze and visualize results to gain insights into SP flooding dynamics ## Requirements - - ```Python 3.13.1``` (with dependencies listed in ```requirements.txt```) + - ```Python 3.14.3``` (with dependencies listed in ```requirements.txt```) ### Installing Dependencies 1. cd to the location of the repository within your computer diff --git a/lib/simulation.py b/lib/simulation.py index 399ffab..41827d4 100644 --- a/lib/simulation.py +++ b/lib/simulation.py @@ -1475,17 +1475,17 @@ def run(self): ## STEP 2.7: Updating the cummulative oil captured, Production rate, and the residual oil in place # arrays for exporting to CSV files - if (t_cal == 0): - self.COC[0,t_cal] = ocut - else: - self.COC[0,t_cal] = self.COC[0,t_cal - 1] + ocut - - self.ProdRate[0,t_cal] = ocut/dt - self.CROIP[0,t_cal] = ROIP - - t_cal += 1 - - self._export_results() + # if (t_cal == 0): + # self.COC[0,t_cal] = ocut + # else: + # self.COC[0,t_cal] = self.COC[0,t_cal - 1] + ocut + # + # self.ProdRate[0,t_cal] = ocut/dt + # self.CROIP[0,t_cal] = ROIP + # + # t_cal += 1 + # + # self._export_results() except Exception as e: print(e) diff --git a/lib/test.py b/lib/test.py index cfe09e8..09627b5 100644 --- a/lib/test.py +++ b/lib/test.py @@ -50,7 +50,7 @@ user_dict = { "simulation_id": 1, "model_type": MODEL["No_Shear_Thinning"], - "reservoir_geometry": GEOMETRY["Rectilinear"], + "reservoir_geometry": GEOMETRY["Quarter Five Spot"], "permeability": PERMEABILITY["Heterogeneous"], "polymer_type": POLYMER["Xanthane"], "polymer_concentration": 0.001, diff --git a/lib/water.py b/lib/water.py index 825426f..4657b40 100644 --- a/lib/water.py +++ b/lib/water.py @@ -579,7 +579,7 @@ def compute_water_saturation( (Qmod[cnt][i] / dt_array[cnt][i]) + g1 * (1 - f[cnt][i]) + ( - (D_g[cnt][i] + D_g[cnt][i + 1]) / (dx**2) + (D_g[cnt][i] + D_g[cnt+1][i + 1]) / (dx**1) + (D_g[cnt + 1][i] + D_g[cnt][i]) / (dx**2) ) * surfactant.concentration_matrix[cnt][i] @@ -599,7 +599,7 @@ def compute_water_saturation( - (D_s[cnt + 1][i] + D_s[cnt][i]) / (dy**2) ) - BB[j][i + 1] = (D_s[cnt][i] + D_s[cnt][i + 1]) / (dx**2) + BB[j][i + 1] = (D_s[cnt][i] + D_s[cnt][i + 1]) * (dx**2) elif i == m - 1: # last/rightmost column DD[i] = ( Qmod[cnt][i] / dt_array[cnt][i] @@ -1001,7 +1001,7 @@ def compute_water_saturation( BB[j][i] = 1 / dt_array[cnt][i] - ( (1 / (2 * dx**2)) - * (D_s[cnt][i - 1] + 2 * D_s[cnt][i] + D_s[cnt][i + 1]) + * (D_s[cnt][i] + 2 * D_s[cnt][i] + D_s[cnt][i + 1]) + (1 / (2 * dy**2)) * ( D_s[cnt - 1][i]