Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “specific phase transitions”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Insight into pseudo-film-boiling in supercritical water under simplified interfacial conditions

The focus of this paper is the physics of supercritical water film boiling through the simulation of natural convection of supercritical water in a rectangular cavity. It has been widely believed that above a fluid's critical point, only one phase exists. However, recent research has indicated the existence of distinct gas-like and liquid-like regions separated by the Widom line, which is the locus of the maxima of the specific heat. Along the Widom line, density decreases sixfold, viscosity drops by a factor of two, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film-boiling, denoted here as pseudo-two-phase flow dynamics. When the cavity is heated from the bottom, a pseudo-film of gas-like fluid is formed at the heated wall, and the flow complexity evolves, characterized by the formation of multiple vortices. Results demonstrate the ability of modern conjugate heat transfer computational fluid dynamics to capture the pseudo two-phase regime, with the formation of a pseudo-film at the heated wall. Further, the boiling curve is studied, showing an increase in heat flux with an increase in the temperature differential between the wall and the pseudo-critical temperature. Furthermore, results also show a previously uncharacterized instability, where the convection cells in the cavity oscillate and can flip directions under certain conditions.

Computational fluid dynamics↗
DOE OSTI2025

End-of-Century Changes in Orographic Precipitation with the Intermediate Complexity Atmospheric Research Model over the Western United States

Abstract Downscaled precipitation projections were created using the Intermediate Complexity Atmospheric Research (ICAR) model over the western United States to increase the physical realism in orographic precipitation changes. End-of-century simulations from eight models in phase 5 of the Coupled Model Intercomparison Project (CMIP5) were downscaled with ICAR and compared to the widely utilized statistically downscaled dataset, localized constructed analogs (LOCAs), to understand where and why projections of cool-season (September–May) precipitation differed. ICAR and LOCA precipitation projections were similar, but their sign differed in hydrologically relevant regions likely due to ICAR’s simulation of microphysics and mesoscale dynamics with high-resolution topography (6 km). In the Pacific Northwest, cool-season precipitation projections from ICAR showed an increase on the windward side of the Cascades and no significant change within the lee. This difference between the windward and leeward side was attributed to reduced zonal wind speeds, allowing more time for microphysical processes within ICAR. This contrast is enhanced by rain’s faster fall speed compared to snow, limiting transport into the lee. Meanwhile, LOCA projected an increase in precipitation across the Cascades. In the Upper Colorado River basin, LOCA projected an increase in precipitation in high elevation regions (>3000 m), but ICAR projected no significant change or a decrease in precipitation. High elevation differences were most evident in the spring and fall and were also attributed to a snow-to-rain transition and dynamical processes that impacted orographic enhancement within ICAR. Idealized, controlled studies are needed to better isolate individual processes, but these results underscore the importance of including microphysics and mesoscale dynamics within regional-scale precipitation projections. Significance Statement A set of global climate model simulations was downscaled using an atmospheric model that contains key physical equations, referred to as Intermediate Complexity Atmospheric Research (ICAR). ICAR was used to examine projected changes in end-of-century cool-season precipitation over mountains in the western United States. Precipitation projections from ICAR were similar to projections that used statistical relationships to downscale climate projections. However, projections differed between ICAR and statistically downscaled datasets in whether they increased, decreased, or stayed the same in specific, hydrologically relevant regions such as the eastern Cascades and high elevation areas of the Upper Colorado River basin. These differences were attributed to the simulation of physical processes in ICAR. The results highlight the importance of kilometer-scale atmospheric processes in regional climate projections.

