Search NASA⌕ Search

SEARCH · Search NASA

Results for “wave physics”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A gas dynamic perspective on particle lifting in electrostatic discharge-like devices using multiphase particle-in-cell approach

Although particle–laden electrostatic discharges are widely used in laboratories as well as in industrial applications, the mechanism of particle lifting for particles initially at rest in such highly unsteady systems is not well understood. A multiphase gas–particle solver is developed using the multiphase particle-in-cell (MP-PIC) approach to emulate the interaction of a compressible shock-dominated gas phase with the dense particle phase. First, the two-way coupled solver is initially used to study the interaction of a planar traveling shock with a vertical curtain of particulates. The gas and particle phase evolution was found to be in good agreement with a similar experimental study in Ling et al. [Phys. Fluids 24, 113301 (2012)]. Second, the MP-PIC code is used to study the interaction of an expanding blast wave with a thick bed of particles. The simulation considered forces such as quasi-steady drag, pressure-gradient, added-mass, Saffman, and Magnus forces. We observe that the vertical liftoff particles close to the shock impingement point in this configuration are associated with the quasi-steady drag, pressure gradient, and added-mass forces. Also, the Saffman lift and Magnus forces contribute to lifting particles located radially farther away from the shock impingement point. In addition, the study finds a decrease in particle lifting efficiency with decreasing plasma kernel length and shock strength.

42 ENGINEERING↗
DOE OSTI2020

Simultaneous Observations of Localized and Global Drift Resonance

In this study, we present Van Allen Probe observations showing that seed (hundreds of keV) and core ( ≳ 1 MeV) electrons can resonate with ultra-low-frequency (ULF) wave modes with distinctive m values simultaneously. An unusual electron energy spectrogram with double-banded resonant structure was recorded by energetic particle, composition, and thermal plasma (ECT)-magnetic electron ion spectrometer (MagEIS) and, meanwhile, boomerang stripes in pitch angle spectrogram appeared at the lower energy band. A localized drift resonance with m = 10 wave component was responsible for the resonant band peaked at ~200 keV while a global drift resonance with m = 3 component gave rise to the upper band resonance peaked at ~1 MeV. Time-Of-Flight on boomerang stripes suggested that the localized drift resonance with ~200 keV electrons was confined within the plasmaspheric plume. Electron flux modulations were reproduced by numerical simulations in good consistency with the observations, supporting the scenario that localized and global drift resonance could coexist in the outer belt electron dynamics simultaneously.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Suppressing the enhancement of stimulated Raman scattering in inhomogeneous plasmas by tuning the modulation frequency of a broadband laser

The stimulated Raman scattering (SRS) instability can inhibit the performance of laser driven inertial confinement fusion (ICF) implosions by scattering light into unwanted directions or by generating hot electrons that preheat the target fuel. In principle, ICF target designs can avoid parameter regimes conducive to large, linear SRS gains. In practice, kinetic inflation—the nonlinear enhancement of SRS due to electron trapping in the excited plasma wave—makes this difficult. Here we show that laser bandwidth in the form of frequency modulation can either decrease or increase the inflationary SRS (iSRS) threshold in inhomogeneous plasmas depending on the maximum chirp of the laser pulse. The threshold, mapped out by a series of particle-in-cell simulations, exhibits a minimum when the frequency change within the pulse cancels the spatial detuning due to density inhomogeneities along the trajectory of the scattered light. In conclusion, by tuning the pump laser parameters away from this minimum, the iSRS threshold can be larger than at zero bandwidth, providing a path to mitigating kinetic inflation in ignition relevant plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Transients following the loss of detonation confinement

We present a theory for the evolution of a one-dimensional, steady-state detonation reaction zone to a two-dimensional reaction zone, when the explosive experiences a sudden loss of side-on confinement as a boundary of the explosive is impulsively withdrawn. Our focus is on condensed-phase explosives, which we describe as having a constant adiabatic gamma equation of state and an irreversible, state-independent reaction rate. We consider two detonation models: (i) the instantaneous reaction heat release Chapman–Jouguet (CJ)-limit and (ii) the spatially resolved reaction heat-release Zel’dovich–von Neumann–Döring (ZND) model, in the limit where only a small fraction of the energy release is resolved (the SRHR-limit). Two competing rarefaction waves are generated by this loss of confinement: (i) a smooth wave coming off the full length of the withdrawn boundary and (ii) a singular fan spreading out from the point where the detonation shock and the withdrawn boundary meet. For the CJ-limit, in all cases the singular rarefaction fan eventually dominates the competition to control the steady-state behaviour. For the SRHR-limit, the spatially resolved heat release moderates this competition. When the withdrawal speed is fast, the rarefaction fan dominates; when the withdrawal speed is slower, the smooth rarefaction eventually dominates, although the flow features a fan at early times. By examining the mathematical properties of the steady two-dimensional fan-based solution, we set down a mechanism for this transition in behaviours.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A numerical study of particle jetting in a dense particle bed driven by an air-blast

