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 649 records · Page 36

Singular spectrum analysis in nonlinear dynamics, with applications to paleoclimatic time series

Two dimensions of a dynamical system given by experimental time series are distinguished. Statistical dimension gives a theoretical upper bound for the minimal number of degrees of freedom required to describe the attractor up to the accuracy of the data, taking into account sampling and noise problems. The dynamical dimension is the intrinsic dimension of the attractor and does not depend on the quality of the data. Singular Spectrum Analysis (SSA) provides estimates of the statistical dimension. SSA also describes the main physical phenomena reflected by the data. It gives adaptive spectral filters associated with the dominant oscillations of the system and clarifies the noise characteristics of the data. SSA is applied to four paleoclimatic records. The principal climatic oscillations and the regime changes in their amplitude are detected. About 10 degrees of freedom are statistically significant in the data. Large noise and insufficient sample length do not allow reliable estimates of the dynamical dimension.

Vautard, R.↗

Investigation of Reduced-Order Modeling for Aircraft Stability and Control Prediction

High-fidelity computational fluid dynamics tools offer the potential to approximate increments for ground-to-flight scaling effects, as well as to augment the dynamic damping derivative data for motion-based flight simulators. Unfortunately, the computational expense is currently prohibitive for populating a complete simulator database. This work investigates an existing surrogate-based, indicial response reduced-order model methodology as a means to efficiently augment a flight simulator database with high-fidelity nonlinear aerodynamic damping derivatives. Creation of the reduced-order model is based on the superposition integrals of the step response with the derivative of its corresponding input signal. Step responses are calculated using a computational grid motion approach that separates the effects of angle of attack and sideslip angle from angular rates, and rates from angle of attack and sideslip. It is demonstrated that the transients produced during the start of a forced-oscillation motion are captured by the reduced-order model to the level of fidelity of a comparable computational solution. Aerodynamic coefficients computed within minutes by the reduced-order model for an aircraft undergoing an 18-second half Lazy-8 maneuver and a 25-second Immelmann turn maneuver are compared with those from full computational flight solutions that required days to complete. Finally, a cost-benefit assessment is included that demonstrates a compelling advantage for this approach. d for maneuvering, flexible vehicles.

Frink, Neal T.↗

The prediction of nonlinear longitudinal combustion instability in liquid propellant rockets

An analytical technique was developed to solve nonlinear longitudinal combustion instability problems. The analysis yields the transient and limit cycle behavior of unstable motors and the threshold amplitude required to trigger a linearly stable motor into unstable operation. The limit cycle waveforms were found to exhibit shock wave characteristics for most unstable engine operating conditions. A method of correlating the analytical solutions with experimental data was developed. Calculated results indicate that a second-order solution adequately describes the behavior of combustion instability oscillations over a broad range of engine operating conditions, but that higher order effects must be accounted for in order to investigate engine triggering.

Lores, M. E.↗

Nonlinear twist-kink instability of a coronal loop

Three-dimensional magnetoinductive particle simulations are used to demonstrate that the mechanical twisting motion applied to a magnetized plasma column induces a current aligned to the external magnetic field direction, pinches the plasma and magnetic fields, and stores the energy in poloidal magnetic fields. As the twist motion continues, the field lines locally begin to wrap around the plasma more than one revolution. A strong MHD instability sets in that is a mixture of kink and ballooning modes, releasing the magnetic energy and causing destruction of coherent column structure and flows of turbulent plasma. A similar episode ensues, exhibiting relaxation oscillations. The buildup of poloidal fields and structure and its sudden release driven by the twist motion may be a model for the solar coronal loop dynamics which exhibits a slow energy buildup with some photospheric motion and a sudden energy release by flares.

Zaidman, E. G.↗

Convective Effects During the Physical Vapor Transport Process. II - Thermosolutal Convection

