Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulated waves”

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 1,117 records · Page 62

Numerical simulation of the generation of electrostatic turbulence in the magnetotail

A two-dimensional plasma model is used to investigate the development of electrostatic turbulence in a magnetized plasma from plasma instabilities. The simulation consists of following the motion of 100,000 ions in their self-consistent electrostatic field. The electrons are treated as a constant neutralizing background. The instabilities modeled are driven by a ring-type velocity distribution and by interpenetrating ion beams in a time-variable magnetic field. Instability growth times are of the order of an ion gyroperiod in the case of the ring distribution and of the order of an ion plasma period in the case of the beam simulation. Maximum potential differences generated are of the order of the ion kinetic energies. These simulations demonstrate the cascade of wave energy to long wavelengths, thus showing the E x B turbulence can be generated from plasma microinstabilities. After the free energy feeding, the instabilities are exhausted, and wave energy at wavelengths less than an ion gyrodiameter decays quickly to equilibrium levels, while longer wavelength modes persist for much longer times. In one model with a time dependent, but spatially uniform, magnetic field the electric field energy at long wavelengths appeared to increase as a result of the increase of the magnetic field.

Swift, D. W.↗

Physics of auroral arc formation; Proceedings of the Chapman Conference on Formation of Auroral Arcs, Fairbanks, AK, July 21-25, 1980

The results of satellite, balloon-borne, and ground-based observations of auroral phenomena are discussed in terms of the morphology of auroral arcs and the behavior and characteristics of auroral electrons and ions. Further attention is given to auroral electric fields and field aligned currents, and to the development of models of auroral potential structures and energization of auroral particles. A simulation of space plasma phenomena is presented, along with numerical simulations of auroral potential structures and related phenomena, including the V-potential double layers and auroral arc deformations, and a simulation of auroral arcs. Plasma waves observed on auroral field lines and in laboratory conditions are reported, and theoretical studies of waves and turbulence in auroral plasmas are described.

Akasofu, S.-I.↗

Transverse acceleration of ions on auroral field lines

This paper examines transverse ion heating on auroral field lines associated with current-driven electrostatic ion cyclotron waves theoretically and by numerical simulations. The auroral plasma is assumed to consist of drifting electrons and stationary hydrogen and oxygen ions. Depending on the ratio of the electron drift speed to the thermal speed and the ratio of hydrogen to oxygen concentrations, preferential heating of either hydrogen or oxygen ions can take place. It is found that unless the oxygen ions are a minority species, oxygen transverse heating generally exceeds that of hydrogen ions. Theory and numerical simulations are in good agreement.

Ashour-Abdalla, M.↗

Applications of finite element and wave envelope element approximations to turbofan engine noise radiation including flight effects

The problem of acoustic radiation from turbofan engine inlets in flow has not lent itself fully to analysis by numerical means because of the large domains and high frequencies involved. The current work has extended the use of finite elements and wave envelope elements, elements which simulate decay and wavelike behaviour in their interpolation functions, from the no-flow case in which they have been proven, to cases incorporating mean flow. By employing an irrotational mean flow assumption, the acoustics problem has been posed in an axisymmetric formulation in terms of acoustic velocity potential, thus minimizing computer solution storage requirements. The results obtained from the numerical procedures agree well with known analytical solutions, static experimental jet engines inflow data, and also with flight test results.

Parrett, A. V.↗

Heating of thermal helium in the equatorial magnetosphere - A simulation study

The heating of He(+) ions observed by the GEOS satellite at L of about seven starting from an initial anisotropic proton distribution is studied. A simulation of electromagnetic ion cyclotron waves in a plasma consisting of hot anisotropic protons, a dominant thermal hydrogen plasma, and a minority species of helium ions is utilized. It is shown that the thermal ions are heated to 100 times their initial temperature due to the growth of large-amplitude ion cyclotron waves. He(+) ions are heated mainly in the perpendicular direction, and they are heated more than cold H(+) ions. The He(+) ions are heated in a two-step process. First, they are set into oscillations by the growing wave until some of the ions reach a parallel velocity on the order of the resonant velocity; then, strong heating occurs.

