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,189 records · Page 66

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Global dynamic modeling of a transmission system

The work performed on global dynamic simulation and noise correlation of gear transmission systems at the University of Akron is outlined. The objective is to develop a comprehensive procedure to simulate the dynamics of the gear transmission system coupled with the effects of gear box vibrations. The developed numerical model is benchmarked with results from experimental tests at NASA Lewis Research Center. The modal synthesis approach is used to develop the global transient vibration analysis procedure used in the model. Modal dynamic characteristics of the rotor-gear-bearing system are calculated by the matrix transfer method while those of the gear box are evaluated by the finite element method (NASTRAN). A three-dimensional, axial-lateral coupled bearing model is used to couple the rotor vibrations with the gear box motion. The vibrations between the individual rotor systems are coupled through the nonlinear gear mesh interactions. The global equations of motion are solved in modal coordinates and the transient vibration of the system is evaluated by a variable time-stepping integration scheme. The relationship between housing vibration and resulting noise of the gear transmission system is generated by linear transfer functions using experimental data. A nonlinear relationship of the noise components to the fundamental mesh frequency is developed using the hypercoherence function. The numerically simulated vibrations and predicted noise of the gear transmission system are compared with the experimental results from the gear noise test rig at NASA Lewis Research Center. Results of the comparison indicate that the global dynamic model developed can accurately simulate the dynamics of a gear transmission system.

Choy, F. K.↗

NASA NTRS2006

Adjoint Method and Predictive Control for 1-D Flow in NASA Ames 11-Foot Transonic Wind Tunnel

This paper describes a modeling method and a new optimal control approach to investigate a Mach number control problem for the NASA Ames 11-Foot Transonic Wind Tunnel. The flow in the wind tunnel is modeled by the 1-D unsteady Euler equations whose boundary conditions prescribe a controlling action by a compressor. The boundary control inputs to the compressor are in turn controlled by a drive motor system and an inlet guide vane system whose dynamics are modeled by ordinary differential equations. The resulting Euler equations are thus coupled to the ordinary differential equations via the boundary conditions. Optimality conditions are established by an adjoint method and are used to develop a model predictive linear-quadratic optimal control for regulating the Mach number due to a test model disturbance during a continuous pitch

Nguyen, Nhan↗

NASA NTRS2009

Modeling and Simulation of Variable Mass, Flexible Structures

The advent of the new Ares I launch vehicle has highlighted the need for advanced dynamic analysis tools for variable mass, flexible structures. This system is composed of interconnected flexible stages or components undergoing rapid mass depletion through the consumption of solid or liquid propellant. In addition to large rigid body configuration changes, the system simultaneously experiences elastic deformations. In most applications, the elastic deformations are compatible with linear strain-displacement relationships and are typically modeled using the assumed modes technique. The deformation of the system is approximated through the linear combination of the products of spatial shape functions and generalized time coordinates. Spatial shape functions are traditionally composed of normal mode shapes of the system or even constraint modes and static deformations derived from finite element models of the system. Equations of motion for systems undergoing coupled large rigid body motion and elastic deformation have previously been derived through a number of techniques [1]. However, in these derivations, the mode shapes or spatial shape functions of the system components were considered constant. But with the Ares I vehicle, the structural characteristics of the system are changing with the mass of the system. Previous approaches to solving this problem involve periodic updates to the spatial shape functions or interpolation between shape functions based on system mass or elapsed mission time. These solutions often introduce misleading or even unstable numerical transients into the system. Plus, interpolation on a shape function is not intuitive. This paper presents an approach in which the shape functions are held constant and operate on the changing mass and stiffness matrices of the vehicle components. Each vehicle stage or component finite element model is broken into dry structure and propellant models. A library of propellant models is used to describe the distribution of mass in the fuel tank or Solid Rocket Booster (SRB) case for various propellant levels. Based on the mass consumed by the liquid engine or SRB, the appropriate propellant model is coupled with the dry structure model for the stage. Then using vehicle configuration data, the integrated vehicle model is assembled and operated on by the constant system shape functions. The system mode shapes and frequencies can then be computed from the resulting generalized mass and stiffness matrices for that mass configuration. The rigid body mass properties of the vehicle are derived from the integrated vehicle model. The coupling terms between the vehicle rigid body motion and elastic deformation are also updated from the constant system shape functions and the integrated vehicle model. This approach was first used to analyze variable mass spinning beams and then prototyped into a generic dynamics simulation engine. The resulting code was tested against Crew Launch Vehicle (CLV-)class problems worked in the TREETOPS simulation package and by Wilson [2]. The Ares I System Integration Laboratory (SIL) is currently being developed at the Marshall Space Flight Center (MSFC) to test vehicle avionics hardware and software in a hardware-in-the-loop (HWIL) environment and certify that the integrated system is prepared for flight. The Ares I SIL utilizes the Ares Real-Time Environment for Modeling, Integration, and Simulation (ARTEMIS) tool to simulate the launch vehicle and stimulate avionics hardware. Due to the presence of vehicle control system filters and the thrust oscillation suppression system, which are tuned to the structural characteristics of the vehicle, ARTEMIS must incorporate accurate structural models of the Ares I launch vehicle. The ARTEMIS core dynamics simulation models the highly coupled nature of the vehicle flexible body dynamics, propellant slosh, and vehicle nozzle inertia effects combined with mass and flexible body properties that vary significant with time during the flight. All forces that act on the vehicle during flight must be simulated, including deflected engine thrust force, spatially distributed aerodynamic forces, gravity, and reaction control jet thrust forces. These forces are used to excite an integrated flexible vehicle, slosh, and nozzle dynamics model for the vehicle stack that simulates large rigid body translations and rotations along with small elastic deformations. Highly effective matrix math operations on a distributed, threaded high-performance simulation node allow ARTEMIS to retain up to 30 modes of flex for real-time simulation. Stage elements that separate from the stack during flight are propagated as independent rigid six degrees of freedom (6DOF) bodies. This paper will present the formulation of the resulting equations of motion, solutions to example problems, and describe the resulting dynamics simulation engine within ARTEMIS.

