Search NASA⌕ Search

SEARCH · Search NASA

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 37 records · Page 2

Propagation of elastic waves in a flexomagnetic solid

Flexomagnetism is the coupling between magnetism and strain gradients and is a technologically relevant phenomenon. We present a theory of elastic wave propagation in a linear elastic flexomagnetic material with microstructure and strain gradient elastic interactions. The expressions of frequency, phase velocity, and group velocity of longitudinal and transverse waves are derived and are shown to depend on the flexomagnetic coefficient and microstructure. We also show that the effect of flexomagnetism and microstructure can lead to some interesting phenomena in wave propagation, which are not observed in classical linear elasticity theory of waves. Specifically, in contrast to classical linear elastic materials, where wave propagation is non-dispersive, flexomagnetic materials with microstructure can exhibit both normal and abnormal dispersion. It is also noteworthy that, in flexomagnetic materials with gradient elasticity, the phase velocities of transverse waves can exceed those of longitudinal waves, which is atypical in classical elasticity. Furthermore, waves can also attenuate for a certain range of wavenumbers that depend on the flexomagnetic coefficient and microstructural parameters. Finally, we explore the possibility of waves exhibiting zero group velocity modes, where waves are non-propagating but have strong local energy confinement, negative group velocity modes, where the wave packet moves in the opposite direction to that of wave propagation, and the phenomenon of wave freezing, where a propagating wave stops in space without diffusing or spreading.

Ghosh, Swarnava [ORNL] (ORCID:0000000338005264)↗

Propagation of EMIC Waves From Shabansky Orbits in the Dayside Magnetosphere

We explore the characteristics of EMIC waves generated in a non-dipole, compressed magnetic field at the minimum of the magnetic field. We conducted 2D full-wave simulations using the Petra-M code, focusing on a compressed magnetic field in the outer dayside magnetosphere for a range of L values (L = 8 - 10). By comparing the simulation results with MMS observations, we aim to understand how the observed wave characteristics are affected by a shifting source region across different L-shells. Our findings indicate that the direction of the Poynting vector systematically changes depending on the local source location of the wave, which is consistent with the observations. EMIC waves propagate along the magnetic field line and reach both the northern and southern hemispheres; however, there is a notable difference in the power of EMIC waves between the two hemispheres, indicating seasonal asymmetries in their occurrence.

58 GEOSCIENCES↗

An integrated approach to optimizing concentration shock wave electrodialysis using 2D multicell simulation and response surface models

Shock wave electrodialysis (SWED) is a highly promising technique for energy-efficient ion separation in the context of a circular economy. This paper presents a approach way of modeling and improving SWED using a two-dimensional multicell model combined with the COMSOL program and response surface methodology. The model integrates the Nernst-Planck equation, Darcy's law, and first-order electroosmosis to examine the local concentration, flux of ionic species, distribution of current, and velocity of flow in SWED cells under various operating conditions. We first illustrate the clear depiction of concentration, velocity, and electric potential distribution through contours which aids in identifying optimal operating conditions and designing scalable SWED systems. The results emphasize the significance of surface charge density and voltage in influencing the features of shock waves for obtaining effective ion separation while optimizing energy consumption and improving current efficiency by controlling the retention time of feed flow. Here, this study defines two crucial characteristics of shock waves, namely the length of the flat depletion zone of a fully developed shock wave (shock wave height) and the distance of shock wave propagation (shock wave length). These properties significantly impact separation performance, as determined by the simulation results. Additionally, the response surface methodology is incorporated with the COMSOL models to develop predictive models and graph responses, enabling a more comprehensive understanding of the interactions between parameters and performance indicators, such as removal ratio, energy consumption, and water recovery. Finally, this work suggests design tactics for expanding SWED processes and outlines potential areas for further research. This research provides valuable insights into the prospective applications, design optimization, and scalability of SWED in the field of electrokinetic separation technologies for green chemistry and a circular economy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Applicability of 1D site response analysis to shallow sedimentary basins: A critical evaluation through physics-based 3D ground motion simulations

