SEARCH · Engineering Papers

Results for “Transport phenomena”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,513 records · Page 84

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Characterization of the inner edge of the neutron star crust

The poorly known crustal equation of state plays a critical role in many observational phenomena associated with a neutron star. Here using semiclassical Monte Carlo simulations, we explore the possible configurations of the inner edge of the neutron star crust for a variety of baryon densities and proton fractions. Applying the Kirkwood-Buff theory to these two-component systems, we observe how the isothermal compressibility reaches a maximum when isolated nonsymmetric (or quasispherical) clusters are formed in an extremely dilute neutron gas. To determine the neutron fraction, we suggest a geometrical model based on the behavior of the proton-neutron pair correlation function. Accordingly, the equation of state of the inner crust is calculated, illustrating that the nuclear energy in β equilibrium follows a power-law behavior with baryon density. As a possible astrophysical outcome of this study, our results could help refine the mass-radius relation. Finally, our results pave the way towards further investigations of the impact of the proton-neutron pair correlation function on transport properties within the neutron star crust.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Probing Electron Correlations in 1D Electronic Materials Using Single Quantum Channels (Final Technical Report)

This research focuses on the investigation of strong electron correlations in quasi-one-dimensional oxide nanostructures. The archetypical system consisting of heterostructures formed from an ultrathin layer of LaAlO 3 and TiO 2 -terminated SrTiO 3 offers a unique laboratory for probing strong correlations due to intrinsic richness of the physical system and its ability to be modulated at extreme nanoscale dimensions. Under suitable conditions, SrTiO 3 is known to exhibit a variety of important properties associated with strong electronic correlations, especially superconductivity, magnetism. Understanding the nature of these electronic correlations at fundamental scales, and in reduced spatial dimensions, will provide the scientific basis for new energy technologies, information technologies, and materials. The research program is motivated by recent discoveries and technical advances in the PIs' lab, specifically, (1) pioneering work to create conductive nanostructures to be created in LaAlO 3 /SrTiO 3 heterointerfaces with 2 nm spatial resolution, (2) the discovery of a novel correlated electronic phases in which electrons pair without forming a superconductor, and (3) the demonstrated ability to form high-quality quantum point contacts (QPCs) and achieve ballistic electron transport over micrometer-scales in LaAlO 3 /SrTiO 3 nanowires. The program has three major components. First, we have press forward to develop a full microscopic understanding of the mechanism for electron pairing in LaAlO 3 /SrTiO 3 nanostructures. Understanding the “glue” that leads to pairing and superconductivity in a doped semiconductor will help to form a scientific basis for the correlated oxide nanoelectronics platform being developed. Second, we have used their ability to create well-defined quantum channels in nanowires allows to allow novel “scattering” experiments to be performed that probe fundamentally new quasiparticle excitations using an unprecedented combination of experimental methodologies. Experiments already underway with artificially defined superlattices (5 nm period) demonstrate the feasibility of this approach. Third, we have investigated 1D correlations in proximally coupled nanowires, which are predicted exhibit Majorana-like zero modes in the absence of superconductivity or spin-orbit coupling. These studies provide insight not only into new and emergent phenomena in the rich class of complex oxides, but also fundamental mechanisms of superconductivity and magnetism. Thematically, the research aims to combine some of the most important and interesting challenges in the physics of semiconductor nanostructures and correlated electronic materials. The increased understanding of electron correlations can pave the way to the development of new materials with outstanding or novel properties, which could be relevant for energy-related technologies such as low-power high performance computing, quantum computation, or the development of novel materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

A Comparison of GPU-Accelerated Multiphase CFD Solvers on the Polaris Supercomputer: Part 1

