Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Include porosity in the pore volume boundary layer #45

Merged
merged 3 commits into from
Apr 2, 2024
+15 −11
Merged

Include porosity in the pore volume boundary layer #45

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
4d50dbf
Include porosity in the porevolume boundary layer
Mar 15, 2024
efb22be
fix code formating
Mar 15, 2024
a64699a
restrict timesteps to make the test pass
Apr 2, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
24 changes: 14 additions & 10 deletions src/pyopmspe11/utils/mapproperties.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -173,7 +173,9 @@ def structured_handling_spe11bc(dic):
dic["satnum"].append(dic["ids_gmsh"][fgl][0])
dic = boxes(dic, dic["xmx_center"][i], z_c, i, dic["satnum"][-1])
dic["permx"].append(dic["rock"][int(dic["ids_gmsh"][fgl][0]) - 1][0])
dic["poro"].append(dic["rock"][int(dic["ids_gmsh"][fgl][0]) - 1][1])
poro = dic["rock"][int(dic["ids_gmsh"][fgl][0]) - 1][1]
dic["poro"].append(poro)
pv = float(poro) * (dic["pvAdded"] + dic["widthBuffer"])
dic["thconr"].append(
f"{dic['rockCond'][int(dic['ids_gmsh'][fgl][0])-1][0]}"
)
Expand All @@ -184,13 +186,13 @@ def structured_handling_spe11bc(dic):
int(dic["ids_gmsh"][fgl][0]) != 1 and int(dic["ids_gmsh"][fgl][0]) != 7
):
dic["porv"].append(
f"PORV {dic['pvAdded']*dic['dy'][0]*dic['dz'][k]} 1 1 1 1 {k+1} {k+1} /"
f"PORV {pv*dic['dy'][0]*dic['dz'][k]} 1 1 1 1 {k+1} {k+1} /"
)
elif i == dic["noCells"][0] - 1 and (
int(dic["ids_gmsh"][fgl][0]) != 1 and int(dic["ids_gmsh"][fgl][0]) != 7
):
dic["porv"].append(
f"PORV {dic['pvAdded']*dic['dy'][0]*dic['dz'][k]} {dic['noCells'][0]} "
f"PORV {pv*dic['dy'][0]*dic['dz'][k]} {dic['noCells'][0]} "
+ f"{dic['noCells'][0]} 1 1 {k+1} {k+1} /"
)
centers.append(
Expand Down Expand Up @@ -235,15 +237,15 @@ def structured_handling_spe11bc(dic):
and int(dic["satnum"][-dic["noCells"][0] + i_i]) != 7
):
dic["porv"].append(
f"PORV {dic['pvAdded']*dic['dy'][j+1]*dic['dz'][k]} 1 1 "
f"PORV {pv*dic['dy'][j+1]*dic['dz'][k]} 1 1 "
+ f"{j+2} {j+2} {k+1} {k+1} /"
)
elif i_i == dic["noCells"][0] - 1 and (
int(dic["satnum"][-dic["noCells"][0] + i_i]) != 1
and int(dic["satnum"][-dic["noCells"][0] + i_i]) != 7
):
dic["porv"].append(
f"PORV {dic['pvAdded']*dic['dy'][j+1]*dic['dz'][k]} {dic['noCells'][0]} "
f"PORV {pv*dic['dy'][j+1]*dic['dz'][k]} {dic['noCells'][0]} "
+ f"{dic['noCells'][0]} {j+2} {j+2} {k+1} {k+1} /"
)
dic["pop1"] = pd.Series(sensor1).argmin()
Expand Down Expand Up @@ -435,21 +437,23 @@ def corner_point_handling_spe11bc(dic):
dic["satnum"].append(dic["ids_gmsh"][fgl][0])
dic = boxes(dic, dic["xyz"][0], z_c, dic["ijk"][0], dic["satnum"][-1])
dic["permx"].append(dic["rock"][int(dic["ids_gmsh"][fgl][0]) - 1][0])
dic["poro"].append(dic["rock"][int(dic["ids_gmsh"][fgl][0]) - 1][1])
poro = dic["rock"][int(dic["ids_gmsh"][fgl][0]) - 1][1]
dic["poro"].append(poro)
pv = float(poro) * (dic["pvAdded"] + dic["widthBuffer"])
dic["disperc"].append(f"{dic['dispersion'][int(dic['ids_gmsh'][fgl][0])-1]}")
dic["thconr"].append(f"{dic['rockCond'][int(dic['ids_gmsh'][fgl][0])-1][0]}")
if dic["ijk"][0] == 0 and (
int(dic["ids_gmsh"][fgl][0]) != 1 and int(dic["ids_gmsh"][fgl][0]) != 7
):
dic["porv"].append(
f"PORV {dic['pvAdded']*dic['d_y'][0]*dic['d_z'][i]} 1 1 1 1 "
f"PORV { pv*dic['d_y'][0]*dic['d_z'][i]} 1 1 1 1 "
+ f"{dic['ijk'][2]+1} {dic['ijk'][2]+1} /"
)
elif dic["ijk"][0] == dic["noCells"][0] - 1 and (
int(dic["ids_gmsh"][fgl][0]) != 1 and int(dic["ids_gmsh"][fgl][0]) != 7
):
dic["porv"].append(
f"PORV {dic['pvAdded']*dic['d_y'][0]*dic['d_z'][i]} {dic['noCells'][0]} "
f"PORV {pv*dic['d_y'][0]*dic['d_z'][i]} {dic['noCells'][0]} "
+ f"{dic['noCells'][0]} 1 1 {dic['ijk'][2]+1} {dic['ijk'][2]+1} /"
)
centers.append(str([dic["xyz"][0], dic["ymy_center"][0], dic["xyz"][2]])[1:-1])
Expand All @@ -476,15 +480,15 @@ def corner_point_handling_spe11bc(dic):
dic["d_zl"] = dic["d_z"][-dic["noCells"][0] + 1 + i]
dic["porv"].append(
"PORV "
+ f"{dic['pvAdded']*dic['d_y'][j+1]*dic['d_zl']} 1 1 "
+ f"{pv*dic['d_y'][j+1]*dic['d_zl']} 1 1 "
+ f"{j+2} {j+2} {dic['ijk'][2]+1} {dic['ijk'][2]+1} /"
)
elif i_i == dic["noCells"][0] - 1 and (
int(dic["satnum"][-dic["noCells"][0] + i_i]) != 1
and int(dic["satnum"][-dic["noCells"][0] + i_i]) != 7
):
dic["porv"].append(
f"PORV {dic['pvAdded']*dic['d_y'][j+1]*dic['d_z'][i]} "
f"PORV {pv*dic['d_y'][j+1]*dic['d_z'][i]} "
+ f"{dic['noCells'][0]} {dic['noCells'][0]} {j+2} {j+2} "
+ f"{dic['ijk'][2]+1} {dic['ijk'][2]+1} /"
)
Expand Down
2 changes: 1 addition & 1 deletion tests/configs/input.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -49,4 +49,4 @@ PERM7 1e-5 PORO7 1e-6 DISP7 0 THCONR7 2.00

"""Define the injection values ([hours] for spe11a; [years] for spe11b/c)"""
"""injection time, time step size to write results, maximum solver time step, injected fluid (0 water, 1 co2) (well1), injection rate [kg/s] (well1), temperature [C] (well1), injected fluid (0 water, 1 co2) (well2), ..."""
25 5 1 1 0.035 10 1 0 10
25 5 .1 1 0.035 10 1 0 10
Loading