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,189 records · Page 66

Lab on a chip for a low-carbon future

Transitioning our society to a sustainable future, with low or net-zero carbon emissions to the atmosphere, will require a wide-spread transformation of energy and environmental technologies. In this perspective article, we describe how lab-on-a-chip (LoC) systems can help address this challenge by providing insight into the fundamental physical and geochemical processes underlying new technologies critical to this transition, and developing the new processes and materials required. We focus on six areas: (I) subsurface carbon sequestration, (II) subsurface hydrogen storage, (III) geothermal energy extraction, (IV) bioenergy, (V) recovering critical materials, and (VI) water filtration and remediation. Here, we hope to engage the LoC community in the many opportunities within the transition ahead, and highlight the potential of LoC approaches to the broader community of researchers, industry experts, and policy makers working toward a low-carbon future.

42 ENGINEERING↗

Techno-Economic Analysis of Repurposing Natural Gas Transmission Pipeline Networks to Accommodate Hydrogen Blends

Blending hydrogen into natural gas infrastructure could supplement natural gas supply and increase resilience for applications such as ammonia production, peaking and load-following power plants, and heating. The United States has an extensive network of natural gas pipelines, but the feasibility of employing this infrastructure to transport hydrogen is unclear. We analyze the costs associated with repurposing three distinct natural gas transmission pipelines in different locations within the United States to carry blends of hydrogen up to 100% via three different pipeline network modification methods and compare against the cost of building a new dedicated hydrogen pipeline. We conduct a sensitivity analysis on the hoop stress limit of the existing pipe, techno-economic parameters, emissions, and relative capacity. The results show that the capital costs required to upgrade a pipeline can vary from tens of millions to billions of dollars depending on the length and capacity of the existing pipeline section, whether the existing pipeline already operates at or below its maximum allowable operating pressure, whether future demand is expected to increase or decrease, the network modification method selected, and pipe material costs. The delivered cost of energy to end users is impacted less by the levelized cost of transporting hydrogen blends than by the cost of the natural gas and hydrogen fuels being transported. The emissions impact of blending hydrogen into natural gas transmission networks scales proportionally with the amount of energy displaced with low-emission hydrogen (such as from natural gas with carbon capture and sequestration or electrolysis powered by nuclear, renewable, or geothermal electricity), therefore low blend ratios (e.g., < 20% vol. hydrogen) will result in low emissions impacts. Factors such as permitting and right-of-way costs, the proximity of the pipeline to hydrogen demand and production, and the compatibility of and/or retrofitting costs of end-use gas-fueled technologies will likely be greater drivers in determining whether converting a particular natural gas pipeline to carry hydrogen makes economic sense.

08 HYDROGEN↗

Stimulated geologic hydrogen: from mechanistic control to engineered rock transformation

Geologic hydrogen (GeoH 2 ) generated from subsurface iron-rich rock–water reaction (i.e., serpentinization) is emerging as a promising candidate for the next primary energy source. Yet, accelerating GeoH 2 production from geological to human timescales via enhanced serpentinization remains a formidable scientific and technical challenge. Here, in this Review, we decipher the mechanistic control and explore strategies for accelerating in situ, engineered iron-rich rock transformation into carbon-free GeoH 2 by orders of magnitude. Serpentinization rate is hindered by low porosity and permeability of source rocks, suboptimal temperatures, unfavorable water chemistry, inefficient Fe 2+ -to-Fe 3+ conversion, thermodynamic constraints, and low reactive surface area. While closed-system experiments provide valuable mechanistic insights, open-system conditions with fluid circulation are more crucial for economically viable GeoH 2 production. We assess stimulation techniques from enhanced hydrocarbon and geothermal recovery as tools to be adopted or adapted for increasing reactive surface areas for stimulated GeoH2 production. We estimate that 7.40 × 10 5 to 1.73 × 10 6 million metric tons (Mt) of hydrogen could be engineered over 20 to 50 years from about 10% iron-rich rocks within 10 km depth of continental crust. Enabling GeoH 2 as a viable energy source requires not only advancing scientific frontiers but also forming a global GeoH 2 research network and innovation ecosystem to address the critical scientific, technical, societal, economic, and policy challenges.

