Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “boulder”

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,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Outstanding Questions and Future Research on Magnetic Reconnection

This short article highlights unsolved problems of magnetic reconnection in collisionless plasma. Advanced in-situ plasma measurements and simulations have enabled scientists to gain a novel understanding of magnetic reconnection. Nevertheless, outstanding questions remain concerning the complex dynamics and structures in the diffusion region, cross-scale and regional couplings, the onset of magnetic reconnection, and the details of particle energization. We discuss future directions for magnetic reconnection research, including new observations, new simulations, and interdisciplinary approaches.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

End-to-end microgrid protection using distributed data-driven methods

This paper introduces an end-to-end microgrid protection framework that offers real-time system monitoring, fault-related decision making, and circuit breaker control. This is achieved through the design of distributed data-driven techniques based on the support vector machine method, where each relay is responsible for distributed data collection, fault detection, fault localization, and fault isolation. Local communication is established among neighboring relays, fostering cooperative fault localization and isolation. This decentralized design not only reduces the computational and communication requirements but also enables the adaptability of each relay under varying operational dynamics. The proposed end-to-end protection framework was validated using MATLAB/Simulink simulations on a 100% renewable microgrid, achieving an accuracy of 93.1% with response time of 0.0523 s, in protecting against a range of fault scenarios that are characterized by various types, locations, impedances, load conditions, photovoltaic power levels, and microgrid operating modes.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Theoretically Informed Kinetics (ThInK): Establishing a modern C 0 -C 3 mechanism for combustion modeling

In contrast to the adage “Models are to be used, not believed”, combustion kinetics models have been intended to be predictive in nature. Theoretical chemical kinetics is now understood to provide a firm foundation for the reaction parameters, thereby facilitating predictive simulations of chemical reactivity, even in regimes that are poorly characterized by chemical kinetic and/or combustion experiments. In this study, we describe a theory-informed kinetics model (ThInK) for small molecule combustion chemistry (H 2 and C 1 – C 3 species) that is based on the prodigious use of theoretical predictions for reaction rate coefficients, thermochemistry, and transport parameters. The distinct features of this kinetics model, which was developed over the course of several decades, are illustrated through simulations of flame propagation and auto-ignition.

Energy transfer↗
DOE OSTI2026

Achieving high rate performance in hybrid pristine-recycled cathodes using model-informed electrode designs

Direct recycling lithium-ion battery cathodes, a process that retains the engineered oxide structures from end-of-life materials, presents a cost-effective and energy-efficient alternative to other battery recycling methods. However, while direct-recycled cathodes have demonstrated performance comparable to that of pristine materials at low cycling rates, their high-rate performance remains uncertain. Morphology changes in cathode particles, a main mode of degradation, directly impact rate performance by limiting surface kinetics and solid-phase diffusion. If direct recycling processes do not sufficiently restore pristine-like morphologies, the recycled materials may retain structural defects that hinder high-rate performance. The present work uses a physics-based pseudo-2D model to simulate hybrid electrodes with pristine and artificially “aged/recycled” NMC materials to investigate potential impacts of incorporating performance-limited aged cathode materials into cells. The study highlights how differences in transport and kinetic properties can influence rate capabilities in mixed electrodes — particularly in high-loading cells in high-demand applications. However, model results also reveal a possible mitigation strategy via dual-layer electrode architectures with lower-performing materials positioned near the current collector. Simulations of 4.0 mAh cm −2 cells cycled at 4C using a dual-layer architecture provided approximately 5%–30% more capacity in constant-current protocols compared to homogeneously blended electrode architectures with the same loadings and mixed-material compositions. These findings highlight the importance of strategic electrode design in minimizing potential performance losses and facilitating the integration of recycled materials into high-performance batteries, advancing sustainable and cost-effective battery manufacturing.

25 ENERGY STORAGE↗
DOE OSTI2024

Optimal operation of multi-plant steam district heating systems for enhanced efficiency and sustainability

