Search NASA⌕ Search

SEARCH · Search NASA

Results for “RUNGE-KUTTA INTEGRATION”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Apollo experience report: Onboard navigational and alignment software

The onboard navigational and alignment routines used during the nonthrusting phases of an Apollo mission are discussed as to their limitations, and alternate approaches that have more desirable capabilities are presented. A more efficient procedure for solving Kepler's equation, which is used in the calculation of Kepler's problem and Lambert's problem is included, and a sixth-order predictor scheme with a Runge-Kutta starter is recommended for numerical integration. The extension of the rendezvous navigation state to include angle biases and the use of a fixed coordinate system is also evaluated.

Savely, R. T.↗
NASA NTRS1977

An optimal Space Shuttle ascent trajectory for the first orbital flight test

An optimal solution of the ascent trajectory of the Space Shuttle for the first orbital flight test is presented; the optimization is a minimum propellant, four-control problem in yaw angle, roll angle, pitch angle and vacuum thrust of each Space Shuttle main engine. Piecewise linear segments with juncture points treated as parameters are employed to model the controls. Equations of motion for a three-dimensional flight with pitch plane moment balance about an oblate are integrated numerically with a fourth-order Runge-Kutta method; two- and one-dimensional cubic spline function curve fits of aerodynamic coefficients are used during the first and second stages, respectively. The constraint minimization problem is solved with the Davidon-Fletcher-Powell function method.

Johnson, I. L., Jr.↗
NASA NTRS1979

Shuttle program. Onorbit navigation integrator results for typical shuttle orbits

Three types of navigation onorbit numerical integrators were evaluated: (1) power integrators with no delta-V incorporation, just coasting (using Taylor series expansion integrators); (2) coasting integrators using the Cowell method of special perturbations; and (3) coasting integrator using the Pines variation of parameter perturbation method. Results show that the super G integrator is a very simple and effective for 2 and 4 second time steps. Since IMU delta-V data can be easily incorporated in the integration scheme, its use as the standard onorbit navigation propagator for the maintenance of the current state was implemented in the onboard navigation software. The Pines formulation method with a Runge-Kutta-Gill fourth-order integrator method produces excellent results up to 300 second time steps. On orbit prediction with this method was implemented in the onboard onorbit navigation scheme. The Runge-Kutta third order, using Cowell's method, is an excellent general purpose determination integrator for time steps up to a 60 second duration.

Olszewski, O. W.↗
NASA NTRS1972

Computer program for investigating effects of nonlinear suspension-system elastic properties on parachute inflation loads and motions

A computer program is presented by which the effects of nonlinear suspension-system elastic characteristics on parachute inflation loads and motions can be investigated. A mathematical elastic model of suspension-system geometry is coupled to the planar equations of motion of a general vehicle and canopy. Canopy geometry and aerodynamic drag characteristics and suspension-system elastic properties are tabular inputs. The equations of motion are numerically integrated by use of an equivalent fifth-order Runge-Kutta technique.

Poole, L. R.↗
NASA NTRS1981

Elasto visco-plastic flow with special attention to boundary conditions

A simple but nontrivial steady-state creeping elasto visco-plastic (Maxwell fluid) radial flow problem is analyzed, with special attention given to the effects of the boundary conditions. Solutions are obtained through integration of a governing equation on stress using the Runge-Kutta method for initial value problems and finite differences for boundary value problems. A more general approach through the finite element method, an approach that solves for the velocity field rather than the stress field and that is applicable to a wide range of problems, is presented and tested using the radial flow example. It is found that steady-state flows of elasto visco-plastic materials are strongly influenced by the state of stress of material as it enters the region of interest. The importance of this boundary or initial condition in analyses involving materials coming into control volumes from unusual stress environments is emphasized.

Shimazaki, Y.↗
NASA NTRS1982

Scaled Runge-Kutta algorithms for treating the problem of dense output

A set of scaled Runge-Kutta algorithms for the third- through fifth-orders are developed to determine the solution at any point within the integration step at a relatively small increase in computing time. Each scaled algorithm is designed to be used with an existing Runge-Kutta formula, using the derivative evaluations of the defining algorithm along with an additional derivative evaluation (or two). Third-order, scaled algorithms are embedded within the existing formulas at no additional derivative expense. Such algorithms can easily be adopted to generate interpolating polynomials (or dependent variable stops) efficiently.

Horn, M. K.↗
NASA NTRS1972

Effect of support flexibilty and damping on the dynamic response of a single mass flexible rotor in elastic bearings