One-dimensional site response analysis (1D SRA) remains the standard practice in considering the effect of local soil deposits and predicting site-specific ground motions, although its range of applicability to realistic seismic wavefields is still in question. In this 1D approach, horizontal and vertical ground shaking are assumed to be induced by vertically propagating shear and compressional waves, respectively. A recent study based on analytical two-dimensional (2D) plane waves and simple point source earthquake simulations has shown two mechanistic limitations in this 1D modelling technique for general inclined seismic waves, that is, systematic over-prediction of the vertical motion and wave trapping in the 1D soil column. In this article, we evaluate in detail the applicability of this 1D modelling approach to realistic three-dimensional (3D) simulated seismic wavefields in shallow sedimentary basins. Linear-viscoelastic 1D SRA predictions using two types of input motions that are commonly used in practice—rock outcrop and in-column motions, are compared with the reference true site response results from 3D earthquake simulations in terms of various measures in the frequency and time domain. It is shown that the horizontal motion in the 3D seismic wavefield exhibits dominant shear wave propagation phenomenon, while the vertical motion is a combined effect of compressional and shear waves and can be over-predicted by the 1D approach when the incident seismic waves are inclined. Direct evidence of the wave refraction process that leads to the vertical motion over-prediction is provided. 1D SRA with in-column inputs can yield motions that have significantly longer duration compared to the true 3D site response solution due to trapped waves, casting in doubt the frequent need for increased soil damping in existing site studies to compensate for wave attenuation due to scattering alone. Sensitivity investigation on the increase of soil profile damping by a multiplier D mul shows D mul values compatible with those found in the literature for both horizontal and vertical motions. It is shown that the level of D mul optimized for a best match of the spectral acceleration is dependent on the characteristic of the input motion and a larger D mul is typically required for the vertical component. In contrast, 1D SRA with outcrop motions predicts motions with shorter significant duration due to its inability to capture the basin-edge generated surface waves. A suite of ground motion simulations was performed to assess the sensitivity of the observations to the basin geologic structure including the velocity gradient, rock-basin impedance contrast and basin depth. The analysis results show that the accuracy of the simplified 1D procedure is dependent on the wavefield composition of both the input motions and the true 3D site response solution. While the horizontal motions in shallow sedimentary basins can, to the first order, be reasonably captured by the simplified 1D approach, 1D SRA for the vertical component is in general not reliable and contributions from inclined shear waves should be accounted for in site-specific evaluation of the vertical design ground motion.

1D assumption↗

High-Fidelity Numerical Wave Tank Verification & Validation Study: Wave Generation Through Paddle Motion: Preprint

This paper presents a numerical benchmark study of wave propagation due to a paddle motion using different high-fidelity numerical models, which are capable of replicating the nearly actual physical wave tank testing. A full time series of the measured wave generation paddle motion which was used to generate wave propagation in the physical wave tank will be utilized in each of the models contributed by IEA OES Task 10's participants, which includes both computational fluid dynamics (CFD) and smooth hydrodynamic particle (SPH). The high-fidelity simulations of the physical wave testcase will allow for the evaluation of the initial transient effects from wave ramp-up and its evolution in the wave tank over time for two representative regular waves with varying levels of nonlinearity. A couple of interesting metrics like the predicted wave surface elevation at select wave probes, wave period, and phase-shift in time will be assessed to evaluate the relative accuracy of numerical models versus experimental data within specified time intervals. These models will serve as a guide for modelers in the wave energy community and provide a base case to allow further and more detailed numerical modeling of the fixed Kramer Sphere Cases under wave excitation force wave tank testing.

HYDRO ENERGY,TIDAL AND WAVE POWER↗

Ultrasonic characterization of material heterogeneities in stainless steel components produced by laser powder bed fusion