08 HYDROGEN↗

Underwater plasma breakdown characteristics with respect to highly pressurized drilling applications

Deep earth drilling is a key technique to extract oil, gas, and geothermal heat from the earth. Many complex energy focusing methods have been explored as an alternative approach to reach these resources but most of them require high energy. However, by utilizing short time span liquid plasma discharges, energy focusing can be achieved within traditional drilling systems. These discharges induce a rapid expansion process and a resulting shockwave. It is believed that this focused energy will lower the required cutting force to progress through the rock. Lowering the required cutting force will allow for lower drill bit wear, quicker rate of penetration, and an overall cost savings of the project. Plasma breakdown characteristics at drilling relevant pressures, ranging from 1 to 350 atm, were studied. A resistance–capacitance circuit with an air gapped spark switch was utilized to generate pulsed plasma between the pressurized electrodes. It was found that the required breakdown voltage increases as the pressure increases. It was also found that a plasma channel formation and an associated breakdown may or may not occur between the electrodes at different pressures due to variation in required breakdown voltages. Breakdown time-lag in the dielectric medium (tap water, 780 μS cm –1 ) increased as the pressure was increased, which indicated a higher voltage drop at higher pressures (>100 atm). The plasma generated cavitation bubble with an associated shockwave occurred as pressures were increased. However, the bubble radius and the bubble duration decreased as the pressure was increased. The plasma generated shockwave speeds fall within the expected speed of sound in water. Lastly, preliminary rock cracking tests were performed on granite at high pressures (340 and 272 atm) and it was found that plasma is able to create cracks in the rock.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

AnisONet: A deep neural operator-based anisotropic permeability upscaler from pore to Darcy scale

Directional permeability variations, which govern directional fluid flow in porous media with anisotropy, are important to accurately predict flow behavior, reactive transport, and fluid–solid interactions for various processes such as enhanced geothermal systems, energy storage devices, and biological systems. However, the intricate architecture of porous media makes it difficult to predict directional permeabilities. In this work, we present a novel machine learning (ML) framework, AnisONet, built upon an integration of a convolutional neural network, Swin transformer, and the deep operator network architecture, designed to predict anisotropic permeability and upscale predictions to larger spatial domains. First, AnisONet was evaluated with three classes of two-dimensional (2D) porous media, including synthetic circular and elliptical grains and natural sandstone grains from micro-computed tomography images. A lattice Boltzmann model (LBM) was used to calculate directional permeabilities at every 10° angle, producing 19 data points per image of porous media. AnisONet is then trained to predict permeability as a function of rotation angle. AnisONet showed strong predictive capability of directional permeability. Second, we tested our model for five upscaling cases with a large image size in the finite-element method (FEM) for 2D Darcy flow with various permeability tensor construction methods. Overall, upscaled permeability tensors in FEM simulations produce a reasonably good match with LBM results, highlighting the importance of selecting appropriate tensor formation strategies for accurate permeability upscaling. AnisONet, as a directional permeability estimator, could be further developed for more complex geometries, with the potential to develop a foundational ML model for various applications in porous media.

42 ENGINEERING↗

Evidence for the stability of ultrahydrous stishovite in Earth’s lower mantle

The distribution and transportation of water in Earth’s interior depends on the stability of water-bearing phases. The transition zone in Earth’s mantle is generally accepted as an important potential water reservoir because its main constituents, wadsleyite and ringwoodite, can incorporate weight percent levels of H 2 O in their structures at mantle temperatures. The extent to which water can be transported beyond the transition zone deeper into the mantle depends on the water carrying capacity of minerals stable in subducted lithosphere. Stishovite is one of the major mineral components in subducting oceanic crust, yet the capacity of stishovite to incorporate water beyond at lower mantle conditions remains speculative. Here, we combine in situ laser heating with synchrotron X-ray diffraction to show that the unit cell volume of stishovite synthesized under hydrous conditions is ~2.3 to 5.0% greater than that of anhydrous stishovite at pressures of ~27 to 58 GPa and temperatures of 1,240 to 1,835 K. Our results indicate that stishovite, even at temperatures along a mantle geotherm, can potentially incorporate weight percent levels of H 2 O in its crystal structure and has the potential to be a key phase for transporting and storing water in the lower mantle.