This report is in support of the Innovative and Novel Computational Impact on Theory and Experiment (INCITE) program sponsored by the U.S. Department of Energy (USDOE). With INCITE-level resources, one project, titled BubblyFlow, was granted computational resources for the 2025 calendar year on the Polaris supercomputer at the Argonne Leadership Computing Facility (ALCF). The project aims to conduct simulations to understand the fundamental characteristics of turbulent bubbly flow phenomena in nature. Staff at the ALCF and Argonne’s Computational Science division, along with collaborators at the City College of New York and University of Illinois at Chicago, helped a summer student to assess the accuracy and performance of two high performance computing (HPC) codes. Both codes, ImExLBM and FluTAS, are fundamentally different in their mathematical and numerical modeling. However, both may be used to solve the same physical problem. The collaboration sought to better understand the differences between both codes in terms of accuracy and efficiency. This would ultimately help the BubblyFlow project better utilize resources and establish a knowledge-base of code capabilities in future simulation campaigns. We compare ImExLBM and FluTAS, two high-performance multiphase computational fluid dynamics (CFD) solvers, in terms of physical fidelity, time-to-solution, and parallel efficiency. We validate ImExLBM (Implicit-Explicit Lattice Boltzmann Method) against a canonical benchmark and assess it’s performance relative to FluTAS (Fluid Transport Accelerated Solver), a well-established open-source CFD code.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Berry curvature memory through electrically driven stacking transitions

In two-dimensional layered quantum materials, the interlayer stacking order determines both crystalline symmetry and quantum electronic properties such as Berry curvature, topology and electron correlation. Electrical stimuli can strongly influence quasi-particle interactions and the free energy landscape, thus making it possible to access hidden stacking orders with novel quantum properties and enabling dynamic engineering of these attributes. In this paper, we demonstrate electrically driven stacking transitions and a new type of nonvolatile memory based on Berry curvature in few-layer WTe 2 . The interplay of out-of-plane electric fields and electrostatic doping controls in-plane interlayer sliding and creates multiple polar and centrosymmetric stacking orders. In-situ nonlinear Hall transport reveals such stacking rearrangements result in a layer-parity-selective Berry curvature memory in momentum space, where the sign reversal of the Berry curvature only occurs in odd layer crystals. Our findings open an avenue towards exploring coupling between topology, electron correlations, and ferroelectricity in hidden stacking orders and demonstrate a new low energy cost, electrically-controlled topological memory in the atomically thin limit.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Universality of Collective Density Fluctuations in Liquids at and away from Equilibrium – An Integrated Neutron Scattering, Theoretical, and Computational Study (Final Technical Report)

Liquids, ubiquitous on earth, are prototypical disordered condensed matter. As the mediate phase between non-structured gases and ordered solids, liquids p ossess intriguing complexities that are absent on either extreme counter phase, especially when they are driven away from equilibrium. Despite decades of studies, the physics of liquids still keeps surprising us and challenges our understanding of condensed matter, emergent phenomena, and complex systems. Among the many interesting properties, the nature of collective dynamics in liquids remains elusive because of both the presence of strong interactions and the absence of translational invariance. This project studied the universal features of collective density fluctuations (or collective modes), which include both phonon-like excitations and collective relaxations, in the Q-dependent generalized hydrodynamic regime in liquids and liquid-like matter both at and away from equilibrium. We performed systematic studies of the fast collective dynamics (including both phonon-like excitations and collective relaxations) of four classes of liquids characterized by different fragilities and interactions using synergistically integrated coherent neutron scattering experiments (Inelastic Neutron Scattering (INS), Quasi-Elastic Neutron Scattering (QENS), and Neutron Spin Echo (NSE)), a ViscoElastic Hydrodynamics theory we have been developing, andMolecular Dynamics (MD) simulations and our analysis package LiquidLib. The transformative knowledge of the collective dynamics in liquids is important to understand and ultimately control the transport of mass, energy, and charge in liquids, both at and away from thermal equilibrium, which may will likely impact a wide range of disciplines.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Two dimensional electron gas in the $\delta$-doped iridates with strong spin–orbit coupling: La$_\delta$Sr 2 IrO 4

Iridates are of significant current interest because of the strong spin–orbit coupling that leads to a variety of new phenomena. Using density-functional studies, we suggest the formation of a spin-orbital entangled two dimensional electron gas (2DEG) in the -doped iridate La$_\delta$Sr 2 IrO 4 , where a single SrO layer is replaced by a LaO layer. The extra La electron resides close to the$\delta$-doped layer, partially occupying the J eff = 1/2 upper Hubbard band and thereby making the interface metallic. The magnetic structure of the bulk is destroyed near the interface, with the Ir 0 layer closest to the interface becoming non-magnetic, while the next layer (Ir 1 ) continues to maintain the AFM structure of the bulk, but with a reduced magnetic moment. The Fermi surface consists of a hole pocket and an electron pocket, located in two different Ir layers (Ir 0 and Ir 1 ), with both carriers derived from the J eff = 1/2 upper Hubbard band. The presence of both electrons and holes at the $\delta$-doped interface suggests unusual transport properties, leading to possible device applications.

