Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear 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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

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

Increasing Accuracy in Computed Inviscid Boundary Conditions

A technique has been devised to increase the accuracy of computational simulations of flows of inviscid fluids by increasing the accuracy with which surface boundary conditions are represented. This technique is expected to be especially beneficial for computational aeroacoustics, wherein it enables proper accounting, not only for acoustic waves, but also for vorticity and entropy waves, at surfaces. Heretofore, inviscid nonlinear surface boundary conditions have been limited to third-order accuracy in time for stationary surfaces and to first-order accuracy in time for moving surfaces. For steady-state calculations, it may be possible to achieve higher accuracy in space, but high accuracy in time is needed for efficient simulation of multiscale unsteady flow phenomena. The present technique is the first surface treatment that provides the needed high accuracy through proper accounting of higher-order time derivatives. The present technique is founded on a method known in art as the Hermitian modified solution approximation (MESA) scheme. This is because high time accuracy at a surface depends upon, among other things, correction of the spatial cross-derivatives of flow variables, and many of these cross-derivatives are included explicitly on the computational grid in the MESA scheme. (Alternatively, a related method other than the MESA scheme could be used, as long as the method involves consistent application of the effects of the cross-derivatives.) While the mathematical derivation of the present technique is too lengthy and complex to fit within the space available for this article, the technique itself can be characterized in relatively simple terms: The technique involves correction of surface-normal spatial pressure derivatives at a boundary surface to satisfy the governing equations and the boundary conditions and thereby achieve arbitrarily high orders of time accuracy in special cases. The boundary conditions can now include a potentially infinite number of time derivatives of surface-normal velocity (consistent with no flow through the boundary) up to arbitrarily high order. The corrections for the first-order spatial derivatives of pressure are calculated by use of the first-order time derivative velocity. The corrected first-order spatial derivatives are used to calculate the second- order time derivatives of velocity, which, in turn, are used to calculate the corrections for the second-order pressure derivatives. The process as described is repeated, progressing through increasing orders of derivatives, until the desired accuracy is attained.

Dyson, Roger↗
NASA NTRS1996

A Simulated Spectrum of Convectively Generated Gravity Waves: Propagation from the Tropopause to the Mesopause and Effects on the Middle Atmosphere

This work evaluates the interaction of a simulated spectrum of convectively generated gravity waves with realistic middle atmosphere mean winds. The wave spectrum is derived from the nonlinear convection model described by Alexander et al. [1995] that simulated a two-dimensional midlatitude squall line. This spectrum becomes input to a linear ray tracing model for evaluation of wave propagation as a function of height through climatological background wind and buoyancy frequency profiles. The energy defined by the spectrum as a function of wavenumber and frequency is distributed spatially and temporally into wave packets for the purpose of estimating wave amplitudes at the lower boundary of the ray tracing model. A wavelet analysis provides an estimate of these wave packet widths in space and time. Without this redistribution of energies into wave packets the Fourier analysis alone inaccurately assumes the energy is evenly distributed throughout the storm model domain. The growth with height of wave amplitudes is derived from wave action flux conservation coupled to a convective instability saturation condition. Mean flow accelerations and wave energy dissipation profiles are derived from this analysis and compared to parameterized estimates of gravity wave forcing, providing a measure of the importance of the storm source to global gravity wave forcing. The results suggest that a single large convective storm system like the simulated squall line could provide a significant fraction of the zonal mean gravity wave forcing at some levels, particularly in the mesosphere. The vertical distributions of mean flow acceleration and energy dissipation do not much resemble the parameterized profiles in form because of the peculiarities of the spectral properties of the waves from the storm source. The ray tracing model developed herein provides a tool for examining the role of convectively generated waves in middle atmosphere physics.

Alexander, M. Joan↗
NASA NTRS1961

Constant Density Approximations for the Flow Behind Axisymmetric Shock Waves

The incompressible rotational flow equations are used to obtain solutions for the flow behind axisymmetric shock waves with conic longitudinal sections. The nonlinear part of the term due to rotation is retained in the analysis. Numerical results for standoff distance and stagnation point velocity gradient are presented for the case in which the shock wave is a paraboloid, a sphere, or an oblate or prolate ellipsoid. A similarity parameter is proposed which correlates approximately the flow behind geometrically similar shock waves at different free-stream conditions.

Munson, Albert G.↗
NASA NTRS1985

The evolution of Ertel's potential vorticity during stratospheric sudden warmings

