SEARCH · Engineering Papers

Results for “wave propagation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 1,009 records · Page 56

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Acoustic Waves From a Distant Explosion Recorded on a Continuously Ascending Balloon in the Middle Stratosphere

Abstract A helium‐filled mylar balloon carrying a smartphone and infrasound sensors ascended to a stratospheric height of 35 km over the surface detonation of a chemical explosive, with a total acoustic propagation distance of 127 km. The smartphone was configured to collect multi‐modal data at high rates from internal sensors. Analysis of the data shows successful collection of the explosion signal by both the smartphone's microphone and its accelerometers, the first from an ascending balloon. Comparison of the acoustic signal with that collected by other infrasound sensors, both airborne and ground‐based, provides insight into the possibilities and limitations of collecting acoustic data from the stratosphere.

47 OTHER INSTRUMENTATION↗

DOE OSTI2024

Direct solution of Minkowski-space Bethe-Salpeter equation in the massive Wick-Cutkosky model

Here in order to solve the Bethe-Salpeter equation (BSE) in the Minkowski space, we first introduce the Nakanishi integral representations of the Bethe-Salpeter amplitude and the Bethe-Salpeter wave function. We then derive the explicit integral equations for the corresponding spectral functions from the BSE for states of two scalar particles bound by a scalar-particle exchange interaction, where the propagators of constituents are allowed to be fully dressed. These integral equations are subsequently solved numerically in the variation of the Wick-Cutkosky model with massive exchange particles, where an algorithm of adaptive mash grid is proposed. The equations and algorithm we develop here serve as the foundation of Minkowski-space formulation of BSE for bound states of fermions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Post-detonation fireball modeling: Validation of freeze out approximations

A numerical simulation study was performed to examine the post-detonation reaction processes produced by the detonation of a 12 mm diameter hemispherical pentaerythritol tetranitrate (PETN) explosive charge. The simulations used a finite rate detailed chemical reaction model consisting of 59 species and 368 reactions to capture post-detonation reaction processes including air dissociation from Mach 19+ shock waves that initially break out of the PETN charge, reactions within the detonation products during expansion, and afterburning when the detonation products mix with the shock heated air. The multi-species and thermodynamically complete Becker-Kistiakowsky-Wilson real-gas equation of state is used for the gaseous phase to allow for the mixing of reactive species. A recent simplified reactive burn model is used to propagate the detonation through the charge and allow for detailed post-detonation reaction processes. The computed blast, shock structures, and mole fractions of species within the detonation products agree well with experimental measurements. A comparison of the simulation results to equilibrium calculations indicates that the assumption of a local equilibrium is fairly accurate until the detonation products rapidly cool to temperatures in the range of 1500-1900 K by expansion waves. Below this range, the computed results show mole fractions that are nearly chemically frozen within the detonation products for a significant portion of expansion. In conclusion, these results are consistent with the freeze out approximation used in the blast modeling community.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Machine learning enhanced predictions of ICRF heating: Overcoming numerical limitations via data curation

In this work, we present the development of robust surrogate models for Ion Cyclotron Range of Frequencies (ICRF) and High-Harmonic Fast Wave (HHFW) heating predictions in fusion plasmas. Building upon our previous efforts to achieve real-time capable models, we identify the cause of the outliers found using TORIC in certain HHFW heating scenarios. The outliers are observed to be spurious ion Bernstein wave (IBW)-like modes caused by a wavelength control algorithm designed to address challenging scenarios with high perpendicular wavenumbers. The effect arises from the modulation in the perpendicular susceptibility, which can induce sign reversal and IBW-like propagation for scenarios featuring normalized ion Larmor radius λ i ≫ 1. We use TORIC with this algorithm disabled to generate a novel HHFW-NSTX database that is free of outliers. Surrogate models trained on this database, including Random Forest Regressor (RFR), Multi-Layer Perceptrons, and Gaussian Process Regressors (GPR), demonstrate the ability to accurately predict HHFW heating profiles, with regression scores of R 2 ∈[0.93−0.99]. Additionally we demonstrate that it is possible to generalize predictions beyond training data by the use of both RFR and GPR models, enabling the prediction of scenarios previously limited to the original model. GPR models also provide uncertainty quantification, offering insights into model confidence. This work introduces a comprehensive Verification, Validation, and Uncertainty Quantification methodology for surrogate modeling, applicable not only to ICRF heating but also to other RF heating challenges and fusion physics problems. Beyond accelerated inference, these models show effective extrapolation capabilities, providing an alternative for addressing numerical challenges.

