Search NASA⌕ Search

SEARCH · Search NASA

Results for “BOUNDARY LAYER - EFFECT OF COMPRESSIBILITY”

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.

36 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

HarDWR - Water Management Area (WMA) Shapefiles

A dataset within the Harmonized Database of Western U.S. Water Rights (HarDWR). For a detailed description of the database, please see the meta-record v2.0. Changelog v2.0 - No changes v1.0 - Initial public release Description Borders of all Water Management Areas (WMAs) across the 11 western-most states of the coterminous United States are available filtered through a single source. The legal name for this set of boundaries varies state-by-state. The data is provided as two compressed shapefiles. One, stateWMAs, contains data for all 11 states. For 10 of those states, Arizona being the exception, the polygons represent the legal management boundaries used by those states to manage their surface and groundwater resources respectively. WMAs refer to the set of boundaries a particular state uses to manage its water resources. Each set of boundaries was collected from the states individually, and then merged into one spatial layer. The merging process included renaming some columns to enable merging with all other source layers, as well as removing columns deemed not required for followup analysis. The retained columns for each boundary are: basinNum - the state provided unique numerical ID; basinName - the state provided English name of the area, where applicable; state - the state name; and uniID - a unique identifier we created by concatenating the state name, and underscore, and the state numerical ID. Arizona is unique within this collection of states in that surface and groundwater resources are managed using two separate sets of boundaries. During our followup analysis (Grogan et al., in review) we decided to focus on one set of boundaries, those for surface water. This is due to the recommendation of our hydrologists that the surface water boundary set is a more realist representation of how water moves across the landscape, as a few of the groundwater boundaries are based on political and/or economic considerations. Therefore, the Arizona surface WMAs are included within stateWMAs. The Arizona groundwater WMAs are provided as a separate file, azGroundWMAs, as a companion to the first file for completeness and general reference. WMA spatial boundary data sources by state: Arizona: Arizona Surface Water Watersheds; Collected February, 2020; https://gisdata2016-11-18t150447874z-azwater.opendata.arcgis.com/datasets/surface-watershed/explore?location=34.158174%2C-111.970823%2C7.50 Arizona: Arizona Ground Water Basins; Collected February, 2020; https://gisdata2016-11-18t150447874z-azwater.opendata.arcgis.com/datasets/groundwater-basin-2/explore?location=34.158174%2C-111.970823%2C7.50 California: California CalWater 2.2.1; Collected February, 2020; https://www.mlml.calstate.edu/mpsl-mlml/data-center/data-entry-tools/data-tools/gis-shapefile-layers/ Colorado: Colorado Water District Boundaries; Collected February, 2020; https://www.colorado.gov/pacific/cdss/gis-data-category Idaho: Idaho Department of Water Resources (IDWR) Administrative Basins; Collected November, 2015; https://data-idwr.opendata.arcgis.com/datasets/fb0df7d688a04074bad92ca8ef74cc26_4/explore?location=45.018686%2C-113.862284%2C6.93 Montana: Collected June, 2019; Directly contacted Montana Department of Natural Resources and Conservation (DNRC) Office of Information Technology (OIT) Nevada: Nevada State Engineer Admin Basin Boundaries; Collected April, 2020 https://ndwr.maps.arcgis.com/apps/mapviewer/index.html?layers=1364d0c3a0284fa1bcd90f952b2b9f1c New Mexico: New Mexico Office of the State Engineer (OSE) Declared Groundwater Basins; Collected April, 2020 https://geospatialdata-ose.opendata.arcgis.com/datasets/ose-declared-groundwater-basins/explore?location=34.179783%2C-105.996542%2C7.51 Oregon: Oregon Water Resources Department (OWRD) Administrative Basins; Collected February, 2020; https://www.oregon.gov/OWRD/access_Data/Pages/Data.aspx Utah: Utah Adjudication Books; Collected April, 2020; https://opendata.gis.utah.gov/datasets/utahDNR::utah-adjudication-books/explore?location=39.497165%2C-111.587782%2C-1.00 Washington: Washington Water Resource Inventory Areas (WRIA); Collected June, 2017; https://ecology.wa.gov/Research-Data/Data-resources/Geographic-Information-Systems-GIS/Data Wyoming: Wyoming State Engineer's Office Board of Control Water Districts; Collected June, 2019; Directly contacted Wyoming State Engineer's Office

