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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

The attenuation of sound in a randomly lined duct

The kinetic theory of wave propagation in random media is applied to assess the relative merits of the peripherally segmented, axially segmented and 'checkerboard' configured acoustic duct treatments proposed by Mani (1976). It is argued that, in a first approximation, a multisegmented liner involving many discrete and distinct sections may be modelled as a liner whose impedance varies randomly with position. This results in a relatively simple analytical description of the acoustic properties of the liner in terms of a system of linear, first order ordinary differential equations, whose solution gives the distribution with axial distance of the acoustic modal intensities. The theory is expected to be valid provided that the variation in the impedance from segment to segment is small, so that the conclusions obtained here are of a preliminary nature, and must be regarded as giving a general indication of features which should be incorporated into the design of efficient acoustic liners. The principal such conclusion is that a checkerboard liner offers the best means of increasing the attenuation over and above that of a uniformly lined duct. On the other hand, only minimal gains are to be derived from the use of strictly peripherally segmented liners.

Howe, M. S.↗
NASA NTRS1984

Nonlinear control of aircraft

Transformations of nonlinear systems were used to design automatic flight controllers for vertical and short takeoff aircraft. Under the assumption that a nonlinear system can be mapped to a controllable linear system, a method using partial differential equations was constructed to approximate transformations in cases where exact ones cannot be found. An application of the design theory to a rotorcraft, the UH-1H helicopter, was presented.

Hunt, L. R.↗
NASA NTRS1996

First-Order System Least Squares for Velocity-Vorticity-Pressure Form of the Stokes Equations, with Application to Linear Elasticity

In this paper, we study the least-squares method for the generalized Stokes equations (including linear elasticity) based on the velocity-vorticity-pressure formulation in d = 2 or 3 dimensions. The least squares functional is defined in terms of the sum of the L(exp 2)- and H(exp -1)-norms of the residual equations, which is weighted appropriately by by the Reynolds number. Our approach for establishing ellipticity of the functional does not use ADN theory, but is founded more on basic principles. We also analyze the case where the H(exp -1)-norm in the functional is replaced by a discrete functional to make the computation feasible. We show that the resulting algebraic equations can be uniformly preconditioned by well-known techniques.

Cai, Zhiqiang↗
NASA NTRS1999

The Vertical Linear Fractional Initialization Problem

This paper presents a solution to the initialization problem for a system of linear fractional-order differential equations. The scalar problem is considered first, and solutions are obtained both generally and for a specific initialization. Next the vector fractional order differential equation is considered. In this case, the solution is obtained in the form of matrix F-functions. Some control implications of the vector case are discussed. The suggested method of problem solution is shown via an example.

Lorenzo, Carl F.↗
NASA NTRS1994

Steady Secondary Momentum and Enthalpy Streaming in the Pulse Tube Refrigerator

Our study investigates the steady secondary streaming that occurs in the pulse tube refrigerator through application of the conservation laws and support by flow visualization experiments. A series expansion of the mass, momentum and energy conservation laws for an axisymmetric system yields a set of linear differential equations. The series solution is expanded in the small parameter epsilon = 1/Str where Str is the Strouhal number. The O(1) basic flow problem is solved for the case of a compressible ideal gas oscillating within a cylindrical tube. The phase and amplitude of the boundary conditions on axial velocity are treated as independent parameters. The basic solution for the gas domain is solved for the case of temperature and heat flux coupling to the tube wall. The temperature coupled solution shows how the Fourier number of the tube wall affects the gas temperature, phase angle, and time-averaged enthalpy flow. The basic solution is also examined for the time-averaged Reynolds stresses which are shown to be the driving force that produces steady secondary momentum streaming. Calculated solutions to the steady momentum streaming are shown to be in qualitative agreement with flow visualization experiments.

Lee, J. M.↗
NASA NTRS1976

Solution of linear systems by a singular perturbation technique

An approximate solution is obtained for a singularly perturbed system of initial valued, time invariant, linear differential equations with multiple boundary layers. Conditions are stated under which the approximate solution converges uniformly to the exact solution as the perturbation parameter tends to zero. The solution is obtained by the method of matched asymptotic expansions. Use of the results for obtaining approximate solutions of general linear systems is discussed. An example is considered to illustrate the method and it is shown that the formulas derived give a readily computed uniform approximation.

Ardema, M. D.↗
NASA NTRS1989

DYSCO: A software system for modeling general dynamic systems

