Search NASA⌕ Search

SEARCH · Search NASA

Results for “Speed of sound”

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

Monolithic AlScN/SiC phononic waveguides for scalable acoustoelectric and quantum devices

Unlike conventional surface acoustic wave devices, phononic waveguide systems enable higher circuit density and stronger strain and piezoelectric fields, making them promising for advanced acoustoelectric and quantum applications. One such material system for generating and guiding phonons at gigahertz frequencies is AlScN on SiC, which can be synthesized by sputter depositing AlScN directly onto SiC wafers. The AlScN on the SiC platform allows for tightly vertically-confined acoustic modes with high electromechanical coupling, high speed of sound, and simple fabrication of strip and rib waveguides. Until now, this system has only been studied as a slab waveguide platform, i.e., without any lateral waveguiding. Here, we demonstrate a two-dimensionally confined phononic architecture in AlScN on SiC that supports guided modes at 2.95 and 4.05 GHz. These modes exhibit strong electromechanical coupling coefficients (k 2 = 4.27%) and propagation losses on the order of 10 dB/mm. Furthermore, this architecture is well-suited for phononic routing and power-efficient active or nonlinear devices such as amplifiers, mixers, and oscillators, and is compatible with the integration of quantum systems, including vacancy centers, charge carriers, photons, and spins, either embedded in SiC or heterogeneously integrated on the surface.

Electrical components↗
DOE OSTI2024

Air-coupled tsunamis generated from impacts and airbursts: Our understanding before Hunga-Tonga Hunga-Ha'apai

The effort to prevent or mitigate the effects of an impact on Earth is known as planetary defense. A significant component of planetary defense research involves risk assessment. Much of our understanding of the risk from near-Earth objects comes from the geologic record in the form of impact craters, but not all asteroid impacts are crater-forming events. Small asteroids explode before reaching the surface, generating an airburst, and most impacts into the ocean do not penetrate the water to form a crater in the sea floor. The risk from these non-crater-forming ocean impacts and airbursts is difficult to quantify and represents a significant uncertainty in our assessment of the overall threat. We are currently working to better understand impact scenarios that can generate dangerous tsunamis. One of the suggested mechanisms for the production of asteroid–generated tsunamis is by direct coupling of the pressure wave to the water, analogous to the means by which a moving weather front can generate a meteotsunami. To test this hypothesis, we ran a series of airburst simulations and provided time-resolved pressure and wind profiles to use as source functions for tsunami models. We used the CTH hydrocode to model the various airburst scenarios to compare to the results of other simulations and provide time dependent boundary conditions as input to shallow-water wave propagation codes. The strongest and most destructive meteotsunamis are generated by atmospheric pressure oscillations with amplitudes of only a few hPa1 (mbar), corresponding to changes in sea level of a few cm. The resulting wave is strongest when there is a resonance between the ocean and the atmospheric forcing. A Proudman resonance takes place when the atmospheric disturbance’s translational speed (U) equals the longwave phase speed $\sqrt{gh}$ of shallow water wave. Coupling is strongest when the Froude number (Fr=U/c) is unity. A weather front propagates much slower than the speed of sound, so meteotsunamis are most common and dangerous in shallow bodies of water such as the Mediterranean Sea or Lake Michigan. By contrast, the blast wave from an airburst or crater-forming impact propagates at a speed faster than a tsunami in the deepest ocean, and a Proudman resonance cannot be achieved even though the overpressures are orders of magnitude greater. However, blast wave profiles are N-waves in which a sharp shock wave leading to overpressure is followed by a more gradual rarefaction to a much longer-duration underpressure phase. Even though the blast outruns the water wave it is forcing, the tsunami should continue to be driven by the out-of-resonance gradient associated with the suction phase, which may depend strongly on the details of the airburst or impact scenario. The open question is whether there are any conditions under which such an airburst-driven tsunami can be dangerous enough to contribute to the overall impact risk. We have also identified other potential mechanisms for airburst-generated tsunamis: 1) reaction force at the surface from the plume ejected into space, which carries significant momentum, 2) expanding toroidal vortices at the surface, which travel more slowly than the shock wave and can generate a Proudman resonance in relatively shallow ocean (such as continental shelf), and 3) steam explosion from seawater ablation by a “Type II” (Libyan Desert Glass-type) airburst in which the hot vapor jet descends to the surface. On January 15, 2022, the Hunga-Tonga Hunga-Ha’apai volcano, located approximately 60 km north of Tongatapu, the main island of Tonga, violently erupted with a powerful explosion, culminating the period of volcanic activity that started in December of 2021. This event and resulting tsunamis provided an existence proof for the air pressure wave coupling mechanism we proposed. It also suggests that it can be stronger and more significant over much greater distances than we contemplated, leading to global tsunamis associated with impact events on land as well as in the water. Large atmospheric explosions generate global Lamb waves with larger amplitudes, longer periods, and slower speeds than the local and regional blast waves we modeled prior to that event. This paper reviews our analysis and modeling of airburst-driven tsunamis prior to the 2022 Hunga-Tonga Hunga-Ha’apai tsunami, which was the subject of two presentations at the 2023 Planetary Defense Conference and is the subject of another paper currently in preparation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Aero-Optics of Hypersonic Turbulent Boundary Layers

