Engineering Papers⌕ Search

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.

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

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Using Dense Array Waveform Correlations to Build a Velocity Model with Stochastic Variability

In an ongoing effort to improve 3D seismic-wave propagation modeling for frequencies up to 10 Hz, we used cross correlations between vertical-component waveforms from an underground chemical explosion to estimate the statistical properties of small-scale velocity heterogeneities. The waveforms were recorded by a dense 2D seismic array deployed during the Source Physics Experiments for event number 5 (SPE-5) in a series of six underground chemical explosions, conducted at the Nevada National Security Site. The array consisted of 996 geophones with a 50–100 m grid spacing, deployed at the SPE site at the north end of the Yucca Flat basin. The SPE were conducted to investigate the generation and propagation of seismic and acoustic waves from underground explosions. Comparisons of decay rates of waveform cross correlations as function of interstation distance, computed for observed and synthetic seismograms from the SPE-5 chemical explosion, were used to constrain statistical properties of correlated stochastic velocity perturbations representing small-scale heterogeneities added to a geology-based velocity model of the Yucca Flat basin. Using comparisons between recorded and simulated waveform cross correlations, we were able to recover sets of statistical properties of small-scale velocity perturbations in the velocity model that produce the best-fit between the recorded and simulated ground motion. The stochastic velocity fluctuations in the velocity model that produced the smallest misfits have a horizontal correlation distance of between 400 and 800 m, a vertical correlation distance between 100 and 200 m, and a standard deviation of 10% from the nominal model velocity in the alluvium basin layers. They also have a horizontal correlation distance of 1000 m, a vertical correlation distance of 250 m, and a standard deviation of 6% in the underlying and consolidated sedimentary layers, up to a depth of 4 km. Comparisons between observed and simulated wavefields were used to assess the proposed small-scale heterogeneity enhancements to the Yucca Flat basin model. We found that adding a depth-resolved stochastic variability to the geology-based velocity model improves the overall performance of ground-motion simulations of an SPE-5 explosion in the modeled frequency range up to 10 Hz. The results may be applicable to other similar basins.

58 GEOSCIENCES↗
DOE OSTI2020

Study on Intelligent Compaction-Equipment Logistics Scheduling and Propagation Characteristics of Vibration Wave in Nonlinear Systems with Multistability Based on Field Test

The equipment scheduling and propagation characteristics of vibration wave from vibratory roller ⟶ filling material nonlinear systems with multistability are the core problems of subgrade intelligent construction technology, and the logistics scheduling of the equipment is directly related to the construction efficiency. Aiming at the shortages, one typical subgrade located at the Gu’an station of Beijing-Xiong’an city railway is selected to research and finish the field tests; some findings are shown as follows: first, some valuable suggestions about the logistics scheduling of intelligent equipment are proposed, which can break the barriers between the organizations and improve construction efficiency; second, when the vibration wave propagates from the vibratory roller ⟶ surface of filling material ⟶ different buried depths of filling material, the peak acceleration of vibration wave gradually decreases and is hyperbolic distribution approximately. At the same time, the sensitive of attenuation is shown as follows: Z R 2 >0.9), and the concrete functional relationship among different amplitudes of harmonics can be summarized as y = Ae −BX ; fourth, the vibration energy is mainly concentrated near 10–30 Hz in the vibratory roller, but when the vibration wave propagates from vibratory roller⟶filling material, the vibration energy gradually decreases with the increase in depth, and the marginal spectrum gradually changes from one peak to two peaks, that is, 30–50 Hz and 50–100 Hz; fifth, the vibration energy in the vibrational wheel is distributed averagely in the compaction process, and the effective compaction time is two seconds, which will be helpful for revealing the propagation characteristics of vibration wave, optimizing the compaction quality control models and providing some support for the development of intelligent compaction theory of railway subgrade.

Han, Yixuan↗
DOE OSTI2022

A high accuracy/resolution spectral element/Fourier–Galerkin method for the simulation of shoaling non-linear internal waves and turbulence in long domains with variable bathymetry