In the winter stratosphere of the Northern Hemisphere, the disruption of the westerly vortex and associated warming of polar latitudes is a well known phenomenon. It has become apparent that some important dynamical processes in the stratosphere are highly nonlinear and are best thought of locally rather than in terms of the interaction between the zonal-mean flow and harmonic waves around latitude circles. The importance of nonlinear processes was suggested by McIntyre and Palmer (1983, 1984) who used isentropic maps of Ertel's Potential Vorticity to show that during disturbed episodes material lines may become strongly and irreversibly deformed in certain places. They adopted the term planetary wave breaking to describe this process. Isentropic maps of Q are used to follow the evolution of a Canadian warming in November - December 1981 and a particularly strong warming in January 1982. The advection of Q over large distances on isentropic surfaces was a striking feature of the flow during each event. This could be identified because of our ability to follow the movement of material lines due to the approximate conservation of Q over several days. The advection of Q was a nonlinear process because its changing distribution affected the advecting wind field. The Canadian warming did not lead to a permanent change in the structure of the westerly vortex, as defined by the coarse-grain field of Q, whereas the January event was accompanied by a substantial loss of resolved Q which was never fully recovered.

Fairlie, T. D. A.↗
NASA NTRS1995

High amplitude waves in the expanding solar wind plasma

We simulated the 1-D nonlinear time-evolution of high-amplitude Alfven, slow and fast magnetoacoustic waves in the solar wind propagating outward at different angles to the mean magnetic (spiral) field, using the expanding box model. The simulation results for Alfven waves and fast magnetoacustic waves fit the observational constraints in the solar wind best, showing decreasing trends for energies and other rms-quantities due to expansion and the appearance of inward propagating waves as minor species in the wind. Inward propagating waves are generated by reflection of Alfven waves propagating at large angles to the magnetic field or they coincide with the occurrence of compressible fluctuations. In our simulations, fast and slow magnetoacoustic waves seem to have a level in the density-fluctuations which is too high when we compare with the observations. Furthermore, the evolution of energies for slow magnetoacoustic waves differs strongly from the evolution of fluctuation energies in situ.

Schmidt, J. M.↗
NASA NTRS1994

Pondermotive versus mirror force in creation of the filamentary cavities in auroral plasma

Recently rocket observations on spikelets of lower-hybrid waves along with strong density cavities and transversely heated ions were reported. The observed thin filamentary cavities oriented along the magnetic field in the auroral plasma have density depletions up to several tens of percent. These observations have been interpreted in terms of a theory for lower-hybrid wave condensation and collapse. The modulational instability leading to the wave consensation of the lower-hybrid waves yields only weak density perturbations, which cannot explain the above strong density depletions. The wave collapse theory is based on the nonlinear pondermotive force in a homogeneous ambient plasma and the density depletion is determined by the balance between the wave pressure (pondermotive force) and the plasma pressure. In the auroral plasma, the balance is achieved in a time tau(sub wc) equal to or less than 1 ms. It is shown here that the mirror force, acting on the transversely heated ions at a relatively long time scale, is an effective mechanism for creating the strong plasma cavities. We suggest that the process of wave condensation, through the pondermotive force causing generation of short wavelength waves from relatively long wavelength waves, is a dominant process until the former waves evolve and become effective in the transverse heating of ions. As soon as this happens, mirror force on ions becomes an important factor in the creation of the density cavities, which may further trap and enhance the waves. Results from a model of cavity formation by transverse ion heating show that the observed depletions in the density cavities can be produced by the heating rates determined by the observed wave amplitudes near the lower-hybrid frequency. It is found that the creation of a strong density cavity takes a few minutes.

Singh, Nagendra↗
NASA NTRS1990

Steepening of kinetic magnetosonic waves into shocklets - Simulations and consequences for planetary shocks and comets

The generation and the nonlinear evolution of oblique low-frequency electromagnetic (kinetic magnetosonic) waves which were observed upstream of planetary bow shocks and at the Giacobini-Zinner comet, and referred to as shocklets, were investigated using an electromagnetic hybrid code. The observations show that the waves, which have a sinusoidal form when their amplitude is small, become steepened and linearly polarized as they grow in amplitude. The results of simulations show the original small-amplitude elliptically polarized wave grows and steepens, so that its polarization changes and becomes somewhat linear. The steepening process is associated with the coherent generation of a broad spectrum of waves on the magnetosonic whistler branch, which propagate at various phase and group velocities. It is shown that the presence of shocklets upstream of a planetary bow shock can modify its local structure by changing the solar wind Mach number and temperature, or by colliding with the shock.

Omidi, N.↗
NASA NTRS2000

Applications of the Space-Time Conservation Element and Solution Element (CE/SE) Method to Computational Aeroacoustic Benchmark Problems

