SEARCH · Search NASA

Results for “Linear and nonlinear operators”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

NASA NTRS1991

A linear control design structure to maintain loop properties during limit operation in a multi-nozzle turbofan engine

The implementation of multi-variable control systems on turbofan engines requires the use of limit protection to maintain safe engine operation. Since a turbofan engine typically encounters limits during transient operation, the use of a limit protection scheme that modifies the feedback loop may void the desired 'guarantees' associated with linear multi-variable control design methods, necessitating considerable simulation to validate the control with limit protection. An alternative control design structure is proposed that maintains the desired linear feedback properties when certain safety limits are encountered by moving the limit protection scheme outside of the feedback loop. This proposed structure is compared to a structure with a limit protection scheme that modifies the feedback loop properties. The two design structures are compared using both linear and nonlinear simulations. The evaluation emphasizes responses where the fan surge margin limit is encountered.

Mattern, Duane↗

NASA NTRS1993

A comparison of two multi-variable integrator windup protection schemes

Two methods are examined for limit and integrator wind-up protection for multi-input, multi-output linear controllers subject to actuator constraints. The methods begin with an existing linear controller that satisfies the specifications for the nominal, small perturbation, linear model of the plant. The controllers are formulated to include an additional contribution to the state derivative calculations. The first method to be examined is the multi-variable version of the single-input, single-output, high gain, Conventional Anti-Windup (CAW) scheme. Except for the actuator limits, the CAW scheme is linear. The second scheme to be examined, denoted the Modified Anti-Windup (MAW) scheme, uses a scalar to modify the magnitude of the controller output vector while maintaining the vector direction. The calculation of the scalar modifier is a nonlinear function of the controller outputs and the actuator limits. In both cases the constrained actuator is tracked. These two integrator windup protection methods are demonstrated on a turbofan engine control system with five measurements, four control variables, and four actuators. The closed-loop responses of the two schemes are compared and contrasted during limit operation. The issue of maintaining the direction of the controller output vector using the Modified Anti-Windup scheme is discussed and the advantages and disadvantages of both of the IWP methods are presented.

Mattern, Duane↗

NASA NTRS1961

Study of Systems Using Inertia Wheels for Precise Attitude Control of a Satellite

Systems using inertia wheels are evaluated in this report to determine their suitability for precise attitude control of a satellite and to select superior system configurations. Various possible inertia wheel system configurations are first discussed in a general manner. Three of these systems which appear more promising than the others are analyzed in detail, using the Orbiting Astronomical Observatory as an example. The three systems differ from each other only by the method of damping, which is provided by either a rate gyro, an error-rate network, or a tachometer in series with a high-pass filter. An analytical investigation which consists of a generalized linear analysis, a nonlinear analysis using the switching-time method, and an analog computer study shows that all three systems are theoretically capable of producing adequate response and also of maintaining the required pointing accuracy for the Orbiting Astronomical Observatory of plus or minus 0.1 second of arc. Practical considerations and an experimental investigation show, however, that the system which uses an error-rate network to provide damping is superior to the other two systems. The system which uses a rate gyro is shown to be inferior because the threshold level causes a significant amount of limit-cycle operation, and the system which uses a tachometer with a filter is shown to be inferior because a device with the required dynamic range of operation does not appear to be available. The experimental laboratory apparatus used to investigate the dynamic performance of the systems is described, and experimental results are included to show that under laboratory conditions with relatively large extraneous disturbances, a dynamic tracking error of less than plus or minus 0.5 second of arc was obtained.

White, John S.↗

NASA NTRS2005

High-Order Tunable Filters Based on a Chain of Coupled Crystalline Whispering Gallery-Mode Resonators

We demonstrate experimentally a tunable third-order optical filter fabricated from the three voltage-controlled lithium niobate whispering gallery-mode resonators. The filter operates at 1550 nm with 30-MHz bandwidth and can be electrooptically tuned by 12 GHz in the linear regime with approximately 80-MHz/V tuning rate. With this filter, we have demonstrated 6-dB fiber-to-fiber insertion loss and 30-ns tuning speed, limited by the resonator buildup time.

optical resonators↗

NASA NTRS1976

A precision star tracker utilizing advanced techniques and materials

An image dissector star tracker has been developed which operates in the photon counting mode making it possible to utilize all digital electronics. A unique pulse processing circuit allows bright stars to be tracked as well as dim stars. Thermal mechanical stability has been greatly enhanced by fabricating a housing with graphite/epoxy composite material with a linear coefficient of thermal expansion near zero. Test results indicate the +10 Mv stars can be acquired and tracked, while position variation with star intensity is less than 2 arc seconds from 2.5 Mv to +10 Mv. The noise equivalent angle for a +8 Mv star is 3 arc seconds. Polynominal correction for remaining cross-coupling and nonlinearity reduced error over total 1 x 1 deg field to 1.5 arc seconds.

