recordMolz1983

Parameter values are taken from Molz et al 1983 and Buscheck et al 1983

Extends from Records.Subsurface_Basic (basic subsurface parameters with default values).

Information

1. Purpose of model

Subsurface parameter record with values taken from the well doublet field experiment described in Molz et al. (1983). Intended as a reference parameter set for validation studies of high-temperature aquifer thermal energy storage (HT-ATES) models.

2. Level of detail, physical effects considered, and physical insight

(Purely technical component without physical modeling.)

3. Limits of validity

- Values represent a shallow sand aquifer at the Auburn University test site (Alabama, USA).

- Initial aquifer temperature is 293.15 K (20°C).

4. Interfaces

(no elements)

5. Nomenclature

\lambda_l lambda_l fluid thermal conductivity = 0.61 W/(m·K)
\rho_s rho_s solid density = 2600 kg/m³
c_{p,s} Cp_s specific heat capacity of solid = 696.15 J/(kg·K)
\lambda_s lambda_s solid thermal conductivity = 2.85 W/(m·K)
n n aquifer porosity = 0.25 [-]
k_v k_v vertical permeability = 9×10-12 m²
k k radial permeability = 6.3×10-11 m²
T_\mathrm{init} T_initial initial aquifer temperature = 293.15 K (20°C)

6. Governing Equations

(no equations)

7. Remarks for Usage

Use this record in conjunction with Val_homogenLayer to reproduce the Molz et al. (1983) well doublet experiment.

8. Validation

Tested in check model "ATES.Validierung.Buoyancy.Val_homogenLayer"

9. References

Molz, F. J., Melville, J. G., Parr, A. D., King, D. A., & Hopf, M. T. (1983): Aquifer thermal energy storage: a well doublet experiment at increased temperatures. Water Resources Research, 19(1), 149–160.

10. Version History

Model created by Markus Gillner (markus.gillner@tuhh.de), May 2026

Parameters

TypeNameDefaultDescription
Property values › fluid properties
SI.Densityrho_l (from Subsurface_Basic)988.04Density of fluid, default: density of water at 50°C
SI.DynamicViscosityeta (from Subsurface_Basic)0.466*10^(-3)Dynamic Viscosity of fluid, default: dynamic viscosity of water at 50°C
SI.SpecificHeatCapacityCp_l (from Subsurface_Basic)4181.5Specific Heat Capacity of fluid, default: capacity of water at 50 °C
SI.ThermalConductivitylambda_l (from Subsurface_Basic)0.64Heat Conductivity of the fluid, default: conductivity of water at 50°C
RealC_w (from Subsurface_Basic)4.03*10^6volumetric Heat Capacity of fluid, default: capacity of water at 50°C
Property values › matrix properties in aquifer
SI.Densityrho_s (from Subsurface_Basic)2000Density of solid
SI.SpecificHeatCapacityCp_s (from Subsurface_Basic)900Specific Heat Capacity of solid
SI.ThermalConductivitylambda_s (from Subsurface_Basic)2.56Heat Conductivity of the solid
Realn (from Subsurface_Basic)0.25Porosity of the aquifer
RealC_a (from Subsurface_Basic)1.81*10^6volumetric Heat Capacity of aquifermatrix
SI.Areak_v (from Subsurface_Basic)0.5*10^(-12)Vertical permeability of aquifer, mostly kv=0.1*k
SI.Areak (from Subsurface_Basic)25*10^(-12)Radial permeability of aquifer
Property values › property values of solid in confining layers
Realn_c (from Subsurface_Basic)0.35Porosity of the confining layers
SI.Densityrho_sc (from Subsurface_Basic)2000Density of solid in confing layers
SI.SpecificHeatCapacityCp_sc (from Subsurface_Basic)900Specific Heat Capacity of solid in confining layers
SI.ThermalConductivitylambda_sc (from Subsurface_Basic)2.56Heat Conductivity of the solid in confining layers
Property values › property values of whole aquifer and confining layers
SI.Densityrho (from Subsurface_Basic)n*rho_l + (1 - n)*rho_sDensity of aquifer
SI.Densityrho_c (from Subsurface_Basic)n_c*rho_l + (1 - n_c)*rho_scDensity of aquitard
SI.SpecificHeatCapacityCp (from Subsurface_Basic)n*Cp_l + (1 - n)*Cp_sSpecific Heat Capacity of aquifer
SI.SpecificHeatCapacityCp_c (from Subsurface_Basic)n_c*Cp_l + (1 - n_c)*Cp_scSpecific Heat Capacity of aquitard
RealC (from Subsurface_Basic)n*C_w + (1 - n)*C_aVolumetric heat capacity of aquifer
SI.ThermalConductivitylambda_c (from Subsurface_Basic)n_c*lambda_l + (1 - n_c)*lambda_scHeat Conductivity of the aquitard
SI.ThermalConductivitylambda_a (from Subsurface_Basic)n*lambda_l + (1 - n)*lambda_sHeat Conductivity of the aquifer
SI.Lengthd_h (from Subsurface_Basic)0.7horizontal dispersion coefficient
SI.Lengthd_v (from Subsurface_Basic)0.05vertical dispersion coefficient
Property values › Initial values
SI.Pressurep_initial (from Subsurface_Basic)1000000inital pressure in aquifer at example 100 bar at the lowest point of the well
SI.TemperatureT_initial (from Subsurface_Basic)311.15initial temperature in aquifer at t=0s