58 GEOSCIENCES↗

Fluid inertia controls mineral precipitation and clogging in pore to network-scale flows

Mineral precipitation caused by fluid mixing presents complex control and predictability challenges in a variety of natural and engineering processes, including carbon mineralization, geothermal energy, and microfluidics. Precipitation dynamics, particularly under the influence of fluid flow, remain poorly understood. Combining microfluidic experiments and three-dimensional reactive transport simulations, we demonstrate that fluid inertia controls mineral precipitation and clogging at flow intersections, even in laminar flows. We observe distinct precipitation regimes as a function of Reynolds number (Re). At low Reynolds numbers (Re < 10), precipitates form a thin, dense layer along the mixing interface, which shuts precipitation off, while at high Reynolds numbers (Re > 50), strong three-dimensional flows significantly enhance precipitation over the entire intersection, resulting in rapid clogging. When injection rates from two inlets are uneven, flow symmetry-breaking leads to unexpected flow bifurcation phenomena, which result in enhanced concurrent precipitation in both downstream channels. Finally, we extend our findings to rough channel networks and demonstrate that the identified inertial effects on precipitation at the intersection scale are also present and even more dramatic at the network scale. This study sheds light on the fundamental mechanisms underlying mixing-induced mineral precipitation and provides a framework for designing and optimizing processes involving mineral precipitation.

Science & Technology - Other Topics↗

Inertia-induced mixing and reaction maximization in laminar porous media flows

Solute transport and biogeochemical reactions in porous and fractured media flows are controlled by mixing, as are subsurface engineering operations such as contaminant remediation, geothermal energy production, and carbon sequestration. Porous media flows are generally regarded as slow, so the effects of fluid inertia on mixing and reaction are typically ignored. Here, we demonstrate through microfluidic experiments and numerical simulations of mixing-induced reaction that inertial recirculating flows readily emerge in laminar porous media flows and dramatically alter mixing and reaction dynamics. An optimal Reynolds number that maximizes the reaction rate is observed for individual pore throats of different sizes. This reaction maximization is attributed to the effects of recirculation flows on reactant availability, mixing, and reaction completion, which depend on the topology of recirculation relative to the boundary of the reactants or mixing interface. Recirculation enhances mixing and reactant availability, but a further increase in flow velocity reduces the residence time in recirculation, leading to a decrease in reaction rate. The reaction maximization is also confirmed in a flow channel with grain inclusions and randomized porous media. Interestingly, the domain-wide reaction rate shows a dramatic increase with increasing Re in the randomized porous media case. This is because fluid inertia induces complex three-dimensional flows in randomized porous media, which significantly increases transverse spreading and mixing. This study shows how inertial flows control reaction dynamics at the pore scale and beyond, thus having major implications for a wide range of environmental systems.

Chen, Michael A. (ORCID:0000000294614710)↗

Three-chamber electrochemical reactor for selective lithium extraction from brine

Efficient lithium recovery from geothermal brines is crucial for the battery industry. Current electrochemical separation methods struggle with the simultaneous presence of Na + , K + , Mg 2+ , and Ca 2+ because these cations are similar to Li + , making it challenging to separate effectively. We address these challenges with a three-chamber reactor featuring a polymer porous solid electrolyte in the middle layer. This design improves the transference number of Li + (t Li+ ) by 2.1 times compared to the two-chamber reactor and also reduces the chlorine evolution reaction, a common side reaction in electrochemical lithium extraction, to only 6.4% in Faradaic Efficiency. Employing a lithium-ion conductive glass ceramic (LICGC) membrane, the reactor achieved high t Li+ of 97.5% in LiOH production from simulated brine, while the concentrations of Na + K + , Mg 2+ , and Ca 2+ are below the detection limit. Electrochemical experiments and surface analysis elucidated the cation transport mechanism, highlighting the impact of Na + on Li + migration at the LICGC interface.

