Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “high-frequency 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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

Electron precipitation in the postmidnight sector of the auroral zones

Measurements of the angular distributions and energy spectra of electron intensities within the energy range 50 eV to 15 keV with electrostatic analyzer arrays on board the low-altitude satellite Injun 5 are reported for the postmidnight sector of the auroral zones during the high-intensity events accompanying magnetic substorms. Precipitation features on closed terrestrial field lines well equatorward of the trapping boundary for electrons with energies greater than 45 keV are examined. Precipitation of low-energy electron intensities was characterized by isotropy for all pitch angles outside the atmospheric backscatter cone. The region of electron precipitation observed is associated with the diffuse aurora and with pulsating aurora in the postmidnight sector. Similar variations of the energetic electron intensities with energies above 45 keV were observed in the regions of fluctuating energy fluxes of low-energy electrons associated with auroral luminosity. The increases of energetic electron intensities were not coincident with those of the principal energy fluxes into the atmosphere, except when the average electron energy for the energy fluxes was unusually high, i.e., in the 10-keV range. Precipitation of electron intensities within these energy ranges is consistent with strong pitch angle diffusion of electron intensities near or at the magnetic equator by high-frequency wave turbulence, the effectiveness of which is modulated by perturbations attributable to micropulsations.

Frank, L. A.↗
NASA NTRS1985

Direct numerical simulations of a reacting mixing layer with chemical heat release

In order to study the coupling between chemical heat release and fluid dynamics, direct numerical simulations of a chemically reacting mixing layer with heat release are performed. The fully compressible equations as well as an approximate set of equations that is asymptotically valid for low-Mach-number flows are treated. These latter equations have the computational advantage that high-frequency acoustic waves have been filtered out, allowing much larger time steps to be taken in the numerical solution procedure. A detailed derivation of these equations along with an outline of the numerical solution technique is given. Simulation results indicate that the rate of chemical product formed, the thickness of the mixing layer, and the amount of mass entrained into the layer all decrease with increasing rates of heat release.

Mcmurtry, P. A.↗
NASA NTRS1986

On uniformly valid high-frequency far-field asymptotic solutions of the Helmholtz equation

An asymptotic, large wave number approximation for the Helmholtz equation is derived. The theory is an extension of the geometric acoustic theory, and provides corrections to that theory in the form of multiplicative functions which satisfy parabolic equations. A simple example is used both to illustrate failure of the geometric theory for large propagation distances, and to show the improvement obtained by use of the new theory.

Mcaninch, G. L.↗
NASA NTRS1972

High-frequency sound waves to eliminate a horizon in the mixmaster universe.

From the linear wave equation for small-amplitude sound waves in a curved space-time, there is derived a geodesiclike differential equation for sound rays to describe the motion of wave packets. These equations are applied in the generic, nonrotating, homogeneous closed-model universe (the 'mixmaster universe,' Bianchi type IX). As for light rays described by Doroshkevich and Novikov (DN), these sound rays can circumnavigate the universe near the singularity to remove particle horizons only for a small class of these models and in special directions. Although these results parallel those of DN, different Hamiltonian methods are used for treating the Einstein equations.

Chitre, D. M.↗
NASA NTRS1973

Geodesic synchrotron radiation

This paper presents the results and methods of computing the high-frequency radiation emitted by freely falling particles moving in circular geodesic orbits in a spherically symmetric gravitational field. The high-frequency radiation, to which the methods of this paper apply, is the principal part of radiated energy only in the case of a particle moving in a highly relativistic, and therefore unstable, circular geodesic. The geodesic synchrotron radiation emitted in this case shows excitation of high-frequency harmonics and a narrow angular distribution. A Green's-function solution of the scalar wave equation is obtained using WKB methods. For application to relativistic circular orbits, a parabolic WKB approximation is required and yields solutions in terms of parabolic cylinder functions.

Breuer, R. A.↗
NASA NTRS1986

Nonlinear analysis of a relativistic beam-plasma cyclotron instability

A self-consistent set of nonlinear and relativistic wave-particle equations are derived for a magnetized beam-plasma system interacting with electromagnetic cyclotron waves. In particular, the high-frequency cyclotron mode interacting with a streaming and gyrating electron beam within a background plasma is considered in some detail. This interaction mode may possibly find application as a high-power source of coherent short-wavelength radiation for laboratory devices. The background plasma, although passive, plays a central role in this mechanism by modifying the dielectric properties in which the magnetized electron beam propagates. For a particular choice of the transverse beam velocity (i.e., the speed of light divided by the relativistic mass factor), the interaction frequency equals the nonrelativistic electron cyclotron frequency times the relativistic mass factor. For this choice of transverse beam velocity the detrimental effects of a longitudinal beam velocity spread is virtually removed. Power conversion efficiencies in excess of 18 percent are both analytically calculated and obtained through numerical simulations of the wave-particle equations. The quality of the electron beam, degree of energy and pitch angle spread, and its effect on the beam-plasma cyclotron instability is studied.

