Modeling, simulation, and optimization of geothermal energy production from hot sedimentary aquifers

Geothermal district heating development has been gaining momentum in Europe with numerous deep geothermal installations and projects currently under development. With the increasing density of geothermal wells, questions related to the optimal and sustainable reservoir exploitation become more and more important. A quantitative understanding of the complex thermo-hydraulic interaction between tightly deployed geothermal wells in heterogeneous temperature and permeability fields is key for a maximum sustainable use of geothermal resources. Motivated by the geological settings of the Upper Jurassic aquifer in the Greater Munich region, we develop a computational model based on finite element analysis and gradient-free optimization to simulate groundwater flow and heat transport in hot sedimentary aquifers, and investigate numerically the optimal positioning and spacing of multi-well systems. Based on our numerical simulations, net energy production from deep geothermal reservoirs in sedimentary basins by smart geothermal multi-well arrangements provides significant amounts of energy to meet heat demand in highly urbanized regions. Our results show that taking into account heterogeneous permeability structures and variable reservoir temperature may drastically affect the results in the optimal configuration. We demonstrate that the proposed numerical framework is able to efficiently handle generic geometrical and geologocal configurations, and can be thus flexibly used in the context of multi-variable optimization problems. Hence, this numerical framework can be used to assess the extractable geothermal energy from heterogeneous deep geothermal reservoirs by the optimized deployment of smart multi-well systems

Modeling, simulation, and optimization of geothermal energy production from hot sedimentary aquifers

Geothermal district heating development has been gaining momentum in Europe with numerous deep geothermal installations and projects currently under development. With the increasing density of geothermal wells, questions related to the optimal and sustainable reservoir exploitation become more and more important. A quantitative understanding of the complex thermo-hydraulic interaction between tightly deployed geothermal wells in heterogeneous temperature and permeability fields is key for a maximum sustainable use of geothermal resources. Motivated by the geological settings of the Upper Jurassic aquifer in the Greater Munich region, we develop a computational model based on finite element analysis and gradient-free optimization to simulate groundwater flow and heat transport in hot sedimentary aquifers, and investigate numerically the optimal positioning and spacing of multi-well systems. Based on our numerical simulations, net energy production from deep geothermal reservoirs in sedimentary basins by smart geothermal multi-well arrangements provides significant amounts of energy to meet heat demand in highly urbanized regions. Our results show that taking into account heterogeneous permeability structures and variable reservoir temperature may drastically affect the results in the optimal configuration. We demonstrate that the proposed numerical framework is able to efficiently handle generic geometrical and geologocal configurations, and can be thus flexibly used in the context of multi-variable optimization problems. Hence, this numerical framework can be used to assess the extractable geothermal energy from heterogeneous deep geothermal reservoirs by the optimized deployment of smart multi-well systems

Analysis of a stabilized penalty-free Nitsche method for the Brinkman, Stokes, and Darcy problems

International audienceIn this paper we study the Brinkman model as a unified framework to allow the transition between the Darcy and the Stokes problems. We propose an unconditionally stable low-order finite element approach, which is robust with respect to the whole range of physical parameters, and is based on the combination of stabilized equal-order finite elements with a non-symmetric penalty-free Nitsche method for the weak imposition of essential boundary conditions. In particular, we study the properties of the penalty-free Nitsche formulation for the Brinkman setting, extending a recently reported analysis for the case of incompressible elasticity (Boiveau and Burman, IMA J. Numer. Anal. 36 (2016) 770-795). Focusing on the two-dimensional case, we obtain optimal a priori error estimates in a mesh-dependent norm, which, converging to natural norms in the cases of Stokes or Darcy ows, allows to extend the results also to these limits. Moreover, we show that, in order to obtain robust estimates also in the Darcy limit, the formulation shall be equipped with a Grad-Div stabilization and an additional stabilization to control the discontinuities of the normal velocity along the boundary. The conclusions of the analysis are supported by numerical simulations

Impact of Alkyl Chain Length on the Formation of Regular- and Reverse-Graded Quasi-2D Perovskite Thin Films

