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,009 records · Page 56

Determination of the regeneration energy of direct air capture solvents/sorbents using calorimetric methods

Regeneration energy and recyclability are important parameters that dictate the viability of a CO 2 direct air capture (DAC) process. In this study, differential scanning calorimetry (DSC) is used in conjunction with thermogravimetric analysis (TGA) and Fourier Transform Infrared Spectroscopy (FTIR) to evaluate the regeneration energy required for CO 2 -loaded 3 M aqueous potassium sarcosinate (K-SAR) solvent and crystalline methyl-glyoxal-bisiminoguanidine (MGBIG) sorbent. Based on calorimetric measurements, for aqueous K-SAR, the sensible heat amounts to 1.5 GJ/tCO 2 and the enthalpy of desorption is estimated at 3.68 GJ/tCO 2 . The total energy of regeneration will also include the enthalpy of solvent vaporization, which will depend on the thermodynamic conditions and the level of regeneration. For MGBIG, the total regeneration energy required was measured at approximately 7.0 GJ/tCO 2 . Further, FTIR measurements were used to observe CO 2 and H 2 O release during the thermal regeneration process. For K-SAR, CO 2 release is seen to take place at temperatures greater than 80 °C and, for MGBIG, the CO 2 release takes place at temperatures greater than 100 °C. The moderate temperatures required for K-SAR and MGBIG regenerations compared to other DAC solvents/sorbents suggest that low-cost sources of heat, such as geothermal heat, could potentially be used for regeneration to drive down the cost of direct air capture.

42 ENGINEERING↗

A Novel Zwitterionic Chromatography Approach to Separate Lithium from Unconventional Resources

Lithium (Li) is a key element for clean energy technologies, and, accordingly, the global lithium demand has been increasing rapidly. Therefore, to meet the Li demand and maintain supply chain stability, it is critical to develop efficient lithium extraction technologies that allow exploitation of unconventional lithium resources, such as geothermal brines and inland brine streams. However, the recovery of Li from these resources is challenging due to low Li concentration, low ratios of Li/Na, Li/Mg, or Li/Ca, and complex feed compositions. To address this, we introduced a new Direct Lithium Extraction (DLE) process using Zwitterionic Chromatography (ZIC) to separate Li from other salts. Since salts are partitioned on ZIC under water elution, no reagent chemicals are needed, and the Li separation is not limited by the adsorption capacity. We prepared 13 different zwitterionic (ZI) resins to investigate the salt retention on various ZI groups and then screened out promising sorbents for efficient Li separation. It was found that salt retention was synergistically affected by the pore size and ZI configurations. Using carboxybetaine (QAC3CA) sorbents, multicomponent separations showed that Li can be partitioned from divalent salts or Na with selectivities of 1.8 or 1.9, respectively. Although the selectivity is relatively low, in real brine tests, Li was separated from Ca and Mg with 79.2 % yield, showing the potential for a continuous process to achieve high productivity and high yield. Simulation studies suggest the salt elution mechanism is related to the hydration reaction energy and the effective hydrated radius of cations.

critical minerals↗

Nitrogen fixation in Pampinifervens , a new species-rich genus of Aquificaceae that inhabits a wide pH range in terrestrial hot springs

Members of the Aquificota play important roles in hydrogen, sulfur, and carbon cycles in geothermal systems as ubiquitous and abundant chemolithoautotrophs; however, their roles in the nitrogen cycle are poorly defined. Here, we show that Aquificaceae isolate T-2 T from Tengchong, China, fixes nitrogen based on growth with dinitrogen as the sole nitrogen source and incorporation of 15 N atoms when grown with 15 N 2 . We further show evidence that suggests its nitrogen fixation (nif) genes are transcribed in situ in a nearby spring. We propose the name Pampinifervens diazotrophicum gen. et sp. nov., with strain T-2 T (= JCM 35475 T = DSM 116324 T ) as the nomenclatural type. Its close relative designated Pampinifervens florentissimum sp. nov. T-8 T (= CGMC 1.5214 T = JCM 33569 T ) does not fix nitrogen but respires nitrate. Analysis of 104 metagenome-assembled genomes from 84 metagenomes allowed us to circumscribe the genus Pampinifervens and identify at least 16 sympatric Pampinifervens species that dominate the Aquificaceae populations in diverse hot springs in southwest China, with distinct species globally. Analysis of these Pampinifervens genomes revealed variable presence of nif and denitrification genes and frequent gene loss and horizontal transfer. Genes for pH homeostasis were also present in Pampinifervens genomes, with a proliferation in species from Tengchong, consistent with the wide pH growth range of both Pampinifervens isolates and the prevalence of Pampinifervens in both alkaline and acidic springs. The dynamic evolution of nif and denitrification genes, along with pH homeostasis genes, may facilitate the diversification of Pampinifervens into diverse springs with different nitrogen availability.

