Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Composite models”

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

Risk-based area of review estimation in overpressured reservoirs to support injection well storage facility permit requirements for CO 2 storage projects

This paper by the Energy & Environmental Research Center presents a workflow and modeling approach for delineating a risk-based area of review (AOR) to support a U.S. Environmental Protection Agency (EPA) Class VI permit for a carbon dioxide (CO 2 ) storage project. The approach combines semianalytical solutions for estimating formation fluid leakage through a hypothetical leaky wellbore with the results of numerical reservoir simulations to define the AOR. The modeling utilizes 1) semianalytical solutions from the peer-reviewed literature for formation fluid leakage through abandoned wellbores by Raven (1990) and Avci (1994), 2) a FORTRAN model compiled and described in Cihan et al. (2011, 2012) called ASLMA (Analytical Solution for Leakage in Multilayered Aquifers), and 3) a computational framework for estimating a risk-based AOR first proposed by Oldenburg et al. (2014, 2016). Therefore, the approach builds upon well-established research and underlying hydrogeological principles that have been upheld for nearly three decades. Moreover, the ASLMA model has been broadly applied to an array of storage projects. The work presented herein extends these earlier works using a custom wrapper written in the software environment, R (R Core Team, 2020), which was developed to perform multiple runs of the ASLMA model using given ranges for one or more input parameters. In addition, the current work simulates the pressure buildup within the storage reservoir in response to CO 2 injection using a compositional simulator to better accommodate the temporospatial evolution of pressure buildup within the storage reservoir that is more accurately modeled using a heterogeneous geologic model and a compositional simulator that accounts for the multiphase interactions. The workflow is demonstrated using a case study for a 180,000-metric-ton-per-year storage project located in the PCOR (Plains CO 2 Reduction) Partnership region. For the storage project evaluated here, under the scenario where the leaky wellbore is open to a saline aquifer (thief zone) between the overlying seal (cap rock) and the underground sources of drinking water (USDW), the risk-based AOR essentially collapses to the areal extent of the CO 2 plume in the storage reservoir because the pressure buildup in the storage reservoir beyond the CO 2 plume is insufficient to drive formation fluids up a hypothetical leaky wellbore into the USDW. However, even under the conservative assumption that the leaky wellbore is not open to a thief zone, beyond the areal extent of the CO 2 plume, the incremental leakage is less than 400 m 3 over 20 years, which represents ~0.0001% or less of the total volume of water contained within the USDW rock volume. As discussed in the text, the threshold criterion for defining the risk-based AOR is site-specific and should be informed by the results of the sensitivity analysis and available site characterization data. The approach outlined in this paper is designed to be protective of USDWs and, therefore, comply with the Safe Drinking Water Act requirements and provisions for the U.S. EPA Class VI Underground Injection Control (UIC) Program (Class VI Rule) and North Dakota Administrative Code Chapter 43-05-01.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Bending and shear improvements in 3D-printed core sandwich composites through modification of resin uptake in the skin/core interphase region

Fused deposition modeling (FDM) printed polymers are rarely used as a structural material due to anisotropic and low mechanical properties compared with conventional composites. In recent years, greater need has been expressed for recycling of materials, such as recyclable FDM, at the end of service life to reduce environmental pollution and manufacture cost. However, how the amount of resin uptake in the skin and skin/core interphase affects the bending and shear performance of the sandwich composites when replacing the low strength and ductile core (conventional core) with a high strength and brittle core (FDM printed PLA (polylactic acid) core) still remains unclear. A new manufacturing routine is needed to improve the incorporation of FDM printed polymers in composite structures. In this work, FDM printed PLA was used as core material and sandwiched between two unidirectional glass fiber reinforced polymer (GFRP) skins to form a sandwich composite by compression-molding (CM) process, which provides a good manufacturing strategy for skin/core interphase modification. The significance of the CM process is proved by investigating the effect of resin uptake on bending and in-plane/out-of-plane shear performances. Current first order shear deformation (FSDT) theory lacks a direct connection between the in-plane shear stress and out-of-shear stress in the core region of sandwich composites. With the help of DIC, a connection between the in-plane shear and the out-of-plane shear strain was built and in-plane shear properties can acquire through out-of-plane shear properties, hence reducing the redundancy of sample preparation or the need for simulation. A significant improvement was found compared with the optimized resin uptake (Optimized resin uptake range: 20.43%–22.86 wt%) 3D-printed PLA core sandwich composite and lowest performance sandwich composite (Improvement: in-plane shear strength (~34%)/modulus (~29%), out-of-plane shear strength (~25%)/modulus (~31%), specific peak bending load (~19%)). Finally, compared with balsa core sandwich composites, the 3D-printed cores are suitable for use in composite sandwich structures in many applications with a satisfactory strength-to-weight ratio.

