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.
A model describing the evaportion of an isolated drop of a multicomponent fuel containing hundreds of species has been developed. The model is based on Continuous Thermodynamics concepts wherein the composition of a fuel is statistically described using a Probability Distribution Function (PDF).
Combining multiple optical resonators or engineering dispersion of complex media has provided an effective method for demonstrating topological physics controlling photons in unprecedented ways such as unidirectional light propagation and spatially localized modes between an interface or on a corner. Further, adding nonlinear responses to those topological photonic systems has enabled achieving diverse phases of photons in both space and time, allowing for more functionalities in photonic devices that provide a new playground for studying dynamic features of nonlinear topological systems. However, most methods for describing nonlinear topological photonic systems rely on linear topological theories, making it challenging to accurately characterize the topology of nonlinear systems. Thus, substantial efforts have focused on rigorously describing nonlinear topological phases and developing effective tools to analyze nonlinear topological effects. Meanwhile, coupled multimode optical waveguides with nonlinear dynamic responses provide an excellent platform for the statistical description of photons, opening a new paradigm called “optical thermodynamics”. This review will introduce the basic concepts of nonlinear topological photonics and the recent development of theoretical approaches focusing on data-driven approaches for creating phase diagrams as well as the spectral localizer framework and the pseudospectrum method for understanding optical nonlinearities in topological systems. In addition, the new concept of optical thermodynamics will be introduced with some recent theoretical works.
Due to the failure of thermodynamics for low temperature near-extremal black holes, it has long been conjectured that a 'thermodynamic mass gap' exists between an extremal black hole and the lightest near-extremal state. For non-supersymmetric near-extremal black holes in Einstein gravity with an AdS2 throat, no such gap was found. Rather, at that energy scale, the spectrum exhibits a continuum of states, up to non-perturbative corrections. In this paper, we compute the partition function of near-BPS black holes in supergravity where the emergent, broken, symmetry is PSU(1, 1|2). To reliably compute this partition function, we show that the gravitational path integral can be reduced to that of a $\mathcal{N}=4$ supersymmetric extension of the Schwarzian theory, which we define and exactly quantize. In contrast to the non-supersymmetric case, we find that black holes in supergravity have a mass gap and a large extremal black hole degeneracy consistent with the Bekenstein–Hawking area. Furthermore, our results verify a plethora of string theory conjectures, concerning the scale of the mass gap and the counting of extremal micro-states.
The emergence of scanning probe and electron beam imaging techniques has allowed quantitative studies of atomic structure and minute details of electronic and vibrational structure on the level of individual atomic units. These microscopic descriptors, in turn, can be associated with local symmetry breaking phenomena, representing the stochastic manifestation of the underpinning generative physical model. In this work, we explore the reconstruction of exchange integrals in the Hamiltonian for a lattice model with two competing interactions from observations of microscopic degrees of freedom and establish the uncertainties and reliability of such analysis in a broad parameter-temperature space. In contrast to other approaches, we specifically specify a loss function inherent to thermodynamic systems and utilize it to estimate uncertainty in simulated realizations of different models. As an ancillary task, we develop a machine learning approach based on histogram clustering to predict phase diagrams efficiently using a reduced descriptor space. We further demonstrate that reconstruction is possible well above the phase transition and in the regions of parameter space when the macroscopic ground state of the system is poorly defined due to frustrated interactions. This suggests that this approach can be applied to the traditionally complex problems of condensed matter physics such as ferroelectric relaxors and morphotropic phase boundary systems, spin and cluster glasses, and quantum systems once the local descriptors linked to the relevant physical behaviors are known.
Here we present a mechanistic model of carbon nucleation and growth from a fluid at elevated temperature (T) and pressure conditions, typical of those found in the shallow Earth’s lithosphere. Our model uses a replica exchange reactive molecular dynamics framework in which molecular configurations are swapped between adjacent T replica at regular intervals according to underlying statistical mechanics. This framework allows predicting complex molecular structures and thermodynamics while remaining computationally efficient. Here we simulate the reactivity of an unstable mixture of CO 2 and CH 4 at 1000 K and 1 GPa. We find that the path to thermodynamic equilibrium is initially entropy-driven, producing a diversity of short-lived species, including various alcohols with intermediate carbon oxidation states. Cyclic and polycyclic radicals that are sometimes resonance-stabilized form next and set the stage for carbon nucleation. The carbon exsolution process releases abundant water, is exothermic and starts with the nucleation of a large aggregate of hydrogenated graphene flakes from covalently bonded polycyclic units. The carbon backbone of this nucleus subsequently grows into a hydrogen-depleted fullerene-like structure, before evolving toward a partially bilayered graphene layer. Overall, our results show that the mechanism of graphitic C formation is certainly not bimolecular, and that it may involve a combination of key condensation and radical chain reactions. This will help understand the isotopic, and reactive characteristics of carbon-bearing fluids during their upward transit through the Earth’s mantle and crust. Moreover, the mechanistic insights outlined here present intriguing similarities with the process of soot and interstellar dust formation, which suggests that the widespread distribution of abiotic polyaromatic and graphitic material on Earth and beyond may reflect the prevalence of a fundamental chemical pathway.
The Jarzynski equality allows the calculation of free-energy differences using values of work measured from nonequilibrium trajectories. The number of trajectories required to accurately estimate free-energy differences in this way grows sharply with the size of work fluctuations, motivating the search for protocols that perform desired transformations with minimum work. However, protocols of this nature can involve varying temperature, to which the Jarzynski equality does not apply. Here, we derive a variant of the Jarzynski equality that applies to varying-temperature protocols, and show that it can have better convergence properties than the standard version of the equality. We derive this modified equality and the associated fluctuation relation within the framework of Markovian stochastic dynamics, complementing related derivations done within the framework of Hamiltonian dynamics.
The Community Climate Model (CCM) is a comprehensive three-dimensional global atmospheric model assembled at the National Center for Atmospheric Research. Detailed comparisons are made between the climate simulated by a seasonal version of CCM, namely CCM1, at 12 level, R15 spectral resolution, and that revealed by ECMWF operational analyses over 1980-86 truncated to a similar resolution. A variety of circulation statistics are presented in order to reveal the spatial character and seasonality of CCM1 biases in temperature, winds, and wave flux quantities. The biases revealed for CCM1 are similar in many respects to those found in other climate models run at similar resolution. Results of this study are then presented, discussing temperature biases and thermodynamic balances, wind fields, zonal mean wave diagnostics, stationary waves, and local eddy statistics. The major findings are then summarized.
The Global Reference Atmospheric Model is used along with the revised perturbation statistics to evaluate and computer graph various atmospheric statistics along a space shuttle reference mission and abort trajectory. The trajectory plots are height vs. ground range, with height from ground level to 155 km and ground range along the reentry trajectory. Cross sectional plots, height vs. latitude or longitude, are also generated for 80 deg longitude, with heights from 30 km to 90 km and latitude from -90 deg to +90 deg, and for 45 deg latitude, with heights from 30 km to 90 km and longitudes from 180 deg E to 180 deg W. The variables plotted are monthly average pressure, density, temperature, wind components, and wind speed and standard deviations and 99th inter-percentile range for each of these variables.
Explore the source record for details and available documents.
Abstract Gaining continued insights into the impact of global warming on the occurrence of hurricane-associated intense record downpours is essential for building climate resilient communities. This study investigates projected future changes in extreme rainfall over the Northeast United States, as represented by extreme daily amounts during Hurricane Ida in 2021. We used historical control simulations of Weather Research and Forecasting (WRF) Model generated from 40 years of weather events (1980–2014, 12 km) forced by the fifth generation European Centre for Medium-Range Weather Forecasts atmospheric reanalysis. These simulations are thermodynamically modified (2060–2100) via an imposed warming for the high-emission scenario of shared socioeconomic pathway (SSP585) from a range of general circulation models. Ground observations from the Global Historical Climatology Network (1950–2014) and WRF simulations (historical, 1980–2014, and future, 2060–2100) are integrated into a nonstationary generalized extreme value (GEV) framework to assess the frequency of Ida’s heaviest daily rain rates under the SSP585 scenario. Results show that Ida’s daily maximum rainfall recorded at different observation locations was higher than the single highest September daily maximum observed (1950–2014) for 5 out of 17 stations (∼30% of the stations). Ida-like extreme daily rain rates are projected to be, on average, more than 2 times more likely to occur at the end of the century in the simulations (with some regions as high as 5 times). This work demonstrates that integrating a high-resolution atmospheric model’s present-day and thermodynamically modified future simulations along with ground observations, within a nonstationary statistical framework, is crucial for understanding changing characteristics of extreme weather events. Significance Statement Daily scale extreme precipitation is expected to become more frequent and severe, as evidenced by observations and model simulations. While it is important to investigate how these intensifying heavy rainfall events affect current engineering standards, fewer studies have contextualized how warming impacts the most extreme rainfall from a single storm event relative to historical heavy downpours. In this study, we focused on the daily extreme rainfall associated with the extratropical transition of Hurricane Ida (2021), particularly over the northeastern United States—some of which exceeded the commonly used hydrologic design criteria for a 100-yr storm. Using a high-resolution atmospheric model simulation, we investigated how continued warming may influence the frequency of such daily rain rates. Under a high-emission scenario, these events are projected to become up to 5 times more likely at the end of the twenty-first century.
In general, there are two broad scientific objectives when using cloud resolving models (CRMs or cloud ensemble models-CEMs) to study tropical convection. The first one is to use them as a physics resolving models to understand the dynamic and microphysical processes associated with the tropical water and energy cycles and their role in the climate system. The second approach is to use the CRMs to improve the representation of moist processes and their interaction with radiation in large-scale models. In order to improve the credibility of the CRMs and achieve the above goals, CRMs using identical initial conditions and large-scale influences need to produce very similar results. Two CRMs produced different statistical equilibrium (SE) states even though both used the same initial thermodynamic and wind conditions. Sensitivity tests to identify the major physical processes that determine the SE states for the different CRM simulations were performed. Their results indicated that atmospheric horizontal wind is treated quite differently in these two CRMs. The model that had stronger surface winds and consequently larger latent and sensible heat fluxes from the ocean produced a warmer and more humid modeled thermodynamic SE state. In addition, the domain mean thermodynamic state is more unstable for those experiments that produced a warmer and more humid SE state. Their simulated wet (warm and humid) SE states are thermally more stable in the lower troposphere (from the surface to 4-5 km in altitude). The large-scale horizontal advective effects on temperature and water vapor mixing ratio are needed when using CRMs to perform long-term integrations to study convective feedback under specified large-scale environments. In addition, it is suggested that the dry and cold SE state simulated was caused by enhanced precipitation but not enough surface evaporation. We find some problems with the interpretation of these three phenomena.
Metallic glasses derive their properties from the statistics of local atomic motifs rather than from long-range order, yet a quantitative, chemistry-specific link between motif populations and the underlying glassy state has remained elusive. In this work we combine large-scale molecular dynamics, Voronoi tessellation, deep neural networks, and SHapley Additive exPlanations (SHAP) to identify which local structural motifs define the glassy state of Cu—Zr metallic glasses. A dataset of 17,180 atomistic configurations spanning ten compositions (Cu 20 Zr 80 –Cu 80 Zr 20 ) and four quench rates (10 9 –10 12 K/s) is used to train a feed-forward neural network that regresses temperature across the 50–2000 K liquid–supercooled–glass range, achieving a mean absolute error of 19.89 K and R 2 = 0.9974, confirming that the local structural state is faithfully encoded in motif-level structure. SHAP analysis then reveals that a tightly coupled near-icosahedral family of motifs (coordination numbers (CN) 11–13, including the full icosahedron 001200 and its single-atom-perturbation sibling 10930) collectively encodes the thermodynamic state of the system across the full glass-forming range. The CN = 11–13 ordered members carry negative SHAP values at high populations, tracking the most deeply-quenched configurations, while 10930 shows the reversed signature consistent with its role as a soft-spot host whose population shrinks as the icosahedral network deepens. The analysis demonstrates that explainable machine learning can isolate the minimal motif vocabulary defining the glassy state and recovers the near-icosahedral building blocks previously identified by data-driven analyses of Cu—Zr. The approach provides a general, chemistry-specific route for characterizing the structural state of disordered materials.
Calculating thermodynamic potentials and observables efficiently and accurately is key for the application of statistical mechanics simulations to materials science. However, naive Monte Carlo approaches, on which such calculations are often dependent, struggle to scale to complex materials in many state-of-the-art disciplines such as the design of high entropy alloys or multi-component catalysts. To address this issue, we adapt sampling tools built upon machine learning-based generative modeling to the materials space by transforming them into the semi-grand canonical ensemble. Furthermore, we show that the resulting models are transferable across wide ranges of thermodynamic conditions and can be implemented with any internal energy model U, allowing integration into many existing materials workflows. We demonstrate the applicability of this approach to the simulation of benchmark systems (AgPd, CuAu) that exhibit diverse thermodynamic behavior in their phase diagrams. Finally, we discuss remaining challenges in model development and promising research directions for future improvements.
In internal combustion engine research, cylinder pressure measurements provide valuable information about the underlying thermodynamic and combustion processes, and are typically collected in ensembles of several 100 traces. Although in some particular fields of combustion research all traces are analyzed, in most cases only one trace is studied because analyzing all the traces is impractical due to the large number of collected samples. Instead, an ensemble-averaged pressure trace is commonly calculated and used for analysis. However, this pressure trace is highly smoothed and dynamic information is lost during the averaging process. With the average trace, pressure rise rates are lower and pressure oscillations such as the ones resulting from combustion knock are lost. In this work, a statistical method was developed to determine the “most representative cycle,” which is the cycle from the ensemble that has the pressure trace most representative of the engine operating condition. Eleven characteristic parameters are computed from each pressure trace and probabilistic distributions are obtained for each of the parameters using all the traces in the ensemble. Finally, the most representative cycle is selected by means of a cost function minimization. The benefits of this method are illustrated using experimental data from four very different engine platforms, under four different combustion modes and over a range of operating conditions.
While classical size effects usually lead to a reduced effective thermal conductivity, we report here that nonthermal phonon populations produced by a micro/nanoscale heat source can lead to enhanced heat conduction, exceeding the prediction from Fourier's law. We study nondiffusive thermal transport by phonons at small distances within the framework of the Boltzmann transport equation (BTE) and demonstrate that the transport is significantly affected by the distribution of phonons emitted by the source. We discuss analytical solutions of the steady-state BTE for a source with a sinusoidal spatial profile, as well as for a three-dimensional Gaussian “hot spot,” and provide numerical results for single crystal silicon at room temperature. If a micro/nanoscale heat source produces a thermal phonon distribution, it gets hotter than that predicted by the heat diffusion equation; however, if the source predominantly produces low-frequency acoustic phonons with long mean free paths, it may get significantly cooler than that predicted by the heat equation, yielding an enhanced heat transport beyond bulk heat conduction.
Review of the literature on the nature and extent of hydrogen interactions with metals and the role of hydrogen in metal failure. The classification of hydrogen-containing systems is discussed, including such categories as covalent hydrides, volatile hydrides, polymeric hydrides, and transition metal hydride complexes. The use of electronegativity as a correlating parameter in determining hydride type is evaluated. A detailed study is made of the thermodynamics of metal-hydrogen systems, touching upon such aspects as hydrogen solubility, the positions occupied by hydrogen atoms within the solvent metal lattice, the derivation of thermodynamic functions of solid solutions from solubility data, and the construction of statistical models for hydrogen-metal solutions. A number of theories of hydrogen-metal bonding are reviewed, including the rigid-band model, the screened-proton model, and an approach employing the augmented plane wave method to solve the one-electron energy band problem. Finally, the mechanism of hydrogen embrittlement is investigated on the basis of literature data concerning stress effects and the kinetics of hydrogen transport to critical sites.
The transformation of cold neutral intergalactic hydrogen into a highly ionized warm plasma marks the end of the cosmic dark ages and the beginning of the age of galaxies. The details of this process reflect the nature of the early sources of radiation and heat, the statistical characteristics of the large-scale structure of the Universe, the thermodynamics and chemistry of cosmic baryons, and the histories of star formation and black hole accretion. A number of massive data sets from new ground- and space-based instruments and facilities over the next decade are poised to revolutionize our understanding of primeval galaxies, the reionization photon budget, the physics of the intergalactic medium (IGM), and the fine-grained properties of hydrogen gas in the "cosmic web". In this review we survey the physics and key aspects of reionization-era modeling and describe the diverse range of computational techniques and tools currently available in this field.