Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear oscillation”

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 613 records · Page 34

The attractor dimension of solar decimetric radio pulsations

The temporal characteristics of decimetric pulsations and related radio emissions during solar flares are analyzed using statistical methods recently developed for nonlinear dynamic systems. The results of the analysis is consistent with earlier reports on low-dimensional attractors of such events and yield a quantitative description of their temporal characteristics and hidden order. The estimated dimensions of typical decimetric pulsations are generally in the range of 3.0 + or - 0.5. Quasi-periodic oscillations and sudden reductions may have dimensions as low as 2. Pulsations of decimetric type IV continua have typically a dimension of about 4.

Kurths, J.↗

Nonlinear attitude control of elastic spacecraft-antenna system

An approach to large angle rotational maneuvers of a spacecraft-beam-tip body (an antenna or a reflector) configuration based on nonlinear invertibility and linear feedback stabilization is presented. A control law u sub d is derived to obtain independent decoupled control of attitude angles, lateral elastic deflections, slopes due to bending and angular deflection due to torsion at the tip of the beam using torquers and force actuators. For the stabilization of the elastic oscillations, a linear feedback control law, u sub s, is obtained based on a linearized model about the terminal state augmented with a servo-compensator. Simulation results obtained for single axis control, for simplicity, show that large slewing and elastic mode stabilization can be accomplished in spite of uncertainty in the system using the total control u = u sub d + u sub s.

Singh, S. N.↗

On some numerical schemes for transonic flow problems

New second-order-accurate finite-difference approximations for a class of nonlinear PDEs of mixed type, which includes the two-dimensional low-frequency transonic small-disturbance equation (TSD) and full-potential equation (FP), are presented. For the TSD equation, the scheme is implemented via a time-splitting algorithm; the inclusion of flux limiters keeps the total variation nonincreasing and eliminates spurious oscillations near shocks. Global-linear-stability, TVD, and entropy-stability results are proved. Numerical results for the flow over a thin airfoil are presented. Current techniques used to solve the TSD equation may easily be extended to second-order accuracy by this method. For the FP equation, the new scheme requires no subsonic/supersonic switching and no numerical flux biasing. Global linear stability for all values of the Mach number is proved.

Mostrel, Marco Mosche↗

Higher-Order Spectral Analysis of F-18 Flight Flutter Data

Royal Australian Air Force (RAAF) F/A-18 flight flutter test data is presented and analyzed using various techniques. The data includes high-quality measurements of forced responses and limit cycle oscillation (LCO) phenomena. Standard correlation and power spectral density (PSD) techniques are applied to the data and presented. Novel applications of experimentally-identified impulse responses and higher-order spectral techniques are also applied to the data and presented. The goal of this research is to develop methods that can identify the onset of nonlinear aeroelastic phenomena, such as LCO, during flutter testing.

Silva, Walter A.↗

Quantum-well diode frequency multipliers - Varistor case

Local oscillators for heterodyne receivers at submillimeter wavelengths are typically made using a fundamental source followed by a harmonic frequency multiplier. An investigation of the required circuit embedding conditions for a possible new harmonic generator, the quantum-well resonant-tunneling diode, is summarized. A low-frequency multiplier has been tested that employs the resistive nonlinearity of the device as opposed to the reactive nonlinearity. The results show good agreement between practice and theory.

Batelaan, Paul D.↗

Nonlinear behavior of solar gravity modes driven by He-3 in the core. I - Bifurcation analysis

The nonlinear development of solar gravity modes driven by He-3 burning in the solar core is investigated by means of an idealized dynamical model. Possible outcomes that have been suggested in the literature include the triggering of subcritical direct convection, leading to core mixing, and the saturation of the excitation processes, leading to sustained finite-amplitude oscillations. The present simple model suggests that the latter is the more likely. The limiting amplitude of the oscillations is estimated, ignoring possible resonances with other gravity modes, to be of order 10 km/s at the solar surface. Such oscillations would be easily observable. That large-amplitude gravity modes have not been observed suggests either that these modes are not unstable in the present era or that they are limited to much smaller amplitudes by resonant coupling.

Merryfield, William J.↗

In-Flight Aeroelastic Stability of the Thermal Protection System on the NASA HIAD, Part II: Nonlinear Theory and Extended Aerodynamics

