Description
Mesh, properties, initial conditions, injection/withdrawal rates for modeling thermal, hydrological, and mechanical effects of fluid injection to and withdrawal from ground for Stockton University reservoir cooling system (aquifer storage cooling system), Galloway, New Jersey, on large scale grid, with some results. First simulation of 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, from which Figures 1-5 pertain.
Resources
| Name | Format | Description | Link |
|---|---|---|---|
| displacement.068.tec.zip | 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. | Open resource |
| File Overview.txt | Text | An overview of the files included in this submission. | Open resource |
| flowdata068.tec.zip | 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. | Open resource |
| flow.stockton.8lay.6years.inp.inp | Unknown | 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. | Open resource |
| flowvector.068.tec.zip | 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. | Open resource |
| GENER.multiW.monthlywspike | Unknown | 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. | Open resource |
| initialstress_0.5029sand_0.7003clay_1.zip | ZIP | (zipped) TReactMech initial stress file (a text file). Format described in file TReactMech_mech_3D_inp_format.txt. Units are Pa. | Open resource |
| mech_3D.stockton.7lay.notensclay.updated.inp | Unknown | 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. | Open resource |
| MESH.stockton.constPTBC.psvado.cgrav.zip | ZIP | (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). | Open resource |
| Plasticity.tec | Unknown | 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. | Open resource |
| SAVE.1e6yrs.zip | ZIP | (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). | Open resource |
| stress_strain.068.tec.zip | 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. | Open resource |
| TReactMech_mech_3D_inp_format.txt | Text | Program TReactMech mech_3D.inp input file. The file includes keyword format descriptions for inputs. | Open resource |
| T-H-M Modeling of Stockton University Reservoir Cooling System | HTML | 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. | Open resource |
| TOUGH2 User's Guide Version 2 | K. Pruess, C. Oldenburg, and G. Moridis, 2012. TOUGH2 User's Guide, Version 2, Lawrence Berkeley National Laboratory, LBNL-43134. | Open resource |