Aero-optics refers to optical distortions due to index-of-refraction gradients that are induced by aerodynamic density gradients. At hypersonic flow conditions, the bulk velocity is many times the speed of sound and density gradients may originate from shock waves, compressible turbulent structures, acoustic waves, thermal variations, etc. Due to the combination of these factors, aero-optic distortions are expected to differ from those common to sub-sonic and lower super-sonic speeds. This report summarizes the results from a 2019-2022 Laboratory Directed Research and Development (LDRD) project led by Sandia National Laboratories in collaboration with the University of Notre Dame, New Mexico State University, and the Georgia Institute of Technology. Efforts extended experimental and simulation methodologies for the study of turbulent hypersonic boundary layers. Notable experimental advancements include development of spectral de-aliasing techniques for highspeed wavefront measurements, a Spatially Selective Wavefront Sensor (SSWFS) technique, new experimental data at Mach 8 and 14, a Quadrature Fringe Imaging Interferometer (QFII) technique for time-resolved index-of-refraction measures, and application of QFII to shock-heated air. At the same time, model advancements include aero-optic analysis of several Direct Numerical Simulation (DNS) datasets from Mach 0.5 to 14 and development of wall-modeled Large Eddy Simulation (LES) techniques for aero-optic predictions. At Mach 8 measured and predicted root mean square Optical Path Differences agree within confidence bounds but are higher than semi-empirical trends extrapolated from lower Mach conditions. Overall, results show that aero-optic effects in the hypersonic flow regime are not simple extensions from prior knowledge at lower speeds and instead reflect the added complexity of compressible hypersonic flow physics.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Decoding THz‐Driven Dynamic Fingerprints of Ferroelectric Nanotwin Networks

Ultrafast polarization dynamics in ferroelectrics are of considerable interest for high-speed tunable dielectrics and electro-optics. Extended domain wall networks formed in ferroelectric twin nanodomains can support collective dynamics in the terahertz regime but require techniques that track polarization and strain evolution driven by ultrafast stimulus. Here, we use multi-modal probing of THz-pulse-driven excitations in PbTiO 3 /SrTiO 3 superlattices by combining X-ray free electron laser measurements that directly tracks lattice changes, with optical second harmonic generation that tracks the electronic potential coupled with the lattice potential. Dynamical phase-field modeling enables fingerprinting of these collective modes as superpositions of domain “breathing” through wall oscillations and polarization “rotations” with still walls. Ultrafast domain wall motion at 0.1–0.5 THz is observed at practical fields of 100 kV/cm with wall velocities of >4000 m/s, approaching typical speed of sound in PbTiO 3 . A unique “charging” mode is discovered that can electrically charge and discharge domain walls on ∼4 ps time scale thus dynamically tuning wall conductivity. Integrated experimental and theoretical fingerprinting of the dynamical landscape presented here enables ultrafast control of ferroics for high-speed microelectronics and optical applications.

THz dynamics↗
DOE OSTI2020

