Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part I: Model Formulation

Here, this paper formulates a new particle-in-cell method for the Vlasov–Maxwell system. Under the Lorenz gauge condition, Maxwell’s equations for the electromagnetic fields can be written as a collection of scalar and vector wave equations. The use of potentials for the fields motivates the adoption of a Hamiltonian formulation for particles that employs the generalized (conjugate) momentum. A notable advantage offered by the Hamiltonian formulation is the elimination of time derivatives in the Lorenz gauge formulation that are required by the standard Newton–Lorentz treatment of the particles. This allows the fields to retain the full time-accuracy guaranteed by the field solver. The resulting updates for particles require only knowledge of the fields and their spatial derivatives. An analytical method for constructing these spatial derivatives is presented that exploits the underlying integral solution used in the field solver for the wave equations. Moreover, these derivatives are demonstrated to converge at the same rate as the fields in both time and space. The Method of Lines Transpose field solver we consider in this work is globally first-order accurate in time and high-order accurate in space (e.g., fourth- and fifth-order) and belongs to a larger class of methods which are unconditionally stable, can address geometry, and leverage $\mathcal {O}(N)$ fast summation methods for efficiency. We demonstrate the method on several well-established benchmark problems on bounded domains, including a plasma sheath as well as a relativistic particle beam. The efficacy of the proposed formulation is established by comparing with a second-order accurate finite-difference time-domain method that employs a leapfrog time advance for particles and a charge conserving map suitable for bounded domains. The new method shows mesh-independent numerical heating properties even in cases where the plasma Debye length is smaller than the grid spacing. This is an important feature of the new method for problems defined on bounded domains, because it permits the use of coarser grids in space in the representation of the fields. Such a capability has significant implications for the simulation of plasmas in bounded domains with complex geometry, where the ratio between the largest and smallest cells can vary significantly. The use of high-order spatial approximations in the new method also means that fewer grid points are required in order to achieve a fixed accuracy. Our results also suggest that the new method can be used with fewer simulation particles per cell compared to the benchmark explicit method, which permits further computational savings.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1983

Analytical study of the twin-jet shielding

The development of the analytical model of twin-jet shielding is summarized. The models consist of a point noise source impinging on a cylinder of heated flow in which the temperature and flow velocity are uniform cross the cross section. In the formulation of the model, the wave equations are written for the regions outside the flow and within the flow cylinder. The solutions to the wave equations are matched at the jet boundary under the conditions of continuity of pressure and continuity of the vortex sheet. The solution reduces to an indefinite integral involving Bessel functions. The integral is solved approximately using a saddle point method.

Gerhold, C. H.↗
NASA NTRS1990

Physical interpretation and application of principles of ultrasonic nondestructive evaluation of high-performance materials

An ultrasonic measurement system employed in the experimental interrogation of the anisotropic properties (through the measurement of the elastic stiffness constants) of the uniaxial graphite-epoxy composites is presented. The continuing effort for the development of improved visualization techniques for physical parameters is discussed. The background is set for the understanding and visualization of the relationship between the phase and energy/group velocity for propagation in high-performance anisotropic materials by investigating the general requirements imposed by the classical wave equation. The consequences are considered when the physical parameters of the anisotropic material are inserted into the classical wave equation by a linear elastic model. The relationship is described between the phase velocity and the energy/group velocity three dimensional surfaces through graphical techniques.

Miller, James G.↗
DOE OSTI2016

Photoacoustic effect generated by moving optical sources: Motion in one dimension

Although the photoacoustic effect is typically generated by pulsed or amplitude modulated optical beams, it is clear from examination of the wave equation for pressure that motion of an optical source in space will result in the production of sound as well. Here, the properties of the photoacoustic effect generated by moving sources in one dimension are investigated. The cases of a moving Gaussian beam, an oscillating delta function source, and an accelerating Gaussian optical sources are reported. The salient feature of one-dimensional sources in the linear acoustic limit is that the amplitude of the beam increases in time without bound.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1973

Aerodynamic noise emission from turbulent shear layers.

The Phillips (1960) convected wave equation is employed in this paper to study aerodynamic noise emission processes in subsonic and supersonic shear layers. The wave equation in three spatial dimensions is first reduced to an ordinary differential equation by Fourier transformation and then solved via the WKBJ method. Three typical solutions are required for discussions in this paper. The current results are different from the classical conclusions. The effects of refraction, convection, Mach-number dependence and temperature dependence of turbulent noise emission are analyzed in the light of solutions to the Phillips equation.

Pao, S. P.↗
NASA NTRS1989

The Helmholtz-Kirchhoff integral relation as a framework for developing algorithms for sound propagation through inhomogeneous moving media