2DEG↗

DOE OSTI2022

Capturing the ground state of uranium dioxide from first principles: Crystal distortion, magnetic structure, and phonons

Uranium dioxide (UO 2 ) remains a formidable challenge for first-principles approaches due to the complex interplay among spin-orbit coupling, Mott physics, magnetic ordering, and crystal distortions. Here we use DFT+U to explore UO 2 at zero temperature, incorporating all the aforementioned phenomena. The technical challenge is to navigate the many metastable electronic states produced by DFT+U, which is accomplished using f-orbital occupation matrix control to search for the ground state. We restrict our search to the high-symmetry ferromagnetic phase, including spin-orbit coupling, which produces a previously unreported occupation matrix. This newfound occupation matrix is then used as an initialization to explore the broken symmetry phases. We find the oxygen cage distortion of the 3k antiferromagnetic state to be in excellent agreement with experiments, and both the spin-orbit coupling and the Hubbard U are critical ingredients. Further, we demonstrate that only select phonon modes have a strong dependence on the Hubbard U, whereas magnetic ordering has only a small influence overall. We perform measurements of the phonon dispersion curves using inelastic neutron scattering, and our calculations show good agreement when using reasonable values of U. The quantitative success of DFT+U warrants exploration of thermal transport and other observables within this level of theory.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Anatomy of Continuous Mars SEIS and Pressure Data from Unsupervised Learning

The seismic noise recorded by the Interior Exploration using Seismic Investigations, Geodesy, and Heat Transport (InSight) seismometer (Seismic Experiment for Interior Structure [SEIS]) has a strong daily quasi-periodicity and numerous transient microevents, associated mostly with an active Martian environment with wind bursts, pressure drops, in addition to thermally induced lander and instrument cracks. That noise is far from the Earth’s microseismic noise. Quantifying the importance of nonstochasticity and identifying these microevents is mandatory for improving continuous data quality and noise analysis techniques, including autocorrelation. Cataloging these events has so far been made with specific algorithms and operator’s visual inspection. We investigate here the continuous data with an unsupervised deep-learning approach built on a deep scattering network. This leads to the successful detection and clustering of these microevents as well as better determination of daily cycles associated with changes in the intensity and color of the background noise. We first provide a description of our approach, and then present the learned clusters followed by a study of their origin and associated physical phenomena. We show that the clustering is robust over several Martian days, showing distinct types of glitches that repeat at a rate of several tens per sol with stable time differences. We show that the clustering and detection efficiency for pressure drops and glitches is comparable to or better than manual or targeted detection techniques proposed to date, noticeably with an unsupervised approach. Finally, here we discuss the origin of other clusters found, especially glitch sequences with stable time offsets that might generate artifacts in autocorrelation analyses. We conclude with presenting the potential of unsupervised learning for long-term space mission operations, in particular, for geophysical and environmental observatories.

58 GEOSCIENCES↗

DOE OSTI2023

Theory and Simulation of Metal–Insulator–Semiconductor (MIS) Photoelectrodes

A metal–insulator–semiconductor (MIS) structure is an attractive photoelectrode-catalyst architecture for promoting photoelectrochemical reactions, such as the formation of H 2 by proton reduction. The metal catalyzes the generation of H 2 using electrons generated by photon absorption and charge separation in the semiconductor. The insulator layer between the metal and the semiconductor protects the latter element from photo-corrosion and, also, significantly impacts the photovoltage at the metal surface. Understanding how the insulator layer determines the photovoltage and what properties lead to high photovoltages is critical to the development of MIS structures for solar-to-chemical energy conversion. Herein, we present a continuum model for charge-carrier transport from the semiconductor to the metal with an emphasis on mechanisms of charge transport across the insulator. The polarization curves and photovoltages predicted by this model for a Pt/HfO 2 /p-Si MIS structure at different HfO 2 thicknesses agree well with experimentally measured data. The simulations reveal how insulator properties (i.e., thickness and band structure) affect band bending near the semiconductor/insulator interface and how tuning them can lead to operation closer to the maximally attainable photovoltage, the flat-band potential. This phenomenon is understood by considering the change in tunneling resistance with insulator properties. The model shows that the best MIS performance is attained with highly symmetric semiconductor/insulator band offsets (e.g., BeO, MgO, SiO 2 , HfO 2 , or ZrO 2 deposited on Si) and a low to moderate insulator thickness (e.g., between 0.8 and 1.5 nm). Beyond 1.5 nm, the density of filled interfacial trap sites is high and significantly limits the photovoltage and the solar-to-chemical conversion rate. These conclusions are true for photocathodes and photoanodes. This understanding provides critical insight into the phenomena enhancing and limiting photoelectrode performance and how this phenomenon is influenced by insulator properties. The study gives guidance toward the development of next-generation insulators for MIS structures that achieve high performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Impact of organic solutes on capillary phenomena in water-CO2-quartz systems

