SEARCH · Engineering Papers
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.
We seek to answer the following questions: What interactions across spatial and temporal scales drive multivariate compound and cascading water cycle extremes, and how can we leverage understanding of such interactions to improve predictive understanding of these events and quantification of their impacts?
Plutonium metal exhibits an anomalously large softening of its bulk modulus at elevated temperatures that is made all the more extraordinary by the finding that it occurs irrespective of whether the thermal expansion coefficient is positive, negative, or zero—representing an extreme departure from conventional Grüneisen scaling. We show here that the cause of this softening is the compressibility of plutonium’s thermally excited electronic configurations, which has thus far not been considered in thermodynamic models. We show that when compressible electronic configurations are thermally activated, they invariably give rise to a softening of the bulk modulus regardless of the sign of their contribution to the thermal expansion. The electronically driven softening of the bulk modulus is shown to be in good agreement with elastic moduli measurements performed on the gallium-stabilized δ phase of plutonium over a range of temperatures and compositions and is shown to grow rapidly at small concentrations of gallium and at high temperatures, where it becomes extremely sensitive to hydrostatic pressure.
Extreme heat has been associated with adverse health outcomes, yet its impact on asthma exacerbations remains understudied. This is, in part, due to data limitations: research that relies on weather station records and aggregated health statistics cannot resolve fine-scale differences in heat impacts. This study investigates the association between heat wave definitions and summertime asthma-related emergency department visits in Baltimore, Maryland from 2016 to 2022, including 819 adult and 695 pediatric exacerbations. Using geocoded electronic health records and air temperature measurements at several spatial resolutions, we applied a case-crossover design with conditional logistic regressions at the census block group and tract levels. We found strong associations between asthma exacerbations and nighttime heat wave definitions based on relative thresholds of minimum temperatures when census block group or tract level temperature estimates were used. These relationships were significant for both age groups and showed elevated risks in socially vulnerable areas. In contrast, heat wave definitions derived from the city's primary National Weather Service synoptic weather station show associations between asthma and daytime heat extremes, suggesting that the character of the heat hazard depends on the scale at which it is defined. The extreme heat event definition used by Baltimore City's Code Red system showed no significant association with exacerbations. These findings highlight the importance of data resolution in shaping health inferences related to extreme heat in urban environments. Further, this study demonstrates that, regardless of spatial scale, extreme heat is associated with asthma exacerbations in both age groups.
Modeling effective transport properties of 3D porous media, such as permeability, at multiple scales is challenging as a result of the combined complexity of the pore structures and fluid physics—in particular, confinement effects which vary across the nanoscale to the microscale. While numerical simulation is possible, the computational cost is prohibitive for realistic domains, which are large and complex. Although machine learning (ML) models have been proposed to circumvent simulation, none so far has simultaneously accounted for heterogeneous 3D structures, fluid confinement effects, and multiple simulation resolutions. By utilizing numerous computer science techniques to improve the scalability of training, we have for the first time developed a general flow model that accounts for the pore-structure and corresponding physical phenomena at scales from Angstrom to the micrometer. Using synthetic computational domains for training, our ML model exhibits strong performance (R 2 = 0.9) when tested on extremely diverse real domains at multiple scales.
Modeling effective transport properties of 3D porous media, such as permeability, at multiple scales is challenging as a result of the combined complexity of the pore structures and fluid physics—in particular, confinement effects which vary across the nanoscale to the microscale. While numerical simulation is possible, the computational cost is prohibitive for realistic domains, which are large and complex. Although machine learning (ML) models have been proposed to circumvent simulation, none so far has simultaneously accounted for heterogeneous 3D structures, fluid confinement effects, and multiple simulation resolutions. By utilizing numerous computer science techniques to improve the scalability of training, we have for the first time developed a general flow model that accounts for the pore-structure and corresponding physical phenomena at scales from Angstrom to the micrometer. Using synthetic computational domains for training, our ML model exhibits strong performance (R 2 = 0.9) when tested on extremely diverse real domains at multiple scales.
Resolving the most fundamental questions in cosmology requires simulations that match the scale, fidelity, and physical complexity demanded by next-generation sky surveys. To achieve the realism needed for this critical scientific partnership, detailed gas dynamics must be treated self-consistently with gravity for end-to-end modeling of structure formation. Exascale computing enables simulations that span survey-scale volumes while incorporating key astrophysical processes that shape complex cosmic structures. We present results from CRK-HACC, a cosmological hydrodynamics code built for extreme scalability. Using separation-of-scale techniques, GPU-resident tree solvers, in situ analysis pipelines, and multi-tiered I/O, CRK-HACCexecuted Frontier-E: a four trillion particle full-sky simulation, over an order of magnitude larger than previous efforts. The run achieved 513.1 PFLOPs peak performance, processing 46.6 billion particles per second and writing more than 100 PB of data in just over one week of runtime. Frontier-E marks a significant advance in predictive modeling for next-generation cosmological science.
A carbon fiber rope thermal barrier is being evaluated as a replacement for the conventional room temperature vulcanizing (RTV) thermal barrier that is currently used to protect o-rings in Reusable Solid Rocket Motor (RSRM) nozzle joints. Performance requirements include its ability to cool any incoming, hot propellant gases that fill and pressurize the nozzle joints, filter slag and particulates, and to perform adequately in various joint assembly conditions as well as dynamic flight motion. Modified National Aeronautics and Space Administration (MNASA) motors, with their inherent and unique ability to replicate select RSRM internal environment features, were an integral step in the development path leading to full scale RSRM static test demonstration of the carbon fiber rope (CFR) joint concept. These 1/4 scale RSRM motors serve to bridge the gap between the other classes of subscale test motors (extremely small and moderate duration, or small scale and short duration) and the critical asset RSRM static test motors. A series of MNASA tests have been used to demonstrate carbon fiber rope performance and have provided rationale for implementation into a full-scale static motor and flight qualification.
Assessing energy resources (e.g., solar, wind, and hydro) under future scenarios requires datasets with sufficient spatial and temporal detail to capture variability and extreme events. While global-scale Earth System Model (ESM) projections are widely used, their coarse resolution limits direct application to regional energy system analyses. Dynamical downscaling offers a robust approach to generate physically consistent, fine-scale datasets that better represent local atmospheric processes impacting energy resources. In this work, we present a two-stage approach for producing high-resolution historical and future projections over the contiguous United States (CONUS). First, we optimize the Weather Research and Forecasting (WRF) model configuration for energy-relevant variables - solar irradiance, wind speed, and precipitation - by conducting ERA5-driven simulations at 8-km and 28-km resolution. Multiple physics schemes and model configurations within the WRF are evaluated against observational datasets including the National Solar Radiation Database (NSRDB), the Parameter-elevation Regressions on Independent Slopes Model (PRISM), and the Stage IV multi-radar/multi-sensor precipitation product for the CONUS domain. Using the best-performing configuration, we dynamically downscale MPI-ESM1-2-HR simulations for 2000-2060 under SSP2-4.5 and SSP5-8.5 scenarios at 4-km spatial and hourly temporal resolution. This presentation will provide a comprehensive analysis of the results from multiple numerical experiments and high-resolution ESM projections. In addition, we will discuss potential applications of our high-resolution datasets within the energy sector and outline future research avenues dedicated to evaluating how extreme weather events influence system performance and resilience.
Abstract Extreme rainfall events drive the amount and spatial distribution of rainfall in the Amazon and are a key driver of forest dynamics across the basin. This study investigates how the 3-hourly predictions in the High Resolution Model Intercomparison Project (HighResMIP, a component of the recent Coupled Model Intercomparison Project, CMIP6) represent extreme rainfall events at annual, seasonal, and sub-daily time scales. TRMM 3B42 (Tropical Rainfall Measuring Mission) 3 h data were used as observations. Our results showed that eleven out of seventeen HighResMIP models showed the observed association between rainfall and number of extreme events at the annual and seasonal scales. Two models captured the spatial pattern of number of extreme events at the seasonal and annual scales better (higher correlation) than the other models. None of the models captured the sub-daily timing of extreme rainfall, though some reproduced daily totals. Our results suggest that higher model resolution is a crucial factor for capturing extreme rainfall events in the Amazon, but it might not be the sole factor. Improving the representation of Amazon extreme rainfall events in HighResMIP models can help reduce model rainfall biases and uncertainties and enable more reliable assessments of the water cycle and forest dynamics in the Amazon.
It is well known that extreme or prolonged rainfall is the dominant trigger of landslides; however, there remain large uncertainties in characterizing the distribution of these hazards and meteorological triggers at the global scale. Researchers have evaluated the spatiotemporal distribution of extreme rainfall and landslides at local and regional scale primarily using in situ data, yet few studies have mapped rainfall-triggered landslide distribution globally due to the dearth of landslide data and consistent precipitation information. This research uses a newly developed Global Landslide Catalog (GLC) and a 13-year satellite-based precipitation record from Tropical Rainfall Measuring Mission (TRMM) data. For the first time, these two unique products provide the foundation to quantitatively evaluate the co-occurence of precipitation and rainfall-triggered landslides globally. The GLC, available from 2007 to the present, contains information on reported rainfall-triggered landslide events around the world using online media reports, disaster databases, etc. When evaluating this database, we observed that 2010 had a large number of high-impact landslide events relative to previous years. This study considers how variations in extreme and prolonged satellite-based rainfall are related to the distribution of landslides over the same time scales for three active landslide areas: Central America, the Himalayan Arc, and central-eastern China. Several test statistics confirm that TRMM rainfall generally scales with the observed increase in landslide reports and fatal events for 2010 and previous years over each region. These findings suggest that the co-occurrence of satellite precipitation and landslide reports may serve as a valuable indicator for characterizing the spatiotemporal distribution of landslide-prone areas in order to establish a global rainfall-triggered landslide climatology. This research also considers the sources for this extreme rainfall, citing teleconnections from ENSO as likely contributors to regional precipitation variability. This work demonstrates the potential for using satellite-based precipitation estimates to identify potentially active landslide areas at the global scale in order to improve landslide cataloging and quantify landslide triggering at daily, monthly and yearly time scales.
Energy equipartition is a powerful theoretical tool for understanding astrophysical plasmas. It is invoked, for example, to measure magnetic fields in the interstellar medium (ISM), as evidence for small-scale turbulent dynamo action, and, in general, to estimate the energy budget of star-forming molecular clouds. Here, in this study, we motivate and explore the role of the volume-averaged root-mean-squared (rms) magnetic coupling term between the turbulent, δB , and large-scale, B 0 , fields, $\langle$(δB·B 0 ) 2 $\rangle$$_{V}^{1/2}$. By considering the second moments of the energy balance equations we show that the rms coupling term is in energy equipartition with the volume-averaged turbulent kinetic energy for turbulence with a sub-Alfvénic large-scale field. Under the assumption of exact energy equipartition between these terms, we derive relations for the magnetic and coupling term fluctuations, which provide excellent, parameter-free agreement with time-averaged data from 280 numerical simulations of compressible magnetohydrodynamic (MHD) turbulence. Furthermore, we explore the relation between the turbulent mean field and total Alfvén Mach numbers, and demonstrate that sub-Alfvénic turbulence can only be developed through a strong, large-scale magnetic field, which supports an extremely super-Alfvénic turbulent magnetic field. This means that the magnetic field fluctuations are significantly subdominant to the velocity fluctuations in the sub-Alfvénic large-scale field regime. Throughout our study, we broadly discuss the implications for observations of magnetic fields and understanding the dynamics in the magnetized ISM.
In many applications, physical domains are geometrically complex making it challenging to perform coarse-scale approximation. A defeaturing process is often used to simplify the domain in preparation for approximation and analysis at the coarse scale. In this paper, a methodology is presented for constructing a coarse-scale reproducing basis on geometrically complex domains given an initial fine-scale mesh of the fully featured domain. The initial fine-scale mesh can be of poor quality and extremely refined. The construction of the basis functions begins with a coarse-scale covering of the domain and generation of weighting functions with local support. Manifold geodesics are used to define distances within the local support for general applicability to non-convex domains. Conventional moving least squares is used to construct the coarse-scale reproducing basis. Applications in quasi-interpolation and linear elasticity are presented.
Origami offers a distinct approach for designing and engineering new material structures and properties. The folding and stacking of atomically thin van der Waals (vdW) materials, for example, can lead to intriguing new physical properties including bandgap tuning, Van Hove singularity, and superconductivity. On the other hand, achieving well-controlled folding of vdW materials with high spatial precision has been extremely challenging and difficult to scale toward large areas. In this paper, a deterministic technique is reported to fold vdW materials at a defined position and direction using microfluidic forces. Electron beam lithography (EBL) is utilized to define the folding area, which allows precise control of the folding geometry, direction, and position beyond 100 nm resolution. Using this technique, single-atomic-layer vdW materials or their heterostructures can be folded without the need for any external supporting layers in the final folded structure. In addition, arrays of patterns can be folded across a large area using this technique and electronic devices that can reconfigure device functionalities through folding are also demonstrated. Such scalable formation of folded vdW material structures with high precision can lead to the creation of new atomic-scale materials and superlattices as well as opening the door to realizing foldable and reconfigurable electronics.
Technologies based on water adsorption such as water harvesting from air have tremendous potential in mitigating important global crises such as water scarcity. An important challenge to the deployment of such technologies is finding optimal adsorbent materials. Given the large materials space of available adsorbents, large-scale computational screening can be extremely helpful for this task. This work explores the methods and details associated with such screening procedures and recommends best practices. We also shed light on the limitations of traditionally used and inexpensive to compute prescreening approaches involving geometric and energetic features to predict water adsorption behavior of porous materials. Such approaches can provide general trends to predict adsorption behavior but may lead to the overlook of potentially important structures due to the complex nature of water adsorption. Finally, this study offers insights for future water adsorption simulations to facilitate the development of optimal water adsorbents.
An analytical framework is presented to predict the effects of sample size on the maximum value distribution (MVD) of the driving forces for fatigue crack formation in metals and alloys. The distribution of the maximum driving force for fatigue crack formation over the domain follows the generalized extreme value theory in the limit as the domain size increases to infinity. Here, a simulation-based analysis of microstructure influences on fatigue resistance for polycrystalline metals and alloys is very costly, and reaching those limits is intractable. This work models the MVD of Fatigue Indicator Parameters (FIPs), which serve as surrogate measures for the driving force for fatigue crack formation, at finite sample sizes prior to their convergence to a limiting extreme value distribution. Large-scale crystal plasticity finite element (CPFE) simulations of FCC Al 7075-T6 with microstructure realizations of various sizes are incorporated to calibrate and evaluate the developed framework, and a total of ∼6.5 million grains of Al 7075-T6 are examined. The calibrated analytical solution agrees well with the brute force Monte Carlo simulation results extracted from the CPFE simulations. Furthermore, the developed formulation can predict the MVD of FIPs for different sample sizes using a size-dependent parameter, and it is capable of accurately extrapolating the MVD of FIPs for much larger microstructure sample sizes than the size used for its calibration.
Using what has become a celebrated catchphrase, Philip W. Anderson once wrote that “more is different” (Science, Vol. 177, Issue 4047, pp. 393–396, 1972). First formulated in the context of condensed matter, this statement carries far beyond the sole limits of solid-state physics. It emphasizes that collective behavior can be more than the mere sum of what happens for elementary constituents or the mere collation of the evolution of each degree of freedom. Said otherwise, complex phenomena can arise out of the interplay between multiple sub-phenomena each of which can be relatively simple. The process of particle resuspension, in which discrete particles adhering on a surface are pulled off and carried away by a fluid flow, is another example involving a web of phenomena pertaining to fluid mechanics, particle dynamics and interface chemistry whose cross-effects create an intricate topic. The purpose of this review is to analyze the physics at play in particle resuspension in order to bring insights into the rich complexity of this common but challenging concern. Following the more-is-different vision, this is performed by starting from a range of practical observations and experimental data. We then work our way through the investigation of the key mechanisms which play a role in the overall process. In turn, these mechanisms reveal an array of fundamental interactions, such as particle–fluid, particle–particle and particle–surface, whose combined effects create the tapestry of current applications. At the core of this analysis are descriptions of these physical phenomena and the different ways through which they are intertwined to build up various models used to provide quantitative assessment of particle resuspension. The physics of particle resuspension implies to hold together processes occurring at extremely different space and time scales and models are key in providing a single vehicle to lead us through such multiscale journeys. This raises questions on what makes up a model and one objective of the present work is to clarify the essence of a modeling approach. In spite of its ubiquitous nature, particle resuspension is still at the early stages of developments. Many extensions need to be worked out and revisiting the art of modeling is not a moot point. The need to consider more complex objects than small and spherical particles and, moreover, to come up with unified descriptions of mono- and multilayer resuspension put the emphasis on solid model foundations if we are to go beyond current limits. Further, this is very much modeling in the making and new ideas are proposed to stimulate interest into this everyday but challenging issue in physics.
Understanding Li distribution in layered lithium transition-metal oxide (LiTMO) cathodes in Li-ion batteries has been a major challenge at the atomic scale and nanoscale. Li is extremely difficult to study by transmission electron microscopy (TEM) because the high-energy electrons impart significant energy and cause massive migration. Here, we directly map the intrinsic spatial distribution and bonding of Li in LiNiO 2 -layered cathode materials using low-dose and low-loss electron energy loss spectroscopy (EELS). EELS spectra of the Li–K edge are measured simultaneously with O–K and Ni–L, M 3,2 edges from layered, cation-mixed, and rock-salt phases and directly matched with atomic-resolution scanning TEM images to correlate the changes in peak intensities and positions to the stoichiometry changes with continual loss of Li and O. Changes in the Li content in the LiNiO 2 particles as a function of electron beam dose are studied by sequential Li spectroscopic mapping. We show that the “intrinsic” Li distribution can be observed using a total dose of less than ~1.5 × 10 8 e – nm –2 at an accelerating voltage of 80 kV. The method of nanoscale mapping of Li distribution introduced in this study is applicable to high-Ni LiTMO cathode materials (>89% of Ni) as well as LiNiO 2 . Further study on the extra peaks of the Li–K edge reveals that the peak at ~59 eV is from the Li ions intercalated in between NiO 2 layers barely interacting with each other with less Li K shell electrons pulled to the L shell electrons in NiO 2 . Furthermore, the results shown here provide improved low-loss TEM characterization approaches that can be used to understand the intrinsic fundamental behaviors in Li-ion batteries.
The 237 Np Mössbauer isomer shift and quadrupole splitting (QS) are powerful probes for the metal–ligand bonding of neptunium, a 5f-element of vital importance in the nuclear fuel cycle. A large set of Np compounds with different oxidation states (III) to (VII) is studied to investigate, by first-principles calculations, isomer shifts and the QS trends in relation to the Np oxidation state. Natural Bond Orbital analysis reveals that in addition to donation bonding to the 5f shell, participation of the 6d and 7s neptunium shells in covalent (donation) bonding substantially impacts the isomer shifts. The isomer shift cannot be interpreted solely by the 5f shell electron count. The isomer shift for Np(II) compounds is estimated to be in the range of 31–34 mm/s, less positive than for Np(III) compounds. For the QS, density functional calculations fail to reproduce the quadrupole splitting for some Np(VI) ionic solids. A multiconfigurational wave function approach reproduces the observed QS trends. The calculations give a semiquantitative interpretation of the trends for Np oxidation states (V) to (VII). Furthermore, the contrasting QS for standard and “reverse” neptunyl(VI), at the opposite extremes of the observed QS scale, arises predominantly from the different crystal environments.