SEARCH · Search NASA

Results for “Ordinary Differential Equation”

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 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS1982

Spline approximations for nonlinear hereditary control systems

A sline-based approximation scheme is discussed for optimal control problems governed by nonlinear nonautonomous delay differential equations. The approximating framework reduces the original control problem to a sequence of optimization problems governed by ordinary differential equations. Convergence proofs, which appeal directly to dissipative-type estimates for the underlying nonlinear operator, are given and numerical findings are summarized.

Daniel, P. L.↗

NASA NTRS1977

On the stability of Galerkin methods for initial-boundary value problems for hyperbolic systems

The stability of approximating the solution of mixed initial-boundary value problems for hyperbolic systems by semidiscrete Galerkin methods is studied. It is shown that a particular straightforward Galerkin method yields an unstable approximation, and that this numerical instability is caused by an improper treatment of the boundary. Stable schemes are then presented, one of which differs from the unstable scheme only insofar as the treatment of the boundary is concerned. These stable schemes make use of a particular matrix which symmetrizes the differential system. It is therefore shown that the use of this matrix is crucial to the stability of the computations as well as for obtaining a priori bounds on the energy of the continuous system. This symmetrizing matrix is also related to the diagonalizing matrix for the system of hyperbolic equations and to the Liapunov matrix for the system of ordinary differential equations resulting from the application of Galerkin's method.

Gunzburger, M. D.↗

NASA NTRS1990

An equation for the evolution of solar and stellar flare loops

An ordinary differential equation describing the evolution of a coronal loop subjected to a spatially uniform but time-varying heating rate is discussed. It is assumed that the duration of heating is long compared to the sound transit time through the loop, which is assumed to have uniform cross section area. The form of the equation changes as the loop evolves through three states: 'strong evaporation', 'scaling law behavior', and 'strong condensation'. Solutions to the equation may be used to compute the time dependence of the average coronal temperature and emission measure for an assumed temporal variation of the flare heating rate. The results computed from the model agree reasonably well with recent published numerical simulations and may be obtained with far less computational effort. The model is then used to study the May 21, 1980, solar flare observed by SMM and the giant April 12, 1985, flare observed on the star AD Leo.

Fisher, George H.↗

NASA NTRS1977

Initial values for the integration scheme to compute the eigenvalues for propagation in ducts

A scheme for the calculation of eigenvalues in the problem of acoustic propagation in a two-dimensional duct is described. The computation method involves changing the coupled transcendental nonlinear algebraic equations into an initial value problem involving a nonlinear ordinary differential equation. The simplest approach is to use as initial values the hardwall eigenvalues and to integrate away from these values as the admittance varies from zero to its actual value with a linear variation. The approach leads to a powerful root finding routine capable of computing the transverse and axial wave numbers for two-dimensional ducts for any frequency, lining, admittance and Mach number without requiring initial guesses or starting points.

Eversman, W.↗

NASA NTRS1993

Inclusion of transverse shear deformation in the exact buckling and vibration analysis of composite plate assemblies

The problem considered is the development of the necessary plate stiffnesses for use in the general purpose program VICONOPT for buckling and vibration of composite plate assemblies. The required stiffnesses include the effects of transverse shear deformation and are for sinusoidal response along the plate length as required in VICONOPT. The method is based on the exact solution of the plate differential equations for a composite laminate having fully populated A, B, and D stiffness matrices which leads to an ordinary differential equation of tenth order.

Anderson, Melvin S.↗

NASA NTRS1993

A numerical procedure for analysis of finite rate reacting flows

Combustion processes in rocket propulsion systems are characterized by the existence of multiple, vastly differing time and length scales, as well as flow-speeds at wide variation of Mach numbers. The chemical kinetics processes in the highly active reaction zone are characterized by much smaller scales compared to fluid convective and diffusive time scales. An operator splitting procedure for transient finite rate chemistry problems has been developed using a pressure based method, which can be applied to all speed flows without difficulties. The splitting of chemical kinetics terms formed the fluid-mechanical terms of the species equation ameliorated the difficulties associated with the disparate time scales and stiffness in the set of equations which describes highly exothermic combustion. A combined efficient ordinary differential equations (ODE) solver was used to integrate the effective chemical source terms over the residence time at each grid cell. One and two dimensional reacting flow situations were carried out to demonstrate and verify the current procedure. Different chemical kinetics with different degrees of nonlinearity have also been incorporated to test the robustness and generality of the proposed method.

Shang, H. M.↗

