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,693 records · Page 94

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Effect Of Noise In The Ideal State Reconstructor

More measurements yield better estimate. Report discusses effects of measurement noise on system including linear, time-invariant plant of measurable inputs and known parameters and deterministic digital control subsystem governed by algorithm called "ideal state reconstructor." So named because in absence of noise, it exactly reconstructs vector represting state of plant, even without knowledge of initial state. Ideal state reconstructor adds no new states or eigenvalues to system; affects measurement equation only.

Polites, Michael E.↗

NASA NTRS1989

Numerical studies of identification in nonlinear distributed parameter systems

An abstract approximation framework and convergence theory for the identification of first and second order nonlinear distributed parameter systems developed previously by the authors and reported on in detail elsewhere are summarized and discussed. The theory is based upon results for systems whose dynamics can be described by monotone operators in Hilbert space and an abstract approximation theorem for the resulting nonlinear evolution system. The application of the theory together with numerical evidence demonstrating the feasibility of the general approach are discussed in the context of the identification of a first order quasi-linear parabolic model for one dimensional heat conduction/mass transport and the identification of a nonlinear dissipation mechanism (i.e., damping) in a second order one dimensional wave equation. Computational and implementational considerations, in particular, with regard to supercomputing, are addressed.

Banks, H. T.↗

NASA NTRS1993

Annular suspension and pointing system with controlled DC electromagnets

The Annular Suspension and Pointing System (ASPS) developed by the Flight System division of Sperry Corporation is a six-degree of freedom payload pointing system designed for use with the space shuttle. This magnetic suspension and pointing system provides precise controlled pointing in six-degrees of freedom, isolation of payload-carrier disturbances, and end mount controlled pointing. Those are great advantages over the traditional mechanical joints for space applications. In this design, we first analyzed the assumed model of the single degree ASPS bearing actuator and obtained the plant dynamics equations. By linearizing the plant dynamics equations, we designed the cascade and feedback compensators such that a stable and satisfied result was obtained. The specified feedback compensator was computer simulated with the nonlinearized plant dynamics equations. The results indicated that an unstable output occurred. In other words, the designed feedback compensator failed. The failure of the design is due to the Taylor's series expansion not converging.

Vu, Josephine Lynn↗

NASA NTRS2002

Field-Programmable Gate Array Computer in Structural Analysis: An Initial Exploration

This paper reports on an initial assessment of using a Field-Programmable Gate Array (FPGA) computational device as a new tool for solving structural mechanics problems. A FPGA is an assemblage of binary gates arranged in logical blocks that are interconnected via software in a manner dependent on the algorithm being implemented and can be reprogrammed thousands of times per second. In effect, this creates a computer specialized for the problem that automatically exploits all the potential for parallel computing intrinsic in an algorithm. This inherent parallelism is the most important feature of the FPGA computational environment. It is therefore important that if a problem offers a choice of different solution algorithms, an algorithm of a higher degree of inherent parallelism should be selected. It is found that in structural analysis, an 'analog computer' style of programming, which solves problems by direct simulation of the terms in the governing differential equations, yields a more favorable solution algorithm than current solution methods. This style of programming is facilitated by a 'drag-and-drop' graphic programming language that is supplied with the particular type of FPGA computer reported in this paper. Simple examples in structural dynamics and statics illustrate the solution approach used. The FPGA system also allows linear scalability in computing capability. As the problem grows, the number of FPGA chips can be increased with no loss of computing efficiency due to data flow or algorithmic latency that occurs when a single problem is distributed among many conventional processors that operate in parallel. This initial assessment finds the FPGA hardware and software to be in their infancy in regard to the user conveniences; however, they have enormous potential for shrinking the elapsed time of structural analysis solutions if programmed with algorithms that exhibit inherent parallelism and linear scalability. This potential warrants further development of FPGA-tailored algorithms for structural analysis.

Singleterry, Robert C., Jr.↗

NASA NTRS1949

On the Theory of the Laval Nozzle

In the present paper, the motion of a gas in a plane-parallel Laval nozzle in the neighborhood of the transition from subsonic to supersonic velocities is studied. In a recently published paper, F. I. Frankl, applying the holograph method of Chaplygin, undertook a detailed investigation of the character of the flow near the line of transition from subsonic to supersonic velocities. From the results of Tricomi's investigation on the theory of differential equations of the mixed elliptic-hyperbolic type, Frankl introduced as one of the independent variables in place of the modulus of the velocity, a certain specially chosen function of this modulus. He thereby succeeded in explaining the character of the flow at the point of intersection of the transition line and the axis of symmetry (center of the nozzle) and in studying the behavior of the stream function in the neighborhood of this point by separating out the principal term having, together with its derivatives, the maximum value as compared with the corresponding corrections. This principal term is represented in Frankl's paper in the form of a linear combination of two hypergeometric functions. In order to find this linear combination, it is necessary to solve a number of boundary problems, which results in a complex analysis. In the investigation of the flow with which this paper is concerned, a second method is applied. This method is based on the transformation of the equations of motion to a form that may be called canonical for the system of differential equations of the mixed elliptic-hyperbolic type to which the system of equations of the motion of an ideal compressible fluid refers. By studying the behavior of the integrals of this system in the neighborhood of the parabolic line, the principal term of the solution is easily separated out in the form of a polynomial of the third degree. As a result, the computation of the transitional part of the nozzle is considerably simplified.

