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 415 records · Page 23

Discrete integration of continuous Kalman filtering equations for time invariant second-order structural systems

A general form for the first-order representation of the continuous second-order linear structural-dynamics equations is introduced to derive a corresponding form of first-order continuous Kalman filtering equations. Time integration of the resulting equations is carried out via a set of linear multistep integration formulas. It is shown that a judicious combined selection of computational paths and the undetermined matrices introduced in the general form of the first-order linear structural systems leads to a class of second-order discrete Kalman filtering equations involving only symmetric sparse N x N solution matrices.

Park, K. C.↗

Time domain desensitized specific optimal system design

A digital computer algorithm for desensitized specific optimal system design for general time-invariant nonlinear systems is presented. The algorithm utilizes a technique of partitioning and uncoupling of the original specific closed-loop system state and sensitivity equations into linear, nonlinear, and control subsystems. Results of the application of the algorithm in the design of a Saturn V Launch Vehicle analog attitude control system and an LST image motion compensation digital controller are given.

Henson, T. F.↗

A Method for Longitudinal and Lateral Range Control for a High-Drag Low-Lift Vehicle Entering the Atmosphere of a Rotating Earth

A study has been made of a method for controlling the trajectory of a high-drag low-lift entry vehicle to a desired longitude and latitude on the surface of a rotating earth. By use of this control technique the vehicle can be guided to the desired point when the present position and heading of the vehicle are known and the desired longitude and latitude are specified. The present study makes use of a single reference trajectory and an estimate of the lift and side-force capabilities of the vehicle. This information is stored in a control-logic system and used with linear control equations to guide the vehicle to the desired destination. Results are presented of a number of trajectory studies which describe the operation of the control system and illustrate its ability to control the vehicle trajectory to the desired landing area.

Young, John W.↗

Some estimation formulae for continuous time-invariant linear systems

In this brief paper we examine a Riccati equation decomposition due to Reid and Lainiotis and apply the result to the continuous time-invariant linear filtering problem. Exploitation of the time-invariant structure leads to integration-free covariance recursions which are of use in covariance analyses and in filter implementations. A super-linearly convergent iterative solution to the algebraic Riccati equation (ARE) is developed. The resulting algorithm, arranged in a square-root form, is thought to be numerically stable and competitive with other ARE solution methods. Certain covariance relations that are relevant to the fixed-point and fixed-lag smoothing problems are also discussed.

Bierman, G. J.↗

A Scalable Multigrid Reduction Framework for Multiphase Poromechanics of Heterogeneous Media

Simulation of multiphase poromechanics involves solving a multiphysics problem in which multiphase flow and transport are tightly coupled with the porous medium deformation. To capture this dynamic interplay, fully implicit methods, also known as monolithic approaches, are usually preferred. The main bottleneck of a monolithic approach is that it requires solution of large linear systems that result from the discretization and linearization of the governing balance equations. Because such systems are nonsymmetric, indefinite, and highly ill-conditioned, preconditioning is critical for fast convergence. Recently, most efforts in designing efficient preconditioners for multiphase poromechanics have been dominated by physics-based strategies. Current state-of-the-art “black-box” solvers such as algebraic multigrid (AMG) are ineffective because they cannot effectively capture the strong coupling between the mechanics and the flow subproblems, as well as the coupling inherent in the multiphase flow and transport process. In this work, we develop an algebraic framework based on multigrid reduction (MGR) that is suited for tightly coupled systems of PDEs. Using this framework, the decoupling between the equations is done algebraically through defining appropriate interpolation and restriction operators. One can then employ existing solvers for each of the decoupled blocks or design a new solver based on knowledge of the physics. We demonstrate the applicability of our framework when used as a “black-box” solver for multiphase poromechanics. Here, we show that the framework is flexible to accommodate a wide range of scenarios, as well as efficient and scalable for large problems.

97 MATHEMATICS AND COMPUTING↗

An Optimized Multicolor Point-Implicit Solver for Unstructured Grid Applications on Graphics Processing Units