Currier, William Ryan [NOAA/Physical Sciences Labo↗
DOE OSTI2025

Carbonation and Sulfidation of Mg- and Ni-Containing Solutions: Implications for Carbon Mineralization and Critical Element Recovery

Carbonation of alkaline earth metals ( e.g. , magnesium (Mg)) and sulfidation of nickel (Ni) are promising methods to achieve concurrent carbon mineralization and selective Ni recovery. However, the coexistence of alkaline earth metals and Ni from silicate ores or mining wastewater complicates the carbonation and sulfidation owing to cation coprecipitation. To better understand simultaneous metal carbonation and Ni-sulfide formation, we used Mg- and Ni-containing solutions and systematically investigated the Mg and Ni coprecipitates’ phase transformation during sequential/concurrent carbonation and sulfidation. During a single carbonation process, hydromagnesite dehydrated and formed magnesite over time. Nickel bicarbonate formed and became a Mg–Ni carbonate solid solution because of their similar ionic radii. During a single sulfidation process, the pH did not affect Ni-sulfide formation, but it controlled Mg behavior. Specifically, at pH 9.6, brucite formed, while at pH 7.8, Mg 2+ remained in the solution. For the sequential carbonation–sulfidation process, Ni-carbonate formed during carbonation converted to Ni-sulfide because of the low Ni-sulfide K sp . For the sulfidation–carbonation process, Ni-sulfide remained the same even after carbonation and Mg-carbonate precipitates. For the concurrent carbonation and sulfidation process, Mg-carbonate and Ni-sulfide formed simultaneously. Finally, this study develops a scientific foundation of carbonation and sulfidation processes, benefiting coupled CO 2 storage and sulfide-enabled resource recovery.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Lattice vacancy migration barriers in Fe-Ni alloys, and an indication as to why Ni atoms diffuse slowly: A first-principles study

Lattice vacancy migration barriers in ferromagnetic Fe 𝑥 ⁢Ni 1−𝑥 alloys (0.4 ≤ 𝑥 ≤ 0.6) are accurately quantified within the framework of ab initio electronic structure calculations using the nudged elastic band (NEB) method. Both the atomically disordered (A1) fcc phase, as well as the atomically ordered, tetragonal L⁢1 0 phase—which is under consideration as a material for a rare-earth-free gap magnet for advanced engineering applications—are investigated. Across an ensemble of NEB calculations performed on supercell configurations spanning a range of compositions and containing disordered, partially ordered, and fully ordered structures, we find that Ni-vacancy interchanges encounter significantly higher energetic barriers than do Fe-vacancy interchanges. We contend that this aspect is a key factor in determining the differences in mobility between Fe and Ni atoms in this ferromagnetic alloy. Moreover, we are able to interpret these findings in terms of the ferromagnetic alloy's underlying spin-polarized electronic structure. Specifically, we report a coupling between the size of local lattice distortions and the magnitude of the local electronic spin polarization around vacancies. This causes Fe atoms to relax into lattice vacancies, while Ni atoms remain rigidly fixed to their original lattice positions. These results give atomic-scale insight into the longstanding experimental observation that Ni exhibits remarkably slow atomic diffusion in Fe-Ni alloys.

density functional theory↗
DOE OSTI2022

Ferromagnetic MnBi 4 Te 7 obtained with low-concentration Sb doping: A promising platform for exploring topological quantum states

The tuning of the magnetic phase, chemical potential, and structure is crucial to observe diverse exotic topological quantum states in Mn Bi 2 Te 4 (Bi 2 Te 3 ) m (m = 0–3). Here we show a ferromagnetic (FM) phase with a chiral crystal structure in Mn (Bi 1–x Sb x ) 4 Te 7 , obtained via tuning the growth conditions and Sb concentration. Unlike previously reported Mn (Bi 1–x Sb x ) 4 Te 7 , which exhibits FM transitions only at high Sb doping levels, our samples show FM transitions (T C = 13.5 K) at 15%–27% doping levels. Furthermore, our single-crystal x-ray-diffraction structure refinements find Sb doping leads to a chiral structure with the space group of P3, contrasted with the centrosymmetric $P\bar{3}m1$ crystal structure of the parent compound MnBi 4 Te 7 . Through angle-resolved photoemission spectroscopy measurements, we also demonstrated that the nontrivial band topology is preserved in the Sb-doped FM samples. Given that the nontrivial band topology of this system remains robust for low Sb doping levels, our success in making FM Mn(Bi 1–x Sb x ) 4 Te 7 with x = 0.15, 0.175, 0.2, and 0.27 paves the way for realizing the predicted topological quantum states, such as the axion insulator and Weyl semimetals. Additionally, we also observed magnetic glassy behavior in both antiferromagnetic MnBi 4 Te 7 and FM Mn (Bi 1–x Sb x )4 Te 7 samples, which we believe originates from cluster spin-glass phases coexisting with long-range antiferromagnetic/FM orders. Further, we have also discussed how the antisite Mn ions impact the interlayer magnetic coupling and how FM interlayer coupling is stabilized in this system.

36 MATERIALS SCIENCE↗
NASA NTRS2018

Robust Exploration and Commercial Missions to the Moon Using Nuclear Thermal Rocket Propulsion and Lunar Liquid Oxygen Derived from FeO-Rich Pyroclasitc Deposits

The nuclear thermal rocket (NTR) has frequently been identified as a key space asset required for the human exploration of Mars. This proven technology can also provide the affordable access through cislunar space necessary for commercial development and sustained human presence on the Moon. It is a demonstrated technology capable of generating both high thrust and high specific impulse (I(sub sp) approx. 900 s) twice that of today's best chemical rockets. Nuclear lunar transfer vehicles-consisting of a propulsion stage using three approx. 16.5-klb(sub f) small nuclear rocket engines (SNREs), an in-line propellant tank, plus the payload-are reusable, enabling a variety of lunar missions. These include cargo delivery and crewed lunar landing missions. Even weeklong ''tourism'' missions carrying passengers into lunar orbit for a day of sightseeing and picture taking are possible. The NTR can play an important role in the next phase of lunar exploration and development by providing a robust in-space lunar transportation system (LTS) that can allow initial outposts to evolve into settlements supported by a variety of commercial activities such as in-situ propellant production used to supply strategically located propellant depots and transportation nodes. The use of lunar liquid oxygen (LLO2) derived from iron oxide (FeO)-rich volcanic glass beads, found in numerous pyroclastic deposits on the Moon, can significantly reduce the launch mass requirements from Earth by enabling reusable, surface-based lunar landing vehicles (LLVs)that use liquid oxygen and hydrogen (LO2/LH2) chemical rocket engines. Afterwards, a LO2/LH2 propellant depot can be established in lunar equatorial orbit to supply the LTS. At this point a modified version of the conventional NTR-called the LO2-augmented NTR, or LANTR-is introduced into the LTS allowing bipropellant operation and leveraging the mission benefits of refueling with lunar-derived propellants for Earth return. The bipropellant LANTR engine utilizes the large divergent section of its nozzle as an ''afterburner'' into which oxygen is injected and supersonically combusted with nuclear preheated hydrogen emerging from the engine's choked sonic throat-essentially ''scramjet propulsion in reverse.'' By varying the oxygen-to-hydrogen mixture ratio, LANTR engines can operate over a range of thrust and I(sub sp) values while the reactor core power level remains relatively constant. A LANTR-based LTS offers unique mission capabilities including short-transit-time crewed cargo transports. Even a ''commuter'' shuttle service may be possible allowing ''one-way'' trip times to and from the Moon on the order of 36 hours or less. If only 1% of the extracted LLO2 propellant from identified resource sites were available for use in lunar orbit, such a supply could support daily commuter flights to the Moon for many thousands of years! This report outlines an evolutionary architecture and examines a variety of mission types and transfer vehicle designs, along with the increasing demands on LLO2 production as mission complexity and velocity change delta V requirements increase. A comparison of vehicle features and engine operating characteristics, for both NTR and LANTR engines, is also provided along with a discussion of the propellant production and mining requirements associated with using FeO-rich volcanic glass as source material.

rocket↗
NASA NTRSDate not supplied

Spaceflight Medical Evacuation Risk Assessment Principles - A Qualitative Investigation

BACKGROUND Future human space exploration beyond Low Earth Orbit (LEO) will require innovative solutions in many areas, primarily those that provide direct medical support to crew on long-duration missions and optimize their health and performance. The associated challenges therein will be numerous, including but not necessarily limited to extended one-way or "asynchronous" communication delays, minimal to non-existent resupply, and a prolonged transit time to "definitive care" ranging from 3-days to 9-months. At the same time, long-duration exploration spacecraft and crew will face restrictions on mass, power, volume, and data far more significant than that seen in current LEO settings. Given these limitations, medical risk assessment is of primary importance, especially evaluating the implications of a medical evacuation of an ill or injured crewmember. Such evacuations are complicated, potentially dangerous, and well may be impossible in certain phases of the mission. Regardless, such issues must be weighed against the risks of the injured crew remaining aboard a spacecraft with limited medical resources. OBJECTIVE This qualitative study drew from the experiences of subject matter experts (SMEs) in spaceflight and appropriate analog environments (i.e., military, disaster, and extreme environment fields) to identify unique principles common amongst medical evacuation considerations helpful in informing future risk assessment tools. Appropriate analog environments included austere operational settings where multiple factors (weather, logistics/limited resupply, denied/extreme environments, and patient condition) resulted in a limited ability to provide definitive local medical care. The fundamental principles in question revolved around scenarios where evacuation became a complicating yet necessary consideration and where life-threatening medical concerns had to be weighed against critical mission objective(s). The primary authors collected semi-structured data gathered through in-depth interviews with 16 subject matter experts (SMEs). Interview questions investigated how these SMEs consider and weigh the attendant risks present in medical evacuation scenarios. Among the critical questions posed were those that sought to understand how the SME balanced the challenges and requirements of medical evacuation (or keeping an injured patient "on-site" aka: "prolonged field care") against the evacuation operation's risks on impacting overarching mission success. The team analyzed interview transcripts for common themes and principles using the qualitative methods of thematic analysis based on consensus, co-occurrence, and comparison. As a result, nine primary risk consideration themes and nine contributing factor themes emerged, all of which will ideally inform future medical evacuation decision-making tools and operational decision-making for exploration class missions. Specific aims for this study included: 1. Identification of common principles used to assess risks and benefits of medical evacuations in extreme environments 2. Identification of common points of friction or complication and challenges in extreme environment evacuations

A T Almand↗
NASA NTRS2024

Development of a High Frequency Class-A Amplifier for High Temperature Wireless Microsystems based on a Silicon Carbide Static Induction Transistor

Next generation gas turbine engine health monitoring for high performance jet engine-powered aircraft will incorporate integrated microelectronic sensors specifically designed to enhance aircraft functionality, engine efficiency, and safety. Temperature, air flow and pressure associated with the gas path as well as emissions from the engine correlate to the state of the combustor and other engine health conditions. Monitoring these parameters accurately requires positioning electronic sensors in the hot zone of the gas turbine engine, which exposes them to temperatures in excess of 400ºC. At these temperatures, signal amplification is needed at the transducer level to distinguish the desired measured signal from the significant electronic noise generated by the harsh environment. In this presentation, we report on the simulation, development, fabrication, and testing on a high temperature Class-A amplifier which utilizes a 4H-SiC Static Induction Transistor (SIT) as the active device and input/output matching and DC bias networks comprised of thin-film spiral inductors, metal-insulator-metal capacitors, and thick film chip resistors. A small signal model that emulates the operation of the 4H-SiC SIT from 25 to 400ºC, with an emphasis on operation at 400οC, was utilized. Measurements were performed from 25 to 400ºC to generate current-voltage curves, capacitive transistor characteristics and high frequency scattering parameters (S-parameters). The measured data was used to extrapolate the transconductance, gm, as a function of temperature for model development. Circuit simulation tools were used to generate S-parameters, which were compared to the measured values. At 400οC, a maximum difference between measured and simulated S-parameters for frequencies from 20 to 100 MHz were 3.84%, 0.68%, 10.61%, and 3.26% for S21, S12, S11, and S22, respectively. The average transit frequency, ft, was calculated from measured values to be 197.8 MHz, while the simulated value from the model was found to be 200 MHz. The Amplifier’s S-parameters were recorded at 20-100 MHz over a temperature range of 25-400ºC and show a gain of approximately 15.8 and 5.80 dBm at 25 and 400ºC, respectively. The input and output reflection coefficients at 50 MHz and 400ºC were -18.5 and -15.2 dB, respectively. The noise figure and phase noise were measured over the temperature range of 25-400ºC and recorded. The noise figure increased 21% at 50 MHz over the temperature range, while the 1 kHz offset of the phase noise remained below -110 dB. The stability factor, K, calculated with measured and simulated data demonstrates unconditional stability over the frequency range at 400ºC. Lastly, the 1-dB compression point was measured at 50 MHz and 400ºC with an approximated output of 9.5 dBm.

SiC↗
NASA NTRSDate not supplied

Machine Learning Approaches for Rare-Earth Silicate Environmental Barrier Coating Thermochemical and Thermomechanical Property Predictions

Environmental barrier coatings (EBCs) are a necessary enabling technology for the transition from superalloys to silicon carbide (SiC) ceramic matrix composites (CMCs) in gas turbine engines for increased efficiency and decreased fuel costs. SiC-based CMCs are prone to oxidation-based degradation in the engine hot section, and rare-earth (RE) silicates are promising candidates for EBCs due to their close thermal expansion match to the composite substrate and oxidation resistance. However, the design of EBCs is hindered by the large chemical space of candidate materials and the difficulty in obtaining material properties for engineering optimization. This is especially difficult as research continues into mixed-cation or “high-entropy” RE silicates. First-principles computational methods such as density functional theory (DFT) are highly effective at calculating material properties to guide coating design but are limited by their computational cost. Atomistic simulations have the potential to both accelerate property calculations and expand the properties able to be calculated due to their lower computational compared to DFT. However, they require interatomic potentials (IAPs) specific to the material system of interest, and, to our knowledge, there are no suitable IAPs for RE silicates. Machine learning (ML) is a promising technique to accelerate material property predictions indirectly by generating IAPs for atomistic simulations or via direct prediction. In this work, we present two ML approaches to accelerate the calculation of RE silicate properties relevant to EBC design: 1) a ML-derived interatomic potential (IAP) for atomistic simulations of yttrium disilicate (Y2Si2O7) from DFT training data, and 2) a neural network (NN) model to directly predict thermochemical properties of RE silicates and oxides directly from easily obtainable unit cell parameters. Classical MD simulations using the IAP yield lattice properties and bond lengths in good agreement with both DFT and experimental results from x-ray diffraction. Thermodynamic properties calculated using the finite-displacement phonon method and quasi-harmonic approximation were orders of magnitude faster than DFT with good agreement to the DFT results. The IAP was also used to calculate properties such as coefficient of thermal expansion (CTE) that require large simulation supercells and are therefore difficult with DFT. The IAP correctly predicted the anisotropic nature of the CTE in three different phases of Y2Si2O7. The NN model predicts constant pressure heat capacity, Cp, orders of magnitude faster than DFT calculations, which can enable its use as a surrogate model for multiscale simulations. The two methods presented in this work demonstrate the utility of ML for accelerating the prediction of RE silicate properties, which can in turn accelerate EBC design and optimization.