The Internal Propagation problems, Fan Noise problem, and Turbomachinery Noise problems are solved using the space-time conservation element and solution element (CE/SE) method. The problems in internal propagation problems address the propagation of sound waves through a nozzle. Both the nonlinear and linear quasi 1D Euler equations are solved. Numerical solutions are presented and compared with the analytical solution. The fan noise problem concerns the effect of the sweep angle on the acoustic field generated by the interaction of a convected gust with a cascade of 3D flat plates. A parallel version of the 3D CE/SE Euler solver is developed and employed to obtain numerical solutions for a family of swept flat plates. Numerical solutions for sweep angles of 0, 5, 10, and 15 deg are presented. The turbomachinery problems describe the interaction of a 2D vortical gust with a cascade of flat-plate airfoils with/without a downstream moving grid. The 2D nonlinear Euler Equations are solved and the converged numerical solutions are presented and compared with the corresponding analytical solution. All the comparisons demonstrate that the CE/SE method is capable of solving aeroacoustic problems with/without shock waves in a simple and efficient manner. Furthermore, the simple non-reflecting boundary condition used in the CE/SE method which is not based on the characteristic theory works very well in 1D, 2D and 3D problems.

Wang, Xiao-Yen↗
NASA NTRS1988

Nonlinear evolution of the Kelvin-Helmholtz instability in the high-latitude ionosphere

The first numerical simulations of the nonlinear evolution of the electrostatic Kelvin-Helmholtz (K-H) instability with ionospheric Perderson conductivity coupling are presented. It is found that the K-H instability develops in a distinctly different manner in the nonlinear regime with Pedersen coupling than without it. Pedersen coupling effects, in conjunction with a neutral wind and density gradient, are shown to result in an increased time scale for K-H instability wave growth, to inhibit K-H vortex formation, to lead to nonlinear structures which can be described as 'breaking waves', and to generate, in the nonlinear regime, small-scale turbulence by means of secondary instabilities growing on primary waves. The spatial power spectra of the electrostatic potential and density fluctuations are computed, and differences with and without Pedersen effects are reported.

Keskinen, M. J.↗
NASA NTRS1966

Quasi-linear theory of plasma waves.

Plasma wave quasi-linear theory, noting nonlinear effect treatment via distribution function, electric field, Fourier component decay, etc

DISTRIBUTION FUNCTION↗
NASA NTRS1995

Application of essentially nonoscillatory methods to aeroacoustic flow problems

A finite-difference essentially nonoscillatory (ENO) method has been applied to several of the problems prescribed for the workshop sponsored jointly by the Institute for Computer Applications in Science and Engineering and by NASA Langley Research Center entitled 'Benchmark Problems in Computational Aeroacoustics'. The workshop focused on computational challenges specific to aeroacoustics. Among these are long-distance propagation of a short-wavelength disturbance, propagation of small-amplitude disturbances, and nonreflective boundary conditions. The shock capturing-capability inherent to the ENO method effectively eliminates oscillations near shock waves without the need to add and tune dissipation or filter terms. The method-of-lines approach allows the temporal and spatial operators to be chosen separately in accordance with the demands of a particular problem. The ENO method was robust and accurate for all problems in which the propagating wave was resolved with 8 or more points per wavelength. The finite-wave-model boundary condition, a local nonlinear acoustic boundary condition, performed well for the one-dimensional problems. The buffer-domain approach performed well for the two-dimensional test problem. The amplitude of nonphysical reflections were less than 1 percent of the exiting wave's amplitude.

Atkins, Harold L.↗
NASA NTRS1989

Extended inertial range phenomenology of magnetohydrodynamic turbulence

A phenomenological treatment of the inertial range of isotropic statistically steady magnetohydrodynamic turbulence is presented, extending the theory of Kraichnan (1965). The role of Alfven wave propagation is treated on equal footing with nonlinear convection, leading to a simple generalization of the relations between the times characteristic of wave propagation, convection, energy transfer, and decay of triple correlations. The theory leads to a closed-form steady inertial range spectral law that reduces to the Kraichnan and Kolmogorov laws in appropriate limits. The Kraichnan constant is found to be related in a simple way to the Kolmogorov constant; for typical values of the latter constant, the former has values in the range 1.22-1.87. Estimates of the time scale associated with spectral transfer of energy also emerge from the new approach, generalizing previously presented 'golden rules' for relating the spectral transfer time scale to the Alfven and eddy-turnover time scales.

Matthaeus, William H.↗
NASA NTRS1993

Spiral density waves resonantly excited by a rapidly rotating bar

Recent observations at millimeter wavelengths have revealed bar-spiral structures in the central regions of several nearby spiral galaxies. Moreover, new infrared data and even COBE data seem to support the idea that there is a bar in the central regions of the Milky Way. All of these new findings have rekindled the interest in bar-driven spirals again. First, the mechanism that is responsible for exciting the spiral density waves near the resonances is reviewed, and then a theory to account for the highly nonlinear behavior of the gas associated with the waves is presented. The new results have again demonstrated that the 3 kpc arm phenomenon can be reproduced by a rapidly rotating bar in the galactic center. Similarly, in terms of the resonance excitation mechanism, the recent high-resolution radio observations of the bar-spiral structure in the central regions of nearby galaxies can be explained.