Omura, Y.↗

The generation of low-frequency electrostatic waves in association with auroral kilometric radiation

Observations of the electron distribution in the source region of auroral kilometric radiation (AKR) indicate that it can have positive gradients with respect to both v-perpendicular and v-parallel. It is shown that this type of distribution is unstable to both the electron cyclotron maser instability, which is responsible for AKR, and the bump-in-tail instability which generates electron acoustic waves in the hiss band. Simulations indicate that the two instabilities do not develop independently but compete for the available free energy due to the reduction in the positive gradients of the distribution caused by the quasi-linear diffusion associated with each instability. The dominant instability depends on the ratio of the density of the energetic electrons n(E) to that of the background electrons n(b). For n(E)/n(b) greater than about 1, the maser instability dominates, and the bump-in-tail instability is suppressed, whereas the reverse is true if n(E)/n(b) is smaller than about 1. Comparison with observations indicates that probably n(E)/n(b) is greater than about 1 in the source region of AKR.

Winglee, R. M.↗

Cell averaging Chebyshev methods for hyperbolic problems

A cell averaging method for the Chebyshev approximations of first order hyperbolic equations in conservation form is described. Formulas are presented for transforming between pointwise data at the collocation points and cell averaged quantities, and vice-versa. This step, trivial for the finite difference and Fourier methods, is nontrivial for the global polynomials used in spectral methods. The cell averaging methods presented are proven stable for linear scalar hyperbolic equations and present numerical simulations of shock-density wave interaction using the new cell averaging Chebyshev methods.

Wei, Cai↗

Testing continuum descriptions of low-Mach-number shock structures