Conical shell theory and a supersonic potential flow aerodynamic theory are used to study the nonlinear pressure buckling and aeroelastic limit cycle behavior of the thermal protection system for NASA's Hypersonic Inflatable Aerodynamic Decelerator. The structural model of the thermal protection system consists of an orthotropic conical shell of the Donnell type, resting on several circumferential elastic supports. Classical Piston Theory is used initially for the aerodynamic pressure, but was found to be insufficient at low supersonic Mach numbers. Transform methods are applied to the convected wave equation for potential flow, and a time-dependent aerodynamic pressure correction factor is obtained. The Lagrangian of the shell system is formulated in terms of the generalized coordinates for all displacements and the Rayleigh-Ritz method is used to derive the governing differential-algebraic equations of motion. Aeroelastic limit cycle oscillations and buckling deformations are calculated in the time domain using a Runge-Kutta method in MATLAB. Three conical shell geometries were considered in the present analysis: a 3-meter diameter 70 deg. cone, a 3.7-meter 70 deg. cone, and a 6-meter diameter 70 deg. cone. The 6-meter configuration was loaded statically and the results were compared with an experimental load test of a 6-meter HIAD. Though agreement between theoretical and experimental strains was poor, the circumferential wrinkling phenomena observed during the experiments was captured by the theory and axial deformations were qualitatively similar in shape. With Piston Theory aerodynamics, the nonlinear flutter dynamic pressures of the 3-meter configuration were in agreement with the values calculated using linear theory, and the limit cycle amplitudes were generally on the order of the shell thickness. The effect of axial tension was studied for this configuration, and increasing tension was found to decrease the limit cycle amplitudes when the circumferential elastic supports were neglected, but resulted in more complex behavior when the supports were included. The nominal flutter dynamic pressure of the 3.7-meter configuration was significantly lower than that of the 3-meter, and it was found that two sets of natural modes coalesce to flutter modes near the same dynamic pressure. This resulted in a significant drop in the limit cycle frequencies at higher dynamic pressures, where the flutter mode with the lower frequency becomes more critical. Pre-buckling pressure loads and the aerodynamic pressure correction factor were studied for all geometries, and these effects resulted in significantly lower flutter boundaries compared with Piston Theory alone. The maximum dynamic pressure predicted by aerodynamic simulations of a proposed 3.7-meter HIAD vehicle was still lower than any of the calculated flutter dynamic pressures, suggesting that aeroelastic effects for this vehicle are of little concern.

Goldman, Benjamin D.↗

Computed and Experimental Flutter/LCO Onset for the Boeing Truss-Braced Wing Wind-Tunnel Model

This paper presents high fidelity Navier-Stokes simulations of the Boeing Subsonic Ultra Green Aircraft Research truss-braced wing wind-tunnel model and compares the results to linear MSC. Nastran flutter analysis and preliminary data from a recent wind-tunnel test of that model at the NASA Langley Research Center Transonic Dynamics Tunnel. The simulated conditions under consideration are zero angle of attack, so that structural nonlinearity can be neglected. It is found that, for Mach number greater than 0.78, the linear flutter analysis predicts flutter onset dynamic pressure below the wind-tunnel test and that predicted by the Navier-Stokes analysis. Furthermore, the wind-tunnel test revealed that the majority of the high structural dynamics cases were wing limit cycle oscillation (LCO) rather than flutter. Most Navier-Stokes simulated cases were also LCO rather than hard flutter. There is dip in the wind-tunnel test flutter/LCO onset in the Mach 0.76-0.80 range. Conditions tested above that Mach number exhibited no aeroelastic instability at the dynamic pressures reached in the tunnel. The linear flutter analyses do not show a flutter/LCO dip. The Navier-Stokes simulations also do not reveal a dip; however, the flutter/LCO onset is at a significantly higher dynamic pressure at Mach 0.90 than at lower Mach numbers. The Navier-Stokes simulations indicate a mild LCO onset at Mach 0.82, then a more rapidly growing instability at Mach 0.86 and 0.90. Finally, the modeling issues and their solution related to the use of a beam and pod finite element model to generate the Navier-Stokes structure mode shapes are discussed.

Bartels, Robert E.↗

Nonlinear optical processes for the mid-IR region

The 2 - 5 micron (mid-IR) region is of interest for a number of applications. Efficient up-conversion and down-conversion techniques are being developed to obtain optical sources at mid-IR wavelengths. These techniques are reviewed and recent results using AgGaSe2 are reported. Gain in AgGaSe2 as high as 13 has been observed with a pump wavelength of 1.73 microns and a signal wavelength of 3.39 microns. Optical parametric oscillation in the 2.9 - 6.8 micron region has also been demonstrated.

Schepler, K. L.↗

A theory of rotating stall of multistage axial compressors