Science & Technology - Other Topics↗

Cement sensors with acoustic bandgaps using carbon nanotubes

Cement is widely used in wellbores to stabilize the steel casing used in wellbore operations for oil and gas production, enhanced geothermal systems and carbon sequestration, and to limit fluid movement between sub-surface strata. Flaws such as microcracks in wellbore cement can lead to leakage along the wellbore compromising wellbore integrity. There is an increasing need for methods to monitor cement crack propagation in wellbore environments. In this study, we develop and report the first cementitious sensors capable of exhibiting high frequency acoustic bandgaps (ABGs) using carbon nanotubes (CNTs). Computational simulations of a sensor unit cell are used to design cement-multi walled carbon nanotubes (MWCNTs) sensors that show a wide bandgap. When the cement-MWCNTs sensors is embedded in cement specimens, bandgaps were measured experimentally under 300 kHz and under 600 kHz, consistent with the computationally predicted bandgaps in the range of 290–360 kHz, 410–460 kHz and 515–585 kHz. These bandgap features were absent in homogeneous cement specimens. X-ray tomographic reconstructions showed microscopic debonding at cement-MWCNTs sensor interface. Frequency response analysis of a three-dimensional computational model indicated a shift of frequency of minimum transmission due to the interface debonding, but no perturbation of bandgap response was observed. Here, the cement-MWCNTs sensors developed in this study show the potential of a packed CNT inclusion material in cementitious matrix to create ABGs in a cement matrix.

36 MATERIALS SCIENCE↗

Environmental risks and opportunities of orphaned oil and gas wells in the United States

Abstract Hundreds of thousands of documented and undocumented orphaned oil and gas wells exist in the United States (U.S.). These wells have the potential to contaminate water supplies, degrade ecosystems, and emit methane and other air pollutants. Thus, orphaned wells present risks to climate stability and to environmental and human health, which can be reduced by plugging. To quantify environmental risks and opportunities of well plugging at the national level, we analyze data on 81 857 documented orphaned wells across the U.S. We find that > 4.6 million people live within 1 km of a documented orphaned well. 35% of the documented orphaned wells are located within 1 km of a domestic groundwater well, yet only 8% of the wells have groundwater quality data within a 1 km radius. Methane emissions from the documented orphaned wells represent approximately 3%–6% of total U.S. methane emissions from abandoned oil and gas wells, but this estimate is based on measurements at < 0.03% of U.S. abandoned wells. 91% of the documented orphaned wells overlie formations favorable for geologic storage of carbon dioxide and hydrogen, meaning that orphaned well plugging can reduce leakage risks from future storage projects. Finally, we estimate plugging costs for documented orphaned wells to exceed the $4.7 billion federal funding by 30%–80%, emphasizing the importance of prioritizing federal spending on wells with large remediation benefits. Overall, environmental monitoring data are not extensive enough to quantify risks, especially those related to air and water quality and human health. Plugging orphaned wells can provide opportunities for geologic storage of carbon dioxide and hydrogen and geothermal energy development, thereby facilitating efforts to transition to net-zero energy systems. Our analysis on environmental risks and opportunities of orphaned wells provides a framework that can be used to manage the millions of documented and undocumented orphaned wells in the U.S. and abroad.

54 ENVIRONMENTAL SCIENCES↗

High-resolution Bayesian spatial autocorrelation (SPAC) quasi-3-D Vs model of Utah FORGE site with a dense geophone array

SUMMARY We expand the application of spatial autocorrelation (SPAC) from typical 1-D Vs profiles to quasi-3-D imaging via Bayesian Monte Carlo inversion (BMCI) using a dense nodal array (49 nodes) located at the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) site. Combinations of 4 and 9 geophones in subarrays provide for 36 and 25 1-D Vs profiles, respectively. Profiles with error bars are determined by calculating coherency functions that fit observations in a frequency range of 0.2–5 Hz. Thus, a high-resolution quasi-3-D Vs model from the surface to 2.0 km depth is derived and shows that surface-parallel sedimentary strata deepen to the west, consistent with a 3-D seismic reflection survey. Moreover, the resulting Vs profile is consistent with a Vs profile derived from distributed acoustic sensing (DAS) data located in a borehole at the FORGE site. The quasi-3-D velocity model shows that the base of the basin dips ∼22° to the west and topography on the basement interface coincident with the Mag Lee Wash suggests that the bedrock interface is an unconformity.

