Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differential algebra”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Response of joint dominated space structures

An efficient linearization method is presented for calculating the transient response of nonlinear systems due to initial disturbances. The method is an extension of the describing function approach in which the steady state response of the system is calculated by representing the nonlinear element, typically joints in the case of space structures, by impedances which are functions of the amplitude of response. Thus, the problem of solving the differential equation for the steady state response becomes one of solving a set of nonlinear algebraic equations involving the steady state amplitudes and phases of the system. It is shown that for the transient case the steady state impedances can be averaged over the range of responses in order to provide equivalent values of stiffness and damping that, for a given set of initial displacements, may be treated as being constant for purposes of calculating system response. Single degree of freedom system are used to demonstrate the method and to develop an approach for optimizing the joint's characteristics so as to minimize transient response times. The use of this method for response estimation and optimization in multiple degree of freedom systems is investigated.

Source record↗
NASA NTRS1988

The ATOMFT integrator - Using Taylor series to solve ordinary differential equations

This paper discusses the application of ATOMFT, an integration package based on Taylor series solution with a sophisticated user interface. ATOMFT has the capabilities to allow the implementation of user defined functions and the solution of stiff and algebraic equations. Detailed examples, including the solutions to several astrodynamics problems, are presented. Comparisons with its predecessor ATOMCC and other modern integrators indicate that ATOMFT is a fast, accurate, and easy method to use to solve many differential equation problems.

Berryman, Kenneth W.↗
NASA NTRS1975

Some properties of extremal solutions of zero-sum differential games

Necessary conditions for infinite-time linear-quadratic two-player differential games with perfect state information are examined in the light of Willems' work on the algebraic Riccati equation. Necessary conditions for the nonzero-sum game do not always reduce precisely to those of the zero-sum game, when the weighting parameters of the cost index are taken as those of a zero-sum game. Formulation of zero-sum game problems in extended L2-spaces reveals the existence of nonstabilizing but infimal controls, and explicit frequency-domain conditions for existence of solutions and stabilizing solutions are established through the algebraic Riccati equation (ARE).

Johnson, T. L.↗
NASA NTRS1975

Fixed-range optimum trajectories for short-haul aircraft

An algorithm, based on the energy-state method, is derived for calculating optimum trajectories with a range constraint. The basis of the algorithm is the assumption that optimum trajectories consist of, at most, three segments: an increasing energy segment (climb); a constant energy segment (cruise); and a decreasing energy segment (descent). This assumption allows energy to be used as the independent variable in the increasing and decreasing energy segments, thereby eliminating the integration of a separate adjoint differential equation and simplifying the calculus of variations problem to one requiring only pointwise extremization of algebraic functions. The algorithm is used to compute minimum fuel, minimum time, and minimum direct-operating-cost trajectories, with range as a parameter, for an in-service CTOL aircraft and for an advanced STOL aircraft. For the CTOL aircraft and the minimum-fuel performance function, the optimum controls, consisting of air-speed and engine power setting, are continuous functions of the energy in both climb and descent as well as near the maximum or cruise energy. This is also true for the STOL aircraft except in the descent where at one energy level a nearly constant energy dive segment occurs, yielding a discontinuity in the airspeed at that energy. The reason for this segment appears to be the relatively high fuel flow at idle power of the engines used by this STOL aircraft. Use of a simplified trajectory which eliminates the dive increases the fuel consumption of the total descent trajectory by about 10 percent and the time to fly the descent by about 19 percent compared to the optimum.

Erzberger, H.↗
NASA NTRS2007

Analytic Guidance for the First Entry in a Skip Atmospheric Entry

