Search NASA⌕ Search

SEARCH · Search NASA

Results for “WAVE EQUATIONS”

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 109 records · Page 6

Finite difference identification of noisy distributed systems using scanning measurements

Most of the present-day literature concerned with identification theory and techniques is directed toward lumped parameter systems, and many comprehensive surveys of the field are available. Relatively little has appeared in the literature concerning distributed identification, and even more noticeable is the scarcity of papers dealing with systems described by the one-dimensional wave equation. Perdeauville and Goodson were perhaps the first researchers with a workable but time consuming method for the identification of coefficients of the wave equation. Fairman and Shen, also considering the wave equation, used the technique of finite differencing to approximate spatial derivatives, and Poisson filter chains to approximate temporal derivatives.

Hughes, R. O.↗

A Fast Method of Deriving the Kirchhoff Formula for Moving Surfaces

The Kirchhoff formula for a moving surface is very useful in many wave propagation problems, particularly in the prediction of noise from rotating machinery. Several publications in the last two decades have presented derivations of the Kirchhoff formula for moving surfaces in both time and frequency domains. Here we present a method originally developed by Farassat and Myers in time domain that is both simple and direct. It is based on generalized function theory and the useful concept of imbedding the problem in the unbounded three-dimensional space. We derive an inhomogeneous wave equation with the source terms that involve Dirac delta functions with their supports on the moving data surface. This wave equation is then solved using the simple free space Green's function of the wave equation resulting in the Kirchhoff formula. The algebraic manipulations are minimal and simple. We do not need the Green's theorem in four dimensions and there is no ambiguity in the interpretation of any terms in the final formulas. Furthermore, this method also gives the simplest derivation of the classical Kirchhoff formula which has a fairly lengthy derivation in physics and applied mathematics books. The Farassat-Myers method can be used easily in frequency domain.

Farassat, F.↗

Acoustic propagation in a thermally stratified atmosphere

Acoustic propagation in an atmosphere with a specific form of temperature profile has been investigated by analytical means. The temperature profile used is representative of an actual atmospheric profile and contains three free parameters. Both lapse and inversion cases have been considered. Although ray solution have been considered the primary emphasis has been on solutions of the acoustic wave equation with point force where the sound speed varies with height above the ground corresponding to the assumed temperature profile. The method used to obtain the solution of the wave equation is based on Hankel transformation of the wave equation, approximate solution of the transformed equation for wavelength small compared to the scale of the temperature (or sound speed) profile, and approximate or numerical inversion of the Hankel transformed solution. The solution displays the characteristics found in experimental data but extensive comparison between the models and experimental data has not been carried out.

Vanmoorhem, W. K.↗

Acoustic propagation in a thermally stratified atmosphere

Acoustic propagation in an atmosphere with a specific form of a temperature profile has been investigated by analytical means. The temperature profile used is representative of an actual atmospheric profile and contains three free parameters. Both lapse and inversion cases have been considered. Although ray solutions have been considered, the primary emphasis has been on solutions of the acoustic wave equation with point source where the sound speed varies with height above the ground corresponding to the assumed temperature profile. The method used to obtain the solution of the wave equation is based on Hankel transformation of the wave equation, approximate solution of the transformed equation for wavelength small compared to the scale of the temperature (or sound speed) profile, and approximate or numerical inversion of the Hankel transformed solution. The solution displays the characteristics found in experimental data but extensive comparison between the models and experimental data has not been carried out.

Vanmoorhem, W. K.↗

The Extended Parabolic Equation Method and Implication of Results for Atmospheric Millimeter-Wave and Optical Propagation

The extended wide-angle parabolic wave equation applied to electromagnetic wave propagation in random media is considered. A general operator equation is derived which gives the statistical moments of an electric field of a propagating wave. This expression is used to obtain the first and second order moments of the wave field and solutions are found that transcend those which incorporate the full paraxial approximation at the outset. Although these equations can be applied to any propagation scenario that satisfies the conditions of application of the extended parabolic wave equation, the example of propagation through atmospheric turbulence is used. It is shown that in the case of atmospheric wave propagation and under the Markov approximation (i.e., the -correlation of the fluctuations in the direction of propagation), the usual parabolic equation in the paraxial approximation is accurate even at millimeter wavelengths. The methodology developed here can be applied to any qualifying situation involving random propagation through turbid or plasma environments that can be represented by a spectral density of permittivity fluctuations.

Manning, Robert M.↗

Analysis of acoustic radiation in a jet flow environment

