Explore the vast range of titles available.

With their high storage capacity and energy efficiency as well as the compatibilities with renewable energy sources, high‐temperature aquifer thermal energy storage (HT‐ATES) systems are frequently the target today in the design of temporally and spatially balanced and continuous energy supply systems. The inherent density‐driven buoyancy flow is of greater importance with HT‐ATES, which may lead to a lower thermal recovery efficiency than conventional low‐temperature ATES. In this study, the governing equations for HT‐ATES considering buoyancy flow are nondimensionalized, and four key dimensionless parameters regarding thermal recovery efficiency are determined. Then, using numerical simulations, recovery efficiency for a sweep of the key dimensionless parameters for multiple cycles and storage volumes is examined. Ranges of the key dimensionless parameters for the three displacement regimes, that is, a buoyancy‐dominated regime, a conduction‐dominated regime, and a transition regime, are identified. In the buoyancy‐dominated regime, recovery efficiency is mainly correlated to the ratio between the Rayleigh number and the Peclet number. In the conduction‐dominated regime, recovery efficiency is mainly correlated to the product of a material‐related parameter and the Peclet number. Multivariable regression functions are provided to estimate recovery efficiency using the dimensionless parameters. The recovery efficiency estimated by the regression function shows good agreement with the simulation results. Additionally, well screen designs for optimizing recovery efficiency at various degrees of intensity of buoyancy flow are investigated. The findings of this study can be used for a quick assessment and characterization of the potential HT‐ATES systems based on the geological and operational parameters. Key Points Four key dimensionless parameters for the high‐temperature aquifer thermal energy storage systems are identified The displacement processes are classified into a buoyancy‐dominated regime, a conduction‐dominated regime, and a transition regime Multivariable regression functions are demonstrated for the estimation of thermal recovery efficiency

A novel reactive smart tracer method, termed the kinetic interface-sensitive (KIS) tracer test, has been demonstrated in laboratory column experiments to enable measurement of the specific capillary-associated fluid–fluid interfacial area in dynamic two-phase flow displacement processes in porous media. Development of the tracer method towards effective application in real field conditions requires investigation of the influence of the porous media heterogeneity on the front size and the specific interfacial area, and, consequently, in how far a kinetic interface-sensitive tracer experiment, and the corresponding breakthrough curves, are affected. This study employs a two-dimensional Darcy-scale two-phase flow reactive transport model to investigate numerically the KIS tracer transport in heterogeneous porous media. Simulations were carried out for the primary drainage process in a domain formed of fine and coarse porous media. Various heterogeneity patterns, having different numbers of inclusions and different geometrical distributions, were studied. It is shown that the shape of the breakthrough curves can be used as an indicator for quantifying the displacement front roughness, the specific interfacial area in the domain, and the domain heterogeneity, e.g., the existence of preferential flow pathways inside the porous media. The results indicate that when the displacement front roughness is small, the concentration breakthrough curves exhibit a linear increase. The slopes of the breakthrough curves linearly depend on the fraction of the bulk volume occupied by the low-permeability sand inclusions. The volume-averaged specific interfacial area and the size of the transition zone can be determined from the slopes of the breakthrough curves.

In this article we consider the Dirichlet problem on a bounded domain Ω ⊂ R d with respect to a second-order elliptic differential operator in divergence form. We do not assume a divergence condition as in the pioneering work by Stampacchia, but merely assume that 0 is not a Dirichlet eigenvalue. The purpose of this article is to define and investigate a solution of the Dirichlet problem, which we call Perron solution, in a setting where no maximum principle is available. We characterise this solution in different ways: by approximating the domain by smooth domains from the interior, by variational properties, by the pointwise boundary behaviour at regular boundary points and by using the approximative trace. We also investigate for which boundary data the Perron solution has finite energy. Finally we show that the Perron solution is obtained as an H 0 1 -perturbation of a continuous function on  Ω ¯ . This is new even for the Laplacian and solves an open problem.

