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 1,045 records · Page 58

A theoretical study of the acoustic impedance of orifices in the presence of a steady grazing flow

An analysis of the oscillatory fluid flow in the vicinity of a circular orifice with a steady grazing flow is presented. The study is similar to that of Hersh and Rogers but with the addition of the grazing flow. Starting from the momentum and continuity equations, a considerably simplified system of partial differential equations is developed with the assumption that the flow can be described by an oscillatory motion superimposed upon the known steady flow. The equations are seen to be linear in the region where the grazing flow effects are dominant, and a solution and the resulting orifice impedance are presented for this region. The nonlinearity appears to be unimportant for the usual conditions found in aircraft noise suppressors. Some preliminary conclusions of the study are that orifice resistance is directly proportional to grazing flow velocity (known previously from experimental data) and that the orifice inductive (mass reactance) end correction is not a function of grazing flow. This latter conclusion is contrary to the widely held notion that grazing flow removes the effect of the orifice inductive end correction. This conclusion also implies that the experimentally observed total inductance reduction with grazing flow might be in the flow within the orifice rather than in the end correction.

Rice, E. J.↗

Euclidean controllability of linear delay systems with limited controls

For linear delay systems with limited controls the notion of a proper control system is introduced. If the uncontrolled system x(t) = Ax(t) + Bx(t - h) is uniformly asymptotically stable and the control equation x(t) = Ax(t) + Bx(t-) + Cu(t) is proper, then the control system is Euclidean null-controllable. It is noted that controllability is equivalent to the system being proper for delay systems with unlimited power.

Chukwu, E. N.↗

The well-posedness of the Kuramoto-Sivashinsky equation

The Kuramoto-Sivashinsky equation arises in a variety of applications, among which are modeling reaction diffusion systems, flame propagation and viscous flow problems. It is considered here, as a prototype to the larger class of generalized Burgers equations: those consist of a quadratic nonlinearity and an arbitrary linear parabolic part. It is shown that such equations are well posed, thus admitting a unique smooth solution, continuously dependent on its initial data. As an attractive alternative to standard energy methods, existence and stability are derived in this case, by patching in the large short time solutions without loss of derivatives.

Tadmor, E.↗

Characteristics pertaining to a stiffness cross-coupled Jeffcott model

Rotordynamic studies of complex systems utilizing multiple degree-of-freedom analysis have been performed to understand response, loads, and stability. In order to understand the fundamental nature of rotordynamic response, the Jeffcott rotor model has received wide attention. The purpose of this paper is to provide a generic rotordynamic analysis of a stiffness cross-coupled Jeffcott rotor model to illustrate characteristics of a second order stiffness-coupled linear system. The particular characteristics investigated were forced response, force vector diagrams, response orbits, and stability. Numerical results were achieved through a fourth order Runge-Kutta method for solving differential equations and the Routh Hurwitz stability criterion. The numerical results were verified to an exact mathematical solution for the steady state response.

Spanyer, K. L.↗

The well-posedness of the Kuramoto-Sivashinsky equation

The Kuramoto-Sivashinsky equation arises in a variety of applications, among which are modeling reaction diffusion systems, flame propagation and viscous flow problems. It is considered here, as a prototype to the larger class of generalized Burgers equations: those consist of a quadratic nonlinearity and an arbitrary linear parabolic part. It is shown that such equations are well posed, thus admitting a unique smooth solution, continuously dependent on its initial data. As an attractive alternative to standard energy methods, existence and stability are derived in this case, by patching in the large short time solutions without 'loss of derivatives'.

Tadmor, E.↗

Piecewise linear approximation for hereditary control problems

Finite dimensional approximations are presented for linear retarded functional differential equations by use of discontinuous piecewise linear functions. The approximation scheme is applied to optimal control problems when a quadratic cost integral has to be minimized subject to the controlled retarded system. It is shown that the approximate optimal feedback operators converge to the true ones both in case the cost integral ranges over a finite time interval as well as in the case it ranges over an infinite time interval. The arguments in the latter case rely on the fact that the piecewise linear approximations to stable systems are stable in a uniform sense. This feature is established using a vector-component stability criterion in the state space R(n) x L(2) and the favorable eigenvalue behavior of the piecewise linear approximations.

Propst, Georg↗

The response of a laminar boundary layer in supersonic flow to small amplitude progressive waves

The effect of a small amplitude progressive wave on the laminar boundary layer on a semi-infinite flat plate, due to a uniform supersonic freestream flow, is considered. The perturbation to the flow divides into two streamwise zones. In the first, relatively close to the leading edge of the plate, on a transverse scale comparable to the boundary layer thickness, the perturbation flow is described by a form of the unsteady linearized compressible boundary layer equations. In the freestream, this component of flow is governed by the wave equation, the solution of which provides the outer velocity conditions for the boundary layer. This system is solved numerically, and also the asymptotic structure in the far downstream limit is studied. This reveals a breakdown and a subsequent second streamwise zone, where the flow disturbance is predominantly inviscid. The two zones are shown to match in a proper asymptotic sense.

