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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

A study of equation solvers for linear and non-linear finite element analysis on parallel processing computers

Concurrent computing environments provide the means to achieve very high performance for finite element analysis of systems, provided the algorithms take advantage of multiple processors. The authors have examined several algorithms for both linear and nonlinear finite element analysis. The performance of these algorithms on an Alliant FX/80 parallel supercomputer has been studied. For single load case linear analysis, the optimal solution algorithm is strongly problem dependent. For multiple load cases or nonlinear analysis through a modified Newton-Raphson method, decomposition algorithms are shown to have a decided advantage over element-by-element preconditioned conjugate gradient algorithms.

Watson, Brian C.↗
NASA NTRS1971

A dynamic model of the human postural control system.

Description of a digital simulation of the pitch axis dynamics of a stick man. The difficulties encountered in linearizing the equations of motion are discussed; the conclusion reached is that a completely linear simulation is of such restricted validity that only a nonlinear simulation is of any practical use. Typical simulation results obtained from the full nonlinear model are illustrated.

Hill, J. C.↗
NASA NTRS1973

The three dimensional motion and stability of a rotating space station-cable - Counterweight configuration

The three dimensional equations of motion for a cable connected space station - counterweight system are developed using a Lagrangian formulation. The system model employed allows for cable and end body damping and restoring effects. The equations are then linearized about the equilibrium motion and nondimensionalized. To first degree, the out-of-plane equations uncouple from the in-plane equations. Therefore, the characteristic polynomials for the in-plane and out-of-plane equations are developed and treated separately. From the general in-plane characteristic equation, necessary conditions for stability are obtained. The Routh-Hurwitz necessary and sufficient conditions for stability are derived for the general out-of-plane characteristic equation. Special cases of the in-plane and out-of-plane equations (such as identical end masses, and when the cable is attached to the centers of mass of the two end bodies) are then examined for stability criteria.

Bainum, P. M.↗
NASA NTRS1974

The three dimensional motion and stability of a rotating space station: Cable-counterweight configuration

The three dimensional equations of motion for a cable connected space station--counterweight system are developed using a Lagrangian formulation. The system model employed allows for cable and end body damping and restoring effects. The equations are then linearized about the equilibrium motion and nondimensionalized. To first degree, the out-of-plane equations uncouple from the inplane equations. Therefore, the characteristic polynomials for the in-plane and out-of-plane equations are developed and treated separately. From the general in-plane characteristic equation, necessary conditions for stability are obtained. The Routh-Hurwitz necessary and sufficient conditions for stability are derived for the general out-of-plane characteristic equation. Special cases of the in-plane and out-of-plane equations (such as identical end masses, and when the cable is attached to the centers of mass of the two end bodies) are then examined for stability criteria.

Bainum, P. M.↗
NASA NTRS1982

Stability of mean flows over an infinite flat plate

The present investigation is mainly concerned with a stability analysis for the linearized Navier-Stokes equations for parallel and nonparallel mean flows over an infinite flat plate. The system of equations for parallel flows is presented. The system is viewed as a generalized Orr-Sommerfeld equation. Attention is given to an explicit criterion characterizing the case when the stability of all physically reasonable solutions is determined by the eigenvalues. The proof given in the investigation is applicable to both the generalized Orr-Somerfeld equations and the modified equations for nonparallel flow. The fact that the criterion is independent of the completeness or incompleteness of eigenfunctions is contrary to some expectations.

Miklavcic, M.↗
NASA NTRS1983

Canonical forms for nonlinear systems

Necessary and sufficient conditions for transforming a nonlinear system to a controllable linear system have been established, and this theory has been applied to the automatic flight control of aircraft. These transformations show that the nonlinearities in a system are often not intrinsic, but are the result of unfortunate choices of coordinates in both state and control variables. Given a nonlinear system (that may not be transformable to a linear system), we construct a canonical form in which much of the nonlinearity is removed from the system. If a system is not transformable to a linear one, then the obstructions to the transformation are obvious in canonical form. If the system can be transformed (it is called a linear equivalent), then the canonical form is a usual one for a controllable linear system. Thus our theory of canonical forms generalizes the earlier transformation (to linear systems) results. Our canonical form is not unique, except up to solutions of certain partial differential equations we discuss. In fact, the important aspect of this paper is the constructive procedure we introduce to reach the canonical form. As is the case in many areas of mathematics, it is often easier to work with the canonical form than in arbitrary coordinate variables.

Su, R.↗
NASA NTRS1988

A Mathematical Formulation of the SCOLE Control Problem. Part 2: Optimal Compensator Design

The study initiated in Part 1 of this report is concluded and optimal feedback control (compensator) design for stability augmentation is considered, following the mathematical formulation developed in Part 1. Co-located (rate) sensors and (force and moment) actuators are assumed, and allowing for both sensor and actuator noise, stabilization is formulated as a stochastic regulator problem. Specializing the general theory developed by the author, a complete, closed form solution (believed to be new with this report) is obtained, taking advantage of the fact that the inherent structural damping is light. In particular, it is possible to solve in closed form the associated infinite-dimensional steady-state Riccati equations. The SCOLE model involves associated partial differential equations in a single space variable, but the compensator design theory developed is far more general since it is given in the abstract wave equation formulation. The results thus hold for any multibody system so long as the basic model is linear.