A high-order hybrid continuous-Galerkin numerical method, designed for the simulation of non-linear, non -hydrostatic internal waves and turbulence in long computational domains with complex bathymetry, is presented. The spatial discretization in the non-periodic wave-propagating directions, utilizes the nodal spectral element method. Such a high-order element-based discretization allows the highly accurate representation of complex domain geometry along with the flexibility of concentrating resolution in areas of interest. Under the assumption of the normal-to-isobath propagation of non-linear internal waves, a third periodic direction is incorporated via a Fourier-Galerkin discretization. The distinct non-hydrostatic nature of non-linear internal waves and, any instabilities and turbulence therein, necessitates the numerically challenging solution of the pressure Poisson problem. A defining feature of this work is the application of a domain decomposition approach, combined with block-Jacobi/deflation-based preconditioning to the pressure Poisson problem. Such a combined approach is particularly suitable for the long high aspect-ratio complex domains of interest and enables the efficient high-accuracy reproduction of the non-hydrostatic dynamics of non-linear internal waves. Implementation details are also described in the context of the stability of the solver and its parallelization strategy. A series of benchmarks of increasing complexity demonstrate the robustness of the flow solver. The benchmarks culminate with the three-dimensional simulation of a convectively breaking mode-one non-linear internal wave over a realistic South-China-Sea bathymetric transect and background current/stratification profiles.

Deflation↗
DOE OSTI2022

Effects of polarization-reversed EMIC waves on the ring current dynamics

Electromagnetic ion cyclotron (EMIC) waves are widely believed to play an important role in influencing the radiation belt and ring current dynamics. Most studies have investigated the effects or characteristics of EMIC waves by assuming their left-handed polarization. However, recent studies have found that the reversal of polarization, which occurs at higher latitudes along the wave propagation path, can change the wave-induced pitch angle diffusion coefficients. Whether such a polarization reversal can influence the global ring current dynamics remains unknown. In this study, we investigate the ring current dynamics and proton precipitation loss in association with polarization-reversed EMIC waves by using the ring current–atmosphere interactions model (RAM). The results indicate that the polarization reversal of H-band EMIC waves can truly decrease the scattering rates of protons of 10 to 50 keV or >100 keV in comparison with the scenario in which the EMIC waves are considered purely left-handed polarized. Additionally, the global ring current intensity and proton precipitation may be slightly affected by the polarization reversal, especially during prestorm time and the recovery phase, but the effects are not large during the main phase. This is probably because the H-band EMIC waves contribute to the proton scattering loss primarily at E < 10 keV, an energy range that is not strongly affected by the polarization reversal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Visualization of Thermal and Strain Dynamics in X-ray Optics by Time-Resolved Rocking Curve Imaging

The absorption of intense X-ray pulses from free-electron lasers by X-ray crystal optics triggers rapid thermal and mechanical responses, including lattice compression, expan- sion, and strain-wave propagation, which dynamically modify diffraction conditions and wavefront properties. These effects pose significant challenges for next-generation high- repetition-rate X-ray sources, particularly for crystal monochromators used in XFEL beam- lines and self-seeding. Here, we demonstrate time-resolved rocking curve imaging (TrRCI) of X-ray Bragg optics with high spatial resolution (6.4 μm) and sub-1-μrad angular resolu- tion, enabling direct visualization of lattice dynamics from nanoseconds to microseconds. By combining a high-flux beamline at SPring-8, a scintillator-coupled sCMOS detector, and a synchronized femtosecond laser system, we capture lattice dynamics induced by tran- sient laser heating with a strain sensitivity on the order of 10−6 . Our results reveal transient lattice deformation and the propagation of laser-induced strain waves over millimeter-scale distances. This approach provides experimental benchmarks for understanding heat-load effects in high-repetition-rate X-ray optics.

Source record↗
DOE OSTI2023

Ion cyclotron heating at high plasma density in Proto-MPEX