machine learning↗
DOE OSTI2024

Machine Learning for Joint Quality Control

The use of lightweight material combinations has been highly demanded in manufacturing automotive structures. However, making robust dissimilar material joints of such lightweight materials is still challenging. A significant barrier to achieving high-quality and repeatable joint performance is a deficient understanding of the relationship between the welding process, joint attributes, and joint performance. In this context, welding factors refer to material, equipment, environment, and process parameters, while joint features comprise specific microstructural attributes of the weld such as nugget size, heat affected zone (HAZ) topology, intermetallic layer thickness, and sheet thickness reduction. Joint performance is quantified in terms of strength (e.g., tensile shear, coach peel, cross-tension), weld size, and hardness, among other factors. While there have been many attempts to establish this process-structure-property relationship by developing a model derived from the associated physics and first principles, the complexity of the joining processes compounded by the complex interactions with different materials in an automotive assembly line environment, has hindered the usefulness of such attempts. The complexity is further exacerbated using different stacking materials, especially comprising dissimilar material combinations. In practice, the common approach has been the laborious process of creating welds, characterizing them, and then physically testing them through experimentation. With the emergence of artificial intelligence (AI) methods, an alternative pathway to eliciting the desired process-structure-property relationship at an accelerated pace is to use a data-driven approach by employing machine-learning (ML) techniques. This approach is benefitted by the availability of large streams of data, generated through years of research and testing by original equipment manufacturers, in the form of material, process, environmental, equipment, microstructural, and bulk-scale performance information from multimodal, multiscale sensors making measurements from laboratory-scale to production-scale processes. During Phase I efforts, which ended in fiscal year (FY) 2021, the Oak Ridge National Laboratory and Pacific Northwest National Laboratory (ORNL/PNNL) team demonstrated the effectiveness of different ML/AI frameworks in modeling complex relationships between resistance spot welding (RSW) process parameters, weld attributes, and joint properties using a subset of data from General Motors (GM). In FY 2022, the project team further refined and expanded their respective ML models to analyze additional welds with new weld stack-ups and materials to enhance the ML model predictive capability. ORNL extended its unified deep neural networks (DNN) ML training and prediction framework with new data streams of process parameters, and PNNL extended its model describing RSW process parameters’ associations with weld attributes. In FY 2023, the project team completed the development of the AI/ML architecture for analyzing aluminum/steel joints manufactured by GM via RSW and transitioned into the inline welding quality monitoring task for steel/steel RSW joints provided by GM.