In the field of computational fluid dynamics, the Navier-Stokes equations are often solved using an unstructuredgrid approach to accommodate geometric complexity. Implicit solution methodologies for such spatial discretizations generally require frequent solution of large tightly-coupled systems of block-sparse linear equations. The multicolor point-implicit solver used in the current work typically requires a significant fraction of the overall application run time. In this work, an efficient implementation of the solver for graphics processing units is proposed. Several factors present unique challenges to achieving an efficient implementation in this environment. These include the variable amount of parallelism available in different kernel calls, indirect memory access patterns, low arithmetic intensity, and the requirement to support variable block sizes. In this work, the solver is reformulated to use standard sparse and dense Basic Linear Algebra Subprograms (BLAS) functions. However, numerical experiments show that the performance of the BLAS functions available in existing CUDA libraries is suboptimal for matrices representative of those encountered in actual simulations. Instead, optimized versions of these functions are developed. Depending on block size, the new implementations show performance gains of up to 7x over the existing CUDA library functions.

Zubair, Mohammad↗

The transonic aerofoil problem with embedded shocks.

The integral equation approach to the mixed flow problem of infinite wings at high subsonic speeds is adopted for non-circulatory and circulatory (lifting) flows. The solutions are determined from a system of non-linear algebraic equations and, to ensure always unique solutions, the method of differentiation with respect to a parameter has been applied. The resulting Cauchy problem is then solved with the linearised flow solution as the initial value vector. For the case of embedded shocks in the flow field, the method of steepest descent has been added to the calculation scheme. Results for subcritical and supercritical flows past aerofoils are given and compared with solutions obtained by finite-difference techniques.

Norstrud, H.↗

A high-order explicit Runge-Kutta approximation technique for the shallow water equations

Here, we introduce a high-order space–time approximation of the Shallow Water Equations with sources that is invariant-domain preserving (IDP), well-balanced with respect to rest states, and employs a novel explicit Runge–Kutta (ERK) introduced in Ern and Guermond (SIAM J. Sci. Comput. 44(5), A3366–A3392, 2022) for systems of non-linear conservation equations. The resulting method is then numerically illustrated through verification and validation.

97 MATHEMATICS AND COMPUTING↗

Influence of wind turbine design parameters on linearized physics-based models in OpenFAST

Abstract. While most physics involved in wind energy are nonlinear, linearization of the underlying nonlinear wind system equations is often important for understanding the system response and exploiting well-established methods and tools for analyzing linear systems. Linearized models are important for eigenanalysis (to derive structural natural frequencies, damping ratios, and mode shapes), controls design (based on linear state-space models), etc. In controls co-design, wherein methods often rely on linearized time-domain models of the physics, the physical structure (often called the plant) and controller are designed and optimized concurrently, so it is important to understand how changes to the physical design affect the linearized system. This work summarizes efforts done to understand the impact of design parameter variations in the physical system (e.g., mass, stiffness, geometry, and aerodynamic and hydrodynamic coefficients) on the linearized system using OpenFAST.

17 WIND ENERGY↗

ORACLS: A system for linear-quadratic-Gaussian control law design

A modern control theory design package (ORACLS) for constructing controllers and optimal filters for systems modeled by linear time-invariant differential or difference equations is described. Numerical linear-algebra procedures are used to implement the linear-quadratic-Gaussian (LQG) methodology of modern control theory. Algorithms are included for computing eigensystems of real matrices, the relative stability of a matrix, factored forms for nonnegative definite matrices, the solutions and least squares approximations to the solutions of certain linear matrix algebraic equations, the controllability properties of a linear time-invariant system, and the steady state covariance matrix of an open-loop stable system forced by white noise. Subroutines are provided for solving both the continuous and discrete optimal linear regulator problems with noise free measurements and the sampled-data optimal linear regulator problem. For measurement noise, duality theory and the optimal regulator algorithms are used to solve the continuous and discrete Kalman-Bucy filter problems. Subroutines are also included which give control laws causing the output of a system to track the output of a prescribed model.

Armstrong, E. S.↗

A variable-step-size robust delta modulator.

Description of an analytically obtained optimum adaptive delta modulator-demodulator configuration. The device utilizes two past samples to obtain a step size which minimizes the mean square error for a Markov-Gaussian source. The optimum system is compared, using computer simulations, with a linear delta modulator and an enhanced Abate delta modulator. In addition, the performance is compared to the rate distortion bound for a Markov source. It is shown that the optimum delta modulator is neither quantization nor slope-overload limited. The highly nonlinear equations obtained for the optimum transmitter and receiver are approximated by piecewise-linear equations in order to obtain system equations which can be transformed into hardware. The derivation of the experimental system is presented.

