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 37 records · Page 2

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.↗

On the receptivity and nonparallel stability of traveling disturbances in rotating-disk flow

The generation and evolution of small amplitude long wavelength traveling disturbances in rotating-disk flow are discussed. The steady rotational speed of the disk is perturbed so as to introduce high-frequency oscillations in the flow field. Secondly, surface imperfections on the disk such as roughness elements are introduced. The interaction of these two disturbances will generate the instability waves whose evolution is governed by parabolic partial differential equations which are solved numerically. It is found that, for the class of disturbances considered here (wavelength on the order of Reynolds number), eigensolutions exist which decay or grow algebraically in the radial direction. However, these solutions grow only for frequencies larger than 4.58 times the steady rotational speed of the disk. The computed receptivity coefficient shows that there is an optimum size of roughness for which these modes are preferentially excited.

Balakumar, P.↗

A new interpretation of plasma-line overshoot phenomena

The temporal evolution of Langmuir waves excited by high-power, high-frequency (HF) radio waves in the ionosphere is studied theoretically. This study is motivated by past observations made with the 450 MHz radar at Arecibo Observatory in Puerto Rico. Two kinds of nonlinear damping to the parametric decay instability are considered in the derivation of the rate equation for the spectral intensity of enhanced Langmuir waves. They are Langmuir wave cascading caused by nonlinear Landau damping and cross-field electron diffusion. The first damping process leads to the saturation of individual unstable Langmuir wave. The second process, which results from the incoherent scattering of electron orbits by the total excited Langmuir waves, yields anomalous damping that applies to each Langmuir wave. Consequently, Langmuir waves with smaller growth rates will be suppressed by those with larger growth rates. Such a mode competition process may cause the overshoot of the HF-enhanced plasma line observed with the Arecibo 430 MHz radar. Favorable agreement is obtained between theory and the Arecibo observations.

Kuo, S. P.↗

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↗

The possible role of high-frequency waves in heating solar coronal loops

We investigate the role of high-frequency waves in the heating of solar active region coronal loops. We assume a uniform background magnetic field, and we introduce a density stratification in a direction perpendicular to this field. We focus on ion compressive viscosity as the damping mechanism of the waves. We incorporate viscosity self-consistently into the equations, and we derive a dispersion relation by adopting a slab model, where the density inside the slab is greater than that outside. Such a configuration supports two types of modes: surface waves and trapped body waves. In order to determine under what conditions these waves may contribute to the heating of active regions, we solve our dispersion relation for a range of densities, temperatures, magnetic field strengths, density ratios, wavevector magnitudes, wavevector ratios, and slab widths. We find that surface waves exhibit very small damping, but body waves can potentially damp at rates needed to balance radiative losses. However, the required frequencies of these body waves are very high. For example, the wave frequency must be at least 5.0/s for a slab density of 10(exp 9,5)/cc, a slab temperature of 10(exp 6,5) K, a field strength of 100 G, and a density ratio of 5. For a slab density of 10(exp 10)/cc, this frequency increases to 8.8/s. Although these frequencies are very high, there in no observational evidence to rule out their existence, and they may be generated both below the corona and at magnetic reconnection sites in the corona. However, we do find that, for slab densities of 10(exp 10)/cc or less, the dissipation of high-frequency waves will be insufficient to balance the radiative losses if the magnetic field strength exceeds roughly 200 G. Because the magnetic field is known to exceed 200 G in many active region loops, particularly low-lying loops and loops emanating from sunspots, it is unlikely that high-frequency waves can provide sufficient heating in these regions.

Porter, Lisa J.↗

Effect of artificial viscosity on three-dimensional flow solutions

Artificial viscosity is added either implicity or explicitly in practically every numerical scheme for suppressing spurious oscillations in the solution of fluid-dynamics equations. In the present central-difference scheme, artificial viscosity is added explicitly for suppressing high-frequency oscillations and achieving good convergence properties. The amount of artificial viscosity added is controlled through the use of preselected coefficients. In the standard scheme, scalar coefficients based on the spectral radii of the Jacobian of the convective fluxes are used. However, this can add too much viscosity to the slower waves. Hence, the use of matrix-valued coefficients, which give appropriate viscosity for each wave component, is suggested. With the matrix-valued coefficients, the central-difference scheme produces more accurate solutions on a given grid, particularly in the vicinity of shocks and boundary layers, while still maintaining good convergence properties.

Turkel, Eli↗