Tobbe, Patrick A.↗

NASA NTRS1966

Theory and Practical Design of Phase-locked Receivers, Volume I

This is Volume I of a two-volume work on the theory, of phase-locked receivers, with pertinent reference material on practical receiver design. Volume I is primarily devoted to the performance of carrier-tracking loops, including a rigorous treatment of narrow-band systems having IF limiters. The bulk of the work is based on a theoretical linear model, but a nonlinear method is also presented to predict behavior near the threshold. The emphasis throughout the work is toward completeness, simplicity, and internal consistency of the material assembled. Part I of this Volume is an exposition of the theory, the resulting equations, and the design philosophy that have enabled the phase-lock concept to evolve into the basic principle underlying the most sensitive receivers il, the world today. Part II is a condensed version of Part I, intended as a quick reference to formulas, definitions, and salient design considerations.

RECEIVER↗

NASA NTRS1993

Global organization of tectonic deformation on Venus

The geographic organization of surface deformation on Venus as on Earth is a key to understanding the global tectonic system. To date we have mapped the distribution of three unambiguous tectonic land forms on Venus: (1) linear foldbelts analogous to those at plate margins of the Earth; (2) linear rift zones, analogous to continental rifts on the Earth; and (3) distributed plains deformation in the form of wrinkle ridges and extensional faults and fractures. The linear foldbelts are the dominant structural style in the Northern Hemisphere; ninety percent of the planet's foldbelts lie above the equator. In contrast, compressive deformation in the Southern Hemisphere is dominated by two large, sweeping patterns of wrinkle ridges. The two hemispheres are divided by an equatorial region that is largely covered by rift zones and several large tessera blocks. A tectonic model of generally poleward convergence of the Northern Hemisphere explains the distribution of foldbelts and rift zones. In our model, a northern hemispherical plate (or system of plates) moves poleward and deforms along discrete, predominately longitudinal bands. We recognize four types of foldbelts based on their relationships to other large-scale tectonic features on Venus. There are foldbelts that lie within the low plains, foldbelts associated with coronae, novae and chasmata, foldbelts that lie at the margins of poly-deformed tessera plateaus, and the folded mountain belts around Lakshmi Planum. We see a geometric increase in the area of fold belts when normalized to percent area at a given latitude. This increase is consistent with our model of poleward convergence. Also, the orientations of most foldbelts are either approximately north-south or parallel to lines of latitude in the northern hemisphere. This observation is also consistent with the model in that the longitudinal bands are the result of the decreasing area of the sphere as the plate moves poleward and the latitudinal belts are the direct result of poleward compression. The trends of wrinkle ridges have been mapped over the planet and several large, sweeping patterns evidently reflect long-wavelength topography. Using wrinkle ridges as paleostress indicators, we have developed local and regional stress trajectory maps.