Aquificae↗

Interfacial and Confinement-Mediated Organization of Gas Hydrates, Water, Organic Fluids, and Nanoparticles for the Utilization of Subsurface Energy and Geological Resources

Harnessing the subsurface geologic environments in an efficient and environmentally sustainable manner is challenged by uncertainties associated with predicting the fate of fluids and sustaining porosity and permeability in subsurface geologic environments. Some of these uncertainties arise from confined and interfacially induced structures of fluids in subsurface geologic environments. The formation of gas hydrates, phase transitions of confined fluids, assembly and deposition of heavy hydrocarbons, and the agglomeration and fate of nanoparticles in confined environments are summarized in this review. Nanoscale confinement contributes to anisotropic structures and dynamics of fluids, which is the basis for anomalous phase transition thermodynamics, reactivity, transport, and geomechanical behavior. In this review, we discuss the structures of confined fluids and deviation in observed properties from bulk fluids. The factors influencing the structures of confined fluids can be generally divided into two groups: (a) pore characteristics including pore size, pore surface chemistry, and pore geometry and (b) confined fluid/solid characteristics such as molecular structure, concentrations, charges, pore filling, and presence of additives. Scientific advancements and knowledge gaps in our understanding of the structures of confined fluids and the associated differences in observed properties compared to bulk fluids are discussed. Here, the phenomena discussed in this review are of particular relevance to our efforts in harnessing the subsurface environments for a low carbon future by increasing the utilization of geothermal energy, using CO 2 as a working fluid, and storing CO 2 in subsurface geologic environments.

42 ENGINEERING↗

Quantification of the Impact of Acidified Brine on Fracture-Matrix Transport in a Naturally Fractured Shale Using in Situ Imaging and Modeling

Understanding flow, transport, chemical reactions, and hydromechanical processes in fractured geologic materials is key for optimizing a range of subsurface processes including carbon dioxide and hydrogen storage, unconventional energy resource extraction, and geothermal energy recovery. Flow and transport processes in naturally fractured shale rocks have been challenging to characterize due to experimental complexity and the multiscale nature of quantifying continuum scale descriptions of mass exchange between micrometer-scale fractures and nanometer-scale pores. In this study, we use positron emission tomography (PET) to image the transport of a conservative tracer in a naturally fractured Wolfcamp shale core before and after the core was exposed to low pH brine conditions. Image-based experimental observations are interpreted by fitting an analytical transport model to fracture-containing voxels in the core. Results of this analysis indicate subtle increases in matrix diffusivity and a slightly more uniform fracture velocity distribution following exposure to low pH conditions. These observations are compared with a multicomponent one-dimensional reactive transport model that indicates the capacity for a 10% increase in porosity at the fracture-matrix interface as a result of the low pH brine exposure. This porosity change is the result of the dissolution of carbonate minerals in the shale matrix to low pH conditions. Further, this image-based workflow represents a new approach for quantifying spatially resolved fracture-matrix transport processes and provides a foundation for future work to better understand the role of coupled transport, reaction, and mechanical processes in naturally fractured rocks.

58 GEOSCIENCES↗

Integrated Carbon Capture and Storage in Hydrogen Production: A Combined Techno-Economic and Life Cycle Assessment

This paper presents a coupled techno-economic and life cycle assessment of “blue” hydrogen to be produced at a hydrogen facility through steam methane reforming (SMR) equipped with carbon capture and storage (CCS). Blue hydrogen was modeled in ChemCAD, while an integrated asset model represented the carbon capture and storage chain. An unabated carbon dioxide (CO 2 ) configuration release 11.99 kgCO 2 -eq/ kgH 2 . Capturing ≥95% of the CO 2 stream lowers the carbon footprint to 6.59 kgCO 2 -eq/kgH2 but raises the levelized cost of hydrogen (LCOH) from $\$$1.82/kgH 2 (no CO 2 capture) to $\$$3.22/kgH 2 ; the U.S. 45Q tax credit reduces it to $\$$2.59/kgH 2 . Incorporating CCS reduces the levelized net present value from $\$$0.87/kgH 2 to $\$$0.74/kgH 2 , owing to additional capture, transport, and storage costs. Supplying SMR with low-carbon electricity, especially nuclear, wind, or hydro, delivers the lowest carbon footprint relative to geothermal or grid mixes. Sensitivity analysis identifies that hydrogen sales price, internal rate of return, and CCS cost as the strongest economic levers, while electricity demand dominates residual lifecycle emissions. The results underscore a clear trade-off; substantial CO 2 reductions are achievable, but only with higher production costs, making supportive policy instruments, access to clean power, and robust hydrogen markets essential for large-scale deployment of blue hydrogen.