Economics↗
DOE OSTI2026

Scaling kinetic Monte-Carlo simulations of grain growth with combined convolutional and graph neural networks

Graph neural networks (GNN) have emerged as a promising machine learning method for microstructure simulations such as grain growth. However, accurate modeling of realistic grain boundary networks requires large simulation cells, which GNN has difficulty scaling up to. To alleviate the computational costs and memory footprint of GNN, we suggest a hybrid architecture combining a convolutional neural network (CNN) based bijective autoencoder to compress the spatial dimensions, and a GNN that evolves the microstructure in the latent space of reduced spatial sizes. Our results demonstrate that the new design significantly reduces computational costs with using fewer message passing layer (from 12 down to 3) compared with GNN alone. The reduction in computational cost becomes more pronounced as the spatial size increases, indicating strong computational scalability. For the largest mesh evaluated (160 3 ), our method reduces memory usage and runtime in inference by 117× and 115×, respectively, compared with GNN-only baseline. More importantly, it shows higher accuracy and stronger spatiotemporal capability than the GNN-only baseline, especially in long-term testing. Such combination of scalability and accuracy is essential for simulating realistic material microstructures over extended time scales. The improvements can be attributed to the bijective autoencoder’s ability to compress information losslessly from spatial domain into a high dimensional feature space, thereby producing more expressive latent features for the GNN to learn from, while also contributing its own spatiotemporal modeling capability. Training data are generated from stochastic grain growth simulations, providing realistic variability for learning robust microstructure evolution. Comprehensive system validation confirms that the model is accurate, robust, and scalable.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The Flux-Differencing Discontinuous Galerkin Method Applied to an Idealized Fully Compressible Nonhydrostatic Dry Atmosphere

Dynamical cores used to study the circulation of the atmosphere employ various numerical methods ranging from finite-volume, spectral element, global spectral, and hybrid methods. In this work, we explore the use of Flux-Differencing Discontinuous Galerkin (FDDG) methods to simulate a fully compressible dry atmosphere at various resolutions. We show that the method offers a judicious compromise between high-order accuracy and stability for large-eddy simulations and simulations of the atmospheric general circulation. In particular, filters, divergence damping, diffusion, hyperdiffusion, or sponge-layers are not required to ensure stability; only the numerical dissipation naturally afforded by FDDG is necessary. We apply the method to the simulation of dry convection in an atmospheric boundary layer and in a global atmospheric dynamical core in the standard benchmark of Held and Suarez

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Data-driven closure modeling for hypersonic turbulent flows

The Reynolds-averaged Navier–Stokes (RANS) equations remain a workhorse technology for simulating compressible fluid flows of practical interest. Due to model-form errors, however, RANS models can yield erroneous predictions that preclude their use on mission-critical problems. This report summarizes work performed from FY22-FY24 focused on improving RANS models for hypersonic flows using data-driven modeling and scientific machine learning. In this work we: 1. Investigate the current capabilities of RANS models in Sandia’s parallel aerodynamics and re-entry code (SPARC) for hypersonic flows with a focus on shock boundary layer interactions (SBLIs), 2. Assess several established corrections that exist in the literature aimed at improving predictions for SBLIs, 3. Develop improved models for the Reynolds stress tensor using tensor-basis neural networks, 4. Develop a neural-network-based variable turbulent Prandtl number model to reduce errors in wall heating in SBLIs. 5. Begin future investigations including employing the LIFE framework to improve wall heating predictions in SBLIs as well as the ensemble Kalman filter. We find that current RANS models in SPARC are deficient for complex SBLI flows. In particular, no current model jointly predicts wall heat flux, wall shear stress, and wall pressure with reasonable accuracy. Existing corrections help, but do not alleviate this issue altogether. The development of improved models for the Reynolds stress tensor via tensor-basis neural networks results in more predictive RANS models across a suite of low-speed and high-speed cases. For hypersonic boundary layers, the inclusion of the wall-normal Reynolds stress via TBNNs has an appreciable impact on the wall-normal momentum balance and wall quantities. However, we find that improvements to the Reynolds stress tensor do not address the over-prediction in wall heat flux in SBLIs. We find that a neural-network-based variable turbulent Prandtl number model systematically and substantially improves wall heating predictions for a range of SBLI cases.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Low Mach number fluctuating hydrodynamics model for ionic liquids