A steady state analysis of the shaft and the bearing housing motion was made by assuming synchronous precession of the system. The conditions under which the support system would act as a dynamic vibration absorber at the rotor critical speed were studied; plots of the rotor and support amplitudes, phase angles, and forces transmitted were evaluated by the computer, and the performance curves were automatically plotted by a CalComp plotter unit. Curves are presented on the optimization of the support housing characteristics to attenuate the rotor unbalance response over the entire rotor speed range. The complete transient motion including rotor unbalance was examined by integrating the equations of motion numerically using a modified fourth order Runge-Kutta procedure, and the resulting whirl orbits were plotted by the CalComp plotter unit. The results of the transient analysis are discussed with regards to the design optimization procedure derived from the steady-state analysis.

Kirk, R. G.↗
NASA NTRS1974

Comparison of numerical integration techniques for orbital applications

The present work gives a brief comparison of the performance of programs for integrating differential equations for orbital applications. The evaluation criteria and the method of testing are described, and the results of the test problem set are included. Integration methods that were chosen for comparison include high-order Runge-Kutta methods; a rational extrapolation method (Bulirsch and Stoer, 1966); a variable step, variable order, multistep method (Krogh, 1969); classical multistep methods of Adams and Cowell; and modified multistep methods. The high-order Runge-Kutta methods used in the comparison include RKF 7(8) and RKF 8(9) (Fehlberg, 1968), and RKS 8-10 (Shanks, 1966).

Moore, H.↗
NASA NTRS1979

Nonlinear wave choked inlets

The quasi-one dimensional flow program was modified in two ways. The Runge-Kutta subroutine was replaced with a subroutine which used a modified divided difference form of the Adams Pece method and the matrix inversion routine was replaced with a pseudo inverse routine. Calculations were run using both the original and modified programs. Comparison of the calculations showed that the original Runge-Kutta routine could not detect singularity near the throat and was integrating across it. The modified version was able to detect the singularity and therefore gave more valid calculations.

Source record↗
NASA NTRS1971

Numerical integration of second order differential equations

Performance characteristics of higher order approximations of Runge-Kutta type are analyzed, and performance predictors for time required on machine and for error size are developed. Technique is useful in evaluating system performance, analyzing material characteristics, and designing inertial guidance and nuclear instrumentation and materials.

Shanks, E. B.↗
NASA NTRS1976

Conservation properties of numerical integration methods for systems of ordinary differential equations

If a system of ordinary differential equations represents a property conserving system that can be expressed linearly (e.g., conservation of mass), it is then desirable that the numerical integration method used conserve the same quantity. It is shown that both linear multistep methods and Runge-Kutta methods are 'conservative' and that Newton-type methods used to solve the implicit equations preserve the inherent conservation of the numerical method. It is further shown that a method used by several authors is not conservative.

Rosenbaum, J. S.↗
NASA NTRS1981

Modified Runge-Kutta methods for solving ODES

A class of Runge-Kutta formulas is examined which permit the calculation of an accurate solution anywhere in the interval of integration. This is used in a code which seldom has to reject a step; rather it takes a reduced step if the estimated error is too large. The absolute stability implications of this are examined.

Vanvu, T.↗
NASA NTRS1977

An ILLIAC program for the numerical simulation of homogeneous incompressible turbulence

An algorithm and ILLIAC computer program, developed for the simulation of homogeneous incompressible turbulence in the presence of an applied mean strain, are described. The turbulence field is represented spatially by a truncated triple Fourier series (spectral method) and followed in time using a fourth-order Runge-Kutta algorithm. These include: (1) transformation of variables suggested by Taylor's sudden-distortion theory; (2) implicit viscous diffusion by use of an integrating factor; (3) implicit pressure calculation suggested by Taylor's sudden-distortion theory, and (4) inexpensive control of aliasing by random and phased coordinate shifts.

Rogallo, R. S.↗
NASA NTRS1980

On time discretizations for spectral methods

New methods are introduced for the time integration of the Fourier and Chebyshev methods of solution for dynamic differential equations. These methods are unconditionally stable, even though no matrix inversions are required. Time steps are chosen by accuracy requirements alone. For the Fourier method both leapfrog and Runge-Kutta methods are considered. For the Chebyshev method only Runge-Kutta schemes are tested. Numerical calculations are presented to verify the analytic results. Applications to the shallow water equations are presented.

Gottlieb, D.↗
NASA NTRS1972

On the stability of numerical integration routines.

Numerical integration methods for the solution of initial value problems for ordinary vector differential equations may be modelled as discrete time feedback systems. The stability criteria discovered in modern control theory are applied to these systems and criteria involving the routine, the step size and the differential equation are derived. Linear multistep, Runge-Kutta, and predictor-corrector methods are all investigated.

Glover, K.↗
NASA NTRS1973

On the stability of numerical integration routines for ordinary differential equations.

Numerical integration methods for the solution of initial value problems for ordinary vector differential equations may be modelled as discrete time feedback systems. The stability criteria discovered in modern control theory are applied to these systems and criteria involving the routine, the step size and the differential equation are derived. Linear multistep, Runge-Kutta, and predictor-corrector methods are all investigated.

Glover, K.↗