36 MATERIALS SCIENCE↗
NASA NTRS2013

Characterizing Uranus with an Ice giant Planetary Origins Probe (Ice-POP)

We now know from studies of planetary transits and microlensing that Neptune-mass planets are ubitquitous and may be the most common class of planets in the Galaxy. As such it is crucial that we understand the formation and evolution of the ice giant planets in our own solar system so that we can better understand planet formation throughout the galaxy. An entry probe mission to Uranus would help accomplish this goal. In fact the Planetary Decadal Survey recommended a Uranus orbiter with entry probe but did not explore in detail the specifications for the entry probe. NASA Ames is currently studying thermal protection system requirements for such a mission and this has led to questions regarding the minimum interesting science payload of such an entry probe. The single most important in-situ measurement for an ice giant entry probe is a measurement of atmospheric composition. For Uranus this would specifically include the methane and noble gas abundances. An in situ measurement of the methane abundance, from below the methane cloud, would constrain the atmospheric carbon abundance, which is believed to be roughly 30 to 50 times solar. There are hints from the transiting planets that extrasolar ice giants show comparable or even greater enhancements of heavy elements compared to their primary stars. However the origin of this carbon enhancement is controversial. Is Uranus a "failed core" of a larger gas giant or was the atmosphere enhanced by accretion of icy planetesimals' Constraining atmospheric abundances of C and perhaps S or even N from below 5 bars would provide badly needed data to address such issues. A measurement of the N abundance would provide clues on the origin of the planetesimals that formed Uranus. Low N-abundance indicates planetesimals from 'warmer' regions where N was mainly in form of NH3, whereas a strong enrichment could indicate planetesimals / cometary material from the colder outer regions of the nebula. Furthermore CO and HCN have been detected in Neptune but not in Uranus. A measurement of the abundance of either would constrain the source mechanisms for these molecules (exogenic or internal). A major surprise from the Galileo Entry Probe was that the heavier noble gases Ar, Kr, and Xe are enhanced in Jupiter's atmosphere at a level comparable to what was seen for the chemically active volatiles N, C, and S. It had been generally expected that Ar, Kr, and Xe would be present in solar abundances, as all were expected to accrete with hydrogen during the gravitational capture of nebular gases. Enhanced abundances of Ar, Kr, and Xe is equivalent to saying that these noble gases have been separated from hydrogen. There are several mechanisms that could accomplish this but these hypotheses require further testing. Measurement of noble gas abundances in an ice giant would constrain the planetary formation and nebular mechanisms responsible for this enhancement. Standard three-layer models of Uranus find that the outer, predominantly H/He layer of Uranus does not reach pressures high enough (approximately 1 Mbar) for H2 to transition to liquid metallic hydrogen. However, valid models can also be constructed with a smaller intermediate water-rich layer, with hydrogen then reaching the metallic hydrogen phase. If this occurs, He should phase separate from the hydrogen and ``rain out," taking along a substantial abundance of Ne, as suggested for Jupiter (and likely also for Saturn). Hence He and Ne depletions could be probes of the planet's structure in the much deeper interior. A determination of Uranus' atmospheric abundances, particularly of the noble gasses, is thus critical to understanding the formation of Uranus, and giant planets in general. These measurements can only be performed with an entry probe. The second key measurement would be a temperature-pressure sounding to provide ground truth for remote measurements of atmospheric temperature and composition and to constrain the internal heat flow. This would also establish that the methane abundance measurements have indeed been made below any possible methane cloud. Finally an ultra stable oscillator would measure wind speeds and constrain atmospheric dynamics. In our presentation we will discuss the importance of all of these measurements and argue that an entry probe is a crucial component of any ice giant mission.

