Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear wave propagation”

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 307 records · Page 17

Convective amplification of Type 1 irregularities in the equatorial electrojet

Wave propagation and refraction of Type 1 irregularities in the equatorial electrojet were investigated. Quantitative calculation of wave refraction in a model electrojet showed that the direction of wave refraction must change sign at one altitude. Waves propagating with the electrons rotate their wave vectors upwards in the upper electrojet and downwards in the lower electrojet during the day, and vice versa at night. Furthermore, the altitude region of largest linear growth rate is also the one with the weakest refraction rate. Consequently, computations of the ray-path integrated wave growth shows that this region would dominate the backscatter spectrum from the electrojet if linear theory were valid, and it is further noted that the maximum amplitude wave should have phase velocities exceeding the ion acoustic speed. It was concluded that propagation alone, without inclusion of nonlinear effects, cannot explain backscatter observations of a constant Doppler frequency shift given by the ion acoustic speed.

Lee, K.↗

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

Heat Entrapment Effects Within Liquid Acquisition Devices

We introduce a model problem to address heat entrapment effects or the local accumulation of thermal energy within liquid acquisition devices. We show that the parametric space consists of six parameters, namely the Rayleigh and Prandtl numbers, the aspect ratio, and heat flux ratios for the bottom, side, and top boundaries of the enclosure. For the range of Ra considered 1 to 10(sup 9), beyond Ra on the order of 10(sup 5), convective instability is the dominant mode of convection in comparison to natural convection. The flow field transitions to asymmetric modes at Ra on the order of 10(sup 7). Direct numerical simulation of a large geometric length scale prototype for Ra on the order of 10(sup 9) shows that the flow field evolves from small wavelength instability which gives rise to nonlinear growth of thermals, propagation of the instability occurs via growth of secondary and tertiary modes, and a travelling wave mode occurs prior to asymmetry. The effect of a large aspect ratio is to increase the number of modes in the vertical direction. Due to the slow diffusion of heat in the prototype, asymptotic states are not readily attained, we show that dynamical similarity can be used for a model which allows the attainment of asymptotic states and that transition to a chaotic state occurs for Ra on the order of 10(sup 9) via a broadband power spectrum. These dynamical events show that for the baseline condition in which heat is absorbed from background laboratory environment, higher heat flux is absorbed at the top and bottom boundaries of the enclosure than a nominal value of 34.9 ergs per square centimeter -second.

Duval, W. M. B.↗

Nonlinear Internal Waves in the South China Sea During ASIAEX

Internal wave distribution maps have been compiled from more than one hundred ERS-1/2, RADARSAT and Space Shuttle SAR images in the South China Sea (SCS) from 1993 to 2000. Based on these distribution maps, most of internal waves in the northeast part of SCS were propagating westward. The wave crest can be as long as 200 km with amplitude of 100 m, due to strong current from the Kuroshio branching out into the SCS. In recent Asian Seas International Acoustics Experiment (ASIAEX), moorings have been deployed in April 2000 and May 2001. Simultaneous RADARSAT ScanSAR images have been collected during the field test to integrate with the model and the in-situ measurements in the SCS. During ASIAEX in May 2001, many large internal waves were observed at the test area and were the major features for acoustic volume interaction. The environmental parameters have been calculated based on extensive CTD castings and mooring data. Nonlinear internal wave models have been applied to simulate the wave evolution on the continental shelf and the results compare reasonably with mooring measurements. The evolution and dissipation of huge internal waves on the shelf break, mode-two waves, elevation waves, and wave-wave interaction are very important issues for acoustic propagation. The implication of internal wave effects on acoustic propagation will also be discussed.

ACOUSTIC PROPAGATION↗

Modified radiative transfer theory for a two-layer random medium

Modified radiative transfer (MRT) equations appropriate for electromagnetic wave propagation in bounded random media are derived from the Bethe-Salpeter equation with the ladder approximation and the Dyson equation with the nonlinear approximation. The MRT equations are then solved with the first-order renormalization approximation to obtain analytical results for the backscattering cross sections of a two-layer random medium with arbitrary three-dimensional correlation functions. The coherent effects of the MRT theory are illustrated and comparisons are made with backscattering cross sections obtained with the first Born approximation to the wave equation.

Zuniga, M. A.↗

Time dependent simulation of cosmic-ray shocks including Alfven transport