Artificial neural networks↗

DOE OSTI2021

Generation of isolated soft x-ray attosecond pulses with mid-infrared driving lasers via transient phase-matching gating

Abstract Isolated attosecond pulses (IAPs) in the soft x-ray (SXR) region are highly desirable for attosecond time-resolved experiments. Here we identify a transient phase matching gating method for the generation of such IAPs with mid-infrared (MIR) lasers. This gating method works when a loosely focused, long-duration MIR Gaussian driving beam is slightly reshaped during its propagation in the gas medium under the ‘critical’ ionization condition. Quantitatively, the calculated coherence length of high harmonic is used to analyze the mechanism of the gating method, by using one-dimensional plane-wave beams and by comparing 2000 nm and 800 nm lasers. The robustness of the generation method is checked by varying laser’s carrier-envelope-phase and gas pressure. This gating method provides with an alternative route to efficiently produce tabletop ultrashort attosecond SXR light sources with the emerging MIR lasers.

Physics↗

DOE OSTI2023

Emergent quantum mechanics at the boundary of a local classical lattice model

We formulate a model in which quantum mechanics emerges from classical mechanics. Given a local Hamiltonian H acting on n qubits, we define a local classical model with an additional spatial dimension whose boundary dynamics is approximately—but to arbitrary precision—described by Schrödinger's equation and H. The bulk consists of a lattice of classical bits that propagate towards the boundary through a circuit of stochastic matrices. The bits reaching the boundary are governed by a probability distribution whose deviation from the uniform distribution can be interpreted as the quantum-mechanical wave function. Bell nonlocality is achieved because information can move through the bulk much faster than the boundary speed of light. Finally, we analytically estimate how much the model deviates from quantum mechanics, and we validate these estimates using computer simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Rupture Model of the 5 April 2024 Tewksbury, New Jersey, Earthquake Based on Regional Lg -Wave Data

On 5 April 2024, an earthquake of magnitude 4.8 occurred in Tewksbury, New Jersey. It was the largest instrumentally recorded event since 1900 in New Jersey and southern New York. Millions of people around New York City, ~65 km east–northeast of Tewksbury, felt the shaking from the mainshock, but the epicentral area experienced no known significant property damages. We determine the focal mechanism, which is oblique faulting, and retrieve the Lg–wave relative source time functions (RSTFs) from the stations at regional distances to understand rupture processes and ground motions. Our fault–slip models well explain azimuthal variations of the RSTFs. The models show the rupture propagating toward the east–northeast (~50° to 60°), not along the fault strike. The slip distribution on the nodal plane striking north and dipping to the east shows a slip area of 1.1 km radius with the rupture propagating down–dip. The down–dip rupture may account for the observed lack of strong shaking in the epicentral area.

58 GEOSCIENCES↗

DOE OSTI2022

An improved framework for the dynamic likelihood filtering approach to data assimilation

Here, we propose improvements to the Dynamic Likelihood Filter (DLF), a Bayesian data assimilation filtering approach, specifically tailored to wave problems. The DLF approach was developed to address the common challenge in the application of data assimilation to hyperbolic problems in the geosciences and in engineering, where observation systems are sparse in space and time. When these observations have low uncertainties, as compared to model uncertainties, the DLF exploits the inherent nature of information and uncertainties to propagate along characteristics to produce estimates that are phase aware as well as amplitude aware, as would be the case in the traditional data assimilation approach. Along characteristics, the stochastic partial differential equations underlying the linear or nonlinear stochastic dynamics are differential equations. This study focuses on developing the explicit challenges of relating dynamics and uncertainties in the Eulerian and Lagrangian frames via dynamic Gaussian processes. It also implements the approach using the ensemble Kalman filter (EnKF) and compares the DLF approach to the conventional one with respect to wave amplitude and phase estimates in linear and nonlinear wave problems. Numerical comparisons show that the DLF/EnKF outperforms the EnKF estimates, when applied to linear and nonlinear wave problems. This advantage is particularly noticeable when sparse, low uncertainty observations are used.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Altering Terahertz Sound Propagation in a Liquid upon Nanoparticle Immersion