Transient sound propagation in an inhomogeneous moving medium is considered. For circumstances in which the medium is slowly varying over distances of a wavelength but possibly varying substantially over the propagation distance, a derivation is given of a new wave equation, which implicitly allows for diffraction and scattering and which also is consistent with earlier geometrical acoustics formulations. This wave equation is used as a starting point to derive a version of the Helmholtz-Kirchhoff integral relation that applies to inhomogeneous moving medium. It is suggested that a good approximation to the Green's function that appears in this relation is that derived from geometrical acoustics, the approximation becoming progressively better the shorter the distance between surfaces in the computation. It should also be at least as good as conventional ray acoustics, but can account for diffraction effects, such as at caustics.

Pierce, Allan D.↗
NASA NTRS1990

Graphics applications utilizing parallel processing

The results are presented of research conducted to develop a parallel graphic application algorithm to depict the numerical solution of the 1-D wave equation, the vibrating string. The research was conducted on a Flexible Flex/32 multiprocessor and a Sequent Balance 21000 multiprocessor. The wave equation is implemented using the finite difference method. The synchronization issues that arose from the parallel implementation and the strategies used to alleviate the effects of the synchronization overhead are discussed.

Rice, John R.↗
NASA NTRS1995

Pseudoclassical description of the Dirac Oscillator

In this paper we discuss the Dirac Oscillator wave equation in terms of pseudoclassical language, using Grassmann variables to describe the internal degrees of freedom of the oscillator. Regarding the original wave equation as a classical constraint, we use the theory of constrained systems, to develop a reparameterization invariant lagrangian, which is the pseudoclassical equivalent of the quantum case. The consistency of the Hamiltonian formalism and the quantization procedure are also analyzed.

Delsolmesa, Antonio↗
NASA NTRS1991

A numerical method for unsteady aerodynamics via acoustics

Formal solutions to the wave equation may be conveniently described within the framework of generalized function theory. A generalized function theory is used to yield a formulation and formal solution of a wave equation describing oscillation of a flat plate from which a numerical method may be derived.

Hodge, Steve↗
DOE OSTI2016

Erratum: Nonlinear Dirac equation solitary waves in external fields [Phys. Rev. E 86, 046602 (2012)]

In Sec. IV of our original paper, we assumed a particular conservation law Eq. (4.6), which was true in the absence of external potentials, to derive some particular potentials for which we obtained solutions to the nonlinear Dirac equation (NLDE). Because the conservation law of Eq. (4.6) for the component T 11 of the energy-momentum tensor is not true in the presence of these external potentials, the solutions we found do not satisfy the NLDEs in the presence of these potentials. Thus all the equations from Eq. (4.6) through Eq. (4.44) are not correct, since the exact solutions that followed in that section presumed Eq. (4.6) was true. Also Eqs. (A3)–(A5) are a restatement of Eq. (4.6) and also are not correct. These latter equations are also not used in Sec. V and beyond. The rest of our original paper (starting with Sec. V) was not concerned with exact solutions, rather it was concerned with how the exact solitary-wave solutions to the NLDE in the absence of an external potential responded to being in the presence of various external potentials. This Erratum corrects this mistake.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1991

Critical problems of computational aeroacoustics

The aeroacoustics of rigid boundaries is discussed. Lighthill gave a formulation of this problem in which he showed that the sources of the acoustic field were quadrupole in nature. We have preferred a different formulation of the problem, in which the quadrupoles are sources for a nonlinear wave equation, as opposed to the linear one used by Lighthill. This is given here in a figure which also gives further details of a solution procedure for the Euler equations appropriate for the aeroelastic problem and motivated by the analysis of Crow. In accordance with our formulation, we contend that an algorithm for accurate solutions of an inhomogeneous nonlinear wave equation is of prime importance.

Rogers, Joel C. W.↗
DOE OSTI2016

Integrated surface/subsurface permafrost thermal hydrology: Model formulation and proof-of-concept simulations

The need to understand potential climate impacts and feedbacks in Arctic regions has prompted recent interest in modeling of permafrost dynamics in a warming climate. A new fine-scale integrated surface/subsurface thermal hydrology modeling capability is described and demonstrated in proof-of-concept simulations. The new modeling capability combines a surface energy balance model with recently developed three-dimensional subsurface thermal hydrology models and new models for nonisothermal surface water flows and snow distribution in the microtopography. Surface water flows are modeled using the diffusion wave equation extended to include energy transport and phase change of ponded water. Variation of snow depth in the microtopography, physically the result of wind scour, is also modeled heuristically with a diffusion wave equation. The multiple surface and subsurface processes are implemented by leveraging highly parallel community software. Fully integrated thermal hydrology simulations on the tilted open book catchment, an important test case for integrated surface/subsurface flow modeling, are presented. Fine-scale 100-year projections of the integrated permafrost thermal hydrological system on an ice wedge polygon at Barrow Alaska in a warming climate are also presented. Finally, these simulations demonstrate the feasibility of microtopography-resolving, process-rich simulations as a tool to help understand possible future evolution of the carbon-rich Arctic tundra in a warming climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Quasi-periodic travelling gravity–capillary waves

