Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ANL”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Cradle-to-Gate greenhouse gas emissions of the production of ethylene from U.S. Corn ethanol and comparison to fossil-derived ethylene production

Conventional ethylene production heavily depends on fossil-derived feedstocks via steam cracking, a very energy- and emission-intensive process. Researchers have been exploring alternatives to reduce CO 2 emissions including producing ethylene from biobased feedstocks. This paper evaluates the cradle-to-gate greenhouse gas (GHG) emissions of bioethylene produced from U.S. corn ethanol. The analysis includes different pathways for the dehydration of corn ethanol to ethylene and co-processing routes via fluid catalytic cracking (FCC) processes. For the FCC co-processing route carbon-14 analysis is used to determine bioethanol yields. A 127% reduction in life cycle GHG emissions of bioethylene is estimated compared to fossil-derived ethylene for the base case. Additional case studies are also discussed to understand the reduction of GHG emissions due to sustainable corn farming and renewable power use, biogenic carbon capture, and fuel switch with biofuels at the ethanol plant, and its impact on bioethylene GHG emissions.

carbon footprint↗
DOE OSTI2025

LABQ3: Bayesian method for quantification of mineral compositions and nano-scale elemental mapping of 3D synchrotron XCT data

Quantitative analysis of mineral compositions is essential in understanding geochemical, mineralogical and environmental processes. Fine-resolution 3D imaging is widely done using synchrotron X-ray computed tomography (XCT), but existing analyses are limited to visualization and segmentation. This paper presents a new method, Linear Attenuation Bayesian Quantitative 3D-mapper (LABQ3), based on the linearity of X-ray attenuation with respect to elemental concentrations. To address the random variability in attenuation measurements, LABQ3 employs Bayesian decision theory to minimize classification error, using reference attenuation distributions from scans of pure mineral standards. To demonstrate LABQ3 and test its performance, we studied precipitated carbonate samples. XCT scans were done at multiple energies using the transmission X-ray microscope (TXM) at beamline 32-ID-C of the Advanced Photon Source at Argonne National Laboratory. The reconstructed 3D images have a voxel size of 20 nm. Analyses revealed rich nano-scale compositional heterogeneity within individual particles. A mixture of calcium and cadmium produced an overall stoichiometric composition of (Ca 0.78 ,Cd 0.22 )CO 3 , with some voxels containing nearly pure CdCO 3 . The addition of zinc led to an overall stoichiometric composition of 33% Ca, 28% Cd, 39% Zn, with a nearly pure CaCO 3 core and compositional zonation through the rim. These compositional gradients are related to temporal sequences of carbonate mineral formation where Cd precipitated at the beginning in (Ca,Cd)CO 3 , while Cd and Zn precipitated at the end in (Ca, Cd,Zn)CO 3 . Results differ from bulk analyses using Inductively Coupled Plasma-Mass Spectrometry (ICP-MS), showing that LABQ3 provides particle-specific insights. LABQ3 distinguishes itself by quantifying chemical compositions along a continuum, making it different from XCT analyses based on segmentation. LABQ3 allows simultaneous acquisition of morphology and chemical composition in 3D, facilitating the interpretation of chemical gradients of trace elements, quantification of solid solution compositions, inferences about temporal sequences of mineral precipitation, and addressing other concerns about solid-phase chemistry.

58 GEOSCIENCES↗
DOE OSTI2025

Evaluation of IrO 2 catalysts doped with Ti and Nb at industrially relevant electrolyzer conditions: A comprehensive study

A series of commercial Oxygen Evolution Reaction (OER) IrO 2 -based materials doped with acid-stable titanium and niobium species were comprehensively characterized by Brunauer-Emmett-Teller (BET), X-ray diffraction analysis (XRD), X-ray photoelectron spectroscopy (XPS), transmission electron microscopy (TEM) with energy dispersive spectroscopy (EDS), and X-ray Scattering. Electrocatalysts were integrated into Membrane Electrode Assembly (MEA) using a fabrication method developed under the US DOE H2NEW consortium. An electrolysis performance in a commercial setup as well as a laboratory screening system was performed at conditions relevant to industrial application. According to the comprehensive characterizations, the studied materials are closer to doped iridium oxides rather than core–shell structures. In an electrolysis cell, the IrO 2 /TiO x catalyst slightly outperforms the IrO 2 /NbO x based on the activity. It was demonstrated that the operation of electrolysis cells at elevated temperatures and the implementation of thinner Nafion-type membranes allows for substantially increased performance, which is consistent with the literature report. Finally, this work provides valuable baselines including characterization and performance for guiding future research in this direction.

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