3D-printed core sandwich composite↗
DOE OSTI2023

Image Characterization of Reactions Generated by an Aeroengine Micromixing Injector for Lean Direct Injection of Hydrogen and Hydrogen/Natural Gas Blend

An aeroengine micromixing injector, originally designed for lean direct injection of jet fuel, was adapted to work with hydrogen, natural gas, and any blend in between. The ultimate goal of the design was to achieve low NOx emissions when operated on pure hydrogen; to better characterize its performance, flame imaging diagnostics was implemented. Three types of cameras mounted parallel to the injector base and pointing towards the quartz combustor cylinder with the same angle were operated simultaneously: Nikon D90, Dynacolor FB-N9-U, and Phantom v7.1, to capture the visible spectrum, OH* chemiluminescence, and dynamic behavior of the flames, respectively. The current work also includes the overlapping of visible spectrum and OH* chemiluminescence, as an effort to qualitatively define the heat release over the flame area, i.e., UV over visible domain. Additionally, a z-type schlieren configuration was employed to reinforce the presence of some small-scale details occurring near the injector outlet ports. Using the airbox pressure drop, preheat temperature, fuel composition and flame temperature as the factors of study, a Box Behnken model was designed yielding into a 16-points matrix for the operability region, at atmospheric pressure. The main responses were extracted from the OH* images, because of its correlation with the flame heat release, these being the flame area and its average brightness, the heat release area and its center of gravity and leading edge. To assess the importance of the main factors and their interactions, an analysis of variance (ANOVA) was performed for each response, considering as significant each contributor with a p-value below 0.05. For these five responses it was found that all four main factors need to be included into the model (most of them because its own significance, others to ensure the model hierarchy), as well as some common interactions as the preheat and flame temperatures or pressure drop and fuel composition, and quadratic terms from fuel composition or flame temperature. Three of the five models presented R2 and coefficient of variance (C.V.) indicators around 0.97 and 5%, respectively; the leading edge ANOVA yielded into R2=0.78 and C.V.=18.88%; and the heat release area showed R2=0.91 and C.V.=24.13%. Lastly, the confinement ratio effect of the combustor over the flame structure was quantified, for three different conditions: the original 80 mm inner diameter and 200 mm length cylinder was compared against a shorter 150 mm tube (same inner diameter) and a narrower 47 mm inner diameter one (same 200 mm length). The length had a smaller effect than the inner diameter, even though both factors were not strong enough to make a significant deviation on most of the flame parameters: the difference was within the repeatability margin of error for the brightness, flame area, and center of gravity and leading edge of the heat release area.

Imaging diagnostics, flame characterization, OH* c↗
DOE OSTI2022

Generalization of thermodynamic Langmuir isotherm for mixed‐gas adsorption equilibria

Abstract This work presents a comprehensive thermodynamic model for both pure component isotherms and mixed‐gas adsorption equilibria. A generalization of thermodynamic Langmuir isotherm, the proposed model assumes competitive adsorption of multiple adsorbates on adsorbent surface for mixed‐gas adsorption equilibria, and it applies an area‐based adsorption nonrandom two‐liquid activity coefficient model in the activity coefficient calculations for the adsorbate phase. The resulting generalized Langmuir isotherm properly captures both surface loading dependence and adsorbate phase composition dependence for mixed‐gas adsorption equilibria. The model is validated with accurate representations of gas adsorption equilibrium data for varieties of unary, binary, and ternary gas systems. The model results are further compared with those calculated from extended Langmuir isotherm and Ideal Adsorbed Solution Theory.

Hamid, Usman↗
DOE OSTI2025

Supply Chain Energy and Greenhouse Gas Analysis Using the Materials Flows Through Industry (MFI) Tool: Examination of Alternative Technology Scenarios for the U.S. Chemical Sector