The analysis of sound fields from arbitrary source distributions in terms of Legendre and spherical Hankel functions is well known. The purpose of this paper is to extend this classical method of analysis to environments such as jet flows where flow and flow gradients are inherently present. The wave-equation governing the radiation of sound in such an environment is derived. The steady state flow and flow gradients in the axial and transverse directions appear as coefficients in the terms of the wave-equation. A semi-numerical method is used to solve the wave-equation in terms of modified spherical harmonics yielding the phase velocities and the directivities of an infinite set of modes. The directivity of each mode is obtained in terms of modified Legendre functions by numerical integration. Some results of these directivity and phase-velocity calculations are presented for a limited number of frequency and flow parameters. Both convective and shear refraction are shown to be important.

Mungur, P.↗

A Kosloff/Basal method, 3D migration program implemented on the CYBER 205 supercomputer

Conventional finite difference migration has relied on approximations to the acoustic wave equation which allow energy to propagate only downwards. Although generally reliable, such approaches usually do not yield an accurate migration for geological structures with strong lateral velocity variations or with steeply dipping reflectors. An earlier study by D. Kosloff and E. Baysal (Migration with the Full Acoustic Wave Equation) examined an alternative approach based on the full acoustic wave equation. The 2D, Fourier type algorithm which was developed was tested by Kosloff and Baysal against synthetic data and against physical model data. The results indicated that such a scheme gives accurate migration for complicated structures. This paper describes the development and testing of a vectorized, 3D migration program for the CYBER 205 using the Kosloff/Baysal method. The program can accept as many as 65,536 zero offset (stacked) traces.

Pyle, L. D.↗

Full wave modeling of radio-frequency beams in tokamaks in the electron cyclotron frequency range

Simulation of full wave, without paraxial approximation, high-resolution solution of wave equations in frequency domain in the electron cyclotron resonance (ECR) frequency range for realistic Tokamak plasma parameters became possible by using recently formulated hybrid iterative algorithm [Svidzinski et al., Phys. Plasmas 25, 082509 (2018)] for numerically solving discretized wave equations. This approach combines time evolution and iterative relaxation techniques into iteration cycles. This algorithm is implemented in 2D code FullWave, solving wave equations in Tokamaks in cold and hot plasma models, and it has been tested in 3D full wave iterative RF beams simulation tool, which is presently being developed to model 3D ECRH RF beams in fusion devices using dynamic grid adaptation. The results of 2D full wave modeling, assuming specified toroidal mode number, of ECRH RF beams in DIII-D plasma, performed in the cold and hot plasma models for outboard and top launch scenarios using FullWave are presented. Nonlocal hot plasma response model, based on accurate numerical solution of linearized Vlasov equation, is used to model beam propagation and absorption in the 2nd electron cyclotron harmonic region. Demonstration of capability of the hybrid iterative algorithm to model ECRH RF beams in 3D is made by simulating a substantial part of realistic beam in DIII-D, launched from outboard side of the machine. All relevant physics of RF beam propagation, most of which is not captured in paraxial approximation, such as beam's divergence, interference between the X and O modes in the beam, X-O mode conversion, beam splitting into the X and O mode beams, transformation of beam's cross section, and absorption at the 2nd electron cyclotron harmonic, is captured in the simulations. A numerical technique to find an optimal beam polarization at the launcher to launch a nearly pure X or O mode beam in plasma is developed and tested.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Radiation boundary conditions for wave-like equations

In the numerical computation of hyperbolic equations it is not practical to use infinite domains; instead, the domain is truncated with an artificial boundary. In the present study, a sequence of radiating boundary conditions is constructed for wave-like equations. It is proved that as the artificial boundary is moved to infinity the solution approaches the solution of the infinite domain as O(r exp -m-1/2) for the m-th boundary condition. Numerical experiments with problems in jet acoustics verify the practical nature of the boundary conditions.

Bayliss, A.↗

High frequency green function for aerodynamic noise in moving media. I - General theory. II - Noise from a spreading jet

It is shown how a high frequency analysis can be made for general problems involving flow-generated noise. In the parallel shear flow problem treated by Balsa (1976) and Goldstein (1982), the equation governing sound propagation in the moving medium could be transformed into a wave equation for a stationary medium with an inhomogeneous index of refraction. It is noted that the procedure of Avila and Keller (1963) was then used to construct a high frequency Green function. This procedure involves matching a solution valid in an inner region around the point source to an outer, ray-acoustics solution. This same procedure is used here to construct the Green function for a source in an arbitrary mean flow. In view of the fact that there is no restriction to parallel flow, the governing equations cannot be transformed into a wave equation; the analysis therefore proceeds from the equations of motion themselves.

Durbin, P. A.↗

Reflection and interference of electromagnetic waves in inhomogeneous media

Solutions were obtained of the wave equation for a plane horizontally polarized electro-magnetic wave incident on a semi infinite two dimensional inhomogeneous medium. Two problems were considered: An inhomogeneous half space, and an inhomogeneous layer of arbitrary thickness. Solutions of the wave equation were obtained in terms of Hankel functions with complex arguments. Numerical calculations were made of the reflection coefficient R at the interface of the homogeneous medium. The startling results show that the reflection coefficient for a complex dielectric constant with gradient, can be less than that of the same medium with zero gradient.