The effect of an inert gas on the diffusive-convective physical vapor transport process is investigated for the case when the temperature gradient is stabilizing, and the concentration gradient destabilizing, for a wide parametric range. When an inert gas is present, the thermal and solutal convection oppose each other. The solutal field is destabilizing while the thermal field and the advective-diffusive flux stabilize the flow field. When the pressure of the inert component is increased, the stabilizing effect of the advective-diffusive flux is decreased; thus, convection becomes more vigorous. The nonlinear dynamics of the flow field here show a transition from quasi-periodic to chaotic state. When both stabilizing mechanisms are present, the flow field shows a transition to a steady state. Toward steady state, growth and amalgamation of rolls occur, which result in an overturning motion. This leads to a superposed flow consisting of one roll and a unidirectional flow. However, when the pressure is increased, the advective-diffusive stability mechanism is decreased. and oscillations of the flow field occur. The low gravity environment is effective at eliminating oscillatory behavior of the flow field and results in uniform temperature and concentration gradients.

Duval, Walter M. B.↗

A Numerical Method of Calculating Propeller Noise Including Acoustic Nonlinear Effects

Using the transonic flow fields(s) generated by the NASPROP-E computer code for an eight blade SR3-series propeller, a theoretical method is investigated to calculate the total noise values and frequency content in the acoustic near and far field without using the Ffowcs Williams - Hawkings equation. The flow field is numerically generated using an implicit three dimensional Euler equation solver in weak conservation law form. Numerical damping is required by the differencing method for stability in three dimensions, and the influence of the damping on the calculated acoustic values is investigated. The acoustic near field is solved by integrating with respect to time the pressure oscillations induced at a stationary observer location. The acoustic far field is calculated from the near field primitive variables as generated by NASPROP-E computer code using a method involving a perturbation velocity potential as suggested by Hawkings in the calculation of the acoustic pressure time-history at a specified far field observed location. the methodologies described are valid for calculating total noise levels and are applicable to any propeller geometry for which a flow field solution is available.

Korkan, K. D.↗

Resonant frequency calculations using a hybrid perturbation-Galerkin technique

A two-step hybrid perturbation Galerkin technique is applied to the problem of determining the resonant frequencies of one or several degree of freedom nonlinear systems involving a parameter. In one step, the Lindstedt-Poincare method is used to determine perturbation solutions which are formally valid about one or more special values of the parameter (e.g., for large or small values of the parameter). In step two, a subset of the perturbation coordinate functions determined in step one is used in Galerkin type approximation. The technique is illustrated for several one degree of freedom systems, including the Duffing and van der Pol oscillators, as well as for the compound pendulum. For all of the examples considered, it is shown that the frequencies obtained by the hybrid technique using only a few terms from the perturbation solutions are significantly more accurate than the perturbation results on which they are based, and they compare very well with frequencies obtained by purely numerical methods.

Geer, James F.↗

Resonant frequency calculations using a hybrid perturbation-Galerkin technique

A two-step hybrid perturbation Galerkin technique is applied to the problem of determining the resonant frequencies of one or several degrees of freedom nonlinear systems involving a parameter. In one step, the Lindstedt-Poincare method is used to determine perturbation solutions which are formally valid about one or more special values of the parameter (e.g., for large or small values of the parameter). In step two, a subset of the perturbation coordinate functions determined in step one is used in Galerkin type approximation. The technique is illustrated for several one degree of freedom systems, including the Duffing and van der Pol oscillators, as well as for the compound pendulum. For all of the examples considered, it is shown that the frequencies obtained by the hybrid technique using only a few terms from the perturbation solutions are significantly more accurate than the perturbation results on which they are based, and they compare very well with frequencies obtained by purely numerical methods.

Geer, James F.↗

Unsteady transonic flow simulation on a full-span-wing-body configuration

The presence of a body influences both the aerodynamic and aeroelastic performance of wings. Such effects are more pronounced in the transonic regime. To accurately account for the effect of the body, particularly when the wings are experiencing asymmetric modal motions, it is necessary to model the full configuration in the nonlinear transonic regime. In this study, full-span-wing-body configurations are simulated for the first time by a theoretical method that uses the unsteady potential equations based on the small-disturbance theory. The body geometry is modeled exactly as the physical shape, instead of as a rectangular box, which has been done in the past. Steady pressure computations for wing-body configurations compare well with the available experimental data. Unsteady pressure computations when the wings are oscillating in asymmetric modes show significant influence of the body.