carbon capture↗

Electrically Powered High-Salinity Brine Separation Using Dimethyl Ether

Dewatering highly saline aqueous streams, from mining and geothermal leachates to industrial wastewater, is essential for effective resource recovery and safe disposal. Membraneless water extraction (MWE) uses a low-polarity solvent to separate water from concentrated aqueous solutions. In this study, we design a new MWE that uses dimethyl ether (DME) to selectively extract water from high-salinity brines, leveraging the volatility of DME to achieve rapid solvent recovery. By separating water and dissolved salts at a liquid–liquid interface, MWE minimizes the deleterious effects of scaling on vulnerable membrane and heat exchanger surfaces, reducing the need for extensive pretreatment and expensive materials. We begin by developing a computational framework for a multistage counterflow liquid–liquid contactor, which extracts water into DME, coupled with a multistage solvent regenerator that uses vapor compression to efficiently separate the desalinated water from the DME extractant. Excess Gibbs free energy and equation of state frameworks are used to model fluid phase equilibria in water–DME–sodium chloride (NaCl) mixtures, with interaction parameters estimated from experimental data. Incorporating equilibrium calculations into a system-scale computational model, we examine the performance of MWE using DME for the first time. Our analysis demonstrates that MWE can concentrate seawater desalination brine (>1.0 mol NaCl kg –1 ) to zero-liquid discharge salinities, with an energy consumption of under 50 kW h per m 3 of water extracted with a solvent recovery ratio greater than 99.9%. We highlight the importance of staging the vapor compression process to simultaneously minimize energy consumption while enabling brine concentration and product water solvent contamination. Finally, the thermodynamic framework developed here allows for the robust evaluation of new MWE solvents and systems for critical brine concentration and fractional precipitation applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of Anions on the Delithiation of [Li–Al] Layered Double Hydroxides: Thermodynamic Insights

The utilization of lithium aluminum [Li–Al] layered double hydroxides (LDHs) is explored for direct lithium extraction (DLE) from geothermal brines and minerals. Following extraction, the Li + ions need to be removed from the LDH structure and converted into LiOH or Li 2 CO 3 products, making them suitable for battery applications. The research investigates the delithiation of [Li–Al-X] LDHs (where X = Cl – , OH – , and SO 4 2– ), which were synthesized and dried under different conditions. Herein the study aims to understand how the choice of anions and drying conditions affects the delithiation process. To determine the stability of these [Li–Al] LDHs, high-temperature oxide melt solution calorimetry is employed. The results reveal that the stability of LDHs varies significantly based on postdrying conditions, anion choice, and the water content in the interlayer. The order of stability, as indicated by the enthalpy of formation values, is found to be OH-LDH-O > OH-LDH-A > SO 4 -LDH-O > SO4-LDH-A > Cl-LDH-O > Cl-LDH-A. This trend is attributed to the interactions among the interlayer species and the metal hydroxide layers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Carbonate-Metal Reactions in the Lower Mantle

Carbonates are important carbon-bearing phases in the mantle. While their role in upper mantle petrologic processes has been well studied, their effect on phase relations, melting, and transport properties in the lower mantle is less understood. The stability of carbonates in the mantle depends on a host of factors, including pressure, temperature, oxygen fugacity, and reactions with surrounding mantle phases. To understand the stability of carbonates in the presence of metal in the lower mantle, carbonate-metal reaction experiments on the Fe–Si–Ca–Mg–C–O system were conducted up to 124 GPa and 3200 K. We find that carbonates react with iron alloys to form silicates, iron carbides, and oxides. However, the temperature at which these reactions occur increases with pressure, indicating that along a geotherm in the lowermost mantle carbonates are the stable carbon-bearing phase. Carbon is found to be less siderophilic at high-pressure compared to silicon.