The DYSCO program has been under development since 1979. It has been funded by Army and Air Force laboratories and by the Kaman Aerospace Corporation. It is presently available at a number of government and nongovernment installations. It has been used to analyze a very broad range of dynamics problems. A principle feature of the software design of DYSCO is the separation of the executive from the technology. The executive, which controls all the operations, is intelligent in the sense that it knows that its function is to assemble differential equations and to prepare them for solution. The technology library contains FORTRAN routines which perform standard functions, such as, computing the equation coefficients of an element (or component) given the local state at any time. The technology library also contains algorithms and procedures for solving the coupled system equations. The system was designed to allow easy additional of technology to the library. Any linear or nonlinear structural entity, control system, or set of ordinary differential equations may be simply coded and added to the library, as well as algorithms for time or frequency domain solution. The program is described with emphasis on its usefulness in easily modeling unusual concepts and configurations, performing analysis of damage, evaluating new algorithms, and simulating dynamic tests.

Berman, Alex↗
NASA NTRS1976

An iterative decoupling solution method for large scale Lyapunov equations

A great deal of attention has been given to the numerical solution of the Lyapunov equation. A useful classification of the variety of solution techniques are the groupings of direct, transformation, and iterative methods. The paper summarizes those methods that are at least partly favorable numerically, giving special attention to two criteria: exploitation of a general sparse system matrix structure and efficiency in resolving the governing linear matrix equation for different matrices. An iterative decoupling solution method is proposed as a promising approach for solving large-scale Lyapunov equation when the system matrix exhibits a general sparse structure. A Fortran computer program that realizes the iterative decoupling algorithm is also discussed.

Athay, T. M.↗
NASA NTRS1981

Determination of the potential in a synodic system

Determination of the potential field in a fixed (inertial) system may be accomplished by the solution of a homogeneous linear partial differential equation when a family of orbits of a body moving in the field is given. This partial differential equation was presented and thoroughly analyzed earlier. The present paper discusses the same problem in a rotating system where the centrifugal and Coriolis effects render the pertinent partial differential equation in general non-homogeneous and non-linear. A linear, though non-homogeneous, partial differential equation for the determination of the synodic potential is obtained only in the special case of iso-energetic families of orbits.

Szebehely, V.↗
NASA NTRS1996

Multiresolution and Explicit Methods for Vector Field Analysis and Visualization

We first report on our current progress in the area of explicit methods for tangent curve computation. The basic idea of this method is to decompose the domain into a collection of triangles (or tetrahedra) and assume linear variation of the vector field over each cell. With this assumption, the equations which define a tangent curve become a system of linear, constant coefficient ODE's which can be solved explicitly. There are five different representation of the solution depending on the eigenvalues of the Jacobian. The analysis of these five cases is somewhat similar to the phase plane analysis often associate with critical point classification within the context of topological methods, but it is not exactly the same. There are some critical differences. Moving from one cell to the next as a tangent curve is tracked, requires the computation of the exit point which is an intersection of the solution of the constant coefficient ODE and the edge of a triangle. There are two possible approaches to this root computation problem. We can express the tangent curve into parametric form and substitute into an implicit form for the edge or we can express the edge in parametric form and substitute in an implicit form of the tangent curve. Normally the solution of a system of ODE's is given in parametric form and so the first approach is the most accessible and straightforward. The second approach requires the 'implicitization' of these parametric curves. The implicitization of parametric curves can often be rather difficult, but in this case we have been successful and have been able to develop algorithms and subsequent computer programs for both approaches. We will give these details along with some comparisons in a forthcoming research paper on this topic.

Source record↗
NASA NTRS1992

On the continuous dependence with respect to sampling of the linear quadratic regulator problem for distributed parameter system

The convergence of solutions to the discrete- or sampled-time linear quadratic regulator problem and associated Riccati equation for infinite-dimensional systems to the solutions to the corresponding continuous time problem and equation, as the length of the sampling interval (the sampling rate) tends toward zero(infinity) is established. Both the finite-and infinite-time horizon problems are studied. In the finite-time horizon case, strong continuity of the operators that define the control system and performance index, together with a stability and consistency condition on the sampling scheme are required. For the infinite-time horizon problem, in addition, the sampled systems must be stabilizable and detectable, uniformly with respect to the sampling rate. Classes of systems for which this condition can be verified are discussed. Results of numerical studies involving the control of a heat/diffusion equation, a hereditary or delay system, and a flexible beam are presented and discussed.

Rosen, I. G.↗
NASA NTRS1973

Some conclusions regarding the aeroelastic stability of hingeless helicopter blades in hover and in forward flight

In this paper results and conclusions obtained from the study of the aeroelastic instability of hingeless helicopter blades are presented. First, the large amplitude coupled flap-lag equations of motion of a hingeless elastic helicopter blade are solved using an asymptotic expansion procedure in multiple time scales. Both hover and forward flight cases are considered. Stability boundaries and amplitudes of nonlinear response are obtained. From these, the importance of the nonlinear coupling and the effect of the periodic coefficients is determined. Next, using a system of linearized coupled flap-lag-pitch equations in hover, various divergence mechanisms for hingeless blades are shown. Finally, the flutter boundaries for coupled flap-lag-pitch are obtained. The effect of the torsional degree of freedom on the flap-lag type of instability is investigated. Similarly the effect of lag on the flap-pitch type of instability is considered. In addition, the effect of various blade parameters on the stability boundaries is shown.