Sprangle, P.↗
NASA NTRS1974

Gravitational Stokes parameters

The electromagnetic and gravitational Stokes parameters are defined in the general theory of relativity. The general-relativistic equation of radiative transfer for polarized radiation is then derived in terms of the Stokes parameters for both high-frequency electromagnetic and gravitational waves. The concept of Stokes parameters is generalized for the most general class of metric theories of gravity, where six (instead of two) independent states of polarization are present.

Anile, A. M.↗
NASA NTRS1981

Functional methods for waves in random media

Some basic ideas in functional methods for waves in random media are illustrated through a simple random differential equation. These methods are then generalized to solve certain random parabolic equations via an exponential representation given by the Feynman-Kac formula. It is shown that these functional methods are applicable to a number of problems in random wave propagation. They include the forward-scattering approximation in Gaussian white-noise media; the solution of the optical beam propagation problem by a phase-integral method; the high-frequency scattering by bounded random media, and a derivation of approximate moment equations from the functional integral representation.

Chow, P. L.↗
NASA NTRS1981

Functional methods for waves in random media

Some basic ideas in functional methods for waves in random media are illustrated through a simple random differential equation. These methods are then generalized to solve certain random parabolic equations via an exponential representation given by the Feynman-Kac formula. It is shown that these functional methods are applicable to a number of problems in random wave propagation. They include the forward-scattering approximation in Gaussian white-noise media; the solution of the optical beam propagation problem by a phase-integral method; the high-frequency scattering by bounded random media; and a derivation of approximate moment equations from the functional integral representation.

Chow, P. L.↗
NASA NTRS1973

Propagation of plane waves in flow through a variable area duct

Duct geometry is an important factor that influences the transmission of sound in a duct. Nonconstant area produces variations in steady flow quantities through the duct, which can cause reflections of acoustic disturbances and the creation of standing wave patterns, and which can attenuate and disperse propagating waves. To study the problem, a second-order-accurate numerical method has been developed. The method of characteristics is used to solve the acoustic equations with the assumption of quasi-one-dimensional flow in the duct. Numerical results are presented over a broad range of frequencies. Low-frequency results are compared with approximate solutions for compact nozzles and high-frequency calculations are compared with results from the short-wave theory. The numerical results are also compared with an analytical solution of the acoustic equations for exponential horn sections with zero mean flow to further establish the validity of the numerical method.

King, L. S.↗
NASA NTRS1973

Acoustic propagation in ducts with varying cross sections

The method of multiple scales is used to derive the equations that describe the spatial and temporal variation of the amplitudes and phases of a wave packet propagating in slowly varying hard-walled or lined ducts. The analysis is carried out for rectangular as well as circular ducts. These equations are statements of the conservation of energy. For large admittance or high-frequency modes, an approximate expression is obtained for the attenuation. This expression shows that all possible acoustic modes are attenuating. The results also show that decreasing the cross sectional area can lead to elimination of some of the acoustic modes.

Nayfeh, A. H.↗
NASA NTRS1986

Taylor-Goertler instabilities of Tollmien-Schlichting waves and other flows governed by the interactive boundary-layer equations

The Taylor-Goertler vortex instability equations are formulated for steady and unsteady interacting boundary-layer flows. The effective Goertler number is shown to be a function of the wall shape in the boundary layer and the possibility of both steady and unsteady Taylor-Goertler modes exists. As an example the steady flow in a symmetrically constricted channel is considered and it is shown that unstable Goertler vortices exist before the boundary layers at the wall develop the Goldstein singularity discussed by Smith and Daniels (1981). As an example of an unsteady spatially varying basic state, it is considered the instability of high-frequency large-amplitude two- and three-dimensional Tollmien-Schlichting waves in a curved channel. It is shown that they are unstable in the first 'Stokes-layer stage' of the hierarchy of nonlinear states discussed by Smith and Burggraf (1985). This instability of Tollmien-Schlichting waves in an internal flow can occur in the presence of either convex or concave curvature. Some discussion of this instability in external flows is given.

Hall, Philip↗
DOE OSTI2025

Cold freeze out of superheavy dark matter and Hubble tension