Accuracy Study of the Space-Time CE/SE Method for Computational Aeroacoustics Problems Involving Shock Waves

The space-time conservation element and solution element(CE/SE) method is used to study the sound-shock interaction problem. The order of accuracy of numerical schemes is investigated. The linear model problem.govemed by the 1-D scalar convection equation, sound-shock interaction problem governed by the 1-D Euler equations, and the 1-D shock-tube problem which involves moving shock waves and contact surfaces are solved to investigate the order of accuracy of numerical schemes. It is concluded that the accuracy of the CE/SE numerical scheme with designed 2nd-order accuracy becomes 1st order when a moving shock wave exists. However, the absolute error in the CE/SE solution downstream of the shock wave is on the same order as that obtained using a fourth-order accurate essentially nonoscillatory (ENO) scheme. No special techniques are used for either high-frequency low-amplitude waves or shock waves.

Wang, Xiao Yen↗

Computational Study of Near-limit Propagation of Detonation in Hydrogen-air Mixtures

A computational investigation of the near-limit propagation of detonation in lean and rich hydrogen-air mixtures is presented. The calculations were carried out over an equivalence ratio range of 0.4 to 5.0, pressures ranging from 0.2 bar to 1.0 bar and ambient initial temperature. The computations involved solution of the one-dimensional Euler equations with detailed finite-rate chemistry. The numerical method is based on a second-order spatially accurate total-variation-diminishing (TVD) scheme, and a point implicit, first-order-accurate, time marching algorithm. The hydrogen-air combustion was modeled with a 9-species, 19-step reaction mechanism. A multi-level, dynamically adaptive grid was utilized in order to resolve the structure of the detonation. The results of the computations indicate that when hydrogen concentrations are reduced below certain levels, the detonation wave switches from a high-frequency, low amplitude oscillation mode to a low frequency mode exhibiting large fluctuations in the detonation wave speed; that is, a 'galloping' propagation mode is established.

Yungster, S.↗

Direct Numerical Simulation of Transition in a Swept-Wing Boundary Layer

Direct numerical simulation (DNS) is performed to examine laminar to turbulent transition due to high-frequency secondary instability of stationary crossflow vortices in a subsonic swept-wing boundary layer for a realistic natural-laminar-flow airfoil configuration. The secondary instability is introduced via inflow forcing derived from a two-dimensional, partial-differential-equation based eigenvalue computation; and the mode selected for forcing corresponds to the most amplified secondary instability mode which, in this case, derives a majority of its growth from energy production mechanisms associated with the wall-normal shear of the stationary basic state. Both the growth of the secondary instability wave and the resulting onset of laminar-turbulent transition are captured within the DNS computations. The growth of the secondary instability wave in the DNS solution compares well with linear secondary instability theory when the amplitude is small; the linear growth is followed by a region of reduced growth resulting from nonlinear effects before an explosive onset of laminar breakdown to turbulence. The peak fluctuations are concentrated near the boundary layer edge during the initial stage of transition, but rapidly propagates towards the surface during the process of laminar breakdown. Both time-averaged statistics and flow visualization based on the DNS reveal a sawtooth transition pattern that is analogous to previously documented surface flow visualizations of transition due to stationary crossflow instability. The memory of the stationary crossflow vortex is found to persist through the transition zone and well beyond the location of the maximum skin friction.

Duan, Lian↗

Structure and Stability of One-Dimensional Detonations in Ethylene-Air Mixtures

The propagation of one-dimensional detonations in ethylene-air mixtures is investigated numerically by solving the one-dimensional Euler equations with detailed finite-rate chemistry. The numerical method is based on a second-order spatially accurate total-variation-diminishing scheme and a point implicit, first-order-accurate, time marching algorithm. The ethylene-air combustion is modeled with a 20-species, 36-step reaction mechanism. A multi-level, dynamically adaptive grid is utilized, in order to resolve the structure of the detonation. Parametric studies over an equivalence ratio range of 0.5 less than phi less than 3 for different initial pressures and degrees of detonation overdrive demonstrate that the detonation is unstable for low degrees of overdrive, but the dynamics of wave propagation varies with fuel-air equivalence ratio. For equivalence ratios less than approximately 1.2 the detonation exhibits a short-period oscillatory mode, characterized by high-frequency, low-amplitude waves. Richer mixtures (phi greater than 1.2) exhibit a low-frequency mode that includes large fluctuations in the detonation wave speed; that is, a galloping propagation mode is established. At high degrees of overdrive, stable detonation wave propagation is obtained. A modified McVey-Toong short-period wave-interaction theory is in excellent agreement with the numerical simulations.

