thermal-hydrological-mechanical-modelling-of-stockton-university-reservoir-cooling-system-large

The National Renewable Energy Laboratory in the US specializes in the research and development of renewable energy, energy efficiency, energy systems integration, and sustainable transportation.

national-renewable-energy-laboratory-nrel

National Renewable Energy Laboratory (NREL)

Program TReactMech mech_3D.inp input file. The file includes keyword format descriptions for inputs.

TReactMech_mech_3D_inp_format.txt

Filename "flow.stockton.8lay.6years". Program TReactMech flow.inp input file (a text file).
Global and local flow parameters.
Local parameters organized by (MESH file) rock type.
Formats for individual items described in
TOUGH2 User's Guide, Version 2, appendix E.
All quantities in S.I. (MKS) units.

flow.stockton.8lay.6years.inp.inp

(Zipped) TReactMech MESH file (a text file).
Grid information in TOUGH2 MESH file format.
Lines after "ELEME" have; element name (characters 1-5),
rock type (char. 16-20), volume (char 21-30),
center x, y, z locations, (char 51-60,61-70,71-80).
Lines after "CONNE" have connection information.
Lines after "+++" have connection element number pairs (reduplicatively)  (Fortran (15i6) format  ).
Elements are in x,z,y order; x varying fastest,
z next, and y last; i.e., in successive vertical
slices normal to the y direction.
Formats described in TOUGH2 User's Guide, appendix E.
Coordinates are approximately x = distance NE from cold well 1, y=  distance NW from cold well 1, z = elevation (topography neglected).
Volume stated for constant P T boundary condition cells is
10^50 times actual volume (e.g., see TOUGH2 Users Guide).

MESH.stockton.constPTBC.psvado.cgrav.zip

Program TReactMech mech_3D.inp input file, (a text file).
Global and local mechanics parameters.
Local parameters organized by (MESH file) rock type.
Structure described in file TReactMech_mech_3D_inp_format.txt
All quantities in S.I. (MKS) units.

mech_3D.stockton.7lay.notensclay.updated.inp

(Zipped) TReactMech INCON file (a text file) (TOUGH2 compatible).
Initial pressure and temperature by element name.
File alternates lines of element name (with porosity and perm_xx, perm_yy, perm_zz),
and pressure (Pa) and temperature (deg. C).

SAVE.1e6yrs.zip

TReactMech GENER file (a text file).
Fluid source and sink terms in TOUGH2 GENER file format.
Sources and sinks are separate entries.
For each source after line with element number;
table of times (s), table of injection rates (kg/s),
table of specific enthalpies (J/kg).
For each sink after line with element number;
table of times (s), table of (negative) injection rates (kg/s).
Simulation injection/widthdrawl schedule time 0
at beginning of a November.
Format described in TOUGH2 User's Guide, appendix E.


GENER.multiW.monthlywspike

(zipped) TReactMech initial stress file (a text file).
Format described in file TReactMech_mech_3D_inp_format.txt.
Units are Pa.

initialstress_0.5029sand_0.7003clay_1.zip

(Zipped) TReactMech flow data output file
(at time beginning of February, sixth year).
(a tecplot file, text file.)
Element center x,y,z, pressure, temperature, gas saturation,
liquid saturation, other, other, assumed capillary pressure,
gas density, liquid density, porosity, xx permeability,
yy permeability, zz permeability, vertical permeability/
starting vertical permeability, assumed gas relative permeability,
assumed liquid relative permeability, specific heat capacity,
thermal conductivity, bulk density, diffusivity.
Units are m, m, m, Pa, deg. C, none, none, none, none, Pa,
kg/m^3, kg/m^3, none, m^2, m^2, m^2, none, none, none,
J/m/deg. C, kg/m^3, m^2 /s.
Points are in x,y,z order, i.e., by horizontal planes,
x varying fastest.


flowdata068.tec.zip