Geochemistry & Geophysics↗

Regional seismic velocity changes following the 2019 M w 7.1 Ridgecrest, California earthquake from autocorrelations and P / S converted waves

SUMMARY We examine regional transient changes of seismic velocities generated by the 2019 Mw 7.1 Ridgecrest earthquake in California, using autocorrelations of moving time windows in continuous waveforms recorded at regional stations. We focus on traveltime differences in a prominent phase generated by an interface around 2 km depth, associated with transmitted Pp waves and converted Ps waves from the ongoing microseismicity. Synthetic tests demonstrate the feasibility of the method for monitoring seismic velocity changes. Taking advantage of the numerous aftershocks in the early period following the main shock, we obtain a temporal resolution of velocity changes up to 20 min in the early post-main-shock period. The results reveal regional coseismic velocity drops in the top 1–3 km with an average value of ∼2 per cent over distances up to 100 km from the Ridgecrest event. These average velocity drops are likely dominated by larger changes in the shallow materials and are followed by rapid recoveries on timescales of days. Around the north end of the Ridgecrest rupture and the nearby Coso geothermal region, the observed coseismic velocity drops are up to ∼8 per cent. The method allows monitoring temporal changes of seismic velocities with high temporal resolution, fast computation and precise spatial mapping of changes. The results suggest that significant temporal changes of seismic velocities of shallow materials are commonly generated on a regional scale by large events.

58 GEOSCIENCES↗

Mineralogy of the deep lower mantle in the presence of H 2 O

Understanding the mineralogy of the Earth's interior is a prerequisite for unravelling the evolution and dynamics of our planet. Here, we conducted high pressure-temperature experiments mimicking the conditions of the deep lower mantle (DLM, 1800–2890 km in depth) and observed surprising mineralogical transformations in the presence of water. Ferropericlase, (Mg, Fe)O, which is the most abundant oxide mineral in Earth, reacts with H 2 O to form a previously unknown (Mg, Fe)O 2 H x (x ≤ 1) phase. The (Mg, Fe)O 2 H x has a pyrite structure and it coexists with the dominant silicate phases, bridgmanite and post-perovskite. Depending on Mg content and geotherm temperatures, the transformation may occur at 1800 km for (Mg 0.6 Fe 0.4 )O or beyond 2300 km for (Mg 0.7 Fe 0.3 )O. The (Mg, Fe)O 2 H x is an oxygen excess phase that stores an excessive amount of oxygen beyond the charge balance of maximum cation valences (Mg 2+ , Fe 3+ and H + ). This important phase has a number of far-reaching implications including extreme redox inhomogeneity, deep-oxygen reservoirs in the DLM and an internal source for modulating oxygen in the atmosphere.

58 GEOSCIENCES↗

Dominant effect of anharmonicity on the equation of state and thermal conductivity of MgO under extreme conditions

Understanding heat transfer in the earth mantle region is of great scientific interest but highly challenging due to dissimilar conditions from our daily lives. Thermal conductivity of highly pressured and thermally activated materials can be completely different from that at ambient condition. In this study we calculate the equation of state and thermal conductivity of cubic crystalline MgO, a major component of the mantle, in such an environment. It is shown that the material properties are not accurately captured unless rigorous treatment of anharmonicity of phonon vibrational motion is explicitly applied. To predict renormalized phonon dispersion relations, interatomic force constants of high-order expansion are required, especially at elevated temperature. The anharmonicity plays an increasingly substantial role for regulating volume expansion as temperature rises, because the Gibbs free energy of MgO considerably deviates from what was predicted by the quasiharmonic approximation. We find that the lattice thermal conductivity of MgO and the heat transfer of the earth's mantle are reliably estimated only when both the temperature-dependent renormalization of phonon frequencies and four-phonon scattering effects are taken into consideration. We map the characteristics at Geotherm into contour diagrams in which equilibrium volumes and the thermal conductivities are described as functions of temperature and pressure. Furthermore, the thermal conductivity values of MgO in this study are smaller than those in previous theoretical studies due to the additional incorporation of four-phonon scatterings in the calculations. Specifically, deviation between the values in previous studies and in this study becomes larger as thermodynamic variables reach those near the core-mantle boundary. For instance, thermal conductivities of 52.15 and 23.67 Wm -1 K -1 are estimated at ambient and extreme conditions (T = 4000K and P = 135 GPa), respectively.