Time evolution of plane, cosmic-ray modified shocks was simulated numerically for the case with parallel magnetic fields. Computations were done in a 'three-fluid' dynamical model incorporating cosmic-ray and Alfven wave energy transport equations. Nonlinear feedback from the cosmic-rays and Alfven waves is included in the equation of motion for the underlying plasma, as is the finite propagation speed and energy dissipation of the Alfven waves. Exploratory results confirm earlier, steady state analyses that found these Alfven transport effects to be potentially important when the upstream Alfven speed and gas sound speeds are comparable. As noted earlier Alfven transport effects tend to reduce the transfer of energy through a shock from gas to energetic particles. These studies show as well that the time scale for modification of the shock is altered in nonlinear ways. It is clear, however, that the consequences of Alfven transport are strongly model dependent and that both advection of cosmic-rays by the waves and dissipation of wave energy in the plasma will be important to model correctly when quantitative results are needed. Comparison is made between simulations based on a constant diffusion coefficient and more realistic diffusion models allowing the diffusion coefficient to vary in response to changes in Alfven wave intensity. No really substantive differences were found between them.

Jones, T. W.↗

Assessment of shock capturing schemes for resonant flows in nonlinear instability analysis

The paper presents computational assessment of advanced numerical schemes for nonlinear acoustic problems related to combustion instabilities in liquid rocket engines. Several time-accurate, shock capturing schemes have been evaluated on a benchmark, closed-end resonant pipe flow problem. It involves the numerical solution of inviscid, compressible gas dynamics equations to predict acoustic wave propagation, wave steepening, formation of shocks, acoustic energy dissipation and wave-wall reflection for several hundred wave cycles. It was demonstrated that high accuracy TVD type schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time. The selected scheme was then applied to analyze the acoustic responses of resonant pipe-resonator, radial acoustic modes and hub-baffle configurations. Interesting observations of wave shape and damping characteristics have been drawn from presented computational studies.

Przekwas, A. J.↗

Strongly nonlinear evolution of low-frequency wave packets in a dispersive plasma

The evolution of strongly nonlinear, strongly modulated wave packets is investigated in a dispersive plasma using a hybrid numerical code. These wave packets have amplitudes exceeding the strength of the external magnetic field, along which they propagate. Alfven (left helicity) wave packets show strong steepening for p < 1, while fast (fight heIicity) wave packets hardly steepen for any beta. Substantial regions of opposite helicity form on the leading side of steepened Alfven wave packets. This behavior differs qualitatively from that exhibited by the solutions to the derivative nonlinear Schrodinger (DNLS) equation.

Vasquez, Bernard J.↗

Effects of ion two-stream instability on auroral ion heating

The ion two-stream instability develops between upflowing H(+) and O(+) in the auroral acceleration region, if the two species are accelerated to different velocities in a semistationary electrostatic potential. Parallel propagating modes are unstable only for small relative drifts, while oblique modes remain unstable at larger drifts and, consequently, higher altitudes. The effects of oblique modes on ion heating are investigated, and it is found that O(+) heats primarily by trapping in the direction of wave propagation during the early, coherent phase of wave growth in the temporal evolution problem. After undergoing trapping along the magnetic field direction, H(+) heats primarily by quasi-linear diffusion when the spectrum broadens nonlinearly. The parallel ion distributions evolve by formation of a high-energy (low-energy) tail on oxygen (hydrogen) along with bulk slowing of hydrogen, which results in more energetic oxygen than hydrogen. Total heating is greatest in the direction of wave propagation, and this effect is greater for hydrogen than for oxygen.

Roth, I.↗

Structure in the nonlinear saturation spectrum of parametric instabilities.

The nonlinear saturation spectrum of the decay instability is obtained in the limit of small spontaneous emission, for comparable ion and electron temperatures, from numerical solutions of a kinetic equation based on an accurate expression for the nonlinearity. The spectral energy occupies several pairs of isolated saucer-shaped regions in wave-vector space. The regions increase in thickness, angular diameter, and number as the pump power is increased. The theory thus predicts the generation of waves propagating in directions which can differ substantially from the direction of the pump field. Ionospheric observations confirm this prediction; they were difficult to reconcile with the predictions of previous theories based on an approximated expression for the nonlinearity. The present work also corrects the results of previous one-dimensional theories that used an accurate expression for the nonlinearity and predicted 'spectral lines' in the limit of vanishing spontaneous emission.

Fejer, J. A.↗

Convective amplification of type I irregularities in the equatorial electrojet.

Wave propagation and refraction of 'type I' irregularities in the equatorial electrojet are investigated. Quantitative calculation of wave refraction in a model electrojet shows that the direction of wave refraction must change sign at one altitude. Waves propagating with the electrons rotate their wave vectors upward in the upper electrojet and downward in the lower electrojet during the day, and vice versa at night. Furthermore, the altitude region of largest linear growth rate is also the one with the weakest refraction rate. Consequently, computations of the ray-path integrated wave growth show that this region would dominate the backscatter spectrum from the electrojet if linear theory were valid, and it is further noted that the maximum amplitude wave should have phase velocities exceeding the ion acoustic speed. It is therefore concluded that propagation alone, without inclusion of nonlinear effects, cannot explain backscatter observations of a constant Doppler frequency shift given by the ion acoustic speed.