NASA NTRS1977

Symbolic computation of recurrence equations for the Chebyshev series solution of linear ODE's

If a linear ordinary differential equation with polynomial coefficients is converted into integrated form then the formal substitution of a Chebyshev series leads to recurrence equations defining the Chebyshev coefficients of the solution function. An explicit formula is presented for the polynomial coefficients of the integrated form in terms of the polynomial coefficients of the differential form. The symmetries arising from multiplication and integration of Chebyshev polynomials are exploited in deriving a general recurrence equation from which can be derived all of the linear equations defining the Chebyshev coefficients. Procedures for deriving the general recurrence equation are specified in a precise algorithmic notation suitable for translation into any of the languages for symbolic computation. The method is algebraic and it can therefore be applied to differential equations containing indeterminates.

Geddes, K. O.↗

NASA NTRS1984

A comparison of the efficiency of numerical methods for integrating chemical kinetic rate equations

The efficiency of several algorithms used for numerical integration of stiff ordinary differential equations was compared. The methods examined included two general purpose codes EPISODE and LSODE and three codes (CHEMEQ, CREK1D and GCKP84) developed specifically to integrate chemical kinetic rate equations. The codes were applied to two test problems drawn from combustion kinetics. The comparisons show that LSODE is the fastest code available for the integration of combustion kinetic rate equations. It is shown that an iterative solution of the algebraic energy conservation equation to compute the temperature can be more efficient then evaluating the temperature by integrating its time-derivative.

Radhakrishnan, K.↗

NASA NTRS1995

On the Dynamics of Some Discretizations of Convection-Diffusion Equations

Numerical discretizations of differential equations which model physical processes can possess dynamics quite different from that of the equations themselves. Recently the emphasis has been on the the dynamics of numerical discretizations for Ordinary Differential Equations (ODEs). For Partial Differential Equations (PDEs) using a method of lines approach the situation is more complex. First, the spatial discretisation may introduce dynamics not present in the original equations; second, the solution of the resulting system of ODEs is open to the modified dynamics of the ODE solver used. These two effects may interact in a complex manner. In this talk we present some results of our recent work on the dynamics of discretizations of convection-diffusion equations, including those produced using Total Variation Diminishing (TVD) schemes and adaptive grid techniques. A more general overview of the area may be found on our accompanying poster presentation.

Sweby, Peter K.↗

NASA NTRS1973

Aerodynamic noise emission from turbulent shear layers.

The Phillips (1960) convected wave equation is employed in this paper to study aerodynamic noise emission processes in subsonic and supersonic shear layers. The wave equation in three spatial dimensions is first reduced to an ordinary differential equation by Fourier transformation and then solved via the WKBJ method. Three typical solutions are required for discussions in this paper. The current results are different from the classical conclusions. The effects of refraction, convection, Mach-number dependence and temperature dependence of turbulent noise emission are analyzed in the light of solutions to the Phillips equation.

Pao, S. P.↗

NASA NTRS1978

Application of the method of integral relations to boundary layer flows over blunt bodies

Calculations of boundary layer flows past blunt bodies at angles of incidence are presented. Using the method of integral relations together with the method of lines, the full three-dimensional boundary layer equations are reduced to a system of first-order ordinary differential equations. The streamwise shear stress function and the cross-flow velocity component are represented as suitable functions of the streamwise velocity component. The role of the zone of dependence is automatically satisfied by the choice of differencing in the method of lines. Solutions correct to the second order are obtained in the positive shear region for flow over an ellipsoid at 30-deg incidence. The results are compared with corresponding finite difference solutions.

Modarress, D.↗

NASA NTRS1980

The primer vector in linear, relative-motion equations

Primer vector theory is used in analyzing a set of linear, relative-motion equations - the Clohessy-Wiltshire equations - to determine the criteria and necessary conditions for an optimal, N-impulse trajectory. Since the state vector for these equations is defined in terms of a linear system of ordinary differential equations, all fundamental relations defining the solution of the state and costate equations, and the necessary conditions for optimality, can be expressed in terms of elementary functions. The analysis develops the analytical criteria for improving a solution by (1) moving any dependent or independent variable in the initial and/or final orbit, and (2) adding intermediate impulses. If these criteria are violated, the theory establishes a sufficient number of analytical equations. The subsequent satisfaction of these equations will result in the optimal position vectors and times of an N-impulse trajectory. The solution is examined for the specific boundary conditions of (1) fixed-end conditions, two-impulse, and time-open transfer; (2) an orbit-to-orbit transfer; and (3) a generalized rendezvous problem. A sequence of rendezvous problems is solved to illustrate the analysis and the computational procedure.