Song, C. L.↗

Identification of general linear mechanical systems

Previous work in identification theory has been concerned with the general first order time derivative form. Linear mechanical systems, a large and important class, naturally have a second order form. This paper utilizes this additional structural information for the purpose of identification. A realization is obtained from input-output data, and then knowledge of the system input, output, and inertia matrices is used to determine a set of linear equations whereby we identify the remaining unknown system matrices. Necessary and sufficient conditions on the number, type and placement of sensors and actuators are given which guarantee identificability, and less stringent conditions are given which guarantee generic identifiability. Both a priori identifiability and a posteriori identifiability are considered, i.e., identifiability being insured prior to obtaining data, and identifiability being assured with a given data set.

Sirlin, S. W.↗

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↗

Linearized frequency domain Landau-Lifshitz-Gilbert equation formulation

We present a general finite element linearized Landau-Lifshitz-Gilbert equation (LLGE) solver for magnetic systems under weak time-harmonic excitation field. The linearized LLGE is obtained by assuming a small deviation around the equilibrium state of the magnetic system. Inserting such expansion into LLGE and keeping only first order terms gives the linearized LLGE, which gives a frequency domain solution for the complex magnetization amplitudes under an external time-harmonic applied field of a given frequency. We solve the linear system with an iterative solver using generalized minimal residual method. We construct a preconditioner matrix to effectively solve the linear system. The validity, effectiveness, speed, and scalability of the linear solver are demonstrated via numerical examples.

36 MATERIALS SCIENCE↗

Attitude stability of a spinning spacecraft during appendage deployment/retraction

The work presented is motivated by the need for a national satellite rescue policy, not the ad hoc policy now in place. In studying different approaches for a national policy, the issue of capture and stabilization of a tumbling spacecraft must be addressed. For a rescue mission involving a tumbling spacecraft, it may be advantageous to have a rescue vehicle which is compact and 'rigid' during the rendezvous/capture phase. After capture, passive stabilization techniques could be utilized as an efficient means of detumbling the resulting system (i.e., both the rescue vehicle and captures spacecraft). Since the rescue vehicle is initially compact and 'rigid,' significant passive stabilization through energy dissipation can only be achieved through the deployment of flexible appendages. Once stabilization is accomplished, retraction of the appendages before maneuvering the system to its final destination may also prove advantageous. It is therefore of paramount interest that we study the effect of appendage deployment/retraction on the attitude stability of a spacecraft. Particular interest should be paid to appendage retraction, since if this process is destabilizing, passive stabilization as proposed may not be useful. Over the past three decades, it has been an 'on-again-off-again affair' with the problem of spacecraft appendage deployment. In most instances, these studies have been numerical simulations of specific spacecraft configurations for which there were specific concerns. The primary focus of these studies was the behavior of the appendage during deployment; the effects of appendage retraction was considered only in one of these studies. What is missing in the literature is a thorough study of the effects of appendage deployment/retraction on the attitude stability of a spacecraft. This paper presents a rigorous analysis of the stability of a spinning spacecraft during the deployment or the retraction of an appendage. The analysis is simplified such that meaningful insights into the problem can be inferred; it is not overly simplified such that critical dynamical behavior is neglected. The system is analyzed assuming that the spacecraft hub is rigid. The appendage deployment mechanism is modeled as a point mass on a massless rod whose length undergoes prescribed changes. Simplified flexibility effects of the appendage are included. The system is examined for stability by linearizing the equations in terms of small deviations from steady, noninterfering coning motion. Routh's procedure for analyzing small deviations from steady motion in dynamical systems is utilized in the analysis. The system of equations are nondimensionalized to facilitate parametric studies. The results are presented in terms of a reduced number of nondimensional parameters so that some general conclusions may be drawn. Verification of the linear analysis is presented through numerical simulations of the complete nonlinear, nonautonomous, coupled equations.

Fitz-Coy, Norman↗

Error magnitude estimation in model-reference adaptive systems

A second order approximation is derived from a linearized error characteristic equation for Lyapunov designed model-reference adaptive systems and is used to estimate the maximum error between the model and plant states, and the time to reach this peak following a plant perturbation. The results are applicable in the analysis of plants containing magnitude-dependent nonlinearities.

Colburn, B. K.↗