Guruswamy, Guru P.↗

Diagnostics and control of wavenumber stability and purity of tunable diode lasers relevant to their use as local oscillators in heterodyne systems

Initial operation of the tunable diode lasers (TDL) showed that it was not possible to adjust the wavenumber to one selected a priori in the TDL tuning range. During operation, the operating point would change by 0.1/cm over the longer term with even larger changes occurring during some thermal cycles. Most changes during thermal cycling required using lower temperatures and higher currents to reach the former wavenumber (when it could be reached). In many cases, an operating point could be selected by changing TDL current and temperature to give both the desired wavenumber and most of the power in a single mode. The selection procedure had to be used after each thermal cycling. Wavenumber nonlinearities of about 10% over a 0.5 cm tuning range were observed. Diagnostics of the single mode selected by a grating monochromator showed wavenumber fine structure under certain operating conditions. The characteristics due to the TDL environment included short term wavenumber stability, the instrument lineshape function, and intermediate term wavenumber stability.

Poultney, S.↗

Nonlinear response of infinitely long circular cylindrical shells to subharmonic radial loads

The nonlinear response of infinitely long circular cylindrical shells (thin circular rings) in the presence of a two-to-one internal (autoparametric) resonance to a subharmonic excitation of order one-half of the higher mode is analyzed with the multiple-scale method. Four autonomous first-order ordinary differential equations are derived for the modulation of the amplitudes and phases of the interacting models. These modulation equations are used to determine the fixed points and their stability. The fixed points correspond to periodic oscillations of the shell, whereas the limit-cycle solutions of the modulation equations correspond to amplitude and phase-modulated oscillations of the shell. The force response curves exhibit saturation, jumps, and Hopf bifurcation. As excitation frequency changes, all limit cycles deform and lose stability through either pitchfork or cyclic-fold (saddle-node) bifurcations. Some of these saddle-node bifurcations cause a transition to chaos. The pitchfork bifurcations break the symmetry of the limit cycles.

Nayfeh, Ali H.↗

Hydrodynamic Stability and Magnetic Reconnection in Disks and Stars

The purpose of this grant is to study parametric instability. The simplest example of parametric instability is a harmonic oscillator with a periodic modulation of the spring constant. If the modulation frequency is close to twice the natural frequency of the oscillator, the amplitude of oscillation tends to grow exponentially. The growth rate is proportional to the strength of the modulation, but it also depends upon the closeness to resonance of the two frequencies, and upon natural damping rate or "Q" of the oscillator. Parametric instabilities are very common in physics. A familiar example is a jogger's ponytail--normally a very strongly damped pendulum, it can be destabilized by the variation in effective gravity during the jogger's stride. Observation confirms that the period of the pendulum is half that of the jogger's vertical motion. In astrophysics, parametric instability may occur by external tidal forcing, or by interaction among eigenmodes. In the latter case, an energetic eigenmode may destabilize modes of half its frequency, provided some weak nonlinearity exists to couple them. Under a previous Astrophysical Theory grant (NAGW-2419), the PI discovered a parametric instability of tidally forced disks such as the accretion disks in cataclysmic variables and X ray binaries [2]. The destabilized modes are tightly-wound, incompressible, three-dimensional waves analogous to g-modes and r-modes in stars. Later work has confirmed our analysis [4]. It was hoped that these modes might provide a source of turbulence and angular momentum transport in accretion disks. However, a follow-up investigation of this instability by local numerical simulations, although confirming the analytically estimated growth rates, found negligible angular momentum flux [3]. Other work, partly supported by the ATP, now strongly indicates that the transport mechanism in such disks is magnetohydrodynamic turbulence [6]. Nevertheless, the parametric mechanism may truncate the outer edges of disks in close binaries [2], and it may be important in disks of very low ionization such as protostellar disks, or even cataclysmic-variable disks in quiescence where the MHD mechanism may be ineffective [5]. All analyses up to 1996 were done in a local approximation where the orbital frequency, shear rate, and tidal field were treated as constants. The locally computed growth rate turns out to depend strongly on radius, and it was unclear how to average these local rates to obtain the correct global rate. This is a critical issue for accretion disks in close binaries, because the local growth rate is comparable to the orbital frequency towards the outer edge of the disk but decreases rapidly inwards. Paper #1 examined this issue in a simplified global model where the destabilizing terms vary with position. We found that the global growth rate is essentially equal to the maximum local rate, provided that the latter is smoothed over a radial range equal to the distance that the destabilized wave propagates at its group speed in one growth time. Thus, in an accretion disk, waves would grow rapidly in the outer parts but would propagate both inwards and outwards at a maximum group speed of order the disk thickness divided by the orbital period.