Duck, Peter W.↗

Computing Linear Mathematical Models Of Aircraft

Derivation and Definition of Linear Aircraft Model (LINEAR) computer program provides user with powerful, and flexible, standard, documented, and verified software tool for linearization of mathematical models of aerodynamics of aircraft. Intended for use in software tool to drive linear analysis of stability and design of control laws for aircraft. Capable of both extracting such linearized engine effects as net thrust, torque, and gyroscopic effects, and including these effects in linear model of system. Designed to provide easy selection of state, control, and observation variables used in particular model. Also provides flexibility of allowing alternate formulations of both state and observation equations. Written in FORTRAN.

Duke, Eugene L.↗

Theory of finite disturbances in a centrifugal compression system with a vaneless radial diffuser

A previous small perturbation analysis of circumferential waves in circumferential compression systems, assuming inviscid flow, is shown to be consistent with observations that narrow diffusers are more stable than wide ones, when boundary layer displacement effect is included. The Moore-Greitzer analysis for finite strength transients containing both surge and rotating stall in axial machines is adapted for a centrifugal compression system. Under certain assumptions, and except for a new second order swirl, the diffuser velocity field, including resonant singularities, can be carried over from the previous inviscid linear analysis. Nonlinear transient equations are derived and applied in a simple example to show that throttling through a resonant value of flow coefficient must occur in a sudden surge-like drop, accompanied by a transient rotating wave. This inner solution is superseded by an outer surge response on a longer time scale. Surge may occur purely as result of circumferential wave resonance. Numerical results are shown for various parametric choices relating to throttle schedule and the characteristic slope. A number of circumferential modes considered simultaneously is briefly discussed.

Moore, F. K.↗

Attitude Determination Error Analysis System (ADEAS) mathematical specifications document

The mathematical specifications of Release 4.0 of the Attitude Determination Error Analysis System (ADEAS), which provides a general-purpose linear error analysis capability for various spacecraft attitude geometries and determination processes, are presented. The analytical basis of the system is presented. The analytical basis of the system is presented, and detailed equations are provided for both three-axis-stabilized and spin-stabilized attitude sensor models.

Nicholson, Mark↗

A homotopy algorithm for digital optimal projection control GASD-HADOC

The linear-quadratic-gaussian (LQG) compensator was developed to facilitate the design of control laws for multi-input, multi-output (MIMO) systems. The compensator is computed by solving two algebraic equations for which standard closed-loop solutions exist. Unfortunately, the minimal dimension of an LQG compensator is almost always equal to the dimension of the plant and can thus often violate practical implementation constraints on controller order. This deficiency is especially highlighted when considering control-design for high-order systems such as flexible space structures. This deficiency motivated the development of techniques that enable the design of optimal controllers whose dimension is less than that of the design plant. A homotopy approach based on the optimal projection equations that characterize the necessary conditions for optimal reduced-order control. Homotopy algorithms have global convergence properties and hence do not require that the initializing reduced-order controller be close to the optimal reduced-order controller to guarantee convergence. However, the homotopy algorithm previously developed for solving the optimal projection equations has sublinear convergence properties and the convergence slows at higher authority levels and may fail. A new homotopy algorithm for synthesizing optimal reduced-order controllers for discrete-time systems is described. Unlike the previous homotopy approach, the new algorithm is a gradient-based, parameter optimization formulation and was implemented in MATLAB. The results reported may offer the foundation for a reliable approach to optimal, reduced-order controller design.

Collins, Emmanuel G., Jr.↗

Projection techniques for iterative solution of A(bar x) equals (bar b) with successive right-hand sides

Two projection techniques for computing approximate solutions to linear systems of the form A(bar x)(sup n) = (bar b)(sup n), for a sequence n = 1, 2, ..., e.g., such as arises from time discretization of a partial differential equation, are presented. The inexpensive approximate solutions can be used as initial guesses for iterative solution of the system, resulting in significantly reduced computational expense. Examples of two- and three-dimensional incompressible Navier-Stokes calculations are presented in which x represents the pressure, and A is a discrete Poisson operator. In flows containing significant dynamic activity, these projection techniques lead to as much as a two-fold reduction in solution time.

Fischer, Paul F.↗

CGRO Guest Investigator Program