Geiger, F. E.↗

Quantum principles and free particles

The quantum principles that establish the energy levels and degeneracies needed to evaluate the partition functions are explored. The uncertainty principle is associated with the dual wave-particle nature of the model used to describe quantized gas particles. The Schroedinger wave equation is presented as a generalization of Maxwell's wave equation; the former applies to all particles while the Maxwell equation applies to the special case of photon particles. The size of the quantum cell in phase space and the representation of momentum as a space derivative operator follow from the uncertainty principle. A consequence of this is that steady-state problems that are space-time dependent for the classical model become only space dependent for the quantum model and are often easier to solve. The partition function is derived for quantized free particles and, at normal conditions, the result is the same as that given by the classical phase integral. The quantum corrections that occur at very low temperatures or high densities are derived. These corrections for the Einstein-Bose gas qualitatively describe the condensation effects that occur in liquid helium, but are unimportant for most practical purposes otherwise. However, the corrections for the Fermi-Dirac gas are important because they quantitatively describe the behavior of high-density conduction electron gases in metals and explain the zero point energy and low specific heat exhibited in this case.

Source record↗

Small-amplitude disturbances in turbomachine flows with swirl

A method is proposed for determining the response of the steady swirling potential flow past a blade row in an axial turbomachine to the most general type of 'nonacoustic' incident disturbance. The method is based on Goldstein's decomposition of the disturbance velocity and only requires solving a linear inhomogeneous wave equation. It is believed that numerical solutions to this wave equation can be obtained more efficiently and will be more accurate than corresponding solutions to the linearized compressible Euler equations.

Wundrow, David W.↗

Progress in modeling atmospheric propagation of sonic booms

The improved simulation of sonic boom propagation through the real atmosphere requires greater understanding of how the transient acoustic pulses popularly termed sonic booms are affected by humidity and turbulence. A realistic atmosphere is invariably somewhat turbulent, and may be characterized by an ambient fluid velocity v and sound speed c that vary from point to point. The absolute humidity will also vary from point to point, although possibly not as irregularly. What is ideally desired is a relatively simple scheme for predicting the probable spreads in key sonic boom signature parameters. Such parameters could be peak amplitudes, rise times, or gross quantities obtainable by signal processing that correlate well with annoyance or damage potential. The practical desire for the prediction scheme is that it require a relatively small amount of knowledge, possibly of a statistical nature, concerning the atmosphere along, the propagation path from the aircraft to the ground. The impact of such a scheme, if developed, implemented, and verified, would be that it would give the persons who make planning decisions a tool for assessing the magnitude of environmental problems that might result from any given overflight or sequence of overflights. The technical approach that has been followed by the author and some of his colleagues is to formulate a hierarchy of simple approximate models based on fundamental physical principles and then to test these models against existing data. For propagation of sonic booms and of other types of acoustic pulses in nonturbulent model atmospheres, there exists a basic overall theoretical model that has evolved as an outgrowth of geometrical acoustics. This theoretical model depicts the sound as propagating within ray tubes in a manner analogous to sound in a waveguide of slowly varying cross-section. Propagation along the ray tube is quasi-one-dimensional, and a wave equation for unidirectional wave propagation is used. A nonlinear term is added to this equation to account for nonlinear steepening, and the formulation has been carried through to allow for spatially varying sound speed, ambient density, and ambient wind velocities. The model intrinsically neglects diffraction, so it cannot take into account what has previously been mentioned in the literature as possibly important mechanisms for turbulence-related distortion. The model as originally developed could predict an idealized N-waveform which often agrees with data in terms of peak amplitude and overall positive phase duration. It is possible, moreover, to develop simple methods based on the physics of relaxation processes for incorporating molecular relaxation into the quasi-one-dimensional model of nonlinear propagation along ray tubes.

Pierce, Allan D.↗

The adiabatic semiclassical perturbation theory for vibrationally inelastic scattering. I - Collinear calculations. II - Three-dimensional treatment

A semiclassical approximation to treat vibrationally inelastic scattering is developed. The vibrational basis set used is adiabatic with respect to a reference potential which is chosen to be as close as possible to the true potential and also gives easily obtainable solutions to the vibrational wave equation. The radial wave functions are obtained using the WKB approximation, and the coupled Schroedinger equations are solved by a first-order perturbation method to yield a phase shift matrix which is exponentiated to give the full scattering matrix. Results were obtained for all the cases computed by Secrest and Johnson and by Clark and Dickinson, and the agreement is better than 10% for half of the cross-sections and rarely off by more than a factor of 2.

Cross, R. J., Jr.↗