Search NASA⌕ Search

SEARCH · Search NASA

Results for “wave equation”

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

Joint PP and PS plane-wave wave-equation migration-velocity analysis

Conventional joint PP and PS velocity analysis is based on ray tomography. We develop a joint PP and PS wave-equation migration-velocity-analysis method using plane-wave common-image gathers (CIGs) to produce accurate P- and S-wave velocity models. The objective function of our new method consists of three terms: The first and second terms penalize the moveout residuals computed from PP and PS plane-wave CIGs, respectively, and the third term constrains the nonzero relative depth shifts between the PP and PS migration images. The moveout of plane-wave CIGs is automatically picked using a semblance analysis method, and the relative depth shifts between the PP and PS images are automatically computed using dynamic warping or manually picking the depths of certain primary reflectors. The moveout residuals and the relative depth shifts are transformed into weighted image perturbations, and they are then projected into the velocity models to update the P- and S-wave velocity models using the scalar-wave equations and their linearized forms. Here, numerical tests with synthetic and multicomponent field data demonstrate that our method can simultaneously invert for accurate P- and S-wave velocity models for elastic migration.

58 GEOSCIENCES↗

Decoupled Fréchet kernels based on a fractional viscoacoustic wave equation

We formulate the Fréchet kernel computation using the adjoint-state method based on a fractional viscoacoustic wave equation. We first numerically prove that both the 1/2- and the 3/2-order fractional Laplacian operators are self-adjoint. Using this property, we show that the adjoint wave propagator preserves the dispersion and compensates the amplitude, while the time-reversed adjoint wave propagator behaves identically as the forward propagator with the same dispersion and dissipation characters. Without introducing rheological mechanisms, this formulation adopts an explicit Q parameterization, which avoids the implicit Q in the conventional viscoacoustic/viscoelastic full waveform inversion (Q-FWI). In addition, because of the decoupling of operators in the wave equation, the viscoacoustic Fréchet kernel is separated into three distinct contributions with clear physical meanings: lossless propagation, dispersion, and dissipation. Here, we find that the lossless propagation kernel dominates the velocity kernel, while the dissipation kernel dominates the attenuation kernel over the dispersion kernel.

58 GEOSCIENCES↗

On the Cauchy Problem for the Wave Equation with Data on the Boundary

The Cauchy problem for the wave equation in Ω × ℝ with data given on some part of the boundary ∂Ω × ℝ is considered. A reconstruction algorithm for this problem based on analytic expressions is given. This result is applicable to the problem of determining a nonstationary wave field arising in geophysics, photoacoustic tomography, tsunami wave source recovery.

97 MATHEMATICS AND COMPUTING↗

Hermite Methods for the Scalar Wave Equation

Arbitrary order dissipative and conservative Hermite methods for the scalar wave equation are presented here. Both methods use $(m+1)^d$ degrees of freedom per node for the displacement in $d$ dimensions; the dissipative and conservative methods achieve orders of accuracy $(2m-1)$ and $2m$, respectively. Stability and error analyses as well as implementation strategies for accelerators are also given.

97 MATHEMATICS AND COMPUTING↗

Stability and accuracy analysis of some fully-discrete algorithms for the one-dimensional second-order wave equation

The present investigation is concerned with some basic results for a predictor-multicorrector algorithm applied to the one-dimensional wave equation, giving particular attention to so-called 2-pass explicit schemes in which both lumped and coupled mass matrices are employed. In an assessment of the accuracy and stability properties of the algorithms, use is made of the one-dimensional, second-order wave equation. The maximum stable time step of the lumped right-hand-side mass, 2-pass explicit algorithm is twice that of the 1-pass explicit algorithm. Improved accuracy is obtained by employing higher-order, or consistent, right-hand-side, mass.

Hughes, T. J. R.↗

Working With the Wave Equation in Aeroacoustics: The Pleasures of Generalized Functions

The theme of this paper is the applications of generalized function (GF) theory to the wave equation in aeroacoustics. We start with a tutorial on GFs with particular emphasis on viewing functions as continuous linear functionals. We next define operations on GFs. The operation of interest to us in this paper is generalized differentiation. We give many applications of generalized differentiation, particularly for the wave equation. We discuss the use of GFs in finding Green s function and some subtleties that only GF theory can clarify without ambiguities. We show how the knowledge of the Green s function of an operator L in a given domain D can allow us to solve a whole range of problems with operator L for domains situated within D by the imbedding method. We will show how we can use the imbedding method to find the Kirchhoff formulas for stationary and moving surfaces with ease and elegance without the use of the four-dimensional Green s theorem, which is commonly done. Other subjects covered are why the derivatives in conservation laws should be viewed as generalized derivatives and what are the consequences of doing this. In particular we show how we can imbed a problem in a larger domain for the identical differential equation for which the Green s function is known. The primary purpose of this paper is to convince the readers that GF theory is absolutely essential in aeroacoustics because of its powerful operational properties. Furthermore, learning the subject and using it can be fun.

Farassat, F.↗

Well-posedness of one-way wave equations and absorbing boundary conditions

The one-way wave equations occurring in geophysics, underwater acoustics, and numerical studies involving absorbing boundary conditions are characterized analytically. The conditions under which such equations are well posed are obtained by examining the rational functions used to reduce them. Results are presented for Chebyshev, Pade, least-squares, and symmetric least-squares approximations.

Trefethen, Lloyd N.↗

Transient difference solutions of the inhomogeneous wave equation: Simulation of the Green's function

A time-dependent finite difference formulation to the inhomogeneous wave equation is derived for plane wave propagation with harmonic noise sources. The difference equation and boundary conditions are developed along with the techniques to simulate the Dirac delta function associated with a concentrated noise source. Example calculations are presented for the Green's function and distributed noise sources. For the example considered, the desired Fourier transformed acoustic pressures are determined from the transient pressures by use of a ramping function and an integration technique, both of which eliminates the nonharmonic pressure associated with the initial transient.