Yuan, Chi↗
NASA NTRS1991

Alfven shock trains

The Cohen-Kulsrud-Burgers equation (CKB) is used to consider the nonlinear evolution of resistive, quasi-parallel Alfven waves subject to a long-wavelength, plane-polarized, monochromatic instability. The instability saturates by nonlinear steepening, which proceeds until the periodic waveform develops an interior scale length comparable to the dissipation length; a fast or an intermediate shock then forms. The result is a periodic train of Alfven shocks of one or the other type. For propagation strictly parallel to the magnetic field, there will be two shocks per instability wavelength. Numerical integration of the time-dependent CKB equation shows that an initial, small-amplitude growing wave asymptotes to a stable, periodic stationary wave whose analytic solution specifies how the type of shock embedded in the shock train, and the amplitude and speed of the shock train, depend on the strength and phase of the instability. Waveforms observed upstream of the earth's bowshock and cometary shocks resemble those calculated here.

Malkov, M. A.↗
NASA NTRS1996

Mixing and Transition Control Studied

Considerable progress in understanding nonlinear phenomena in both unbounded and wallbounded shear flow transition has been made through the use of a combination of high- Reynolds-number asymptotic and numerical methods. The objective of this continuing work is to fully understand the nonlinear dynamics so that ultimately (1) an effective means of mixing and transition control can be developed and (2) the source terms in the aeroacoustic noise problem can be modeled more accurately. Two important aspects of the work are that (1) the disturbances evolve from strictly linear instability waves on weakly nonparallel mean flows so that the proper upstream conditions are applied in the nonlinear or wave-interaction streamwise region and (2) the asymptotic formulations lead to parabolic problems so that the question of proper out-flow boundary conditions--still a research issue for direct numerical simulations of convectively unstable shear flows--does not arise. Composite expansion techniques are used to obtain solutions that account for both mean-flow-evolution and nonlinear effects. A previously derived theory for the amplitude evolution of a two-dimensional instability wave in an incompressible mixing layer (which is in quantitative agreement with available experimental data for the first nonlinear saturation stage for a plane-jet shear layer, a circular-jet shear layer, and a mixing layer behind a splitter plate) have been extended to include a wave-interaction stage with a three-dimensional subharmonic. The ultimate wave interaction effects can either give rise to explosive growth or an equilibrium solution, both of which are intimately associated with the nonlinear self-interaction of the three dimensional component. The extended theory is being evaluated numerically. In contrast to the mixing-layer situation, earlier comparisons of theoretical predictions based on asymptotic methods and experiments in wall-bounded shear-flow transition have been somewhat lacking in one aspect or another. The current work strongly suggests that the main weakness is the underlying asymptotic representation of the linear "part" of the problem and not the explicit modeling of the nonlinear/wave-interaction effects. Consequently, the long-wave-length/high-Reynolds-number asymptotic limit for the Blasius boundary-layer stability problem was reexamined, and a new dispersion relationship for the instability waves that is uniformly valid for both the upper- and lower branch regions to the required order of approximation was obtained. A comparison with numerical results, obtained by solving the Orr-Sommerfeld stability problem, shows that the asymptotic formula provides surprisingly good results, even for values of the frequency parameter usually encountered in experimental investigations. This is particularly evident in the dynamically important upper-branch region, where much of the nonlinear interactions in transition experiments are believed to take place. The result is important in that it can be used to greatly improve the accuracy of weakly nonlinear critical-layer-based theories, and a consistent nonlinear theory is currently under evaluation.

Source record↗
NASA NTRS1955

From linear mechanics to nonlinear mechanics

Consideration is given to the techniques used in telecommunication where a nonlinear system (the modulator) results in a linear transposition of a signal. It is then shown that a similar method permits linearization of electromechanical devices or nonlinear mechanical devices. A sweep function plays the same role as the carrier wave in radio-electricity. The linearizations of certain nonlinear functionals are presented.

RESEARCH TECHNIQUE, MATHEMATICS↗
NASA NTRS1985

Upconversion of lower-hybrid waves by gyrating ion beams in a plasma

A gyrating ion beam, with ring shape distribution in velocity space, supports negative energy modes near the harmonics of beam gyrofrequency. An investigation of the nonlinear interaction of high-frequency lower-hybrid waves with the negative energy beam cyclotron mode is made. A nonlinear dispersion relation is derived for the coupled modes. It is shown that a gyrating ion beam frequency upconverts the lower-hybrid waves separated by harmonics of the beam gyrofrequency. This process of upconversion may thus provide a tool for the diagnostics of the gyrating ion beams. The process may also be important in the saturation of lower-hybrid waves observed in the auroral zone S3-3 satellite data.

Sharma, O. P.↗