Marley, Mark S.↗
DOE OSTI2026

Dynamic compression and spallation behavior of NbTaTiVZr high-entropy alloy

High-entropy alloys are a new class of materials with promising properties for aerospace and defense applications. The unique behavior of these materials is driven by the complex interactions between dissimilar atoms in a crystal and requires combined theoretical and experimental efforts to unlock their full potential. Here, we evaluate the relationship between the microstructure and the dynamic response of the equiatomic NbTaTiVZr alloy. Specifically, the shock Hugoniot, or equation of state, was measured up to a particle velocity of 0.5 mm μs −1 using gas gun plate impacts. Shock wave profile and incipient spallation experiments were used to characterize wave propagation and damage formation with post-mortem recovery experiments. Molecular dynamics simulations confirm the experimental findings and extend them up to 2.1 mm μs −1 particle velocity. Computational thermodynamic calculation of phase diagram simulations explain details of the material microstructure, which explains the measured strength and experimental damage patterns. Overall, this work provides detailed high-strain-rate characterization of a refractory high-entropy alloy, and more importantly, a framework and demonstration of the utility and necessity of a combined theoretical and experimental approach, outlining the importance of considering processing and manufacturing conditions when evaluating the performance of new materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Optimal surface-tension isotropy in the Rothman-Keller color-gradient lattice Boltzmann method for multiphase flow