We present a new mesoscale model for ionic liquids based on a low Mach number fluctuating hydrodynamics formulation for multicomponent charged species. The low Mach number approach eliminates sound waves from the fully compressible equations leading to a computationally efficient incompressible formulation. The model uses a Gibbs free energy functional that includes enthalpy of mixing, interfacial energy, and electrostatic contributions. These lead to a new fourth-order term in the mass equations and a reversible stress in the momentum equations. We calibrate our model using parameters for [DMPI+][F6P-], an extensively-studied room temperature ionic liquid (RTIL), and numerically demonstrate the formation of mesoscopic structuring at equilibrium in two and three dimensions. In simulations with electrode boundaries the measured double layer capacitance decreases with voltage, in agreement with theoretical predictions and experimental measurements for RTILs. Finally, we present a shear electroosmosis example to demonstrate that the methodology can be used to model electrokinetic flows.

42 ENGINEERING↗
DOE OSTI2025

Bulk-scale stress–strain hysteresis in layered crystalline solids: A study of graphite and Ti 3 SiC 2

Polycrystalline graphite and the MAX phase Ti 3 SiC 2 are layered crystalline solids with similar deformation mechanisms, including basal slip, ripplocation boundaries (RBs), kink boundaries (KBs), and cracking. The interplay of these mechanisms, notably in energy dissipation, has been much discussed in the past twenty-five years. This study builds upon previous work, investigating deformation with a renewed emphasis on the bulk-scale and given recent findings concerning RBs. Our investigation compares the evolution of energy dissipation, nonlinear recoverable and irrecoverable strain, and damage upon increasing stress for graphite and Ti 3 SiC 2 . Benitez et al.’s (2016) methodology of compressive cyclic loading and post-mortem electron backscatter diffraction (EBSD) to assess the prevalence of kinking based on low-angle grain boundaries (LAGBs) was used. Strains were measured with digital image correlation and EBSD was conducted on Ti 3 SiC 2 leveraging dictionary indexing, which was necessary herein to identify LAGBs accurately. The stress–strain stages of Ti 3 SiC 2 agree with literature on Ti 2 AlC. Damage and energy dissipation were more accelerated in graphite. No significant difference was observed in the fraction of LAGBs between pristine and unloaded Ti 3 SiC 2 . Trends observed and EBSD evidence that KBs were not dominant suggest that RBs are the primary dissipator of energy in both materials.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Detonation propagation for shock-driven, subsonic and supersonic confiner flow

In this work, we study the compressible flow dynamics of two-dimensional, steady detonation wave propagation in a high explosive (HE) confined by aluminium (Al) or stainless steel (SS), outside of which is an air layer. We examine how the thickness of the confinement affects the subsonic detonation driving zone structure (DDZ) and the detonation speed (D 0 ), demonstrating a strong dependence on whether the oblique shock-driven flow in the confiner is supersonic, as for SS, or subsonic, as for Al. A characteristic path analysis is used to examine the information flow from the material boundaries through the supersonic flow regions in both the HE and confiner that can impact the sonic surfaces bounding the subsonic flow regions. It is shown that the nature of gas-dynamic wave reflection off the SS–air or Al–air boundary can significantly influence the DDZ and D 0 .

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Atomic-Scale Scanning of Domain Network in the Ferroelectric HfO 2 Thin Film