This paper presents an analytic method to generate a reference drag trajectory for the first entry portion of a skip atmospheric entry. The drag reference, expressed as a polynomial function of the velocity, will meet the conditions necessary to fit the requirements of the complete entry phase. The generic method proposed to generate the drag reference profile is further simplified by thinking of the drag and the velocity as density and cumulative distribution functions respectively. With this notion it will be shown that the reference drag profile can be obtained by solving a linear algebraic system of equations. The resulting drag profile is flown using the feedback linearization method of differential geometric control as guidance law with the error dynamics of a second order homogeneous equation in the form of a damped oscillator. This approach was first proposed as a revisited version of the Space Shuttle Orbiter entry guidance. However, this paper will show that it can be used to fly the first entry in a skip entry trajectory. In doing so, the gains in the error dynamics will be changed at a certain point along the trajectory to improve the tracking performance.

Garcia-Llama, Eduardo↗
NASA NTRS1993

A critical analysis of the accuracy of several numerical techniques for combustion kinetic rate equations

A detailed analysis of the accuracy of several techniques recently developed for integrating stiff ordinary differential equations is presented. The techniques include two general-purpose codes EPISODE and LSODE developed for an arbitrary system of ordinary differential equations, and three specialized codes CHEMEQ, CREK1D, and GCKP4 developed specifically to solve chemical kinetic rate equations. The accuracy study is made by application of these codes to two practical combustion kinetics problems. Both problems describe adiabatic, homogeneous, gas-phase chemical reactions at constant pressure, and include all three combustion regimes: induction, heat release, and equilibration. To illustrate the error variation in the different combustion regimes the species are divided into three types (reactants, intermediates, and products), and error versus time plots are presented for each species type and the temperature. These plots show that CHEMEQ is the most accurate code during induction and early heat release. During late heat release and equilibration, however, the other codes are more accurate. A single global quantity, a mean integrated root-mean-square error, that measures the average error incurred in solving the complete problem is used to compare the accuracy of the codes. Among the codes examined, LSODE is the most accurate for solving chemical kinetics problems. It is also the most efficient code, in the sense that it requires the least computational work to attain a specified accuracy level. An important finding is that use of the algebraic enthalpy conservation equation to compute the temperature can be more accurate and efficient than integrating the temperature differential equation.

Radhadrishnan, Krishnan↗
NASA NTRS1973

High speed flow past wings

The analytical solution to the transonic small perturbation equation which describes steady compressible flow past finite wings at subsonic speeds can be expressed as a nonlinear integral equation with the perturbation velocity potential as the unknown function. This known formulation is substituted by a system of nonlinear algebraic equations to which various methods are applicable for its solution. Due to the presence of mathematical discontinuities in the flow solutions, however, a main computational difficulty was to ensure uniqueness of the solutions when local velocities on the wing exceeded the speed of sound. For continuous solutions this was achieved by embedding the algebraic system in an one-parameter operator homotopy in order to apply the method of parametric differentiation. The solution to the initial system of equations appears then as a solution to a Cauchy problem where the initial condition is related to the accompanying incompressible flow solution. In using this technique, however, a continuous dependence of the solution development on the initial data is lost when the solution reaches the minimum bifurcation point. A steepest descent iteration technique was therefore, added to the computational scheme for the calculation of discontinuous flow solutions. Results for purely subsonic flows and supersonic flows with and without compression shocks are given and compared with other available theoretical solutions.

Norstrud, H.↗
NASA NTRS1982

Improved Numerical Differencing Analyzer

SINDA, Systems Improved Numerical Differencing Analyzer, solves differential and algebric equations representing physical systems. SINDA solves numerically almost any set of ordinary differential equations that represent transient behavior of a lumped-parameter system or any set of nonlinear algebraic equations that represents the steady state conditions of a physical system.

Skladany, J. T.↗
NASA NTRS1989

Numerical computation of shock wave-turbulent boundary layer interaction in transonic flow over an axisymmetric curved hill

A control-volume based finite difference computation of a turbulent transonic flow over an axisymmetric curved hill is presented. The numerical method is based on the SIMPLE algorithm, and hence the conservation of mass equation is replaced by a pressure correction equation for compressible flows. The turbulence is described by a k-epsilon turbulence model supplemented by a near-wall turbulence model. In the method, the dissipation rate in the region very close to the wall is obtained from an algebraic equation and that for the rest of the flow domain is obtained by solving a partial differential equation for the dissipation rate. The other flow equations are integrated up to the wall. It is shown that the present turbulence model yields the correct location of the compression shock. The other computational results are also in good agreement with experimental data.

