Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear 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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Ocular torsion on earth and in weightlessness

Otolith function is studied by means of measurements of ocular torsion under various acceleration environments on earth and in weightlessness. Photographic measurements of ocular torsion as indicated by the rotation of landmarks on the iris with respect to head-fixed fiducial marks were obtained in subjects undergoing horizontal linear acceleration in a ground-based version of the space sled, lateral acceleration from weightlessness during pullout from the free-fall portion of parabolic flight, and optokinetic stimulation about the roll axis in the supine position in the laboratory and during weightlessness. The responses of ocular torsion to horizontal acceleration are in agreement with a simple low-order linear system with a dominant time constant of 0.33 sec, with a transfer function fit by a model with a pure delay of 0 to 400 msec and a first-order lag. In the pullout experiment, torsion was not observed in response to the onset of acceleration in the right-ear-down position, although it was present in response to the lateral stimulus. Results of the roll vection experiments indicate the independence of ocular torsion and visually induced tilt. In addition, an automatic video system using a soft contact lens target is presented which has been developed for ocular torsion measurements.

Young, L. R.↗
NASA NTRS1980

Asymptotic eigenstructures

The behavior of the closed loop eigenstructure of a linear system with output feedback is analyzed as a single parameter multiplying the feedback gain is varied. An algorithm is presented that computes the asymptotically infinite eigenstructure, and it is shown how a system with high gain feedback decouples into single input single ouput systems. Then a synthesis algorithm is presented which uses full state feedback to achieve a desired asymptotic eigenstructure.

Thompson, P. M.↗
NASA NTRS1980

Hamiltonian indices and rational spectral densities

Several (global) topological properties of various spaces of linear systems, particularly symmetric, lossless, and Hamiltonian systems, and multivariable spectral densities of fixed McMillan degree are announced. The study is motivated by a result asserting that on a connected but not simply connected manifold, it is not possible to find a vector field having a sink as its only critical point. In the scalar case, this is illustrated by showing that only on the space of McMillan degree = /Cauchy index/ = n, scalar transfer functions can one define a globally convergent vector field. This result holds both in discrete-time and for the nonautonomous case. With these motivations in mind, theorems of Bochner and Fogarty are used in showing that spaces of transfer functions defined by symmetry conditions are, in fact, smooth algebraic manifolds.

Byrnes, C. I.↗
NASA NTRS1983

Calculations of transonic potential flows by a parameter free procedure

A parameter-free computational procedure for transonic potential flows around an airfoil is presented. The method is based on an inexact, Newton-like iterative scheme, and a preconditioned minimal residual algorithm is used to solve the linear system of equations arising at each iterate. This iterative procedure is formulated so that the difficulties associated with the approximate factorization technique, namely the requirement of iteration parameters and the boundary conditions for the intermediate variable are eliminated. Numerical experiments are given for potential flows around NACA 0012 airfoil with different Mach numbers and different angles of attack. These results are then compared with those obtained by the approximate factorization techniques.

Wong, Y. S.↗
NASA NTRS1983

Acceleration of convergence of vector sequences

A general approach to the construction of convergence acceleration methods for vector sequence is proposed. Using this approach, one can generate some known methods, such as the minimal polynomial extrapolation, the reduced rank extrapolation, and the topological epsilon algorithm, and also some new ones. Some of the new methods are easier to implement than the known methods and are observed to have similar numerical properties. The convergence analysis of these new methods is carried out, and it is shown that they are especially suitable for accelerating the convergence of vector sequences that are obtained when one solves linear systems of equations iteratively. A stability analysis is also given, and numerical examples are provided. The convergence and stability properties of the topological epsilon algorithm are likewise given.

Sidi, A.↗
NASA NTRS1984

Preliminary results in implementing a model of the world economy on the CYBER 205: A case of large sparse nonsymmetric linear equations

A brief description of the Model of the World Economy implemented at the Institute for Economic Analysis is presented, together with our experience in converting the software to vector code. For each time period, the model is reduced to a linear system of over 2000 variables. The matrix of coefficients has a bordered block diagonal structure, and we show how some of the matrix operations can be carried out on all diagonal blocks at once.

Szyld, D. B.↗
NASA NTRS1984

Two-dimensional computer simulation of EMVJ and grating solar cells under AMO illumination