We present a unified framework, the "X miracle", in which dark matter consists of superheavy, nonthermal X particles whose relic abundance is determined not by the conventional weak-scale, semi-relativistic ("hot") freeze-out of WIMPs, but by annihilation or decay occurring within the smallest and earliest gravitationally bound objects. Unlike thermal WIMPs, which decouple at velocities of order 0.3c with relic abundance ρ∞ set by weak-scale interactions, X particles are produced nonthermally with an initial overabun dance ρ ini >> ρ ∞ . They become nonrelativistic extremely early, redshift to ultra-cold velocities, allowing collapse into compact bound structures characterized by a novel quantum gravitational scale, r X = 4$\hbar$ 2 $/Gm^3_X$ = 10 −13 m $\hbar$$/m_Xc$, much larger than the Compton wavelength. The framework predicts a particle mass of 10 12 GeV and an enhanced cross section of 10 −21 m 3 /s. Overlapping particle wavefunctions in these compact structures drive annihilation or decay into additional radiation, leading to a "cold" freeze-out that converts most of ρ ini into radiation while leaving a relic density ρ ∞ . Solutions to the Boltzmann equation indicate that an extreme ("big") depletion, with only one particle in a billion surviving, yields an additional radiation contribution $ΔN_{eff}$ ≈ 0.4, which could help alleviate the Hubble tension. For particles of 10 12 GeV, the scenario predicts a dark coupling constant α X = 0.09 that is responsible for an instanton-induced decay process, consistent with current UHECR bounds. Early collapse at 10 −6 s may release binding energy as high-frequency (100kHz) gravitational waves or ultralight GUT-scale axions (10 −9 eV). Superheavy sterile neutrinos provide a natural particle realization, linking dark matter to neutrino mass and baryogenesis. If gravitationally produced, this framework favors high-scale inflation and effi cient reheating. The "X miracle" thus demonstrates that dark matter need not be weak-scale: gravitational dynamics can control freeze-out and evolution, producing multi-messenger observational signatures in UHECRs, axions, gravitational waves, and small-scale structures.

Xu, Zhijie Jay [Pacific Northwest National Laborat↗
DOE OSTI2025

The shear Alfvén continuum of quasisymmetric stellarators

The shear Alfvén wave (SAW) continuum plays a critical role in the stability of energetic particle-driven Alfvén eigenmodes. We develop a theoretical framework to analyze the SAW continuum in three-dimensional (3-D) quasisymmetric magnetic fields, focusing on its implications for stellarator design. By employing a near-axis model and degenerate perturbation theory, the continuum equation is solved, highlighting unique features in 3-D configurations, such as the interactions between spectral gaps. Numerical examples validate the theory, demonstrating the impact of flux-surface shaping and quasisymmetric field properties on continuum structure. The results provide insights into optimizing stellarator configurations to minimize resonance-driven losses of energetic particles. This work establishes a basis for incorporating Alfvénic stability considerations into the stellarator design process, demonstrated through optimization of a quasihelical configuration to avoid high-frequency spectral gaps.

fusion plasma↗
DOE OSTI2025

A Linear-Complexity Tensor Butterfly Algorithm for Compressing High-Dimensional Oscillatory Integral Operators

This paper presents a multilevel tensor compression algorithm called tensor butterfly algorithm for efficiently representing large-scale and high-dimensional oscillatory integral operators, including Green's functions for wave equations and integral transforms such as Radon transforms and Fourier transforms. The proposed algorithm leverages a tensor extension of the so-called complementary low-rank property of existing matrix butterfly algorithms. The algorithm partitions the discretized integral operator tensor into subtensors of multiple levels and factorizes each subtensor at the middle level as a Tucker-type interpolative decomposition, whose factor matrices are formed in a multilevel fashion. For a d-dimensional (d > 1) integral operator discretized into a 2d-mode tensor with n2d entries, the overall CPU time and memory requirement scale as O(nd), in stark contrast to the O(nd log n) complexity of existing matrix algorithms such as matrix butterfly algorithms and fast Fourier transforms (FFTs), where n is the number of points per direction. When comparing with other tensor algorithms such as quantized tensor train (QTT), the proposed algorithm also shows superior CPU and memory performance for tensor contraction. Remarkably, the tensor butterfly algorithm can efficiently model high-frequency Green's function interactions between two unit cubes, each spanning 512 wavelengths per direction, which represents problems of scale over 512× larger than that existing butterfly algorithms can handle, with the same amount of computation resources. On the other hand, for a problem representing 64 wavelengths per direction, which is the largest size existing algebraic matrix algorithms can handle, our tensor butterfly algorithm exhibits 200x speedups and 30× memory reduction compared with existing ones. Moreover, the tensor butterfly algorithm also permits O(nd)-complexity FFTs and Radon transforms up to d = 6 dimensions.

Kielstra, P Michael↗
NASA NTRS1976

An experimental and computational investigation of the flow field about a transonic airfoil in supercritical flow with turbulent boundary-layer separation

A combined experimental and computational research program for testing and guiding turbulence modeling within regions of separation induced by shock waves incident on turbulent boundary layers is described. Specifically, studies are made of the separated flow over the rear portion of an 18%-thick circular-arc airfoil at zero angle of attack in high Reynolds number supercritical flow. The measurements include distributions of surface static pressure and local skin friction. The instruments employed include high-frequency response pressure cells and a large array of surface hot-wire skin-friction gages. Computations at the experimental flow conditions are made using time-dependent solutions of ensemble-averaged Navier-Stokes equations, plus additional equations for the turbulence modeling.

Rubesin, M. W.↗