Kim, S.-W.↗
NASA NTRS1989

State-of-the-art surveys on computational mechanics

Topics considered include advances in finite difference techniques for computational fluid dynamics, the spectral element methods for the incompressible Navier-Stokes equations, a review of recent developments in time integration, and advances and trends in element-by-element techniques. Also examined are the algebraic multigrid methods applied to problems in computational structural mechanics, grid generation for the solution of partial differential equations, advances in adaptive improvements, and new computing systems and their impact on computational mechanics.

Noor, Ahmed K.↗
NASA NTRS2008

Simulink Model of the Ares I Upper Stage Main Propulsion System

A numerical model of the Ares I upper stage main propulsion system is formulated based on first principles. Equation's are written as non-linear ordinary differential equations. The GASP fortran code is used to compute thermophysical properties of the working fluids. Complicated algebraic constraints are numerically solved. The model is implemented in Simulink and provides a rudimentary simulation of the time history of important pressures and temperatures during re-pressurization, boost and upper stage firing. The model is validated against an existing reliable code, and typical results are shown.

Burchett, Bradley T.↗
NASA NTRS2010

Analytic Development of a Reference Profile for the First Entry in a Skip Atmospheric Entry

This note shows that a feasible reference drag profile for the first entry portion of a skip entry can be generated as a polynomial expression of the velocity. The coefficients of that polynomial are found through the resolution of a system composed of m + 1 equations, where m is the degree of the drag polynomial. It has been shown that a minimum of five equations (m = 4) are required to establish the range and the initial and final conditions on velocity and flight path angle. It has been shown that at least one constraint on the trajectory can be imposed through the addition of one extra equation in the system, which must be accompanied by the increase in the degree of the drag polynomial. In order to simplify the resolution of the system of equations, the drag was considered as being a probability density function of the velocity, with the velocity as a distribution function of the drag. Combining this notion with the introduction of empirically derived constants, it has been shown that the system of equations required to generate the drag profile can be successfully reduced to a system of linear algebraic equations. For completeness, the resulting drag profiles have been flown using the feedback linearization method of differential geometric control as a guidance law with the error dynamics of a second order homogeneous equation in the form of a damped oscillator. Satisfactory results were achieved when the gains in the error dynamics were changed at a certain point along the trajectory that is dependent on the velocity and the curvature of the drag as a function of the velocity. Future work should study the capacity to update the drag profile in flight when dispersions are introduced. Also, future studies should attempt to link the first entry, as presented and controlled in this note, with a more standard control concept for the second entry, such as the Apollo entry guidance, to try to assess the overall skip entry performance. A guidance law that includes an integral feedback term, as is the case in the actual Space Shuttle entry guidance and as is proposed in Ref 29, could be tried in future studies to assess whether its use results in an improvement of the tracking performance, and to evaluate the design needs when determining the control gains.

Garcia-Llama, Eduardo↗
NASA NTRS2004

Rapid Airplane Parametric Input Design(RAPID)

An efficient methodology is presented for defining a class of airplane configurations. Inclusive in this definition are surface grids, volume grids, and grid sensitivity. A small set of design parameters and grid control parameters govern the process. The general airplane configuration has wing, fuselage, vertical tail, horizontal tail, and canard components. The wing, tail, and canard components are manifested by solving a fourth-order partial differential equation subject to Dirichlet and Neumann boundary conditions. The design variables are incorporated into the boundary conditions, and the solution is expressed as a Fourier series. The fuselage has circular cross section, and the radius is an algebraic function of four design parameters and an independent computational variable. Volume grids are obtained through an application of the Control Point Form method. Grid sensitivity is obtained by applying the automatic differentiation precompiler ADIFOR to software for the grid generation. The computed surface grids, volume grids, and sensitivity derivatives are suitable for a wide range of Computational Fluid Dynamics simulation and configuration optimizations.