The physics of ion cyclotron heating (ICH) relevant to the steady-state linear machine MPEX (Material Plasma Exposure eXperiment) has been explored in its predecessor, short-pulse device: Proto-MPEX. MPEX will utilize fundamental ICH to increase heat flux at the target and produce ion temperatures and velocity distributions with improved fidelity to those found in a tokamak divertor region, in comparison to those produced by substrate biasing. In the experiments on Proto-MPEX described here, bulk ion temperatures up to ~15 eV have been achieved with 20 kW net ICH power at 6.5 MHz, using ICH heating of a deuterium plasma produced by a helicon plasma source. The heat flux at the target has been observed to increase throughout the plasma cross section, including in the core region. Core $T_i$ and target heat flux are observed to scale linearly with injected ICH power. Further, measurements of plasma loading and target heat flux as a function of the magnetic field strength at the antenna, together with modeling of the wave propagation from the antenna to the ion cyclotron resonance using the ANTENA and COMSOL codes with a warm plasma dielectric tensor, indicate that power is coupled to the core plasma via fast wave excitation of a kinetic Alfvén wave.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Validation of seismic velocity models in southern California with full-waveform simulations

SUMMARY Crustal seismic velocity models provide essential information for many applications including earthquake source properties, simulations of ground motion and related derivative products. We present a systematic workflow for assessing the accuracy of velocity models with full-waveform simulations. The framework is applied to four regional seismic velocity models for southern California: CVM-H15.11, CVM-S4.26, CVM-S4.26.M01 that includes a shallow geotechnical layer, and the model of Berg et al. For each model, we perform 3-D viscoelastic wave propagation simulations for 48 virtual seismic noise sources (down to 2 s) and 44 moderate-magnitude earthquakes (down to 2 s generally and 0.5 s for some cases) assuming a minimum shear wave velocity of 200 m s–1. The synthetic waveforms are compared with observations associated with both earthquake records and noise cross-correlation data sets. We measure, at multiple period bands for well-isolated seismic phases, traveltime delays and normalized zero-lag cross-correlation coefficients between the synthetic and observed data. The obtained measurements are summarized using the mean absolute derivation of time delay and the mean correlation coefficient. These two metrics provide reliable statistical representations of model quality with consistent results in all data sets. In addition to assessing the overall (average) performance of different models in the entire study area, we examine spatial variations of the models’ quality. All examined models show good phase and waveform agreements for surface waves at periods longer than 5 s, and discrepancies at shorter periods reflecting small-scale heterogeneities and near-surface structures. The model performing best overall is CVM-S4.26.M01. The largest misfits for both body and surface waves are in basin structures and around large fault zones. Inaccuracies generated in these areas may affect tomography and model simulation results at other regions. The seismic velocity models for southern California can be improved by adding better resolved structural representations of the shallow crust and volumes around the main faults.

58 GEOSCIENCES↗
DOE OSTI2022

Visualization and quantum control of light-accelerated condensates by terahertz multi-dimensional coherent spectroscopy

Abstract Characterizing and controlling high-order correlation of quantum systems is key for developing quantum devices and switching technologies. Although conventional static and ultrafast spectroscopy gives access to collective excitations characterizing quantum states, more exotic correlations cannot be easily separated from other contributions. Here we develop density matrix simulations to show that seventh-order-wave-mixing peaks with distinct temperature and field dependences in two-dimensional terahertz nonlinear spectra reveal light-induced correlations in non-equilibrium superconducting states. Above critical terahertz driving, these emerging peaks split from conventional peaks along the second axis introduced by pump-probe relative phase in two-dimensional frequency space. They are photo-generated by correlations between two-photon fluctuations and interacting quasi-particle and quasi-particle/Higgs superconductor excitations. By photo-inducing persistent symmetry breaking via light-wave propagation, we also demonstrate seventh-order-wave-mixing sensing of Higgs collective modes. Our theory suggests to use multi-dimensional spectroscopy for quantum sensing of light-driven superconductivity and paves a path for quantum operations by few-cycle-THz-periodic photocurrent modulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A description of structured waves in shock compressed particulate composites