Baumeiste, K. J.↗

Transient difference solutions of the inhomogeneous wave equation - Simulation of the Green's function

A time-dependent finite difference formulation to the inhomogeneous wave equation is derived for plane wave propagation with harmonic noise sources. The difference equation and boundary conditions are developed along with the techniques to simulate the Dirac delta function associated with a concentrated noise source. Example calculations are presented for the Green's function and distributed noise sources. For the example considered, the desired Fourier transformed acoustic pressures are determined from the transient pressures by use of a ramping function and an integration technique, both of which eliminates the nonharmonic pressure associated with the initial transient.

Baumeister, K. J.↗

Higher order parabolic approximations of the reduced wave equation

Asymptotic solutions of order k to the nth are developed for the reduced wave equation. Here k is a dimensionless wave number and n is the arbitrary order of the approximation. These approximations are an extension of geometric acoustics theory, and provide corrections to that theory in the form of multiplicative functions which satisfy parabolic partial differential equations. These corrections account for the diffraction effects caused by variation of the field normal to the ray path and the interaction of these transverse variations with the variation of the field along the ray. The theory is applied to the example of radiation from a piston, and it is demonstrated that the higher order approximations are more accurate for decreasing values of k.

Mcaninch, G. L.↗

The Mean Curvature of the Influence Surface of Wave Equation With Sources on a Moving Surface

The mean curvature of the influence surface of the space-time point (x, t) appears in linear supersonic propeller noise theory and in the Kirchhoff formula for a supersonic surface. Both these problems are governed by the linear wave equation with sources on a moving surface. The influence surface is also called the Sigma - surface in the aeroacoustic literature. This surface is the locus, in a frame fixed to the quiescent medium, of all the points of a radiating surface f(x, t) = 0 whose acoustic signals arrive simultaneously to an observer at position x and at the time t. Mathematically, the Sigma- surface is produced by the intersection of the characteristic conoid of the space-time point (x, t) and the moving surface. In this paper, we derive the expression for the local mean curvature of the Sigma - space of the space-time point for a moving rigid or deformable surface f(x, t) = 0. This expression is a complicated function of the geometric and kinematic parameters of the surface f(x, t) = 0. Using the results of this paper, the solution of the governing wave equation of high speed propeller noise radiation as well as the Kirchhoff formula for a supersonic surface can be written as very compact analytic expression.

Farassat, F.↗

Fréchet kernels based on a fractional viscoacoustic wave equation

Incorporating the seismic attenuation into the waveform inversion framework could not only improve the accuracy of the velocity model but also provide an additional Q model. Recently, we proposed a viscoacoustic wave equation assisted by the fractional Laplacian operators to accurately model the wave propagation in heterogeneous attenuating media with computational efficiency. The explicit presence of Q as a coefficient in this equation suggests the potential to conveniently develop its full waveform inversion scheme. In this study, we utilize the adjoint-state method to formulate the computation of the Fréchet kernels with respect to both velocity and attenuation based on this wave equation. These kernels will play a fundamental role in the viscoacoustic multiparameter waveform inversion.

Xing, Guangchi↗

The Parabolic Wave Equation Generalized to a Spatially Varying Plasma in a Vertical Magnetic Field

This report walks through the process of deriving the wave equation for a vertical wave being propagated through a spatially varying plasma in a vertical magnetic field. Additionally, this report shows how to derive the index of refraction of the plasma corresponding to three different modes of the vertical magnetic field. Finally, this paper briefly explores the next steps planned to successfully simulate and test the equations in a pre-existing phase screen scintillation code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The nonlinear wave equation for higher harmonics in free-electron lasers

The nonlinear wave equation and self-consistent pendulum equation are generalized to describe free-electron laser operation in higher harmonics; this can significantly extend their tunable range to shorter wavelengths. The dynamics of the laser field's amplitude and phase are explored for a wide range of parameters using families of normalized gain curves applicable to both the fundamental and harmonics. The electron phase-space displays the fundamental physics driving the wave, and this picture is used to distinguish between the effects of high gain and Coulomb forces.

Colson, W. B.↗

An implementation of a high-order generalized finite difference method for solving the time-harmonic cold plasma wave equation in toroidal geometry

A high-order physics-informed meshless finite difference numerical technique is introduced for solving the time-harmonic cold plasma wave equation in toroidal geometries, presenting a novel application of the generalized finite difference (GFD) method to plasma wave simulations. The algorithm employs an irregular distribution of computational points, with local point density informed by the shortest wavelength derived from the cold plasma dispersion relation. Numerical stability and robustness are addressed using regularization techniques. The algorithm, implemented for two spatial dimensions, solves for the wave electric field and is demonstrated to achieve convergence rates of $\mathcal{O}$($\mathcal{h}$ $\mathcal{P}$ )⁠. Verification tests reproduce plane wave solutions, and example simulations of ion cyclotron resonance heating and electron cyclotron resonance heating demonstrate its capability, approaching realistic tokamak plasma scenarios. This work contributes to laying a foundation for the GFD method to be used in more sophisticated, optimized, and physically realistic full-wave simulations in time-harmonic plasma wave research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Exact finite difference schemes for the non-linear unidirectional wave equation

Attention is given to the construction of exact finite difference schemes for the nonlinear unidirectional wave equation that describes the nonlinear propagation of a wave motion in the positive x-direction. The schemes constructed for these equations are compared with those obtained by using the usual procedures of numerical analysis. It is noted that the order of the exact finite difference models is equal to the order of the differential equation.

Mickens, R. E.↗