Knowledge of soil moisture dynamics in the unsaturated soil zone provides valuable information on the temporal and spatial variability of groundwater recharge. This is especially true for the Mediterranean region, where a substantial fraction of long-term groundwater recharge is expected to occur during high magnitude precipitation events of above-average wet winters. To elucidate process understanding of infiltration processes during these extreme events, a monitoring network of precipitation gauges, meteorological stations, and soil moisture plots was installed in an area with a steep climatic gradient in the Jordan Valley region. In three soil moisture plots, Hydrus-1D was used to simulate water movement in the unsaturated soil zone with soil hydraulic parameters estimated by the Shuffled Complex Evolution Metropolis algorithm. To generalize our results, we modified soil depth and rainfall input to simulate the effect of the pronounced climatic gradient and soil depth variability on percolation fluxes and applied the calibrated model to a time series with 62 years of meteorological data. Soil moisture measurements showed a pronounced seasonality and suggested rapid infiltration during heavy rainstorms. Hydrus-1D successfully simulated short and long-term soil moisture patterns, with the majority of simulated deep percolation occurring during a few intensive rainfall events. Temperature drops in a nearby groundwater well were observed synchronously with simulated percolation pulses, indicating rapid groundwater recharge mechanisms. The 62-year model run yielded annual percolation fluxes of up to 66% of precipitation depths during wet years and of 0% during dry years. Furthermore, a dependence of recharge on the temporal rainfall distribution could be shown. Strong correlations between depth of recharge and soil depth were also observed.

Karst aquifers can be particularly vulnerable to human activities and climate change due to their relatively high degree of connection with the surface. This study utilized an ensemble of event-based recharge calculation methods to address the problem of structural uncertainty for the example of the Western Mountain Aquifer (WMA), a Mediterranean karst aquifer located in Israel and the West Bank. Spatially distributed recharge estimates derived from the Soil and Water Assessment Tool (SWAT) and the process-based infiltration model (PIM) were compared to site-specific, empirical regression models. The SWAT and PIM mean annual recharge estimates ranged from 32–34.6% of precipitation, almost equating to the results of empirical regression models (32–36%). Future recharge predictions under the influence of climate change were quantified by parameterizing the SWAT and PIM methods with a downscaled regional climate model of Israel. SWAT predicts a 23% decrease in recharge by 2051–2070 relative to 1981–2001. In contrast, PIM shows a 9% decrease, possibly due to the representation of infiltration through preferential flow pathways and exclusion of surface runoff processes. These divergent projections underline key methodological differences in the representation of hydrological processes. Nevertheless, both methods effectively provided good estimates of groundwater recharge. The recharge rates estimated from the various methods were integrated into MODFLOW to assess their relative impacts on groundwater storage dynamics. The ensemble of MODFLOW projected groundwater storage outputs can provide guidance for sustainable groundwater management in the region.

The objective of this work is the simulation of saturated and unsaturated flow in a karstified aquifer using a double continuum approach. The HydroGeoSphere code (Therrien et al., 2006) is employed to simulate spring discharge with the Richards equations and van Genuchten parameters to represent flow in the (1) fractured matrix and (2) conduit continuum coupled by a linear exchange term. Rapid vertical small-scale flow processes in the unsaturated conduit continuum are accounted for by applying recharge boundary conditions at the bottom of the saturated model domain. An extensive sensitivity analysis is performed on single parameters as well as parameter combinations. The transient hydraulic response of the karst spring is strongly controlled by the matrix porosity as well as the van Genuchten parameters of the unsaturated matrix, which determine the head dependent inter-continuum water transfer when the conduits are draining the matrix. Sensitivities of parameter combinations partially reveal a non-linear dependence over the parameter space. This can be observed for parameters not belonging to the same continuum as well as combinations, which involve the exchange parameter, showing that results of the double continuum model may depict a certain degree of ambiguity. The application of van Genuchten parameters for simulation of unsaturated flow in karst systems is critically discussed.