Chemical manufacturing is a large and diverse sector of the U.S. economy, with products, fuels, and a wide assortment of materials used daily by both the public and businesses. Currently, several of the largest volume chemicals produced in the United States rely on fossil fuels as a feedstock, energy source, or both. The list of chemicals includes steam cracking products such as ethylene, propylene, benzene, and xylenes as well as products such as ammonia and methanol. The focus for this work is on platform chemicals that are both produced in the largest volume and have a high potential for subsequent processing into more specialized products. In this study, we explore several new pathways that reduce the overall energy consumption and greenhouse gas (GHG) emissions for each product. These pathways include energy efficiency measures applied to existing production methods, the use of bio‐based fuels and/or feedstocks as new production methods, and electrification of high‐energy‐input stages within current production methods. Scenarios for energy demand and GHG reduction were conducted with the National Renewable Energy Laboratory's Materials Flows through Industry tool. Projections of the energy demand and GHG emissions in 2030 and 2050 are included, using grid composition projections from the NREL ReEDS model. The alternative scenarios selected showcase the effect of realistic changes the industry could make, focusing on technologies with a high level of technical readiness.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Optimizing the Composite Cathode Microstructure in All‐Solid‐State Batteries by Structure‐Resolved Simulations

Abstract All‐solid‐state batteries are considered as an enabler for applications requiring high energy and power density. However, they still fall short of their theoretical potential due to various limitations. One issue is poor charge transport kinetics resulting from both material inherit limitations and non‐optimized design. Therefore, a better understanding of the relevant properties of the cathode microstructure is necessary to improve cell performance. In this article, we identify optimization potentials of the composite cathode by structure‐resolved electrochemical 3D‐simulations. In our simulation study, we investigate the influence of cathode active material fraction, density, particle size, and active material properties on cell performance. Special focus is set on the impact of grain boundaries on the cathode design. Based on our simulation results, we can predict target values for cell manufacturing and reveal promising optimization strategies for an improved cathode design.

25 ENERGY STORAGE↗
DOE OSTI2020

Co-optimizing water-alternating-carbon dioxide injection projects using a machine learning assisted computational framework

In this report, a robust machine-learning-based computational framework that couples multi-layer neural network (MLNN) proxies and a multi-objective particle swarm optimizer (MOPSO) to design water-alternating-carbon dioxide injection (CO 2 -WAG) projects is presented. The proposed optimization protocol considers various objectives, including oil recovery and CO 2 storage volume. Expert MLNN systems are trained and employed as surrogate models of the high-fidelity compositional simulator in the optimization workflow. When multiple objective functions are considered, two approaches are employed to treat the objectives: the weighted sum method and the Pareto-front-based scheme. A field-scale implementation focusing on tertiary recovery in the Morrow B formation at Farnsworth Unit (FWU) is presented. The developed Pareto-optimal solutions indicate the maximal available oil production can be 1.64 × 10 7 barrels and maximal carbon storage can achieve 2.35 × 10 7 tons. Trade-offs factor is defined to divide the constructed Pareto front into 4 sections with the trade-off factors’ value ranges from 0.35 to 49.9. This work also compares the optimum solution found by the aggregative objective function and the solution repository covered by the Pareto front that considers the physical and operational constraints and reduces uncertainties involved by the multi-objective optimization process. Our comparison indicates multiple solutions exist to satisfy the objective criteria of the WAG design, and these results cannot be found using the traditional weighted sum method. The Pareto front solution can provide more options for project designers, but decisions regarding necessary trade-offs must be made using the solution repository to balance the project economics and CO 2 storage amount.

10 SYNTHETIC FUELS↗
DOE OSTI2025

In-service corrosion and grain boundary oxidation in neutron-irradiated 316 stainless steel baffle-former bolts

Reactor core internal components such as baffle-former bolts (BFBs) are subjected to significant mechanical stress, corrosive environment, and neutron irradiation from the reactor core during the plant operation. Over the long operation period, these conditions lead to potential degradation and of the bolts. In this work, characterization was performed on the oxidized surface of stainless steel BFBs harvested from a commercial pressurized water reactor (PWR) after 40 years of operation. The analysis shows that a complex multilayered surface oxide with six identified layers formed that is different from 2-layer structure commonly observed in model experiments. The oxide varies by composition – predominantly Fe, Cr, and Ni, grain size, and phase, and has features resembling both unirradiated and radiation/ corrosion experiments likely due to the low radiation flux compared to ion-irradiation or the test reactor radiation. In addition, grain boundary oxidative attack featured a pathway for Fe and other elements to move from the metal matrix to the outermost oxide. In conclusion, the results help assess PWR lifetime extension, put into context previous experimental studies, and provide input for designing experiments combining radiation and corrosion effects.