Despite their crucial role in supplying heat and power to universities, industries, and healthcare facilities, many steam-based district heating systems rely on outdated control methods. Among these, multi-central plant districts are particularly challenging due to the complexities of coordinating multiple plants, optimizing load distributions, and managing system downtime. In response, new operational strategies are developed to enhance the efficiency and sustainability of steam districts while utilizing existing resources. These strategies include reducing plant operational pressure without compromising the reliable supply to buildings and optimizing load allocation across multiple plants. The load allocation considers boiler part-load efficiency, runtime, network losses, and building pressure set points, and is compared with traditional multi-boiler controls. To support this exploration, new dynamic Modelica models are developed. In addition, methods to reduce modeling complexities are incorporated, enhancing their suitability for practical applications. A holistic district-wide analysis using a real university case study demonstrates a 4.7% fuel savings by lowering boiler operational pressure from 900 kPa to 600 kPa, along with a 13.3% reduction in condensation losses across the distribution network. Furthermore, the load allocation approach results in a 13.1% reduction in fuel consumption during peak winter periods and 15.3% during shoulder periods, with corresponding decreases in carbon emissions and fuel costs. This approach can also save maintenance costs by reducing the boiler runtime by 49.6%. In conclusion, this research underscores the benefits of retrofitting aging steam district heating systems, offering immediate operational improvements by enhancing efficiency, meeting regulatory compliance, and extending infrastructure lifespans while delaying costly overhauls.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Development of a River Dynamical Core for E3SM to simulate compound flooding on Exascale-class heterogeneous supercomputers

Flooding events pose significant risk to human life, property, and infrastructure. Physically-consistent quantification of altered flood risks in global models requires hyper-resolution (~1 km) or fine flood simulations using two-dimensional (2D) physics schemes, both of which are unavailable in the current generation Earth System Models. Here, in this work, we have developed the River Dynamical Core (RDycore), which is an open-source, 2D shallow water equation (SWE) library for the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM). RDycore uses PETSc and libCEED libraries that allows it to run efficiently on CPUs and GPUs, as well as select a time-integration algorithm at runtime without requiring any code modifications. RDycore achieves spatial error convergence rates for problems with analytical and manufactured solutions similar to those reported previously in the literature, or consistent with the implemented first-order spatial discretization scheme. RDycore's accuracy in predicting flooding for a well-studied dam break problem is comparable to existing SWE models. For a problem with 471 million grid cells, RDycore achieves a speedup of 6.6x and 7.6x on GPUs compared to CPUs when using 320 compute nodes on DOE's Perlmutter and Frontier supercomputers, respectively. The one-way coupling of the RDycore library within E3SM is demonstrated by performing multiple 5-day flooding simulations during Hurricane Harvey driven by five precipitation datasets. The E3SM--RDycore simulations at 30 m spatial resolution accurately simulate maximum water height during the hurricane when benchmarked against a previously published study and achieve a speedup of 15x (Perlmutter) and 21x (Frontier) on GPUs relative to CPUs. The work presented here is the foundational step in providing hardware and algorithmic portability framework for simulating kilometer-scale river dynamics within E3SM.

Flood Simulation↗
DOE OSTI2025

Earthquake-triggered submarine canyon flushing transfers young terrestrial and marine organic carbon into the deep sea