Gates, R. F.↗

NASA NTRS2003

Numerical and Experimental Determination of the Geometric Far Field for Round Jets

To reduce ambiguity in the reporting of far field jet noise, three round jets operating at subsonic conditions have recently been studied at the NASA Glenn Research Center. The goal of the investigation was to determine the location of the geometric far field both numerically and experimentally. The combination of the WIND Reynolds-Averaged Navier-Stokes solver and the MGBK jet noise prediction code was used for the computations, and the experimental data was collected in the Aeroacoustic Propulsion Laboratory. While noise sources are distributed throughout the jet plume, at great distances from the nozzle the noise will appear to be emanating from a point source and the assumption of linear propagation is valid. Closer to the jet, nonlinear propagation may be a problem, along with the known geometric issues. By comparing sound spectra at different distances from the jet, both from computational methods that assume linear propagation, and from experiments, the contributions of geometry and nonlinearity can be separately ascertained and the required measurement distance for valid experiments can be established. It is found that while the shortest arc considered here (approx. 8D) was already in the geometric far field for the high frequency sound (St greater than 2.0), the low frequency noise due to its extended source distribution reached the geometric far field at or about 50D. It is also found that sound spectra at far downstream angles does not strictly scale on Strouhal number, an observation that current modeling does not capture.

Koch, L. Danielle↗

NASA NTRS1986

Aircraft control surface failure detection and isolation using the OSGLR test

The performance of the orthogonal series generalized likelihood ratio (OSGLR) test in detecting and isolating commercial aircraft control surface and actuator failures is evaluated. A modification to incorporate age-weighting which significantly reduces the sensitivity of the algorithm to modeling errors is presented. The steady-state implementation of the algorithm based on a single linear model valid for a cruise flight condition is tested using a nonlinear aircraft simulation. A number of off-nominal no-failure flight conditions including maneuvers, nonzero flap deflections, different turbulence levels and steady winds were tested. Based on the no-failure decision functions produced by off-nominal flight conditions, the failure detection and isolation performance at the nominal flight condition was determined. The extension of the algorithm to a wider flight envelope by scheduling on dynamic pressure and flap deflection is examined. Based on this testing, the OSGLR algorithm should be capable of detecting control surface failures that would affect the safe operation of a commercial aircraft. Isolation may be difficult if there are several surfaces which produce similar effects on the aircraft. Extending the algorithm over the entire operating envelope of a commercial aircraft appears feasible.

Bonnice, W. F.↗

NASA NTRS2002

Reentry Vehicle Flight Controls Design Guidelines: Dynamic Inversion

This report addresses issues in developing a flight control design for vehicles operating across a broad flight regime and with highly nonlinear physical descriptions of motion. Specifically it addresses the need for reentry vehicles that could operate through reentry from space to controlled touchdown on Earth. The latter part of controlled descent is achieved by parachute or paraglider - or by all automatic or a human-controlled landing similar to that of the Orbiter. Since this report addresses the specific needs of human-carrying (not necessarily piloted) reentry vehicles, it deals with highly nonlinear equations of motion, and then-generated control systems must be robust across a very wide range of physics. Thus, this report deals almost exclusively with some form of dynamic inversion (DI). Two vital aspects of control theory - noninteracting control laws and the transformation of nonlinear systems into equivalent linear systems - are embodied in DI. Though there is no doubt that the mathematical tools and underlying theory are widely available, there are open issues as to the practicality of using DI as the only or primary design approach for reentry articles. This report provides a set of guidelines that can be used to determine the practical usefulness of the technique.

Ito, Daigoro↗

NASA NTRS1992

Correction for nonlinear photon counting effects in lidar systems

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.

Donovan, D. P.↗

NASA NTRS2017

Novel Estimation of Pilot Performance Characteristics

Two mechanisms internal to the pilot that affect performance during a tracking task are: 1) Pilot equalization (i.e. lead/lag); and 2) Pilot gain (i.e. sensitivity to the error signal). For some applications McRuer's Crossover Model can be used to anticipate what equalization will be employed to control a vehicle's dynamics. McRuer also established approximate time delays associated with different types of equalization - the more cognitive processing that is required due to equalization difficulty, the larger the time delay. However, the Crossover Model does not predict what the pilot gain will be. A nonlinear pilot control technique, observed and coined by the authors as 'amplitude clipping', is shown to improve stability, performance, and reduce workload when employed with vehicle dynamics that require high lead compensation by the pilot. Combining linear and nonlinear methods a novel approach is used to measure the pilot control parameters when amplitude clipping is present, allowing precise measurement in real time of key pilot control parameters. Based on the results of an experiment which was designed to probe workload primary drivers, a method is developed that estimates pilot spare capacity from readily observable measures and is tested for generality using multi-axis flight data. This paper documents the initial steps to developing a novel, simple objective metric for assessing pilot workload and its variation over time across a wide variety of tasks. Additionally, it offers a tangible, easily implementable methodology for anticipating a pilot's operating parameters and workload, and an effective design tool. The model shows promise in being able to precisely predict the actual pilot settings and workload, and observed tolerance of pilot parameter variation over the course of operation. Finally, an approach is proposed for generating Cooper-Harper ratings based on the workload and parameter estimation methodology.