Ferroelectric HfO 2 -based thin films have attracted much interest in the utilization of ferroelectricity at the nanoscale for next-generation electronic devices. However, the structural origin and stabilization mechanism of the ferroelectric phase are not understood because the film is typically nanocrystalline with active yet stochastic ferroelectric domains. Here, in this study, electron microscopy is used to map the in-plane domain network structures of epitaxially grown ferroelectric Y:HfO 2 films in atomic resolution. The ferroelectricity is confirmed in free-standing Y:HfO 2 films, allowing for investigating the structural origin for their ferroelectricity by 4D-STEM, high-resolution STEM, and iDPC-STEM. At the grain boundaries of <111>-oriented Pca2 1 orthorhombic grains, a high-symmetry mixed-(R3m, Pnm2 1 ) phase is induced, exhibiting enhanced polarization due to in-plane compressive strain. Nanoscale Pca2 1 orthorhombic grains and their grain boundaries with mixed-(R3m, Pnm2 1 ) phases of higher symmetry cooperatively determine the ferroelectricity of the Y:HfO 2 film. It is also found that such ferroelectric domain networks emerge when the film thickness is beyond a finite value. Furthermore, in-plane mapping of oxygen positions overlaid on ferroelectric domains discloses that polarization is suppressed at vertical domain walls, while it is active when domains are aligned horizontally with subangstrom domain walls. In addition, randomly distributed 180° charged domain walls are confined by spacer layers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Thermal Transport and Mechanical Stress Mapping of a Compression Bonded GaN/Diamond Interface for Vertical Power Devices

Bonding diamond to the back side of gallium nitride (GaN) electronics has been shown to improve thermal management in lateral devices; however, engineering challenges remain with the bonding process and characterizing the bond quality for vertical device architectures. Here, in this study, integration of these two materials is achieved by room-temperature compression bonding centimeter-scale GaN and a diamond die via an intermetallic bonding layer of Ti/Au. Recent attempts at GaN/diamond bonding have utilized a modified surface activation bonding (SAB) method, which requires Ar fast atom bombardment immediately followed by bonding within the same tool under ultrahigh vacuum (UHV) conditions. The method presented here does not require a dedicated SAB tool yet still achieves bonding via a room-temperature metal–metal compression process. Imaging of the buried interface and the total bonding area is achieved via transmission electron microscopy (TEM) and confocal acoustic scanning microscopy (C-SAM), respectively. The thermal transport quality of the bond is extracted from spatially resolved frequency-domain thermoreflectance (FDTR) with the bonded areas boasting a thermal boundary conductance of >100 MW/m 2 ·K. Additionally, Raman maps of GaN near the GaN–diamond interface reveal a low level of compressive stress, <80 MPa, in well-bonded regions. FDTR and Raman were coutilized to map these buried interfaces and revealed some poor thermally bonded areas bordered by high-stress regions, highlighting the importance of spatial sampling for a complete picture of bond quality. Overall, this work demonstrates a novel method for thermal management in vertical GaN devices that maintains low intrinsic stresses while boasting high thermal boundary conductances.

Raman↗
DOE OSTI2025

High temporal frequency data from a four turbine, blade-resolved wind farm simulation with ExaWind

