Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Pressure Core”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

High Flux Isotope Reactor (HFIR)

Oak Ridge National Laboratory’s High Flux Isotope Reactor (HFIR) is the highest-flux reactor-based neutron source in the United States. HFIR’s intense neutron flux and state-of-the-art facilities result in world-class capabilities, including neutron scattering, radioisotope production, materials irradiation, and neutron activation analysis. HFIR’s multimission capabilities are attributed to its versatile high-power-density core design consisting of a series of concentric regions, including a flux trap target region, an inner fuel element, an outer fuel element, a control element region, and a beryllium reflector. The pressurized, light-water-cooled research reactor operates at 85 MW and is fueled by 9.4 kg 235U.

Chandler, David↗
DOE OSTI2022

Enhancing reactant selectivity for Ni/Mg reforming catalysts using silicalite-1 shells: A modeling study

The ability of a zeolite shell to enhance the selective conversion of hydrocarbons through diffusional limitations was investigated using a multi-scale model of a fixed-bed reactor. The impact of shell thickness and molecule/pore size on the catalytic performance of silicalite-1 zeolite encapsulated nickel catalyst pellets for steam reforming of C1-C7 hydrocarbons is reported. A reaction–diffusion model using kinetic expressions established in literature was employed. The model was verified through comparison with reported experimental results for steam reforming data over a temperature range of 748 – 1113 K and pressure of 1 – 10 bar. Comparisons are also made against experimental data for steam reforming in the presence of a zeolite shell. Evaluation of the Weisz-Prater criterion for both the core and encapsulated catalyst confirmed mass transfer limitation induced by the utilization of a zeolite shell. The model was used to suggest an optimal thickness that balances diffusional limitations imposed by the zeolite layer on methane versus that of the heavier hydrocarbons. The optimum thickness varied as a function of hydrocarbon size and shape which determined the diffusion rates. For toluene and heptane, a 50 nm thick shell was sufficient to wholly prevent reaction. Hydrocarbons like propane and butane required a shell 7.5 and 5 μm thick. Increasing the gas-hourly-space-velocity from 10,000 to 60,000 h –1 caused a decrease in the optimum shell thickness. Furthermore, this approach can be modified for application to other mixed hydrocarbon systems to predict optimal catalyst design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Utah FORGE: Documentation on Discrete Fracture Network and Fracture Propagation Modelling

This dataset includes reports and a slide presentation on discrete fracture network (DFN) generation and hydraulic fracture modeling at the Utah FORGE site. It details the characterization of natural fractures using well log and core data, as well as stochastic modeling techniques. The reports describe simulations of hydraulic fracture propagation, fluid-mechanical interactions, and induced microseismicity. The dataset also includes history-matching of net pressure and analyses of fracture growth in naturally fractured geothermal reservoirs. The slides summarize key findings and future research directions.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Adrian Portales Coop Internship Report

Using the thesis of W.R. Nolan (cite) as a guide, Cobalt Iron (CoFe) powders were reacted with 0.1 wt.% and 0.2 wt.% phosphoric acid in a 20:1 ratio of acetone to phosphoric acid. The powders were then dried at room temperature. The resulting phosphate coated CoFe was then mixed with 0.75 wt.% of the lubricant N,N' ethylene bis-stearamide (trade name: Acrawax C) and hot pressed to form a consolidated soft magnetic material referred to as CoFeP. With an avenue of synthesis for CoFeP determined, a proper amount of stock was synthesized for continuous “brick” production. While under current optimization, these 1x1 mm magnetic bricks will ultimately be placed and secured along the inside wall of each MK Magnetics transformer core by an appropriate CoFeP dispersed epoxy. As of now each brick has been produced though a pressing and annealing process via square 2x2 cm die. Before a brick is made a pressure calculation is run to ensure the dies maximum operating pressure is not exceeded. Figure 1. ensures the user’s safety by showing that the tons-on-ram required for a 2x2 cm square die to reach 760 MPa is below the point of die failure.

36 MATERIALS SCIENCE↗
DOE OSTI2022

FY22 Status Report on the ART-GCR CMVB and CNWG International Collaborations