Baffle-former bolt↗
DOE OSTI2020

Assessment of CO 2 trapping mechanisms in partially depleted oil-bearing sands

The goal of this work is to evaluate various CO 2 sequestration mechanisms occurring in the Morrow B Sandstone in the Farnsworth Unit. A history-matched numerical simulation model was created using extensive geological, petrophysical, and operational data collected from the field. The numerical model is competent to investigate the impact of residual, structural-stratigraphic, solubility, and mineral trapping mechanisms on the fluid transportation dynamics and petrophysical property variations. The model forecasts the field response of 20 years of WAG injections. Afterward, all wells were shut-in, and the reservoir was allowed to evolve for 1000 years to investigate the fate of injected CO 2 . In this paper, we assess the impacts of various trapping mechanisms on oil recovery and CO 2 storage efficacy. By analyzing the results reported from the numerical simulation model, the in-situ fluid composition and mineralogy changes are also investigated. More importantly, we seek to confirm the petrophysical property variations due to the CO 2 injection with observations from laboratory measurements. The experiences gained from this study provide valuable insights regarding physiochemical storage induced by the CO 2 injection activities and serve as a benchmark case for future CO 2 enhanced oil recovery (EOR) projects involving reactive solute transport.

42 ENGINEERING↗
DOE OSTI2021

Thermal diffusivity and thermal conductivity of SiC composite tubes: the effects of microstructure and irradiation

Cladding thermal conductivity is an important physical property in assessing the performance of silicon carbide (SiC)-cladded fuels for nuclear reactors. However, there is a significant lack of reliable data, particularly for irradiated materials, because the geometry complicates the measurement. This study investigates the thermal diffusivity of coupons with a curvature, machined from SiC fiber–reinforced SiC matrix composite tubes, with and without neutron irradiation under light water reactor–relevant temperature and dose conditions. The tested materials included full composite and duplex SiC composite tubes. The measurements were conducted using a modern flash diffusivity apparatus. The analyzed area on the specimen during diffusivity testing was reduced for improved measurement accuracy due to sample curvature. Post-irradiation measurements showed that the effects of neutron irradiation on thermal conductivity (e.g. thermal defect resistivity) are different between SiC composite plates versus tubes. The difference was explained by higher matrix density of the tube than the plate. This study provides reliable thermal properties of prototypic SiC composite tubes useful for fuel performance modeling of SiC-based cladding.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Physics-informed machine learning of redox flow battery based on a two-dimensional unit cell model

In this paper, we present a physics-informed neural network (PINN) approach for predicting the performance of an all-vanadium redox flow battery, with its physics constraints enforced by a two-dimensional (2D) mathematical model. The 2D model, which includes 6 governing equations and 24 boundary conditions, provides a detailed representation of the electrochemical reactions, mass transport and hydrodynamics occurring inside the redox flow battery. To solve the 2D model with the PINN approach, a composite neural network is employed to approximate species concentration and potentials; the input and output are normalized according to prior knowledge of the battery system; the governing equations and boundary conditions are first scaled to an order of magnitude around 1, and then further balanced with a self-weighting method. Our numerical results show that the PINN is able to predict cell voltage correctly, but the prediction of potentials shows a constant-like shift. To fix the shift, the PINN is enhanced by further constrains derived from the current collector boundary. Finally, we show that the enhanced PINN can be even further improved if a small number of labeled data is available.

25 ENERGY STORAGE↗
DOE OSTI2023

Scalable aqueous-phase fabrication of reduced graphene oxide nanofiltration membranes by an integrated roll-to-roll (R2R) process