36 MATERIALS SCIENCE↗

Anharmonic thermodynamic properties and phase boundary across the postperovskite transition in MgSi O 3

To address the effects of lattice anharmonicity across the perovskite to postperovskite transition in MgSiO 3 , we conduct calculations using the phonon quasiparticle (PHQ) approach. The PHQ is based on abinitio molecular dynamics and, in principle, captures full anharmonicity. Free energies in the thermodynamic limit (N → ∞) are computed using temperature-dependent quasiparticle dispersions within the phonon gas model. Systematic results on anharmonic thermodynamic properties and phase boundary are reported. Both the local density approximation and the generalized gradient approximation calculations are performed to provide confident constraints on these properties. Anharmonic effects are demonstrated by comparing results with those obtained using the quasiharmonic approximation (QHA). The inadequacy of the QHA is indicated by its overestimation of thermal expansivity and thermodynamic Grüneisen parameter and its converged isochoric heat capacity in the high-temperature limit. The PHQ phase boundary has a Clapeyron slope (dP/dT) that increases with temperature. This result contrasts with the nearly zero curvature of the QHA phase boundary. Anharmonicity bends the phase boundary to lower temperatures at high pressures. Implications for the double-crossing of the phase boundary by the mantle geotherm are discussed.

36 MATERIALS SCIENCE↗

A machine learning degradation model for electrochemical capacitors operated at high temperature

Electrochemical capacitors (ECs) have only recently been considered as an alternative power source for telemetry sensors of drilling equipment for geothermal or oil and gas exploration. The lifecycle analysis and modelling of ECs is underrepresented in literature in comparison to other storage devices e.g. Li-ion batteries. This paper investigates the degradation of ECs when cycled outside the manufacturer-specified operating temperature envelope and proposes a machine learning-based approach for modelling the degradation. Experimental results show that end of life, defined as a 30% decrease in capacitance, occurs at 1,000 cycles when the environmental temperature exceeds the maximum operating temperature by 30%. The life-cycle test data is then used as an input to a Gaussian process regression (GPR) algorithm to predict the capacitance fade trend. The GPR is validated on a total of nine commercial cells from two different manufacturers, achieving an average root mean squared percent error of less than 2% and a mean calibration score of 93% when referenced to a 95% confidence interval. The model can be utilized to determine the EC degradation rate at a range of operating temperature values.

42 ENGINEERING↗

Model-Based Reconstruction for Collimated Beam Ultrasound Systems

Collimated beam ultrasound systems are a novel technology for imaging inside multi-layered structures such as geothermal wells. Such systems include a transmitter and multiple receivers to capture reflected signals. Common algorithms for ultrasound reconstruction use delay-and-sum (DAS) approaches; these have low computational complexity but produce inaccurate images in the presence of complex structures and specialized geometries such as collimated beams.In this paper, we propose a multi-layer, ultrasonic, model-based iterative reconstruction algorithm designed for collimated beam systems. We introduce a physics-based forward model to accurately ac-count for the propagation of a collimated ultrasonic beam in multi-layer media and describe an efficient implementation using binary search. We model direct arrival signals, detector noise, and a spatially varying image prior, then cast the reconstruction as a maximum a posteriori estimation problem. Using simulated and experimental data we obtain significantly fewer artifacts relative to DAS while running in near real time using commodity compute resources.

Alanazi, Abdulrahman M.↗