Hypothesis: The migration of supercritical CO2 (scCO2) injected into underground reservoirs as part of carbon capture and storage is influenced by organic contamination affecting mineral wettability. Molecular dynamics (MD) simulations of relevant systems that incorporate representative organic solutes allow detailed investigation of changes in fundamental interfacial and capillary properties. Experiments: We use MD simulations to explore the effects of four organic solutes (quinoline, decanoic acid, coronene, sorgoleone) on the wettability of quartz by water in the presence of scCO2. We examine the impacts of polar, alkyl, and aromatic moieties as well as fluid flow velocity at elevated temperatures and pressures. Findings: Organic molecules accumulate at the water-CO2 interface, where they distribute according to their size and functional groups. Certain organics penetrate the adsorbed water film at the quartz-CO2 interface, revealing two modes of hydrogen bonding between polar organic functional group, water, and quartz surface –OH groups. Interfacial energies and contact angles are minimally impacted by organic adsorption at the water-CO2 interface, possibly due to simultaneous CO2 desorption. Non- equilibrium MD simulations reveal that flow-induced redistribution of organic compounds modulates the radii of curvature of the advancing and receding water-CO2 interfaces. Our results indicate that organic adsorption on water surfaces is likely ubiquitous during multi-phase flow in soils and sedimen- tary rocks, with implications for the mobilization and transport of organic compounds.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Estimation of stress and stress-induced permeability change in a geological nuclear waste repository in a thermo-hydrologically coupled simulation

Geologic disposal is a promising solution for a safe permanent isolation of accumulated high-level nuclear waste from nuclear powerplants. The behavior of host rock is highly coupled thermally, hydromechanically, and chemically. Numerical simulations of such coupled phenomena for the extremely long term (>100,000 years) and large length scale (>1 km) of geologic disposal remain to be computationally challenging. In this study, a methodology has been developed to approximate stress and stress-induced permeability change in host rock using only thermo-hydrological (TH) variables. A coupled thermo-hydromechanical (THM) simulation is carried out using TOUGH-FLAC simulator to model THM behaviors of a generic nuclear waste repository, in order to evaluate the performance of the developed methodology, which is implemented in a coupled TH simulation. Results show that stress and permeability change estimated by the developed methodology in the TH-coupled simulation match those calculated in the THM-coupled simulation over the simulated timeframe of over 10,000 years. Details about errors in stress and permeability estimates accrued by the developed methodology are discussed here. The developed methodology will help incorporate stress-induced permeability change into existing TH simulators for the long-term radionuclide transport in geologic disposal.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS2015

Deep Chandra Observation and Numerical Studies of the Nearest Cluster Cold Front in the Sky