Development of Methods for In-line Monitoring of Yield Stress During the Transfer of Radioactive Waste - 20435

Radioactive waste slurries at the Savannah River Site (SRS) and Hanford often behave as non-Newtonian fluids with a yield stress. The yield stress increases the size of equipment needed to transport or mix these slurries and can increase the risk of process upsets. This property is often measured by collection of samples and transportation to a laboratory, which can increase the risk of personnel exposure and even change the rheology. In addition, slurry rheology has been observed to change with time. Real-time in-line monitoring (RTIM) of rheological properties of slurries can have substantial benefits to critical processes which are currently in operation or will be part of future plans in Department of Energy (DOE) complex processes. The need for an appropriate RTIM technology that can potentially minimize or eliminate numerous sampling operations was recently identified. This research at Florida International University focuses on addressing the need for a reliable RTIM technology at US Department of Energy sites. The objective of this work is to investigate and develop appropriate and reliable technologies that can perform accurate and easy measurements of yield stress with minimal changes to the flow conditions at US Department of Energy sites. Use of ultrasound in longitudinal and shear modes was considered as a method for monitoring of yield stress variations with no disturbance to the flow. Measurement of shear wave velocity (in shear mode) and monitoring of sound velocity variations with start of flow from rest (longitudinal mode) are the underlying principles. Measurement of pressure loss and liquid rise (due to static differential pressure) during flow of slurries were two other methods that were considered and investigated. In addition, analysis of stability depth of a penetrometer in slurry medium (penetrometer method) revealed a potential suitability. This paper presents efforts associated with investigation of different methods mentioned earlier for kaolin-water mixtures at different concentrations. Theoretical analysis showed that all methods under investigation were capable of monitoring yield stress in a desirable range of 0 to 70 pascal. A bench scale test setup was constructed for investigation of pressure loss and liquid rise methods. Simulants at different concentrations were created for initial testing campaigns. For ultrasonic testing, effects of yield stress on sound propagation speed in longitudinal and shear modes were investigated for simulants placed inside PVC pipes. Parameters such as simulant concentration and pipe diameter were varied to find optimal values. It was found that spacing (or medium thickness) of 1/8 inch or less was necessary to observe echoes with sufficient strength. In addition, a significant effect of yield stress on sound propagation velocity was observed. For the liquid rise method, kaolin-water simulant at 1.17 specific gravity was pumped through a 1/4 inch pipe and liquid height was measured in a vertical riser branching off a three-way connection. Static gauge pressure was varied within [0.3 to 1.1] psi range by increase of pumping pressure and flow. Results showed consistent yield stress data within 10% deviation from values reported by a rheometer. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Ultrasonic Measurements of Temperature Profile and Heat Fluxes in Coal-Fired Power Plants (Final Report)