Submarine canyons transfer substantial amounts of sediment and organic carbon (OC) into the deep ocean, nourishing deep-sea ecosystems and contributing to the global carbon cycle through OC burial and sequestration. Tracking lateral OC transport through submarine canyon systems is challenged by the deep-ocean setting, difficulties with constraining episodic depositional events, and the need to assess the composition and age of marine and terrestrial organic matter. We apply innovative parallel ramped pyrolysis oxidation accelerator mass spectrometry and pyrolysis-gas chromatography-mass spectrometry with isotope analyses to track OC age and sources in the 2016 Kaikoura earthquake-triggered, canyon-flushing event that deposited along >1300 km of a submarine canyon-channel system, offshore Aotearoa New Zealand. Specifically, these techniques allow us to determine the ages, sources, and partitioning of OC within the Kaikoura turbidite deposit and test hypotheses of how submarine canyon systems contribute to lateral OC flux and burial. Our results show that, despite considerable canyon floor erosion, substantial amounts of young OC were flushed into the deep sea, with relatively little (~2%) pre-Holocene OC contributions. Even without a direct connection between rivers and submarine canyons, most (~55%) of the OC in the Kaikoura event bed is from terrestrial sources. However, the deposit also contains substantial amounts (~22%) of marine-derived OC and ~23% of the material is of unassignable origin. Particle sorting imparts variability on the age and composition of OC within turbidite deposits and along the turbidity current flow path. Terrestrial-derived OC is preferentially older than marine-derived OC and concentrated in coarser particle sizes found more commonly at the deposit base and in proximal settings. Young, marine-derived OC is concentrated at the surface of the deposits and tends to be enriched in finer particle sizes. Such OC partitioning in turbidites supports the relevance of depositional models for predicting and quantifying distribution of OC in deep-sea deposits. Earthquake-triggered, canyon flushing events and resulting turbidites enhance OC burial efficiency and can sequester OC effectively, contributing an important carbon sink to the sedimentary carbon cycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Temporal relationships among lunar crustal rocks

Temporal relationships among the three most common suites of lunar crustal rocks have been investigated by obtaining new high precision ages on Felsic/Alkali-suite Quartz monzodiorite Clast B from breccia 15405 and Magnesian-suite norite 78235/6/8/55/56 and comparing them to previously dated ferroan anorthosite sample 60025. The weighted average age of 4337.19 ± 0.49 Ma of 15405 Clast B is defined by zircon U-Pb and Pb-Pb ages as well as mineral isochron Sm-Nd and Nd-Nd ages. It is identical to the weighted average age for Apollo 17 norite 78235/6/8/55/56 of 4334.1 ± 3.5 Ma which is defined by Pb-Pb ages measured on baddeleyites in this investigation and less precise Pb-Pb and Sm-Nd ages reported in the literature. Both ages are ∼ 25 Ma younger than the weighted average of Sm-Nd and Pb-Pb ages reported in the literature on ferroan anorthosite 60025 of 4359.3 ± 2.3 Ma. The fact that ages of all three samples are defined by multiple U-Pb, Pb-Pb, Sm-Nd, and 142 Nd- 143 Nd chronometers provide confidence that they record the igneous crystallization history of the samples and do not represent disturbances or mixing lines with no temporal significance. Here, the extent to which these three ages represent broader scale magmatism is difficult to evaluate. Nevertheless, the age defined for 15405 Clast B, 78235/6/8/55/56, and 60025 are contemporaneous with the peak of ages observed in detrital zircons from the Apollo 12, 14, 15, and 17 landing sites (4340 ± 20 Ma), a Mg-suite Sm-Nd whole rock isochron defined by samples from Apollo 14, 15, 16, and 17 landing sites (4348 ± 25 Ma), and a Ferroan Anorthosite-suite Sm-Nd whole rock isochron defined by samples from the Apollo 15 and 16 landing sites (4354 ± 29 Ma). This implies that Ferroan Anorthosite-suite magmatism is temporally distinct and earlier than magmatism associated with the Mg-suite and the Felsic/Alkali-suite, as predicted by the lunar magma ocean model of lunar differentiation. The short 35 ± 10 Ma interval between primary ferroan anorthosite magmatism and secondary magmatism suggests that the lunar crust formed over a limited period of time. Although heat from decay of long-lived isotopes, large impacts, tidal heating associated with interactions between the Earth and Moon, and density driven overturn of the magma ocean have all been invoked to explain production of ancient secondary crustal magmatism, only tidal heating and cumulate overturn are consistent with the apparent short duration of secondary crustal magmatism and the great depth of crystallization implied for some Mg-suite samples. The initial ε 143 Nd values derived from the 15405 Clast B and 78238 Mg-suite norite isochrons, as well as a Mg-suite whole rock isochron are −0.23 ± 0.11, −0.27 ± 0.74, and −0.25 ± 0.09, respectively. They are identical within uncertainty indicating that Mg-suite and Felsic/Alkali-suite magmas were derived from materials that had the same time averaged Sm/Nd ratios since the formation of the solar system. This, combined with the contemporaneous nature of 15405 Clast B and 78235/6/8/55/56 Mg-suite norite, is consistent with evolution of both samples, and likely both magma suites, from a common source through closed system fractional crystallization or partial melting processes.

