SEARCH · Engineering Papers
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.
Control-theory design package, called Optimal Regulator Algorithms for Control of Linear Systems (ORACLS), aids in design of controllers and optimal filters for systems modeled by linear, time invariant differential and difference equations. ORACLS is particularly attractive rigorous tool for dealing with multi-input and multi-output dynamic systems in both continuous and discrete forms.
A method is presented which transforms certain nonlinear differential equations of dynamics into linear equations by introducing an independent variable and utilizing the integrals of motion. As examples of special interest, the linearizations of unperturbed and perturbed Keplerian motions are discussed.
A set of equations known as Chandrasekhar equations arising in the linear quadratic optimal control problem is considered. In this paper, we consider the linear time-invariant system defined in Hilbert spaces involving unbounded input and output operators. For a general class of such systems, the Chandrasekhar equations are derived and the existence, uniqueness, and regularity of the results of their solutions established.
The advent of affordable computing, low-cost sensor hardware, and high-speed and reliable communications have spurred installation of ubiquitous sensors in complex engineered systems. However, ensuring reliable data quality remains a challenge. Exploitation of redundancy among sensor signals can help improving the precision of measured variables, detecting the presence of gross errors, and identifying faulty sensors. The cost of sensor ownership, maintenance efforts in particular, can still be cost-prohibitive however. Maximizing the ability to assess and control data quality while minimizing the cost of ownership thus requires a careful sensor placement. To solve this challenge, in this work we develop a generally applicable method to solve the multi-objective sensor placement problem in systems governed by linear and bilinear balance equations. Importantly, the method computes all Pareto-optimal sensor layouts with conventional computational resources and requires no information about the expected sensor quality.
The current state of the EUTERPE code is described with emphasis on the implemented models and their numerical implementation. The code solves the multi-species electromagnetic gyrokinetic equations in the full volume of a three-dimensional domain. Noise reduction of the particle-in-cell method is achieved by using a δf-method and Fourier filters. The field equations are discretized with B-splines and the resulting system of equations is solved iteratively. For linear simulations a phase-factor transformation is applied in order to strongly reduce the necessary grid resolution. Apart from the full gyrokinetic model, other numerically less expensive hybrid models are also implemented. They are mainly tailored for comparison with fluid theory and for studying the interaction of the bulk plasma with fast particles. The code is parallelized for CPUs by particle and domain decomposition. Good scalability up to several thousand nodes is demonstrated.
Runge-kutta integration to approximate system of nonlinear equations by series of linear equations
Perturbation of solutions of system of nonlinear differential equations under general linear boundary conditions
Stochastic approximation for identification of distributed parameter system solutions described by linear partial differential equations
Modal analysis computer program package for finding frequencies and mode shapes of any linear discrete system governed by generalized eigenvalue equation
Equations of motion simulating the landing approach case for the supersonic cruise vehicle were programmed and exercised using a fixed base simulation facility. Unaugmented and augmented system comparisons using this facility are provided. Refinements necessary for system performance improvement are made. The unaugmented longitudinal responses to elevator commands are slow and sluggish, requiring augmentation to increase the speed of the response. In the lateral directional case, the Dutch roll is highly underdamped and requires an augmentation system to increase this damping and provide satisfactory flying qualities. The status of this fixed base study is that the longitudinal equations, updated with recent wind tunnel data, were evaluated on the simulator and the system found to be satisfactory. The lateral axis equations are linearized and are not yet updated to large excursion capability; consequently, only limited, preliminary findings on this system are available. The basic results so far indicate augmentation systems are required to provide a satisfactory longitudinal system, and that additional study and evaluation of the lateral directional case are necessary before a more complete assessment can be made.
Numerical weather prediction (NWP) is an initial value problem for a system of nonlinear partial differential equations in which the initial values are known only incompletely and inaccurately. Data at initial time can be supplemented, however, by observations of the system distributed over a time interval preceding it. Estimation theory was successful in approaching such problems for models governed by systems of ordinary differential equations and of linear PDEs. Estimation-theoretic methods for NWP are developed. A model exhibiting many features of large scale atmospheric flow important in NWP is the one governed by the shallow fluid equations. The estimation problem for a linearized formulation of these equations is studied. A finite difference version of the equations is used as a forecast model to simulate the numerical models used in NWP.
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.
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.
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.
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.
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.
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.
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.