58 GEOSCIENCES↗

Membranes for Lithium Recovery From Conventional and Unconventional Sources

Lithium has been deemed a critical mineral of national importance that finds uses in a wide range of applications, and its demand has been rising significantly in recent years. The urgency of meeting this demand requires lithium extraction from various aqueous sources such as continental brines, geothermal brines, seawater, produced water, and battery waste. While direct lithium extraction (DLE) technologies such as adsorption, ion exchange, and solvent extraction have emerged as possible solutions, membrane technologies are also being investigated for various sources and at different stages of the recovery process. Here, we analyze the application of membranes for pretreatment of lithium source waters, bring management, lithium/magnesium separation, lithium/sodium separation, and lithium hydroxide conversion, and evaluate performance metrics for critical lithium separations from the literature. We explore the potential of membranes at every stage of the recovery process and describe their current status and future prospects. We describe hypothetical process trains with integrated membrane technologies for each source type and address their feasibility and challenges. The potential energy and water impacts of membrane-integrated and conventional DLE processes are also critically considered alongside performance and selectivity metrics, and this is illustrated using examples and calculated from published technical reports. This paper thus provides a comprehensive overview of the application of membranes along every stage of the lithium recovery process, emphasizing the versatility and potential of membrane technologies for critical mineral recovery.

36 MATERIALS SCIENCE↗

Extreme Temperature Cryptography Based On Nitrogen-Incorporated Ultrananocrystalline Diamond

Physical entropy sources that remain stable under extreme temperatures are essential for cryptography in emerging technological frontiers in deep space exploration, geothermal energy harvesting, and nuclear energy. However, conventional semiconductor platforms fail to generate stable and reliable cryptographic keys above 200 degrees C due to performance degradation. Here, we report a diamond-based cryptographic primitive that exploits the defect-rich sp 2 -bonded grain boundary network in nitrogen-incorporated ultrananocrystalline diamond (n-UNCD) film as a robust entropy source to generate cryptographic keys that remain operationally stable even after enduring extreme temperatures of 700 degrees C for 54 h while also surviving thermal cycling between room temperature and 700 degrees C for 48 h. The strength of the generated keys is assessed through several cryptographic metrics such as bit uniformity, entropy, hamming distances, and correlation coefficients, all of which are found to be near their respective ideal values. Moreover, the generated keys pass the NIST SP 800 and SP 800-90B tests and are also resilient to supply bias variations and a regression-based machine learning attack model based on the Fourier series. The robustness of the keys is attributed to the better thermal stability and chemical inertness of the n-UNCD film. This is supported by high-resolution energy-dispersive X-ray spectroscopy (EDS), which shows no significant lateral diffusion of metal atoms into the n-UNCD layer, and by Raman spectroscopy, which reveals no significant changes in the bonding configuration of the n-UNCD structure. Our findings highlight the remarkable potential of n-UNCD film for extreme environment cryptography by expanding the operational limits of conventional hardware security platforms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Integration of the Biot–Gassmann Fluid Substitution Method and Machine Learning-Based Velocity–Stress Relationship for Estimating In Situ Stresses

Recent advancements have shown that in situ stresses can be reliably estimated through an integrated machine/deep learning (ML/DL)-based framework, which relies on models trained and validated using true triaxial ultrasonic velocity (TUV) experimental data that involve measurements of ultrasonic velocity in saturated rocks under varying stress configurations. However, when the goal is to interpret lower frequency measurements, it may be more appropriate to run experiments on dry rocks and then obtain Biot–Gassmann-derived equivalent saturated velocities (low-frequency approximation) and employ these quantities for training ML/DL models to predict in situ stress. Whether the dispersion effect of frequency on the velocity–stress relationship substantially impacts in situ stress prediction is an important and unresolved question. This work presents an enhancement of ML/DL-based workflow by training and implementing ML/DL models using equivalent saturated acoustic velocities (low-frequency) obtained by applying Biot–Gassmann fluid substitution on the ultrasonic velocities of dry cores. The models were trained on TUV data sets derived from three subsurface cores extracted from the geothermal well 16B(78)-32 at the Utah FORGE site. Each core was subjected to 75 unique stress configurations for velocity measurement in the dry state. The ML/DL trained on the TUV data set with equivalent saturated velocities demonstrated promising performance to predict in situ stress in subsurface geological rocks using velocity–stress relationships with R 2 of 0.86, 0.971, and 0.975 and root mean squared error (RMSE) of 2.59, 1.92, and 1.80 for validation/testing phases of vertical, minimum horizontal, and maximum horizontal stress models, respectively. Additionally, interpretation and explanation by Shapley additive explanations (SHAP) analysis further improved scientific validation and model reliability for estimating in situ stresses.