Falkovich, S. V.↗

NASA NTRS2005

Adaptive Failure Compensation for Aircraft Flight Control Using Engine Differentials: Regulation

The problem of using engine thrust differentials to compensate for rudder and aileron failures in aircraft flight control is addressed in this paper in a new framework. A nonlinear aircraft model that incorporates engine di erentials in the dynamic equations is employed and linearized to describe the aircraft s longitudinal and lateral motion. In this model two engine thrusts of an aircraft can be adjusted independently so as to provide the control flexibility for rudder or aileron failure compensation. A direct adaptive compensation scheme for asymptotic regulation is developed to handle uncertain actuator failures in the linearized system. A design condition is specified to characterize the system redundancy needed for failure compensation. The adaptive regulation control scheme is applied to the linearized model of a large transport aircraft in which the longitudinal and lateral motions are coupled as the result of using engine thrust differentials. Simulation results are presented to demonstrate the effectiveness of the adaptive compensation scheme.

Yu, Liu↗

NASA NTRS2006

3D Relativistic Magnetohydrodynamic Simulations of Magnetized Spine-Sheath Relativistic Jets

Numerical simulations of weakly magnetized and strongly magnetized relativistic jets embedded in a weakly magnetized and strongly magnetized stationary or weakly relativistic (v = c/2) sheath have been performed. A magnetic field parallel to the flow is used in these simulations performed by the new GRMHD numerical code RAISHIN used in its RMHD configuration. In the numerical simulations the Lorentz factor gamma = 2.5 jet is precessed to break the initial equilibrium configuration. In the simulations sound speeds are less than or equal to c/the square root of 3 in the weakly magnetized simulations and less than or equal to 0.56 c in the strongly magnetized simulations. The Alfven wave speed is less than or equal to 0.07 c in the weakly magnetized simulations and less than or equal to 0.56 c in the strongly magnetized simulations. The results of the numerical simulations are compared to theoretical predictions from a normal mode analysis of the linearized relativistic magnetohydrodynamic (RMHD) equations capable of describing a uniform axially magnetized cylindrical relativistic jet embedded in a uniform axially magnetized relativistically moving sheath. The theoretical dispersion relation allows investigation of effects associated with maximum possible sound speeds, Alfven wave speeds near light speed and relativistic sheath speeds. The prediction of increased stability of the weakly magnetized system resulting from c/2 sheath speeds and the stabilization of the strongly magnetized system resulting from c/2 sheath speeds is verified by the numerical simulation results.

Mizuno, Yosuke↗

NASA NTRS1983

Orbit structure of Hamiltonian systems arising from Lie transformation group actions

This paper associates the Riccati group and its group action on linear-quadratic optimal control problems to the action of a Lie transformation group on a set of Hamiltonian matrices. In this Lie theoretic setting results are presented concerning the associated orbit structure and the structure of the group itself. These results are of importance in understanding the solution structure of matrix Riccati differential equations, and thus also of importance in linear-quadratic optimal control.

Garzia, M. R.↗

DOE OSTI2021

Implicit Multirate GARK Methods

This article considers multirate generalized-structure additively partitioned Runge–Kutta methods for solving stiff systems of ordinary differential equations with multiple time scales. These methods treat different partitions of the system with different timesteps for a more targeted and efficient solution compared to monolithic single rate approaches. With implicit methods used across all partitions, methods must find a balance between stability and the cost of solving nonlinear equations for the stages. In order to characterize this important trade-off, we explore multirate coupling strategies, problems for assessing linear stability, and techniques to efficiently implement Newton iterations for stage equations. Unlike much of the existing multirate stability analysis which is limited in scope to particular methods, we present general statements on stability and describe fundamental limitations for certain types of multirate schemes. New implicit multirate methods up to fourth order are derived, and their accuracy and efficiency properties are verified with numerical tests.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Plasma Wave and Particle Dynamics During Interchange Events in the Jovian Magnetosphere Using Juno Observations

