Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear equation systems”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

A computational algorithm for spacecraft control and momentum management

Developments in the area of nonlinear control theory have shown how coordinate changes in the state and input spaces of a dynamical system can be used to transform certain nonlinear differential equations into equivalent linear equations. These techniques are applied to the control of a spacecraft equipped with momentum exchange devices. An optimal control problem is formulated that incorporates a nonlinear spacecraft model. An algorithm is developed for solving the optimization problem using feedback linearization to transform to an equivalent problem involving a linear dynamical constraint and a functional approximation technique to solve for the linear dynamics in terms of the control. The original problem is transformed into an unconstrained nonlinear quadratic program that yields an approximate solution to the original problem. Two examples are presented to illustrate the results.

Dzielski, John↗
NASA NTRS1993

Optimal fixed-finite-dimensional compensator for Burgers' equation with unbounded input/output operators

The problem of using reduced order dynamic compensators to control a class of nonlinear parabolic distributed parameter systems was considered. Concentration was on a system with unbounded input and output operators governed by Burgers' equation. A linearized model was used to compute low-order-finite-dimensional control laws by minimizing certain energy functionals. Then these laws were applied to the nonlinear model. Standard approaches to this problem employ model/controller reduction techniques in conjunction with linear quadratic Gaussian (LQG) theory. The approach used is based on the finite dimensional Bernstein/Hyland optimal projection theory which yields a fixed-finite-order controller.

Burns, John A.↗
NASA NTRS1994

Magnetotail dynamics under isobaric constraints

Using linear theory and nonlinear MHD simulations, we investigate the resistive and ideal MHD stability of two-dimensional plasma configurations under the isobaric constraint dP/dt = 0, which in ideal MHD is equivalent to conserving the pressure function P = P(A), where A denotes the magnetic flux. This constraint is satisfied for incompressible modes, such as Alfven waves, and for systems undergoing energy losses. The linear stability analysis leads to a Schroedinger equation, which can be investigated by standard quantum mechanics procedures. We present an application to a typical stretched magnetotail configuration. For a one-dimensional sheet equilibrium characteristic properties of tearing instability are rediscovered. However, the maximum growth rate scales with the 1/7 power of the resistivity, which implies much faster growth than for the standard tearing mode (assuming that the resistivity is small). The same basic eigen-mode is found also for weakly two-dimensional equilibria, even in the ideal MHD limit. In this case the growth rate scales with the 1/4 power of the normal magnetic field. The results of the linear stability analysis are confirmed qualitatively by nonlinear dynamic MHD simulations. These results suggest the interesting possibility that substorm onset, or the thinning in the late growth phase, is caused by the release of a thermodynamic constraint without the (immediate) necessity of releasing the ideal MHD constraint. In the nonlinear regime the resistive and ideal developments differ in that the ideal mode does not lead to neutral line formation without the further release of the ideal MHD constraint; instead a thin current sheet forms. The isobaric constraint is critically discussed. Under perhaps more realistic adiabatic conditions the ideal mode appears to be stable but could be driven by external perturbations and thus generate the thin current sheet in the late growth phase, before a nonideal instability sets in.

Birn, Joachim↗
NASA NTRS2000

Bending Boundary Layers in Laminated-Composite Circular Cylindrical Shells

An analytical, parametric study of the attenuation of bending boundary layers or edge effects in balanced and unbalanced, symmetrically and unsymmetrically laminated thin cylindrical shells is presented for nine contemporary material systems. The analysis is based on the linear Sanders-Koiter shell equations and specializations to the Love-Kirchhoff shell equations and Donnell's equations are included. Two nondimensional parameters are identified that characterize and quantify the effects of laminate orthotropy and laminate anisotropy on the bending boundary-layer decay length in a very general and encompassing manner. A substantial number of structural design technology results are presented for a wide range of laminated-composite cylinders. For all the laminate constructions considered, the results show that the differences between results that were obtained with the Sanders-Koiter shell equations, the Love-Kirchhoff shell equations, and Donnell's equations are negligible. The results also show that the effect of anisotropy in the form of coupling between pure bending and twisting has a negligible effect on the size of the bending boundary-layer decay length of the balanced, symmetrically laminated cylinders considered. Moreover, the results show that coupling between the various types of shell anisotropies has a negligible effect on the calculation of the bending boundary-layer decay length in most cases. The results also show that in some cases neglecting the shell anisotropy results in underestimating the bending boundary-layer decay length and in other cases it results in an overestimation.