Lee, K.↗

Numerical simulation of MHD shock waves in the solar wind

The effects of the interplanetary magnetic field on the propagation speed of shock waves through an ambient solar wind are examined by numerical solutions of the time-dependent nonlinear equations of motion. The magnetic field always increases the velocity of strong shocks. Although the field may temporarily slow down weak shocks inside 1 AU, it eventually also causes weak shocks to travel faster than they would without the magnetic field at larger distances. Consistent with the increase in the shock velocity, the gas pressure ratio across a shock is reduced considerably in the presence of the magnetic field. The numerical method is used to simulate (starting at 0.3 AU) the large deceleration of a shock observed in the lower corona by ground-based radio instrumentation and the more gradual deceleration of the shock in the solar wind observed by the Pioneer 9 and Pioneer 10 spacecraft.

Steinolfson, R. S.↗

Excitation of turbulence by density waves

A nonlinear system describes the microdynamical state of turbulence that is excited by density waves. It consists of an equation of propagation and a master equation. A group-scaling generates the scaled equations of many interacting groups of distribution functions. The two leading groups govern the transport processes of evolution and eddy diffusivity. The remaining sub-groups represent the relaxation for the approach of diffusivity to equilibrium. In strong turbulence, the sub-groups disperse themselves and the ensemble acts like a medium that offers an effective damping to close the hierarchy. The kinetic equation of turbulence is derived. It calculates the eddy viscosity and identifies the effective damping of the assumed medium self-consistently. It formulates the coupling mechanism for the intensification of the turbulent energy at the expense of the wave energy, and the transfer mechanism for the cascade. The spectra of velocity and density fluctuations find the power law k sup-2 and k sup-4, respectively.

Tichen, C. M.↗

Probabilistic finite elements for transient analysis in nonlinear continua

The probabilistic finite element method (PFEM), which is a combination of finite element methods and second-moment analysis, is formulated for linear and nonlinear continua with inhomogeneous random fields. Analogous to the discretization of the displacement field in finite element methods, the random field is also discretized. The formulation is simplified by transforming the correlated variables to a set of uncorrelated variables through an eigenvalue orthogonalization. Furthermore, it is shown that a reduced set of the uncorrelated variables is sufficient for the second-moment analysis. Based on the linear formulation of the PFEM, the method is then extended to transient analysis in nonlinear continua. The accuracy and efficiency of the method is demonstrated by application to a one-dimensional, elastic/plastic wave propagation problem. The moments calculated compare favorably with those obtained by Monte Carlo simulation. Also, the procedure is amenable to implementation in deterministic FEM based computer programs.

Liu, W. K.↗

Excitation of gravity waves in common envelopes

We study the excitation of gravity waves by a low-mass companion orbiting inside the envelope of a giant star, concentrating on brown dwarfs inside the envelope of asymptotic giant branch stars. Efficient g-wave excitations occur only after the brown dwarf has spiraled-in to the radiative zone, well inside the envelope, of the asymptotic giant branch star. The brown dwarf excites g-waves when its orbital radius is about 3-10 solar radii. At this stage of the evolution the envelope mass is below 0.1 solar mass. The g-waves propagate inward from the secondary orbit, carrying angular momentum and energy. We find that the angular momentum transport leads to an efficient spin-up of the inner envelopes. The differential rotation between the envelope and core and nonlinear wave effects, can cause a mixing of heavy elements from the core to the envelope.

Soker, Noam↗

Acoustoelastic characterization of plates using zero group velocity Lamb modes

Acoustoelasticity, a characteristic of material anharmonicity, gives rise to a link between wave propagation velocity and the stress state in materials. Ultrasonic techniques to monitor this coupling, particularly with high sensitivity and in a noncontact manner, can have widespread application both in the quantification of applied and residual stress and in the characterization of nonlinear material behavior through measurement of higher order elastic constants. Here, we use a laser ultrasonic technique to excite and detect zero group velocity (ZGV) Lamb wave resonances in aluminum plates under uniaxial loading. A laser line source is used to excite these resonances at different orientations with respect to the applied load, and the signals are detected using an interferometer. The effects of stress and source orientation on ZGV resonance frequencies are validated using the theory of acoustoelastic Lamb wave propagation. In addition, a model-based inversion technique is used to extract Murnaghan's third-order elastic constants from measurements of the stress dependence of the first two ZGV modes generated parallel and perpendicular to the applied load. Laser generation and detection of ZGV resonances is shown to be an effective and powerful approach for the noncontact and nondestructive acoustoelastic characterization of elastic waveguides.