A computer program, SCAP2D (Solar Cell Analysis Program in 2-Dimensions), is used to evaluate the Etched Multiple Vertical Junction (EMVJ) and grating solar cells. The aim is to demonstrate how SCAP2D can be used to evaluate cell designs. The cell designs studied are by no means optimal designs. The SCAP2D program solves the three coupled, nonlinear partial differential equations, Poisson's Equation and the hole and electron continuity equations, simultaneously in two-dimensions using finite differences to discretize the equations and Newton's Method to linearize them. The variables solved for are the electrostatic potential and the hole and electron concentrations. Each linear system of equations is solved directly by Gaussian Elimination. Convergence of the Newton Iteration is assumed when the largest correction to the electrostatic potential or hole or electron quasi-potential is less than some predetermined error. A typical problem involves 2000 nodes with a Jacobi matrix of order 6000 and a bandwidth of 243.

Gray, J. L.↗
NASA NTRS1984

Investigation and appreciation of optimal output feedback. Volume 1: A convergent algorithm for the stochastic infinite-time discrete optimal output feedback problem

The stochastic, infinite time, discrete output feedback problem for time invariant linear systems is examined. Two sets of sufficient conditions for the existence of a stable, globally optimal solution are presented. An expression for the total change in the cost function due to a change in the feedback gain is obtained. This expression is used to show that a sequence of gains can be obtained by an algorithm, so that the corresponding cost sequence is monotonically decreasing and the corresponding sequence of the cost gradient converges to zero. The algorithm is guaranteed to obtain a critical point of the cost function. The computational steps necessary to implement the algorithm on a computer are presented. The results are applied to a digital outer loop flight control problem. The numerical results for this 13th order problem indicate a rate of convergence considerably faster than two other algorithms used for comparison.

Halyo, N.↗
NASA NTRS1983

Direct and implicit optical matrix-vector algorithms

New direct and implicit algorithms for optical matrix-vector and systolic array processors are considered. Direct rather than indirect algorithms to solve linear systems and implicit rather than explicit solutions to solve second-order partial differential equations are discussed. In many cases, such approaches more properly utilize the advantageous features of optical systolic array processors. The matrix-decomposition operation (rather than solution of the simplified matrix-vector equation that results) is recognized as the computationally burdensome aspect of such problems that should be computed on an optical system. The Householder QR matrix-decomposition algorithm is considered as a specific example of a direct solution. Extensions to eigenvalue computation and formation of matrices of special structure are also noted.

Casasent, D.↗
NASA NTRS1984

Differential rotation in a solar-driven quasi-axisymmetric circulation

The concept of a quasi-axisymmetric circulation is used to explore the global scale dynamics of planetary atmospheres. A numerical circulation model applicable to Jupiter is presented, and an analytical study is performed elucidating the conditions leading to differential rotation in an atmosphere which is convectively unstable. A linear system forced by solar differential heating is considered, with nonlinear effects arising from advection being represented in the form of eddy diffusion. An empirical, latitudinal spectrum of the observed zonal wind field on Jupiter is discussed. Numerical solutions are presented which reveal banded wind fields with alternating and equatorial zonal jets and a multicellular Ferrel-Thomson meridional circulation consistent with the observed cloud striations on Jupiter. The vertical derivatives are parameterized to construct a simplified one-layer model.

Mayr, H. G.↗
NASA NTRS1984

Implicit Numerical Methods in Meteorology

The development of a fully implicit finite-difference model, whose time step is chosen solely to resolve accurately the physical flow of interest is discussed. The method is based on an operator factorization which reduces the dimensionality of the implicit approach: at each time step only (spatially) one-dimensional block-tridiagonal linear systems must be solved. The scheme uses two time levels and is second-order accurate in time. Compact implicit spatial differences are used, yielding fourth-order accuracy both vertically and horizontally. In addition, the development of a fully interactive computer code is discussed. With this code the user will have a choice of models, with various levels of accuracy and sophistication, which are imbedded, as subsets of the fully implicit 3D code.

Augenbaum, J.↗
NASA NTRS1984

The effects of time delay in man-machine control systems: Implications for design of flight simulator Visual-Display-Delay compensation

When human operators are performing precision tracking tasks, their dynamic response can often be modeled by quasilinear describing functions. That fact permits analysis of the effects of delay in certain man machine control systems using linear control system analysis techniques. The analysis indicates that a reduction in system stability is the immediate effect of additional control system delay, and that system characteristics moderate or exaggerate the importance of the delay. A selection of data (simulator and flight test) consistent with the analysis is reviewed. Flight simulator visual-display delay compensation, designed to restore pilot aircraft system stability, was evaluated in several studies which are reviewed here. The studies range from single-axis, tracking-task experiments (with sufficient subjects and trials to establish the statistical significance of the results) to a brief evaluation of compensation of a computer generated imagery (CGI) visual display system in a full six degree of freedom simulation. The compensation was effective, improvements in pilot performance and workload or aircraft handling qualities rating (HQR) were observed. Results from recent aircraft handling qualities research literature, which support the compensation design approach, are also reviewed.