The Rothman-Keller color-gradient (CG) lattice Boltzmann method is a popular method to simulate two-phase flow because of its ability to deal with fluids with large viscosity contrasts and a wide range of interfacial tensions. Here, two fluids are labeled red and blue, and the gradient in the color difference is used to compute the effect of interfacial tension. It is well known that finite-difference errors in the color-gradient calculation lead to anisotropy of interfacial tension and errors such as spurious currents. Here, we investigate the accuracy of the CG calculation for interfaces between fluids with several radii of curvature and find that the standard CG calculations lead to significant inaccuracy. Specifically, we observe significant anisotropy of the color gradient of order 7% for high curvature of an interface such as when a pinchout occurs. We derive a second order accurate color gradient and find that the diagonal nearest neighbors can be weighted differently than in the usual color-gradient calculation such that anisotropy is minimized to a fraction of a percent. The optimal weights that minimize anisotropy for the smallest radius of curvature interface are found to be w = (0.298, 0.284, 0.275) for diagonal nearest neighbors for the cases of the interface smoothing parameter β = (0.5, 0.7, 0.99), somewhat higher than the w = 0.25 value derived by Leclaire et al. [Leclaire, Reggio, and Trepanier, Computers and Fluids 48, 98 (2011)] based on obtaining isotropic errors to second order. We find that use of these optimal w values yields over a factor of 10 decrease in anisotropy and over a factor of 30 decrease in mean anisotropy relative to using the standard w = 1 value. And we find a factor of about 2 decrease in the anisotropic error and up to factor 15 decrease in mean anisotropic error relative to the choice of w = 0.25 for small radius of curvature interfaces. The improved CG calculations will allow the method to be more reliably applied to studies of phenomenology and pore scale processes such as viscous and capillary fingering, and droplet formation where surface-tension isotropy of narrow fingers and small droplets plays a crucial role in correctly capturing phenomenology. We present an example illustrating how different phenomena can be captured using the improved color-gradient method. Namely, we present simulations of a wetting fluid invading a fluid filled pipe where the viscosity ratio of fluids is unity in which droplets form at the transition to fingering using the improved CG calculations that are not captured using the standard CG calculations. We present an explanation of why this is so which relates to anisotropy of the surface tension, which inhibits the pinchouts needed to form droplets.

58 GEOSCIENCES↗
DOE OSTI2020

Ferromagnetic cluster-glass phase in Ca(Co 1-x Ir x ) 2-y As 2 crystals

In this study, single crystals of Ca(Co 1-x Ir x ) 2-y As 2 with 0≤x≤0.35 and 0.10≤y≤0.14 have been grown using the self-flux technique and characterized by single-crystal x-ray diffraction (XRD), energy-dispersive x-ray spectroscopy, magnetization M, and magnetic susceptibility χ measurements versus temperature T, magnetic field H, and time t, and heat-capacity C p (H,T) measurements. The XRD refinements reveal that all the Ir-substituted crystals crystallize in a collapsed-tetragonal structure as does the parent CaCo 2-y As 2 compound. A small 3.3% Ir substitution for Co in CaCo 1.86 As 2 drastically lowers the A-type antiferromagnetic (AFM) transition temperature T N from 52 to 23 K with a significant enhancement of the Sommerfeld electronic heat-capacity coefficient. The A-type AFM structure consists of ab-plane layers of spins ferromagnetically aligned along the c axis with AFM alignment of the spins in adjacent layers along this axis. The positive Weiss temperatures obtained from Curie-Weiss fits to the χ(T>T N ) data indicate that the dominant magnetic interactions are ferromagnetic (FM) for all x. A magnetic phase boundary is inferred to be present between x=0.14 and x=0.17 from a discontinuity in the x dependencies of the effective moment and Weiss temperature in the Curie-Weiss fits. FM fluctuations that strongly increase with increasing x are also revealed from the χ(T) data. Furthermore, the magnetic ground state for x≥0.17 is a spin glass as indicated by hysteresis in χ(T) between field-cooled and zero-field-cooled measurements and from the relaxation of M in a small field that exhibits a stretched-exponential time dependence. The spin glass has a small FM component to the ordering and is hence inferred to be comprised of small FM clusters. The competing AFM and FM interactions along with crystallographic disorder associated with Ir substitution are inferred to be responsible for the development of a FM cluster-glass phase. A logarithmic T dependence of Cp at low T for x=0.14 is consistent with the presence of significant FM quantum fluctuations. This composition is near the T=0 boundary at x≈0.16 between the A-type AFM phase containing ferromagnetically-aligned layers of spins and the FM cluster-glass phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Electrochemically-Formed Disordered Rock Salt ω-Li x V 9 Mo 6 O 40 as a Fast-Charging Li-Ion Electrode Material

Electrochemically-formed disordered rock salt compounds are an emerging class of Li-ion electrode materials for fast-charging energy storage. However, the specific factors that govern the formation process and the resulting charge storage performance are not well understood. Here, we characterize the transformation mechanism and charge storage properties of an electrochemically-formed disordered rock salt from V 9 Mo 6 O 40 (VMO). The crystal structure of VMO has similar motifs to that of α-V 2 O 5 , a well-studied analogue, but VMO has less mechanical flexibility due to additional corner-sharing octahedra in its structure. As a result, VMO undergoes a single-step transformation pathway, which we characterize through operando X-ray diffraction, and forms an unusual highly distorted lamellar microstructure, as we show with high-resolution transmission electron microscopy. The resulting Li x VMO material shows fast charging and other electrochemical characteristics and performance typical of many nanomaterials, even though the material is composed of relatively large particles.