Goodman, Jeremy↗

The Investigation of the Effects of Gravity on Single Bubble Sonoluminescence

In single bubble following it's rapid collapse each cycle of oscillation of an ultrasonic field. Since widely varying length and time scales affect the bubble dynamics and optical emission processes, it is difficult to anticipate the importance of the effects of gravity present for observations on earth. Our bubble is driven in an acoustically resonating cavity at it's first harmonic mode. The acoustical radiation pressure (Bjerknes force) will then keep it suspended in the center near the pressure antinode. When driven in a region where the diffusive processes balance the bubble it acts in a nonlinear but regular way, emitting a short (approx. 200ps) burst of light each acoustic cycle. Balancing the Bjerknes force with buoyancy, as in, we can see that the bubble should be displaced from the velocity node approximately 20m at normal gravity. Therefore, water flows past the bubble at the time of collapse. Gravitation also changes the ambient pressure at the bubble's location, as Delta.P = rho.g.h this gives a change of approximately -0.5% in our experiment when going from 1.8g to 0g. Studies of ambient pressure changes were also done in order to assess these effects. Inside a pressure sealed chamber a spherical glass cell is filled with distilled water which has been degassed to 120mmHg. A bubble is then trapped in the center and driven by a piezoelectric transducer at 32.2kHz attached to the side of the cell. An optical system is then set up to take strobbed video images along and light emission data simultaneously. Temperature, pressure, drive voltage, and listener voltage are also monitored. PMT output in Volts The radii of the bubbles for both experiment s are fit using the Rayleigh-Plesset equation and the acoustic drive amplitude and the ambient bubble radius are found. There is little change in the acoustic drive amplitude as we expect, since we are not varying the drive voltage. However. the ambient bubble radius goes up considerably. These changes (increased light output, increased maximum bubble radius, and increased ambient bubble radius) are also observed when the ambient pressure is varied in the laboratory by an amount similar to that due to gravitation. The changes in the ambient bubble radius and light output with a change in ambient pressure are predicted by the "dissociation hypothesis" and have been observed by other groups in the laboratory. It seems clear that buoyancy's effect on light output and bubble radius, are at best on the same order as the effects of ambient pressure.

Dzikowicz, Ben↗

Experimental and Computational Analysis of Unidirectional Flow Through Stirling Engine Heater Head

A high efficiency Stirling Radioisotope Generator (SRG) is being developed for possible use in long-duration space science missions. NASA s advanced technology goals for next generation Stirling convertors include increasing the Carnot efficiency and percent of Carnot efficiency. To help achieve these goals, a multi-dimensional Computational Fluid Dynamics (CFD) code is being developed to numerically model unsteady fluid flow and heat transfer phenomena of the oscillating working gas inside Stirling convertors. In the absence of transient pressure drop data for the zero mean oscillating multi-dimensional flows present in the Technology Demonstration Convertors on test at NASA Glenn Research Center, unidirectional flow pressure drop test data is used to compare against 2D and 3D computational solutions. This study focuses on tracking pressure drop and mass flow rate data for unidirectional flow though a Stirling heater head using a commercial CFD code (CFD-ACE). The commercial CFD code uses a porous-media model which is dependent on permeability and the inertial coefficient present in the linear and nonlinear terms of the Darcy-Forchheimer equation. Permeability and inertial coefficient were calculated from unidirectional flow test data. CFD simulations of the unidirectional flow test were validated using the porous-media model input parameters which increased simulation accuracy by 14 percent on average.

Wilson, Scott D.↗