In this work, the particle jetting behavior in a blast-driven dense particle bed is studied at early times. Four-way coupled Euler–Lagrange simulations are performed using a high-order discontinuous Galerkin spectral element solver coupled with a high-order Lagrangian particle solver, wherein the inter-particle collisions are resolved using a discrete element method collision model. Following the experiments of Rodriguez et al. [“Formation of particle jetting in a cylindrical shock tube,” Shock Waves 23(6), 619–634 (2013)] and the simulations of Osnes et al. [“Numerical simulation of particle jet formation induced by shock wave acceleration in a Hele-Shaw cell,” Shock Waves 28(3), 451–461 (2018)], the simulations are performed in a quasi-two-dimensional cylindrical geometry (Hele-Shaw cell). Parametric studies are carried out to assess the impact of the coefficient of restitution and the strength of the incident shock on the particle jetting behavior. The deposition of vorticity through a multiphase (gas–particle) analog of Richtmyer–Meshkov instability is observed to play a crucial role in channeling the particles into well-defined jets at the outer edge of the particle bed. This is confirmed by the presence of vortex pairs around the outer jets. Furthermore, the effect of the relaxation of the relative velocity between the two phases on the vorticity generation is explored by analyzing the correlation between the radial velocity of particles and the radial velocity of the gas at the particle location.

42 ENGINEERING↗
DOE OSTI2024

Small-angle scattering interferometry with neutron orbital angular momentum states

Methods to prepare and characterize neutron helical waves carrying orbital angular momentum (OAM) were recently demonstrated at small-angle neutron scattering (SANS) facilities. These methods enable access to the neutron orbital degree of freedom which provides new avenues of exploration in fundamental science experiments as well as in material characterization applications. However, it remains a challenge to recover phase profiles from SANS measurements. We introduce and demonstrate a novel neutron interferometry technique for extracting phase information that is typically lost in SANS measurements. An array of reference beams, with complementary structured phase profiles, are put into a coherent superposition with the array of object beams, thereby manifesting the phase information in the far-field intensity profile. We demonstrate this by resolving petal-structure signatures of helical wave interference for the first time: an implementation of the long-sought recovery of phase information from small-angle scattering measurements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Numerical investigation of laser-driven shock interaction with a deformable particle

A laser-driven shock propagating through an isolated particle embedded in a plastic (CH) target was studied using the radiation-hydrodynamic code FLASH. Preliminary simulations using IONMIX equations of state (EOS) showed significant differences in the shock Hugoniot of aluminum compared to experimental data in the low-pressure regime [O(10) GPa], resulting in higher streamwise compression and deformation of an aluminum particle. Hence, a simple modification to the ideal gas EOS was developed and employed to describe the target materials and examine the particle dynamics. The evolution of the pressure field demonstrated a complex wave interaction, resulting in a highly unsteady particle drag which featured two drag minima due to shock focusing at the rear end of the particle and rarefaction stretching due to laser shut-off. Although ~30% lateral expansion and ~25% streamwise compression were observed, the aluminum particle maintained considerable integrity without significant distortion. Additional simulations examined the particle response for a range of particle densities, sizes, and acoustic impedances. The data revealed that lighter particles such as aluminum gained significant momentum, reaching up to ~96% of the shocked CH's speed, compared to ~29% for the heavier tungsten particles. Despite the differences seen in the early stage of shock interaction, particles with varying acoustic impedances ultimately reached the same peak velocity. This identified particle-to-host density ratio is an important factor in determining the inviscid terminal velocity of the particle. In addition, the modified EOS model presented in this study could be used to approximate solid materials in hydrocodes that lack material strength models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Kinetic study of strong shock waves in fully ionized plasmas

Fully kinetic simulations are used in this work to study strong shock waves in fully ionized plasmas for charge state 𝑍 ranging from 1 to 18. The ion stress tensor amplitude agrees well with classical predictions within the compression layer; however, its spatial structure is noticeably shaped by nonlocal transport effects. The electron stress tensor remains significantly smaller than that of the ions for all values of 𝑍 considered. Ion overheating in the compression layer is observed to scale approximately linearly with 𝑍. Both findings diverge from hydrodynamic models based on classical electron stress formulations. In conclusion, the electrons exhibit strong nonlocal behavior in the compression layer, with Knudsen numbers approaching unity for 𝑍 ≥ 6.