Nemeth, Michael P.↗
NASA NTRS2000

Bending Boundary Layers in Laminated-Composite Circular Cylindrical Shells

A study of the attenuation of bending boundary layers in balanced and unbalanced, symmetrically and unsymmetrically laminated cylindrical shells is presented for nine contemporary material systems. The analysis is based on the linear Sanders-Koiter shell equations and specializations to the Love-Kirchhoff shell equations and Donnell's equations are included. Two nondimensional parameters are identified that characterize the effects of laminate orthotropy and anisotropy on the bending boundary-layer decay length in a very general manner. A substantial number of structural design technology results are presented for a wide range of laminated-composite cylinders. For all laminates considered, the results show that the differences between results obtained with the Sanders-Koiter shell equations, the Love-Kirchhoff shell equations, and Donnell's equations are negligible. The results also show that the effect of anisotropy in the form of coupling between pure bending and twisting has a negligible effect on the size of the bending boundary-layer decay length of the balanced, symmetrically laminated cylinders considered. Moreover, the results show that coupling between the various types of shell anisotropies has a negligible effect on the calculation of the bending boundary-layer decay length in most cases. The results also show that, in some cases, neglecting the shell anisotropy results in underestimating the bending boundary-layer decay length and, in other cases, results in an overestimation.

Nemeth, Michael P.↗
NASA NTRS2001

Multigrid approaches to non-linear diffusion problems on unstructured meshes

The efficiency of three multigrid methods for solving highly non-linear diffusion problems on two-dimensional unstructured meshes is examined. The three multigrid methods differ mainly in the manner in which the nonlinearities of the governing equations are handled. These comprise a non-linear full approximation storage (FAS) multigrid method which is used to solve the non-linear equations directly, a linear multigrid method which is used to solve the linear system arising from a Newton linearization of the non-linear system, and a hybrid scheme which is based on a non-linear FAS multigrid scheme, but employs a linear solver on each level as a smoother. Results indicate that all methods are equally effective at converging the non-linear residual in a given number of grid sweeps, but that the linear solver is more efficient in cpu time due to the lower cost of linear versus non-linear grid sweeps.

Mavriplis, Dimitri J.↗
NASA NTRS2001

An Experimental Evaluation of Generalized Predictive Control for Tiltrotor Aeroelastic Stability Augmentation in Airplane Mode of Flight

The results of a joint NASA/Army/Bell Helicopter Textron wind-tunnel test to assess the potential of Generalized Predictive Control (GPC) for actively controlling the swashplate of tiltrotor aircraft to enhance aeroelastic stability in the airplane mode of flight are presented. GPC is an adaptive time-domain predictive control method that uses a linear difference equation to describe the input-output relationship of the system and to design the controller. The test was conducted in the Langley Transonic Dynamics Tunnel using an unpowered 1/5-scale semispan aeroelastic model of the V-22 that was modified to incorporate a GPC-based multi-input multi-output control algorithm to individually control each of the three swashplate actuators. Wing responses were used for feedback. The GPC-based control system was highly effective in increasing the stability of the critical wing mode for all of the conditions tested, without measurable degradation of the damping in the other modes. The algorithm was also robust with respect to its performance in adjusting to rapid changes in both the rotor speed and the tunnel airspeed.

Kvaternik, Raymond G.↗
NASA NTRS1988

On-orbit parametric identification methodology

On-orbit system identification (ID) of large space systems is essential for various reasons. For example, the complex composite structure of such systems cannot be ground-tested; their structural dynamic characteristics must be known accurately in order to accomplish active control. Furthermore, such capability can be used to characterize/identify various disturbances. The identification process is consisted of four principal elements: (1) modeling, (2) the estimation algorithm, (3) input system, and (4) measurement system. These elements are highly correlated and all togerher determine the success of the identification problem. Accurate modeling of large space systems is the most important element of the identification process. Large flexible structures are non-linear and infinite dimensional systems with highly coupled parameters and low frequency packed modes. In addition, these systems are subject to stochastic and time-varying disturbances, they have structural parameters which can vary due to on-orbit assembly deployment, and operations. These systems are generally; however, represented by constant coefficient, finite order differential equations. The non-linearities, coupling and noise effects are also often neglected. Moreover, identification experiment designs which lead to highly complex optimization problems usually require the simultaneous choice of ID algorithm, sensor, and actuator type and placement. On-orbit bandwidth and power restrictions on excitation, limited data window, and restrictions on sensor/actuator type, placement and number, has led to practical questions of implementations.