Chronology↗
DOE OSTI2026

How deep is your soil? Quantifying and spatially analyzing understudied deep soil in the United States

Deep soil is largely understudied and important in understanding biogeochemical processes in soil. Here, understudied soil is defined as the difference between soil studied to a known depth and the estimated bedrock depth. To understand more about deep soil, the understudied soil in the US was quantified and spatially analyzed using soil survey data and model estimates of bedrock depth. An equation was derived to find understudied soil using the dataset parameters “max lower depth studied”, “depth to bedrock”, and “likelihood of bedrock in the top 200 cm”. The survey data and bedrock model revealed that soil has been studied to an average depth of 1-2 meters, and the average depth to bedrock is 20 meters. Soil data density in the soil surveys was greatest in the West Coast, Midwest, and areas historically managed for agricultural, while the non-contiguous US and interior West were underrepresented. The soil had been studied deeper than the estimated soil depth in 455 out of 56,889 observation points concentrated in Alaska, California, Texas, Florida, Puerto Rico, and the US Virgin Islands. To understand the diversity and any taxonomic bias of the global soil data available, soil order was compared to US-based National Resource Conservation Service percentages and it was found that Oxisols, Alfisols, Ultisols, Andisols, and Histosols were overrepresented while Gelisols, Aridisols, Vertisols, Entisols, and Spodosols are underrepresented. Soil depth is important in exploring the complexity of biogeochemical processes that take place in soil.

Bedrock↗
DOE OSTI2025

Pseudoatom molecular dynamics plasma microfields

Spectral lines are powerful diagnostic tools for both laboratory and astrophysical plasmas, as their shape is sensitive to the plasma environment. The low-frequency component of the electric microfield is an important input for semi-analytic line broadening codes. Here, in this paper, we detail a new method of calculating plasma microfields using configuration-resolved pseudoatom molecular dynamics. This approach accounts for both quantum atomic structure and N-body effects, similar to density functional theory molecular dynamics, but with less computational cost. We present pseudoatom microfields at plasma conditions relevant for recent high energy density laboratory astrophysics experiments conducted at the Sandia Z-Machine, National Ignition Facility, and Linac Coherent Light Source. Compared to established microfield codes we find moderate deviations at solid density conditions and strong agreement at lower plasma densities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Elevated-Mn ChemCam targets illuminating Mn redox cycling and diagenesis in the Bradbury Rise, Gale Crater, Mars

Manganese plays a crucial role as a paleo-environmental and geological indicator due to its sensitivity to redox potential and pH variations in the environment. On Earth, the association between the rise of atmospheric oxygen during the Great Oxidation Event and the presence of Mn in the sedimentary rock record underscores its significance. Here, in this study, we reexamined ChemCam targets from the first 600 sols of the Mars Science Laboratory mission, focusing on identifying instances of above-average Mn within these targets. These elevated-Mn targets were categorized into distinct geologic classes, revealing a pattern linking heightened Mn levels with diagenetically altered materials, such as calcium-sulfate veins and concretions, as well as clay minerals within the same targets, indicating a compelling relationship between Mn enrichment and diagenetic processes. High concentrations of Mn were observed in chemically altered targets, suggesting the occurrence of multiple fluid events: the first to alter the material and the second to deposit Mn. The observed patterns suggest multiple diagenetic events and redox cycling that facilitated the deposition and transport of Mn subsequent to the initial dissolution of basaltic materials. This research sheds light on the complexity of martian diagenetic processes and their implications for the planet’s environmental evolution.