Many industrial processes are inaccessible or inhospitable to characterization by traditional temperature measurement methods, such as thermocouples, especially over prolonged exposure to harsh environments. Ultrasound is an established characterization technology with diverse applications ranging from medical imaging to therapies to flaw detection to nondestructive evaluation. Ultrasound may characterize solid materials and components noninvasively as a nondestructive evaluation modality and obtain internal measurements of material properties. For example, the speed of ultrasound propagation changes with Young’s modulus and Poisson’s ratio, which can be found from its measurements. Traditional ultrasonic characterization assumes all material properties remain constant with the position. When this assumption holds, a property of interest may be measured by relating it to the speed of ultrasound propagation (or a speed of sound, SOS) and measuring the SOS by timing the ultrasound propagation through a known distance. However, when a property of interest is spatially distributed, the propagation time depends on the SOS changing with the position along the ultrasound propagation path. The multiple temperature distributions may lead to an identical time of flight (TOF). Temperature is one property that impacts the speed of ultrasound and often cannot be assumed to remain constant with the position. Previously, in the context of temperature, we addressed the challenge of ultrasonic characterization of spatially distributed properties by developing a method for measuring segmental temperature distributions (MSTD). This method divides the ultrasonic propagation into segments bound by echogenic features. These features provide ultrasonic interfaces where some energy is reflected toward the receiving transducer, and the rest continues through the medium. The time-of-flight between the echoes reflected from echogenic features characterizes the spatial distribution in the properties of interest in the corresponding segment of the ultrasonic propagation path. This project demonstrated the application of the MSTD method in industrial conditions of the coal-fired power plant. We implemented the MSTD using metals and alloys waveguides, which may be the existing structure for which the temperature distribution is characterized or purposefully designed waveguides added to the structure by welding or other means specifically to quantify thermal properties using the MSTD method. Previous iterations of the MSTD method used ceramic and cementitious waveguides, which significantly attenuate ultrasound. On the other hand, low attenuation in metallic waveguides creates interactions between echogenic features which compilates the signal analysis in the segmental TOF measurements. We have established the WG design principles that minimize the interferences between trailing and primary echoes and, in some cases, eliminate them. The waveguides in which echoes do not interfere improve the timing accuracy and the robustness of ultrasonic measurements of the spatial distributions in material properties. Our emphasis remained on the estimation of the temperature distributions. We have developed general recommendations for designing ultrasonically segmented waveguides with the reduced influence of trailing echoes. Two of our waveguide designs were tested in the industry. The first waveguide was designed for insertion into a combustion zone of the utility-scale coal-fired power plant boiler. The second design allows the characterization of temperature distribution in the direction normal to the boiler’s water wall, a large heat exchanger converting the chemical energy released during combustion to the steam driving the electrical power generation turbines. These waveguides were designed to operate within a restrictive space of thermally insulated water wall and incorporate densely located echogenic features while combatting the influence of trailing echoes. The project has successfully demonstrated the feasibility of using the developed method for accurate, continuous, and robust temperature measurements in extreme environments of power generation and other industrial processes. It, therefore, has achieved its overarching goal of advancing the technology readiness level of the novel Ultrasound Measurements of Segmental Temperature Distribution (US-MSTD) method for real-time measurements of the temperature distribution and heat fluxes closer to commercial availability, developing a prototype multipoint measurement system, and validating its performance on coal-fired utility boilers. The success of this project was achieved in collaboration with the power generator, Rocky Mountain Power, and set the stage for the transfer of this technology from the laboratory to the industry.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2025

Reference lattice, sound, stiffness, and magnetic transitions of Ising monolayers

A reference lattice, away from which elastic distortions induced by the spin texturing of two-dimensional (2D) magnets take hold, is motivated from a picture of pairwise Biot-Savart interactions among identical solenoids that either elongate or compress a (“zero-current”) spring lattice. Applied to a paradigmatic CrSiTe 3 monolayer, the reference is given by the average between the atomic positions of ferromagnetic (FM) and Néel antiferromagnetic (AFM) lattices; such an atomic disposition permits understanding structural distortions and elastic energies due to magnetism readily. Furthermore, the anisotropic speed of sound in the magnetic ground state explains an observed anisotropy of vibrational frequencies on similar magnets. Elastic stiffness constants are reported, too. Magnetic energies in four Ising structural configurations were calculated and the strain needed for those 2D magnets to undergo an AFM to FM quantum phase transition was determined as well.

density functional theory↗
DOE OSTI2022

Plasma production and composition from hypervelocity impacts on solar cell cover glass

Satellites are subjected to detrimental effects from the space environment, including impacts by orbital debris and meteoroids that can cause mechanical or electrical damage. While mechanical damage is well studied, electrical effects resulting from hypervelocity impact (HVI) remain poorly understood. In HVIs, where the projectile speed exceeds the speed of sound in the target material, the impact energy ionizes the material near the surface, creating a plasma that can cause radio frequency emission. This emission is the main source of electrical anomalies in satellites, but the mechanism behind its creation is still unknown. Ground-based experiments have been used to characterize HVI plasma, including the empirical power law relations that describe the impact charge produced as a function of impactor mass and speed. Quantifying these power law relations gives us a better understanding of the species formed in the impact which can eventually lead to plasma evolution models and help us first, understand the mechanism behind radio frequency emission and second, ways to prevent its damage to the satellites. In this paper, we interpret the results of the HVI experiments on solar cell coverglass targets and identify the potential chemical or molecular reactions that form the impact ion species. This is done by looking simultaneously at multiple data sets for different species within the same impact configuration, and using a new technique for quantifying the power law relations for each emitted species. We use time of flight analysis (TOF) to detect the charged particle compositions and determine their initial speeds. The results suggest that the charge production power law is species specific, with different formation mechanisms involved in a single HVI event, some of which are associated with dissociation and some with ionization of a molecule. Furthermore, our results also show that electrons have different power laws that depend on the spacecraft surface charge. TOF measurements of ions from impact events on solar cell cover glass shows the presence of SiO + and SiO – 2 as the dominant species present in plasmas on positively and negatively charged targets, respectively, as well as the formation of FeO + .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Transonic dislocation propagation in diamond