This work presents the numerical model for the high temperature test reactor (HTTR) loss of forced cooling (LOFC) experiment with the INL codes Griffin, BISON, and RELAP-7, based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Promising results were obtained, with the overall behavior of the reactor successfully captured. Changes in the heat transfer coupling, as compared to the fiscal year (FY)-21 model, enabled drastic improvement of the steady-state solution, both in terms of computational time and global energy discrepancies. The former was reduced by a factor of roughly 60, whereas the latter decreased from 11% to less than 2%. Furthermore, the discrepancy in the steady-state multiplication factor was improved from +2,300 and +2,900 pcm (for the 30 and 9 MW cases) to -700 and +1,200 pcm, respectively, and now falls well within the large measurement uncertainties stemming from graphite impurities. Validation of the Monte Carlo model used to generate cross sections was also performed against available measurements. Though significant, the discrepancies remain acceptable overall in light of the large uncertainty. Numerous improvements are still needed to better compare with the experiment involving the 9 MW case and to instill greater confidence in the model’s ability to accurately predict the 30 MW behavior. Specifically, the power levels predicted by the 9 MW transient simulation following re-criticality remain low, pointing to an underestimation of the passive cooling of the core. A key aspect of future work will be to better understand the flow pattern during the LOFC event, particularly to determine if natural or forced convection is occurring inside the reactor pressure vessel (RPV). More generally, additional validation data would be immensely useful for further enhancing the numerical model and better matching the experiments. In addition, a more sophisticated thermal-hydraulics model that simulates all the channels as a single system model should be considered to take into account the rest of the primary loop. Finally, even if the results are in better agreement with the experiments, sensitivity analysis and uncertainty quantification will be necessary to evaluate the model uncertainty.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Effects of Interactions Between Produced Formation Fluid and Rock Matrix on Pore Structure of Caney Shale, Southern Oklahoma

ABSTRACT: Rock-fluid interactions change properties of shales during exploitation. To investigate effects of rock-fluid interactions on pore structure of shales matrix after hydraulic fracturing, powder samples from two late Mississippian Caney Shale cores in the Ardmore Basin, southern Oklahoma, were used to react with formation produced fluid from the field in the batch reactor analysis. X-ray diffraction for mineralogy and Low-pressure nitrogen adsorption isotherms for pore structure were measured for original, after-7days, and after-30days samples. Results show that the samples consist mainly of quartz, followed by clay minerals, carbonates, and feldspar. The pore sizes of micropore (<2 nm) and mesopore (2-50 nm) increase 14%-233% due to dissolution of pyrite, feldspar, and carbonates after 7 days. Due to the transformation from smectite to illite and the increase of pore size, the specific surface area (SSA) decreases after 7-days interactions. After 30-days interactions, the micropore volume slightly increases and the mesopore and macropore volume decreases. Due to the decrease of pore size, the SSA of 30-days reacted samples increases correspondingly and is lower (for the clay-rich sample) or higher (for the calcareous sample) than that of the unreacted samples. Findings improve our understanding of dynamic alteration of shale properties during production. 1. INTRODUCTION Energy demand will continuously grow owing to the increasing global population as well as energy consumption (EIA, 2023). On the other hand, shale gas and oil reshaped the energy market in the United States, enabling the United States to become a net-export of natural gas country in 2017 (EIA, 2023). However, shale reservoirs are challenging tight formations that are still poorly understood in the extraction and production of hydrocarbons (Ross and Bustin, 2009; Curtis et al., 2012; Xiong et al., 2015, 2021a; Li Y. et al., 2016; Gong et al., 2019a; Benge et al., 2021; Awejori et al., 2022; Huang et al., 2022). One of the most challenging topics is the rock-fluid interactions post hydraulic fracturing and its subsequent impacts on the pore structures of fractured formation matrix.

Xiong, Fengyang↗
DOE OSTI2022

CTF Improved Drag Model and Flow Regime Transition Criteria