Balakrishnan, A. V.↗
NASA NTRS1989

Autonomous spacecraft attitude control using magnetic torquing only

Magnetic torquing of spacecraft has been an important mechanism for attitude control since the earliest satellites were launched. Typically a magnetic control system has been used for precession/nutation damping for gravity-gradient stabilized satellites, momentum dumping for systems equipped with reaction wheels, or momentum-axis pointing for spinning and momentum-biased spacecraft. Although within the small satellite community there has always been interest in expensive, light-weight, and low-power attitude control systems, completely magnetic control systems have not been used for autonomous three-axis stabilized spacecraft due to the large computational requirements involved. As increasingly more powerful microprocessors have become available, this has become less of an impediment. These facts have motivated consideration of the all-magnetic attitude control system presented here. The problem of controlling spacecraft attitude using only magnetic torquing is cast into the form of the Linear Quadratic Regulator (LQR), resulting in a linear feedback control law. Since the geomagnetic field along a satellite trajectory is not constant, the system equations are time varying. As a result, the optimal feedback gains are time-varying. Orbit geometry is exploited to treat feedback gains as a function of position rather than time, making feasible the onboard solution of the optimal control problem. In simulations performed to date, the control laws have shown themselves to be fairly robust and a good candidate for an onboard attitude control system.

Musser, Keith L.↗
NASA NTRS1990

Output model-following control synthesis for an oblique-wing aircraft

Recent interest in oblique-wing aircraft has focused on the potential aerodynamic performance advantage of a variable-skew oblique wing over a conventional or symmetric sweep wing. Unfortunately, the resulting asymmetric configuration has significant aerodynamic and inertial cross-coupling between the aircraft longitudinal and lateral-directional axes. Presented here is a decoupling control law synthesis technique that integrates stability augmentation, decoupling, and the direct incorporation of desired handling qualities into an output feedback controller. The proposed design technique uses linear quadratic regulator concepts in the framework of explicit model following. The output feedback strategy used is a suboptimal projection from the state space to the output space. Dynamics are then introduced into the controller to improve steady-state performance and increase system robustness. Closed-loop performance is shown by application of the control laws to the linearized equations of motion and nonlinear simulation of an oblique-wing aircraft.

Pahle, Joseph W.↗
NASA NTRS1992

Large-amplitude finite element flutter analysis of composite panels in hypersonic flow

Consideration is given to a finite-element approach for determining the nonlinear flutter characteristics of 3D thin laminated composite panels using the full third-order-piston transverse-loading aerodynamic theory. The unsteady hypersonic aerodynamic theory and the von Karman large-deflection-plate theory are used to formulate the aeroelasticity problem. Nonlinear flutter analyses are performed to assess the influence of the higher-order aerodynamic theory on the structure's limit-cycle amplitude and the dynamic pressure of the flow velocity. A solution procedure is presented to solve the nonlinear panel flutter and large-amplitude free-vibration finite-element equations. Nonlinear flutter analyses are performed for different boundary support-conditions and for various system parameters. Linear finite-element flutter for isotropic and composite panels and large-amplitude isotropic panel flutter results are compared with existing classical solutions. The large-amplitude panel flutter results using the full third-order piston aerodynamic theory are presented to assess the influence of the nonlinear aerodynamic theory.

Gray, Carl E., Jr.↗
NASA NTRS1990

Effects of structural offset, axial shortening, and gravitational torque on the slewing of a flexible spacecraft

The large-angle maneuvering of a Shuttle-beam-reflector spacecraft in the plane of a circular earth orbit is examined by considering the effects of the structural offset connection, the axial shortening, and the gravitational torque on the slewing motion. The offset effect is analyzed by changing the attachment point of the reflector to the beam. As the attachment point is moved away from the mass center of the reflector, the responses of the nonlinear system deviate from those of the linearized system. The axial geometric shortening effect induced by the deformation of the beam contributes to the system equations through second order terms in the modal amplitudes and rates. The gravitational torque effect is relatively small.

Li, Feiyue↗
NASA NTRS1993

Laplace's equation and the Dirichlet-Neumann map in multiply connected domains

A variety of problems in material science and fluid dynamics require the solution of Laplace's equation in multiply connected domains. Integral equation methods are natural candidates for such problems, since they discretize the boundary alone, require no special effort for free boundaries, and achieve superalgebraic convergence rates on sufficiently smooth domains in two space dimensions, regardless of shape. Current integral equation methods for the Dirichlet problem, however, require the solution of M independent problems of dimension N, where M is the number of boundary components and N is the total number of points in the discretization. In this paper, we present a new boundary integral equation approach, valid for both interior and exterior problems, which requires the solution of a single linear system of dimension N + M. We solve this system by making use of an iterative method (GMRES) combined with the last multipole method for the rapid calculation of the necessary matrix vector products. For a two-dimensional system with 200 components and 100 points on each boundary, we gain a speedup of a factor of 100 from the new analytic formulation and a factor of 50 from the fast multipole method. The resulting scheme brings large scale calculations in extremely complex domains within practical reach.