Bilotti, Frank↗

NASA NTRS1974

Power-spectral-density relationship for retarded differential equations

The power spectral density (PSD) relationship between input and output of a set of linear differential-difference equations of the retarded type with real constant coefficients and delays is discussed. The form of the PSD relationship is identical with that applicable to unretarded equations. Since the PSD relationship is useful if and only if the system described by the equations is stable, the stability must be determined before applying the PSD relationship. Since it is sometimes difficult to determine the stability of retarded equations, such equations are often approximated by simpler forms. It is pointed out that some common approximations can lead to erroneous conclusions regarding the stability of a system and, therefore, to the possibility of obtaining PSD results which are not valid.

Barker, L. K.↗

NASA NTRS1976

Scaled equation of state parameters for gases in the critical region

In the light of recent theoretical developments, the paper presents an accurate characterization of anomalous thermodynamic behavior of xenon, helium 4, helium 3, carbon dioxide, steam and oxygen in the critical region. This behavior is associated with long range fluctuations in the system and the physical properties depend primarily on a single variable, namely, the correlation length. A description of the thermodynamic behavior of fluids in terms of scaling laws is formulated, and the two successfully used scaled equations of state (NBS equation and Linear Model parametric equation) are compared. Methods for fitting both equations to experimental equation of state data are developed and formulated, and the optimum fit for each of the two scaled equations of the above gases are presented and the results are compared. By extending the experimental data for the above one-component fluids to partially miscible binary liquids, superfluid liquid helium, ferromagnets and solids exhibiting order-disorder transitions, the principle of universality is concluded. Finally by using this principle, the critical regions for nine additional fluids are described.

Sengers, J. M. H. L.↗

NASA NTRS1979

Space construction base control system

Several approaches for an attitude control system are studied and developed for a large space construction base that is structurally flexible. Digital simulations were obtained using the following techniques: (1) the multivariable Nyquist array method combined with closed loop pole allocation, (2) the linear quadratic regulator method. Equations for the three-axis simulation using the multilevel control method were generated and are presented. Several alternate control approaches are also described. A technique is demonstrated for obtaining the dynamic structural properties of a vehicle which is constructed of two or more submodules of known dynamic characteristics.

Kaczynski, R. F.↗

NASA NTRS1991

Computational attributes of the integral form of the equation of transfer

Difficulties can arise in radiative and neutron transport calculations when a highly anisotropic scattering phase function is present. In the presence of anisotropy, currently used numerical solutions are based on the integro-differential form of the linearized Boltzmann transport equation. This paper, departs from classical thought and presents an alternative numerical approach based on application of the integral form of the transport equation. Use of the integral formalism facilitates the following steps: a reduction in dimensionality of the system prior to discretization, the use of symbolic manipulation to augment the computational procedure, and the direct determination of key physical quantities which are derivable through the various Legendre moments of the intensity. The approach is developed in the context of radiative heat transfer in a plane-parallel geometry, and results are presented and compared with existing benchmark solutions. Encouraging results are presented to illustrate the potential of the integral formalism for computation. The integral formalism appears to possess several computational attributes which are well-suited to radiative and neutron transport calculations.

Frankel, J. I.↗

NASA NTRS1993

Recent developments in high order K-exact reconstruction on unstructured meshes

This paper presents recent improvements in high-order K-exact reconstruction on unstructured meshes. The new reconstruction procedures are incorporated into a basic upwind finite-volume scheme suitable for solving scalar advection-diffusion equations as well as the Euler and Navier-Stokes equations. Numerical calculations are performed comparing the present method with lower order accurate reconstruction procedures (piecewise constant and piecewise linear) and various competing technologies such as the fluctuation splitting method of Roe (1987) and Deconinck et al. (1992) and a system-variant of the streamline diffusion Petrov-Galerkin method developed by Hansbo (1991) and Hansbo and Johnson (1991). Five test problems are used in the numerical comparisons: scalar circular advection, transonic and supersonic Euler flow, laminar boundary-layer flow, and general compressible Navier-Stokes flow.

Barth, Timothy J.↗