The demand for accurate prediction of two-phase flow behavior in a boiling water reactor (BWR) requires a comprehensive understanding of flow regime, void fraction, heat transfer, and pressure drop. The CTF subchannel code, which is used for the Thermal/Hydraulic (T/H) solution in the Consortium for Advanced Simulation of Light Water Reactors (CASL)-developed Virtual Environment for Reactor Application (VERA) core simulator, is being further developed for BWR applications. In support of this goal, the present work highlights some of the two-phase closure model developments towards improving the CTF void fraction prediction, especially for subcooled boiling. The drift-flux approach has been well-developed for upward dispersed two-phase flows and proven to be accurate in predicting void fraction in bubbly and slug flow regimes. In this work, these kinematic constitutive relations for the drift-flux velocity have been implemented into CTF to describe the interfacial drag of bubbly flow as an alternative to the existing model for better void fraction prediction. The success of these constitutive relations also relies on a good flow regime map that accounts for flow conditions and channel geometry. A more reliable flow regime transition criteria that account for the flow condition has also been implemented in this study for modeling the flow regime transition criteria. The newly implemented models are shown to give improved void fraction predictions in comparison to experimental data.

Hizoum, Belgacem↗
DOE OSTI2025

Thermally induced structural competitiveness and metastability of body-centered-cubic iron under nonequilibrium conditions

The structure and stability of iron near melting at multi-megabar pressures are of significant interest in high pressure physics and earth and planetary sciences. While the body-centered cubic (BCC) phase is generally recognized as unstable at lower temperatures, its stability relative to the hexagonal close-packed (HCP) phase at high temperatures (approximately 0.5 eV) in the Earth’s inner core (IC) remains a topic of ongoing theoretical and experimental debate. Motivated by a series of studies attempting to understand the state of the IC and resulting in contradictory findings, we aim to explore the physical factors that may contribute to this confusion in a broader context. Our ab initio calculations show a significant drop in energy, the emergence of a plateau and a local minimum in the potential energy surface, and stabilization of all phonon modes of BCC iron at elevated electron temperatures (>1–1.5 eV). These effects increase the competition among the BCC, HCP, and the face-centered cubic (FCC) phases and lead to the metastability of the BCC structure. Furthermore, at pressure-temperature conditions relevant to the IC, we can expect substantial coexistence of FCC and HCP iron due to their similar free energies and low transition barrier (∼0.2 eV/atom), and the thermodynamic stability of BCC iron is enhanced by its substantial lattice vibration entropy. In conclusion, these findings provide a clear theoretical framework for understanding iron phase relations and solidification processes, both experimentally and in the IC.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Risk-Informed Safety Analysis for Accident Tolerant Fuels