Hadaegh, Fred Y.↗
NASA NTRS1990

Linear iterative solvers for implicit ODE methods

The numerical solution of stiff initial value problems, which lead to the problem of solving large systems of mildly nonlinear equations are considered. For many problems derived from engineering and science, a solution is possible only with methods derived from iterative linear equation solvers. A common approach to solving the nonlinear equations is to employ an approximate solution obtained from an explicit method. The error is examined to determine how it is distributed among the stiff and non-stiff components, which bears on the choice of an iterative method. The conclusion is that error is (roughly) uniformly distributed, a fact that suggests the Chebyshev method (and the accompanying Manteuffel adaptive parameter algorithm). This method is described, also commenting on Richardson's method and its advantages for large problems. Richardson's method and the Chebyshev method with the Mantueffel algorithm are applied to the solution of the nonlinear equations by Newton's method.

Saylor, Paul E.↗
NASA NTRS1984

Discrete approximation methods for parameter identification in delay systems

Approximation schemes for parameter identification problems in which the governing state equation is a linear functional differential equation of retarded type are constructed. The basis of the schemes is the replacement of the parameter identification problem having an infinite dimensional state equation by a sequence of approximating parameter identification problems in which the states are given by finite dimensional discrete difference equations. The difference equations are constructed using linear semigroup theory and rational function approximations to the exponential. Sufficient conditions are given for the convergence of solutions to the approximating problems, which can be obtained using conventional methods, to solutions to the original parameter identification problem. Finite difference and spline based schemes using Paderational function approximations to the exponential are constructed, and shown to satisfy the sufficient conditions for convergence. A discussion and analysis of numerical results obtained through the application of the schemes to several examples is included.

Rosen, I. G.↗
NASA NTRS1961

An Analytical Approach to the Design of an Automatic Discontinuous Control System

An analytical investigation was conducted concerning the design of an attitude-stabilization system for stabilizing a vehicle experiencing negligible external moments. The system studied was an automatic discontinuous control system employing a linear switching function including effects of pure time delays, rise and decay time, and neutral zone. Equations were developed which generalize the transient and limit-cycle performance of the control system. understanding of how the physical constants of the system affect its performance, the equations enable the optimization of the system with regard to most considerations that can be expressed mathematically. Design charts are presented which enable rapid calculation of the best thrust level and switching-function coefficient for minimizing power required to stabilize the vehicle within some amplitude, minimizing attitude error, and minimizing angular-velocity error after a period of operation.

Taylor, Lawrence W., Jr.↗
NASA NTRS1991

The RAPT control method for nonlinear systems

The Rapid Advancement Preview Tracking (RAPT) control method uses both a nonlinear system model and a linearized system model to provide a form of feedback/preview control for nonlinear systems. It uses measurements of the states to calculate the current preview control signal from a local linearization of the state equations, and then iteratively calculates a perturbation correction for that control using the nonlinear equations. The corrected control signal is then applied to the system, driving the end-effector to the next position along a specified trajectory. Computer simulations of tracking control for two-link planar manipulators indicate that this algorithm may be fast enough for real-time control.

Long, T. W.↗
NASA NTRS2003

Control Law Design in a Computational Aeroelasticity Environment

A methodology for designing active control laws in a computational aeroelasticity environment is given. The methodology involves employing a systems identification technique to develop an explicit state-space model for control law design from the output of a computational aeroelasticity code. The particular computational aeroelasticity code employed in this paper solves the transonic small disturbance aerodynamic equation using a time-accurate, finite-difference scheme. Linear structural dynamics equations are integrated simultaneously with the computational fluid dynamics equations to determine the time responses of the structure. These structural responses are employed as the input to a modern systems identification technique that determines the Markov parameters of an "equivalent linear system". The Eigensystem Realization Algorithm is then employed to develop an explicit state-space model of the equivalent linear system. The Linear Quadratic Guassian control law design technique is employed to design a control law. The computational aeroelasticity code is modified to accept control laws and perform closed-loop simulations. Flutter control of a rectangular wing model is chosen to demonstrate the methodology. Various cases are used to illustrate the usefulness of the methodology as the nonlinearity of the aeroelastic system is increased through increased angle-of-attack changes.

Newsom, Jerry R.↗
NASA NTRS2011

Development and Verification of the Charring, Ablating Thermal Protection Implicit System Simulator

