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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,135 records · Page 63

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Model-following control for an oblique-wing aircraft

A variable-skew oblique wing offers a substantial aerodynamic performance advantage for aircraft missions that require both high efficiency in subsonic flight and supersonic dash or cruise. The most obvious characteristic of the oblique-wing concept is the asymmetry associated with wing-skew angle which results in significant aerodynamic and inertial cross-coupling between the aircraft longitudinal and lateral-directional axes. This paper presents a technique for synthesizing a decoupling controller while providing the desired stability augmentation. The proposed synthesis procedure uses the concept of explicit model following. Linear quadratic optimization techniques are used to design the linear feedback system. The effectiveness of the control laws developed in achieving the desired decoupling is illustrated for a given flight condition by application to linearized equations of motion, and also to the nonlinear equations of six degrees of freedom of motion with nonlinear aerodynamic data.

Alag, G. S.↗

NASA NTRS1989

The theory of magnetohydrodynamic wave generation by localized sources. II - Collisionless dissipation of wave packets

The dispersion equation of Barnes (1966) is used to study the dissipation of asymptotic wave packets generated by localized periodic sources. The solutions of the equation are linear waves, damped by Landau and transit-time processes, in a collisionless warm plasma. For the case of an ideal MHD system, most of the waves emitted from a source are shown to cancel asympotically through destructive interference. The modes transporting significant flux to asymptotic distances are found to be Alfven waves and fast waves with theta (the angle between the magnetic field and the characteristics of the far-field waves) of about 0 and about pi/2.

Collins, William↗

NASA NTRS1991

A controls engineering approach for analyzing airplane input-output characteristics

An engineering approach for analyzing airplane control and output characteristics is presented. State-space matrix equations describing the linear perturbation dynamics are transformed from physical coordinates into scaled coordinates. The scaling is accomplished by applying various transformations to the system to employ prior engineering knowledge of the airplane physics. Two different analysis techniques are then explained. Modal analysis techniques calculate the influence of each system input on each fundamental mode of motion and the distribution of each mode among the system outputs. The optimal steady state response technique computes the blending of steady state control inputs that optimize the steady state response of selected system outputs. Analysis of an example airplane model is presented to demonstrate the described engineering approach.

Arbuckle, P. Douglas↗

NASA NTRS1991

Nonlinear dynamics near the stability margin in rotating pipe flow

The nonlinear evolution of marginally unstable wave packets in rotating pipe flow is studied. These flows depend on two control parameters, which may be taken to be the axial Reynolds number R and a Rossby number, q. Marginal stability is realized on a curve in the (R, q)-plane, and the entire marginal stability boundary is explored. As the flow passes through any point on the marginal stability curve, it undergoes a supercritical Hopf bifurcation and the steady base flow is replaced by a traveling wave. The envelope of the wave system is governed by a complex Ginzburg-Landau equation. The Ginzburg-Landau equation admits Stokes waves, which correspond to standing modulations of the linear traveling wavetrain, as well as traveling wave modulations of the linear wavetrain. Bands of wavenumbers are identified in which the nonlinear modulated waves are subject to a sideband instability.

Yang, Z.↗

NASA NTRS1992

Adaptive attitude control and momentum management for large-angle spacecraft maneuvers

The fully coupled equations of motion are systematically linearized around an equilibrium point of a gravity gradient stabilized spacecraft, controlled by momentum exchange devices. These equations are then used for attitude control system design of an early Space Station Freedom flight configuration, demonstrating the errors caused by the improper approximation of the spacecraft dynamics. A full state feedback controller, incorporating gain-scheduled adaptation of the attitude gains, is developed for use during spacecraft on-orbit assembly or operations characterized by significant mass properties variations. The feasibility of the gain adaptation is demonstrated via a Space Station Freedom assembly sequence case study. The attitude controller stability robustness and transient performance during gain adaptation appear satisfactory.

Parlos, Alexander G.↗