Accident Tolerant Fuels (ATF) are being tested by different nuclear vendors and research organizations and their introduction in the US light water reactors fleet is planned for the second half of the 2020`s. In the framework of the US-DOE Light-Water Reactor Sustainability program, Risk-Informed Systems Analysis pathway, researches are being conducted at the Idaho National Laboratory for developing tools and methods that can help the industry quantifying the benefits from the adoption of ATF technology. In this paper we describe the developed risk-informed methodology including the safety analysis codes improvements, and we present some results for selected accident scenarios. The developed methodology combines INL state-of-the-art deterministic Best Estimate code RELAP5-3D, and Probabilistic Risk Analysis tools RAVEN and SAPHIRE. The analyses are performed on a three-loops pressurized water reactor, simulating Station Blackout and LB-LOCA accidents and considering near-term ATFs, or FeCrAl and Chromium-coated clads. Finally, we show how, applying our methodology, the new Core Damage Frequency can be assessed. The results indicated that the main benefits in introducing near-term ATFs concern a significant reduction in the hydrogen production during accident conditions. No significant Core Damage Frequency reduction was found.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Updates to Nek5000: RANS wall functions, NRC support, and documentation

As a key part of the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation program (NEAMS), the next generation of tools for nuclear reactor design are being developed. A primary metric for success of this program is the adoption of those tools by industry stakeholders. To help meet this metric for Nek5000, development efforts this year have focused on the expansion of the available Reynolds-averaged Navier-Stokes (RANS) turbulence models to include wall function models, continued support of the U.S. Nuclear Regulatory Commission in simulating a hydrogen mitigation benchmark experiment, and expansion of the Nek5000 documentation to enhance the code’s usability. The implementation of wall functions is intended to significantly decrease time to-solution for a wide range of problems and it has been identified as an important feature by stakeholders across the nuclear community. The initial implementation of standard and pressure corrected wall functions has been tested in Nek5000 on the classical problems of channel flow and a backward facing step. It has also been tested on a 2D molten salt fast reactor core and a T-junction. The recommended formulation of choosing y + = 30 as a boundary value combined with a Neumann-Neumann formulation of the boundary condition for k and τ produces the most consistent results compared to a wall resolved approach. The collaboration with the NRC is a continuing exercise that has been ongoing for multiple years. This year we have concluded investigations into appropriate inlet conditions, showing that fully developed turbulent conditions are adequate for the full PANDA domain. We have also concluded simulations for an unobstructed jet and are now refocusing to the original obstructed jet case. The expansion of the code documentation has also been an ongoing effort over the past few years to address feedback from the various stakeholders. The documentation now more fully covers many of the available features and closes identified SQA gaps. Finally we report on three training sessions that were offered over the past year in an effort to expand the Nek5000 user base.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Experimental Characterization of an Additively Manufactured Inconel 718 Heat Exchanger for High-Temperature Applications

This work presents the experimental results of a novel, air-to-air, additively manufactured manifold-microchannel heat exchanger with straight fins on both sides. The heat exchanger was made of Inconel 718 using a direct metal laser sintering technique. The overall core size of the heat exchanger was 94 mm × 87.6 mm × 94.4 mm, with a fin thickness of 0.220 mm on both the hot and cold sides. The heat exchanger was tested with pressurized nitrogen gas at 300 °C and 340 kPa for the hot side, while air at an ambient condition was used for the cold side. An overall heat transfer of 276 W/m2K was obtained for Reynolds number values of 132 and 79 for the cold and hot sides, respectively. A gravimetric heat transfer density (Q/m∆T) of 4.7–6.7 W/kgK and a volumetric heat transfer density (Q/V∆T) of 6.9–9.8 kW/m3K were recorded for this heat exchanger with a coefficient of performance value that varied from 42 to 52 over the operating conditions studied here. The experimental pressure drop results were within 10% of the numerical values, while the corresponding heat transfer results were within 17% of the numerical results, mainly due to imperfections in the fabrication process. Despite this penalty, the performance of the tested heat exchanger was superior to the conventional plate-fin heat exchangers: more than 60% of improvements in both gravimetric and volumetric heat transfer densities were recorded for the entire range of experimental data.

42 ENGINEERING↗
DOE OSTI2025

Morphology and property tuning in ZnO–Ni hybrid metamaterials in vertically aligned nanocomposite (VAN) form

ZnO thin films have attracted significant interest in the past decades owing to their unique wide band gap properties, piezoelectric properties, non-linearity and plasmonic properties. Recent efforts have been made in coupling ZnO with secondary phases to enhance its functionalities, such as Au–ZnO nanocomposite thin films with tunable optical and plasmonic properties. In this work, magnetic nanostructures of Ni are incorporated in ZnO thin films in a vertically aligned nanocomposite (VAN) form to couple magnetic and plasmonic response in a complex hybrid metamaterial system. Nickel (Ni) is of interest due to its ferromagnetic and plasmonic properties along with gold (Au) which is also plasmonic. Therefore, two approaches, namely, tuning of the deposition pressure and use of a ZnO–Au seeding layer have been attempted to achieve unique Ni nanostructures in addition to tuning of the microstructure. Together, both approaches demonstrate a range of microstructures such as core–shell, nanodisk, nanocup, and nanocube-like morphologies not previously attempted. Additionally, there is effective tuning of properties. Specifically, the seeding layer thickness causes hyperbolic behavior as well as redshift in the surface plasmon resonance (SPR) wavelength. The addition of the ZnO–Au seeding layer directly influences the optical properties. Plus, regardless of the different approaches, the films demonstrate magnetic anisotropy based on the composition and microstructure of the film which impacted the saturation magnetization and coercivity. This study demonstrates the potential of ZnO-based complex hybrid metamaterials with coupled electro-magneto-optical properties for integrated photonic devices.

Bhatt, Nirali A. [Purdue Univ., West Lafayette, IN↗
DOE OSTI2023

Simulation of the NuScale SMR and Investigation of the Effect of Load-Following on Component Lifetimes

The NuScale SMR has been modelled using the Virtual Environment for Reactor Applications (VERA) multiphysics environment and the results compared with the publicly reported data in the Design Certification Application (DCA). The results show an excellent agreement for the compared axial and radial power distributions, temperature coefficients of reactivity, boron and control rod worths, and fast neutron flux. This NuScale model is then used to investigate the effect of different operational modes on reactor components to determine how flexible load-following operation may affect control rod and reactor pressure vessel (RPV) lifetimes. The control rod degradation is confirmed to primarily affect the silver-indium-cadmium (AIC) rod tip. The degradation rate is observed to follow a non-linear function of core power level where the increase in degradation decreases with insertion depth. For the variation in core power levels expected with current load-following schemes, the total control rod degradation is found to be mild, at 5-10% of usable life per cycle for a reactor operating at <80% power. Nonetheless, this enables load following strategies to be confirmed and/or modified to ensure that control rods do not need to be replaced during the 60+ year life of the reactor. The RPV degradation was found to be almost directly proportional to the core power level and was not overly sensitive to flux shape perturbations. Future work is planned using these damage functions to optimize operation over multiple NuScale SMR units and develop strategies for prognostics and health management.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Simulation of metal nanoparticles growth in methane atmosphere of arc discharge: comparison to experiment

A direct current arc discharge in a methane atmosphere is a scalable and sustainable method to produce metal-carbon core–shell nanoparticles and single-walled carbon nanotubes, where a metal catalyst can be continuously supplied through evaporation of an anode made from the catalyst material. The size of catalyst particles is of critical importance as it can affect the synthesis yield and properties of nanotubes and core–shell nanoparticles. This study presents a numerical model describing the formation and growth of metal particles for the conditions representative of the arc discharge with an evaporating iron anode at near-atmospheric pressure of a methane-rich atmosphere. The model incorporates carbon adsorption to the metal surface and explains the limiting effect of carbon coverage on the size of metal nanoparticles. The predicted particle sizes are compared with experimental observations. The model also predicts higher concentrations of metal particles with the increasing partial pressure of methane.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

MHD simulations of cold bubble formation from 2/1 tearing mode during massive gas injection in a tokamak

Massive gas injection (MGI) experiments have been carried out in many tokamaks to study disruption dynamics and mitigation schemes. Two events often observed in those experiments are the excitation of the m = 2, n = 1 magnetohydrodynamic mode, and the formation of cold bubble structure in the temperature distribution before the thermal quench (TQ). Here m is the poloidal mode number, n the toroidal mode number. The physics mechanisms underlying those phenomena, however, have not been entirely clear. In this work, our recent NIMROD simulations of the MGI process in a tokamak have reproduced the main features of both events, which has allowed us to examine and establish the causal relation between them. In these simulations, the 3/1 and 2/1 islands are found to form successively after the arrival of impurity ion cold front at the corresponding q = 3 and q = 2 rational surfaces. At the interface between impurity and plasma, a local thin current sheet forms due to an enhanced local pressure gradient and moves inward following the gas cold front, this may contribute to the formation of a dominant 2/1 mode. Following the growth of the 2/1 tearing mode, the impurity penetration into the core region inside the q = 2 surface gives rise to the formation of the cold bubble temperature structure and initiates the final TQ. Here, a subdominant 1/1 mode developed earlier near the q = 1 surface alone does not cause such a cold bubble formation, however, the exact manner of the preceding impurity penetration depends on the nature of the 1/1 mode: kink-tearing or quasi-interchange.

2/1 tearing mode↗
DOE OSTI2020

Evaluation of the Four-Dimensional Ensemble-Variational Hybrid Data Assimilation with Self-Consistent Regional Background Error Covariance for Improved Hurricane Intensity Forecasts

The feasibility of a hurricane initialization framework based on the Gridpoint Statistical Interpolation (GSI)-based four-dimensional ensemble-variational (GSI-4DEnVar) hybrid data assimilation system for the Hurricane Weather Research and Forecasting model (HWRF) model is evaluated in this study. The system considers the temporal evolution of error covariances via the use of four-dimensional ensemble perturbations that are provided by high-resolution, self-consistent HWRF ensemble forecasts. It is different from the configuration of the GSI-based three-dimensional ensemble-variational (GSI-3DEnVar) hybrid data assimilation system, similar to that used in the operational HWRF, which employs background error covariances provided by coarser-resolution global ensembles from the National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS) ensemble Kalman filtering data assimilation system. In addition, our proposed initialization framework discards the empirical intensity correction in the vortex initialization package that is employed by the GSI-3DEnVar initialization framework in operational HWRF. Data assimilation and numerical simulation experiments for Hurricanes Joaquin (2015), Patricia (2015), and Matthew (2016) are conducted during their intensity changes. The impacts of two initialization frameworks on the HWRF analyses and forecasts are compared. It is found that GSI-4DEnVar leads to a reduction in track, minimum sea level pressure (MSLP), and maximum surface wind (MSW) forecast errors in all of the HWRF simulations, compared with the GSI-3DEnVar initialization framework. With assimilating high-resolution observations within the hurricane inner-core region, GSI-4DEnVar can produce the initial hurricane intensity reasonably well without the empirical vortex intensity correction. Further diagnoses with Hurricane Joaquin indicate that GSI-4DEnVar can significantly alleviate the imbalances in the initial conditions and enhance the performance of the data assimilation and subsequent hurricane intensity and precipitation forecasts

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

DCMIP2016: the tropical cyclone test case

This paper describes and analyzes the Reed–Jablonowski (RJ) tropical cyclone (TC) test case used in the 2016 Dynamical Core Model Intercomparison Project (DCMIP2016). This intermediate-complexity test case analyzes the evolution of a weak vortex into a TC in an idealized tropical environment. Reference solutions from nine general circulation models (GCMs) with identical simplified physics parameterization packages that participated in DCMIP2016 are analyzed in this study at 50 km horizontal grid spacing, with five of these models also providing solutions at 25 km grid spacing. Evolution of minimum surface pressure (MSP) and maximum 1 km azimuthally averaged wind speed (MWS), the wind–pressure relationship, radial profiles of wind speed and surface pressure, and wind composites are presented for all participating GCMs at both horizontal grid spacings. While all TCs undergo a similar evolution process, some reach significantly higher intensities than others, ultimately impacting their horizontal and vertical structures. TCs simulated at 25 km grid spacings retain these differences but reach higher intensities and are more compact than their 50 km counterparts. These results indicate that dynamical core choice is an essential factor in GCM development, and future work should be conducted to explore how specific differences within the dynamical core affect TC behavior in GCMs.

Willson, Justin L.↗
DOE OSTI2024

Liquid Structure of Iron and Iron–Nitrogen–Carbon Alloys Within the Cores of Small Terrestrial Bodies

Abstract Nitrogen has been proposed to be stored within planetary cores, but its effects on the structure and density of molten Fe–alloys have not been explored experimentally. Using energy‐dispersive X‐ray diffraction, we determined the structure of Fe–N(–C) liquids at core conditions (1–7 GPa and 1700–1900°C) within a Paris‐Edinburgh press. Variation of N up to 7 wt.% and C up to 1.5 wt.% results in near‐linear changes in Fe–Fe atom distances and structure factor with increasing light element content. We did not observe a significant pressure‐driven structural transition in Fe–N(–C) liquids. We model the expansion of the Fe–Fe bonds using a modified Birch‐Murnaghan equation of state. With this model, we demonstrate that N or C contamination could lead to an overestimation of the Fe–Fe distances of pure Fe. We observe that the incorporation of 1 wt.% N or C into Fe results in a change in Fe–Fe distances that is twice as significant as the effect of 1 GPa. By approximating the change in volume, we infer that N and C incorporated in liquid iron could contribute to the density deficit observed in the cores of terrestrial bodies.

Pease, Allison [Department of Earth and Environmen↗