Dynamic compression of composite materials is of scientific interest because the mechanical mismatch between internal phases challenges continuum theories. Typical assumptions about steady wave propagation and quasi-instantaneous state changes require reexamination along with the need for time-dependent models. To that end, data and models are presented here for the shock compression of an idealized particulate composite. To serve as a generic representative of this material class, a polymer matrix was filled with tungsten particles, ranging from 1 to 50 vol. %. This creates a simple microstructure containing randomly scattered particles with an extreme impedance mismatch to the binding matrix. These materials were parallel plate impact loaded by Al flyers traveling at 1.8–5.0 km/s. Velocimetry provided records of the equilibrium state and the compression wave structure for each case with trends quantified by an empirical fit. The same quantities were also studied as a function of the wave's propagation distance. A homogenized viscoelastic model then made it possible to progress from cataloging trends to making predictions. Starting from a Mie–Gruneisen equation of state, additional time varying terms were added to capture the transient response. After calibration, accurate predictions of the steady wave structure were possible.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Unitary Qubit Lattice Algorithms for Plasma Physics

This final technical report summarizes research conducted under DOE Award DE-SC0021653 to develop unitary Quantum Lattice Algorithms for modeling electromagnetic wave propagation and scattering in complex media, including plasmas. The project developed and validated quantum-inspired formulations of Maxwell's equations that preserve unitary evolution and can be evaluated on classical high-performance computing systems while providing a foundation for future quantum-computing implementations. Major accomplishments include the development of two- and three-dimensional algorithms for electromagnetic scattering; scalable, distributed-memory implementations demonstrated on the Perlmutter supercomputer; formulations for nonlinear lossless fluid dynamics and cold, lossless, inhomogeneous magnetized plasmas; and an explicit quantum algorithm for a time-discretized Lorenz model. Simulations reproduced a range of characteristic wave phenomena, including transient effects that are not readily apparent in conventional frequency-domain studies, demonstrating the effectiveness of the proposed approach for modeling complex electromagnetic and plasma systems. The work establishes a unified theoretical and computational framework for quantum and quantum-inspired simulation and provides a foundation for future implementation on fault-tolerant quantum systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Dependence of the kinetic energy absorption capacity of bistable mechanical metamaterials on impactor mass and velocity

Using an alternative mechanism to dissipation or scattering, bistable structures and mechanical metamaterials have shown promise for mitigating the detrimental effects of impact by reversibly locking energy into strained material. Herein, we extend prior works on impact absorption via bistable metamaterials to computationally explore the dependence of kinetic energy transmission on the velocity and mass of the impactor, with strain rates exceeding 10 2 s –1 . We observe a large dependence on both impactor parameters, ranging from significantly better to worse performance than a comparative linear material. We then correlate the variability in performance to solitary wave formation in the system and give analytical estimates of idealized energy absorption capacity under dynamic loading. In addition, we find a significant dependence on damping accompanied by a qualitative difference in solitary wave propagation within the system. The complex dynamics revealed in this study offer potential future guidance for the application of bistable metamaterials to applications including human and engineered system shock and impact protection devices.

42 ENGINEERING↗
DOE OSTI2021

A generalized machine learning workflow to visualize mechanical discontinuity

Accurate detection and mapping of mechanical discontinuity in materials has widespread industrial and research applications. Herein, we developed a generalized machine-learning framework for visualizing single mechanical discontinuity embedded in material of any composition, velocity, density, porosity, and size with limited data. The proposed visualization of discontinuity requires accurate estimations of the length, location, and orientation of the embedded discontinuity by processing multipoint wave-transmission measurements. k-Wave simulator is used to create a large dataset of elastic waveforms recorded during multi-point wave-transmission measurements through materials containing single mechanical discontinuity. k-Wave simulator considers the wave attenuation, dispersion, and mode conversion in wave motion. Discrete wavelet transform (DWT) and statistical feature extraction are essential for data preprocessing prior to the data-driven model development. DWT also minimizes the effect of noise. Using hyper-parameter tuning and cross validation, gradient boosting regression can visualize the mechanical discontinuity with an accuracy of 0.85, in terms of coefficient of determination. A double-layered neural network-based regression has better performance with an accuracy of 0.95. Use of convolutional neural network converts the predictive task from a waveform processing to an image processing problem. Convolutional neural network achieved a generalization performance of 0.91. The proposed generalized workflow requires robust simulation of wave propagation, signal processing, feature engineering, and model evaluation. Sensors closest to the source and those located opposite the source are the most significant for the desired visualization. Notably, the sensors closest to the source capture the non-linear associations, whereas the sensor on the border opposite to the source capture the linear associations between the measured waveforms and the properties of the mechanical discontinuity.