Current activated reactive ultrafast joining (CARUJ) of silicon carbide

In this work we propose and demonstrate a novel approach for rapid fabrication of ceramic-ceramic joint assemblies, using sintered Silicon Carbide (SiC) as an initial example. Current Activated Reactive Ultrafast Joining (CARUJ) utilizes resistive heating of carbon-based materials to apply localized heat at or around the joint zone at heating rates of 101– 103 °C/min. CARUJ is used to fabricate SiC-SiC joints with a Si-SiC interface using minimal pressure (1–2 MPa) at time scales considerably shorter (several minutes opposed to hours) than conventional approaches for similar systems. The reaction of an interfacial precursor based on elemental silicon and carbon leads to in-situ SiC formation to produce a continuous and dense bond in a single step. Measured average joint strengths of roughly 15 MPa are achieved when tested in single lap offset (SLO) compressive shear. Optional additions of refractory metals such as molybdenum can be utilized to introduce secondary inclusions such as MoSi 2 within the joint interface. We further demonstrate SiC-SiC joining using Active Brazing Alloys (ABA) and proof of concept joining of tubular geometries. The localized and rapid heat application realizes a versatile material joining technique that could be extended for joining components at the plant site.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Upcycling mixed cathode materials to high-energy-density LiFe 0.75 Mn 0.25 PO 4

To address the demand for next-generation cathode materials with high energy density, upcycling LiFePO 4 into LiFe 0.75 Mn 0.25 PO 4 has attracted considerable attention. Nevertheless, existing strategies have yet to achieve both morphology and full elemental recovery under mild ambient conditions. Here, we report an upcycling route that can address this issue by combining leaching and a high-temperature treatment process. The upcycled LiFe 0.75 Mn 0.25 PO 4 exhibits enlarged lattice spacing and a high discharge plateau, and it delivers an energy density of 563.7 Wh/kg, 40.3 Wh/kg higher than recycled LiFePO 4 , which indicates the high value of the proposed upcycling strategy. At 1 C, LiFe 0.75 Mn 0.25 PO 4 also exhibits excellent cycling stability of 91% over 700 cycles. Techno-economic analysis also indicates impressive economic and environmental benefits, including 10.4% less raw materials usage and 12.2% less energy consumption and wastewater generation. This work demonstrates a scalable and economic upcycling strategy and provides a promising pathway for sustainable battery upcycling compatible with industrial conditions.

LMFP↗
DOE OSTI2025

A discontinuous Galerkin spectral element method for compressible reacting flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large-eddy simulations because of their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reacting Navier-Stokes equations. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of the DG approach. The framework, implemented in the spectral element code Nek5000, is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. An entropy-residual based artificial viscosity is added to smooth shocked regions of flow, and a positivity-preserving limiter is implemented to suppress non-physical oscillations. These enhancements support the numerical stability of the hydrodynamic sub-step, which is decoupled from the chemistry integration through a second-order operator splitting method. Here, a series of smooth and discontinuous validation cases are presented in increasing physical and computational complexity for both inviscid and viscous flows. In particular, simulations of canonical one-dimensional and two-dimensional detonations are performed, and the high-order numerical results are validated against available literature data. Additional validation studies are carried out for classical three-dimensional numerical simulations of incompressible and compressible turbulent flows.

Compressible reacting flows↗
DOE OSTI2024

Salt matters: How ionic strength and electrolytes impact redox polymer reactivity and dynamics for energy storage