25 ENERGY STORAGE↗
DOE OSTI2020

Drive-pressure optimization in ramp-wave compression experiments through differential evolution

Ramp-wave dynamic-compression experiments are used to examine quasi-isentropic loading paths in materials. The gradual and continuous increase in pressure created by ramp waves make these types of experiments ideal for studying nonequilibrium material behavior, such as solidification kinetics. In ramp-wave compression experiments, the input drive pressure to the experimental setup may be exerted through one of a number of different mechanisms (e.g., magnetic fields, gas-gun-driven impactors, or high-energy lasers) and is generally required for simulating such experiments. Yet, regardless of the specific mechanism, this drive pressure cannot be measured directly (measurements are generally taken at a location near the back of the experimental setup through a transparent window), leading to an inverse problem where one must determine the drive pressure at the front of the experimental setup (i.e., the input) that corresponds to the particle velocity (the output) measured near the back of the experimental setup. Furthermore, we solve this inverse problem using a heuristic optimization algorithm, known as differential evolution, coupled with a multiphysics, hydrodynamics code that simulates the compression of the experimental setup. By running many rounds of forward simulations of the experimental setup, our optimization process iteratively searches for a drive pressure that is optimized to closely reproduce the experimentally measured particle velocity near the back of the experimental setup. While our optimization methodology requires a significant number of hydrodynamics simulations to be conducted, many of these can be performed in parallel, which greatly reduces the time cost of our methodology. One novel aspect of our method for determining the drive pressure is that it does not require physical modeling of the drive mechanism and can thus be broadly applied to many types of ramp-compression experiments, regardless of the drive mechanism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Multi-objective Bayesian optimization of ferroelectric materials with interfacial control for memory and energy storage applications