We present the results of a very deep (500 ks) Chandra observation, along with tailored numerical simulations, of the nearest, best resolved cluster cold front in the sky, which lies 90 kpc (19 arcmin) to the north-west of M87. The northern part of the front appears the sharpest, with a width smaller than 2.5 kpc (1.5 Coulomb mean free paths; at 99 per cent confidence). Everywhere along the front, the temperature discontinuity is narrower than 4-8 kpc and the metallicity gradient is narrower than 6 kpc, indicating that diffusion, conduction and mixing are suppressed across the interface. Such transport processes can be naturally suppressed by magnetic fields aligned with the cold front. Interestingly, comparison to magnetohydrodynamic simulations indicates that in order to maintain the observed sharp density and temperature discontinuities, conduction must also be suppressed along the magnetic field lines. However, the northwestern part of the cold front is observed to have a non-zero width. While other explanations are possible, the broadening is consistent with the presence of Kelvin-Helmholtz instabilities (KHI) on length-scales of a few kpc. Based on comparison with simulations, the presence of KHI would imply that the effective viscosity of the intracluster medium is suppressed by more than an order of magnitude with respect to the isotropic Spitzer-like temperature dependent viscosity. Underneath the cold front, we observe quasi-linear features that are approximately 10 per cent brighter than the surrounding gas and are separated by approximately 15 kpc from each other in projection. Comparison to tailored numerical simulations suggests that the observed phenomena may be due to the amplification of magnetic fields by gas sloshing in wide layers below the cold front, where the magnetic pressure reaches approximately 5-10 per cent of the thermal pressure, reducing the gas density between the bright features.

hydrodynamics↗

NASA NTRS2001

Experiments and Model Development for the Investigation of Sooting and Radiation Effects in Microgravity Droplet Combustion

Today, despite efforts to develop and utilize natural gas and renewable energy sources, nearly 97% of the energy used for transportation is derived from combustion of liquid fuels, principally derived from petroleum. While society continues to rely on liquid petroleum-based fuels as a major energy source in spite of their finite supply, it is of paramount importance to maximize the efficiency and minimize the environmental impact of the devices that burn these fuels. The development of improved energy conversion systems, having higher efficiencies and lower emissions, is central to meeting both local and regional air quality standards. This development requires improvements in computational design tools for applied energy conversion systems, which in turn requires more robust sub-model components for combustion chemistry, transport, energy transport (including radiation), and pollutant emissions (soot formation and burnout). The study of isolated droplet burning as a unidimensional, time dependent model diffusion flame system facilitates extensions of these mechanisms to include fuel molecular sizes and pollutants typical of conventional and alternative liquid fuels used in the transportation sector. Because of the simplified geometry, sub-model components from the most detailed to those reduced to sizes compatible for use in multi-dimensional, time dependent applied models can be developed, compared and validated against experimental diffusion flame processes, and tested against one another. Based on observations in microgravity experiments on droplet combustion, it appears that the formation and lingering presence of soot within the fuel-rich region of isolated droplets can modify the burning rate, flame structure and extinction, soot aerosol properties, and the effective thermophysical properties. These observations led to the belief that perhaps one of the most important outstanding contributions of microgravity droplet combustion is the observation that in the absence of asymmetrical forced and natural convection, a soot shell is formed between the droplet surface and the flame, exerting an influence on the droplet combustion response far greater than previously recognized. The effects of soot on droplet burning parameters, including burning rate, soot shell dynamics, flame structure, and extinction phenomena provide significant testing parameters for studying the structure and coupling of soot models with other sub-model components.

Choi, Mun Young↗

NASA NTRS2001

Flight Test Analysis of the Forces and Moments Imparted on a B737-100 Airplane During Wake Vortex Encounters

