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.
Three-dimensional Fréchet sensitivity kernels for electromagnetic wave propagation
Electromagnetic imaging methods are useful tools for monitoring subsurface changes in pore-fluid content and the associated changes in electrical permittivity and conductivity. The most common method for georadar tomography uses a high frequency ray-theoretic approximation that is valid when material variations are sufficiently small relative to the wavelength of the propagating wave. Georadar methods, however, often utilize electromagnetic waves that propagate within heterogeneous media at frequencies where ray theory may not be applicable. In this paper we describe the 3-D Fréchet sensitivity kernels for EM wave propagation. Various data functional types are formulated that consider all three components of the electric wavefield and incorporate near-, intermediate-, and far-field contributions. We show that EM waves exhibit substantial variations for different relative source-receiver component orientations. The 3-D sensitivities also illustrate out-of-plane effects that are not captured in 2-D sensitivity kernels and can influence results obtained using 2-D inversion methods to image structures that are in reality 3-D.
Rapid calculation of paraxial wave propagation for cylindrically symmetric optics
Not Available
Obliquely propagating waves in the magnetized strongly coupled one-component plasma
The quasi-localized charge approximation is used to calculate the wave spectrum of the magnetized three-dimensional strongly coupled one-component plasma at arbitrary angles θ between the wave vector and the magnetic field axis. Three frequency branches are identified whose interplay is strongly determined by β=ω{sub c}/ω{sub p}, the ratio of the cyclotron frequency ω{sub c}, and the plasma frequency ω{sub p}. The frequency dispersion relations for the three principal modes along the magnetic field cross in the case β<1, which strongly affects the transition from parallel to perpendicular wave propagation. For β>1, the frequencies of the different branches are well separated, and the long-wavelength dispersion in the intermediate and upper branch changes sign as θ is varied from 0 to π/2. In addition to the frequencies, we also investigate the waves' polarization properties.
Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media
It is important to develop fast yet accurate numerical methods for seismic wave propagation to characterize complex geological structures and oil and gas reservoirs. However, the computational cost of conventional numerical modeling methods, such as finite-difference method and finite-element method, becomes prohibitively expensive when applied to very large models. We propose a Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media, where we construct basis functions from multiple local problems for both the boundaries and interior of a coarse node support or coarse element. The application of multiscale basis functions can capture the fine scale medium property variations, and allows us to greatly reduce the degrees of freedom that are required to implement the modeling compared with conventional finite-element method for wave equation, while restricting the error to low values. We formulate the continuous Galerkin and discontinuous Galerkin formulation of the multiscale method, both of which have pros and cons. Applications of the multiscale method to three heterogeneous models show that our multiscale method can effectively model the elastic wave propagation in anisotropic media with a significant reduction in the degrees of freedom in the modeling system.
The detection of upwardly propagating waves channeling energy from the chromosphere to the low corona
There have been ubiquitous observations of wave-like motions in the solar atmosphere for decades. Recent improvements to space- and ground-based observatories have allowed the focus to shift to smaller magnetic structures on the solar surface. In this paper, high-resolution ground-based data taken using the Swedish 1 m Solar Telescope is combined with co-spatial and co-temporal data from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO) satellite to analyze running penumbral waves (RPWs). RPWs have always been thought to be radial wave propagation that occurs within sunspots. Recent research has suggested that they are in fact upwardly propagating field-aligned waves (UPWs). Here, RPWs within a solar pore are observed for the first time and are interpreted as UPWs due to the lack of a penumbra that is required to support RPWs. These UPWs are also observed co-spatially and co-temporally within several SDO/AIA elemental lines that sample the transition region and low corona. The observed UPWs are traveling at a horizontal velocity of around 17 ± 0.5 km s{sup –1} and a minimum vertical velocity of 42 ± 21 km s{sup –1}. The estimated energy of the waves is around 150 W m{sup –2}, which is on the lower bound required to heat the quiet-Sun corona. This is a new, yet unconsidered source of wave energy within the solar chromosphere and low corona.
A staggered-grid finite-difference scheme optimized in the time–space domain for modeling scalar-wave propagation in geophysical problems
For modeling scalar-wave propagation in geophysical problems using finite-difference schemes, optimizing the coefficients of the finite-difference operators can reduce numerical dispersion. Most optimized finite-difference schemes for modeling seismic-wave propagation suppress only spatial but not temporal dispersion errors. We develop a novel optimized finite-difference scheme for numerical scalar-wave modeling to control dispersion errors not only in space but also in time. Our optimized scheme is based on a new stencil that contains a few more grid points than the standard stencil. We design an objective function for minimizing relative errors of phase velocities of waves propagating in all directions within a given range of wavenumbers. Dispersion analysis and numerical examples demonstrate that our optimized finite-difference scheme is computationally up to 2.5 times faster than the optimized schemes using the standard stencil to achieve the similar modeling accuracy for a given 2D or 3D problem. Compared with the high-order finite-difference scheme using the same new stencil, our optimized scheme reduces 50 percent of the computational cost to achieve the similar modeling accuracy. This new optimized finite-difference scheme is particularly useful for large-scale 3D scalar-wave modeling and inversion.
Numerical investigation of fast-wave propagation and radio-frequency sheath interaction with a shaped tokamak wall
Interactions between propagating fast waves and radio-frequency (RF) sheaths in the ion cyclotron range of frequencies are numerically investigated based on a cold fluid plasma model coupled with a sheath boundary condition. In this two-dimensional study, the capability of the finite element code rfSOL, which was developed in previous numerical work, is extended to analyze self-consistent RF sheath-plasma interaction problems in a tokamak with a non-circular cross-section. It is found that a large sheath voltage is generated near the edges of the limiter-shaped deformation as a result of the conversion from fast to slow waves on the sheaths. The sheath voltage associated with this conversion is particularly significant in the localized region where the contact angle between the magnetic field line and the conducting wall varies rapidly along the curved sheath surface, which is consistent with the results in previous one-dimensional theoretical work. The dependences of the RF sheaths on various parameters in plasma such as the toroidal wavenumber, edge plasma density, and the degree of the RF wave absorption in the core region are also examined in detail.
Frequency-selective plasmonic wave propagation through the overmoded waveguide with photonic-band-gap slab arrays
Confined propagation of guided waves through the periodically corrugated channel sandwiched between two staggered dielectric photonic-band-gap slab arrays is investigated with the band-response analysis. Numerical simulations show that longitudinally polarized evanescent waves within the band gap propagate with insertion loss of {approx}-0.2 to 1 dB (-0.05 to 0.4 dB/mm at G-band) in the hybrid band filter. This structure significantly suppresses low energy modes and higher-order-modes beyond the band-gap, including background noises, down to {approx}-45 dB. This would enable the single-mode propagation in the heavily over-moded waveguide (TEM-type), minimizing abnormal excitation probability of trapped modes. This band filter could be integrated with active and passive RF components for electron beam and optoelectronic devices.
Compact representations of partially coherent undulator radiation suitable for wave propagation
Not Available
A Variable Grid FDM Technique for Treating Solid-Liquid and Solid-Air Boundaries in Wave Propagation Modeling of Elastic Waves. Task 2, Computer Program
In Task1 we explained the layout of the staggered grid adopted in the UL0_2.4 a finite difference computer program developed by Pitarka (1996,1999). The objective of the second task is to modify the UL0_2.4 program in order to extend its application to computing elastic wave propagation in heterogeneous structures with irregular surface topography and solid-liquid boundaries. The new program UL0_2.4W preserves the structure and variable-grid finite-difference technique used in the original computer program. The list of computer programs for running the FD code, and programs for post processing the output data is given in the Appendix.
Experimental realization of an additively manufactured monatomic lattice for studying wave propagation
Increasing interest in wave propagation in phononic systems and metamaterials motivates the development of experimental designs, measurement techniques, and fabrication methods for use in basic research and classroom demonstrations. The simplest phononic system, the monatomic chain, exhibits rich physics such as dispersion and frequency-domain filtering. However, a limited number of experimental studies showcase monatomic chains for macroscale observation of phonons. Herein, we discuss the design, fabrication, and testing of monatomic lattices as enabled by three-dimensional (3D) printing. Using this widely available technology, we provide design guidelines for realization of a monatomic chain composed of 3D printed serpentine springs and press-fitted cylindrical masses. We also present measurement techniques that record propagating waves and algorithms for the experimental determination of dispersion behavior.
Serpentine: Finite Difference Methods for Wave Propagation in Second Order Formulation
Wave propagation phenomena are important in many DOE applications such as nuclear explosion monitoring, geophysical exploration, estimating ground motion hazards and damage due to earthquakes, non-destructive testing, underground facilities detection, and acoustic noise propagation. There are also future applications that would benefit from simulating wave propagation, such as geothermal energy applications and monitoring sites for carbon storage via seismic reflection techniques. In acoustics and seismology, it is of great interest to increase the frequency bandwidth in simulations. In seismic exploration, greater frequency resolution enables shorter wave lengths to be included in the simulations, allowing for better resolution in the seismic imaging. In nuclear explosion monitoring, higher frequency seismic waves are essential for accurate discrimination between explosions and earthquakes. When simulating earthquake induced motion of large structures, such as nuclear power plants or dams, increased frequency resolution is essential for realistic damage predictions. Another example is simulations of micro-seismic activity near geothermal energy plants. Here, hydro-fracturing induces many small earthquakes and the time scale of each event is proportional to the square root of the moment magnitude. As a result, the motion is dominated by higher frequencies for smaller seismic events. The above wave propagation problems are all governed by systems of hyperbolic partial differential equations in second order differential form, i.e., they contain second order partial derivatives of the dependent variables. Our general research theme in this project has been to develop numerical methods that directly discretize the wave equations in second order differential form. The obvious advantage of working with hyperbolic systems in second order differential form, as opposed to rewriting them as first order hyperbolic systems, is that the number of differential equations in the second order system is significantly smaller. Another issue with re-writing a second order system into first order form is that compatibility conditions often must be imposed on the first order form. These (Saint-Venant) conditions ensure that the solution of the first order system also satisfies the original second order system. However, such conditions can be difficult to enforce on the discretized equations, without introducing additional modeling errors. This project has previously developed robust and memory efficient algorithms for wave propagation including effects of curved boundaries, heterogeneous isotropic, and viscoelastic materials. Partially supported by internal funding from Lawrence Livermore National Laboratory, many of these methods have been implemented in the open source software WPP, which is geared towards 3-D seismic wave propagation applications. This code has shown excellent scaling on up to 32,768 processors and has enabled seismic wave calculations with up to 26 Billion grid points. The WPP calculations have resulted in several publications in the field of computational seismology, e.g.. All of our current methods are second order accurate in both space and time. The benefits of higher order accurate schemes for wave propagation have been known for a long time, but have mostly been developed for first order hyperbolic systems. For second order hyperbolic systems, it has not been known how to make finite difference schemes stable with free surface boundary conditions, heterogeneous material properties, and curvilinear coordinates. The importance of higher order accurate methods is not necessarily to make the numerical solution more accurate, but to reduce the computational cost for obtaining a solution within an acceptable error tolerance. This is because the accuracy in the solution can always be improved by reducing the grid size h. However, in practice, the available computational resources might not be large enough to solve the problem with a low order method.
Hybrid SAFE-GMM approach for predictive modeling of guided wave propagation in layered media
Not provided.
Identification of the dominant ULF wave mode and generation mechanism for obliquely propagating waves in the Earth's foreshock
Abstract We discuss mechanisms of the generation of ultralow frequency (ULF) upstream waves in the terrestrial foreshock that are essential for the acceleration of ions in space plasmas. The analysis is based on global hybrid kinetic simulations of the magnetosphere that provide realistic environment for the growth of the ULF waves in a quasi‐radial configuration of the interplanetary magnetic field. We focus on a long‐debated problem of the generation mechanism of oblique and parallel ULF waves and provide quantitative arguments in favor of the ion/ion cyclotron resonant instability. We also show that parallel propagating waves are predominantly generated in this configuration, but geometrical effects related to the phase space density in wave vector space lead to apparent predominance of obliquely propagating waves. Correspondence between the results outlined above and previously published experimental claims is thoroughly discussed and our results are shown to be consistent with spacecraft measurements.
GLOBAL AND LOCAL CUTOFF FREQUENCIES FOR TRANSVERSE WAVES PROPAGATING ALONG SOLAR MAGNETIC FLUX TUBES
It is a well-established result that the propagation of linear transverse waves along a thin but isothermal magnetic flux tube is affected by the existence of the global cutoff frequency, which separates the propagating and non-propagating waves. In this paper, the wave propagation along a thin and non-isothermal flux tube is considered and a local cutoff frequency is derived. The effects of different temperature profiles on this local cutoff frequency are studied by considering different power-law temperature distributions, as well as the semi-empirical VAL C model of the solar atmosphere. The obtained results show that the conditions for wave propagation strongly depend on the temperature gradients. Moreover, the local cutoff frequency calculated for the VAL C model gives constraints on the range of wave frequencies that are propagating in different parts of the solar atmosphere. These theoretically predicted constraints are compared to observational data and are used to discuss the role played by transverse tube waves in the atmospheric heating and dynamics, and in the excitation of solar atmospheric oscillations.
Evanescent waves propagation along a periodically corrugated surface and their amplification by relativistic electron beam (quasi-optical theory)
By using a quasi-optical approach, we study propagation of evanescent waves along a periodically corrugated surface and their excitation by relativistic electron beams. Under assumption of a shallow (in the scale of period) corrugation, the dispersion equation for normal waves is derived and two particular cases are studied. In the first case, the wave frequency is far from the Bragg resonance; therefore, the evanescent wave propagation can be described by using the impedance approximation with deceleration of the zeroth spatial harmonic. The second case takes place at the frequencies close to the Bragg resonance. There, the field can be represented as two counter-propagating quasi-optical wave beams, which are coupled on the corrugated surface and form an evanescent normal wave. With regard to the interaction with an electron beam, the first case corresponds to the convective instability that can be used for amplification of radiation, while the second case corresponds to the absolute instability used in surface-wave oscillators. This paper is focused on studying main features of amplifier schemes, such as the increments, electron efficiency, and formation of a self-consistent spatial structure of the radiated field. For practical applications, the feasibility of realization of relativistic surface-wave amplifiers in the submillimeter wavelength range is estimated.
Modeling Frequency–Independent Q Viscoacoustic Wave Propagation in Heterogeneous Media
Quantifying the attenuation of seismic waves propagating in the Earth interior is critical to study the subsurface structure. Previous studies have proposed fractional anelastic wave equations to model the frequency–independent Q seismic wave propagation. Such wave equations involve fractional derivatives that pose computational challenges for the numerical schemes in terms of accuracy and efficiency when dealing with heterogeneous Earth media. To tackle these challenges, here we derive a new viscoacoustic wave equation, where the power terms of the fractional Laplacian operators are spatially independent, thus accurate and efficient methods (e.g., the Fourier pseudospectral method) can be adopted. Our derivation enables the resultant equation to capture both amplitude and phase signatures of the anelastic wave propagation by matching the complex wave numbers for all the frequencies of interest. We verify the derivation by comparing the dispersion curves of both the attenuation factor and the phase velocity produced by the new wave equation with their theoretical values as well as the Pierre Shale in situ measurements. Following that, we use a synthetic attenuating gas chimney model to demonstrate the attenuation effects on seismic waveforms and then construct the Q–compensated reverse time migration to undo these effects for seismic image enhancement. Finally, we find that our forward modeling results can characterize the spatiotemporal attenuation effects revealed in the Frio–II CO2 injection time–lapse seismic monitoring data. Here, we expect this proposed equation to be useful to quantify the attenuation in seismic data to push the resolution limits of seismic imaging and inversion.