We introduce pulse-echo ultrasound as a method for characterizing the impact of powder bed fusion parameters on the properties of additively manufactured stainless-steel components, their material anisotropy, and location-dependent heterogeneity. Our results indicate that accurate characterization requires careful selection of ultrasonic propagation paths, which must consider the direction of additive layering, variations in processing parameters, and the component's geometry. We employed two distinct methods to estimate material properties from ultrasonic data: One assumes isotropy, while the other accounts for anisotropic interactions during the propagation of elastic waves. When applied to samples fabricated with laser energy densities ranging from 24 to 42 J/mm³ , these methods revealed transverse isotropy and weak anisotropy (quantified by small Thomsen parameters, ε = 0.0651 and γ = 0.0092) and less than a ∼6 % change in acoustic impedance. The assumption of isotropy, in this case, leads to small errors (less than 4 % or 1 % for Young's modulus in the build or transverse directions) when estimating orthotropic material properties using ultrasonic data measured along just two orthogonal directions, one of which must align with the build direction. By comparing ultrasonic measurements — which aggregate the spatial variability in material properties along the length of elastic wave propagation into a single value — with localized measurements obtained from surface nanoindentation, we uncovered and spatially profiled significant differences between the surface and interior properties. Specifically, the surface Young's modulus decreased from approximately 210 GPa to 180 GPa within a depth of about 3 mm. We attribute this surface-localized heterogeneity in PBF-fabricated components to distinct thermal histories experienced by the surface and interior regions. Collectively, the results of this study establish a framework for the ultrasonic characterization of material heterogeneity and anisotropy in material properties and demonstrate its application in additively manufactured metal components.

36 MATERIALS SCIENCE↗

Tailored ultrasound propagation in microscale metamaterials via inertia design

The quasi-static properties of micro-architected (meta)materials have been extensively studied over the past decade, but their dynamic responses, especially in acoustic metamaterials with engineered wave propagation behavior, represent a new frontier. However, challenges in miniaturizing and characterizing acoustic metamaterials in high-frequency (megahertz) regimes have hindered progress toward experimentally implementing ultrasonic-wave control. Here, we present an inertia design framework based on positioning microspheres to tune responses of 3D microscale metamaterials. We demonstrate tunable quasi-static stiffness by up to 75% and dynamic longitudinal-wave velocities by up to 25% while maintaining identical material density. Using noncontact laser-based dynamic experiments of tunable elastodynamic properties and numerical demonstrations of spatio-temporal ultrasound wave propagation, we explore the tunable static and elastodynamic property relation. This design framework expands the quasi-static and dynamic metamaterial property space through simple geometric changes, enabling facile design and fabrication of metamaterials for applications in medical ultrasound and analog computing.

42 ENGINEERING↗

Mantaray: A Rust Package for Ray Tracing Ocean Surface Gravity Waves

Ocean surface gravity waves are an important component of air-sea interaction, influencing energy, momentum, and gas exchanges across the ocean-atmosphere interface. In specific applications such as refraction by ocean currents or bathymetry, ray tracing provides a computationally efficient way to gain insight into wave propagation. In this paper, we introduce Mantaray, an open-source software package implemented in Rust, with a Python interface, that solves the ray equations for ocean surface gravity waves. Mantaray is designed for performance, robustness, and ease of use. The package is modular to facilitate further development and can currently be applied to both idealized and realistic wave propagation problems (Fig. 1).

16 TIDAL AND WAVE POWER↗

Final Technical Report - Center for Simulation of Fusion Relevant RF Actuators

We have developed a suite of 3D electromagnetic field solvers, both FEM and FDTD based, that account for the RF antenna and vacuum vessel geometries with unprecedented accuracy. Workflows were developed that make it possible to translate CAD models for the antenna and vacuum vessel to physics meshes for RF wave simulation. Nonlinear RF sheath formation has been incorporated self-consistently as a boundary condition in these solvers. We have also carried out extensive studies of the impact of RF sheaths on the ion energy angle distribution at plasma-material interfaces, using high fidelity particle-in-cell codes. Comprehensive simulation models were developed to assess the impact of blob-like edge turbulence on RF wave propagation and the impact of the RF ponderomotive force on the plasma scrape-off layer (SOL). A fluid transport solver for the far-SOL was also developed which accounts for the high parallel to perpendicular heat anisotropy on an unstructured mesh, thus making it possible to precisely represent an antenna structure in the presence of edge transport. Finally we have developed a hierarchy of core wave propagation and absorption models that self-consistently combine continuum Fokker Planck and Monte Carlo treatments of fast ion evolution with ICRF full-wave field solvers and continuum Fokker Planck treatments of fast electron evolution with both full-wave and ray tracing models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation of Physics-Based 0-10Hz Strong Motion Using High Performance Computing Supporting Refinements to Regional Ground Motion Models for the Central Eastern US