Human pilot models↗

NASA NTRS1983

Formulation and solution of rotary-wing aeroelastic stability and response problems

The state of the art in the formulation and solution of rotary-wing aeroelastic stability and response problems is reviewed in detail. The approximations used in the structural, inertia and aerodynamic operators are discussed. The important role of geometric nonlinearities, due to moderate deflections, and aerodynamic stall in the aeroelastic stability and response problem are identified. It is also shown that geometric nonlinearities are of primary importance in aeroelastic stability calculations, and have a more limited, though important, role in response calculations. Next, formulation of coupled rotor/fuselage problems is described, for both air and ground resonance type problems. Both topics, the isolated blade problem and the coupled rotor/fuselage problem, are treated for both hover and forward flight. Solution of aeroelastic stability and response problems proceeds in two stages. First, the spatial dependence is eliminated by using Galerkin's method, or by using the finite element method. Next the nonlinear, or linear, ordinary differential equation with periodic coefficients have to be solved for stability or response. Efficient numerical methods for accomplishing these objectives are presented in a comprehensive manner. The paper contains a number of illustrative numerical results which are intended to clarify various aspects of the modeling process and serve as representative results for both aeroelastic stability and response calculations for a variety of blade and rotor configurations.

Friedmann, P. P.↗

NASA NTRS1998

Parallel Implicit Algorithms for CFD

The main goal of this project was efficient distributed parallel and workstation cluster implementations of Newton-Krylov-Schwarz (NKS) solvers for implicit Computational Fluid Dynamics (CFD.) "Newton" refers to a quadratically convergent nonlinear iteration using gradient information based on the true residual, "Krylov" to an inner linear iteration that accesses the Jacobian matrix only through highly parallelizable sparse matrix-vector products, and "Schwarz" to a domain decomposition form of preconditioning the inner Krylov iterations with primarily neighbor-only exchange of data between the processors. Prior experience has established that Newton-Krylov methods are competitive solvers in the CFD context and that Krylov-Schwarz methods port well to distributed memory computers. The combination of the techniques into Newton-Krylov-Schwarz was implemented on 2D and 3D unstructured Euler codes on the parallel testbeds that used to be at LaRC and on several other parallel computers operated by other agencies or made available by the vendors. Early implementations were made directly in Massively Parallel Integration (MPI) with parallel solvers we adapted from legacy NASA codes and enhanced for full NKS functionality. Later implementations were made in the framework of the PETSC library from Argonne National Laboratory, which now includes pseudo-transient continuation Newton-Krylov-Schwarz solver capability (as a result of demands we made upon PETSC during our early porting experiences). A secondary project pursued with funding from this contract was parallel implicit solvers in acoustics, specifically in the Helmholtz formulation. A 2D acoustic inverse problem has been solved in parallel within the PETSC framework.

Keyes, David E.↗

NASA NTRS1983

Transformations of the perturbed two-body problem to unperturbed harmonic oscillators

Singular, nonlinear, and Liapunov unstable equations are made regular and linear through transformations that change the perturbed planar problem of two bodies into unperturbed and undamped harmonic oscillators with constant coefficients, so that the stable solution may be immediately written in terms of the new variables. The use of arbitrary and special functions for the transformations allows the systematic discussion of previously introduced and novel anomalies. For the case of the unperturbed two-body problem, it is proved that if transformations are power functions of the radial variable, only the eccentric and the true anomalies (with the corresponding transformations of the radial variable) will result in harmonic oscillators. The present method significantly reduces computation requirements in autonomous space operations.

Szebehely, V.↗

NASA NTRS1999

A Hysteresis Model for Piezoceramic Materials

This paper addresses the modeling of nonlinear constitutive relations and hysteresis inherent to piezoceramic materials at moderate to high drive levels. Such models are, necessary to realize the, full potential of the materials in high performance control applications, and a necessary prerequisite is the development of techniques which permit control implementation. The approach employed here is based on the qualification of reversible and irreversible domain wall motion in response to applied electric fields. A comparison with experimental data illustrates that because the resulting ODE model is physics-based, it can be employed for both characterization and prediction of polarization levels throughout the range of actuator operation. Finally, the ODE formulation is amenable to inversion which facilitates the development of an inverse compensator for linear control design.