42 ENGINEERING↗
DOE OSTI2021

Transverse impact by RCCs on S-glass and Kevlar® FRC strips

This study aims to isolate interactions between plies and tows and reveal fundamental physics involved in transverse impact on fiber-reinforced composite (FRC) structures. Composite strips were sectioned from large panels and characterized by optical photography, three-dimensional synchrotron X-ray computed tomography, and scanning electron microscopy (SEM). Each strip was impacted perpendicularly by a right circular cylinder (RCC) projectile at a velocity ranging from ~ 150 to 600 m/s. The global strip behavior, as well as localized deformation and failure of the strip near the projectile corner, were both captured by high-speed optical imaging. S-Glass FRC strips were observed to fail in tension ahead of the RCC projectiles’ flat surfaces while Kevlar® FRC strips fractured at the projectile corners. The concept of critical velocity region previously used for impact on yarns was introduced to define different failure modes of each composite strip type. The strip damage extent was found to increase with the impact velocity and reach the maximum at the upper limit of the critical velocity region. Above the critical velocity region, the damage extent decreased with impact velocity. Wave propagations and load histories in the composite strips during impact were quantified and compared with Smith’s theory. Finally, critical velocities of single fibers, yarns, and composite strips and ballistic limits of single-ply and multi-ply composite panels were compared to provide insight into the design of impact-resistant fabrics and composites.

42 ENGINEERING↗
DOE OSTI2020

Global adjoint tomography—model GLAD-M25

SUMMARY Building on global adjoint tomography model GLAD-M15, we present transversely isotropic global model GLAD-M25, which is the result of 10 quasi-Newton tomographic iterations with an earthquake database consisting of 1480 events in the magnitude range 5.5 ≤ Mw ≤ 7.2, an almost sixfold increase over the first-generation model. We calculated fully 3-D synthetic seismograms with a shortest period of 17 s based on a GPU-accelerated spectral-element wave propagation solver which accommodates effects due to 3-D anelastic crust and mantle structure, topography and bathymetry, the ocean load, ellipticity, rotation and self-gravitation. We used an adjoint-state method to calculate Fréchet derivatives in 3-D anelastic Earth models facilitated by a parsimonious storage algorithm. The simulations were performed on the Cray XK7 ‘Titan’ and the IBM Power 9 ‘Summit’ at the Oak Ridge Leadership Computing Facility. We quantitatively evaluated GLAD-M25 by assessing misfit reductions and traveltime anomaly histograms in 12 measurement categories. We performed similar assessments for a held-out data set consisting of 360 earthquakes, with results comparable to the actual inversion. We highlight the new model for a variety of plumes and subduction zones.

58 GEOSCIENCES↗
DOE OSTI2026

Optimizing infrasound observations for sample return capsule re-entry: Insights from OSIRIS-REx and Hayabusa2