Economic fractionation of multicomponent biomass-derived streams is a key challenge in biobased fuels and chemicals production. Kraft black liquor (BL) is generated at ~1 billion tons/yr globally from biomass pulping and contains ~15wt% total solids including lignin, hemicellulose fragments, and inorganics. Membrane-based BL concentration is attractive but challenged by low solute rejections and poor stability in BL, which combines alkaline pH (~13), high dissolved solids content (15+ wt%) and high temperature (70–85°C). Here we discuss the scaled fabrication and characterization of reduced graphene oxide (rGO) nanofiltration membranes by slot die coating on a roll-to-roll (R2R) with integrated vacuum filtration. Here, we obtain high-quality 90×30 cm rGO membranes (referred to as “R2R-rGO membranes” in this work) with ~100 nm average thickness supported on porous poly(ethersulfone) (PES) sheets, and an effective area of 2700cm 2 . We present characterizations of their microstructure and uniformity as mapped over large areas by scanning electron microscopy (SEM), ellipsometry, adhesion test and X-ray photoelectron spectroscopy (XPS). The constructed ellipsometry model for the rGO membrane/PES substrate composite was successfully applied to 48 different locations of the R2R-rGO membranes. The R2R-rGO membranes showed >98% lignin rejection with a steady state flux of ~12 LMH at 50bar in a 15.7wt% total solids (TS) kraft BL stream at 72°C and 2–3 gal/min crossflow rates. Slot die coating on a R2R platform with integrated vacuum filtration can enable rapid fabrication of high-quality GO membranes at scale without organic solvents or volatile organic compounds.

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

Investigation of the thermal conductivity of molten LiF-NaF-KF with experiments, theory, and equilibrium molecular dynamics

Molten salts are being proposed for numerous advanced energy applications, including advanced nuclear reactors, concentrating solar power plants, thermal energy storage, and fusion reactors. Accurate knowledge of the thermophysical properties of molten salts directly impact the performance of these energy systems and are essential for design and safety analyses. Thermal conductivity data for fluoride molten salts and mixtures are especially lacking. In this work, experimental measurements of thermal conductivity using the steady-state variable gap technique were performed on eutectic LiF-NaF-KF from 834 to 1195 K. The experiment accounts for radiative, convective, and conductive heat losses. In addition, theoretical and molecular dynamics models are used, from 750 K up to 1300 K, to estimate the thermal conductivity for comparison with the experimental results. The results of experiments show a weak negative deviation of thermal conductivity with temperature, unlike previous experimental results in the literature. The measured thermal conductivity magnitudes agree with the theoretical and molecular dynamics predictions, aside from the data above 1100 K, where heat losses and radiative errors are the most significant, having a 16% maximum deviation from theory. These experimental results provide new thermal conductivity data for the LiF-NaF-KF system and further validation of the predictive models. The theoretical model was used to map the composition and temperature dependent thermal conductivity of LiF-NaF-KF and the mapping’s deviation from a linear additivity estimation of thermal conductivity. Additionally, this mapping showed the highest deviations from linearity for KF-LiF rich mixtures and increasing deviation with temperature. Notably, the deviation from linearity near the LiF-NaF-KF eutectic composition was around 25%.

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

William A. Bardeen: A life in physics and the legacy of the chiral anomaly

William Allan Bardeen (September 15, 1941 − November 18, 2025) was an American theoretical physicist who worked at the Fermi National Accelerator Laboratory. He is renowned for his foundational work on the chiral anomaly, the Adler-Bardeen theorem, the non-Abelian anomaly and gravitational anomalies. He was instrumental in the development of quantum chromodynamics and its applications, such as semileptonic decays and the Λ $\overline{MS}$ scheme frequently used in perturbative analysis of high energy processes involving strong interactions. Bardeen also played a major role in developing a theory of dynamical breaking of electroweak symmetry via top quark condensates, leading to one of the first composite Brout-Englert-Higgs boson models. His work on the chiral symmetry dynamics of heavy-light quark bound states correctly predicted abnormally long-lived resonances which are chiral symmetry partners of the ground state.