Crystallization of low-dimensional perovskites is a complex process that leads to multidimensional films comprising two-dimensional (2D), quasi-2D, and three-dimensional (3D) phases. Most quasi-2D perovskite films possess a regular gradient with 2D phases located at the bottom of the film and 3D phases at the top. Recently, multiple studies have reported reverse-graded perovskite films, where the location of the 2D and 3D structures is inverted. The underlying reasons for such a peculiar phase distribution are unclear. While crystallization of regular-graded quasi-2D perovskites has been described as starting with 3D phases from the liquid-air interface, the film formation of reverse-graded films has not been investigated yet. Here, we examine the impact of the alkyl chain length on the formation of regular- and reverse-graded perovskites using n-alkylammonium ions. We find that long alkyl chains reverse the phase distribution gradient. By combining photoluminescence spectroscopy with in situ optical absorption measurements, we demonstrate that crystallization starts at the liquid-N2 interface, though as 3D phases for short-chain n-alkylammonium ions and as quasi-2D phases for long chains. We link this behavior to enhanced van der Waals interactions between long-chain n-alkylammonium ions in polar solvents and their tendency to accumulate at the liquid-N2 interface, creating a concentration gradient along the film thickness.</p

ParMooN - a modernized program package based on mapped finite elements

{\sc ParMooN} is a program package for the numerical solution of elliptic and parabolic partial differential equations. It inherits the distinct features of its predecessor {\sc MooNMD} \cite{JM04}: strict decoupling of geometry and finite element spaces, implementation of mapped finite elements as their definition can be found in textbooks, and a geometric multigrid preconditioner with the option to use different finite element spaces on different levels of the multigrid hierarchy. After having presented some thoughts about in-house research codes, this paper focuses on aspects of the parallelization for a distributed memory environment, which is the main novelty of {\sc ParMooN}. Numerical studies, performed on compute servers, assess the efficiency of the parallelized geometric multigrid preconditioner in comparison with some parallel solvers that are available in the library {\sc PETSc}. The results of these studies give a first indication whether the cumbersome implementation of the parallelized geometric multigrid method was worthwhile or not.Comment: partly supported by European Union (EU), Horizon 2020, Marie Sk{\l}odowska-Curie Innovative Training Networks (ITN-EID), MIMESIS, grant number 67571

Volcanic Fluxes Over the Last Millennium as Recorded in the Gv7 Ice Core (Northern Victoria Land, Antarctica)

Major explosive volcanic eruptions may significantly alter the global atmosphere for about 2-3 years. During that period, volcanic products (mainly H2SO4) with high residence time, stored in the stratosphere or, for shorter times, in the troposphere are gradually deposited onto polar ice caps. Antarctic snow may thus record acidic signals providing a history of past volcanic events. The high resolution sulphate concentration profile along a 197 m long ice core drilled at GV7 (Northern Victoria land) was obtained by Ion Chromatography on around 3500 discrete samples. The relatively high accumulation rate (241 +/- 13 mm we yr (-1)) and the 5-cm sampling resolution allowed a preliminary counted age scale. The obtained stratigraphy covers roughly the last millennium and 24 major volcanic eruptions were identified, dated, and tentatively ascribed to a source volcano. The deposition flux of volcanic sulphate was calculated for each signature and the results were compared with data from other Antarctic ice cores at regional and continental scale. Our results show that the regional variability is of the same order of magnitude as the continental one

Characterization of free L- and D-amino acids in size-segregated background aerosols over the Ross Sea, Antarctica

The study of airborne chemical markers is crucial for identifying sources of aerosols, and their atmospheric processes of transport and transformation. The investigation of free amino acids and their differentiation between the L-and D- en-antiomers are even more important to understand their sources and atmospheric fate. Aerosol samples were collected with a high-volume sampler with cascade impactor at Mario Zucchelli Station (MZS) on the coast of the Ross Sea (Antarctica) for two summer campaigns (2018/19 and 2019/20). The total mean concentration of free amino acids in PM10 was 4 +/- 2 pmol m-3 for both campaigns and most of free amino acids were distributed in fine particles. The coarse mode of airborne D-Alanine and dimethylsufoniopropionate in seawater showed a similar trend during both Antarctic campaigns. Thus, the study of D/L Ala ratio in fine, coarse and PM10 fractions indicated the microlayer as the local source. This paper demonstrated that free amino acids follow the trend of DMS and MSA release occurred in the Ross Sea, confirming their applicability as markers for phytoplankton bloom also in paleoclimatic studies