As the global demand for sustainable energy grows, redox-active polymers (RAPs) have emerged as promising materials for batteries due to their advantages in stability, ease of preparation, and low-cost processability. Despite factors traditionally known to impact polymer dynamics (e.g., temperature, viscosity, and structure), we posit that investigating the effect of ionic strength and/or supporting electrolyte types on the electrochemical performance of RAP systems is crucial, both in aqueous and nonaqueous systems. Here, we first highlight recent findings on RAP-electrolyte interactions, elucidating how their polyelectrolyte nature determines their redox activity. Then, we focus on strategies to enhance RAP performance for energy storage through ionic strength optimization and tailored electrolyte composition. These insights into the modulation of RAP reactivity provide a foundation for improving battery performance in both flow and stationary configurations, thus facilitating progress toward next-generation energy storage solutions.

25 ENERGY STORAGE↗
DOE OSTI2025

Turbulent burning velocity of lean premixed hydrogen/air flames at engine conditions: Effects of turbulence intensity and length scale

For turbulent lean premixed hydrogen flames with strong thermodiffusively instabilities, most previous studies have focused on the influence of turbulence intensity, whereas the role of turbulence length scale is less well understood. Here, this study addresses this gap by conducting direct numerical simulations (DNS) of statistically planar turbulent premixed flames for a lean (ϕ=0.35) hydrogen/air mixture under independently varied turbulence intensity (u') and length scale (l T ) at engine-relevant thermodynamics conditions. Results show that as u' increases, the flame front becomes increasingly wrinkled, forming smaller cellular structures. In contrast, l T variations do not significantly alter the size of these structures. For the turbulent burning velocity (S T ), the normalized S T (i.e., S T /S L , where S L is the laminar flame speed) increases linearly with u', driven by both enhanced flame surface wrinkling (i.e., increased A T /A L ) and enhanced local burning rate (i.e., increased I 0 ). However, increasing l T reduces I 0 , despite a continued increase in A T /A L , resulting in only a marginal increase in S T /S L . To reveal the underlying mechanisms, especially the decreasing trend of I 0 with l T , local flame dynamics analyses are performed. It is found that as l T increases, the interaction between thermodiffusive effects and turbulence weakens due to the reduced tangential strain rate, while the flame curvature remains largely unchanged. This suppresses local reactivity enhancement and thus decreases I 0 , In contrast, an increase in u' enhances the interaction by amplifying both curvature fluctuation and tangential strain rate, leading to increased local reactivity (increased I 0 ). Finally, based on the DNS data, several new scaling models are proposed for the three global properties, S T /S L , A T /A L , and I 0 , and showed improvements compared to existing models. These findings provide new insights into the flame-turbulence interactions in thermodiffusively unstable hydrogen flames. The DNS dataset is also useful for the development of turbulent combustion models applicable to practical engine simulations.

Engine-relevant condition↗
DOE OSTI2025

HFBTHO-AD: Differentiation of a nuclear energy density functional code

The HFBTHO code implements a nuclear energy density functional solver to model the structure of atomic nuclei. HFBTHO has previously been used to calibrate energy functionals and perform sensitivity analysis by using derivative-free methods. To enable derivative-based optimization and uncertainty quantification approaches, we must compute the derivatives of HFBTHO outputs with respect to the parameters of the energy functional, which are a subset of all input parameters of the code. Here, we use the algorithmic/automatic differentiation (AD) tool Tapenade to differentiate HFBTHO. We compare the derivatives obtained using AD against finite-difference approximation and examine the performance of the derivative computation.

Algorithmic differentiation↗
DOE OSTI2025

The value of hydropower flexibility for electricity system decarbonization

Hydropower is an abundant, dispatchable, clean energy resource that will play an important role in supporting the clean energy transition. In particular, dispatchable hydropower can provide the operational flexibility that will be required in future systems with high variable renewable energy penetrations. However, the theoretical operational flexibility of hydropower can be restricted in practice by various non-power constraints. In this paper, we quantify how increasing the operational flexibility of dispatchable hydropower resources with reservoirs impacts least-cost generation portfolios and supports power system decarbonization. Specifically, we conduct a capacity expansion analysis of a two-zone system: a hydro-dominated region and a neighboring region with aggressive decarbonization targets that are represented by the United States Pacific Northwest and California respectively. We then introduce a quantifiable index for characterizing the operational flexibility of reservoir hydropower and assess how changes in this metric impact the system-optimal generation portfolio. We find that increasing hydropower flexibility leads to more investment in wind generation, less investment in natural gas generation, lower system costs, and lower system emissions. We further demonstrate a substitution effect between the grid services provided by flexible hydropower operation, increased transmission capacity on a congested line, and energy storage resources. Finally, we show that increasing the operational flexibility of hydropower increases the effective load carrying capability of both hydropower and wind resources. This research supports a more nuanced understanding of how hydropower can support electricity system decarbonization and may motivate reassessing the cost-benefit tradeoffs of non-power constraints that restrict operational flexibility.

Capacity expansion modeling↗
DOE OSTI2025

Electrodeposition of near-equiatomic CoCuFeNi multi-principal element alloys from an acidic glycine-citrate-triton X100 aqueous electrolyte

Understanding the composition and morphology control of electrodeposited CoCuFeNi is the first step to finding a general strategy for developing electrodeposition processes of unconventional alloys with large redox potential differences and complicated deposition mechanisms. Here, in this work, we have successfully synthesized the near-equiatomic (<5 at% error) CoCuFeNi films with ~200nm thickness by electrodeposition from glycine-citrate-Triton X-100 acidic electrolytes. This system generally follows the Principle II in Brenner’s paradigm on alloy composition control in electrodeposition. X-ray diffraction (XRD) profiles show that the films only consist of one crystalline phase, different from the deposit from ammonia-citrate-boric acid electrolytes and the equilibrium phases predicted by CALPHAD. The near-equiatomic deposits at Ru substrates were successfully annealed at 400 ° C without significant intermixing between the substrate and the deposit, in which no phase separation of the crystalline phase was observed in its XRD profile. Xray photoelectron spectroscopy (XPS) reveals that some depositing elements (mainly Fe) exist as metal oxides. Cu is dissolved in the crystalline phases, stabilized by small crystalline domain size and the amorphous metal oxides inside the CoCuFeNi deposits.

CoCuFeNi↗
DOE OSTI2025

An update to the Sandia method for creating Typical Meteorological Years from a limited pool of calendar years

Typical Meteorological Years (TMYs) are essential for the efficient evaluation of energy system performance. Ideally, 30 years of weather data are required to generate TMYs, but significantly fewer years are typically available due to practical limitations. To address this issue, an update to the Sandia method was developed, referred to as the Argonne method, to create TMYs from a limited number of years. Furthermore, this method enhances candidate diversity by systematically shifting original candidate months forward or backward by specific days, creating an expanded pool of candidates. The effectiveness of the Argonne method was validated through statistical testing, comparison of monthly average weather parameters, and numerical simulations. The results demonstrate a high probability of identifying at least one shifted month whose cumulative distribution functions of weather parameters closely align with long-term distributions. In 67 % of all comparisons, the monthly average weather parameters in TMYs generated using the Argonne method exhibit better agreement with long-term averages than TMY3. Moreover, in 74 % of the 318 building simulation cases, the Argonne method outperforms TMY3 in estimating long-term average building heating and cooling demands. Therefore, the Argonne method effectively diversifies the candidate pool and produces typical years that provide more accurate estimations of long-term averages compared to TMY3 when only a limited pool of calendar years (10 years or fewer) is available.

Building energy modeling↗
DOE OSTI2025

Thermal decomposition pathways of bulk electrolytes on vanadium oxide nanocrystals

The thermal stability of electrolytes at an elevated temperature induced by battery charge-discharge cycling is critical for the long cycling performance of a rechargeable battery. For many multivalent systems, such as rechargeable magnesium batteries, which offer great potential for high energy density and utilize earth-abundant resources, electrolyte instability and electrode surface passivation, arising from electrolyte decomposition, remain as major roadblocks. Understanding the electrolyte decomposition pathways at the electrode-electrolyte interface is essential to provide guidance in overcoming this challenge. In this work, in situ 13 C magic angle spinning nuclear magnetic resonance (MAS NMR) and first-principles calculations were used to investigate the thermal decomposition of the electrolyte in a system consisting of MgV 2 O 4 , a novel cathode for magnesium batteries, mixed with a bulk electrolyte consisting of magnesium bis(trifluoromethanesulfonyl)imide (Mg(TFSI) 2 ) in diglyme (G2). We show that significant electrolyte decomposition is observed in bulk 1.0 M Mg(TFSI) 2 in G2 mixed with nanometer sized MgV 2 O 4 powder at elevated temperatures. This observation is to mimic the possible thermal decomposition that might happen during battery cycling. We demonstrate that the MgV 2 O 4 surface is covered by a layer of decomposed G2 products. We conclude that the dominant reaction pathway for electrolyte decomposition is the thermal decomposition of the pure electrolytes at elevated temperatures, followed by adsorption of G2 decomposition products to the MgV 2 O 4 surface. The activation energy for the major decomposition pathway is obtained. In conclusion, this work highlights the importance of studying thermal decomposition of electrolytes for overall system stability and explores electrolyte stability at significantly elevated temperatures.

25 ENERGY STORAGE↗
DOE OSTI2025

Upcycling Mixed Spent Ni-Lean Cathodes into Ni-Rich Polycrystalline Cathodes

Sustainable battery recycling is vital for conserving resources and reducing environmental impacts. Current open- and closed-loop recycling strategies often focus on recovering individual components, making the reuse of mixed cathode materials a complex challenge. Meanwhile, the research on upcycling has been limited to using pristine cathode feedstocks and virgin materials for synthesis. Here, to address this issue, we present an upcycling approach for spent Ni-lean mixed cathode materials that integrate an upcycling hydrometallurgical recycling process with traditional hydrometallurgical methods. This strategy achieves a utilization of 92.31 mol % of recycled materials, enabling the regeneration of Ni-rich cathode materials while significantly reducing the reliance on virgin resources. The regenerated 83Ni cathode materials demonstrate physical properties comparable to those produced from virgin materials. Electrochemical evaluations using single-layer pouch cells show that both recycled and virgin cathodes exhibit initial specific capacity close to 201.1 mAh/g and maintain approximately 88 % capacity retention after 500 cycles. Additionally, 2Ah cells confirmed these findings, delivering 85 % capacity retention after about 900 cycles. Techno-economic analysis demonstrates notable environmental benefits, including reductions in greenhouse gas emissions and energy consumption, achieving 232.75 MJ/kg of product, which is 8.6 % lower than traditional methods and comparable to direct upcycling. Furthermore, the upcycling hydrometallurgical recycling process generates the highest profit, proving its economic viability. This scalable and versatile process is adaptable to varying transition metal compositions, facilitating a closed-loop recycling system that bridges mixed spent cathodes with next-generation cathode materials, and offers a sustainable solution for managing waste battery materials.

Hydrometallurgical recycling↗
DOE OSTI2025

Thermodynamic modeling of aqueous acetic acid, butyric acid, and lactic acid solutions

Based on the activity coefficient – fugacity coefficient approach, a rigorous thermodynamic modeling study is presented for accurate correlation of vapor-liquid equilibrium data of aqueous solutions of acetic acid (293 to 391 K), butyric acid (325 to 436 K), lactic acid (378 to 409 K), and acetic acid + butyric acid binary mixture (358 to 421 K). In addition, the pH data of the three aqueous, single carboxylic acid solutions were measured at 298 to 328 K and successfully correlated. Given that these aqueous carboxylic acid solutions exhibit various degrees of association behavior in both vapor and liquid phases, the thermodynamic models considered for this study include the Redlich-Kwong equation of state (RK-EoS) and the Hayden-O’Connell equation of state (HOC-EoS) for the vapor phase fugacity coefficients and the electrolyte non-random two-liquid model (eNRTL) and the association electrolyte non-random two-liquid model (AeNRTL) for the liquid phase activity coefficients. The combination of the HOC-EoS for the vapor phase and the AeNRTL model for the liquid phase is found to provide the best correlation results, consistent with the fact that the HOC-EoS and the AeNRTL model explicitly account for association behaviors in the vapor phase and the liquid phase, respectively.

09 BIOMASS FUELS↗
DOE OSTI2025

CFD modeling of near-wall combustion and unburned methane prediction in natural gas spark ignition engines

Natural gas-powered engines play a critical role in gas drilling, compression, and transmission sectors, but methane (CH 4 ) from engine combustion slip can be significant over their lifespan, contributing to atmospheric pollution and signaling reduced engine efficiency. Here, to address this challenge, computational fluid dynamics (CFD) simulations offer valuable insights into the in-cylinder combustion process, enabling the optimization of combustion strategies and engine designs to minimize unburned CH 4 slip. This study aims to evaluate and improve combustion models for simulating the combustion process and predicting unburned CH 4 concentrations in natural gas spark-ignition (SI) engines, including engines that are part of combined reformer-engine systems. Specifically, the performance of two flamelet-based combustion models—the Extended Coherent Flame Model (ECFM) and the G-equation model—was assessed using experimental engine data collected under varying excess-air ratio (λ) conditions and fuel compositions, including natural gas and syngas blends. In addition, to enhance the predictive capabilities of the G-equation model, a flame-wall interaction (FWI) sub-model was integrated into its framework. The effects of its model parameters, such as quenching and influence distance, on combustion behavior and unburned methane predictions were analyzed in detail. The ECFM tended to predict delayed combustion phasing under diluted mixture conditions, resulting in overprediction of unburned CH 4 concentrations. In contrast, the G-equation model provided reasonable predictions of combustion pressure, while representing higher the CH 4 reduction rate across the operating condition compared to experimental data. Incorporating the FWI sub-model—with the quenching distance calculated based on a pressure-dependent relation (P -0.48 ) and a fixed influence distance of 1.5 mm—further improved the G-equation model’s accuracy in predicting CH 4 reduction rates without compromising its ability to simulate the combustion process.

Combustion model↗
DOE OSTI2025

Computing material volume fractions on a superimposed mesh as applied to Monte Carlo particle transport simulations

Here, we present a newly implemented ray tracing algorithm in OpenMC for efficiently computing material volume fractions on superimposed meshes in complex geometries. By firing rays along each coordinate direction through the geometry, the approach accumulates track-length data in each mesh element, thereby determining the fractional composition of each material. Scaling studies on three different models—a random tetrahedra configuration, the Frascati Neutron Generator ITER dose rate benchmark, and a stellarator design—show excellent parallel performance, with nearly linear speedup on modern multi-threaded and distributed-memory systems. An analysis of the residual error relative to high-resolution reference solutions demonstrated that under optimal conditions it decreases as 1/R, where R is the number of rays fired, making it straightforward to achieve user-prescribed accuracy. This new functionality enables practical, mesh-based approaches for detailed nuclear analyses in production Monte Carlo workflows without resorting to expensive, fully conformal or unstructured meshing.

Monte Carlo↗
DOE OSTI2025

Quantifying spatial and vertical variations in soil C:N relationships in permafrost-affected landscapes

Permafrost regions are experiencing rapid changes that affect carbon (C) and nitrogen (N) cycles, with implications for vegetation dynamics and gas exchanges with the atmosphere. Soil C:N ratio is a key indicator of organic matter quality, yet spatial estimates of N stocks and C:N ratios lag behind those for C. We used quantile regression forests to compare direct and indirect digital soil mapping approaches for predicting soil C:N ratios at 0–30, 30–60, and 60–100 cm depths across a latitudinal transect in Alaska. The indirect approach – deriving C:N from separately predicted C and N stocks – outperformed direct mapping for the surface layer (0–30 cm), while direct mapping was marginally better at greater depths. However, prediction accuracy decreased with depth for both methods. Temperature and topography were the most important predictors. Both approaches overestimated low and underestimated high C:N ratios, with direct mapping showing greater bias. Our results underscore the challenges of modeling C:N ratios in heterogeneous, data-sparse permafrost soils, but also suggest that indirect mapping holds promise if supported by more extensive datasets.

54 ENVIRONMENTAL SCIENCES↗