Smith, Ralph C.↗

NASA NTRS2004

A Study of Penalty Function Methods for Constraint Handling with Genetic Algorithm

COMETBOARDS (Comparative Evaluation Testbed of Optimization and Analysis Routines for Design of Structures) is a design optimization test bed that can evaluate the performance of several different optimization algorithms. A few of these optimization algorithms are the sequence of unconstrained minimization techniques (SUMT), sequential linear programming (SLP) and the sequential quadratic programming techniques (SQP). A genetic algorithm (GA) is a search technique that is based on the principles of natural selection or "survival of the fittest". Instead of using gradient information, the GA uses the objective function directly in the search. The GA searches the solution space by maintaining a population of potential solutions. Then, using evolving operations such as recombination, mutation and selection, the GA creates successive generations of solutions that will evolve and take on the positive characteristics of their parents and thus gradually approach optimal or near-optimal solutions. By using the objective function directly in the search, genetic algorithms can be effectively applied in non-convex, highly nonlinear, complex problems. The genetic algorithm is not guaranteed to find the global optimum, but it is less likely to get trapped at a local optimum than traditional gradient-based search methods when the objective function is not smooth and generally well behaved. The purpose of this research is to assist in the integration of genetic algorithm (GA) into COMETBOARDS. COMETBOARDS cast the design of structures as a constrained nonlinear optimization problem. One method used to solve constrained optimization problem with a GA to convert the constrained optimization problem into an unconstrained optimization problem by developing a penalty function that penalizes infeasible solutions. There have been several suggested penalty function in the literature each with there own strengths and weaknesses. A statistical analysis of some suggested penalty functions is performed in this study. Also, a response surface approach to robust design is used to develop a new penalty function approach. This new penalty function approach is then compared with the other existing penalty functions.

Ortiz, Francisco↗

NASA NTRS2017

Parallel Computation of the Jacobian Matrix for Nonlinear Equation Solvers Using MATLAB

Demonstrating speedup for parallel code on a multicore shared memory PC can be challenging in MATLAB due to underlying parallel operations that are often opaque to the user. This can limit potential for improvement of serial code even for the so-called embarrassingly parallel applications. One such application is the computation of the Jacobian matrix inherent to most nonlinear equation solvers. Computation of this matrix represents the primary bottleneck in nonlinear solver speed such that commercial finite element (FE) and multi-body-dynamic (MBD) codes attempt to minimize computations. A timing study using MATLAB's Parallel Computing Toolbox was performed for numerical computation of the Jacobian. Several approaches for implementing parallel code were investigated while only the single program multiple data (spmd) method using composite objects provided positive results. Parallel code speedup is demonstrated but the goal of linear speedup through the addition of processors was not achieved due to PC architecture.

Rose, Geoffrey K.↗

NASA NTRS1998

Efficient and Optimal Attitude Determination Using Recursive Global Positioning System Signal Operations

In this paper, a new and efficient algorithm is developed for attitude determination from Global Positioning System signals. The new algorithm is derived from a generalized nonlinear predictive filter for nonlinear systems. This uses a one time-step ahead approach to propagate a simple kinematics model for attitude determination. The advantages of the new algorithm over previously developed methods include: it provides optimal attitudes even for coplanar baseline configurations; it guarantees convergence even for poor initial conditions; it is a non-iterative algorithm; and it is computationally efficient. These advantages clearly make the new algorithm well suited to on-board applications. The performance of the new algorithm is tested on a dynamic hardware simulator. Results indicate that the new algorithm accurately estimates the attitude of a moving vehicle, and provides robust attitude estimates even when other methods, such as a linearized least-squares approach, fail due to poor initial starting conditions.

Crassidis, John L.↗

NASA NTRS1985

Shortwave infrared 512 x 2 line sensor for earth resources applications

As part of the NASA remote-sensing Multispectral Linear Array Program, an edge-buttable 512 x 2 IRCCD line image sensor with 30-micron Pd2Si Schottky-barrier detectors is developed for operation with passive cooling at 120 K in the 1.1-2.5 micron short infrared band. On-chip CCD multiplexers provide one video output for each 512 detector band. The monolithic silicon line imager performance at a 4-ms optical integration time includes a signal-to-noise ratio of 241 for irradiance of 7.2 microwatts/sq cm at 1.65 microns wavelength, a 5000 dynamic range, a modulation transfer function, greater than 60 percent at the Nyquist frequency, and an 18-milliwatt imager chip total power dissipation. Blemish-free images with three percent nonuniformity under illumination and nonlinearity of 1.25 percent are obtained. A five SWIR imager hybrid focal plane was constructed, demonstrating the feasibility of arrays with only a two-detector loss at each joint.

Tower, J. R.↗