Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geothermal”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Reactive transport modeling of the Aquifer Thermal Energy Storage (ATES) system at Stockton University, New Jersey during seasonal operations

Hydrogeochemical processes associated with Aquifer Thermal Energy Storage (ATES) operations can often impact the system performance owing to mineral precipitation either at the wellbore or in the aquifer owing to changes in temperature and fluid disequilibria. Although failure of ATES systems due to mineral precipitation ("fouling") is common, predictive reactive-transport models have rarely been applied to plan their design and operation. Here, the objective of this study is to develop a reactive-transport model by coupling thermal, hydrological, and chemical (THC) processes to evaluate effects of introduced atmospheric oxygen on water chemistry, mineral precipitation/dissolution, porosity, and permeability changes associated with an ATES system at Stockton University (New Jersey, USA). The THC model builds on a Thermal-Hydrological-Mechanical (THM) model of the site that evaluated system failure owing to possible fracturing in the caprock or around the wellbore. The causes of the system failure are not known – potential causes include hydraulic fracturing owing to elevated pump pressures that took place, a flow pathway created by one of the boreholes, or a pre-existing natural hydrologic connection between the upper unconfined aquifer and the ATES aquifer, any of which could have led to oxygenated water entering the reservoir and causing the observed Fe-oxide fouling on well screens. The THC model is used to evaluate some of the hypotheses and observations regarding system failure owing to geochemical processes. The reactive-transport code TOUGHREACT V4 was used to model the THC processes during seasonal heating and cooling operations at the Stockton ATES site over 6 years of operation. In the THC simulations, the primary effects on geochemistry were observed when the injection water is saturated with atmospheric oxygen. Simulations show greater precipitation of goethite near the cold wells as compared to the warm wells. Although volume fractions of Fe-hydroxides were relatively small, the model was aimed at processes in the aquifer at the scale of meters and larger rather than at the scale of mm or cm (i.e., a well screen). Kaolinite is the dominant precipitating phase, also around the cold wells. Illite dissolves near the cold wells and precipitates near the warm wells. There is a net decrease in the porosity near the cold wells and increase near the warm wells, although a slight amount of thermal contraction near the cold wells and expansion near the warm wells is responsible for a significant proportion of the porosity change. Owing to the coarse discretization of the numerical grid near the wells (compared to the screen thickness) the magnitude of permeability changes at the wellbore are likely underestimated. The reactive transport model in this study can be used for characterization of aquifers, optimizing the operational parameters (temperature, pressure, pH etc.), and planning of mitigation strategies for ATES systems.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Assessing the impacts of air-sealing on the sizing, operation, and economic feasibility of ground-source heat pumps for electrifying single-family houses in the US

According to recent studies and reports, in single-family houses (SFHs), air-sealing can significantly lower the thermal loads for space heating and cooling. Thus, air-sealing in SFHs could reduce the required size and cost of ground source heat pump (GSHP) systems for electrifying SFHs. Here, this study investigated the costs and benefits of integrating air-sealing with GSHPs for retrofitting existing SFHs when compared with air-source heat pumps. A whole building energy simulation tool integrated with an advanced design tool for modeling ground heat exchangers was used to calculate changes in required GSHP capacity, total borehole length, and building energy consumption with and without air-sealing in SFHs in 16 US climatic regions. The results from this study showed that reducing outdoor air infiltration from 0.8 air changes per hour (ACH) to the minimum ventilation requirement (0.35 ACH) can significantly reduce borehole length (up to 55 %), GSHP capacity (up to 48 %), and total heating electricity reduction, especially in cold climates (up to 44 %). The results also showed that for airtight homes (0.03 ACH infiltration) with a direct outdoor air system, the minimum required borehole length, GSHP capacity, and total heating electricity consumption can be reduced up to 70 %, 68 %, and 67 %, respectively, when compared with SFHs with 0.8 ACH infiltration. Moreover, the life cycle cost analysis showed that air-sealing in conjunction with a GSHP is more profitable than replacing the existing system with an air-source heat pump, even without any incentives for most climatic regions in the US (except for some hot regions).

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Multiscale spatial analysis of fracture nodes in two dimensions

Spatial arrangement of fractures as a function of scale is an important component of fracture quantification for inferential and predictive modeling. Available methods that analyze fracture spatial arrangement are based on one-dimensional spacing data; therefore, they are limited to semi-parallel fractures. Such methods cannot be applied to fracture networks in higher dimensions, particularly when fractures have different orientations. Here to characterize fracture arrangements in two dimensions, we propose using Ripley’s K-function, as a method of point pattern analysis, to quantify spatial arrangement of fracture nodes. Fracture nodes, such as barycenters, intersection points, and tips, are point-based representations of fracture locations and connectivity within the fracture network. We introduce formulations for isotropic as well as directional analyses of spatial arrangement. In addition, we derive formulations for edge correction in circular and rectangular study domains. Finally, we demonstrate applications of Ripley’s K-function on two natural fracture datasets. Our proposed method supports quantification and characterization of fracture spatial arrangements that allow practitioners to build representative models of fractures in the subsurface.