In collaboration with the U.S. Nuclear Regulatory Commission (NRC) the LLNL has developed a computationally efficient simulation platform designed to perform physics-based ground motion simulations for crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using high-performance computing. The main objective of the earthquake simulations was to use synthetic ground motion to provide constrains to refinements of existing ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-10Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model (1994). In simulations performed during the first stage of this project the Saikia’s velocity model demonstrated better performance in modelling high frequency regional wave propagation for the CEUS region recorded during the Mw5.0 November 7, 2016, Cushing Oklahoma (Taylor et al., 2017), and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes. The proposed regional 3D model includes random perturbations to the 1D background model using the stochastic scheme of Pitarka and Mellors (2021). In addition, validation analysis of the rupture generator and regional wave propagation models, using comparisons with different GMMs for Mw6.5 and Mw7.0 scenario earthquakes in the CEUS region resulted in a very good match between the simulated and empirical ground motion models. For the purposes of seismic hazard assessment at the existing and planned nuclear power plants, NRC is interested in studies aimed at improving the current ground motion models (GMM) for both Stable Continental Regions (SCR) in the Central and Eastern US and Active Crustal Regions (ACR) in the Western US. Due to lack of recorded data, these improvements require synthetic data for short fault distances and large magnitude earthquakes for which the existing recorded data is not enough to uniquely constrain the GMMs. The need for simulations and strong motion data is especially critical for the CEUS region where we do not have recorded data from potentially large damaging earthquakes with moment magnitudes 6.0 and higher. In this the project, we focused on 10Hz simulations of Mw7.0 scenario earthquakes with strike slip and thrust faulting mechanisms. We used more than 50 Mw7.0 earthquake rupture scenarios to investigate the ground motion uncertainty due to unknown earthquake rupture parameters, in particular, the slip distribution, rupture velocity, and faulting mechanism, and their implication on ground motion amplification due to forward rupture directivity effects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Experimental Study of Instabilities in Hydrogen-Air Fueled Rotating Detonation Combustion Presentation

Conventional gas turbine engines rely on an idealized constant pressure combustion process that in reality produces a pressure decrease as a result of viscous and other non-reversible losses. An alternative approach is rotating detonation combustion (RDC) which is a form of pressure gain combustion in which one or more detonation waves propagate an annular channel resulting in an increase in pressure across, subsequently providing greater work availability compared to deflagration ultimately leading to opportunities for greater thermodynamic efficiency when used in gas turbine engines that conventionally relies on constant. Modern gas turbine engines often rely on pre-mixed reactants to limit NOx emissions, although this may result in greater susceptibility to instabilities such as flashback and thermoacoustic oscillation, particularly for applications that utilize hydrogen as the fuel. Research in RDC has focused on non-premixed reactants thus limiting the occurrence of flashback, and high frequency detonation wave propagation (kHz) may interfere with the occurrence of thermoacoustic oscillations. Thermal NOx emissions are controlled through rapid combustion and sudden expansion of the working fluid. Although RDC may not be susceptible to instabilities encountered in conventional state of the art gas turbine engine combustion, there may be other mechanisms occurring that support instabilities that could be detrimental to performance.

Weber, Justin↗

Performance evaluation of the USGS velocity model for the San Francisco Bay Area