Aircraft travel has become a major form of transportation. Several of our major airports are operating near their capacity limit, increasing congestion and delays for travelers. As a result, the National Aeronautics and Space Administration (NASA) has been working in conjunction with the Federal Aviation Administration (FAA), airline operators, and the airline industry to increase airport capacity without sacrificing public safety. One solution to the problem is to increase the number of airports and build new. runways; yet, this solution is becoming increasingly difficult due to limited space. A better solution is to increase the production per runway. This solution increases the possibility that one aircraft will encounter the trailing wake of another aircraft. Hazardous wake vortex encounters occur when an aircraft encounters the wake produced by a heavier aircraft. This heavy-load aircraft produces high-intensity wake turbulence that redistributes the aerodynamic loads of trailing smaller aircraft. This situation is particularly hazardous for smaller aircraft during takeoffs and landings. In order to gain a better understanding of the wake-vortex/aircraft encounter phenomena, NASA Langley Research Center conducted a series of flight tests from 1995 through 1997. These tests were designed to gather data for the development a wake encounter and wake-measurement data set with the accompanying atmospheric state information. This data set is being compiled into a database that can be used by wake vortex researchers to compare with experimental and computational results. The purpose of this research is to derive and implement a procedure for calculating the wake-vortex/aircraft interaction portion of that database by using the data recorded during those flight tests. There were three objectives to this research. Initially, the wake-induced forces and moments from each flight were analyzed based on varying flap deflection angles. The flap setting alternated between 15 and 30 degrees while the separation distance remained constant. This examination was performed to determine if increases in flap deflection would increase or decrease the effects of the wake-induced forces and moments. Next, the wake-induced forces and moments from each flight were analyzed based on separation distances of 1-3 nautical miles. In this comparison, flap deflection was held constant at 30 degrees. The purpose of this study was to determine if increased separation distances reduced the effects of the wake vortex on the aircraft. The last objective compared the wake-induced forces and moments of each flight as it executed a series of maneuvers through the wake-vortex. This analysis was conducted to examine the impact of the wake on the B737 as it traversed the wake horizontally and vertically. Results from the first analysis indicated that there was no difference in wake effect at flap deflections of 15 and 30 degrees. This conclusion is evidenced in the cases of the wake-induced sideforce, rolling moment, and yawing moment. The wake-induced lift, drag, and pitching moment cases yielded less conclusive results. The second analysis compared the wake-induced forces and moments at separation distances of 1-3 nautical miles. Results indicated that there was no significant difference in the wake-induced lift, drag, sideforce, or yawing moment coefficients. The analysis compared the wake-induced forces and moments based on different flight maneuvers. It was found that the wake-induced forces and moments had the greatest impact on out-to-in and in-to-out maneuvers.

Roberts, Chistopher L.↗

NASA NTRS2002

Solar Supergranulation Revealed as a Superposition of Traveling Waves

40 years ago two new solar phenomena were described: supergranulation and the five-minute solar oscillations. While the oscillations have since been explained and exploited to determine the properties of the solar interior, the supergranulation has remained unexplained. The supergranules, appearing as convective-like cellular patterns of horizontal outward flow with a characteristic diameter of 30 Mm and an apparent lifetime of 1 day, have puzzling properties, including their apparent superrotation and the minute temperature variations over the cells. Using a 60-day sequence of data from the MDI (Michelson-Doppler Imager) instrument onboard the SOHO (Solar and Heliospheric Observatory) spacecraft, we show that the supergranulation pattern is formed by a superposition of traveling waves with periods of 5-10 days. The wave power is anisotropic with excess power in the direction of rotation and toward the equator, leading to spurious rotation rates and north-south flows as derived from correlation analyses. These newly discovered waves could play an important role in maintaining differential rotation in the upper convection zone by transporting angular momentum towards the equator.

Gizon, L.↗

DOE OSTI2023

Simulations with machine learning potentials identify the ion conduction mechanism mediating non-Arrhenius behavior in LGPS

Abstract Li 10 Ge(PS 6 ) 2 (LGPS) is a highly concentrated solid electrolyte, in which Coulombic repulsion between neighboring cations is hypothesized as the underlying reason for concerted ion hopping, a mechanism common among superionic conductors such as Li 7 La 3 Zr 2 O 12 (LLZO) and Li 1.3 Al 0.3 Ti 1.7 (PO 4 ) 3 (LATP). While first principles simulations using molecular dynamics (MD) provide insight into the Li + transport mechanism, historically, there has been a gap in the temperature ranges studied in simulations and experiments. Here, we used a neural network potential trained on density functional theory (DFT) simulations, to run up to 40-nanosecond long MD simulations at DFT-like accuracy to characterize the ion conduction mechanisms across a range of temperatures that includes previous simulations and experimental studies. We have confirmed a Li + sublattice phase transition in LGPS around 400 K, below which the ab -plane diffusivity D a b ∗ is drastically reduced. Concomitant with the sublattice phase transition near 400 K, there is less cation-cation (cross) correlation, as characterized by Haven ratios closer to 1, and the vibrations in the system are more harmonic at lower temperature. Intuitively, at high temperature, the collection of vibrational modes may be sufficient to drive concerted ion hops. However, near room temperature, the vibrational modes available may be insufficient to overcome electrostatic repulsion, thus resulting in less correlated ion motion and comparatively slower ion conduction. Such phenomena of a sublattice phase transition, below which concerted hopping plays a less significant role, may be extended to other highly concentrated solid electrolytes such as LLZO and LATP.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