NASA NTRS1995

The Duffin-Kemmer-Petiau oscillator

In view of current interest in relativistic spin-one systems and the recent work on the Dirac Oscillator, we introduce the Duffin-Kemmer-Petiau (DKP) equation obtained by using an external potential linear in r. Since, in the non-relativistic limit, the spin 1 representation leads to a harmonic oscillator with a spin-orbit coupling of the Thomas form, we call the equation the DKP oscillator. This oscillator is a relativistic generalization of the quantum harmonic oscillator for scalar and vector bosons. We show that it conserves total angular momentum and that it is exactly solvable. We calculate and discuss the eigenspectrum of the DKP oscillator in the spin 1 representation.

Nedjadi, Youcef↗

NASA NTRS1999

1-D Numerical Analysis of ABCC Engine Performance

ABCC engine combines air breathing and rocket engine into a single engine to increase the specific impulse over an entire flight trajectory. Except for the heat source, the basic operation of the ABCC is similar to the basic operation of the RBCC engine. The ABCC is intended to have a higher specific impulse than the RBCC for single stage Earth to orbit vehicle. Computational fluid dynamics (CFD) is a useful tool for the analysis of complex transport processes in various components in ABCC propulsion system. The objective of the present research was to develop a transient 1-D numerical model using conservation of mass, linear momentum, and energy equations that could be used to predict flow behavior throughout a generic ABCC engine following a flight path. At specific points during the development of the 1-D numerical model a myriad of tests were performed to prove the program produced consistent, realistic numbers that follow compressible flow theory for various inlet conditions.

Holden, Richard↗

NASA NTRS2004

A-Posteriori Error Estimation for Hyperbolic Conservation Laws with Constraint

This lecture considers a-posteriori error estimates for the numerical solution of conservation laws with time invariant constraints such as those arising in magnetohydrodynamics (MHD) and gravitational physics. Using standard duality arguments, a-posteriori error estimates for the discontinuous Galerkin finite element method are then presented for MHD with solenoidal constraint. From these estimates, a procedure for adaptive discretization is outlined. A taxonomy of Green's functions for the linearized MHD operator is given which characterizes the domain of dependence for pointwise errors. The extension to other constrained systems such as the Einstein equations of gravitational physics are then considered. Finally, future directions and open problems are discussed.

Barth, Timothy↗

NASA NTRS2010

A General Algorithm for Reusing Krylov Subspace Information. I. Unsteady Navier-Stokes

A general algorithm is developed that reuses available information to accelerate the iterative convergence of linear systems with multiple right-hand sides A x = b (sup i), which are commonly encountered in steady or unsteady simulations of nonlinear equations. The algorithm is based on the classical GMRES algorithm with eigenvector enrichment but also includes a Galerkin projection preprocessing step and several novel Krylov subspace reuse strategies. The new approach is applied to a set of test problems, including an unsteady turbulent airfoil, and is shown in some cases to provide significant improvement in computational efficiency relative to baseline approaches.

Carpenter, Mark H.↗

DOE OSTI2022

Frequency-domain vs time-domain SP{sub N} equations to simulate neutron noise

Inside the reactor core mechanical vibrations of fuel assemblies can produce high fluctuations around a steady-state configuration, known as neutron noise. This effect can cause the triggering of power reduction measures. Classically, diffusion theory has been used to simulate this behavior. However, this equation has some limitations if the materials of the reactor have strong variations. In this work, we use the diffusive time-dependent simplified spherical harmonics equations that improve the previous results without the necessity of using high computational requirements. In particular, two types of analyses with these equations (SP{sub 3}) are made: a frequency-domain and a time-domain. A numerical neutron noise benchmark tests the methodology and compare both formulations. First, numerical results show a good agreement between the amplitudes and phases of the SP3 equations computed with the frequency-domain and time-domain. Therefore, as the frequency-domain computation only requires to solve a linear system, it is a recommendable option for neutron noise computations. Second, one can conclude that for this type of nuclear systems, where the assemblies are not homogenized, the SP{sub 3} approximation results improve considerably the accuracy of the diffusion theory.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS2018