Physics - Plasma physics↗
DOE OSTI2024

A three-point velocity estimation method for two-dimensional coarse-grained imaging data

Time delay and velocity estimation methods have been widely studied subjects in the context of signal processing, with applications in many different fields of physics. The velocity of waves or coherent fluctuation structures is commonly estimated as the distance between two measurement points divided by the time lag that maximizes the cross correlation function between the measured signals, but this is demonstrated to result in erroneous estimates for two spatial dimensions. We present an improved method to accurately estimate both components of the velocity vector, relying on three non-aligned measurement points. We introduce a stochastic process describing the fluctuations as a superposition of uncorrelated pulses moving in two dimensions. Using this model, we show that the three-point velocity estimation method, using time delays calculated through cross correlations, yields the exact velocity components when all pulses have the same velocity. The two- and three-point methods are tested on synthetic data generated from realizations of such processes for which the underlying velocity components are known. The results reveal the superiority of the three-point technique. Finally, we demonstrate the applicability of the velocity estimation on gas puff imaging data of strongly intermittent plasma fluctuations due to the radial motion of coherent, blob-like structures at the boundary of the Alcator C-Mod tokamak.

Materials Science↗
DOE OSTI2020

GROWTH on S190814bv: Deep Synoptic Limits on the Optical/Near-infrared Counterpart to a Neutron Star–Black Hole Merger