On Problems Associated with Modeling Wing-Tail Configurations from Wind Tunnel Data

This paper considers factors that contribute to poor identification of unsteady aerodynamics from wind tunnel data for an airliner configuration. One approach to modeling a wing-tail configuration is considered and applied to both steady and large-amplitude forced pitch oscillation wind tunnel data taken over a wide range of angles of attack but a limited range of amplitude and frequencies. The identified models fit the measured data well but in some cases with inaccurate parameters. Only limited conclusions can be drawn from analysis of the current data set until further experiments can be performed to resolve the identification issues. The analysis of measured and simulated data provides some insights and guidance on how an effective experiment may be designed for wing-tail configurations with nonlinear unsteady aerodynamics.

Murphy, Patrick C.↗

Pulse control of single degree-of-freedom system

The possibility of using active control for structural and mechanical systems is examined. Open loop adaptive control of structural systems subjected to deterministic and stochastic excitations is considered by means of an algorithm for a single-degree-of-freedom oscillator. The algorithm is heuristic and requires continuous monitoring of the state variables, with a pulse control applied when some prespecified level of control is exceeded. The control algorithm is based on the principle that the gradual rhythmic increase of the structural response can be eliminated by applying a pulse of suitable magnitude in the proper direction, thereby avoiding the application of large forces over sustained periods of time. The optimum pulse magnitude is determined as the minimum rms response of the system. Application to linear and nonlinear systems, as well as in chatter suppression, are outlined.

Udwadia, F. E.↗

Finite element nonlinear panel flutter with arbitrary temperatures in supersonic flow

A finite element frequency domain method for predicting nonlinear flutter response of panels with temperature effects is presented. By using the principle of virtual work, the element nonlinear stiffness formulation for a panel under a combined thermal and aerodynamic loads is derived on the bases of von Karman's large deflection plate theory, the first-order piston theory aerodynamics and the quasi-steady thermal stress theory. The system equations of motion can be mathematically separated into two sets of equations and then solved in sequence. The first set of equations yields the panel thermal-aerodynamic equilibrium and the second set of equations of motion leads to the flutter limit-cycle oscillations. Stability and flutter boundaries can also be obtained from the two sets of system equations. Finite element large amplitude limit-cycle flutter results at different uniform temperatures are obtained for a simply supported square panel and are compared with existing Galerkin/time integration and other finite element solutions. Effects of nonuniform temperature distributions, panel length-to-width ratios, and boundary conditions on flutter responses of rectangular and triangular panels are presented.

Xue, David Y.↗

High Order Difference Method for Low Mach Number Aeroacoustics

A high order finite difference method with improved accuracy and stability properties for computational aeroacoustics (CAA) at low Mach numbers is proposed. The Euler equations are split into a conservative and a symmetric non- conservative portion to allow the derivation of a generalized energy estimate. Since the symmetrization is based on entropy variables, that splitting of the flux derivatives is referred to as entropy splitting. Its discretization by high order central differences was found to need less numerical dissipation than conventional conservative schemes. Owing to the large disparity of acoustic and stagnation quantities in low Mach number aeroacoustics, the split Euler equations are formulated in perturbation form. The unknowns are the small changes of the conservative variables with respect to their large stagnation values. All nonlinearities and the conservation form of the conservative portion of the split flux derivatives can be retained, while cancellation errors are avoided with its discretization opposed to the conventional conservative form. The finite difference method is third-order accurate at the boundary and the conventional central sixth-order accurate stencil in the interior. The difference operator satisfies the summation by parts property analogous to the integration by parts in the continuous energy estimate. Thus, strict stability of the difference method follows automatically. Spurious high frequency oscillations are suppressed by a characteristic-based filter similar to but without limiter. The time derivative is approximated by a 4-stage low-storage second-order explicit Runge-Kutta method. The method has been applied to simulate vortex sound at low Mach numbers. We consider the Kirchhoff vortex, which is an elliptical patch of constant vorticity rotating with constant angular frequency in irrotational flow. The acoustic pressure generated by the Kirchhoff vortex is governed by the 2D Helmholtz equation, which can be solved analytically using separation of variables.

Mueller, B.↗