Friedmann, P.↗
NASA NTRS1993

Theory of Metastable State Relaxation for Non-Critical Binary Systems with Non-Conserved Order Parameter

A new mathematical ansatz for a solution of the time-dependent Ginzburg-Landau non-linear partial differential equation is developed for non-critical systems such as non-critical binary solutions (solute + solvent) described by the non-conserved scalar order parameter. It is demonstrated that in such systems metastability initiates heterogeneous solute redistribution which results in formation of the non-equilibrium singly-periodic spatial solute structure. It is found how the time-dependent period of this structure evolves in time. In addition, the critical radius r(sub c) for solute embryo of the new solute rich phase together with the metastable state lifetime t(sub c) are determined analytically and analyzed.

Izmailov, Alexander↗
NASA NTRS1982

A computer package for the design and eigenproblem solution of damped linear multidegree of freedom systems

Systems described by the matrix differental equation are considered. An interactive design routine is presented for positive definite mass, damping, and stiffness matrices. Designing is accomplished by adjusting the mass, damping, and stiffness matrices to obtain a desired oscillation behavior. The algorithm also features interactively modifying the physical structure of the system, obtaining the matrix structure and a number of other system properties. In case of a general system, where the M, C, and K matrices lack any special properties, a routine for the eigenproblem solution of the system is developed. The latent roots are obtained by computing the characteristic polynomial of the system and solving for its roots. The above routines are prepared in FORTRAN IV and prove to be usable for the machines with low core memory.

Ahmadian, M.↗
DOE OSTI2022

A symbolic framework to obtain mid-fidelity models of flexible multibody systems with application to horizontal-axis wind turbines

Abstract. The article presents a symbolic framework (also called computer algebra program) that is used to obtain, in symbolic mathematical form, the linear and nonlinear equations of motion of a mid-fidelity multibody system including rigid and flexible bodies. Our approach is based on Kane's method and a nonlinear shape function representation for flexible bodies. The shape function approach does not represent the state of the art for flexible multibody dynamics but is an effective trade-off to obtain mid-fidelity models with few degrees of freedom, taking advantage of the separation of space and time. The method yields compact symbolic equations of motion with implicit account of the constraints. The general and automatic framework facilitates the creation and manipulation of models with various levels of complexity by adding or removing degrees of freedom. The symbolic treatment allows for analytical gradients and linearized equations of motion. The linear and nonlinear equations can be exported to Python code or dedicated software. There are multiple applications, such as time domain simulation, stability analyses, frequency domain analyses, advanced controller design, state observers, and digital twins. In this article, we describe the method we used to systematically generate the equations of motion of multibody systems and present the implementation of the framework using the Python package SymPy. We apply the framework to generate illustrative land-based and offshore wind turbine models. We compare our results with OpenFAST simulations and discuss the advantages and limitations of the method. The Python implementation is provided as an open-source project.

Branlard, Emmanuel (ORCID:0000000277506128)↗
DOE OSTI2024

Jacobian-based Model Diagnostics and Application to Equation Oriented Modeling of a Carbon Capture System

Equation-oriented (EO) modeling has the potential to enable the effective design and optimization of the operation of advanced energy systems. However, advanced modeling of energy systems results in a large number of variables and non-linear equations, and it can be difficult to search through these to identify the culprit(s) responsible for convergence issues. The Institute for the Design of Advanced Energy Systems Integrated Platform (IDAES-IP) contains a tool to identify poorly scaled constraints and variables by searching for rows and columns of the Jacobian matrix with small L2-norms so they can be rescaled. A further singular value decomposition can be per-formed to identify degenerate sets of equations and remaining scaling issues. This work presents an EO model of a flowsheet developed for post-combustion carbon capture using a monoethanolamine (MEA) solvent system as a case study. The IDAES diagnostics tools were successfully applied to this flowsheet to identify problems to improve model robustness and enable the optimization of process design and operating conditions of a carbon capture system.

Allan, Douglas↗
NASA NTRS1990

On the Decadal Modes of Oscillation of an Idealized Ocean-atmosphere System

Axially-symmetric, linear, free modes of global, primitive equation, ocean-atmosphere models are examined to see if they contain decadal (10 to 30 years) oscillation time scale modes. A two-layer ocean model and a two-level atmospheric model are linearized around axially-symmetric basic states containing mean meridional circulations in the ocean and the atmosphere. Uncoupled and coupled, axially-symmetric modes of oscillation of the ocean-atmosphere system are calculated. The main conclusion is that linearized, uncoupled and coupled, ocean-atmosphere systems can contain axially-symmetric, free modes of variability on decadal time scales. These results have important implications for externally-forced decadal climate variability.

Mehta, Vikram M.↗