The data was generated with ExaWind (https://github.com/Exawind) which couples AMR-Wind (https://github.com/Exawind/amr-wind/), Nalu-Wind (https://github.com/Exawind/nalu-wind), TIOGA (https://github.com/Exawind/tioga), and OpenFAST (https://github.com/OpenFAST/openfast). This is a large-scale simulation of a blade-resolved wind farm using the ExaWind software stack. ExaWind couples together a background flow solver, AMR-Wind, and a near-body solver, Nalu-Wind, through an overset technique from the TIOGA application. Another application, OpenFAST, handles the structural dynamics of the turbine blades and towers, which informs the fluid-structure interaction of the wind turbines with the flow solvers. This particular simulation includes four blade-resolved wind turbines operating in a turbulent atmospheric boundary layer. The AMR-Wind solver uses 500 million cells and is being solved on 256 AMD GPUs of the Oakridge Leadership Computing Facility Frontier supercomputer. Each turbine is assigned its own Nalu-Wind solver with over 13 million elements per turbine and solved using 448 CPU cores, for a total of 1792 CPU cores. For each node, 56 cores contain Nalu-Wind, while 8 cores correspond to AMR-Wind operations on the GPUs. Consequently, ExaWind is entirely utilizing the CPUs and the GPUs of the nodes concurrently. The data used in the visualization is full flow field data output from the simulation. It is lossy-compressed to a specific accuracy using ZFP and written to disk every 16 time-steps to enable real-time flow visualization. The flow fields are sampled at a high temporal frequency to enable real-time, 24fps visualization. The flow fields are sampled every 12 simulation time steps (every 0.04132s).

17 WIND ENERGY↗
DOE OSTI2022

High-temperature compressive creep tests of U 3 Si 2 with spark plasma sintering: Experiments and finite element modeling

This paper reports a systematic high-temperature creep behavior of dense U 3 Si 2 pellets using a spark plasma sintering (SPS) apparatus at elevated temperatures and pressures under vacuum conditions. The stress exponent was subsequently derived to be 3.21 at 1173 K and 2.17 at 1223 K, respectively, indicating a grain boundary sliding creep mechanism. The creep activation energy was determined to be 203.6 ± 19.0 kJ/mol, which agrees well with the literature. Finite element modeling was performed using the creep parameters fitted from the strain-time plot. Additionally, the results suggest an excellent match with the experimental data, confirming the validity of the experiments. Microstructure characterizations indicate that the main phase of the specimens after creep tests remains to be U 3 Si 2 , with a 4 μm thick layer of nano-sized particles induced from the diffusion between U 3 Si 2 and alumina disc used to avoid electric current passing through the sample. The successful conduct of creep experiments demonstrates the great potential of SPS to perform high-temperature mechanical testing of nuclear fuels under vacuum conditions. The subsequent finite element modeling exhibits excellent capabilities for accurately predicting material performance in the creep tests and provides a practical tool in evaluating nuclear fuels’ performance for a much-extended time scale.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The shape of convection in 2D and 3D global simulations of stellar interiors

Context. Theoretical descriptions of convective overshooting in stellar interiors often rely on a basic one-dimensional parameterization of the flow called the filling factor for convection. Several different definitions of the filling factor have been developed for this purpose, based on: (1) the percentage of the volume, (2) the mass flux, and (3) the convective flux that moves through the boundary. Aims. We examine these definitions of the filling factor with the goal of establishing their ability to explain differences between 2D and 3D global simulations of stellar interiors that include fully compressible hydrodynamics and realistic microphysics for stars. Methods. We study convection and overshooting in pairs of identical two-dimensional (2D) and three-dimensional (3D) global simulations of stars produced with MUSIC , a fully compressible, time-implicit hydrodynamics code. We examine pairs of simulations for (1) a 3 M ⊙ red giant star near the first dredge-up point, (2) a 1 M ⊙ pre-main-sequence star with a large convection zone, (3) the current sun, and (4) a 20 M ⊙ main-sequence star with a large convective core. Results. Our calculations of the filling factor based on the volume percentage and the mass flux indicate asymmetrical convection near the surface for each star with an outer convection zone. However, near the convective boundary, convective flows achieve inward-outward symmetry for each star that we study; for 2D and 3D simulations, these filling factors are indistinguishable. A filling factor based on the convective flux is contaminated by boundary-layer-like flows, making a theoretical interpretation difficult. We present two possible new alternatives to these frequently used definitions of a filling factor, which instead compare flows at two different radial points. The first alternative is the penetration parameter of Anders et al. (2022, ApJ, 926, 169). The second alternative is a new statistic that we call the plume interaction parameter. We demonstrate that both of these parameters captures systematic differences between 2D and 3D simulations around the convective boundary.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Perovskite superlattices with efficient carrier dynamics

Compared with their three-dimensional (3D) counterparts, low-dimensional metal halide perovskites (2D and quasi-2D; B 2 A n-1 M n X 3n+1 , such as B = R-NH 3 + , A = HC(NH 2 ) 2 + , Cs + ; M = Pb 2+ , Sn 2+ ; X = Cl - , Br - , I - ) with periodic inorganic-organic structures have shown promising stability and hysteresis-free electrical performance. However, their unique multiple-quantum-well structure limits the device efficiencies because of the grain boundaries and randomly oriented quantum wells in polycrystals. In single crystals, the carrier transport through the thickness direction is hindered by the layered insulating organic spacers. Furthermore, the strong quantum confinement from the organic spacers limits the generation and transport of free carriers. Also, lead-free metal halide perovskites have been developed but their device performance is limited by their low crystallinity and structural instability 11 . Here, in this paper, we report a low-dimensional metal halide perovskite BA 2 MA n-1 Sn n I 3n+1 (BA, butylammonium; MA, methylammonium; n = 1, 3, 5) superlattice by chemical epitaxy. The inorganic slabs are aligned vertical to the substrate and interconnected in a criss-cross 2D network parallel to the substrate, leading to efficient carrier transport in three dimensions. A lattice-mismatched substrate compresses the organic spacers, which weakens the quantum confinement. The performance of a superlattice solar cell has been certified under the quasi-steady state, showing a stable 12.36% photoelectric conversion efficiency. Moreover, an intraband exciton relaxation process may have yielded an unusually high open-circuit voltage (V OC ).

36 MATERIALS SCIENCE↗
DOE OSTI2023

A study of convective core overshooting as a function of stellar mass based on two-dimensional hydrodynamical simulations

ABSTRACT We perform two-dimensional (2D) numerical simulations of core convection for zero-age main-sequence stars covering a mass range from 3 to 20 M⊙. The simulations are performed with the fully compressible time-implicit code music. We study the efficiency of overshooting, which describes the ballistic process of convective flows crossing a convective boundary, as a function of stellar mass and luminosity. We also study the impact of artificially increasing the stellar luminosity for 3 M⊙ models. The simulations cover hundreds to thousands of convective turnover time-scales. Applying the framework of extreme plume events previously developed for convective envelopes, we derive overshooting lengths as a function of stellar masses. We find that the overshooting distance (dov) scales with the stellar luminosity (L) and the convective core radius (rconv). We derive a scaling law $d_{\rm ov} \propto L^{1/3} r_{\rm conv}^{1/2}$, which is implemented in a one-dimensional stellar evolution code and the resulting stellar models are compared to observations. The scaling predicts values for the overshooting distance that significantly increase with stellar mass, in qualitative agreement with observations. Quantitatively, however, the predicted values are underestimated for masses ≳10 M⊙. Our 2D simulations show the formation of a nearly adiabatic layer just above the Schwarzschild boundary of the convective core, as exhibited in recent three-dimensional simulations of convection. The most luminous models show a growth in size with time of the nearly adiabatic layer. This growth seems to slow down as the upper edge of the nearly adiabatic layer gets closer to the maximum overshooting length and as the simulation time exceeds the typical thermal diffusive time-scale in the overshooting layer.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Numerical and Experimental Investigations on the Ignition Behavior of OME

On the path towards climate-neutral future mobility, the usage of synthetic fuels derived from renewable power sources, so-called e-fuels, will be necessary. Oxygenated e-fuels, which contain oxygen in their chemical structure, not only have the potential to realize a climate-neutral powertrain, but also to burn more cleanly in terms of soot formation. Polyoxymethylene dimethyl ethers (PODE or OMEs) are a frequently discussed representative of such combustibles. However, to operate compression ignition engines with these fuels achieving maximum efficiency and minimum emissions, the physical-chemical behavior of OMEs needs to be understood and quantified. Especially the detailed characterization of physical and chemical properties of the spray is of utmost importance for the optimization of the injection and the mixture formation process. The presented work aimed to develop a comprehensive CFD model to specify the differences between OMEs and dodecane, which served as a reference diesel-like fuel, with regards to spray atomization, mixing and auto-ignition for single- and multi-injection patterns. The simulation results were validated against experimental data from a high-temperature and high-pressure combustion vessel. The sprays’ liquid and vapor phase penetration were measured with Mie-scattering and schlieren-imaging as well as diffuse back illumination and Rayleigh-scattering for both fuels. To characterize the ignition process and the flame propagation, measurements of the OH* chemiluminescence of the flame were carried out. Significant differences in the ignition behavior between OMEs and dodecane could be identified in both experiments and CFD simulations. Liquid penetration as well as flame lift-off length are shown to be consistently longer for OMEs. Zones of high reaction activity differ substantially for the two fuels: Along the spray center axis for OMEs and at the shear boundary layers of fuel and ambient air for dodecane. Additionally, the transient behavior of high temperature reactions for OME is predicted to be much faster.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Sound velocity and compressibility of melts along the hedenbergite (CaFeSi2O6)-diopside (CaMgSi2O6) join at high pressure : Implications for stability and seismic signature of Fe-rich melts in the mantle

Iron-rich silicate melts play an important role in the magmatic history of the Earth and the 16 Moon. However, their elastic properties at high pressures, especially the sound velocity, are poorly 17 understood. Here we determined the ultrasonic sound velocity for the first time of a hedenbergite 18 (Hd, CaFeSi2O6) melt and a melt mixture of 50 mol% hedenbergite + 50 mol% diopside (Hd50Di50) 19 at high pressure and temperature conditions up to 6 GPa and 2329 K, using high-pressure ultrasonic 20 technique combined with synchrotron radiation in a multi-anvil apparatus. Our results show that 21 Fe can significantly reduce the sound velocity while increase the density of silicate melts. 22 Comparing the Di, Hd, and Hd50Di50 melts, we find that the sound velocity does not mix linearly 23 2 for melts in the Hd-Di join, whereas the density for Hd-Di melts at high pressures can be well-24 described by linear mixing. Combined with melt geometry and melt compaction models, we 25 applied our results to study the stability and seismic signature of Fe-rich silicate melts in the Earth’s 26 upper mantle. For the low-velocity zone (LVZ) in mantle asthenosphere, although the degree of 27 seismic velocity reduction can be explained by the presence of a small amount of partial melt 28 distributed in film/band geometry along grain boundaries, Fe-rich melts formed at this depth are 29 unlikely to be gravitationally stable, but may be dynamically unextractable if the melt supply is 30 continuous. For the low-velocity layer (LVL) above the mantle transition zone, the presence of Fe-31 rich melts (with FeO>~10 wt%) distributed in textural equilibrium with the ambient mantle is a 32 plausible explanation.

diopside↗
DOE OSTI2025

Quantum Information for Fusion Energy Sciences (Final Technical Report)

The simulation of plasma dynamics is a critical area of Fusion Energy Sciences (FES) due to it’s usefulness in predicting, controlling, and confining plasmas in the context of potential fusion reactors. The simulation of plasmas is a computationally difficult problem in both classical and quantum physics, motivating investigation into the potential of quantum computers to simulate these systems. This project took several concrete steps towards this goal by developing tools for improving the control, characterization, and calibration of quantum gates on a superconducting quantum computer, developing error suppression and mitigation tools to reduce errors on the quantum computer, and utilizing these advancements to simulate reduced models of plasma dynamics on the quantum computer. In order to efficiently simulate plasma physics, an optimal control method which synthesizes, directly at the pulse level, any quantum gate on qubit and qutrit systems was developed. Using four superconducting transmon quantum processors at Rigetti and LLNL, it was demonstrated that any arbitrary quantum gate on qubits and qutrits could be implemented with high fidelity, leading to a significantly reduced length of a gate sequence. A problem of interest in FES is the nonlinear optical process of laser pulse compression within a plasma. Since quantum physics is linear, simulating nonlinear operations is not naturally feasible on a quantum computer, however it is possible to simulated a quantized version of the nonlinear process. A quantization approach to convert nonlinear wave-wave interaction problems to Hamiltonian simulation problems was developed and demonstrated using two qubits on a Rigetti device. In this experiment, a number of error suppression and mitigation techniques were investigated to determine how best to utilize the finite quantum resources. This study provides an example of how plasma problems may be solved on near-term, noisy quantum computing platforms and identified a promising set of techniques. Building on the insights of these experiments, the investigation turned to linear electron-plasma wave physics. A connection was identified between a local one-dimensional lattice spin model and linear wave phenomena, allowing a plasma physics problem to be efficiently mapped to the quantum computer. In this framework, reflection and transmission of plasma waves at a sharp boundary was studied, as well as the propagation of waves through an inhomogeneous plasma medium. In addition to the suite of error suppression and mitigation techniques developed, this experiment introduced the use of a digital-analog gate scheme designed to efficiently simulate the plasma Hamiltonian. With hardware available at the conclusion of the project, simulation at the scale of 9 qubits and 15 timesteps (60 entangling layers) was achieved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10 −3 to 10 −1 , which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.

79 ASTRONOMY AND ASTROPHYSICS↗