(Zipped) TReactMech flowvector output file
(at time beginning of February, sixth year).
(a zipped tecplot file, text file.)
Variables are
element center x,y,z,
liquid flux magnitude, x liquid flux, y liquid flux, z liquid flux,
liquid pore velocity, liquid x pore velocity, liquid y pore velocity, liquid z pore velocity,
gas flux magnitude, x gas flux, y gas flux, z gas flux,
gas pore velocity, gas x pore velocity, gas y pore velocity, gas z pore velocity,
heat flux magnitude, x heat flux, y heat flux, z heat flux.
Units are m, m, m,
kg/m^2 /s, kg/m^2 /s, kg/m^2 /s, kg/m^2 /s,
m/s, m/s, m/s, m/s,
kg/m^2 /s, kg/m^2 /s, kg/m^2 /s, kg/m^2 /s,
kg/s, kg/s, kg/s, kg/s,
J /s, J /s, J /s, J /s.
Points are in x,y,z order, i.e., by horizontal planes,
x varying fastest.


flowvector.068.tec.zip

(Zipped) TReactMech mechanics displacement output file
(at time beginning of February, sixth year).
(a tecplot file, text file.)
element corner x,y,z (m),
x displacement (m), y displacement (m), z displacement (m).
Output points are at element corners.
Points are in x,y,z order, i.e., by horizontal planes,
x varying fastest.

displacement.068.tec.zip

(Zipped) TReactMech mechanics stress and strain output file
(at time beginning of February, sixth year).
(a tecplot file, text file.)
"X(m)",    "Y(m)",     "Z(m)",
"Sigma_xx",   "Sigma_yy",   "Sigma_zz",
"Sigma_yz",   "Sigma_xz",   "Sigma_xy",
"Strain_xx",  "Strain_yy",  "Strain_zz",
"Strain_yz",  "Strain_xz",  "Strain_xy",
"E_fail_xx",  "E_fail_yy",  "E_fail_zz",
"2*E_fail_yz","2*E_fail_xz","2*E_fail_xy","E_fail_vol",
"Poro_elast",  "Press".
(center x,y,z position (m), total stress (Pa),
total strain (non-dimen.), failure strain (non-dimen.),
poroelastic porosity (non-dimen.), fluid pressure (Pa).)
Output points are at element centers,
stresses and strains are averaged over Gauss points.
Points are in x,y,z order, i.e., by horizontal planes.

stress_strain.068.tec.zip

TReactMech mechanical failure output file, at all times.
(a text file.)
At each time, a header, then
X(m),  Y(m),  Z(m), ielem, N_shear, N_tensl,  N_Surf,
E_fail_xx, E_fail_yy, E_fail_zz,
2*E_fail_yz, 2*E_fail_xz, 2*E_fail_xy, E_fail_vol,
at each mechanically failed (fractured) element.
Output points are at element centers.
Variables are;
center x,y,z position (m), element number, maximum number of independent
planes of shear failure, maximum number of independent planes of
tensile failure, maximum number of independent failure planes,
failure strain, volumetric failure strain.
Failure strains are averaged over the 8 Gauss points of an element.
Maximum numbers are over the 8 Gauss points.
Only elements that have failed up to the reported time step are reported.


Plasticity.tec

J.T. Smith, E. Sonnenthal, P. Dobson, P. Nico,
and M. Worthington, 2021.
Thermal-hydrological-mechanical modeling of Stockton University reservoir cooling system,
Proceedings of the 46th Workshop on Geothermal Reservoir Engineering,
Stanford University, SGP-TR-218.

T-H-M Modeling of Stockton University Reservoir Cooling System

K. Pruess, C. Oldenburg, and G. Moridis, 2012.
TOUGH2 User's Guide, Version 2,
Lawrence Berkeley National Laboratory, LBNL-43134.

TOUGH2 User's Guide Version 2

An overview of the files included in this submission.