ParMooN - a modernized program package based on mapped finite elements

PARMOON is a program package for the numerical solution of elliptic and parabolic partial differential equations. It inherits the distinct features of its predecessor MOONMD [28]: strict decoupling of geometry and finite element spaces, implementation of mapped finite elements as their definition can be found in textbooks, and a geometric multigrid preconditioner with the option to use different finite element spaces on different levels of the multigrid hierarchy. After having presented some thoughts about in-house research codes, this paper focuses on aspects of the parallelization, which is the main novelty of PARMOON. Numerical studies, performed on compute servers, assess the efficiency of the parallelized geometric multigrid preconditioner in comparison with parallel solvers that are available in external libraries. The results of these studies give a first indication whether the cumbersome implementation of the parallelized geometric multigrid method was worthwhile or not

Lmx1a-Dependent Activation of miR-204/211 Controls the Timing of Nurr1-Mediated Dopaminergic Differentiation

The development of midbrain dopaminergic (DA) neurons requires a fine temporal and spatial regulation of a very specific gene expression program. Here, we report that during mouse brain development, the microRNA (miR-) 204/211 is present at a high level in a subset of DA precursors expressing the transcription factor Lmx1a, an early determinant for DA-commitment, but not in more mature neurons expressing Th or Pitx3. By combining different in vitro model systems of DA differentiation, we show that the levels of Lmx1a influence the expression of miR-204/211. Using published transcriptomic data, we found a significant enrichment of miR-204/211 target genes in midbrain dopaminergic neurons where Lmx1a was selectively deleted at embryonic stages. We further demonstrated that miR-204/211 controls the timing of the DA differentiation by directly downregulating the expression of Nurr1, a late DA differentiation master gene. Thus, our data indicate the Lmx1a-miR-204/211-Nurr1 axis as a key component in the cascade of events that ultimately lead to mature midbrain dopaminergic neurons differentiation and point to miR-204/211 as the molecular switch regulating the timing of Nurr1 expression

Biogenic Aerosol in the Artic from Eight Years of MSA Data from Ny Ålesund (Svalbard Islands) and Thule (Greenland)

In remote marine areas, biogenic productivity and atmospheric particulate are coupled through dimethylsulfide (DMS) emission by phytoplankton. Once in the atmosphere, the gaseous DMS is oxidized to produce H2SO4 and methanesulfonic acid (MSA); both species can a ect the formation of cloud condensation nuclei. This study analyses eight years of biogenic aerosol evolution and variability at two Arctic sites: Thule (76.5 N, 68.8 W) and Ny Ålesund (78.9 N, 11.9 E). Sea ice plays a key role in determining the MSA concentration in polar regions. At the beginning of the melting season, in April, up to June, the biogenic aerosol concentration appears inversely correlated with sea ice extent and area, and positively correlated with the extent of the ice-free area in the marginal ice zone (IF-MIZ). The upper ocean stratification induced by sea ice melting might have a role in these correlations, since the springtime formation of this surface layer regulates the accumulation of phytoplankton and nutrients, allowing the DMS to escape from the sea to the atmosphere. The multiyear analysis reveals a progressive decrease in MSA concentration in May at Thule and an increase in July August at Ny Ålesund. Therefore, while the MSA seasonal evolution is mainly related with the sea ice retreat in April, May, and June, the IF-MIZ extent appears as the main factor a ecting the longer-term behavior of MSA.MIUR PRIN 2007 and PRIN 2009, PNRA 2010-2012, PNRA 2015-2016, PNRA 2016-2018Publishedid 3495A. Ricerche polari e paleoclimaJCR Journa