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,783 records · Page 99

Uniform accuracy of implicit-explicit Runge-Kutta (IMEX-RK) schemes for hyperbolic systems with relaxation

Implicit-explicit Runge-Kutta (IMEX-RK) schemes are popular methods to treat multiscale equations that contain a stiff part and a non-stiff part, where the stiff part is characterized by a small parameter. Here, in this work, we prove rigorously the uniform stability and uniform accuracy of a class of IMEX-RK schemes for a linear hyperbolic system with stiff relaxation. The result we obtain is optimal in the sense that it holds regardless of the value of and the order of accuracy is the same as the design order of the original scheme, i.e., there is no order reduction.

97 MATHEMATICS AND COMPUTING↗

NASA-Ames three-dimensional potential flow analysis system (POTFAN) equation solver code (SOLN) version 1

A computer program known as SOLN was developed as an independent segment of the NASA-Ames three-dimensional potential flow analysis systems of linear algebraic equations. Methods used include: LU decomposition, Householder's method, a partitioning scheme, and a block successive relaxation method. Due to the independent modular nature of the program, it may be used by itself and not necessarily in conjunction with other segments of the POTFAN system.

Davis, J. E.↗

Image-gathering system design for information and fidelity

Image gathering and processing are assessed in terms of information and fidelity, and the relationship between these two figures of merit is examined. It is assumed that the system is linear and isoplanatic and that the signal and noise amplitudes are Gaussian, wide-sense stationary, and statistically independent. Within these constraints, it is found that the combined process of image gathering and reconstruction (which is intended to reproduce the output of the image-gathering system) behaves as optical, or photographic, image formation in that the informationally optimized design of the image-gathering system ordinarily does not maximize the fidelity of the reconstructed image. The combined process of image gathering and restoration (which is intended to reproduce the input of the image-gathering system) behaves more as a communication channel in that the informationally optimized design of the image-gathering system tends to maximize the fidelity of optimally restored representations of the input.

Huck, Friedrich O.↗

AESOP: An interactive computer program for the design of linear quadratic regulators and Kalman filters

AESOP is a computer program for use in designing feedback controls and state estimators for linear multivariable systems. AESOP is meant to be used in an interactive manner. Each design task that the program performs is assigned a "function" number. The user accesses these functions either (1) by inputting a list of desired function numbers or (2) by inputting a single function number. In the latter case the choice of the function will in general depend on the results obtained by the previously executed function. The most important of the AESOP functions are those that design,linear quadratic regulators and Kalman filters. The user interacts with the program when using these design functions by inputting design weighting parameters and by viewing graphic displays of designed system responses. Supporting functions are provided that obtain system transient and frequency responses, transfer functions, and covariance matrices. The program can also compute open-loop system information such as stability (eigenvalues), eigenvectors, controllability, and observability. The program is written in ANSI-66 FORTRAN for use on an IBM 3033 using TSS 370. Descriptions of all subroutines and results of two test cases are included in the appendixes.

Lehtinen, B.↗

Subsonic 3-D surface panel method for rapid analysis of multiple geometry perturbations

For any baseline aircraft configuration, a matrix of partial derivatives of surface velocity potential with respect to geometry coordinates is calculated. Linear extrapolation is then used to analyze the subcritical potential flow corresponding to a series of arbitrary small geometry perturbations. Each perturbation analysis is more than an order of magnitude more efficient than a conventional panel method solution because no influence coefficients are calculated and no large system of linear algebraic equations is solved. Wing and wing-fuselage examples are presented to demonstrate that the predicted pressure distributions are nearly exact for large changes to wing camber, thickness, and leading edge radius.

Bristow, D. R.↗

Optimal Prediction of Clocks from Finite Data

This talk is about optimal linear prediction of processes with stationary dth increments, which serve as a class of models for random clock disturbances. The predictor is obtained by orthogonal projection on the affine space of estimators whose errors are invariant to additive polynomials of degree < d. The projection conditions give a system of linear equations thatcan be solved straightforwardly for the regression coefficients. If the data are equally spaced, then the predictor can be obtained by an extension of Levinson's algorithm.

stationary increments↗

Numerical stability in problems of linear algebra.

Mathematical problems are introduced as mappings from the space of input data to that of the desired output information. Then a numerical process is defined as a prescribed recurrence of elementary operations creating the mapping of the underlying mathematical problem. The ratio of the error committed by executing the operations of the numerical process (the roundoff errors) to the error introduced by perturbations of the input data (initial error) gives rise to the concept of lambda-stability. As examples, several processes are analyzed from this point of view, including, especially, old and new processes for solving systems of linear algebraic equations with tridiagonal matrices. In particular, it is shown how such a priori information can be utilized as, for instance, a knowledge of the row sums of the matrix. Information of this type is frequently available where the system arises in connection with the numerical solution of differential equations.