We present a numerical study of spatially quasi-periodic travelling waves on the surface of an ideal fluid of infinite depth. This is a generalization of the classic Wilton ripple problem to the case when the ratio of wavenumbers satisfying the dispersion relation is irrational. We propose a conformal mapping formulation of the water wave equations that employs a quasi-periodic variant of the Hilbert transform to compute the normal velocity of the fluid from its velocity potential on the free surface. We develop a Fourier pseudo-spectral discretization of the travelling water wave equations in which one-dimensional quasi-periodic functions are represented by two-dimensional periodic functions on the torus. This leads to an overdetermined nonlinear least-squares problem that we solve using a variant of the Levenberg–Marquardt method. We investigate various properties of quasi-periodic travelling waves, including Fourier resonances, time evolution in conformal space on the torus, asymmetric wave crests, capillary wave patterns that change from one gravity wave trough to the next without repeating and the dependence of wave speed and surface tension on the amplitude parameters that describe a two-parameter family of waves.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Reducing Frequency Bias of Fourier Neural Operators in 3D Seismic Wavefield Simulations Through Multistage Training

The recent development of neural operator (NeurOp) learning for solutions to the elastic wave equation shows promising results and provides the basis for fast large-scale simulations for different seismological applications. In this article, we use the Fourier neural operator (FNO) model to directly solve the 3D Helmholtz wave equation for fast seismic ground-motion simulations on different frequencies and show the frequency bias of the FNO model, that is, it learns the lower frequencies better comparing to the higher frequencies. To reduce the frequency bias, we adopt the multistage FNO training, that is, after training a stage 1 FNO model for estimating the ground motion, we use a second FNO model as the stage 2 to learn from the residual, which greatly reduced the errors on the higher frequencies. By adopting this multistage training, the FNO models show reduced biases on higher frequencies, which enhanced the overall results of the ground-motion simulations. Thus the multistage training FNO improves the accuracy and realism of the ground-motion simulations.

earthquakes↗
DOE OSTI2019

Moving photoacoustic sources: Acoustic waveforms in one, two, and three dimensions and application to trace gas detection

The photoacoustic effect is governed by a wave equation with a source term proportional to the time derivative of the optical heat deposition per unit volume and time. Although the typical configuration for generation of the photoacoustic effect makes use of pulsed or amplitude modulated optical beams, the form of the source term in the wave equation indicates that a continuous optical source moving in an absorbing medium is capable of sound generation as well. Here, the properties of simple sources moving in one, two, and three space dimensions are reviewed. The salient feature of sources moving in one-dimension at sound speed is that the amplitude of the acoustic wave increases with time without bound according to linear acoustics. Finally, two schemes, one in the time-domain and the other in the frequency-domain, that take advantage of this principle for increasing the sensitivity of trace gas detection are reviewed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1973

Acoustic wave propagation in an axisymmetric swirling jet.

An analysis has been developed to study the acoustic wave propagation in an axisymmetric swirling subsonic jet flow. The governing convected wave equation derived in the spherical coordinates includes mean shears, shear gradients and pressure gradients. The directivity patterns for various spinning and non-spinning modes due to the influence of the mean jet swirl were obtained by numerically integrating the governing wave equation. The mean flow field used in the computation was that obtained semiempirically for subsonic swirling turbulent jet and is completely specified once the degree of swirl is known. The dependence of sound directivity on jet Mach number, swirl ratio and frequency are discussed.

Yu, J. C.↗
NASA NTRS1989

Oxygen acoustic solitary waves in a magnetized plasma

Ion-acoustic solitary waves in a magnetized plasma containing an arbitrary mixture of H(+) and O(+) ions are studied. A nonlinear wave equation has been derived from the Poisson-Vlasov equations, including a uniform magnetic field and dissipation due to reflected electrons. When dissipation is ignored, the equation has soliton solutions associated with both oxygen and hydrogen acoustic modes, which can be either rarefactive or compressive depending on the ion concentrations and the electron/ion temperature ratio and, more weakly, on the bulk drifts of the species. If electron reflection is included, the solitary wave can be intensified. Under somewhat restrictive conditions the oxygen solitary wave is rarefactive and propagates with a velocity comparable to that observed by the Viking satellite. The three-dimensional solitons obey a relation of scales parallel to the magnetic field and in the transverse direction. Computer simulations of one-dimensional versions of the nonlinear wave equation are presented.

Qian, S.↗