colloids↗

Composition of Magma and Characteristics of the Hydrothermal System of Newberry Volcano, Oregon, From Magnetotellurics

Abstract We use 3‐D magnetotellurics to improve our understanding of the structure and magma composition of Newberry Volcano in Oregon, USA. Newberry is a broad shield volcano with a summit caldera and strongly bimodal magmatism. Newberry has long been the subject of geothermal exploration research, but that work has focused on the volcano's west flank, leaving the caldera largely unstudied with geophysical methods until recently. Our modeling shows a relatively resistive magma reservoir of approximately 50 Ωm. Our work builds upon recent seismic models and petrological analysis to interpret Newberry's magma reservoir as a dry rhyolite with 8–11% partial melt, which matches the seismically determined melt fraction. Finding the conditions within the magma reservoir that allow the resistivity and seismic analysis to come to the same melt fraction helps us narrow down the magma temperature and composition. From this we infer a dry rhyolitic magma at 850 °C. We also image a prominent vertical conductive anomaly along the south rim below the vent that produced the most recent eruption. The anomaly extends from magma reservoir depths of 3 km to 1 km below the surface where it disperses into the caldera fill. We interpret this as the main conduit for magmatic fluids to reach Newberry's hydrothermal system. Other features that our model shows include higher conductivity along the caldera rim and a resistive older pluton on the west flank.

Bowles‐Martinez, E.↗

Collapse of Reacted Fracture Surface Decreases Permeability and Frictional Strength

Abstract Geochemical and geomechanical perturbations of the subsurface caused by the injection of fluids present risks of leakage and seismicity. This study investigated how acidic fluid flow affects hydraulic and frictional properties of fractures using experiments with 3.8‐cm‐long specimens of Eagle Ford shale, a laminated shale with carbonate‐rich strata. In low‐pressure flow cells, one set of samples was exposed to acidic brine and another set was exposed to neutral brine. X‐ray computed tomography and energy‐dispersive X‐ray spectroscopy revealed that samples exposed to acidic brine were calcite‐depleted and had developed a porous altered layer, while the other set showed no evidence of alteration. After reaction, samples were compressed and sheared in a triaxial cell that supplied normal stress and differential pore pressure at prescribed sliding velocities, independently measuring friction and permeability. During the initial compression, the porous altered layer collapsed into fine particles that filled the fracture. This effectively impeded flow and sealed the fracture, resulting in fracture permeability to decrease 1 to 2 orders of magnitude relative to the unaltered fractures. This is a favorable outcome in subsurface applications where the goal is to reduce leakage risks. However, during shear the reacted fracture had lower frictional strength because the fine‐grained particles in the collapsed layer prevented the formation of interlocking microasperities. Therefore, coupled geochemical and geomechanical processes that could favorably seal fractures could also increase the likelihood of induced seismicity. These findings have important implications for geological carbon sequestration, pressurized fluid energy storage, geothermal energy, and other subsurface technologies.

58 GEOSCIENCES↗

Time–series Analysis of Volume Change at Brady Hot Springs, Nevada, USA using Geodetic Data from 2003 – 2018

Brady Hot Springs geothermal field has exhibited subsidence, as measured by interferometric synthetic aperture radar (InSAR). Previous studies have examined both the temporal evolution of the deformation from 2004 through 2016 and the spatial extent of the deformation, directly relating the observed subsidence to volumetric changes below the surface. We extend the modeling at Brady to analyze a data set of interferometric pairs spanning from the end of 2003 through 2018. We examine spatial and temporal trends in the observed deformation by time–series analysis of each of the 1656 cubic voxels in a parameterized elastic dislocation model to identify areas where the subsurface volume changes as a function of time. Joint time–series analysis of Global Positioning System and InSAR pairs confirm significant changes in rates of volume change during time intervals when well operations were varied. Here, the rate of subsidence increases with increased injection, consistent with the identification of thermal contraction of the rock matrix as the dominant driving mechanism. Conversely, the modeled volume increases when pumping ceases, suggesting thermal expansion of the rock matrix.

58 GEOSCIENCES↗

Microbial Community Composition in Deep‐Subsurface Reservoir Fluids Reveals Natural Interwell Connectivity