02 PETROLEUM↗
DOE OSTI2022

Electrodialysis through nafion membranes coated with polyelectrolyte multilayers yields >99% pure monovalent ions at high recoveries

Separation of monovalent and multivalent ions is important in processes ranging from water softening to recycling of Li. Previous studies show that cation-exchange membranes coated with multilayer polyelectrolyte films exhibit monovalent/divalent ion selectivities as high as 1000 in electrodialysis (ED) under ideal conditions with low ion recoveries. This work shows that even at 60-80% monovalent-ion recoveries, ED through Nafion membranes coated with protonated poly(allylamine)/poly(4-styrenesulfonate) multilayers can separate equimolar K + and Mg 2+ or Li + and Mg 2+ to give 99.5% pure monovalent ion. Moreover, the use of a low-volume receiving phase allows ED to concentrate the monovalent ion. Typical current efficiencies are around 70% for these separations. Selectivities appear to increase when the membrane-coating procedure uses a higher polyelectrolyte concentration in the deposition solution. Because the limiting current declines with ED time due to depletion of the monovalent ion in the source phase, use of a gradually decreasing applied current enhances current efficiency compared to ED at a constant high current. Further, the membranes also show Li + /Mg 2+ selectivities around 100 when Mg 2+ is in a 20-fold excess in the source phase at an ionic strength of 0.6 M, but selectivity decreases at even higher ionic strengths.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Numerical study of multi-component flow and mixing in a scaled fission product venting system

Numerical simulation of isothermal multi-component gas flow was conducted to study fluid flow and mixing in a scaled low-pressure, low-temperature test facility for the fission product venting system (FPVS) of a gas cooled fast reactor (GFR). The geometry of FPVS test facility was an open loop including gradual expansion coupling and two 90°pipe elbows. First, the validation of large eddy simulation (LES) wall-adapting local eddy-viscosity (WALE) turbulent model was performed for transitional flow in a circular pipe with gradual expansion. The reactingFOAM solver in OpenFOAM v8 was employed. By introducing appropriate turbulence disturbance at the flow inlet, such as turbulence intensity, the numerical results showed good agreement in (1) the velocity profile downstream of the pipe expansion measured by particle image velocimetry (PIV) in this study, and (2) the reattachment length reported in the literature. Using this validated model, the numerical simulation of flow and component mixing was performed for the FPVS with an inlet flowrate corresponding to a Reynolds number of 2,400 to investigate the flow behavior and component mixing. Standard deviation of the mass fraction of component, referred to as absolute mixing index (AMI), was calculated to quantify the mixing of components, showing that the component mixing length determined by AMI is related to the reattachment length downstream of the expansion. Additionally, this study was able to identify the best location in the FPVS test facility to measure the velocity and concentration profies where multicomponent flow is well mixed.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Machine-learning-assisted high-temperature reservoir thermal energy storage optimization

High-temperature reservoir thermal energy storage (HT-RTES) has the potential to become an indispensable component in achieving the goal of the net-zero carbon economy, given its capability to balance the intermittent nature of renewable energy generation. In this study, a machine-learning-assisted computational framework is presented to co-optimize the performance metrics of HT-RTES by combining physics-based simulation with stochastic hydrogeologic formation and thermal energy storage operation parameters, artificial neural network regression of the simulation data, and genetic algorithm-enabled multi-objective optimization. A doublet well configuration with a layered (aquitard-aquifer-aquitard) generic reservoir is simulated for cases of continuous operation and seasonal-cycle operation scenarios. Further, neural network-based surrogate models are developed for the two scenarios and applied to generate the Pareto fronts of the HT-RTES performance for four potential HT-RTES sites. The developed Pareto optimal solutions indicate the performance of HT-RTES is operation-scenario (i.e., fluid cycle) and reservoir-site dependent, and the performance metrics have competing effects for a given site and a given fluid cycle. The developed neural network models can be applied to identify suitable sites for HT-RTES, and the proposed framework sheds light on the design of resilient HT-RTES systems.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Chitosan as a Canvas for Studies of Macromolecular Controls on CaCO 3 Biological Crystallization