The Chaotic Long-Term X-Ray Variability of 4U 1705-44

The low-mass X-ray binary 4U1705-44 exhibits dramatic long-term X-ray time variability with a time-scale of several hundred days. The All-SkyMonitor (ASM) aboard the Rossi X-ray Timing Explorer (RXTE) and the Japanese Monitor of All-sky X-ray Image (MAXI) aboard the International Space Station together have continuously observed the source from 1995 December through 2014 May. The combined ASM-MAXI data provide a continuous time series over 50 times the length of the time-scale of interest. Topological analysis can help us identify 'fingerprints' in the phase space of a system unique to its equations of motion. The Birman-Williams theorem postulates that if such fingerprints are the same between two systems, then their equations of motion must be closely related. The phase-space embedding of the source light curve shows a strong resemblance to the double-welled non-linear Duffing oscillator. We explore a range of parameters for which the Duffing oscillator closely mirrors the time evolution of 4U1705-44. We extract low period, unstable periodic orbits from the 4U1705-44 and Duffing time series and compare their topological information. The Duffing and 4U1705-44 topological properties are identical, providing strong evidence that they share the same underlying template. This suggests that we can look to the Duffing equation to help guide the development of a physical model to describe the long-term X-ray variability of this and other similarly behaved X-ray binary systems.

Rossi X-ray Timing Explorer (RXTE)↗

NASA NTRS1987

Legendre-tau approximation for functional differential equations. II - The linear quadratic optimal control problem

The numerical scheme based on the Legendre-tau approximation is proposed to approximate the feedback solution to the linear quadratic optimal control problem for hereditary differential systems. The convergence property is established using Trotter ideas. The method yields very good approximations at low orders and provides an approximation technique for computing closed-loop eigenvalues of the feedback system. A comparison with existing methods (based on averaging and spline approximations) is made.

Ito, Kazufumi↗

NASA NTRS2018

Design of a Modular Monolithic Implicit Solver for Multi-Physics Applications

The design of a modular multi-physics high-order space-time finite-element framework is presented together with its extension to allow monolithic coupling of different physics. One of the main objectives of the framework is to perform efficient high- fidelity simulations of capsule/parachute systems. This problem requires simulating multiple physics including, but not limited to, the compressible Navier-Stokes equations, the dynamics of a moving body with mesh deformations and adaptation, the linear shell equations, non-re effective boundary conditions and wall modeling. The solver is based on high-order space-time - finite element methods. Continuous, discontinuous and C1-discontinuous Galerkin methods are implemented, allowing one to discretize various physical models. Tangent and adjoint sensitivity analysis are also targeted in order to conduct gradient-based optimization, error estimation, mesh adaptation, and flow control, adding another layer of complexity to the framework. The decisions made to tackle these challenges are presented. The discussion focuses first on the "single-physics" solver and later on its extension to the monolithic coupling of different physics. The implementation of different physics modules, relevant to the capsule/parachute system, are also presented. Finally, examples of coupled computations are presented, paving the way to the simulation of the full capsule/parachute system.

Carton De Wiart, Corentin↗

NASA NTRS1973

The global wind system in the thermosphere.

The upper atmosphere wind distribution is calculated by solving the Navier-Stokes equations with models of upper atmospheric densities and ion distribution. All non-linear terms and viscous terms are included. It is found that the velocities from the non-linear solutions are smaller than those obtained from a linear treatment and that in the equatorial region the magnitude of the average meridional velocity is strongly dependent upon the latitudinal coupling. No appreciable super-rotation is obtained by the inclusion of the nonlinear terms; the small value of the average zonal velocity is mainly due to the ion distribution.

Blum, P. W.↗

NASA NTRS1989

An Entropy-Based Approach to Nonlinear Stability

