Coupling a Borehole Thermal Model and MT3DMS to Simulate Dynamic Ground Source Heat Pump Efficiency

We present a novel integrated two-region model that couples simulation of local heat transfer processes in borehole heat exchangers (BHEs) with field-scale heat transport simulation using MT3DMS to fulfill the dynamic simulation of the borehole geothermal systems. This includes the prediction of subsurface thermal perturbation induced by BHEs, derivation of the U-pipe circulating fluid temperature profile within boreholes, and the evaluation of the ground source heat pump efficiency based on available time series of building heat load. In our approach, MT3DMS is the simulator for the two-dimensional field-scale heat transport, while new Python modules are developed to analytically solve the thermal transfer process within boreholes and interface iteratively with MT3DMS. A Python package for scripting MODFLOW-based code named Flopy is used to establish the MT3DMS numerical model. The proposed model is validated against analytical solutions and we demonstrate the application to more complex test problems with field-scale heterogeneity and pumping. Instructions are provided to access the source code and example problems which are available online.

Distributed Thermal Response Multi-Source Modeling to Evaluate Heterogeneous Subsurface Properties

A thorough assessment of thermal properties in heterogeneous subsurface is necessary in design of low-temperature borehole heat exchangers (BHEs). A distributed thermal response test (DTRT), which combines distributed temperature sensing (DTS) with a conventional thermal response test (TRT), was conducted in a U-bend geothermal loop installed in an open borehole at the University of Illinois at Urbana-Champaign to estimate thermal properties by analyzing the thermal response of different geologic materials while applying a constant heat input rate. Fiber-optic cables in the DTRT were deployed both inside the U-bend geothermal loop and in the center of the borehole to improve the accuracy of calculated heat-loss rates and borehole temperature profile measurements. To assess the subsurface thermal conductivity during the heating phase of the DTRT, a single-source model and a multi-source model, both based on the infinite line source method, were developed using the borehole temperature data and temperatures inside and along the outside of the loop, separately. The two models returned similar thermal conductivity values. The multi-source modeling has the advantage of predicting the thermal conductivity of heterogeneous geologic materials from borehole temperature profiles during the DTRT heating phase. In addition, based on the distributed thermal conductivity measured in the borehole, estimates were made for both radial thermal impacts and the rate of heat loss in the BHE.

Direct Conservative Domain in the Continuous Galerkin Method for Groundwater Models

The continuous Galerkin finite element method is commonly considered locally nonconservative because a single element with fluxes computed directly from its potential distribution is unable to conserve its mass and fluxes across edges that are discontinuous. Some literature sources have demonstrated that the continuous Galerkin method can be locally conservative with postprocessed fluxes. This paper proposes the concept of a direct conservative domain (DCD), which could conserve mass when fluxes are computed directly from the potential distribution. Also presented here is a method for modifying the advection fluxes to obtain different conservative domains from the DCDs. Furthermore, DCDs are used to analyze the local conservation of several postprocessing algorithms, for which DCDs provide the theoretical basis. The local conservation of DCDs and the proposed method are illustrated and verified by using a hypothetical 2-D model.

The effect of groundwater allocation on economic welfare loss

Water scarcity has become a constraint for regional economic development in many cities and regions. Water rationing serves as one instrument to constrain water consumption to persuade users to save water and to moderate their consumption. When the supply of water is unable to satisfy demand, a loss of welfare for the water users will usually occur. This paper conducts an empirical case study on a Chicago suburban county, McHenry County, to evaluate effective water allocation strategies under possible water scarcity scenarios, by specifically taking into consideration of the economic welfare loss under water rationing. It points out the inefficiency of equal rationing and tests a more effective optimal rationing regime which could significantly lower the overall welfare loss for McHenry County. Instead of a conventional watershed-based approach that would provide little advantage for an area that mostly relies on groundwater, this study adopts regional planning/political boundaries as its spatial analytical units. The outcomes suggest that municipality-level water resources management models, powered under economic welfare objective functions, are both possible and practical. The planning strategy drawn under such optimization models suggests a variety of promising approaches to manage groundwater resources at county scales.