NASA NTRS1993

Spectral response of a bilinear oscillator

An approximate analytical procedure is presented to estimate the response power spectral density of a randomly excited spring/mass/damper system having a bilinear spring. The approximate expression for the response spectrum is developed by representing the nonlinear oscillator as a linear system having a natural frequency that depends on the envelope of the random response. This approximate representation of the system leads to estimates of the response spectrum that agree extremely well with those obtained by direct numerical simulation of the governing equation.

Miles, R. N.↗

NASA NTRS1994

Linear and nonlinear dynamic analysis of redundant load path bearingless rotor systems

The goal of this research is to develop the transfer matrix method to treat nonlinear autonomous boundary value problems with multiple branches. The application is the complete nonlinear aeroelastic analysis of multiple-branched rotor blades. Once the development is complete, it can be incorporated into the existing transfer matrix analyses. There are several difficulties to be overcome in reaching this objective. The conventional transfer matrix method is limited in that it is applicable only to linear branch chain-like structures, but consideration of multiple branch modeling is important for bearingless rotors. Also, hingeless and bearingless rotor blade dynamic characteristics (particularly their aeroelasticity problems) are inherently nonlinear. The nonlinear equations of motion and the multiple-branched boundary value problem are treated together using a direct transfer matrix method. First, the formulation is applied to a nonlinear single-branch blade to validate the nonlinear portion of the formulation. The nonlinear system of equations is iteratively solved using a form of Newton-Raphson iteration scheme developed for differential equations of continuous systems. The formulation is then applied to determine the nonlinear steady state trim and aeroelastic stability of a rotor blade in hover with two branches at the root. A comprehensive computer program is developed and is used to obtain numerical results for the (1) free vibration, (2) nonlinearly deformed steady state, (3) free vibration about the nonlinearly deformed steady state, and (4) aeroelastic stability tasks. The numerical results obtained by the present method agree with results from other methods.

Murthy, V. R.↗

NASA NTRS1992

Ballistic limit regression analysis for Space Station Freedom meteoroid and space debris protection system