The atmospheric entry of meteoroids presents a rare and unpredictable phenomenon, posing challenges for systematic observation and detailed characterization. Such events are nonetheless critical for advancing understanding of acoustic wave propagation, atmospheric structure, and entry dynamics. In contrast, sample return capsules (SRCs) from space missions follow well constrained re-entry trajectories, enabling planned observations of shock wave generation and propagation under controlled conditions. This study compares two SRC atmospheric entries, Hayabusa2 in 2020 and OSIRIS-REx in 2023, to assess how different infrasound array configurations influence shock wave detection and trajectory validation. Hayabusa2’s re-entry was monitored using a distributed network of 28 portable infrasound sensors across seven arrays in Woomera, Australia, permitting three-dimensional reconstruction of the trajectory and analysis of wave characteristics. For OSIRIS-REx, a compact four-sensor array deployed near Eureka Airport provided trajectory confirmation through arrival-time differences and back-azimuth estimates. Spectral and waveform analyses revealed differences in signal properties associated with variations in entry angle and velocity. The results illustrate both the strengths and the limitations of the deployed array configurations. The dense and distributed arrays during Hayabusa2’s re-entry enabled detailed trajectory reconstruction, whereas the compact array at Eureka primarily provided confirmation of signal coherence and back-azimuth consistency. These case studies highlight that even relatively small arrays, if located close to the predicted ground track, can still capture useful information on arrival direction and timing. Such insights provide practical guidance for planning future observational campaigns of SRC returns and other controlled atmospheric entries, and may also inform approaches to opportunistic observations of meteoroid events and other atmospheric acoustic phenomena.

celestial mechanics: orbit determination↗
DOE OSTI2021

Seismic waves in medium with poroelastic/elastic interfaces: a two-dimensional P-SV finite-difference modelling

SUMMARY We present a new methodology of the finite-difference (FD) modelling of seismic wave propagation in a strongly heterogeneous medium composed of poroelastic (P) and (strictly) elastic (E) parts. The medium can include P/P, P/E and E/E material interfaces of arbitrary shapes. The poroelastic part can be with (i) zero resistive friction, (ii) non-zero constant resistive friction or (iii) JKD model of the frequency-dependent permeability and resistive friction. Our FD scheme is capable of subcell resolution: a material interface can have an arbitrary position in the spatial grid. The scheme keeps computational efficiency of the scheme for a smoothly and weakly heterogeneous medium (medium without material interfaces). Numerical tests against independent analytical, semi-analytical and spectral-element methods prove the efficiency and accuracy of our FD modelling. In numerical examples, we indicate effect of the P/E interfaces for the poroelastic medium with a constant resistive friction and medium with the JKD model of the frequency-dependent permeability and resistive friction. We address the 2-D P-SV problem. The approach can be readily extended to the 3-D problem.

58 GEOSCIENCES↗
DOE OSTI2022

Spin and charge density waves in quasi-one-dimensional KMn 6 Bi 5

The recent observation that pressure could suppress antiferromagnetic (AFM) order in quasi-one-dimensional AMn 6 Bi 5 Mn-cluster chain materials (A=Na, K, Rb, and Cs) and lead to a superconducting dome offers an alternative Mn-based class of materials with which to study unconventional superconductivity. Using neutron diffraction, we elucidate the exact nature of the previously unknown AFM ground state of KMn 6 Bi 5 and report finding transverse incommensurate spin density waves (SDWs) for the Mn atoms with a propagating direction along the chains. The SDWs have distinct refined amplitudes of ~2.46μ B for the Mn atoms in the pentagons and ~0.29μ B with a large standard deviation for Mn atoms at the center between the pentagons. AFM coupling dominates both the nearest-neighbor Mn-Mn interactions within the pentagon and next-nearest-neighbor Mn-Mn interactions out of the pentagon (along the propagating wave). The SDWs exhibit both local and itinerant characteristics potentially due to cooperative interactions between local magnetic exchange and conduction electrons. Single crystal x-ray diffraction below the AFM transition revealed satellite peaks originating from charge density waves along the chain direction with a q vector twice as large as that of the SDW, pointing to a strong real space coupling between them. Additionally, we report a significant magnetoelastic effect during the AFM transition, especially along the chain direction, observed in temperature-dependent x-ray powder diffraction. Our work not only reveals fascinating intertwined spin, charge, and lattice orders in one-dimensional KMn 6 Bi 5 , but also provides an essential piece of information on its magnetic structure to understand the mechanism of superconductivity in this Mn-based family.

1-dimensional systems↗