Greenbaum, A.↗
NASA NTRS2020

Simulation-Based Analysis and Prediction of Thrust Vector Servoelastic Coupling

A method of analysis and prediction of servoelastic coupling in launch vehicles is presented, surveying the discovery and subsequent resolution of a predicted servoelastic resonance phenomenon affecting the NASA Space Launch System launch vehicle at specific flight conditions. A physics-based linearized multibody mechanization of the governing equations is combined with first principles analysis to demonstrate that antisymmetric bending of the solid rocket motors leads to a reduction of equivalent viscous modal damping through coupling with the thrust vector control actuators. The sensitivity to parameters and the effects of the resonance phenomenon on flight control performance and stability are confirmed through extensive simulation verification in the time and frequency domain. A novel enhancement in model fidelity that accounts for Coriolis effects of fluid flow on bending within the solid rocket motor case and nozzle is shown to add sufficient damping to reduce the risk of adverse control-structure interaction.

Jeb S Orr↗
NASA NTRS2002

Definition of Contravariant Velocity Components

This is an old issue in computational fluid dynamics (CFD). What is the so-called contravariant velocity or contravariant velocity component? In the article, we review the basics of tensor analysis and give the contravariant velocity component a rigorous explanation. For a given coordinate system, there exist two uniquely determined sets of base vector systems - one is the covariant and another is the contravariant base vector system. The two base vector systems are reciprocal. The so-called contravariant velocity component is really the contravariant component of a velocity vector for a time-independent coordinate system, or the contravariant component of a relative velocity between fluid and coordinates, for a time-dependent coordinate system. The contravariant velocity components are not physical quantities of the velocity vector. Their magnitudes, dimensions, and associated directions are controlled by their corresponding covariant base vectors. Several 2-D (two-dimensional) linear examples and 2-D mass-conservation equation are used to illustrate the details of expressing a vector with respect to the covariant and contravariant base vector systems, respectively.

Hung, Ching-Mao↗
NASA NTRS1980

Flutter analysis of an airplane with multiple structural nonlinearities in the control system

Experience has shown that the flutter prediction process for airplanes can be greatly affected by strong concentrated nonlinearities which may be localized in the linking elements of the control mechanism, in the pivot joints of variable-sweep-wing systems, and in the connecting points between wing and pylon-mounted external stores. The principle of equivalent linearization offers an efficent possibility for solving the related nonlinear flutter equations in the frequency domain as a complement to the well-known time domain procedures. Taking as an example an airplane with nonlinear control characteristics, it is demonstrated how the equivalent linearization approach can be extended to rather complicated systems with multiple sets of strongly interacting, concentrated nonlinearities.

Breitbach, E. J.↗
NASA NTRS2017

On High-Order Upwind Methods for Advection

In the fourth installment of the celebrated series of five papers entitled "Towards the ultimate conservative difference scheme", Van Leer (1977) introduced five schemes for advection, the first three are piecewise linear, and the last two, piecewise parabolic. Among the five, scheme I, which is the least accurate, extends with relative ease to systems of equations in multiple dimensions. As a result, it became the most popular and is widely known as the MUSCL scheme (monotone upstream-centered schemes for conservation laws). Schemes III and V have the same accuracy, are the most accurate, and are closely related to current high-order methods. Scheme III uses a piecewise linear approximation that is discontinuous across cells, and can be considered as a precursor of the discontinuous Galerkin methods. Scheme V employs a piecewise quadratic approximation that is, as opposed to the case of scheme III, continuous across cells. This method is the basis for the on-going "active flux scheme" developed by Roe and collaborators. Here, schemes III and V are shown to be equivalent in the sense that they yield identical (reconstructed) solutions, provided the initial condition for scheme III is defined from that of scheme V in a manner dependent on the CFL number. This equivalence is counter intuitive since it is generally believed that piecewise linear and piecewise parabolic methods cannot produce the same solutions due to their different degrees of approximation. The finding also shows a key connection between the approaches of discontinuous and continuous polynomial approximations. In addition to the discussed equivalence, a framework using both projection and interpolation that extends schemes III and V into a single family of high-order schemes is introduced. For these high-order extensions, it is demonstrated via Fourier analysis that schemes with the same number of degrees of freedom 𝐾 per cell, in spite of the different piecewise polynomial degrees, share the same sets of eigenvalues and thus, have the same stability and accuracy. Moreover, these schemes are accurate to order 2𝐾−1, which is higher than the expected order of 𝐾.

high-order methods↗
NASA NTRS1972

Lie theory and control systems defined on spheres

It is shown that in constructing a theory for the most elementary class of control problems defined on spheres, some results from the Lie theory play a natural role. To understand controllability, optimal control, and certain properties of stochastic equations, Lie theoretic ideas are needed. The framework considered here is the most natural departure from the usual linear system/vector space problems which have dominated control systems literature. For this reason results are compared with those previously available for the finite dimensional vector space case.

Brockett, R. W.↗