One of the grand challenges of new generation Condensed Matter physicists is the development of novel devices enabling the control of sound propagation at terahertz frequency. Indeed, phonon excitations in this frequency window are the leading conveyor of heat transfer in insulators. Their manipulation is thus critical to implementing heat management based on the structural design. To explore the possibility of controlling the damping of sound waves, we used high spectral contrast Inelastic X-ray Scattering (IXS) to comparatively study terahertz acoustic damping in a dilute suspension of 50 nm nanospheres in glycerol and on pure glycerol. Bayesian inference-based modeling of measured spectra indicates that, at sufficiently large distances, the spectral contribution of collective modes in the glycerol suspension becomes barely detectable due to the enhanced damping, the weakening, and the slight softening of the dominant acoustic mode.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Comparing theoretical predictions of radiation-free velocities of edge dislocations to molecular dynamics simulations

Transonic defect motion is of interest for high strain-rate plastic deformation as well as for crack propagation. Ever since Eshelby's 1949 prediction [J. D. Eshelby, Proc. Phys. Soc. A 62, 307 (1949)] in the isotropic limit of a “radiation-free” transonic velocity $v$ RF = $\sqrt{2}$c T , where shock waves are absent, there has been speculation about the significance of radiation-free velocities (if they truly exist) for defect mobility. Here, in this work, we argue that they do not play any significant role in dislocation dynamics in metals, based on comparing theoretical predictions of radiation-free velocities for transonic edge dislocations with molecular dynamics simulations for two face-centered cubic metals: Ag, where theory predicts radiation-free states, and Cu, where it does not.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Nonlinear wavepacket dynamics in proximity to a stationary inflection point

A stationary inflection point (SIP) in the Bloch dispersion relation of a periodic waveguide is an exceptional point degeneracy where three Bloch eigenmodes coalesce, forming the so-called frozen mode with a divergent amplitude and vanishing group velocity of its propagating component. For this work, we have developed a theoretical framework to study the time evolution of wavepackets centered at an SIP. Analysis of the evolution of the statistical moment distribution of linear pulses shows a strong deviation from the conventional ballistic wavepacket dynamics in dispersive media. The presence of nonlinear interactions dramatically changes the situation, resulting in a mostly ballistic propagation of nonlinear wavepackets with the speed and even the direction of propagation essentially dependent on the wavepacket amplitude. Such a behavior is unique to nonlinear wavepackets centered at an SIP and can be used for the realization of a novel family of beam power routers for classical waves.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

Evaluation of geophysical and anthropogenic sources of hydroacoustic noise in the Alaskan arctic

Quantifying the ocean soundscape is crucial for ocean-based seismoacoustic monitoring; it sets a baseline for the kinds and sizes of signals that can be detected above the background noise. As sea ice recedes, human activity in and around the Arctic Ocean is increasing, elevating sound levels and heightening the urgency of monitoring. Here, seven years of passive acoustic recordings from a National Oceanographic and Atmospheric Administration hydrophone in the Beaufort Sea are analyzed, focusing on frequencies between 1 and 125 Hz, a band of interest for detecting regional earthquakes and similar events. Sound is related to geophysical and anthropogenic sources, and seasonal and intraseasonal variations in the soundscape are examined. Sea ice emerges as a keystone feature of the Arctic Ocean acoustic environment, controlling or contributing to ambient sound levels at all frequencies studied. During low-ice months, sound levels are dominated by wind and sea surface waves, and in ice-covered months, wind-driven ice noise dominates. Seismic air gun surveys are prominent during low-ice periods, with sound levels decreasing with increasing distance and bathymetric complexity along the propagation path. The implications of this baseline soundscape for event detection in the Alaskan Arctic are discussed.