In this study, we evaluated the performance of the United States Geological Survey velocity model developed for the San Francisco Bay Area (SFBA), version 21.1. The evaluation was performed through high-resolution three-dimensional physics-based ground motion simulations of seven small-magnitude earthquakes (ranging from magnitude 3.8 to 4.4) that occurred on the eastern side of the San Francisco Bay. The simulations were performed in the frequency range from 0 to 5 Hz with a minimum shear-wave velocity of 250 m/s, which allowed the capture of wave propagation effects of the near-surface soft materials that characterize local basins. Based on the direct comparison of Fourier amplitude spectra between recorded and simulated ground motions for more than 250 stations, we found that the velocity model generally performs well in the frequency range of 0.2–5 Hz. The median value of the Fourier amplitude residuals was found to be near zero for all seven earthquakes. The slight over-prediction of 0.2 log-natural units at frequencies above 3 Hz in our simulations was attributed to the potentially inaccurate representation of the source radiation pattern by a double-couple point source model, and simple representation of shallow small-scale underground structural complexity in the velocity model. Maps of spectral amplitude differences between the simulated and recorded data were used to identify areas responsible for systematic ground motion over-predictions or under-predictions. For example, while some sub-domains over soft sediments show over-prediction patterns, the block east of the Hayward fault is prone to exhibit patterns of under-prediction. These maps can be used to guide future refinements of the SFBA velocity model. Since our simulation methodology allows for the decoupling of the source and wave propagation effects, the ground motion data generated by our simulations can also be used to quantify the epistemic uncertainty due to the velocity model, in empirically based ground motion estimates for the SFBA.

58 GEOSCIENCES↗

Reduced-action-integral approach for photon-photon interactions in vacuum

Electromagnetic waves propagating through vacuum can polarize virtual electron–positron pairs; this polarization, in turn, nonlinearly modifies their propagation. A semi-classical nonlinear wave equation describing the propagation is derived from the Euler–Heisenberg Lagrangian density, which captures vacuum polarization effects up to the one-loop level. In this article, we present a reduced-actionintegral approach that enables rapid modeling of nonlinear phenomena arising from the Euler– Heisenberg Lagrangian. Application of the variational principle to the reduced action provides equations of motion for familiar light-pulse parameters, such as spot size, phase, polarization, and phase-front curvature, without requiring full-field simulations. Three examples demonstrate the utility of the approach: phase modulation, birefringence, and frequency mixing.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Thermal modulation of nonlinear ultrasonic waves for nondestructive evaluation of elastic materials

Temperature variation is often considered an undesirable factor in ultrasonic testing, as wave velocity is highly sensitive to thermal fluctuations. However, completely eliminating the temperature effect is difficult, particularly in tests requiring precise velocity measurements. Recently, a method called Thermal Modulation of Nonlinear Ultrasonics (TMNL) has been developed. Instead of eliminating the thermal effect, the TMNL method leverages the temperature variation as a driving force to stimulate the nonlinear response of the medium and modulate the ultrasonic waves propagating in it. These modulated waves can then be used to evaluate the nonlinear behaviours of the test medium. This paper presents a focused review of the TMNL technique, including its theoretical foundations, particularly the conceptual challenges in integrating thermal effects into classical acoustoelastic theory, and recent applications in non-destructive evaluation (NDE). Three case studies are presented to demonstrate TMNL’s application in detecting microcracking in concrete, assessing ageing in polymer materials, and enabling temperature compensation in acoustoelastic tests. In conclusion, the review also summarises related studies, including photothermal crack modulation, and discusses current limitations and future directions of TMNL in elastic media.

NDE↗

Coherent High-Frequency Axial Oscillations in a Partially Magnetized Direct Current Magnetron Discharge

High-frequency oscillations are observed in a neon plasma of a direct current magnetron discharge. At low discharge currents, we see highly coherent 60 MHz fluctuations. Above a distinct current threshold, secondary 5–10 MHz fluctuations emerge in addition to turbulent fluctuations in the 60–100 MHz range. The oscillations in the total discharge current suggest axial wave propagation. A lower-hybrid wave theory is invoked to model the high-frequency oscillations. Here, we attribute the low-frequency modes to a turbulence-driven inverse cascade process, as suggested by recent simulations.

33 ADVANCED PROPULSION SYSTEMS↗

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↗

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↗

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↗