Initial Neutronics and Thermal-Hydraulic Coupling for Spent Nuclear Fuel Canister

This report documents work performed supporting the US Department of Energy (DOE) Nuclear Energy Spent Fuel and Waste Disposition, Spent Fuel and Waste Science and Technology, under work breakdown structure element 1.08.01.03.05, “Direct Disposal of Dual Purpose Canisters.” In particular, this appendix fulfills the M3 milestone, M3SF-19OR010305015, “Multiphysics criticality consequence analysis capability development status report,” as Revision 1 to M3SF-19OR010305016, “Initial neutronic and thermal hydraulic coupling for waste package” within work package SF-19OR01030501, “Direct Disposal of Dual Purpose Canisters–ORNL.” This report presents the initial development status of a multiphysics criticality consequence simulation framework, Terrenus. Terrenus currently couples two physics codes: the Monte Carlo radiation transport code Shift and the thermal-hydraulics code COBRA-SFS. The coupling capability has been demonstrated using a simplified 3 x 3 fuel pin lattice. Terrenus will be developed further to include (1) general geometry package support, allowing modeling of a spent nuclear fuel (SNF) cask with all structural details, (2) a depletion solver to determine the change in nuclide composition at the end of each time step, and (3) a mechanics code to determine any structural impact due to a criticality event. Future research and development works include identifying/developing a modern two-phase thermal-hydraulics code and developing an approach to solve fast neutronic transients. The goal of this multiphysics framework is to determine the feasibility of direct disposal of currently loaded SNF canisters by including or excluding a criticality event from a repository performance analysis framework in terms of consequences. The revision 1 to this report includes extending Terrenus to a 17 x 17 pressurized water reactor (PWR) fuel assembly. We also investigated the phenomena necessary to quench a critical reaction, and investigated the water level necessary to achieve criticality in a typical canister. Finally, we describe the initial implementation of a critical power search capability.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

Transcritical diffuse-interface hydrodynamics of propellants in high-pressure combustors of chemical propulsion systems

Rocket engines and high-power new generations of gas-turbine jet engines and diesel engines oftentimes involve the injection of one or more reactants at subcritical temperatures into combustor environments at high pressures, and more particularly at pressures higher than those corresponding to the critical points of the individual components of the mixture, which typically range from 13 to 50 bars for most propellants. This class of trajectories in the thermodynamic space has been traditionally referred to as transcritical. However, the fundamental understanding of fuel atomization, vaporization, mixing, and combustion processes at such high pressures remains elusive. In particular, whereas fuel sprays are relatively well characterized at normal pressures, analyses of dispersion of fuel in high-pressure combustors are hindered by the limited experimental diagnostics and theoretical formulations available. The description of the thermodynamics of hydrocarbon-fueled mixtures employed in chemical propulsion systems is complex and involves mixing-induced phenomena, including an elevation of the critical point whereby the coexistence region of the mixture extends up to pressures much larger than the critical pressures of the individual components. As a result, interfaces subject to surface-tension forces may persist in multicomponent systems despite the high pressures, and may give rise to unexpected spray-like atomization dynamics that are otherwise absent in monocomponent systems above their critical point. In this article, the current understanding of this phenomenon is reviewed within the context of propulsion systems fueled by heavy hydrocarbons. Emphasis is made on analytical descriptions at mesoscopic scales of interest for computational fluid dynamics. In particular, a set of modifications of the constitutive laws in the Navier–Stokes equations for multicomponent flows, supplemented with a high-pressure equation of state and appropriate redefinitions of the thermodynamic potentials, are introduced in this work based on an extended version of the diffuse-interface theory of van der Waals. The resulting formulation involves revisited forms of the stress tensor and transport fluxes of heat and species, and enables a description of the mesoscopic volumetric effects induced by transcritical interfaces consistently with the thermodynamic phase diagram of the mixture at high pressures. Additionally, applications of the theory are illustrated in canonical problems, including dodecane/nitrogen transcritical interfaces in non-isothermal systems. The results indicate that a transcritical interface is formed between the propellant streams that persists downstream of the injection orifice over distances of the same order as the characteristic thermal-entrance length of the fuel stream. The transcritical interface vanishes at an edge that gives rise to a fully supercritical mixing layer.

02 PETROLEUM↗