Niklasson, Siobhan [Los Alamos National Laboratory↗

DOE OSTI2022

Turbulent Regimes in Collisions of 3D Alfvén-wave Packets

Using three-dimensional gyrofluid simulations, we revisit the problem of Alfvén-wave (AW) collisions as building blocks of the Alfvénic turbulent cascade and their interplay with magnetic reconnection at magnetohydrodynamic (MHD) scales. Depending on the large-scale value of the nonlinearity parameter χ 0 (the ratio between the AW linear propagation time and nonlinear turnover time), different regimes are observed. For strong nonlinearities (χ 0 ~ 1), turbulence is consistent with a dynamically aligned, critically balanced cascade—fluctuations exhibit a scale-dependent alignment $\sin {\theta }_{{k}_{\perp }}\propto {k}_{\perp }^{-1/4}$, resulting in a ${k}_{\perp }^{-3/2}$ spectrum and ${k}_{\parallel }\propto {k}_{\perp }^{1/2}$ spectral anisotropy. At weaker nonlinearities (small χ 0 ), a spectral break marking the transition between a large-scale weak regime and a small-scale ${k}_{\perp }^{-11/5}$ tearing-mediated range emerges, implying that dynamic alignment occurs also for weak nonlinearities. At χ 0 < 1 the alignment angle ${\theta }_{{k}_{\perp }}$ shows a stronger scale dependence than in the χ 0 ~ 1 regime, namely $\sin {\theta }_{{k}_{\perp }}\propto {k}_{\perp }^{-1/2}$ at χ 0 ~ 0.5, and $\sin {\theta }_{{k}_{\perp }}\propto {k}_{\perp }^{-1}$ at χ 0 ~ 0.1. Dynamic alignment in the weak regime also modifies the large-scale spectrum, scaling approximately as ${k}_{\perp }^{-3/2}$ for χ 0 ~ 0.5 and as ${k}_{\perp }^{-1}$ for χ 0 ~ 0.1. A phenomenological theory of dynamically aligned turbulence at weak nonlinearities that can explain these spectra and the transition to the tearing-mediated regime is provided; at small χ 0 , the strong scale dependence of the alignment angle combines with the increased lifetime of turbulent eddies to allow tearing to onset and mediate the cascade at scales that can be larger than those predicted for a critically balanced cascade by several orders of magnitude. Such a transition to tearing-mediated turbulence may even supplant the usual weak-to-strong transition.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Dynamic Rupture Simulations of Caldera Collapse Earthquakes: Effects of Wave Radiation, Magma Viscosity, and Evidence of Complex Nucleation at Kilauea 2018

All instrumented basaltic caldera collapses have generated M w > 5 very long period earthquakes. However, previous studies of source dynamics have been limited to lumped models treating the caldera block as rigid, leaving open questions related to how ruptures initiate and propagate around the ring fault, and the seismic expressions of those dynamics. Here, we present the first 3D numerical model capturing the nucleation and propagation of ring fault rupture, the mechanical coupling to the underlying viscoelastic magma, and the associated seismic wavefield. We demonstrate that seismic radiation, neglected in previous models, acts as a damping mechanism reducing coseismic slip by up to half, with effects most pronounced for large magma chamber volume/ring fault radius or highly compliant crust/compressible magma. Viscosity of basaltic magma has negligible effect on collapse dynamics. In contrast, viscosity of silicic magma significantly reduces ring fault slip. We use the model to simulate the 2018 Kilauea caldera collapse. Three stages of collapse, characterized by ring fault rupture initiation and propagation, deceleration of the downward-moving caldera block and magma column, and post-collapse resonant oscillations, in addition to chamber pressurization, are identified in simulated and observed (unfiltered) near-field seismograms. A detailed comparison of simulated and observed displacement waveforms corresponding to collapse earthquakes with hypocenters at various azimuths of the ring fault reveals a complex nucleation phase for earthquakes initiated on the northwest. Our numerical simulation framework will enhance future efforts to reconcile seismic and geodetic observations of caldera collapse with conceptual models of ring fault and magma chamber dynamics.

58 GEOSCIENCES↗

DOE OSTI2020

Theory of the Tertiary Instability and the Dimits Shift from Reduced Drift-Wave Models

Tertiary modes in electrostatic drift-wave turbulence are localized near extrema of the zonal velocity U(x) with respect to the radial coordinate x. We argue that these modes can be described as quantum harmonic oscillators with complex frequencies, so their spectrum can be readily calculated. The corresponding growth rate $γ$ TI is derived within the modified Hasegawa-Wakatani model. We show that $γ$TI equals the primary-instability growth rate plus a term that depends on the local U''; hence, the instability threshold is shifted compared to that in homogeneous turbulence. This provides a generic explanation of the well-known yet elusive Dimits shift, which we find explicitly in the Terry-Horton limit. Linearly unstable tertiary modes either saturate due to the evolution of the zonal density or generate radially propagating structures when the shear |U'| is sufficiently weakened by viscosity. The Dimits regime ends when such structures are generated continuously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Stopping-power enhancement from discrete particle-wake correlations in high-energy-density plasmas

Three-dimensional (3D) simulations of electron beams propagating in high-energy-density plasmas using the quasistatic Particle-in-Cell (PIC) code QuickPIC demonstrate a significant increase in stopping power when beam electrons mutually interact via their wakes. Each beam electron excites a plasma wave wake of wavelength ~ 2 π c / ω p e , where c is the speed of light and ω p e is the background plasma frequency. Here we show that a discrete collection of electrons undergoes a beam-plasma-like instability caused by mutual particle-wake interactions that causes electrons to bunch in the beam, even for beam densities n b for which fluid theory breaks down. This bunching enhances the beam's stopping power, which we call “correlated stopping,” and the effect increases with the “correlation number” N b ≡ n b ( c / ω p e ) 3 . For example, a beam of monoenergetic 9.7 MeV electrons with N b = 1 / 8 , in a cold background plasma with n e = 10 26 cm – 3 (450 g cm – 3 DT), has a stopping power of 2.28 ± 0.04 times the single-electron value, which increases to 1220 ± 5 for N b = 64 . The beam also experiences transverse filamentation, which eventually limits the stopping enhancement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Analyzing LF/VLF Lightning Waveforms to Estimate D-region Electron Density Profiles

Lightning waveforms in the low frequency (LF; 30-300 kHz) and the very low frequency (VLF; 3-30 kHz) bands can be exploited to produce data-driven ionospheric D-region electron density profile (EDP) estimates with significantly higher spatial and temporal coverage than previously available. The lightning waveforms used in this paper are signals detected in the LF/VLF of negative cloud-to-ground lightning by the Earth Networks Total Lightning Detection Network. Each waveform contains a ground wave and a time-delayed ionospheric reflection. The time delay between the ground wave and ionospheric reflection has previously been used to estimate a single specular reflection altitude, where LF/VLF emissions are reflected by the ionosphere. Here, we expand upon previous methods to include filtering and spectral analysis, and account for oblique propagation to produce higher-order estimates for reflection altitudes and corresponding electron densities. Once estimated, reflection altitudes and corresponding electron densities can be used to derive parameters β and h’, which define an EDP for the D-region. In this study, the lightning waveform (LW) analysis is demonstrated using a single representative 24-hour dataset over the Southeast United States, and then extended to a total of 10 separate datasets with varying locations and ionospheric conditions. The LW-derived D-region EDPs are in agreement with predictions made by the Faraday International Reference Ionosphere model, and the LW EDPs β and h’ values are consistent with previous LF/VLF-derived estimates.

D-region ionosphere↗

DOE OSTI2021

Multipoint Observations of Quasiperiodic Emission Intensification and Effects on Energetic Electron Precipitation

The two Van Allen Probes simultaneously recorded a coherently modulated quasiperiodic (QP) emission that persisted for 3 h. The magnetic field pulsation at the locations of the two satellites showed a substantial difference, and their frequencies were close to but did not exactly match the repetition frequency of QP emissions for most of the time, suggesting that those coherent QP emissions probably originated from a common source, which then propagated over a broad area in the magnetosphere. The QP emissions were amplified by local anisotropic electron distributions, and their large-scale amplitudes were modulated by the plasma density. A novel observation of this event is that chorus waves at frequencies above QP emissions exhibit a strong correlation with QP emissions. Those chorus waves intensified when the QP emissions reach their peak frequency. This indicates that embryonic QP emissions may be critical for its own intensification as well as chorus waves under certain circumstances. A Low-Earth-Orbit POES satellite observed enhanced energetic electron precipitation in conjunction with the Van Allen Probes, providing direct evidence that QP emissions precipitate energetic electrons into the atmosphere. Finally, this scenario is quantitatively confirmed by our quasilinear diffusion simulation results.

79 ASTRONOMY AND ASTROPHYSICS↗