Source record↗

NASA NTRS1985

New integration techniques for chemical kinetic rate equations. 2: Accuracy comparison

A comparison of the accuracy of several techniques recently developed for solving stiff differential equations is presented. The techniques examined include two general purpose codes EEPISODE and LSODE developed for an arbitrary system of ordinary differential equations, and three specialized codes CHEMEQ, CREKID, and GCKP84 developed specifically to solve chemical kinetic rate equations. The accuracy comparisons are made by applying these solution procedures 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. The comparisons show that LSODE is 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 an iterative solution of the algebraic enthalpy conservation equation for the temperature can be more accurate and efficient than computing the temperature by integrating its time derivative.

Radhakrishnan, K.↗

NASA NTRS1985

The Formation of Preplanetary Disks from the Collapse of Rotating Molecular Cloud Cores

Solutions that describe the collapse of a molecular cloud core that is initially in unstable equilibrium, embedded within an envelope of uniform density, and rotating at the same rate as the envelope are given. Hydrodynamic equations, including self gravity, are deduced to a set of ordinary differential equations, which are solved by the method of matched asymptotic expansions. Results of these calculations are: (1) the range of stellar masses derived seems to correspond to realistic ranges of observed stellar masses and interstellar cloud parameters, (2) the proper measure of dissipation rate is the ratio of accretion time to viscous diffusion time, and (3) the pressure distribution on the surface of an accreting protostar is nonuniform in a way that favors the channeling of a stellar wind into a bipolar flow directed along the rotation axis.

Cassen, P.↗

NASA NTRS1986

New integration techniques for chemical kinetic rate equations. II - Accuracy comparison

A comparison of the accuracy of several techniques recently developed for solving stiff differential equations is presented. The techniques examined include two general purpose codes EEPISODE and LSODE developed for an arbitrary system of ordinary differential equations, and three specialized codes CHEMEQ, CREKID, and GCKP84 developed specifically to solve chemical kinetic rate equations. The accuracy comparisons are made by applying these solution procedures 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. The comparisons show that LSODE is 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 an iterative solution of the algebraic enthalpy conservation equation for the temperature can be more accurate and efficient than computing the temperature by integrating its time derivative.

Radhakrishnan, K.↗

NASA NTRS1994

Refraction of high frequency noise in an arbitrary jet flow

Refraction of high frequency noise by mean flow gradients in a jet is studied using the ray-tracing methods of geometrical acoustics. Both the two-dimensional (2D) and three-dimensional (3D) formulations are considered. In the former case, the mean flow is assumed parallel and the governing propagation equations are described by a system of four first order ordinary differential equations. The 3D formulation, on the other hand, accounts for the jet spreading as well as the axial flow development. In this case, a system of six first order differential equations are solved to trace a ray from its source location to an observer in the far field. For subsonic jets with a small spreading angle both methods lead to similar results outside the zone of silence. However, with increasing jet speed the two prediction models diverge to the point where the parallel flow assumption is no longer justified. The Doppler factor of supersonic jets as influenced by the refraction effects is discussed and compared with the conventional modified Doppler factor.

Khavaran, Abbas↗

NASA NTRS1973

Wind-driven currents in a shallow lake or sea.

A technique is given for obtaining exact analytical solutions for wind-driven currents in a shallow lake or sea with a large class of bottom topographies which may or may not contain an island. Welander's second-order partial differential equation for the lake's surface displacement is used for determining the lake circulation. The technique, together with some results based on conformal mapping, is used to reduce the problem of finding solutions to Welander's equation for a closed body of water to the problem of solving an ordinary differential equation. The method is applied to an elliptically shaped lake and a circular lake containing an eccentrically located circular island.

Goldstein, M.↗

NASA NTRS1985

Quasiperiodically forced damped pendula and Schroedinger equations with quasiperiodic potentials - Implications of their equivalence

Certain first-order nonlinear ordinary differential equations exemplified by strongly damped, quasiperiodically driven pendula and Josephson junctions are isomorphic to Schroedinger equations with quasiperiodic potentials. The implications of this equivalence are discussed. In particular, it is shown that the transition to Anderson localization in the Schroedinger problem corresponds to the occurrence of a novel type of strange attractor in the pendulum problem. This transition should be experimentally observable in the frequency spectrum of the pendulum of Josephson junction.

Bondeson, A.↗