In this study the potential of an inversion approach based on hydraulic travel time and hydraulic attenuation tomography was assessed. Both hydraulic travel time and hydraulic attenuation tomography are based on the transformation of the transient groundwater flow equation into the eikonal equation using an asymptotic approach. The eikonal equation allows the calculation of pressure propagation and attenuation along trajectories, which is computationally efficient. The attenuation and travel time‐based inversion approaches are naturally complementary: hydraulic travel times are determined by the hydraulic diffusivity, a combination of hydraulic conductivity and specific storage, whereas the attenuation is determined solely by specific storage. The potential of our hydraulic tomographical approach was investigated at a well‐characterized sand and gravel aquifer located in the Leine River valley near Göttingen, Germany. The database for the hydraulic inversion consists of 392 cross‐well slug interference tests performed between five wells, in which the positions of the sources (injection ports) and the receivers (observation ports), isolated with double packer systems, were varied between tests. The results have shown that the combination of hydraulic travel time and hydraulic attenuation tomography allows the reconstruction of the diffusivity and storage distribution in two and three dimensions with a resolution and accuracy superior to that possible with type curve analysis.

Assessing the hydraulic parameters of karst aquifers is a challenge due to their high degree of heterogeneity. The unknown parameter field generally leads to a high ambiguity for flow and transport calibration in numerical models of karst aquifers. In this study, a distributed numerical model was built for the simulation of groundwater flow and solute transport in a highly heterogeneous karst aquifer in south-western Germany. Therefore, an interface for the simulation of solute transport in one-dimensional pipes was implemented into the software COMSOL Multiphysics® and coupled to the three-dimensional solute transport interface for continuum domains. For reducing model ambiguity, the simulation was matched for steady-state conditions to the hydraulic head distribution in the model area, the spring discharge of several springs and the transport velocities of two tracer tests. Furthermore, other measured parameters such as the hydraulic conductivity of the fissured matrix and the maximal karst conduit volume were available for model calibration. Parameter studies were performed for several karst conduit geometries to analyse the influence of the respective geometric and hydraulic parameters and develop a calibration approach in a large-scale heterogeneous karst system. Results show that it is possible not only to derive a consistent flow and transport model for a 150 km2 karst area but also to combine the use of groundwater flow and transport parameters thereby greatly reducing model ambiguity. The approach provides basic information about the conduit network not accessible for direct geometric measurements. The conduit network volume for the main karst spring in the study area could be narrowed down to approximately 100 000 m3.

Use of the subsurface for energy resources (enhanced geothermal systems, conventional and unconventional hydrocarbons), or for storage of waste (CO 2 , radioactive), requires the prediction of how fluids and the fractured porous rock mass interact. The GREAT cell (Geo-Reservoir Experimental Analogue Technology) is designed to recreate subsurface conditions in the laboratory to a depth of 3.5 km on 200 mm diameter rock samples containing fracture networks, thereby enabling these predictions to be validated. The cell represents an important new development in experimental technology, uniquely creating a truly polyaxial rotatable stress field, facilitating fluid flow through samples, and employing state of the art fibre optic strain sensing, capable of thousands of detailed measurements per hour. The cell’s mechanical and hydraulic operation is demonstrated by applying multiple continuous orientations of principal stress to a homogeneous benchmark sample, and to a fractured sample with a dipole borehole fluid fracture flow experiment, with backpressure. Sample strain for multiple stress orientations is compared to numerical simulations validating the operation of the cell. Fracture permeability as a function of the direction and magnitude of the stress field is presented. Such experiments were not possible to date using current state of the art geotechnical equipment.

Biodegradable hydrogels that promote stem cell differentiation into neurons in three dimensions (3D) are highly desired in biomedical research to study drug neurotoxicity or to yield cell-containing biomaterials for neuronal tissue repair. Here, we demonstrate that oxidized alginate-gelatin-laminin (ADA-GEL-LAM) hydrogels facilitate neuronal differentiation and growth of embedded human induced pluripotent stem cell (hiPSC) derived neurospheres. ADA-GEL and ADA-GEL-LAM hydrogels exhibiting a stiffness close to ~5 kPa at initial cell culture conditions of 37 °C were prepared. Laminin supplemented ADA-GEL promoted an increase in neuronal differentiation in comparison to pristine ADA-GEL, with enhanced neuron migration from the neurospheres to the bulk 3D hydrogel matrix. The presence of laminin in ADA-GEL led to a more than two-fold increase in the number of neurospheres with migrated neurons. Our findings suggest that laminin addition to oxidized alginate—gelatin hydrogel matrices plays a crucial role to tailor oxidized alginate-gelatin hydrogels suitable for 3D neuronal cell culture applications.