Numerical experiments have been performed on normal shock waves with Monte Carlo Direct Simulations (MCDS's) to investigate the validity of continuum theories at very low Mach numbers. Results from the Navier-Stokes and the Burnett equations are compared to MCDS's for both hard-sphere and Maxwell gases. It is found that the maximum-slope shock thicknesses are described equally well (within the MCDS computational scatter) by either of the continuum formulations for Mach numbers smaller than about 1.2. For Mach numbers greater that 1.2, the Burnett predictions are more accurate than the Navier-Stokes results. Temperature-density profile separations are best described by the Burnett equations for Mach numbers greater than about 1.3. At lower Mach numbers the MCDS scatter is too great to differentiate between the two continuum theories. For all Mach numbers above one, the shock shapes are more accurately described by the Burnett equations.

Pham-Van-diep, Gerald C.↗

Magnetospherically reflected whistlers as a source of plasmaspheric hiss

Ray-tracing simulations and estimates of whistler wave damping show that magnetospherically reflected whistlers can persist for about 100 s in a LF band (about 1 kHz). The combined contribution from whistler rays produced by a single lightning flash but entering the magnetosphere at different points form a continuous hisslike signal, as observed at a fixed point. Estimates indicate that the total whistler wave energy input into the magnetosphere from lightning discharges may maintain experimentally observed levels of magnetospheric hiss.

Draganov, A. B.↗

Experimental Investigation of Nozzle/Plume Aerodynamics at Hypersonic Speeds

The work performed by D. W. Bogdanoff and J.-L. Cambier during the period of 1 Feb. - 31 Oct. 1992 is presented. The following topics are discussed: (1) improvement in the operation of the facility; (2) the wedge model; (3) calibration of the new test section; (4) combustor model; (5) hydrogen fuel system for combustor model; (6) three inch calibration/development tunnel; (7) shock tunnel unsteady flow; (8) pulse detonation wave engine; (9) DCAF flow simulation; (10) high temperature shock layer simulation; and (11) the one dimensional Godunov CFD code.

Heinemann, K.↗

A time domain, weighted residual formulation of Maxwell's equations

A finite element model is developed and used to simulate two-dimensional electromagnetic wave propagation and scattering. The spatial discretization of the time-domain electrodynamic equations is accomplished by a Galerkin approach. The semi-discrete equations are solved explicitly using a second-order Runge-Kutta scheme. Both the electric and magnetic fields are discretized using a single grid, with the divergence-free conditions satisfied through a correction approach. Examples depicting the scattering of plane waves in 2D geometries are given to demonstrate the validity of the methodology.

Young, Jeffrey L.↗

Sea surface height variability during 1986-1988 in the tropical Pacific Ocean

A description of the sea surface height variabilities during 1986-1988 in the tropical Pacific Ocean is presented on the basis of Geosat altimeter data and a tropical Pacific Ocean GCM (OGCM) forced with observed winds from Florida State University. Both data sets are validated against in situ sea-level measurements at 12 selected sea-level stations in the equatorial Pacific and were found to be in good agreement with these measurements with correlation coefficients of 0.69 for Geosat and 0.71 for OGCM. The Geosat altimeter data are capable of describing sea surface height variabilities during 1986-1988 in the tropical Pacific, in particular, the seasonal cycle, the 1987 El Nino and 1988 La Nina episodes, and the December 1986 equatorial Kelvin wave event. The OGCM-simulated surface dynamic height is found to be sensitive to different surface wind products.

Chao, YI↗

Computational fluid mechanics

Two papers are included in this progress report. In the first, the compressible Navier-Stokes equations have been used to compute leading edge receptivity of boundary layers over parabolic cylinders. Natural receptivity at the leading edge was simulated and Tollmien-Schlichting waves were observed to develop in response to an acoustic disturbance, applied through the farfield boundary conditions. To facilitate comparison with previous work, all computations were carried out at a free stream Mach number of 0.3. The spatial and temporal behavior of the flowfields are calculated through the use of finite volume algorithms and Runge-Kutta integration. The results are dominated by strong decay of the Tollmien-Schlichting wave due to the presence of the mean flow favorable pressure gradient. The effects of numerical dissipation, forcing frequency, and nose radius are studied. The Strouhal number is shown to have the greatest effect on the unsteady results. In the second paper, a transition model for low-speed flows, previously developed by Young et al., which incorporates first-mode (Tollmien-Schlichting) disturbance information from linear stability theory has been extended to high-speed flow by incorporating the effects of second mode disturbances. The transition model is incorporated into a Reynolds-averaged Navier-Stokes solver with a one-equation turbulence model. Results using a variable turbulent Prandtl number approach demonstrate that the current model accurately reproduces available experimental data for first and second-mode dominated transitional flows. The performance of the present model shows significant improvement over previous transition modeling attempts.

Hassan, H. A.↗

Thermal expansivity, bulk modulus, and melting curve of H2O-ice VII to 20 GPa

Equation of state properties of ice VII and fluid H2O at high pressures and temperatures have been studied experimentally from 6 to 20 GPa and 300-700 K. The techniques involve direct measurements of the unit-cell volume of the solid using synchrotron X-ray diffraction with an externally heated diamond-anvil cell. The pressure dependencies of the volume and bulk modulus of ice VII at room temperature are in good agreement with previous synchrotron X-ray studies. The thermal expansivity was determined as a function of pressure and the results fit to a newly proposed phenomenological relation and to a Mie-Gruneisen equation of state formalism. The onset of melting of ice VII was determined directly by X-ray diffraction at a series of pressures and found to be in accord with previous volumetric determinations. Thermodynamic calculations based on the new data are performed to evaluate the range of validity of previously proposed equations of state for fluid water derived from static and shock-wave compression experiments and from simulations.

Fei, Yingwei↗

Boundary Conditions for Unsteady Compressible Flows

This paper explores solutions to the spherically symmetric Euler equations. Motivated by the work of Hagstrom and Hariharan and Geer and Pope, we modeled the effect of a pulsating sphere in a compressible medium. The literature available on this suggests that an accurate numerical solution requires artificial boundary conditions which simulate the propagation of nonlinear waves in open domains. Until recently, the boundary conditions available were in general linear and based on nonreflection. Exceptions to this are the nonlinear nonreflective conditions of Thompson, and the nonlinear reflective conditions of Hagstrom and Hariharan. The former are based on the rate of change of the incoming characteristics; the latter rely on asymptotic analysis and the method of characteristics and account for the coupling of incoming and outgoing characteristics. Furthermore, Hagstrom and Hariharan have shown that, in a test situation in which the flow would reach a steady state over a long time, Thompson's method could lead to an incorrect steady state. The current study considers periodic flows and includes all possible types and techniques of boundary conditions. The technique recommended by Hagstrom and Hariharan proved superior to all others considered and matched the results of asymptotic methods that are valid for low subsonic Mach numbers.

Hariharan, S. I.↗

Interdecadal Variations in the Alaska Gyre

Climatic dynamic topography variations in the Alaska gyre during the period 1968-1990 are described with an objective analysis of more than 12000 STD and XBT stations, and COADS wind stress data Interannual the dynamic height and SST variations were correlated and were consistent with recently described large-scale climatic shifts in the North Pacific. The gyre was centered more to the east, circulation appeared stronger, and SST was lower during the early to mid-1970s than during the 1980s. The Aleutian low (NP and PNA indices) intensified during the interim, but the response did not appear as a gyre spinup. Instead, the associated wind stress anomalies forced a slowly varying dynamic height anomaly across the eastern and northern part of the gyre through Ekman convergence, which had the effect of displacing the gyre's low somewhat to the WSW in the 1980s. The wind curl spectrum was white, and the slow oceanic response was modeled as stochastic-forced climate variability with a simple first-order Markov autoregression process. Forcing was assumed to be Ekman pumping of the pycnocline, and the damping coefficient was estimated from the data to be approx. 1 yr. A hindcast with observed winds gave estimated dynamic height patterns similar to those observed, with a canonical correlation of 0.79 at 99% confidence. This response was weak in the western half of the gyre, where slow baroclinic variability may have been influenced by long Rossby wave propagation. A simple autoregression simulation using artificial white noise forcing shows the evolution of decadal variations similar in nature to those observed. This result, along with the low frequency correlation between dynamic height and SST, suggests that the upper-ocean climatic variability in this region is primarily wind forced.

Lagerloef, Gary S.E.↗

Minnowbrook I: 1993 Workshop on End-Stage Boundary Layer Transition

This volume contains materials presented at the Minnowbrook I-1993 Workshop on End-Stage Boundary Layer Transition, held at the Syracuse University Minnowbrook Conference Center, New York, from August 15 to 18, 1993. This volume was previously published as a Syracuse University report edited by John E. LaGraff. The workshop organizers were John E. LaGraff (Syracuse University), Terry V. Jones (Oxford University), and J. Paul Gostelow (University of Technology, Sydney). The workshop focused on physical understanding of the late stages of transition from laminar to turbulent flows, with the specific goal of contributing to improving engineering design of turbomachinery and wing airfoils. The workshop participants included academic researchers from the United States and abroad, and representatives from the gas-turbine industry and U.S. government laboratories. To improve interaction and discussions among the participants, no formal papers were required. The physical mechanisms discussed were related to natural and bypass transition, wake-induced transition, effects of freestream turbulence, turbulent spots, hairpin vortices, nonlinear instabilities and breakdown, instability wave interactions, intermittency, turbulence, numerical simulation and modeling of transition, heat transfer in boundary-layer transition, transition in separated flows, laminarization, transition in turbomachinery compressors and turbines, hypersonic boundary-layer transition, and other related topics. This volume contains abstracts and copies of the viewgraphs presented, organized according to the workshop sessions. The workshop summary and the plenary discussion transcript clearly outline future research needs.

LaGraff, John E.↗

The Mock LISA Data Challenge Round 3: New and Improved Sources

The Mock LISA Data Challenges are a program to demonstrate and encourage the development of data-analysis capabilities for LISA. Each round of challenges consists of several data sets containing simulated instrument noise and gravitational waves from sources of undisclosed parameters. Participants are asked to analyze the data sets and report the maximum information they can infer about the source parameters. The challenges are being released in rounds of increasing complexity and realism. Challenge 3. currently in progress, brings new source classes, now including cosmic-string cusps and primordial stochastic backgrounds, and more realistic signal models for supermassive black-hole inspirals and galactic double white dwarf binaries.

Baker, John↗