58 GEOSCIENCES↗
DOE OSTI2026

A decade of progress in understanding and managing legacy well integrity for geologic carbon storage

This study reviews a decade of research progress in legacy well integrity and risk management for geologic carbon storage (GCS) to commemorate the 20 th anniversary of the Intergovernmental Panel on Climate Change’s 2005 Special Report on Carbon Capture and Storage. In the past ten years, legacy well research has benefited from global efforts to constrain emissions from abandoned oil and gas wells, a continued focus on well materials performance in the presence of CO 2 -rich fluids, and practical experience gained through GCS implementation. Field measurements of abandoned well emissions show that leakage is not universal or catastrophic but forms a continuum of low-to-moderate fluxes that depend on isolation integrity and environmental attenuation. Materials research has constrained the conditions under which Portland cements exhibit self-sealing and non-sealing behaviors, and has identified the impact of geomechanical properties, non-uniform pathway apertures, multi-phase flow, and impurities in the CO 2 stream, on leakage pathways as important new areas for investigation. GCS projects at brownfield sites have inspired the creation of new workflows that integrate various tools and technologies to manage legacy well leakage risks. GCS implementation has also motivated a push towards scenario-based well modeling that directly informs permit applications. These advances inspire new research questions for the coming decade, particularly around the level of legacy well leakage risk that is environmentally acceptable and tolerable to stakeholders when sequestering millions of tonnes of CO 2 annually.

Carbon capture and storage↗
DOE OSTI2024

The impact of hot-press conditions on the durability of polymer electrolyte membrane fuel cells

The proton exchange membrane integrity can be compromised during hot-press fabrication of membrane electrode assemblies (MEAs) causing premature cell failures during operation. In this work, infrared (IR) thermography was used as a diagnostic tool to spatially visualize hydrogen (H 2 ) crossover and identify process-induced-membrane irregularities (PIMs). These irregularities were identified as seed locations for MEA failures. Fine tuning of hot-press conditions was used to mitigate premature cell failures informed by accelerated stress testing (AST). The impact of PIMs on the initial performance, high-frequency resistances, open-circuit voltage, and H 2 crossover are reported. Nafion XL and 212 membranes, hot-pressed with a force of 16 kg/cm 2 and temperature of 120°C, were found to be consistently irregularity-free. Irregularity-free MEAs using Nafion 211, 212, and XL membranes demonstrated AST lifetime improvements of 58, 64 and 400%, respectively, compared to those fabricated with non-optimized conditions. In conclusion, this work highlights the importance of fabrication parameters on premature cell failures.

08 HYDROGEN↗
DOE OSTI2025

Multi-plane moment-of-fluid interface reconstruction in 3D

Moment-of-fluid (MOF) methods for interface reconstruction approximate the region occupied by material in each mesh element only through reference to its geometric moments. Here, we present a 3D MOF method that represents the material (POM) in each cell as the convex intersection of the cell and multiple half-spaces, each selected to minimize the least-squares error between computed moments of the approximated material and provided reference moments. This optimization problem is highly non-linear and non-convex, making the numerical result very sensitive to the initial guess. To create an effective initial guess in each cell, we construct an ellipsoid from 0th–2nd order reference moments such that its shape corresponds with that of the POM. Within this ellipsoid we inscribe a polyhedron, and initialize the minimization problem with the half-spaces defined by each of its faces. The inscribed polyhedron has minimally 4 faces, and using up to 3rd order moments permits optimization over up to 20 unknown values. We therefore define MOF methods that utilize 4, 5, or 6 half-spaces, correspondingly initialized with the faces of a single inscribed tetrahedron, triangular prism, or hexahedron. Stability of the non-linear optimization is further improved with a prepossessing step that normalizes the reference moments according to the axes of the reference ellipsoid. Using this approach, the non-linear least-squares solver reliably converges to a near-global minimum from a single initial guess. We demonstrate accuracy and robustness using single-cell and multi-cell examples over a wide spectrum of geometry. In particular, we demonstrate our ability to exactly reproduce several important and complex features defined by up to four half-spaces, such as corners, filaments, filament tips, and embedded material in the cell.