Abstract The identification of natural fractures and the wells they connect is crucial for the development of geological reservoirs because it may have an important impact on reservoir model construction and hydraulic fracture propagation. In this study we investigated the use of a novel data source, the microbial community composition in the reservoir formation fluids, for identification of interwell connectivity caused by natural fractures. We verified this concept at a newly developed mesoscale enhanced geothermal system (EGS) field testbed located 4,850 ft (1,478 m) beneath ground surface at the Sanford Underground Research Facility in Lead, SD. Fluids produced at or near the EGS testbed were sampled and subjected to high‐throughput 16S rRNA gene amplicon sequencing to analyze the microbial community profile therein. Despite the typically substantial heterogeneity across the community profiles of samples spatially distributed (10 m to 1.9 km apart) throughout the site, samples from two wells at the EGS testbed showed highly similar microbial community composition, suggesting the two wells intersected the same natural fracture. This evidence of natural connectivity between the two wells at the EGS testbed was later corroborated by core log analysis and sewer camera surveys into the boreholes. Besides the field case described in this study, microbial community composition as a reservoir diagnostic tool would be applicable in a much broader context such as unconventional hydrocarbon exploitation, groundwater reservoir characterization, and environmental remediation, adding valuable “hard” data capable of pinpointing the origins of fluids unambiguously.

Zhang, Yuran↗

Creation of a Mixed-Mode Fracture Network at Mesoscale Through Hydraulic Fracturing and Shear Stimulation

Enhanced Geothermal Systems could provide a substantial contribution to the global energy demand if their implementation could overcome inherent challenges. Examples are insufficient created permeability, early thermal breakthrough, and unacceptable induced seismicity. In this study we report on the seismic response of a mesoscale hydraulic fracturing experiment performed at 1.5-km depth at the Sanford Underground Research Facility. We have measured the seismic activity by utilizing a 100-kHz, continuous seismic monitoring system deployed in six 60-m length monitoring boreholes surrounding the experimental domain in 3-D. The achieved location uncertainty was on the order of 1 m and limited by the signal-to-noise ratio of detected events. These uncertainties were corroborated by detections of fracture intersections at the monitoring boreholes. Three intervals of the dedicated injection borehole were hydraulically stimulated by water injection at pressures up to 33 MPa and flow rates up to 5 L/min. We located 1,933 seismic events during several injection periods. The recorded seismicity delineates a complex fracture network comprised of multistrand hydraulic fractures and shear-reactivated, preexisting planes of weakness that grew unilaterally from the point of initiation. We find that heterogeneity of stress dictates the seismic outcome of hydraulic stimulations, even when relying on theoretically well-behaved hydraulic fractures. Once hydraulic fractures intersected boreholes, the boreholes acted as a pressure relief and fracture propagation ceased. In order to create an efficient subsurface heat exchanger, production boreholes should not be drilled before the end of hydraulic stimulations.

58 GEOSCIENCES↗

Equation of State Measurements on Iron Near the Melting Curve at Planetary Core Conditions by Shock and Ramp Compressions

The outer core of the Earth is composed primarily of liquid iron, and the inner core boundary is governed by the intersection of the melt line and the geotherm. While there are many studies on the thermodynamic equation of state for solid iron, the equation of state of liquid iron is relatively unexplored. In this work, we use dynamic compression to diagnose the high-pressure liquid equation of state of iron by utilizing the shock-ramp capability at Sandia National Laboratories’ Z-Machine. This technique enables measurements of material states off the Hugoniot by initially shocking samples and subsequently driving a further, shockless compression. Planetary studies benefit greatly from isentropic, off-Hugoniot experiments since they can cover pressure-temperature (P-T) conditions that are close to adiabatic profiles found in planetary interiors. We used this method to drive iron to P-T conditions similar to those of the Earth’s outer-inner core boundary, along an elevated-temperature isentrope in the liquid from 275 GPa to 400 GPa. We derive the equation of state using a hybrid backward integration – forward Lagrangian technique on particle velocity traces to determine the pressure-density history of the sample. Our results are in excellent agreement with SESAME 92141, a previously published equation of state table. With our data and previous experimental data on liquid iron we provide new information on the iron melting line and derive new parameters for a Vinet-based equation of state. The table and our parameterized equation of state are applied to provide an updated means of modeling the pressure, mass, and density of liquid iron cores in exoplanetary interiors.

58 GEOSCIENCES↗