Many numerical methods used in computational fluid dynamics (CFD) incorporate an artificial dissipation term to suppress spurious oscillations and control nonlinear instabilities. The same effect can be accomplished by using upwind techniques, sometimes augmented with limiters to form Total Variation Diminishing (TVD) schemes. An analysis based on numerical satisfaction of the second law of thermodynamics allows many such methods to be compared and improved upon. A nonlinear stability proof is given for discrete scalar equations arising from a conservation law. Solutions to such equations are bounded in the L sub 2 norm if the second law of thermodynamics is satisfied in a global sense over a periodic domain. It is conjectured that an analogous statement is true for discrete equations arising from systems of conservation laws. Analysis and numerical experiments suggest that a more restrictive condition, a positive entropy production rate in each cell, is sufficient to exclude unphysical phenomena such as oscillations and expansion shocks. Construction of schemes which satisfy this condition is demonstrated for linear and nonlinear wave equations and for the one-dimensional Euler equations.

Merriam, Marshal L.↗

NASA NTRS1977

The use of a movable telescoping end mass system for the time-optimal control of spinning spacecraft

The time-optimal control of a spin-stabilized spacecraft with a movable telescoping appendage (boom) is considered analytically and numerically. The motion of a control mass at the end of the boom is determined such that the terminal time will be minimized for two-axis control of a symmetric spacecraft. The equations of rotational motion are linearized about the desired state of spin about the symmetry axis. The equations for the transverse angular velocity components have the form of a coupled two dimensional harmonic oscillator with boom motion as a control force. The control function which brings the system to the desired state is known to be a series of positive and negative pulses. If the initial state is such that the system can be driven to rest in a single switch, the responses, switching and final times, and required boom motion may be determined analytically. Some typical numerical results based on these solutions are discussed.

Bainum, P. M.↗

NASA NTRS1977

Robustness of linear quadratic state feedback designs in the presence of system uncertainty

The paper deals with the problem of expressing the robustness (stability) property of a linear quadratic state feedback (LQSF) design quantitatively in terms of bounds on the perturbations (modeling errors or parameter variations) in the system matrices so that the closed-loop system remains stable. Nonlinear time-varying and linear time-invariant perturbations are considered. The only computation required in obtaining a measure of the robustness of an LQSF design is to determine the eigenvalues of two symmetric matrices determined when solving the algebraic Riccati equation corresponding to the LQSF design problem. Results are applied to a complex dynamic system consisting of the flare control of a STOL aircraft. The design of the flare control is formulated as an LQSF tracking problem.

Patel, R. V.↗

NASA NTRS2006

Stability Limits of a PD Controller for a Flywheel Supported on Rigid Rotor and Magnetic Bearings

Active magnetic bearings are used to provide a long-life, low-loss suspension of a high-speed flywheel rotor. This paper describes a modeling effort used to understand the stability boundaries of the PD controller used to control the active magnetic bearings on a high speed test rig. Limits of stability are described in terms of allowable stiffness and damping values which result in stable levitation of the nonrotating rig. Small signal stability limits for the system is defined as a nongrowth in vibration amplitude of a small disturbance. A simple mass-force model was analyzed. The force resulting from the magnetic bearing was linearized to include negative displacement stiffness and a current stiffness. The current stiffness was then used in a PD controller. The phase lag of the control loop was modeled by a simple time delay. The stability limits and the associated vibration frequencies were measured and compared to the theoretical values. The results show a region on stiffness versus damping plot that have the same qualitative tendencies as experimental measurements. The resulting stability model was then extended to a flywheel system. The rotor dynamics of the flywheel was modeled using a rigid rotor supported on magnetic bearings. The equations of motion were written for the center of mass and a small angle linearization of the rotations about the center of mass. The stability limits and the associated vibration frequencies were found as a function of nondimensional magnetic bearing stiffness and damping and nondimensional parameters of flywheel speed and time delay.

Kascak, Albert F.↗