Here, the motion of line defects (dislocations) has been studied for more than 60 years, but the maximum speed at which they can move is unresolved. Recent models and atomistic simulations predict the existence of a limiting velocity of dislocation motion between the transonic and subsonic ranges at which the self-energy of dislocation diverges, though they do not deny the possibility of the transonic dislocations. We used femtosecond x-ray radiography to track ultrafast dislocation motion in shock-compressed single-crystal diamond. By visualizing stacking faults extending faster than the slowest sound wave speed of diamond, we show the evidence of partial dislocations at their leading edge moving transonically. Understanding the upper limit of dislocation mobility in crystals is essential to accurately model, predict, and control the mechanical properties of materials under extreme conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

An acoustic investigation of the near-surface turbulence on Mars

The Perseverance rover is carrying out an original acoustic experiment on Mars: the SuperCam microphone records the spherical acoustic waves generated by laser sparks at distances from 2 m to more than 8 m. These N-shaped acoustic waves scatter from the multiple local heterogeneities of the turbulent atmosphere. Therefore, large and random fluctuations of sound travel time and intensity develop as the waves cross the medium. The variances of the travel times and the scintillation index (normalized variance of the sound intensity) are studied within the mathematical formalism of the propagation of spherical acoustic waves through thermal turbulence to infer statistical properties of the Mars atmospheric temperature fluctuation field. The comparison with the theory is made by simplifying assumptions that do not include wind fluctuations and diffraction effects. Two Earth years (about one Martian year) of observations acquired during the maximum convective period (10:00–14:00 Mars local time) show a good agreement between the dataset and the formalism: the travel time variance diverges from the linear Chernov solution exactly where the density of occurrence of the first caustic reaches its maximum. Moreover, on average, waves travel faster than the mean speed of sound due to a fast path effect, which is also observed on Earth. To account for the distribution of turbulent eddies, several power spectra are tested and the best match to observation is obtained with a generalized von Karman spectrum with a shallower slope than the Kolmogorov cascade, φ(k)∝(1 + k 2 L 2 ) -4/3 . It is associated with an outer scale of turbulence, L, of 11 cm at 2 m above the surface and a standard deviation of 6 K over 9 s for the temperature. These near-surface atmospheric properties are consistent with a weak to moderate wave scattering regime around noon with little saturation. Overall, this study presents an innovative and promising methodology to probe the near-surface atmospheric turbulence on Mars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A Goldstone boson equivalence for inflation

The effective field theory of single-field inflation characterizes the inflationary epoch in terms of a pattern of symmetry breaking. An operator acquires a time-dependent vacuum expectation value, defining a preferred spatial slicing. In the absence of dynamical gravity, the fluctuations around the time-dependent background are described by the Goldstone boson associated with this symmetry breaking process. With gravity, the Goldstone is eaten by the metric, becoming the scalar metric fluctuation. In this paper, we will show that in general single-field inflation, the statistics of scalar metric fluctuations are given by the statistics of this Goldstone boson decoupled from gravity up to corrections that are controlled as an expansion in slow-roll parameters. This even holds in the presence of additional parameters, like the speed of sound, that naively enhance the impact of the gravitational terms. In the process, we derive expressions for leading and sub-leading gravitational corrections to all-orders in the Goldstone boson.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Hydrodynamics of simply spinning black holes & hydrodynamics for spinning quantum fluids

