SEARCH · Engineering Papers
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.
A simple mathematical jet engine description is presented where the measured transition functions of the engine thrust can be simulated in a quasi-stationary operation as well as in acceleration or deceleration schedules. Because of its simplicity and high fidelity, it is especially suited for representing jet engines in digital simulation programs.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Photomultiplier tubes (PMT's) employed in the photon counting (PC) mode of operation are widely used as detectors in lidar systems. In our laboratory, we have developed a versatile Nd:YAG lidar which is used for measurement of both the middle atmosphere and the troposphere. With this system, we encounter a very wide range of signal levels ranging from the extremely weak signals from the top of the mesosphere to the very strong returns from low level clouds. Although the system is capable of operating the PMT's in either the analog detection or photon counting mode, we find that often when we use photon counting we have portions of our lidar return which contain very useful information but are not within the linear operating regime of the PC system. We report the results of our efforts to explore the extent to which such high intensity PC signals can be quantitatively analyzed. In particular, a useful model relating the mean 'true' count rate and the observed count rate is presented and it's application to our system demonstrated. This model takes into account the variation in height of the PMT output pulses and the effect of the pulse height discrimination threshold.
An analytical study was made of an automatic reaction-control system for the upper stages of a missile to determine limit cycle characteristics, corresponding duty cycles, and the effects of the various system parameters on these quantities. A nonlinear servo analysis (describing function) technique was used to obtain a mathematical representation of the nonlinear components of the system. The results obtained by this analysis are compared with the results obtained from an analog simulation including the jet reaction-control hardware. The good agreement between the results of the two studies tends to indicate the feasibility of using such an analysis technique in the early design phase of a reaction-control system, since by so doing, the parameters for good system performance can be determined quite readily.
Motion resembling that of a pendulum undergoing large-amplitude limit cycle oscillation was observed during a series of flight tests of an unoccupied Orion Capsule Parachute Assembly System (CPAS) comprised of two parachutes and a capsule payload. Large excursions away from vertical by the capsule could cause it to strike the ground or ocean at a large angle with respect to vertical, or at a large horizontal speed. These conditions are undesirable because they would endanger the occupants of the capsule in an actual mission. A simplified planar dynamics model in conjunction with a nonlinear normal force coefficient vs. angle of attack model serves as the basis of an analytical investigation of the fundamental dynamics of this pendulum motion. Output error methodology from system identification theory was used to identify the parameters of the nonlinear aerodynamics model. The identified model yielded excellent comparison with portions of flight test data where the pendulum motion occurred. Due to the inherent nonlinear nature of the pendulum motion limit cycle, traditional nonlinear analysis techniques were applied to gain further insight into the system. Lyapunov’s direct method provided mathematical proof in the absolute stability of the pendulum mode. Describing Function method was used to predict the amplitude and frequency of the limit cycle oscillation. Finally, phase plane analysis allowed easy visualization on the size and shape of the limit cycle with respect to variations in key aerodynamic parameters.
The theory of Hopf bifurcation with symmetry developed by Golubitsky and Stewart (1985) is applied to systems of ODEs having the symmetries of a regular polygon, that is, whose symmetry group is dihedral. The existence and stability of symmetry-breaking branches of periodic solutions are considered. In particular, these results are applied to a general system of n nonlinear oscillators coupled symmetrically in a ring, and the generic oscillation patterns are described. It is found that the symmetry can force some oscillators to have twice the frequency of others. The case of four oscillators has exceptional features.
Nonlinear models and block diagrams of an inner-loop control system for the ATOPS B-737 Research Aircraft are presented. Continuous time linear model representations of the nonlinear inner-loop control systems are derived. Closed-loop aircraft simulations comparing nonlinear and linear dynamic responses to step inputs are used to verify the inner-loop control system models.
An oscillating system with quadratic damping subjected to white noise excitation is replaced by a nonlinear, statistically equivalent system for which the associated Fokker-Planck equation can be exactly solved. The mean square responses are calculated and the optimum damping coefficient is determined with respect to the minimum mean square acceleration criteria. An application of these results to the optimization of automobile suspension damping is given.
The achievement of quality dynamic performance in manipulator systems is difficult using conventional control methods because of both the inherent geometric nonlinearities of these systems and the dependence of the system dynamics on the characteristics of manipulated objects. A model-referenced adaptive control law is developed for maintaining uniformly good performance over a wide range of motions and payloads. The effectiveness of the approach is demonstrated in several simulations and the system stability as a function of input is investigated. Also developed is a 'learning signal' approach designed to minimize initial transients arising from abrupt changes in the inertial payload.
The world coordinate system (WCS) for dispersion calibrated spectra used in the IRAF/NOAO spectroscopy packages is described. In particular, the image header keywords which define the coordinates in an 'image' pixel array. These keywords appear both as part of the IRAF image structure and map directly to FITS format. The types of spectra include multidimensional images with spectral dispersion as one axis, one dimensional images, and a special 'multispec' format consisting of multiple, semi-independent, one dimensional spectra in two or three dimensional images. The types of coordinate systems include multidimensional linear coordinates with individual keywords and one dimensional linear and nonlinear coordinates for the 'multispec' images stored in IRAF WCS attributes. The nonlinear world coordinate systems include polynomial, spline, sampled table, and look-up table function representations.
This paper discusses the calculation of sensitivities. or derivatives, for optimization problems involving systems governed by differential equations and other state relations. The subject is examined from the point of view of nonlinear programming, beginning with the analytical structure of the first and second derivatives associated with such problems and the relation of these derivatives to implicit differentiation and equality constrained optimization. We also outline an error analysis of the analytical formulae and compare the results with similar results for finite-difference estimates of derivatives. We then attend to an investigation of the nature of the adjoint method and the adjoint equations and their relation to directions of steepest descent. We illustrate the points discussed with an optimization problem in which the variables are the coefficients in a differential operator.
An approximate analytic solution of a set of nonlinear model alpha-omega-dynamo equations is obtained. The reaction of the Lorentz force on the velocity shear which stretches and, hence, amplifies the magnetic field is incorporated into the model. To single out the effect of the Lorentz force on the omega-effect, the effect of the Lorentz force on the alpha-effect is neglected in this study. The solution represents a nonlinear oscillation with the amplitude and period determined by the dynamo number N. The amplitude is proportional to N - 1, while the period is almost exactly the same as the dissipation time of the unstable mode (proportional to N).
A "Multimode" or "switched" system is one that switches between various modes of operation. When a switch occurs from one mode to another, a discontinuity may result followed by a smooth evolution under the new regime. Characterizing the switching behavior of these systems is not well understood and, therefore, identification of multimode systems typically requires a preprocessing step to classify the observed data according to a mode of operation. A further consequence of the switched nature of these systems is that data available for parameter estimation of any subsystem may be inadequate. As such, identification and parameter estimation of multimode systems remains an unresolved problem. In this paper, we 1) show that the NARMAX model structure can be used to describe the impulsive-smooth behavior of switched systems, 2) propose a modified extended least squares (MELS) algorithm to estimate the coefficients of such models, and 3) demonstrate its applicability to simulated and real data from the Vestibulo-Ocular Reflex (VOR). The approach will also allow the identification of other nonlinear bio-systems, suspected of containing "hard" nonlinearities.
Estimation of the state of a nonlinear discrete-time system using quantized data is considered. An exact solution for the maximum likelihood estimate is expressed as the solution of a nonlinear two-point boundary-value problem. Approximate recursive solutions for both the maximum likelihood and the conditional-mean estimates are obtained. The results of Monte-Carlo simulations are presented in which the performance of these two algorithms is compared with that of a Kalman filter in which the quantization error is approximated by white noise.-
A numerical method is presented to determine the steady-state nonlinear response of a rotor-support system due to deadband and rubbing using discrete Fourier transformation and inverse discrete Fourier transformation. Damaging subharmonic and superharmonic responses are found to occur in presence of a side force. The calculated results agree with the general trends which have been observed experimentally by other investigators. The effects of selected nondimensionalized parameters on rotor response are studied.
The Cassini spacecraft dynamics-related telemetry during long Main Engine (ME) burns has indicated the presence of stable limit cycles between 0.03-0.04 Hz frequencies. These stable limit cycles cause the spacecraft to possess non-zero oscillating rates for extended periods of time. This indicates that the linear ME guidance control system does not model the complete dynamics of the spacecraft. In this study, we propose that the observed limit cycles in the spacecraft dynamics telemetry appear from a stable interaction between the unmodeled nonlinear elements in the ME guidance control system. Many nonlinearities in the control system emerge from translating the linear engine gimbal actuator (EGA) motion into a spacecraft rotation. One such nonlinearity comes from the gear backlash in the EGA system, which is the focus of this paper. The limit cycle characteristics and behavior can be predicted by modeling this gear backlash nonlinear element via a describing function and studying the interaction of this describing function with the overall dynamics of the spacecraft. The linear ME guidance controller and gear backlash nonlinearity are modeled analytically. The frequency, magnitude, and nature of the limit cycle are obtained from the frequency response of the ME guidance controller and nonlinear element. In addition, the ME guidance controller along with the nonlinearity is simulated. The simulation response contains a limit cycle with similar characterstics as predicted analytically: 0.03-0.04 Hz frequency and stable, sustained oscillations. The analytical and simulated limit cycle responses are compared to the flight telemetry for long burns such as the Saturn Orbit Insertion and Main Engine Orbit Trim Maneuvers. The analytical and simulated limit cycle characteristics compare well with the actual observed limit cycles in the flight telemetry. Both have frequencies between 0.03-0.04 Hz and stable oscillations. This work shows that the stable limit cycles occur due to the interaction between the unmodeled nonlinear elements and linear ME guidance controller.
The equations of horizontal motion of the neutral atmosphere between 120 and 500 km are integrated with the inclusion of all the nonlinear terms of the convective derivative and the viscous forces due to vertical and horizontal velocity gradients. Empirical models of the distribution of neutral and charged particles are assumed to be known. The model of velocities developed is a steady state model. In part 1 the mathematical method used in the integration of the Navier-Stokes equations is described and the various forces are analysed.