The following are highlights from the research supported by this grant: (1) Theory of gamma-ray blazars: We studied the theory of gamma-ray blazars, being among the first investigators to propose that the GeV emission arises from Comptonization of diffuse radiation surrounding the jet, rather than from the synchrotron-self-Compton mechanism. In related work, we uncovered possible connections between the mechanisms of gamma-ray blazars and those of intraday radio variability, and have conducted a general study of the role of Compton radiation drag on the dynamics of relativistic jets. (2) A Nonlinear Monte Carlo code for gamma-ray spectrum formation: We developed, tested, and applied the first Nonlinear Monte Carlo (NLMC) code for simulating gamma-ray production and transfer under much more general (and realistic) conditions than are accessible with other techniques. The present version of the code is designed to simulate conditions thought to be present in active galactic nuclei and certain types of X-ray binaries, and includes the physics needed to model thermal and nonthermal electron-positron pair cascades. Unlike traditional Monte-Carlo techniques, our method can accurately handle highly non-linear systems in which the radiation and particle backgrounds must be determined self-consistently and in which the particle energies span many orders of magnitude. Unlike models based on kinetic equations, our code can handle arbitrary source geometries and relativistic kinematic effects In its first important application following testing, we showed that popular semi-analytic accretion disk corona models for Seyfert spectra are seriously in error, and demonstrated how the spectra can be simulated if the disk is sparsely covered by localized 'flares'.

Begelman, Mitchell C.↗

Active Vibration Control of Elastic Beam by Means of Shape Memory Alloy Layers

The mathematical model of a flexible beam covered with shape memory alloy (SMA) layers is presented. The SMA layers are used as actuators, which are capable of changing their elastic modulus and recovery stress, thus changing the natural frequency of, and adjusting the excitation to, the vibrating beam. The frequency factor variation as a function of SMA Young's modulus, SMA layer thickness and beam thickness is discussed. Also control of the beam employing an optimal linear control law is evaluated. The control results indicate how the system reacts to various levels of excitation input through the non-homogeneous recovery shear term of the governing differential equation.

Chen, Q.↗

System IDentification Programs for AirCraft (SIDPAC)

A collection of computer programs for aircraft system identification is described and demonstrated. The programs, collectively called System IDentification Programs for AirCraft, or SIDPAC, were developed in MATLAB as m-file functions. SIDPAC has been used successfully at NASA Langley Research Center with data from many different flight test programs and wind tunnel experiments. SIDPAC includes routines for experiment design, data conditioning, data compatibility analysis, model structure determination, equation-error and output-error parameter estimation in both the time and frequency domains, real-time and recursive parameter estimation, low order equivalent system identification, estimated parameter error calculation, linear and nonlinear simulation, plotting, and 3-D visualization. An overview of SIDPAC capabilities is provided, along with a demonstration of the use of SIDPAC with real flight test data from the NASA Glenn Twin Otter aircraft. The SIDPAC software is available without charge to U.S. citizens by request to the author, contingent on the requestor completing a NASA software usage agreement.

Morelli, Eugene A.↗

Analysis Method for Inelastic, Adhesively Bonded Joints with Anisotropic Adherends

A one-dimensional analysis method for evaluating adhesively bonded joints composed of anisotropic adherends and adhesives with nonlinear material behavior is presented in the proposed paper. The strain and resulting stress field in a general, bonded joint overlap are determined by using a variable-step, finite-difference solution algorithm to iteratively solve a system of first-order differential equations. Applied loading is given by a system of combined extensional, bending, and shear forces that are applied to the edge of the joint overlap. Adherends are assumed to behave as linear, cylindrically bent plates using classical laminated plate theory that includes the effects of first-order transverse shear deformation. Using the deformation theory of plasticity and a modified von-Mises yield criterion, inelastic material behavior is modeled in the adhesive layer. Results for the proposed method are verified against previous results from the literature and shown to be in excellent agreement. An additional case that highlights the effects of transverse shear deformation between similar adherends is also presented.

Smeltzer, Stanley S., III↗

Reduced Order Methods for Prediction of Thermal-Acoustic Fatigue

The goal of this investigation is to assess the quality of high-cycle-fatigue life estimation via a reduced order method, for structures undergoing random nonlinear vibrations in a presence of thermal loading. Modal reduction is performed with several different suites of basis functions. After numerically solving the reduced order system equations of motion, the physical displacement time history is obtained by an inverse transformation and stresses are recovered. Stress ranges obtained through the rainflow counting procedure are used in a linear damage accumulation method to yield fatigue estimates. Fatigue life estimates obtained using various basis functions in the reduced order method are compared with those obtained from numerical simulation in physical degrees-of-freedom.

Przekop, A.↗

A Nonlinear Reduced Order Method for Prediction of Acoustic Fatigue

The goal of this investigation is to assess the quality of high-cycle-fatigue life estimation via a reduced order method, for structures undergoing geometrically nonlinear random vibrations. Modal reduction is performed with several different suites of basis functions. After numerically solving the reduced order system equations of motion, the physical displacement time history is obtained by an inverse transformation and stresses are recovered. Stress ranges obtained through the rainflow counting procedure are used in a linear damage accumulation method to yield fatigue estimates. Fatigue life estimates obtained using various basis functions in the reduced order method are compared with those obtained from numerical simulation in physical degrees-of-freedom.

Przekop, Adam↗