Yungster, S.↗

On the Behavior of Three-dimensional Wave Packets in Viscously Spreading Mixing Layers

We consider analytically the evolution of a three-dimensional wave packet generated by an impulsive source in a mixing layer. The base flow is assumed to be spreading due to viscous diffusion. The analysis is restricted to small disturbances (linearized theory). A suitable high-frequency ansatz is used to describe the packet; the key elements of this description are a complex phase and a wave action density. It is found that the product of this density and an infinitesimal material volume convecting at the local group velocity is not conserved: there is a continuous interaction between the base flow and the wave action. This interaction is determined by suitable mode-weighted averages of the second and fourth derivatives of the base-flow velocity profile. Although there is some tendency for the dominant wave number in the packet to shift from the most unstable value toward the neutral value, this shift is quite moderate. In practice, wave packets do not become locally neutral in a diverging base flow (as do instability modes), therefore, they are expected to grow more suddenly than pure instability modes and do not develop critical layers. The group velocity is complex; the full significance of this is realized by analytically continuing the equations for the phase and wave action into a complex domain. The implications of this analytic continuation are discussed vis-a-vis the secondary instabilities of the packet: very small scale perturbations on the phase can grow very rapidly initially, but saturate later because most of the energy in these perturbations is convected away by the group velocity. This remark, as well as the one regarding critical layers, has consequences for the nonlinear theories.

Balsa, Thomas F.↗

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.↗

Accretion disk oscillations - A local analysis in a disk of finite thickness

Two types of oscillations are observed to occur in dwarf novae: 'coherent' and 'quasi-periodic' oscillations. These may be associated with the pulsation of the white dwarf or the accretion disk components of the dwarf nova. Here a local (short-wavelength) analysis is utilized to study the oscillation of a self-consistent, two-dimensional model of an accretion disk. The linearized equations describing adiabatic, inviscid, nonaxisymmetric oscillations are used to derive a fifth-order algebraic equation for the (complex) pulsation frequency of the disk. The solutions of this equation for various values of the wavevector k reveal that the disk is capable of supporting (1) a pair of high-frequency acoustic modes (p-modes); (2) a pair of intermediate-frequency modes which may share the characteristics of internal gravity waves (g-modes) and inertial waves; and (3) a mode associated with a dynamical instability (purely imaginary frequency). The role played by the shear in determining the stability or instability of these modes is also considered. Finally, the global oscillation frequencies of the disk are discussed.

Carroll, B. W.↗

A new formulation of electromagnetic wave scattering using an on-surface radiation boundary condition approach

A new formulation of electromagnetic wave scattering by convex, two-dimensional conducting bodies is reported. This formulation, called the on-surface radiation condition (OSRC) approach, is based upon an expansion of the radiation condition applied directly on the surface of a scatterer. It is now shown that application of a suitable radiation condition directly on the surface of a convex conducting scatterer can lead to substantial simplification of the frequency-domain integral equation for the scattered field, which is reduced to just a line integral. For the transverse magnetic case, the integrand is known explicitly. For the transverse electric case, the integrand can be easily constructed by solving an ordinary differential equation around the scatterer surface contour. Examples are provided which show that OSRC yields computed near and far fields which approach the exact results for canonical shapes such as the circular cylinder, square cylinder, and strip. Electrical sizes for the examples are ka = 5 and ka = 10. The new OSRC formulation of scattering may present a useful alternative to present integral equation and uniform high-frequency approaches for convex cylinders larger than ka = 1. Structures with edges or corners can also be analyzed, although more work is needed to incorporate the physics of singular currents at these discontinuities. Convex dielectric structures can also be treated using OSRC.

Kriegsmann, Gregory A.↗

Numerical simulation of shock-induced combustion/detonation in a premixed H2-air mixture using Navier-Stokes equations

A numerical study was conducted to address the structural stability of the oblique detonation wave. An attempt was made to isolate the structural instability predicted by Buckmaster's (1990) linear stability analysis. For this study a detailed viscous flowfield with finite-rate chemistry was computed past a 20-deg wedge. The finite-rate chemistry was modeled by using a seven-species and seven-reaction model. Instability modes were isolated by analyzing the Fourier power spectrum of the H2O mass fraction at selected sample stations. No unstable temporal modes with dominant frequency were isolated. It was concluded that the instability which is predicted by the Buckmaster work was a high-frequency low-amplitude phenomenon. For engineering purposes, the oblique detonation wave is a stable phenomenon as long as sufficient levels of overdrive are present.

