Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear differential equations”

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 109 records · Page 6

Point to Point Control of the Hydrogen Mixer

We continue to build on the theoretical modelling of a liquid hydrogen (LH2) and gas hydrogen (GH2) mixer subsystem. The mixer described in this work is responsible for combining high pressure LH2 and GH2 to produce a hydrogen flow that meets certain thermodynamic properties. The desired properties are maintained by precise control of the LH2 and GH2 flows. The mixer is modelled as a general multi-flow lumped volume for single constituent fluids using density and internal energy as states. The set of nonlinear differential equations is modelled in the SIMULINK environment including a table look-up feature of the fluid thermodynamic properties. A small signal (linear) model is developed based on the nonlinear model and simulated as well. Pulse disturbances are introduced to the valve positions and the quality of the linear model is ascertained by comparing its behavior against the nonlinear model simulations. Valve control strategies that simulate an operator-in-the-loop scenario are then explored demonstrating the need for automatic feedback control. Finally, optimal single-output and multi-output Proportional/Integral controllers are designed based on the linear model and applied to the nonlinear model with excellent results to track simultaneous, constant setpoint changes in desired exit flow, exit temperature, and mixer pressure, as well as to reject unmeasureable but bounded additive step perturbations in the valve positions.

Barbieri, Enrique↗

Varieties of operator manipulation

A change in MACSYMA syntax is proposed to accommodate the operator manipulators necessary to implement direct and indirect methods for the solution of differential equations, calculus of finite differences, and the fractional calculus, as well as their modern counterparts. To illustrate the benefits and convenience of this syntax extension, an example is given to show how MACSYMA's pattern-matching capability can be used to implement a particular set of operator identities which can then be used to obtain exact solutions to nonlinear differential equations.

Doohovskoy, A.↗

Analytic theory of orbit contraction

The motion of a satellite in orbit, subject to atmospheric force and the motion of a reentry vehicle are governed by gravitational and aerodynamic forces. This suggests the derivation of a uniform set of equations applicable to both cases. For the case of satellite motion, by a proper transformation and by the method of averaging, a technique appropriate for long duration flight, the classical nonlinear differential equation describing the contraction of the major axis is derived. A rigorous analytic solution is used to integrate this equation with a high degree of accuracy, using Poincare's method of small parameters and Lagrange's expansion to explicitly express the major axis as a function of the eccentricity. The solution is uniformly valid for moderate and small eccentricities. For highly eccentric orbits, the asymptotic equation is derived directly from the general equation. Numerical solutions were generated to display the accuracy of the analytic theory.

Vinh, N. X.↗

Small Signal Point-to-Point Tracking of a Propellant Mixer

This paper addresses some theoretical modelling and control issues for a mixing chamber used in rocket engine testing at NASA Stennis Space Center. The mixer is responsible for combining high pressure LH2 and GH2 to produce a hydrogen flow that meets certain thermodynamic properties before it is fed into a test article. The desired properties are maintained by precise control of the LH2 and GH2 flows. The mixer is modelled as a general multi-flow lumped volume for single constituent fluids using density and internal energy as states. The set of nonlinear differential equations is modelled in the SIMULINK environment including a table look-up feature of the fluid thermodynamic properties. a small-signal (linear) model is developed based on the nonlinear model and simulated as well. Pulse disturbances are introduced to the valve positions and the quality of the linear model is ascertained by comparing its behavior against the nonlinear model simulations. Valve control strategies that simulate an operator-in-the-loop scenario are then explored demonstrating the need for automatic feedback control. Finally, classical optimal single-output and multi-output Proportional/Integral controllers are designed based on the linear model and applied to the nonlinear model with excellent results to track simultaneous, constant setpoint changes in desired exit flow, exit temperature, and mixer pressure, as well as to reject unmeasurable but bounded additive step perturbations in the valve positions.

Barbieri, Enrique↗

Ultimate boundedness stability and controllability of hereditary systems

By generalizing the Liapunov-Yoshizawa techniques, necessary and sufficient conditions are given for uniform boundedness and uniform ultimate boundedness of a rather general class of nonlinear differential equations of neutral type. Among the applications treated by the methods are the Lienard equation of neutral type and hereditary systems of Lurie type. The absolute stability of this later equation is also investigated. A certain existence result of a solution of a neutral functional differential inclusion with two point boundary values is applied to study the exact function space controllability of a nonlinear neutral functional differential control system. A geometric growth condition is used to characterize both the function space and Euclidean controllability of another nonlinear delay system which has a compact and convex control set. This yields conditions under which perturbed nonlinear delay controllable systems are controllable.