A mechanistic understanding of how macromolecules, typically as an organic matrix, nucleate and grow crystals to produce functional biomineral structures remains elusive. Advances in structural biology indicate that polysaccharides (e.g., chitin) and negatively charged proteoglycans (due to carboxyl, sulfate, and phosphate groups) are ubiquitous in biocrystallization settings and play greater roles than currently recognized. This review highlights studies of CaCO 3 crystallization onto chitinous materials and demonstrates that a broader understanding of macromolecular controls on mineralization has not emerged. With recent advances in biopolymer chemistry, it is now possible to prepare chitosan-based hydrogels with tailored functional group compositions. By deploying these characterized compounds in hypothesis-based studies of nucleation rate, quantitative relationships between energy barrier to crystallization, macromolecule composition, and solvent structuring can be determined. Finally, this foundational knowledge will help researchers understand composition-structure-function controls on mineralization in living systems and tune the designs of new materials for advanced applications.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Potential Link Between 2020 Mentone, West Texas M5 Earthquake and Nearby Wastewater Injection: Implications for Aquifer Mechanical Properties

Abstract The M5 Mentone earthquake that occurred on March 26, 2020, was the largest event recorded over the last 2 decades in West Texas within the Delaware Basin, a U.S. major petroleum‐producing area. Also, numerous hydrofracturing and wastewater disposal wells are spread across this region. Within a 30 km distance to mainshock, eight class‐II injection wells for industrial wastewater disposal target the deep porous Ellenburger aquifer at an average rate of 1.36 × 10 6 barrel (BBL) per month during 2012–2020. Poroelastic models of fluid diffusion show these nearby injectors collectively imparted up to 80.5 kPa of Coulomb stress at the mainshock location, capable of triggering this M5 event. Assuming the Mentone event occurs when pore‐pressure increase is maximum, the time delay between peak injection and the M5 occurrence corresponds with an optimal permeability of 6.76 × 10 ‐14 m 2 for the Ellenburger aquifer layer, in agreement with independent estimates.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Frequency–Magnitude Statistics of Laboratory Foreshocks Vary With Shear Velocity, Fault Slip Rate, and Shear Stress

Understanding the nucleation phase of earthquakes is key for advancing earthquake hazard assessment and improving earthquake early warning systems. However, little progress has been made in this area due to a poor understanding of nucleation processes and incomplete seismic and fault zone measurements. The ability to integrate measured fault zone properties with seismic data could significantly improve our understanding of how earthquakes begin and whether there are systematic variations in seismic properties preceding failure. In this work, we use high-resolution laboratory measurements of fault zone properties along with acoustic emission data to document temporal variations of foreshock properties. Our data show that foreshock size increases with shear stress, loading rate, and fault slip rate. We propose that the preseismic fault slip rate and fault zone thickness (i.e., porosity) work in concert to modulate foreshock properties.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Non-Synchronous Rotation on Europa Driven by Ocean Currents

It has been suggested that the ice shell of Jupiter's moon Europa may drift non-synchronously due to tidal torques. Here we argue that torques applied by the underlying ocean are also important and can result in non-synchronous rotation. The resulting spin rate can be slightly slower than the synchronous angular rate that would have kept the same point of the ice shell facing Jupiter. We develop an ice shell rotation model, driven by ocean stress calculated using a high-resolution state-of-the-art ocean general circulation model, and take into account the viscoelastic deformation of the ice shell. We use the ice shell model results together with observed limits on the ice shell drift speed to constrain ice shell parameters such as effective viscosity, which is currently uncertain by at least four orders of magnitude. Our results suggest, at best, sluggish ice shell convection. Depending on the relaxation time scale of the ice shell and on the ocean currents, the ice shell may exhibit negligible drift, constant drift, or oscillatory drift superimposed on random fluctuations. The expected rotation rate exceeds ~30 m/yr; future spacecraft observations can be used to test these predictions and yield insight into the properties of the ice shell and underlying ocean.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Role of Intersections in Fracture Connectivity

Networks of intersecting fractures often provide the flow paths through subsurface reservoirs. Assessing network connectivity is challenging because fracture intersections compose a vanishingly small fraction of the network void volume. In this paper, motivated by 3D X-ray imaging of the simplest element of fracture network, that is, two orthogonal fractures, we perform a percolation and finite-size scaling analysis to study the connectivity provided by fracture intersections. The conditions when an intersection enhances connectivity across a sample depend on spatial correlations in the fracture aperture distributions, on the stress state, and on the direction of flow. Here we consider three flow directions: (a) across intersections, (b) parallel to intersections and (c) around corners. For (a), intersections provide minimal enhancement of connectivity because they contribute little additional void area. For (b), intersections increase the probability of a connected path near threshold by enabling 3D connected pathways that are not possible in parallel fractures. Flow around corners, (c), is fundamentally the result of the intersection connecting two fractures in series and spatial correlations are broken around corners, suppressing the connectivity relative to (a). When the connected fractures are stressed equally, a joint percolation threshold emerges that continues to have scale invariance. However, when the fractures are stressed unequally, the system has mixed percolation without clearly defined percolation thresholds. In all cases, percolation probabilities are found to be scale dependent which has important consequences for the connectivity of larger fracture networks composed of the fundamental element studied here.

02 PETROLEUM↗