The development and verification of the Charring Ablating Thermal Protection Implicit System Solver (CATPISS) is presented. This work concentrates on the derivation and verification of the stationary grid terms in the equations that govern three-dimensional heat and mass transfer for charring thermal protection systems including pyrolysis gas flow through the porous char layer. The governing equations are discretized according to the Galerkin finite element method (FEM) with first and second order fully implicit time integrators. The governing equations are fully coupled and are solved in parallel via Newton s method, while the linear system is solved via the Generalized Minimum Residual method (GMRES). Verification results from exact solutions and Method of Manufactured Solutions (MMS) are presented to show spatial and temporal orders of accuracy as well as nonlinear convergence rates.

Amar, Adam J.↗
NASA NTRS1989

An optimal momentum management controller for the Space Station

This paper presents a new sequential design procedure for determining an optimal CMG (control moment gyro) momentum management and attitude control system for the Space Station Freedom. First, the Space Station equations of motion are linearized and uncoupled, and the associated state space equations are defined. Next, a new sequential procedure is used for the development of a continuous quadratic regulator with eigenvalue placement in a specified region of the complex plane. The regional pole assignment method is utilized since it is best suited for tradeoffs between eigenvalue locations and robustness with respect to parameter variations, sensor failures, implementation accuracies and gain reductions. The matrix sign function is used for solving the Riccati equations which appear in the design procedure. Simulation results are given which show that the resultant design provides desired system performance.

Sunkel, J. W.↗
NASA NTRS1972

Numerical method for the solution of large systems of differential equations of the boundary layer type

A numerical method for the solution of large systems of nonlinear differential equations of the boundary-layer type is described. The method is a modification of the technique for satisfying asymptotic boundary conditions. The present method employs inverse interpolation instead of the Newton method to adjust the initial conditions of the related initial-value problem. This eliminates the so-called perturbation equations. The elimination of the perturbation equations not only reduces the user's preliminary work in the application of the method, but also reduces the number of time-consuming initial-value problems to be numerically solved at each iteration. For further ease of application, the solution of the overdetermined system for the unknown initial conditions is obtained automatically by applying Golub's linear least-squares algorithm. The relative ease of application of the proposed numerical method increases directly as the order of the differential-equation system increases. Hence, the method is especially attractive for the solution of large-order systems. After the method is described, it is applied to a fifth-order problem from boundary-layer theory.

Green, M. J.↗
NASA NTRS1999

Developing a Fundamental Model for an Integrated GPS/INS State Estimation System with Kalman Filtering

This work will demonstrate the integration of sensor and system dynamic data and their appropriate models using an optimal filter to create a robust, adaptable, easily reconfigurable state (motion) estimation system. This state estimation system will clearly show the application of fundamental modeling and filtering techniques. These techniques are presented at a general, first principles level, that can easily be adapted to specific applications. An example of such an application is demonstrated through the development of an integrated GPS/INS navigation system. This system acquires both global position data and inertial body data, to provide optimal estimates of current position and attitude states. The optimal states are estimated using a Kalman filter. The state estimation system will include appropriate error models for the measurement hardware. The results of this work will lead to the development of a "black-box" state estimation system that supplies current motion information (position and attitude states) that can be used to carry out guidance and control strategies. This black-box state estimation system is developed independent of the vehicle dynamics and therefore is directly applicable to a variety of vehicles. Issues in system modeling and application of Kalman filtering techniques are investigated and presented. These issues include linearized models of equations of state, models of the measurement sensors, and appropriate application and parameter setting (tuning) of the Kalman filter. The general model and subsequent algorithm is developed in Matlab for numerical testing. The results of this system are demonstrated through application to data from the X-33 Michael's 9A8 mission and are presented in plots and simple animations.

Canfield, Stephen↗
NASA NTRS1991

Implementation of a partitioned algorithm for simulation of large CSI problems

The implementation of a partitioned numerical algorithm for determining the dynamic response of coupled structure/controller/estimator finite-dimensional systems is reviewed. The partitioned approach leads to a set of coupled first and second-order linear differential equations which are numerically integrated with extrapolation and implicit step methods. The present software implementation, ACSIS, utilizes parallel processing techniques at various levels to optimize performance on a shared-memory concurrent/vector processing system. A general procedure for the design of controller and filter gains is also implemented, which utilizes the vibration characteristics of the structure to be solved. Also presented are: example problems; a user's guide to the software; the procedures and algorithm scripts; a stability analysis for the algorithm; and the source code for the parallel implementation.

Alvin, Kenneth F.↗