Abstract Interchange instability is known to drive fast radial transport of particles in Jupiter's inner magnetosphere. Magnetic flux tubes associated with the interchange instability often coincide with changes in particle distributions and plasma waves, but further investigations are required to understand their detailed characteristics. We analyze representative interchange events observed by Juno, which exhibit intriguing features of particle distributions and plasma waves, including Z‐mode and whistler‐mode waves. These events occurred at an equatorial radial distance of ∼9 Jovian radii on the nightside, with Z‐mode waves observed at mid‐latitude and whistler‐mode waves near the equator. We calculate the linear growth rate of whistler‐mode and Z‐mode waves based on the observed plasma parameters and electron distributions and find that both waves can be locally generated within the interchanged flux tube. Our findings are important for understanding particle transport and generation of plasma waves in the magnetospheres of Jupiter and other planetary systems.

Geology↗

DOE OSTI2022

Improving Common PV Module Temperature Models by Incorporating Radiative Losses to the Sky

PV module operating temperature is the second-most important factor influencing PV system yield–after irradiance–and a substantial contributor to uncertainty in energy system yield predictions. Models commonly used to predict operating temperature in system simulations are based on a simplified energy balance that lumps together different heat loss mechanisms–including radiation–and assumes an overall linear behavior. Radiative heat loss to the sky is usually substantial, but modeling it accurately requires additional information about down-welling long-wave radiation or sky temperature and increases the complexity of temperature model equations. In this work we show how radiative losses to the sky can be separated into two parts to improve the accuracy of modeling without additional complexity. We also predict and demonstrate the variation of these losses at different tilt angles and show that the effective view factor is reduced by the non- isotropic distribution of down-welling long-wave radiation. Finally, we demonstrate substantial reduction in bias (MBE) and scatter (RMSE) when the new radiative loss term is added to the Faiman model using one year of measurements at Sandia National Labs.

14 SOLAR ENERGY↗

NASA NTRS1982

Coupled flap-lag-torsional dynamics of hingeless rotor blades in forward flight

The aeroelastic stability and response problem of the coupled flap-lag-torsional dynamics of an isolated hingeless rotor blade in forward flight is considered. Linear, quasi-steady aerodynamics below stall is included. The spatial dependence of the partial differential nonliner equations of motion is discretized using a multimodal Galerkin method. The nonlinear time dependent equilibrium position (steady state response) about which the equations are linearized is obtained by solving a sequence of linear periodic response problems (quasilinearization). Results illustrating blade behavior in forward flight for both soft-in-plane and stiff-in-plane designs are presented. Quasilinearization provides a clear indication of the cases when nonlinear terms due to moderate deflections are important. The results indicate that nonlinearities affect system stability much more than system response. This implies that for obtaining blade vibrations and loads, approximate analytical models based on linearized formulations could be employed, whereas for determining blade stability a more accurate nonlinear analysis would be required.

Friedmann, P. P.↗

NASA NTRS1989

Efficient loads analyses of Shuttle-payloads using dynamic models with linear or nonlinear interfaces

An efficient method for the load analysis of Shuttle-payload systems with linear or nonlinear attachment interfaces is presented which allows the kinematics of the interface degrees of freedom at a given time to be evaluated without calculating the combined system modal representation of the Space Shuttle and its payload. For the case of a nonlinear dynamic model, an iterative procedure is employed to converge the nonlinear terms of the equations of motion to reliable values. Results are presented for a Shuttle abort landing event.

Spanos, P. D.↗

NASA NTRS1993

Minimum Hamiltonian Ascent Trajectory Evaluation (MASTRE) program (update to automatic flight trajectory design, performance prediction, and vehicle sizing for support of Shuttle and Shuttle derived vehicles) engineering manual

The Minimum Hamiltonian Ascent Trajectory Evaluation (MASTRE) program and its predecessors, the ROBOT and the RAGMOP programs, have had a long history of supporting MSFC in the simulation of space boosters for the purpose of performance evaluation. The ROBOT program was used in the simulation of the Saturn 1B and Saturn 5 vehicles in the 1960's and provided the first utilization of the minimum Hamiltonian (or min-H) methodology and the steepest ascent technique to solve the optimum trajectory problem. The advent of the Space Shuttle in the 1970's and its complex airplane design required a redesign of the trajectory simulation code since aerodynamic flight and controllability were required for proper simulation. The RAGMOP program was the first attempt to incorporate the complex equations of the Space Shuttle into an optimization tool by using an optimization method based on steepest ascent techniques (but without the min-H methodology). Development of the complex partial derivatives associated with the Space Shuttle configuration and using techniques from the RAGMOP program, the ROBOT program was redesigned to incorporate these additional complexities. This redesign created the MASTRE program, which was referred to as the Minimum Hamiltonian Ascent Shuttle TRajectory Evaluation program at that time. Unique to this program were first-stage (or booster) nonlinear aerodynamics, upper-stage linear aerodynamics, engine control via moment balance, liquid and solid thrust forces, variable liquid throttling to maintain constant acceleration limits, and a total upgrade of the equations used in the forward and backward integration segments of the program. This modification of the MASTRE code has been used to simulate the new space vehicles associated with the National Launch Systems (NLS). Although not as complicated as the Space Shuttle, the simulation and analysis of the NLS vehicles required additional modifications to the MASTRE program in the areas of providing additional flexibility in the use of the program, allowing additional optimization options, and providing special options for the NLS configuration.