3D interface reconstruction↗
DOE OSTI2026

Equilibrium spin polarization arising from chirality

Chirality-induced spin selectivity (CISS) describes how chiral molecules and materials generate spin polarization even at thermal equilibrium. This observation has challenged established principles of microscopic reversibility and Onsager reciprocity. We resolve this paradox by formulating a pseudo-Hermitian quantum framework that separates thermodynamic equilibrium from time-reversal symmetry. Within this approach, structural chirality and electron correlations, irrespective of their microscopic origin, are sufficient to produce CISS observables. Chirality enters through a non-local metric η that couples spin and spatial motion, leading to real spectra, unitary evolution, and thermodynamic consistency. The framework predicts a chirality-induced spin magnetic ordering characterized by a spin-displacement order, which reconciles equilibrium spin polarization with detailed balance and explains the persistence of CISS in materials composed of light elements. We derive generalized Onsager-Casimir relations that respect the observed CISS symmetry, i.e., parity $(\mathscr{P})$-odd and time-reversal $(\mathscr{T})$-odd, but exhibiting $\mathscr{PT}$-even symmetry. This approach establishes a coherent foundation for equilibrium CISS and provides a route to link chemical chirality with measurable spin-to-charge conversion effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Highly charged ion approach to measure nuclear charge radii of Fr, Ra, and Rn isotopes for precision measurements

Conventional electron scattering and muonic atom spectroscopy techniques are challenging to apply to heavy ions unless the element has at least one stable or extremely long-lived isotope. To overcome this limitation, a recently introduced method for determining nuclear charge radii relies on extreme-ultraviolet (EUV) spectroscopy of the D 1 line in highly charged Na-like ions. In this work, we present an experimental approach to measure the nuclear charge radii of isotopes of radioactive elements such as Fr, Ra, and Rn using this method at TRIUMF’s Ion Trap for Atomic and Nuclear Science (TITAN) instrument, located at the ISAC radioactive beam facility. We also explore the potential of using Na-like D 2 lines in the soft x-ray region for future measurements.

Electron beam ion trap↗
DOE OSTI2025

Aerosols are critical for “nowcasting” climate

Here, aerosols are a critical part of the earth’s energy balance and are changing rapidly. Monitoring Earth’s energy budget, clouds and aerosols with a coherent observing system from the ground, air and space, coupled with predictive models is necessary for timely, updated climate ‘nowcasts.’

Gettelman, Andrew [Pacific Northwest National Labo↗
DOE OSTI2025

Flavor Dependence of Charged Pion Fragmentation Functions

We have measured the flavor dependence of multiplicities for π + and π - production in semi-inclusive deep-inelastic scattering (SIDIS) on proton and deuteron to explore a possible charge symmetry violation in fragmentation functions. The experiment used an electron beam with energies of 10.2 and 10.6 GeV at Jefferson Lab and the Hall-C spectrometers. The electron kinematics spanned the range 0.3 < x < 0.6, 2 < Q 2 < 5.5 GeV2, and 2.2 < W < 3.2 GeV. The pion fractional momentum range was 0.3 < z < 0.7, and the transverse momentum range was 0 < p T < 0.25 GeV/c. Assuming factorization and allowing for isospin breaking, the results can be described by two “favored” and two “unfavored” effective low p T fragmentation functions that are flavor-dependent. We find each pair converges to a common flavor-independent fragmentation function at the highest W, where factorization is most applicable.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