Smith, Robert E.↗
NASA NTRS1973

Prediction of free shear flows: A comparison of the performance of six turbulence models

The performance is evaluated of three distinct classes of turbulence model. These classes are: (1) Turbulent-viscosity models in which the length scale of turbulence is found by way of algebraic formulas, (2) turbulent-viscosity models in which the length scale of turbulence is found from a partial differential equation of transport, and (3) models in which the shear stress itself is the dependent variable of a partial differential conservation equation. Two models were examined in each class; thus, six different models were tested. A complete mathematical statement of these models is provided and a brief commentary on the models is included.

Launder, B. E.↗
NASA NTRS1990

Trees, bialgebras and intrinsic numerical algorithms

Preliminary work about intrinsic numerical integrators evolving on groups is described. Fix a finite dimensional Lie group G; let g denote its Lie algebra, and let Y(sub 1),...,Y(sub N) denote a basis of g. A class of numerical algorithms is presented that approximate solutions to differential equations evolving on G of the form: dot-x(t) = F(x(t)), x(0) = p is an element of G. The algorithms depend upon constants c(sub i) and c(sub ij), for i = 1,...,k and j is less than i. The algorithms have the property that if the algorithm starts on the group, then it remains on the group. In addition, they also have the property that if G is the abelian group R(N), then the algorithm becomes the classical Runge-Kutta algorithm. The Cayley algebra generated by labeled, ordered trees is used to generate the equations that the coefficients c(sub i) and c(sub ij) must satisfy in order for the algorithm to yield an rth order numerical integrator and to analyze the resulting algorithms.

Crouch, Peter↗
NASA NTRS1972

A feedback control model for network flow with multiple pure time delays

A control model describing a network flow hindered by multiple pure time (or transport) delays is formulated. Feedbacks connect each desired output with a single control sector situated at the origin. The dynamic formulation invokes the use of differential difference equations. This causes the characteristic equation of the model to consist of transcendental functions instead of a common algebraic polynomial. A general graphical criterion is developed to evaluate the stability of such a problem. A digital computer simulation confirms the validity of such criterion. An optimal decision making process with multiple delays is presented.

Press, J.↗
NASA NTRS1996

Chebyshev polynomials in the spectral Tau method and applications to Eigenvalue problems

Chebyshev Spectral methods have received much attention recently as a technique for the rapid solution of ordinary differential equations. This technique also works well for solving linear eigenvalue problems. Specific detail is given to the properties and algebra of chebyshev polynomials; the use of chebyshev polynomials in spectral methods; and the recurrence relationships that are developed. These formula and equations are then applied to several examples which are worked out in detail. The appendix contains an example FORTRAN program used in solving an eigenvalue problem.

Johnson, Duane↗
NASA NTRS2011

Stellar Mixing: I. Formalism

In this paper we use the Reynolds stress models (RSM) to derive algebraic expressions for the following variables: a) heat fluxes; b) J.l fluxes; and c) momentum fluxes. These relations, which are fully 3D, include: 1) stable and unstable stratification, represented by the Brunt-Vaislila frequency, N(exp 2) =-g/H(sub p_(del - del(sub ad))(1 - RI(sub mu)); 2) double diffusion, salt-fingers, and semi-convection, represented by the density ratio R(sub mu) = del(sub mu)/(del - del(sub ad)); 3) shear (differential rotation), represented by the mean squared shear Sigma(exp 2) or by the Richardson number, Ri =N(exp 2)Sigma(exp -2); 4) radiative losses represented by a Peclet number, Pe; 5) a complete analytical solution of the ID version of the model. In general, the model requires the solution of two differential equations for the eddy kinetic energy K and its rate of dissipation, epsilon. In the local and stationary cases, when production equals dissipation, the model equations are all algebraic.

Canuto, V .M.↗