Explicit Runge-Kutta integration
Fourth order Runge-Kutta method applied to system of linear differential equations
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.
Fourth order Runge-Kutta method applied to system of linear differential equations
A quantitative comparison between the Euler and full potential formulations with respect to speed and accuracy is presented. The robustness of the codes used is tested by a number of transonic airfoil cases. The computed results are from four transonic airfoil computer codes. The full potential codes use fully implicit iteration algorithms. The first Euler code uses a fully implicit ADI iteration scheme. The second Euler code uses an explicit Runge Kutta time stepping algorithm which is enhanced by a multigrid convergence acceleration scheme. Quantitative comparisons are made using various plots of lift coefficient versus the average mesh spacing along the airfoil. Besides yielding an asymptotic limit to the lift coefficient, these results also demonstrate the truncation error behavior of the various codes. Quantitative conclusions regarding the full potential and Euler formulations with respect to accuracy, speed, and robustness can be presented.
The A-contractivity of Runge-Kutta methods with respect to an inner product norm was investigated thoroughly by Butcher and Burrage (who used the term B-stability). Their theory is extended to contractivity in a region bounded by a circle through the origin. The largest possible circle is calculated for many known explicit Runge-Kutta methods. As a rule it is considerably smaller than the stability region, and in several cases it degenerates to a point. It is shown that an explicit Runge-Kutta method cannot be contractive in any circle of this class if it is more than fourth order accurate.
Singular fifth order Runge-Kutta solution conditions for differential system
To identify the fastest algorithm currently available for the numerical integration of chemical kinetic rate equations, several algorithms were examined. Findings to date are summarized. The algorithms examined include two general-purpose codes EPISODE and LSODE and three special-purpose (for chemical kinetic calculations) codes CHEMEQ, CRK1D, and GCKP84. In addition, an explicit Runge-Kutta-Merson differential equation solver (IMSL Routine DASCRU) is used to illustrate the problems associated with integrating chemical kinetic rate equations by a classical method. Algorithms were applied to two test problems drawn from combustion kinetics. These problems included all three combustion regimes: induction, heat release and equilibration. Variations of the temperature and species mole fraction are given with time for test problems 1 and 2, respectively. Both test problems were integrated over a time interval of 1 ms in order to obtain near-equilibration of all species and temperature. Of the codes examined in this study, only CREK1D and GCDP84 were written explicitly for integrating exothermic, non-isothermal combustion rate equations. These therefore have built-in procedures for calculating the temperature.
Three families of classical fifth-order Runge- Kutta formulas including Newton-Cotes, Legendre-Gauss and Lobatto and Randau types
Sixth order Runge-Kutta formula for differential equation system with vector valued functions
The near-term prospects are discussed for calculation of viscous transonic flow fields about realistic configurations at full-scale Reynolds numbers. Three basic algorithms are considered: The central-difference, three-factor ADI method; the central-difference, explicit, multistep Runge-Kutta method with multigrid acceleration; and the relaxation method for the upwind-differenced, flux-split equations. Each method has distinct advantages and disadvantages regarding stability, convergence rate, and vectorizability. It appears that computation times can be 15 to 60 hours on the latest super computers unless 3D algorithms are improved to perform as well as current 2D algorithms.
An explicit multistage Runge-Kutta type time-stepping scheme is used for solving the three-dimensional, compressible, thin-layer Navier-Stokes equations. A finite-volume formulation is employed to facilitate treatment of complex grid topologies encountered in three-dimensional calculations. Convergence to steady state is expedited through usage of acceleration techniques. Further numerical efficiency is achieved through vectorization of the computer code. The accuracy of the overall scheme is evaluated by comparing the computed solutions with the experimental data for a finite wing under different test conditions in the transonic regime. A grid refinement study ir conducted to estimate the grid requirements for adequate resolution of salient features of such flows.
Simulation of high-speed three-dimensional flow about blended wing-body combinations is investigated. A finite-volume explicit scheme with Runge-Kutta time integration is used to solve the compressible Euler equations in order to simulate the flow. The method, augmented by carefully chosen dissipative terms and convergence accelerators such as enthalpy damping and maximum local time-stepping, has been found to be very efficient in solving high-speed flows involving strong shocks. An analytic method for generating the body geometry and an algebraic method for quasi-three-dimensional grid generation are described. Results are presented for various blended wing-body configurations at different Mach numbers and angles of attack.
A numerical algorithm based on a finite-volume explicit scheme with Runge-Kutta time integration of the Euler equations is presented for calculating high-speed three-dimensional flow about complex aerospace configurations. The use of enhancing factors such as artificial dissipative terms, enthalpy damping and local time-stepping are described. An algebraic method for generating quasi-three-dimensional computational grids for realistic aerospace configurations is presented. Computed results for various three-dimensional bodies at different Mach numbers and angles of attack have been obtained using the methods for grid-generation and flow simulation. Comparison of computed and experimental data for an advanced tactical aircraft-like configuration is presented, and a reasonable agreement of the data is noticed.
A new method is derived for solving parabolic partial differential equations arising in transient heat conduction or in boundary-layer flows. The method is based on a combination of the modified differential quadrature (MDQ) method with the rational Runge-Kutta time-integration scheme. It is fully explicit, requires no matrix inversion, and is stable for any time-step for the heat equations. Burgers equation and the one- and two-dimensional heat equations are solved to demonstrate the accuracy and efficiency of the proposed algorithm. The present method is found to be very accurate and efficient when results are compared with analytic solutions.
Tables for elementary weights of Runge-Kutta formulas of first eight orders and for relations of explicit formulas through order seven
A numerical model of eddy diffusive transport of emitted gases from the space shuttle wake, including chemical reactions between the emitted constituents and the ambient atmosphere, has been constructed for 75 km altitude. The numerical methods involve explicit solution of the diffusion equation and Runge-Kutta method for the chemical reactions. The time required to reach background levels of nitric oxide concentration of 7 x 10 to the 7th power molecules/cc has been calculated. This relaxation time depends strongly on atmospheric conditions.
Numerical solutions of the unsteady Euler equations are obtained using the classical fourth order Runge Kutta time marching scheme. This method is fully explicit and is applied to the governing equations in the finite volume, conservation law form. In order to determine the efficiency of this scheme for solving turbomachinery flows, steady blade-to-blade solutions are obtained for compressor and turbine cascades under subsonic and transonic flow conditions. Computed results are compared with other numerical methods and wind tunnel measurements. The present study also focuses on other important numerical aspects influencing the performance of the algorithm and the solution accuracy such as grid types, boundary conditions, and artificial viscosity. For this purpose, H, O, and C type computational grids as well as characteristic and extrapolation type boundary conditions are included in the solution procedure.
Explicit integration of differential equations
The solution of the three dimensional flow field for a flow through nacelle was studied. Both inviscid and viscous inviscid interacting solutions were examined. Inviscid solutions were obtained with two different computational procedures for solving the three dimensional Euler equations. The first procedure employs an alternating direction implicit numerical algorithm, and required the development of a complete computational model for the nacelle problem. The second computational technique employs a fourth order Runge-Kutta numerical algorithm which was modified to fit the nacelle problem. Viscous effects on the flow field were evaluated with a viscous inviscid interacting computational model. This model was constructed by coupling the explicit Euler solution procedure with a flag entrainment boundary layer solution procedure in a global iteration scheme. The computational techniques were used to compute the flow field for a long duct turbofan engine nacelle at free stream Mach numbers of 0.80 and 0.94 and angles of attack of 0 and 4 deg.
A rapid quasi three-dimensional analysis was developed for blade-to-blade flows in turbomachinery. The analysis solves the unsteady Euler or thin layer Navier-Stokes equations in a body-fitted coordinate system. It accounts for the effects of rotation, radius change, and stream-surface thickness. The Baldwin-Lomax eddy-viscosity model is used for turbulent flows. The equations which are solved b a two-stage Runge-Kutta scheme made efficient by use of vectorization, a variable time-step, and a flux-based multigrid scheme, are described. A stability analysis is presented for the two-stage. Results for a flat-plate model problem show the applicability of the method to axial, radial, and rotating geometries. Results for a centrifugal impeller and a radial diffuser show that the quasi three-dimensional viscous analysis can be a practical design tool.