36 MATERIALS SCIENCE↗

Periodic Narrowband Radio Wave Emissions and Inward Plasma Transport at Saturn's Magnetosphere

The abrupt brightening of an Energetic Neutral Atom (ENA) blob or cloud has been interpreted as plasma injection in the Kronian magnetosphere (termed ENA injection herein). Morphologically, there appears to be two types of abrupt ENA cloud brightening: (1) the brightening of a large cloud usually seen at distances >10–12 R {sub s} (R {sub s} ∼ 60,268 km) in the midnight or postmidnight region; (2) the brightening of a smaller cloud usually seen at distances <10–12 R {sub s} around 21–03 magnetic local time. Among many radio waves observed at Saturn, type 2 ENA injections correlate best with the 5 kHz narrowband (NB) waves. Using Cassini Ion and Neutral Camera (INCA) and Radio and Plasma Wave Science (RPWS) data, we examine the periodicities of the type 2 ENA injections and the 5 kHz NB emissions as well as their cross-correlations, which have been previously used to measure the lag times or phase differences. Because correlational analysis can only establish linear relationships, we also use mutual information to establish linear and nonlinear relationships. On average, the peak of the 5 kHz NB emission lags those of the type 2 ENA injection by about a few minutes to 2 hr. The injection of hot plasma to the inner magnetosphere can lead to temperature anisotropy, which can generate electrostatic upper hybrid waves, which upon encountering the high-density gradient at the outer edge of the Enceladus density torus, can mode convert to Z mode and then to O mode. The 5 kHz NB waves commonly propagate in the O mode.

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamics in Layer Models of Solid Flame Propagation

Self-propagating high-temperature synthesis (SHS) is a process in which combustion waves, e.g., "solid flames", which are considered here, are employed to synthesize desired materials. Like many other systems, SHS is a pattern forming system. The problem of describing experimentally observed patterns and of predicting new, as yet unobserved, patterns continues to attract the attention of scientists and mathematicians due to the fundamental significance of the phenomena in combustion in particular, and in nonlinear science in general. Here, we analyze the dynamics of solid flame propagation in a 2D region by considering the region to be composed of parallel, identical layers aligned along the direction of propagation and having thermal contact. Each layer is then described by wave propagation in 1D, with the transverse Laplacian replaced by a term describing heat exchange between neighboring layers. This configuration is the simplest model of a 2D system because it accounts, in a simple way, for the principal feature of the problem, i.e., heat exchange between neighbors in the transverse direction. For simplicity, we describe the situation for two layers. Because the layers are identical, uniformly propagating waves in each layer must be identical, independent of the heat exchange rate alpha. When the Zeldovich number Z exceeds a critical value Z(sub c), which depends on alpha, uniformly propagating waves become unstable. The stability diagram for the two coupled layers reproduces that for the full 2D problem after appropriate identification of parameters in the two problems. Depending on parameter values, we determine three different steady-state dynamical behaviors (though additional behaviors are also expected to occur). The three behaviors are: (i) waves in each layer which pulsate in phase as they propagate, so that together they form a single pulsating propagating wave; (ii) the waves in each layer are no longer identical, and antiphase pulsations occur, with the two waves alternately advancing and receding as they propagate. This wave is the analog of the spinning wave on the surface of a circular cylinder; (iii) finally, there is a region of bistability between the in phase and antiphase waves. with each having its own domain of attraction, so that which of the two behaviors occur depends on the condition of initiation of the wave. The results of our computations indicate a qualitative similarity in the behavior of combustion waves in the layer model and in the full 2D model. Specifically, due to the similarity between the small alpha wave behavior in the layer model and the large diameter behavior in the model of waves on the surface of a cylinder, we are able to predict the behavior of the mean velocity for the waves on the cylinder, where computations of the full problem can be rather difficult. We also compute solutions for three or more layers. The results of our computations prompt us to predict that, while planar uniformly propagating waves are unstable, the wave will be quasiplanar, i.e., the resulting spinning waves have very low amplitude hot spots, and travel with the velocity close to that of the uniformly propagating wave. Such waves may be difficult to distinguish from uniformly propagating waves in experiments. We also find that for both the layer model and full 2D problem, steady-state time-dependent waves, e.g., pulsating and spinning wave, have a conserved quantity H which characterizes the excess energy in the wave, just as in the case of uniformly propagating waves. The quantity H, which is generated by dissipation, does not vary in time and is proportional to the diffusivity and caloricity of the system, and inversely proportional to the mean wave velocity.

Aldushin, A. P.↗