Hill, Christopher T. [Fermi National Accelerator L↗
DOE OSTI2025

Small-Molecule Models of Hydrogen-Evolving MX 2 (M = Mo, W; X = S, Se) Bulk Solids: Composition–Activity Relationships

Triangular metal chalcogenide clusters of the form [M 3 Q 7 L 3 ]An (M = Mo or W; Q = S or Se; L = i Bu 2 NCS 2 – , (CF 3 CH 2 ) 2 NCS 2 – , i Bu 2 NCSe 2 – , or i Bu 2 PS 2 – ; An = Cl – or I – ) have been investigated as molecular analogues of layered metal dichalcogenide (MX 2 ) H 2 -evolution catalysts. These clusters have been evaluated for their relative H 2 -evolving ability under a common photolysis protocol implementing [Ru(bpy) 3 ] 2+ as chromophore and Et 3 N as sacrificial electron donor. With M constant as Mo and with constant supporting ligand, clusters with an all-sulfide core enable greater H 2 -TON than clusters with an all-selenide core. A more active catalyst is produced by [Mo 3 S 7 (S 2 CN i Bu 2 ) 3 ] + I – than its W 3 analogue with the same core sulfide composition and supporting dithiocarbamate ligands. Dichalcogenocarbamate ligands provide more active catalysts than dialkyldithiophosphate ligated clusters, and within the dichalcogenocarbamate set, greater H 2 -turnovers correlate with more-electron-donating ligands (i.e., i Bu 2 NCS 2 – > (CF 3 CH 2 ) 2 NCS 2 – > i Bu 2 NCSe 2 – ). Cluster cations with Cl – as counteranion are very similar in activity H 2 -evolving levels to identical clusters with I – , ruling out any significant interfering effect by I – upon the electron transfer relay between Et 3 N and catalyst. In the aggregate, observations are consistent with a mechanism for H 2 evolution that involves reductive extrusion of H 2 from a metal hydride intermediate.

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

Intercalation of Lithium into Graphite: Insights from First-Principles Simulations

Understanding ion intercalation at electrode–electrolyte interfaces is key to the development of energy storage and water desalination. Here, we investigate Li + kinetics at a prototypical interface between graphite anodes and an organic electrolyte, and we elucidate key factors that determine ion transport, using first-principles methodology coupling ab initio molecular dynamics simulations with a solvation model. We show that surface chemical composition significantly influences the kinetics of ion intercalation from the liquid into graphite. We find that this is partly related to the ion desolvation process, which varies notably for different graphite surface chemical terminations. In addition, interfacial polarization is found to play an important role in determining energy barriers for ion transfer. We also discuss the impact of electrode potentials, which is often neglected in conventional first-principles calculations despite being a key factor in device configurations. Our study provides insights into the coupling of electronic and ionic effects of interfacial chemistry on ion transport at complex electrode–electrolyte interfaces.

25 ENERGY STORAGE↗
DOE OSTI2020

Predicting the Activity and Selectivity of Bimetallic Metal Catalysts for Ethanol Reforming using Machine Learning

Machine learning is ideally suited for the pattern detection in large uniform datasets, but consistent experimental datasets on catalyst studies are often small. Here we demonstrate how a combination of machine learning and first-principles calculations can be used to extract knowledge from a relatively small set of experimental data. The approach is based on combining a complex machine-learning model trained on an extensive computational library of transition-state energies with simple linear regression models of experimental catalytic activities and selectivities from the literature. Using the combined model, we identify the key C–C bond scission reactions involved in ethanol reforming and perform a computational screening for ethanol reforming on monolayer bimetallic catalysts with architectures TM-Pt-Pt(111) and Pt-TM-Pt(111) (TM = 3d transition metals). The model also predicts four promising catalyst compositions for future experimental studies. In conclusion, the approach is not limited to ethanol reforming but is of general use for the interpretation of experimental observations as well as for the computational discovery of novel catalytic materials.

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

Ab Initio Simulation of Medium-Range Ordering in Ionic Glass Electrolytes LiSiPON and LiNaSiPON

Ionic glasses can exhibit a unique combination of optical transparency, electronic resistivity, ionic conductivity, and mechanical ductility. The origin of the relatively high ionic conductivity and ductility is poorly understood. Recently, these ionic glasses were found to have medium-range ordering (MRO) similar to that of metallic glasses. This MRO significantly impacts the properties of metallic glasses and is also expected to significantly impact the properties of ionic glasses. Here, this work used ab initio molecular dynamics (AIMD) simulations to study the effect of ionic glass composition on the MRO. The AIMD models showed that the degree of MRO increased as the temperature of the LiSiPON ionic glass decreased. The modeling also showed that the MRO is suppressed when 50% of Li is substituted with Na. This work lays the foundation for using AIMD to identify clear structure–property relationships in ionic glasses, allowing their properties to be optimized.

Osetsky, Yuri N. [Oak Ridge National Laboratory (O↗