Relationships defining the ballistic limit of Space Station Freedom's (SSF) dual wall protection systems have been determined. These functions were regressed from empirical data found in Marshall Space Flight Center's (MSFC) Hypervelocity Impact Testing Summary (HITS) for the velocity range between three and seven kilometers per second. A stepwise linear least squares regression was used to determine the coefficients of several expressions that define a ballistic limit surface. Using statistical significance indicators and graphical comparisons to other limit curves, a final set of expressions is recommended for potential use in Probability of No Critical Flaw (PNCF) calculations for Space Station. The three equations listed below represent the mean curves for normal, 45 degree, and 65 degree obliquity ballistic limits, respectively, for a dual wall protection system consisting of a thin 6061-T6 aluminum bumper spaced 4.0 inches from a .125 inches thick 2219-T87 rear wall with multiple layer thermal insulation installed between the two walls. Normal obliquity is d(sub c) = 1.0514 v(exp 0.2983 t(sub 1)(exp 0.5228). Forty-five degree obliquity is d(sub c) = 0.8591 v(exp 0.0428) t(sub 1)(exp 0.2063). Sixty-five degree obliquity is d(sub c) = 0.2824 v(exp 0.1986) t(sub 1)(exp -0.3874). Plots of these curves are provided. A sensitivity study on the effects of using these new equations in the probability of no critical flaw analysis indicated a negligible increase in the performance of the dual wall protection system for SSF over the current baseline. The magnitude of the increase was 0.17 percent over 25 years on the MB-7 configuration run with the Bumper II program code.

Jolly, William H.↗

NASA NTRS1971

A numerical solution of the Navier-Stokes equations for supercritical fluid thermodynamic analysis

An explicit numerical solution of the compressible Navier-Stokes equations is applied to the thermodynamic analysis of supercritical oxygen in the Apollo cryogenic storage system. The wave character is retained in the conservation equations which are written in the basic fluid variables for a two-dimensional Cartesian coordinate system. Control-volume cells are employed to simplify imposition of boundary conditions and to ensure strict observance of local and global conservation principles. Non-linear real-gas thermodynamic properties responsible for the pressure collapse phenomonon in supercritical fluids are represented by tabular and empirical functions relating pressure and temperature to density and internal energy. Wall boundary conditions are adjusted at one cell face to emit a prescribed mass flowrate. Scaling principles are invoked to achieve acceptable computer execution times for very low Mach number convection problems. Detailed simulations of thermal stratification and fluid mixing occurring under low acceleration in the Apollo 12 supercritical oxygen tank are presented which model the pressure decay associated with de-stratification induced by an ordinary vehicle maneuver and heater cycle operation.

Heinmiller, P. J.↗

NASA NTRS1972

Stability of a steady, large amplitude whistler wave.

Study of the behavior of weak electrostatic waves in a collisionless magnetoplasma supporting a steady large amplitude whistler wave. All waves are assumed to propagate parallel to a uniform background magnetic field B sub zero. In the presence of the whistler wave fields each particle executes an oscillatory motion parallel to B sub zero, in addition to a translation along B sub zero and transverse motions. This oscillation causes the Landau resonance to be replaced by a series of new resonances between particles and the electrostatic modes. A distribution function for the perturbed plasma is constructed by solving the Vlasov equation, linearized in the electrostatic wave amplitudes. A dispersion relation is obtained and solved approximately for the growth/damping rate of the perturbations. Growing electrostatic modes are found to be approximately uncoupled. Trapped particles have a strong influence on the stability of the system.

Palmadesso, P. J.↗

NASA NTRS1976

Stability and Relative Stability of Linear Systems with Many Constant Time Delays

A method of determining the stability of linear systems with many constant time delays is developed. This technique, an extension of the tau-decomposition method, is used to examine not only the stability but also the relative stability of retarded systems with many delays and a class of neutral equations with one delay. Analytical equations are derived for partitioning the delay space of a retarded system with two time delays. The stability of the system in each of the regions defined by the partitioning curves in the parameter plane is determined using the extended tau-decomposition method. In addition, relative stability boundaries are defined using the extended tau-decompositon method in association with parameter plane techniques. Several applications of the extended tau-decomposition method are presented and compared with stability results obtained from other analyses. In all cases the results obtained using the method outlined herein coincide with and extend those of previous investigations. The extended tau-decomposition method applied to systems with time delays requires less computational effort and yields more complete stability analyses than previous techniques.

Barker, Larry Keith↗

NASA NTRS1979

Comparative study of flare control laws

A digital 3-D automatic control law was developed to achieve an optimal transition of a B-737 aircraft between various initial glid slope conditions and the desired final touchdown condition. A discrete, time-invariant, optimal, closed-loop control law presented for a linear regulator problem, was extended to include a system being acted upon by a constant disturbance. Two forms of control laws were derived to solve this problem. One method utilized the feedback of integral states defined appropriately and augmented with the original system equations. The second method formulated the problem as a control variable constraint, and the control variables were augmented with the original system. The control variable constraint control law yielded a better performance compared to feedback control law for the integral states chosen.

Nadkarni, A. A.↗

NASA NTRS1990

General Rotorcraft Aeromechanical Stability Program (GRASP): Theory manual

The general rotorcraft aeromechanical stability program (GRASP) was developed to calculate aeroelastic stability for rotorcraft in hovering flight, vertical flight, and ground contact conditions. GRASP is described in terms of its capabilities and its philosophy of modeling. The equations of motion that govern the physical system are described, as well as the analytical approximations used to derive them. The equations include the kinematical equation, the element equations, and the constraint equations. In addition, the solution procedures used by GRASP are described. GRASP is capable of treating the nonlinear static and linearized dynamic behavior of structures represented by arbitrary collections of rigid-body and beam elements. These elements may be connected in an arbitrary fashion, and are permitted to have large relative motions. The main limitation of this analysis is that periodic coefficient effects are not treated, restricting rotorcraft flight conditions to hover, axial flight, and ground contact. Instead of following the methods employed in other rotorcraft programs. GRASP is designed to be a hybrid of the finite-element method and the multibody methods used in spacecraft analysis. GRASP differs from traditional finite-element programs by allowing multiple levels of substructure in which the substructures can move and/or rotate relative to others with no small-angle approximations. This capability facilitates the modeling of rotorcraft structures, including the rotating/nonrotating interface and the details of the blade/root kinematics for various types. GRASP differs from traditional multibody programs by considering aeroelastic effects, including inflow dynamics (simple unsteady aerodynamics) and nonlinear aerodynamic coefficients.

Hodges, Dewey H.↗