We find hydrodynamic behavior in large simply spinning five-dimensional Anti-de Sitter black holes. These are dual to spinning quantum fluids through the AdS/CFT correspondence constructed from string theory. Due to the spatial anisotropy introduced by the angular momentum, hydrodynamic transport coefficients are split into groups longitudinal or transverse to the angular momentum, and aligned or anti-aligned with it. Analytic expressions are provided for the two shear viscosities, the longitudinal momentum diffusion coefficient, two speeds of sound, and two sound attenuation coefficients. Known relations between these coefficients are generalized to include dependence on angular momentum. The shear viscosity to entropy density ratio varies between zero and 1/(4π) depending on the direction of the shear. These results can be applied to heavy ion collisions, in which the most vortical fluid was reported recently. In passing, we show that large simply spinning five-dimensional Myers-Perry black holes are perturbatively stable for all angular momenta below extremality.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Equation of State for the Thermodynamic Properties of Trans-1,2-dichloroethene [R-1130(E)]

We present an empirical equation of state in terms of the Helmholtz energy for trans-1,2-dichloroethene [R-1130(E)]. The range of validity is from the triple-point temperature, 223.31 K to 525 K with pressures up to 30 MPa. It may be used to calculate all thermodynamic properties in the fluid phase, including liquid, gas, and supercritical regions. Comparisons are given with existing literature data and estimated uncertainties are provided. In addition, checks were made for correct extrapolation behavior so that the equation behaves in a physically realistic manner when used outside of its range of validity, enabling its use in mixture models. The estimated uncertainties (at a k = 2 or 95 % level of confidence) are based on comparisons with critically assessed data and are 0.25 % for vapor pressure for temperatures in the range 300 K < T < 454 K, rising to 1.5 % as the temperature decreases from 300 K to 265 K. For density in the liquid phase the estimated uncertainty is 0.14 % for temperatures 270 K < T < 410 K and for pressures up to 30 MPa. For the vapor phase the estimated uncertainty in density is 3 %. The uncertainty for liquid-phase heat capacity is 1 % at atmospheric pressure over the temperature range 268 K < T < 309 K, and the uncertainty for the speed of sound in the liquid phase is 0.25 % for temperatures 230 K < T < 420 K and for pressures up to 30 MPa. The uncertainties are larger outside of these specified ranges and in the critical region.

1,2-Dichloroethene↗
DOE OSTI2025

Electronic Melting of Silicon in Nanostructures using X-ray Forbidden Bragg Reflections

Abstract We carried out a short beamtime at the Pohang Accelerator Laboratory x-ray Free Electron Laser to perform a pump-probe (PP) laser excitation diffraction experiment on the silicon (222) forbidden Bragg peak. To limit the x-ray penetration, we used a “device layer” silicon film wafer bonded to a silicon substrate. The sample, specially fabricated by MEMC Electronic Materials, had a Si(100) substrate bonded to a 170 nm Si(100) film rotated at 45° for crystallographic isolation. A second sample was reactive-ion-etched down to 52 nm thickness. In the silicon lattice, the covalent bonds are seen exclusively at the 222 reflection. Upon laser excitation, these electrons are expected to be excited to the valence band on femtosecond electronic time scales. The Si(222) reflection is therefore expected to be extinguished on this fast time scale, while the electron–phonon coupled acoustic response is determined by the lattice dynamics. The latter is determined by the speed of sound over the device thickness, which is in the mid-picosecond range.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Limiting velocities and transonic dislocations in Mg

To accurately predict the mechanical response of materials, especially at high strain rates, it is important to account for dislocation velocities in these regimes. Under these extreme conditions, it has been hypothesized that dislocations can move faster than the speed of sound. However, the presence of such dislocations remains elusive due to challenges associated with measuring these experimentally. In this work, molecular dynamics simulations were used to investigate the dislocation velocities for the basal edge, basal screw, prismatic edge, and prismatic screw dislocations in Mg in the sub-, trans-, and supersonic regimes. Also, our results show that only prismatic edge dislocations achieve supersonic velocities. Furthermore, this work demonstrates that the discrepancy between the theoretical limiting velocity and the MD results for Mg is due to its sensitivity to large hydrostatic stress around the dislocation core, which was not the case for fcc metals such as Cu.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Semi-implicit method to solve compressible multiphase fluid flows without acoustic time step restrictions