Chukwu, E. N.↗

Dynamics of Numerics and CFD

This lecture attempts to illustrate the basic ideas of how the recent advances in nonlinear dynamical systems theory (dynamics) can provide new insights into the understanding of numerical algorithms used in solving nonlinear differential equations (DEs). Examples will be given of the use of dynamics to explain unusual phenomena that occur in numerics. The inadequacy of the use of linearized analysis for the understanding of long time behavior of nonlinear problems will be illustrated, and the role of dynamics in studying the nonlinear stability, accuracy, convergence property and efficiency of using time- dependent approaches to obtaining steady-state numerical solutions in computational fluid dynamics (CFD) will briefly be explained.

Yee, H. C.↗

Theoretical hypervelocity ballistic limit for single or double plates using nonlinear modal analysis

The original research on the use on nonlinear vibration technique to solve for the hypervelocity ballistic limit for double plates is examined. Such structure is commonly found in typical Space Station design where the incoming space or man-made debris would be fragmented upon hitting the outer plate (shield) and the subsequent impact on the main wall would result in a much reduced damge of the space station or spacecraft. The existing few theoretical impact equations do not agree well with each other. The existing computer code bumper used at NASA-Johnson Space Center appears to predict an unconservative ballistic limit when compared with experimental data where the velocity ranges from 3 km/s to 8 km/s. Such unconservative prediction is unacceptable from a practical safe design point of view. The bumper code is based on Wilkinson's (1968) paper and his equations have not been improved nor modified even though they are viewed with suspicion due to lack of agreement with experiments. To make matters worse, there is not another theory which is better than Wilkinson's equation and the designers are forced to use purely empirical Nysmith (1969) or semiempirical equations developed by Cour-Palais in 1969. The Cour-Palais equations were later modified empirically in 1989. The purpose of the present investigation is to examine the many assumptions of the Wilkinson equation. An attempt is made to present design charts based on the modified-Wilkinson equation so that the designer can get a feel of the ranges of the parameters which are of interest and discard a huge range of parameters, thus, significantly reducing the number of test shots required. The analysis is based on a solution of the governing nonlinear differential equations for a plate, assuming axisymmetric behavior using polar coordinates.

Hui, David↗

Numerical derivation of linear perturbation model from nonlinear real-time vehicle simulations - A means of validation and increased credibility

A numerical technique has been developed which generates linear perturbation models from nonlinear real-time aircraft vehicle simulations at Langley Research Center. The technique is very general and can be applied to simulations of any system described by nonlinear differential equations. Linear models about various trim conditions have been obtained from simulations of the Boeing 737, Northrop F5-E, Cessna 180, Sikorsky CH-54, and NASA/Army Rotor Systems Research Aircraft (RSRA). One use of the derived linear models is as a validation and verification tool for the nonlinear simulation.

Dieudonne, J. E.↗

Time scales for molecule formation by ion-molecule reactions

Analytical solutions are obtained for nonlinear differential equations governing the time-dependence of molecular abundances in interstellar clouds. Three gas-phase reaction schemes are considered separately for the regions where each dominates. The particular case of CO, and closely related members of the Oh and CH families of molecules, is studied for given values of temperature, density, and the radiation field. Nonlinear effects and couplings with particular ions are found to be important. The time scales for CO formation range from 100,000 to a few million years, depending on the chemistry and regime. The time required for essentially complete conversion of C(+) to CO in the region where the H3(+) chemistry dominates is several million years. Because this time is longer than or comparable to dynamical time scales for dense interstellar clouds, steady-state abundances may not be observed in such clouds.

Langer, W. D.↗

Analysis and design of an ultrahigh temperature hydrogen-fueled MHD generator

A coupled gas dynamics/radiative heat transfer analysis of partially ionized hydrogen, in local thermodynamic equilibrium, flowing through an ultrahigh temperature (10,000-20,000 K) magnetohydrodynamic (MHD) generator is performed. Gas dynamics are modeled by a set of quasi-one-dimensional, nonlinear differential equations which account for friction, convective and radiative heat transfer, and the interaction between the ionized gas and applied magnetic field. Radiative heat transfer is modeled using nongray, absorbing-emitting 2D and 3D P-1 approximations which permit an arbitrary variation of the spectral absorption coefficient with frequency. Gas dynamics and radiative heat transfer are coupled through the energy equation and through the temperature- and density-dependent absorption coefficient. The resulting nonlinear elliptic problem is solved by iterative methods. Design of such MHD generators as onboard, open-cycle, electric power supplies for a particular advanced airbreathing propulsion concept produced an efficient and compact 128-MWe generator characterized by an extraction ratio of 35.5 percent, a power density of 10,500 MWe/cu m, and a specific (extracted) energy of 324 MJe/kg of hydrogen. The maximum wall heat flux and total wall heat load were 453 MW/sq m and 62 MW, respectively.