Lyons, J. T.↗

NASA NTRS1998

Stability of Capillary Surfaces in Rectangular Containers: The Right Square Cylinder

The linearized governing equations for an ideal fluid are presented for numerical analysis for the stability of free capillary surfaces in rectangular containers against unfavorable disturbances (accelerations,i.e. Rayleigh-Taylor instability). The equations are solved for the case of the right square cylinder. The results are expressed graphically in term of a critical Bond number as a function of system contact angle. A critical wetting phenomena in the corners is shown to significantly alter the region of stability for such containers in contrast to simpler geometries such as the right circular cylinder or the infinite rectangular slot. Such computational results provide additional constraints for the design of fluids systems for space-based applications.

Weislogel, M. M.↗

NASA NTRS1961

A Theoretical Study of the Angular Motions of Spinning Bodies in Space

Euler's dynamic equations were linearized and solved analytically. Analytical expressions which relate angular motions to spin-rate and inertia distributions were obtained and found to be in good agreement with machine solutions of the nonlinear equations for the case of a rectangular-pulse pitching-moment disturbance. Consideration was given to the effects produced by having artificial damping in the system.

Suddath, Jerrold H.↗

NASA NTRS1997

Discretization and Preconditioning Algorithms for the Euler and Navier-Stokes Equations on Unstructured Meshes

Several stabilized discretization procedures for conservation law equations on triangulated domains will be considered. Specifically, numerical schemes based on upwind finite volume, fluctuation splitting, Galerkin least-squares, and space discontinuous Galerkin discretization will be considered in detail. A standard energy analysis for several of these methods will be given via entropy symmetrization. Next, we will present some relatively new theoretical results concerning congruence relationships for left or right symmetrized equations. These results suggest new variants of existing FV, DG, GLS and FS methods which are computationally more efficient while retaining the pleasant theoretical properties achieved by entropy symmetrization. In addition, the task of Jacobian linearization of these schemes for use in Newton's method is greatly simplified owing to exploitation of exact symmetries which exist in the system. These variants have been implemented in the "ELF" library for which example calculations will be shown. The FV, FS and DG schemes also permit discrete maximum principle analysis and enforcement which greatly adds to the robustness of the methods. Some prevalent limiting strategies will be reviewed. Next, we consider embedding these nonlinear space discretizations into exact and inexact Newton solvers which are preconditioned using a nonoverlapping (Schur complement) domain decomposition technique. Elements of nonoverlapping domain decomposition for elliptic problems will be reviewed followed by the present extension to hyperbolic and elliptic-hyperbolic problems. Other issues of practical relevance such the meshing of geometries, code implementation, turbulence modeling, global convergence, etc. will be addressed as needed.

Barth, Timothy↗

NASA NTRS1998

Discretization and Preconditioning Algorithms for the Euler and Navier-Stokes Equations on Unstructured Meshes

Several stabilized demoralization procedures for conservation law equations on triangulated domains will be considered. Specifically, numerical schemes based on upwind finite volume, fluctuation splitting, Galerkin least-squares, and space discontinuous Galerkin demoralization will be considered in detail. A standard energy analysis for several of these methods will be given via entropy symmetrization. Next, we will present some relatively new theoretical results concerning congruence relationships for left or right symmetrized equations. These results suggest new variants of existing FV, DG, GLS, and FS methods which are computationally more efficient while retaining the pleasant theoretical properties achieved by entropy symmetrization. In addition, the task of Jacobean linearization of these schemes for use in Newton's method is greatly simplified owing to exploitation of exact symmetries which exist in the system. The FV, FS and DG schemes also permit discrete maximum principle analysis and enforcement which greatly adds to the robustness of the methods. Discrete maximum principle theory will be presented for general finite volume approximations on unstructured meshes. Next, we consider embedding these nonlinear space discretizations into exact and inexact Newton solvers which are preconditioned using a nonoverlapping (Schur complement) domain decomposition technique. Elements of nonoverlapping domain decomposition for elliptic problems will be reviewed followed by the present extension to hyperbolic and elliptic-hyperbolic problems. Other issues of practical relevance such the meshing of geometries, code implementation, turbulence modeling, global convergence, etc, will. be addressed as needed.

Barth, Timothy J.↗