Babuska, I.↗

Acoustooptic linear algebra processors - Architectures, algorithms, and applications

Architectures, algorithms, and applications for systolic processors are described with attention to the realization of parallel algorithms on various optical systolic array processors. Systolic processors for matrices with special structure and matrices of general structure, and the realization of matrix-vector, matrix-matrix, and triple-matrix products and such architectures are described. Parallel algorithms for direct and indirect solutions to systems of linear algebraic equations and their implementation on optical systolic processors are detailed with attention to the pipelining and flow of data and operations. Parallel algorithms and their optical realization for LU and QR matrix decomposition are specifically detailed. These represent the fundamental operations necessary in the implementation of least squares, eigenvalue, and SVD solutions. Specific applications (e.g., the solution of partial differential equations, adaptive noise cancellation, and optimal control) are described to typify the use of matrix processors in modern advanced signal processing.

Casasent, D.↗

Tracking and disturbance rejection of MIMO nonlinear systems with PI controller

The tracking and disturbance rejection of a class of MIMO nonlinear systems with linear proportional plus integral (PI) compensator is studied. Roughly speaking, it is shown that if the given nonlinear plant is exponentially stable and has a strictly increasing dc steady-state I/O map, then a simple PI compensator can be used to yield a stable unity-feedback closed-loop system which asymptotically tracks reference inputs that tend to constant vectors and asymptotically rejects disturbances that tend to constant vectors.

Desoer, C. A.↗

Transparent Refractory Aerogels for Efficient Spectral Control in High–Temperature Solar Power Generation

Although spectrally selective materials play a key role in existing and emerging solar thermal technologies, temperature-related degradation currently limits their use to below 700 °C in vacuum and even lower temperatures in air. Here a solar-transparent refractory aerogel that offers stable performance up to 800 °C in air is demonstrated, which is significantly greater than its silica counterpart. This improved stability is attributed to the formation of a refractory aluminum silicate phase, which is synthesized using a conformal single cycle of atomic layer deposition within the high-aspect-ratio pores of silica aerogels. Furthermore, based on direct heat loss measurements, the transparent refractory aerogel achieves a receiver efficiency of 75% at 100 suns and an absorber temperature of 700 °C, which is a 5% improvement over the state of the art. Transparent refractory aerogels may find widespread applicability in solar thermal technologies by enabling the use of lower-cost optical focusing systems and eliminating the need for highly evacuated receivers. In particular, a shift to higher operating temperatures while maintaining a high receiver efficiency can enable the use of advanced supercritical CO 2 power cycles and ultimately translate to an ≈10% (absolute) improvement in solar-to-electrical conversion efficiency relative to existing linear concentrating systems.

14 SOLAR ENERGY↗

A first-order Green's function approach to supersonic oscillatory flow: A mixed analytic and numeric treatment

A frequency domain Green's Function Method for unsteady supersonic potential flow around complex aircraft configurations is presented. The focus is on the supersonic range wherein the linear potential flow assumption is valid. In this range the effects of the nonlinear terms in the unsteady supersonic compressible velocity potential equation are negligible and therefore these terms will be omitted. The Green's function method is employed in order to convert the potential flow differential equation into an integral one. This integral equation is then discretized, through standard finite element technique, to yield a linear algebraic system of equations relating the unknown potential to its prescribed co-normalwash (boundary condition) on the surface of the aircraft. The arbitrary complex aircraft configuration (e.g., finite-thickness wing, wing-body-tail) is discretized into hyperboloidal (twisted quadrilateral) panels. The potential and co-normalwash are assumed to vary linearly within each panel. The long range goal is to develop a comprehensive theory for unsteady supersonic potential aerodynamic which is capable of yielding accurate results even in the low supersonic (i.e., high transonic) range.

Freedman, M. I.↗

Elimination on sparse symmetric systems of a special structure.

Consideration of the problem of finding a permutation of rows and columns and an algorithm for solving ordered systems of linear algebraic equations with sparse matrices having a certain regular structure. Two approaches to the solution of this problem, in which the sparsity is used to some extent, are outlined. One of them is a very general approach where optimal (or nearly optimal) ordering is sought and the algorithm for solving the ordered system treats the matrix element by element to perform only necessary operations. The other approach involves the use of band matrices. After comparing these two approaches, a third approach is then suggested which involves the use of pipe matrices, and a means of ordering the rows and columns to obtain this type of matrix is presented. Examples of matrices reordered by the proposed procedure are cited.