A new 3D multiphase numerical capability is presented here for simulating multiphase flow regimes at all Mach numbers (M). The new method is a semi-implicit extension of the finite volume discrete equation method (DEM) of Chinnaya et al. (2004), which originally used explicit time-stepping. The capability is also developed to work with another extension of the DEM to moving grids for arbitrary Lagrangian-Eularian (ALE) methods detailed in Dunn (2011). Rather than solving all phase equations simultaneously, the DEM reduces the equations to a system of single-phase Riemann solves, where each phase has its own velocity and thermodynamic state. Exchanges between the phases are modeled through source terms accounting for the phase interactions. Furthermore, since the original multiphase scheme uses an explicit time-advancement scheme, it has time step restrictions dictated by the speed of sound, which limits the model's ability to simulate weakly compressible flows which typically need to be integrated for longer time periods. Here, we extend the current multiphase formulation by implementing a pressure-correcting step to enable implicit calculations and remove acoustic time step limitations. The new semi-implicit algorithm allows use of relatively large time steps compared to an explicit method. Validation and benefits of the new implicit time-step method are illustrated using several examples including weakly compressible flows and strong shock waves. The scheme presented here is general and may be used for a variety of applications which require capabilities for handling multiphase flow at a wide range of Mach numbers. However, the goal of this research is to simulate all stages of high energy explosions, including the shock formation (high Mach numbers) and evolution of the buoyant cloud (low Mach numbers).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

High pressure/high temperature multiphase simulations of dodecane injection to nitrogen: Application on ECN Spray-A

The present work investigates the complex phenomena associated with pressure/high temperature dodecane injection for the Engine Combustion Network (ECN) Spray-A case, employing more elaborate thermodynamic closures, to avoid well known deficiencies concerning density and speed of sound prediction using traditional cubic models. A tabulated thermodynamic approach is proposed here, based on log 10 (p)-T tables, providing very high accuracy across a large range of pressures, spanning from 0 to 2500 bar, with only a small number of interpolation points. The tabulation approach is directly extensible to any thermodynamic model, existing or to be developed in the future. Here NIST REFPROP properties are used, combined with PC-SAFT Vapor-Liquid-Equilibrium to identify the liquid in mixtures penetration, hence avoiding the use of an arbitrary threshold for mass fraction. Identified liquid and vapour penetration are compared against experimental data from the ECN database showing a good agreement, within approximately 3–8% for axial penetration of liquid, 2% for vapor axial penetration and within experimental uncertainty for radial distribution of mass fraction. Analysis of the vortex evolution indicates that driving mechanisms behind the jet break-up are vortex tilting/stretching, then baroclinic torque, leading to Rayleigh-Taylor instabilities, closely followed by vortex dilation and finally viscous effects.

42 ENGINEERING↗
DOE OSTI2023

Preliminary study of radio frequency waves in hypervelocity impact plasma

Orbital debris and meteoroids can cause mechanical and electrical anomalies by impacting spacecraft. Although mechanical damage tends to occur only from particles with masses >10 –3 g, electrical damage can be a result of hypervelocity (> 10 km/s) impacts (HVI) even with smaller particles. In HVI, the projectile speed exceeds the speed of sound in the material and has enough energy to create plasma by ionizing the material near the surface, which can then emit radio frequency (RF) waves. In this paper, we conduct a preliminary study on the characteristics of impact-generated RF emission with the ultimate goal of building an RF predictive model for risk assessment of HVI electrical effects that can damage spacecraft. First, we introduce wavelet analysis as a method for denoising and quantifying the total RF signal energy. An appropriate mother wavelet and the number of decomposition levels are selected using a synthetic signal. Second, the total energy of the RF emission is calculated using Parseval’s theorem. Results suggest that the total RF signal energy is proportional to the energy flux of the impactor. Lastly, spectrograms of the RF emission are analyzed to characterize features that can be used in the future to predict impacts that can cause electrical damage to satellites.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