Crane, D. F.↗
NASA NTRS1983

Alternatives for jet engine control

The technical progress of researches on alternatives for jet engine control, is reported. The principal new activities involved the initial testing of an input design method for choosing the inputs to a non-linear system to aid the approximation of its tensor parameters, and the beginning of order reduction studies designed to remove unnecessary monomials from tensor models.

Sain, M. K.↗
NASA NTRS1985

A Factored Implicit Scheme for Numerical Weather Prediction with Small Factorization Error

Numerical results show that, for large time steps, the factorization error can be significant, even for the slowly propagating Rossby modes. A new scheme is formulated based on a more accurate factorization of the equations. By grouping separately the terms of the equations which give rise to the fast and slow motion, the equations are factored more accurately. The fast-slow factorization eliminated the factorization error. If each of the fast and slow factors are factored again according to spatial components, the resulting scheme only involves the solution of one dimensional linear systems, and computational efficient. It is shown that the factorization error for the slow made component is negligible for this new scheme.

Augenbaum, J. M.↗
NASA NTRS1985

A fully implicit scheme for the barotropic primitive equations

An efficient implicit finite-difference method is developed and tested for a global barotropic model. The scheme requires, at each time step, the solution of only one-dimensional block-tridiagonal linear systems. This additional computation is offset by the use of a time step chosen independently of the mesh spacing. The method is second-order accurate in time and fourth-order accurate in space. Present experience indicates that this implicit method is practical for numerical simulation on fine meshes.

Cohn, S. E.↗
NASA NTRS1985

Disturbance-accommodating tracking maneuvers of flexible spacecraft

In this paper the problem of maneuvering a flexible spacecraft through a large angle is considered, where the disturbance-accommodating feedback control tracks a desired output state. The desired output state is provided from an open-loop solution for the linear system model. The components of the disturbance vector are assumed to be represented in terms of Fourier series. Closed-form solutions are provided for the Ricati, prefilter, state trajectory, and residual state trajectory equations which define the optimal control. Example maneuvers are presented where control-rate penalties have been included in the performance index for frequency-shaping, in order to smooth both the open- and closed-loop control commands.

Turner, J. D.↗
NASA NTRS1985

GMRES acceleration of computational fluid dynamics codes

The generalized minimal residual algorithm (GMRES) is a conjugate-gradient like method that applies directly to nonsymmetric linear systems of equations. In this paper, GMRES is modified to handle nonlinear equations characteristic of computational fluid dynamics. Attention is devoted to the concept of preconditioning and the role it plays in assuring rapid convergence. A formulation is developed that allows GMRES to be preconditioned by the solution procedures already built into existing computer codes. Examples are provided that demonstrate the ability of GMRES to greatly improve the robustness and rate of convergence of current state-of-the-art fluid dynamics codes. Theoretical aspects of GMRES are presented that explain why it works. Finally, the advantage GMRES enjoys over related methods such as conjugate gradients are discussed.

Wigton, L. B.↗
NASA NTRS1985

The Control of Float Zone Interfaces by the Use of Selected Boundary Conditions

Since the float zone process involves two solid-melt interfaces, possible gas interfaces, heat and mass transfers, various driving forces and complex heating sources, an analysis of the entire process would be very complex. For an initial investigation, a more feasible approach is to examine each component of the process separately. The three principal components are: (1) the shapes of the melt and solid-melt interfaces, (2) the heat and mass transfers, and (3) the heating and cooling sources. This study combined facets of all three components. The purpose of this effort was to study and compute the surface boundary conditions required to give flat (float zone) FZ solid-melt interfaces. This study was undertaken in two phases. The first phase was to investigate the solid zones surface boundary conditions required for flat solid-melt interfaces when given the melt zone surface boundary conditions. The second phase complemented the first and was to investigate the melt zone surface boundary conditions required for flat solid-melt interfaces if given the solid zones surface boundary conditions. Dual integral transform methods were used in both phases; in addition, the use of various numerical methods for differential equations and linear systems of equations were required.

Foster, L. M.↗