Optimization of materials’ performance for specific applications often requires balancing multiple aspects of materials’ functionality. Even for the cases where a generative physical model of material behavior is known and reliable, this often requires search over multidimensional function space to identify low-dimensional manifold corresponding to the required Pareto front. In this work, we introduce the multi-objective Bayesian optimization (MOBO) workflow for the ferroelectric/antiferroelectric performance optimization for memory and energy storage applications based on the numerical solution of the Ginzburg–Landau equation with electrochemical or semiconducting boundary conditions. MOBO is a low computational cost optimization tool for expensive multi-objective functions, where we update posterior surrogate Gaussian process models from prior evaluations and then select future evaluations from maximizing an acquisition function. Using the parameters for a prototype bulk antiferroelectric (PbZrO 3 ), we first develop a physics-driven decision tree of target functions from the loop structures. We further develop a physics-driven MOBO architecture to explore multidimensional parameter space and build Pareto-frontiers by maximizing two target functions jointly—energy storage and loss. This approach allows for rapid initial materials and device parameter selection for a given application and can be further expanded toward the active experiment setting. The associated notebooks provide both the tutorial on MOBO and allow us to reproduce the reported analyses and apply them to other systems (https://github.com/arpanbiswas52/MOBO_AFI_Supplements).

36 MATERIALS SCIENCE↗
NASA NTRS2021

NASA Electric Vertical Takeoff and Landing (eVTOL) Aircraft Technology for Public Services – A White Paper

History has shown that our personal life is highly dependent on the technology that people have developed. A strategic scan of the aerospace environment at the beginning of the 21st century strongly suggests that the world might be approaching a new age of airpower—the era of electrified/hybrid aircraft propulsion. Undeniably, starting from the Montgolfier Brothers balloon flight in 1783, to the Wright Brothers piston engine flight in 1903, and the jet engine of the 1960s, or the space age of today, one can say that leaps in propulsion technology have marked the different ages of human flight. The technological advancements, brought at the beginning of 21st century by the revolution in data exchange, computational power, sensors, wireless communication, internet, and autonomy, contributed to the vision of this new age of propulsion we are approaching. Historically, conventional vertical takeoff and landing (VTOL) aircraft have been equipped with propulsion units relying on complex internal combustion machines (turbines, piston engines, for example), and complex mechanical arrangements (gearboxes, shafts, variable pitch propeller). By contrast, electric VTOL aircraft (eVTOL)1 rely on simpler propulsion units (electric motors and in some cases fixed-pitch propellers). This promotes redundancy and improves tolerance to failures, in turn improving safety. The use of simpler electric propulsion units should also allow significant acquisition and operating cost reductions. Whether full-electric (relying solely on batteries) or hybrid-electric (relying on a combination of batteries, fuel-powered engines, and generators.), eVTOLs are also expected to generate less noise and air pollution than conventional aircraft with similar payloads. According to the 2019 Annual Review of IATA (International Air Transport Association) [ref.1], due to an expected increase in air transport traffic by 5% every year and a doubling of air transport passenger numbers to 8.2 billion by 2037 significant challenges are posed to the aviation industry. Furthermore, this report does not factor in the expected demand for short-range (intra-city) air transportation, which is in development and yet to be operational. The increased demand to fly creates a responsibility to expand in a sustainable manner and an endeavor to develop more environmentally-friendly aircraft. eVTOL aircraft, either piloted or autonomous, is gathering considerable interest worldwide. Modern and novel full-electric or hybrid-electric eVTOL configurations enable a new paradigm shift in air transportation as the aviation industry remains committed to its goals of carbon-neutral growth from 2020 onwards and cutting CO2 emissions to half 2005 levels by 2050. While electric power has been used for decades, recent developments in mobile electric/hybrid propulsion coupled with advanced materials and autonomous systems may create the possibility to transition into the next age of air mobility propelled by electric/hybrid VTOL aircraft technology. Although eVTOL aircraft might seem like an incremental improvement or even a counterintuitive regression with regard to past VTOL development, it has in fact the potential to transform air mobility across a wide range of government applications. Previous transformations in aviation generated dramatic leaps in performance, but the cost was commensurate with performance, limiting quantity produced. This next age appears to take a different approach. Performance may not increase, but at this moment technology is poised for future urban mobility that will spawn commercial passenger drone services, that is, autonomous (pilotless) air taxis and thereby add a new dimension to the urban transportation mix of the future [ref. 2]. Advances in electric propulsion, autonomous flight technology, and 5G communication networks will enable this fast new-growing market to become a reality. It is now time to envision the introduction of electric/hybrid eVTOL aircraft for Public Services2. We believe that in the next decades eVTOL aircraft will have the potential to become an essential tool to Public Service agencies around the world in applications such as firefighting, public safety, search and rescue, disaster relief and law enforcement. This is due to several major factors. • First, with the increasing popularity of small, unmanned aircraft vehicles (UAVs) or drones, many companies today are focusing on the development of passenger UAVs designed to accommodate up to five passengers or equivalent cargo payload. Many such configurations are electric or hybrid-electric designs with VTOL capabilities. Several of these projects have started a flight test program and many more are expected to be in the experimental and development phase in 2020. Such revolutionary vehicles could be in commercial operations by 2030. These eVTOL systems could be ready for selected Public Services missions even sooner. • Second, although these advanced eVTOL vehicles under development still need access to fuel (hybrid) and/or electric charging capability, they can take off and land from almost anywhere. Therefore, such vehicles, both manned and unmanned can be successfully integrated for the critical missions of the Public Services with extra deployment flexibilities. • Third, advancement in electric propulsion systems in the automotive industry together with NASA’s leading efforts in electrification of aircraft propulsion systems, FAA’s ongoing active eVTOL certification programs, and EASA’s proposed framework for the certification of electric/hybrid small category VTOL aircraft in Europe [ref. 3] will help accelerate industry electric propulsion system development and integration. • Finally, eVTOL vehicles could be deployed for Public Services sooner than air taxi or other commercial applications, since Public Services missions may be more easily approved based on specific mission criteria, localized airworthiness authority for public-use aircraft3, and are normally operating under centralized airspace management and control by the theater command. Moreover, public perception and acceptance are generally less of a concern when operations save lives and benefit the wider community. The prioritized introduction of eVTOL aircraft in Public Services is ambitious, but we believe it is achievable in the coming decades if fundamental enablers (people and technologies) are engaged in defining the objectives and needs of these missions. The revolution that is currently taking place in eVTOL aircraft represents an unprecedented opportunity to develop a safer, more affordable, more available and more environmentally friendly future of vertical flight. To ensure that these novel aircraft meet the future expectations of Public Services, it is essential to take a collaborative and multi-disciplinary approach to their development, across engineering disciplines, policy-making, program management, business case development, manufacturing, and flight demonstrations. It should be noted that the term eVTOL (in the near term) used throughout this publication implies aircraft capable of transporting up to 5 persons which may or may not include a pilot if operated fully autonomously, assuming an average of 200 pounds (91 kg) per person or equivalent payload and a range up to 60 miles plus a suitable reserve. Hybrid or hydrogen powered eVTOLs would have greater range. For example, a “3-seat” eVTOL aircraft may only be able to carry two fully equipped firemen, and payload capacity is more relevant when used for the supply mission. Moreover, this paper concentrates on the “last-mile” solutions with a deployment time of no more than 6 hours. Although not specifically discussed in this document, it is understood that the future of Transformative Vertical Flight in general and Public Services, in particular, will also involve smaller UAVs that will undoubtedly play a crucial role in future aerial operations. For example, smaller unmanned aircraft may be used to dispatch medical supplies, portable filtration systems or perform the Search task of future Search and Rescue (SAR) operations. Close collaboration between the aircraft industry, the Civil Aviation Authorities (CAA), e.g., Federal Aviation Administration (FAA), European Aviation Safety Agency (EASA), Transport Canada Civil Aviation (TCCA) and the Department of Defense (DoD) certifiers, will help identify Public Services requirements, define expectations and limit development cost and timescales. Take the US Air Force Agility Prime as an example, the majority of the eVTOL application opportunities and mission elements identified are in line with the NASA TVF WG-4 objectives and use cases. Together, it forms a strong partnership to accelerate the development, certification, and practical deployment for public service missions. The US Air Force Agility Prime has been a collaboration partner on this white paper, and provided valuable input and recommendations. Most of the eVTOL public service mission elements discussed in this paper and additional use cases envisioned by the NASA TVF WG-4 team are shared by the Agility Prime program. The focus and efforts of the Agility Prime in product and system development, industry and government partnership, accelerated certifications as well as early test and deployment are totally in sync with the path forward recommended by this white paper. This kind of collaboration and partnership will help enable the practical use of the eVTOL for public service missions, benefit the eVTOL public acceptance, and accelerate the eVTOL industry revolution.

Johnny T. Doo↗