Segethova, J.↗

Computer aided design of digital controller for radial active magnetic bearings

A five degree of freedom Active Magnetic Bearing (AMB) system is developed which is controlled by digital controllers. The model of the radial AMB system is linearized and the state equation is derived. Based on the state variables feedback theory, digital controllers are designed. The performance of the controllers are evaluated according to experimental results. The Computer Aided Design (CAD) method is used to design controllers for magnetic bearings. The controllers are implemented with a digital signal processing (DSP) system. The control algorithms are realized with real-time programs. It is very easy to change the controller by changing or modifying the programs. In order to identify the dynamic parameters of the controlled magnetic system, a special experiment was carried out. Also, the online Recursive Least Squares (RLS) parameter identification method is studied. It can be realized with the digital controllers. Online parameter identification is essential for the realization of an adaptive controller.

Cai, Zhong↗

An Analytical Calibration Approach for the Polarimetric Airborne C Band Radiometer

Passive microwave remote sensing is sensitive to the quantity and distribution of water in soil and vegetation. During summer 2000, the Microwave Geophysics Group a t the University of Michigan conducted the seventh Radiobrighness Energy Balance Experiment (REBEX-7) over a corn canopy in Michigan. Long time series of brightness temperatures, soil moisture and micrometeorology on the plot were taken. This paper addresses the calibration of the NASA GSFC polarimetric airborne C band microwave radiometer (ACMR) that participated in REBEX-7. These passive polarimeters are typically calibrated using an end-to-end approach based upon a standard artificial target or a well-known geophysical target. Analyzing the major internal functional subsystems offers a different perspective. The primary goal of this approach is to provide a transfer function that not only describes the system in its entire5 but also accounts for the contributions of each subsystem toward the final modified Stokes parameters. This approach does not assume that the radiometric system is linear as it does not take polarization isolation for granted, and it also serves as a realistic instrument simulator, a useful tool for future designs. The ACMR architecture can be partitioned into functional subsystems. The characteristics of each subsystem was extensively measured and the estimated parameters were imported into the overall dosed form system model. Inversion of the model yields a calibration for the modeled Stokes parameters with uncertainties of 0.2 K for the V and H polarizations and 2.4 K for the 3rd and 4th parameters. Application to the full Stokes parameters over a senescent cornfield is presented.

Pham, Hanh↗

Quantum algorithm for the linear Vlasov equation with collisions

The Vlasov equation is a nonlinear partial differential equation that provides a first-principles description of the dynamics of plasmas. Its linear limit is routinely used in plasma physics to investigate plasma oscillations and stability. In this paper, we present a quantum algorithm that simulates the linearized Vlasov equation with and without collisions, in the one-dimensional electrostatic limit. Rather than solving this equation in its native spatial and velocity phase space, we adopt an efficient representation in the dual space yielded by a Fourier-Hermite expansion. For a given simulation time, the Fourier-Hermite representation is exponentially more compact, thus yielding a classical algorithm that can match the performance of a previously proposed quantum algorithm for this problem. Further, this representation results in a system of linear ordinary differential equations (ODEs) which can be solved with well-developed quantum algorithms: a Hamiltonian simulation in the collisionless case, and quantum ODE solvers in the collisional case. In particular, we demonstrate that a quadratic speedup in system size is attainable.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reciprocity principle in duct acoustics

Various reciprocity relations in duct acoustics have been derived on the basis of the spatial reciprocity principle implied in Green's functions for linear waves. The derivation includes the reciprocity relations between mode conversion coefficients for reflection and transmission in nonuniform ducts, and the relation between the radiation of a mode from an arbitrarily terminated duct and the absorption of an externally incident plane wave by the duct. Such relations are well defined as long as the systems remain linear, regardless of acoustic properties of duct nonuniformities which cause the mode conversions.

Cho, Y. C.↗

Reciprocity principle in duct acoustics

Various reciprocity relations in duct acoustics have been derived on the basis of the spatial reciprocity principle implied in Green's functions for linear waves. The derivation includes the reciprocity relations between mode conversion coefficients for reflection and transmission in nonuniform ducts, and the relation between the radiation of a mode from an arbitrarily terminated duct and the absorption of an externally incident plane wave by the duct. Such relations are well defined as long as the systems remain linear, regardless of acoustic properties of duct nonuniformities which cause the mode conversions.

Cho, Y.-C.↗