A theoretical analysis was made of rotating stall in axial compressors of many stages, finding conditions for a permanent, straight-through traveling disturbance, with the steady compressor characteristic assumed known, and with simple lag processes ascribed to the flows in the inlet, blade passages, and exit regions. For weak disturbances, predicted stall propagation speeds agree well with experimental results. For a locally-parabolic compressor characteristic, an exact nonlinear solution is found and discussed. For deep stall, the stall-zone boundary is most abrupt at the trailing edge, as expected. When a complete characteristic having unstalling and reverse-flow features is adopted, limit cycles governed by a Lienard's equation are found. Analysis of these cycles yields predictions of recovery from rotating stall; a relaxation oscillation is found at some limiting flow coefficient, above which no solution exists. Recovery is apparently independent of lag processes in the blade passages, but instead depends on the lags originating in the inlet and exit flows, and also on the shape of the given characteristic diagram. Small external lags and tall diagrams favor early recovery. Implications for future research are discussed.

Moore, F. K.↗

Air-sea interactions and precipitation over the tropical oceans

In this lecture, the author principally discusses air-sea exchanges that are relevant to climate and global problems. The processes of interest are those acting over time scales of months to decades, which in some instances are influenced by smaller-time-scale processes, down to the diurnal time scale. The repsective influence of these processes varies with regions, seasons and scales over which they occur and, because these processes are mostly nonlinear, scale interactions can be quite complex. Owing to the breadth of the topic addressed, the discussion is mostly focused on the tropical regions where air-sea interactions and precipitation processes eventually affect the entire globe. This allows a look in more detail at some air-sea processes, such as those associated with the El Nino southern oscillation (ENSO). This oscillation, which affects the climate of the entire globe, acts over periods of a year or longer and is caused, primarily, by sea surface temperature (SST) variations in the tropical Pacific. As a result, SST variability is often used as an indicator of coupled ocean-atmosphere low-frequency variability. Global or basin scale processes can uniquely be observed from space-born instruments with the coverage required. Space based techniques have been developed during the last decade which can now be used to illustrate the scientific issues presented and the presentation concludes with an overview of some Earth Observing System (EOS) capabilities for addressing air-sea interactions and hydrology issues.

Gautier, C.↗

Low-Dimensional Dynamical Models of Thermal Convection

A low-dimensional dynamic model for transitional buoyancy-driven flow in a differentially heated tall enclosure is presented. The full governing partial differential equations with the associated boundary conditions are solved by a spectral element method for a cavity of aspect ratio A=20. Proper orthogonal decomposition is applied to the oscillatory solution at Prandtl number Pr=P tau (omega) = 0.71 and Grashof number G tau (omega) = 3.2 x 10 (exp 4) to construct empirical eigenfunctions. Using the four most energetic empirical eigenfunctions for the velocity and temperature as basis functions and applying Galerkin's method, a reduced model consisting of eight nonlinear ordinary differential equations is obtained. Close to the 'design' conditions (P tau(omega) G tau(omega)), the low-order model (LOM) predictions are in excellent agreement with the predictions of the full model. In particular, the critical Grashof number at the onset of the first temporal flow instability (Hopf bifurcation) was well as the frequency and amplitude of oscillations at supercritical conditions are in excellent agreement with the predictions of the full model. Far from the 'design' conditions, the LOM predicts the existence of multiple stable steady solutions at large values of G tau, and a unique stable steady solution at small values of G tau, and exhibits hysteretic behavior that is qualitatively similar to that observed in direct numerical simulations based on the full model.

Liakopoulos, Anthony↗

Computation of Feedback Aeroacoustic System by the CE/SE Method

It is well known that due to vortex shedding in high speed flow over cutouts, cavities, and gaps, intense noise may be generated. Strong tonal oscillations occur in a feedback cycle in which the vortices shed from the upstream edge of the cavity convect downstream and impinge on the cavity lip, generating acoustic waves that propagate upstream to excite new vortices. Numerical simulation of such a complicated process requires a scheme that can: (1) resolve acoustic waves with low dispersion and numerical dissipation, (2) handle nonlinear and discontinuous waves (e.g. shocks), and (3) have an effective (near field) nonreflecting boundary condition (NRBC). The new space time conservation element and solution element method, or CE/SE for short, is a numerical method that meets the above requirements.

Loh, Ching Y.↗

On the nonlinear three dimensional instability of Stokes layers and other shear layers to pairs of oblique waves

The nonlinear evolution of a pair of initially oblique waves in a high Reynolds Number Stokes layer is studied. Attention is focused on times when disturbances of amplitude epsilon have O(epsilon(exp 1/3)R) growth rates, where R is the Reynolds number. The development of a pair of oblique waves is then controlled by nonlinear critical-layer effects. Viscous effects are included by studying the distinguished scaling epsilon = O(R(exp -1)). This leads to a complicated modification of the kernel function in the integro-differential amplitude equation. When viscosity is not too large, solutions to the amplitude equation develop a finite-time singularity, indicating that an explosive growth can be introduced by nonlinear effects; we suggest that such explosive growth can lead to the bursts observed in experiments. Increasing the importance of viscosity generally delays the occurrence of the finite-time singularity, and sufficiently large viscosity may lead to the disturbance decaying exponentially. For the special case when the streamwise and spanwise wavenumbers are equal, the solution can evolve into a periodic oscillation. A link between the unsteady critical-layer approach to high-Reynolds-number flow instability, and the wave vortex approach is identified.