Moder, Jeffrey P.↗

Computational alternatives to obtain time optimal jet engine control

Two computational methods to determine an open loop time optimal control sequence for a simple single spool turbojet engine are described by a set of nonlinear differential equations. Both methods are modifications of widely accepted algorithms which can solve fixed time unconstrained optimal control problems with a free right end. Constrained problems to be considered have fixed right ends and free time. Dynamic programming is defined on a standard problem and it yields a successive approximation solution to the time optimal problem of interest. A feedback control law is obtained and it is then used to determine the corresponding open loop control sequence. The Fletcher-Reeves conjugate gradient method has been selected for adaptation to solve a nonlinear optimal control problem with state variable and control constraints.

Basso, R. J.↗

A vectorized algorithm for 3D dynamics of a tethered satellite

Equations of motion characterizing the three dimensional motion of a tethered satellite during the retrieval phase are studied. The mathematical model involves an arbitrary number of point masses connected by weightless cords. Motion occurs in a gravity gradient field. The formulation presented accounts for general functions describing support point motion, rate of tether retrieval, and arbitrary forces applied to the point masses. The matrix oriented program language MATLAB is used to produce an efficient vectorized formulation for computing natural frequencies and mode shapes for small oscillations about the static equilibrium configuration; and for integrating the nonlinear differential equations governing large amplitude motions. An example of time response pertaining to the skip rope effect is investigated.

Wilson, Howard B.↗

The General Equations of Equilibrium of Rotationally Symmetric Membranes and some Static Solutions for Uniform Centrifugal Loading

The general equations of equilibrium of a rotationally symmetrical body undergoing large arbitrary deflections but small strains are derived by means of the principle of virtual work. All deflections are referred to a reference surface which, in general, may be taken to coincide with any intermediate state of the membrane. Static solutions of the general equations for steadily rotating membranes are considered for the cases in which the final or the initial configuration of the membrane is specified. For problems in which the final configuration is specified, conditions are given under which such a configuration is obtainable from a stress-free initial state. Such conditions being satisfied, it is shown that the stresses and radial displacements are determinate while axial displacements may be obtained by quadrature. For problems in which the initial configuration is given, the determination of the stresses and deflections is reduced to the solution of a second-order nonlinear differential equation for a stress function. An approximate method of solution based on boundary-layer considerations is presented for the case of a spinning sphere of constant thickness. Specific examples are given for both classes of problems.

Simmonds, James G.↗

Nonlinear Dynamics of Cantilever-Sample Interactions in Atomic Force Microscopy

The interaction of the cantilever tip of an atomic force microscope (AFM) with the sample surface is obtained by treating the cantilever and sample as independent systems coupled by a nonlinear force acting between the cantilever tip and a volume element of the sample surface. The volume element is subjected to a restoring force from the remainder of the sample that provides dynamical equilibrium for the combined systems. The model accounts for the positions on the cantilever of the cantilever tip, laser probe, and excitation force (if any) via a basis set of set of orthogonal functions that may be generalized to account for arbitrary cantilever shapes. The basis set is extended to include nonlinear cantilever modes. The model leads to a pair of coupled nonlinear differential equations that are solved analytically using a matrix iteration procedure. The effects of oscillatory excitation forces applied either to the cantilever or to the sample surface (or to both) are obtained from the solution set and applied to the to the assessment of phase and amplitude signals generated by various acoustic-atomic force microscope (A-AFM) modalities. The influence of bistable cantilever modes of on AFM signal generation is discussed. The effects on the cantilever-sample surface dynamics of subsurface features embedded in the sample that are perturbed by surface-generated oscillatory excitation forces and carried to the cantilever via wave propagation are accounted by the Bolef-Miller propagating wave model. Expressions pertaining to signal generation and image contrast in A-AFM are obtained and applied to amplitude modulation (intermittent contact) atomic force microscopy and resonant difference-frequency atomic force ultrasonic microscopy (RDF-AFUM). The influence of phase accumulation in A-AFM on image contrast is discussed, as is the effect of hard contact and maximum nonlinearity regimes of A-AFM operation.

Cantrell, John H.↗