On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star–black hole (NSBH) merger, or potentially a low-mass binary BH merger. Due to the low false-alarm rate and the precise localization (23 deg 2 at 90%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to an NSBH merger. To search for potential counterparts, the GROWTH Collaboration performed real-time image subtraction on six nights of public Dark Energy Camera images acquired in the 3 weeks following the merger, covering >98% of the localization probability. Using a worldwide network of follow-up facilities, we systematically undertook spectroscopy and imaging of optical counterpart candidates. Combining these data with a photometric redshift catalog, we ruled out each candidate as the counterpart to S190814bv and placed deep, uniform limits on the optical emission associated with S190814bv. For the nearest consistent GW distance, radiative transfer simulations of NSBH mergers constrain the ejecta mass of S190814bv to be M ej < 0.04 M ⊙ at polar viewing angles, or M ej < 0.03 M ⊙ if the opacity is κ < 2 cm 2 g -1 . Assuming a tidal deformability for the NS at the high end of the range compatible with GW170817 results, our limits would constrain the BH spin component aligned with the orbital momentum to be χ < 0.7 for mass ratios Q < 6, with weaker constraints for more compact NSs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A reduced-order model for nonlinear radiative transfer problems based on moment equations and POD-Petrov-Galerkin projection of the normalized Boltzmann transport equation

A data-driven projection-based reduced-order model (ROM) for nonlinear thermal radiative transfer (TRT) problems is presented. The TRT ROM is formulated by (i) a hierarchy of low-order quasidiffusion (aka variable Eddington factor) equations for moments of the radiation intensity and (ii) the normalized Boltzmann transport equation (BTE). The multilevel system of moment equations is derived by projection of the BTE onto a sequence of subspaces which represent elements of the phase space of the problem. Exact closure for the moment equations is provided by the Eddington tensor. A Petrov-Galerkin (PG) projection of the normalized BTE is formulated using a proper orthogonal decomposition (POD) basis representing the normalized radiation intensity over the whole phase space and time. The Eddington tensor linearly depends on the solution of the normalized BTE. By linear superposition of the POD basis functions, a low-rank expansion of the Eddington tensor is constructed with coefficients defined by the PG projected normalized BTE. The material energy balance (MEB) equation is coupled with the effective gray low-order equations which exist on the same dimensional scale as the MEB equation. The resulting TRT ROM is structure and asymptotic preserving. A detailed analysis of the ROM is performed on the classical Fleck-Cummings (F-C) TRT multigroup test problem in 2D geometry. Numerical results are presented to demonstrate the ROM's effectiveness in the simulation of radiation wave phenomena. Importantly, the ROM is shown to produce solutions with sufficiently high accuracy while using low-rank approximation of the normalized BTE solution. Essential physical characteristics of supersonic radiation wave are preserved in the ROM solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

In-Situ Observation of the Formation of Laser-Induced Periodic Surface Structures with Extreme Spatial and Temporal Resolution

Irradiation of solid surfaces with intense ultrashort laser pulses represents a unique way of depositing energy into materials. It allows to realize states of extreme electronic excitation and/or very high temperature and pressure and to drive materials close to and beyond fundamental stability limits. As a consequence, structural changes and phase transitions often occur along unusual pathways and under strongly nonequilibrium conditions. Furthermore, due to the inherent multiscale nature—both temporally and spatially—of these irreversible processes, their direct experimental observation requires techniques that combine high temporal resolution with the appropriate spatial resolution and the capability to obtain good quality data on a single pulse/event basis. In this respect, fourth-generation light sources, namely, short wavelength and short pulse free electron lasers (FELs), are offering new and fascinating possibilities. As an example, this chapter will discuss the results of scattering experiments carried out at the FLASH free electron laser at DESY (Hamburg, Germany), which allowed us to resolve laser-induced structure formation at surfaces on the nanometer to submicron length scale and in temporal regimes ranging from picoseconds to several nanoseconds with sub-picosecond resolution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Dynamic crushing of metal lattice metamaterials: Shock mode diagrams and transition to topology-independent compaction regime

Additively manufactured lattice metamaterials offer design versatility in strength and energy absorption and provide an additional degree of freedom through the selection of the lattice topology. Under quasistatic loading, the unit cell structure can strongly affect the stiffness, yield, and post-yield behavior, but whether and to what degree the effect of lattice topology persists into dynamic loading scenarios, up to the compaction shock regime, has not been established. LLNL’ s ALE3D hydrocode was used to perform a computational investigation of dynamic loading in multiple lattice types, including the gyroid, octet, Schwarz D, and rhombic dodecahedron, under impact velocities from 0.25 to 2.25 km/s. Shock Hugoniots for each lattice topology are generated and compared, suggesting that above a critical velocity, distinctions between architectures may not persevere and compacted lattices behave similarly. Here, to investigate the transition between topology-dependent quasistatic compression and the topology-independent regime above the critical velocity, a one-dimensional elastic-linear hardening plasticity-densified solid (E-LHP-DS) shock model for lattice materials was developed that relies upon confined compression to link the quasistatic and shock mechanics. Unlike similar works, the model does not assume rigid behavior prior to yield or locking behavior at densification, allowing a richer exploration of lattice mechanics. With only six parameters, the analytical model simultaneously fit quasistatic confined compression simulations for relative densities 0.1 $≤ \bar{ρ} ≤$ 0.9 and predicted dynamic compaction behavior to traverse several distinct shock modes, each defined by a critical impact speed (equivalently, critical stresses). Comparing the numerical results to the one-dimensional E-LHP-DS shock model predictions suggests that the topology-independence under strong shocks is linked to the onset of densification, which can be predicted based on quasistatic confined compression results.

Cellular material↗
DOE OSTI2025

Direct experimental constraints on the spatial extent of a neutrino wavepacket

Despite their high relative abundance in our Universe, neutrinos are the least understood fundamental particles of nature. In fact, the quantum properties of neutrinos emitted in experimentally relevant sources are theoretically contested and the spatial extent of the neutrino wavepacket is only loosely constrained by reactor neutrino oscillation data with a spread of 13 orders of magnitude. Here we present a method to directly access this quantity by precisely measuring the energy width of the recoil daughter nucleus emitted in the radioactive decay of beryllium-7. The final state in the decay process contains a recoiling lithium-7 nucleus, which is entangled with an electron neutrino at creation. The lithium-7 energy spectrum is measured to high precision by directly embedding beryllium-7 radioisotopes into a high-resolution superconducting tunnel junction that is operated as a cryogenic sensor. Under this approach, we set a lower limit on the Heisenberg spatial uncertainty of the recoil daughter of 6.2 pm, which implies that the final-state system is localized at a scale more than a thousand times larger than the nucleus itself. From this measurement, the first, to our knowledge, direct lower limit on the spatial extent of a neutrino wavepacket is extracted. These results may have implications in several areas including the theoretical understanding of neutrino properties, the nature of localization in weak nuclear decays and the interpretation of neutrino physics data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Black hole–neutron star mergers: The first mass gap and kilonovae

Observations of X-ray binaries indicate a dearth of compact objects in the mass range from ~2 –5 M ⊙ . The existence of this (first mass) gap has been used to discriminate between proposed engines behind core-collapse supernovae. From LIGO/Virgo observations of binary compact remnant masses, several candidate first mass gap objects, either neutron stars (NSs) or black holes (BHs), were identified during the O3 science run. Motivated by these new observations, we study the formation of BH-NS mergers in the framework of isolated classical binary evolution, using population synthesis methods to evolve large populations of binary stars (Population I and II) across cosmic time. We present results on the NS to BH mass ratios (q = M NS /M BH ) in merging systems, showing that although systems with a mass ratio as low as q = 0.02 can exist, typically BH-NS systems form with moderate mass ratios q = 0.1 –0.2. If we adopt a delayed supernova engine, we conclude that ~30% of BH-NS mergers may host at least one compact object in the first mass gap (FMG). Even allowing for uncertainties in the processes behind compact object formation, we expect the fraction of BH-NS systems ejecting mass during the merger to be small (from ~0.6 –9%). In our reference model, we assume: (i) the formation of compact objects within the FMG, (ii) natal NS/BH kicks decreased by fallback, (iii) low BH spins due to Tayler-Spruit angular momentum transport in massive stars. We find that ≲1% of BH-NS mergers will have any mass ejection and about the same percentage will produce kilonova bright enough to have a chance of being detected with a large (Subaru-class) 8 m telescope. Interestingly, all these mergers will have both a BH and an NS in the FMG.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Hydrodynamic computations of high-power laser drives generating metal ejecta jets from surface grooves

Understanding dynamic fragmentation in shock-loaded metals and predicting properties of the resulting ejecta are of considerable importance for both basic and applied science. The nature of material ejection has been shown to change drastically when the free surface melts on compression or release. In this work, we present hydrodynamic simulations of laser-driven microjetting from micron-scale grooves on a tin surface. We study microjet formation across a range of shock strengths from drives that leave the target solid after release to drives that induce shock melting in the target. The shock-state particle velocity (Up) varies from 0.3 to 3 km/s and the shock breakout pressure is 3–120 GPa. The microjet tip velocity is 1–8 km/s and the free-surface velocity varies from 0.1 to 5 km/s. Two tin equations of state are examined: a “soft” model (LEOS 501) where the target melts for U p > 1 km/s and a more detailed multiphase model (SESAME 2161) that melts for U p > 1.4 km/s. We use these two models to examine the influence of phase change and the choice of the material model on microjet formation and evolution. We observe in our computational results that jet formation can be classified into three regimes: a low-energy regime where material strength affects jet formation, a moderate-energy regime dominated by the changing phase of tin material, and a high-energy regime where results are insensitive to the material model and jet formation is described by an idealized steady-jet theory. Using an ensemble of 2D simulations, we show that these trends hold across a wide range of drive energies and groove angles.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Low mode implosion symmetry sensitivity in low gas-fill NIF cylindrical hohlraums

Achieving an efficient capsule implosion in National Ignition Facility indirect-drive target experiments requires symmetric hohlraum x-ray drive for the duration of the laser pulse. This is commonly achieved using two-sided two-cone laser irradiation of cylindrical hohlraums that, in principle, can zero the time average of all spherical harmonic asymmetry modes <6 as well as the time dependence of the usually dominant mode 2. In practice, experimental evidence indicates that maintaining symmetric drive becomes limited late in the pulse due to the inward expansion of the hohlraum wall and outward expansion of the capsule ablator plasmas impairing the propagation of the inner-cone laser beams. This effect is enhanced in hohlraums employing low gas-fill, now used almost exclusively as these provide the highest performing implosions and reduce Stimulated Brillouin and Raman backscatter losses, since the gas plasma provides less back pressure to limit blow-in of the hohlraum wall and capsule ablator plasmas. In order to understand this dynamic behavior, we combined multi-keV X-ray imaging of the wall and imploded fuel plasmas as we changed a single parameter at a time: hohlraum gas-fill, laser outer cone picket energy, radius of high density carbon capsules used, and laser beam polar and azimuthal pointing geometry. In this work, we developed a physics-based multi-parameter experimental scaling to explain the results that extend prior scalings and compare those to radiation hydrodynamic simulations to develop a more complete picture of how hohlraum, capsule, and laser parameters affect pole vs equator drive symmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Energy-angle distribution of the ions in the RF sheath of ICRH antennas

Radio-frequency sheaths forming at the surface of ion cyclotron resonance heating (ICRH) actuators in fusion experiments are associated with enhanced impurity sputtering from ICRH plasma-facing components (PFCs). The minimization of impurity fluxes from the ICRH PFC is a vital task for the usage of ICRH systems. Capturing the ion kinetics is important for an accurate description of plasma–material interactions, because the ion dynamics plays a crucial role in RF sheaths. Here, we present a hybrid particle-in-cell (hPIC) model able to capture the kinetic behavior of the ions. We analyze the kinetic ion energy-angle distributions (IEADs) impacting the RF antenna and its dependence on different plasma and RF sheath parameters. In particular, the IEAD dependence on RF frequency and magnetic field alignment is investigated. Using hPIC, we simulated a case emulating the latest experimental campaign from JET. The simulation showed that under specific plasma and RF parameters, the kinetic motion of ions results in a cusp formation in the IEAD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