Wu, Xuesong↗

Lee waves: Benign and malignant

The flow of an incompressible fluid over an obstacle will produce an oscillation in which buoyancy is the restoring force, called a gravity wave. For disturbances of this scale, the atmosphere may be treated as dynamically incompressible, even though there exists a mean static upward density gradient. Even in the linear approximation - i.e., for small disturbances - this model explains a great many of the flow phenomena observed in the lee of mountains. However, nonlinearities do arise importantly, in three ways: (1) through amplification due to the decrease of mean density with height; (2) through the large (scaled) size of the obstacle, such as a mountain range; and (3) from dynamically singular levels in the fluid field. These effects produce a complicated array of phenomena - large departure of the streamlines from their equilibrium levels, high winds, generation of small scales, turbulence, etc. - that present hazards to aircraft and to lee surface areas. The nonlinear disturbances also interact with the larger-scale flow in such a manner as to impact global weather forecasts and the climatological momentum balance. If there is no dynamic barrier, these waves can penetrate vertically into the middle atmosphere (30-100 km), where recent observations show them to be of a length scale that must involve the coriolis force in any modeling. At these altitudes, the amplitude of the waves is very large, and the phenomena associated with these wave dynamics are being studied with a view to their potential impact on high performance aircraft, including the projected National Aerospace Plane (NASP). The presentation shows the results of analysis and of state-of-the-art numerical simulations, validated where possible by observational data, and illustrated with photographs from nature.

Wurtele, M. G.↗

Oscillation Amplitude Growth for a Decelerating Object with Constant Pitch Damping

The equations governing the deceleration and oscillation of a blunt body moving along a planar trajectory are re-expressed in the form of the Euler-Cauchy equation. An analytic solution of this equation describes the oscillation amplitude growth and frequency dilation with time for a statically stable decelerating body with constant pitch damping. The oscillation histories for several constant pitch damping values, predicted by the solution of the Euler-Cauchy equation are compared to POST six degree-of-freedom (6-DoF) trajectory simulations. The simulations use simplified aerodynamic coefficients matching the Euler-Cauchy approximations. Agreement between the model predictions and simulation results are excellent. Euler-Cauchy curves are also fit through nonlinear 6-DoF simulations and ballistic range data to identify static stability and pitch damping coefficients. The model os shown to closely fit through the data points and capture the behavior of the blunt body observed in simulation and experiment. The extracted coefficients are in reasonable agreement with higher fidelity, nonlinear parameter identification results. Finally, a nondimensional version of the Euler-Cauchy equation is presented and shown to be a simple and effective tool for designing dynamically scaled experiments for decelerating blunt capsule flight.

Schoenenberger, Mark↗

A synthesis theory for the externally excited adaptive system /EEAS/

The externally excited adaptive system (EEAS) is a two-degree-of-freedom feedback system with a nonlinearity which is saturated hard by an external periodic signal. Under certain conditions, the EEAS responds quasi-linearly to command and plant disturbance signals, permitting the development of a quantitative synthesis theory for satisfying system tolerances despite large plant uncertainty. The great advantage of the EEAS is its zero sensitivity to plant gain variations, a property it shares with the self-oscillating adaptive system (SOAS). The EEAS is, however, more flexible than the SOAS in satisfying the quasi-linearity constraints. The essential difference is that in the EEAS the loop transmission bandwidth is not rigorously tied to the 'carrier' signal, as it is in the SOAS. There is a class of problems for which the EEAS is superior to the purely linear system, which in turn is superior to the SOAS. The superiority of the EEAS over the SOAS is especially marked in the case of significant plant disturbances, which generally vitiate a SOAS design.

Horowitz, I. M.↗

Viscous boundary layers in rotating fluids driven by periodic flows

The paper analyzes the boundary layers formed in a rotating fluid by an oscillating flow over an infinite half plate, with particular attention paid to the effects of unsteadiness, the critical latitude effect and the structure of the solution to the boundary layer equations at resonance. The Navier-Stokes boundary layer equations are obtained through an asymptotic expansion with the incorporation of the Rossby and Ekman numbers and are analyzed as the sum of a nonlinear steady solution and a linearized unsteady solution. The solution is predominantly composed of two inertial wave vector components, one circularly polarized to the left and the other circularly polarized to the right. The problem considered here has relevance in oceanography and meteorology, with special reference to the unsteady atmospheric boundary layer.

Bergstrom, R. W.↗