Singh, D. J.↗

High-order essentially non-oscillatory methods for computational aeroacoustics

The desire to obtain acoustic information from the numerical solution of a nonlinear system of equations is a demanding proposition for a computational algorithm. High-order accuracy is required for the propagation of high-frequency, low-amplitude waves. In addition, it is desirable to highly resolve discontinuities that can develop in the solutions of the Euler or Navier-Stokes equations. The class of essentially non-oscillatory (ENO) shock-capturing schemes has been designed to have both of these properties. The dual capacity of ENO schemes for high-order accuracy and non-oscillatory shock-capturing is achieved through the use of adaptive stenciling, which makes these schemes highly nonlinear. These schemes are briefly described and referenced herein. A fourth-order algorithm is then applied to the solution of an acoustic wave in a quasi-one-dimensional converging-diverging nozzle.

Casper, Jay↗

Equatorial wave activity derived from fluctuations in observed convection

The spectrum of equatorial wave activity propagating vertically into the stratosphere is calculated from high-resolution imagery of the global convective pattern. Synoptic Global Cloud Imagery (GCI), constructed from six satellites simultaneously observing the earth, is used to diabatically force the linearized primitive equations. Having resolution of 0.5 deg and 3 h, that imagery captures the dominant scales of organized convection, including several harmonics of the diurnal cycle. Its global coverage with high space-time resolution allows the GCI to represent heating variability and dynamical behavior excited by it over a wide range of scales. The dynamical response above the heating is evaluated globally in terms of a space-time spectrum of Hough modes, one which includes planetary-scale Kelvin waves, Rossby waves, and gravity waves down to the resolution of the GCI. The geopotential response, which is indicative of temperature fluctuations observed by satellite, is very red in frequency. Therefore, planetary-scale waves with periods longer than two days dominate the spectrum of geopotential, while high-frequency gravity waves make a comparatively small contribution. Some 80% of the geopotential variance is accounted for by the Kelvin and gravest-symmetric Rossby modes, while the Rossby-gravity mode is comparatively weak. In horizontal eddy motion, the excited wave spectrum is still dominated by planetary-scale components. However, meridional wind fluctuations associated with the Rossby-gravity mode have variance comparable to that of zonal wind fluctuations associated with the Kelvin mode, even though the Rossby-gravity mode is nearly invisible in the geopotential response. Estimates of tropospheric heating lead to amplitudes and propagation characteristics that are broadly consistent with satellite and radiosonde observations of wave activity in the lower stratosphere. The space-time spectrum of EP flux is significantly whiter than the response in either geopotential or motion. Gravity waves of small scale and high frequency carry a large fraction of the upward flux. Although it dominates eastward variance of geopotential and motion, the Kelvin mode carries only about 50% of the eastward EP flux at phase speeds of 20-40 m/s and only 35% of the total eastward flux transmitted to the stratosphere. The remainder is carried by the gravity wave spectrum, which carries nearly all of the westward flux at phase speeds greater than 20 m/s. The gravity wave spectrum also contributes significantly at phase speeds of 10-20 m/s, where only 25% of the flux is accounted for by zonal wavenumbers less than 20. The broad nature of the gravity wave spectrum suggests its absorption at critical levels will be distributed over a deep layer of the middle atmosphere.

Bergman, John W.↗

Using High-Order Methods on Lower-Order Geometries

The desire to obtain acoustic information from the numerical solution of a nonlinear system of equations is a demanding proposition for a computational algorithm. High-order accuracy is required for the propagation of high-frequency, low-amplitude waves. The accuracy of an algorithm can be compromised by low-order errors that naturally occur in the solution of a particular problem. Such errors arise from two sources: the presence of discontinuities in the flow field or because the geometry on which the problem is defined is not everywhere smooth to the order of the scheme. The performance of high-order accurate essentially non-oscillatory (ENO) schemes on piecewise smooth solutions is well documented. Herein, the performance of these methods on smooth solutions defined on piecewise smooth geometries is investigated. The propagation of sound in a quasi-one-dimensional nozzle is considered as a test case. Some of the issues involved in the extension to two spatial dimensions are discussed.

Casper, Jay↗