Search NASA⌕ Search

SEARCH · Search NASA

Results for “Runge-Kutta method”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

Some self starting integrators for x Prime equals f (x, t)

The integration is discussed of the vector differential equation X = F(x, t) from time t sub i to t sub (i = 1) where only the values of x sub i are available for the the integration. No previous values of x or x prime are used. Using an orbit integration problem, comparisons are made between Taylor series integrators and various types and orders of Runge-Kutta integrators. A fourth order Runge-Kutta type integrator for orbital work is presented, and approximate (there may be no exact) fifth order Runge-Kutta integrators are discussed. Also discussed and compared is a self starting integrator ising delta f/delta x. A numerical method for controlling the accuracy of integration is given, and the special equations for accurately integrating accelerometer data are shown.

Lear, W. M.↗
NASA NTRS1989

Numerical modeling of turbulent supersonic reacting coaxial jets

The paper considers the mixing and subsequent combustion within turbulent reacting shear layers. A computer program was developed to solve the axisymmetric Reynolds averaged Navier-Stokes equations. The numerical method integrates the Reynolds averaged Navier-Stokes equations using a finite volume approach while advancing the solution forward in time using a Runge-Kutta scheme. Three separate flowfields are investigated and it is found that no single turbulence model considered could accurately predict the degree of mixing for all three cases.

Eklund, Dean R.↗
NASA NTRS1989

Design of optimally smoothing multi-stage schemes for the Euler equations

In this paper, a method is developed for designing multi-stage schemes that give optimal damping of high-frequencies for a given spatial-differencing operator. The objective of the method is to design schemes that combine well with multi-grid acceleration. The schemes are tested on a nonlinear scalar equation, and compared to Runge-Kutta schemes with the maximum stable time-step. The optimally smoothing schemes perform better than the Runge-Kutta schemes, even on a single grid. The analysis is extended to the Euler equations in one space-dimension by use of 'characteristic time-stepping', which preconditions the equations, removing stiffness due to variations among characteristic speeds. Convergence rates independent of the number of cells in the finest grid are achieved for transonic flow with and without a shock. Characteristic time-stepping is shown to be preferable to local time-stepping, although use of the optimally damping schemes appears to enhance the performance of local time-stepping. The extension of the analysis to the two-dimensional Euler equations is hampered by the lack of a model for characteristic time-stepping in two dimensions. Some results for local time-stepping are presented.

Van Leer, Bram↗
NASA NTRS1986

Computation of three-dimensional, rotational flow through turbomachinery blade rows for improved aerodynamic design studies

The performance of a three dimensional computer code developed for predicting the flowfield in stationary and rotating turbomachinery blade rows is described in this study. The four stage Runge-Kutta numerical integration scheme is used for solving the governing flow equations and yields solution to the full, three dimensional, unsteady Euler equations in cylindrical coordinates. This method is fully explicit and uses the finite volume, time marching procedure. In order to demonstrate the accuracy and efficiency of the code, steady solutions were obtained for several cascade geometries under widely varying flow conditions. Computed flowfield results are presented for a fully subsonic turbine stator and a low aspect ratio, transonic compressor rotor blade under maximum flow and peak efficiency design conditions. Comparisons with Laser Anemometer measurements and other numerical predictions are also provided to illustrate that the present method predicts important flow features with good accuracy and can be used for cost effective aerodynamic design studies.

Subramanian, S. V.↗
NASA NTRS1985

SSME Turbopump Turbine Computations

A two-dimensional viscous code was developed to be used in the prediction of the flow in the SSME high-pressure turbopump blade passages. The rotor viscous code (RVC) employs a four-step Runge-Kutta scheme to solve the two-dimensional, thin-layer Navier-Stokes equations. The Baldwin-Lomax eddy-viscosity model is used for these turbulent flow calculations. A viable method was developed to use the relative exit conditions from an upstream blade row as the inlet conditions to the next blade row. The blade loading diagrams are compared with the meridional values obtained from an in-house quasithree-dimensional inviscid code. Periodic boundary conditions are imposed on a body-fitted C-grid computed by using the GRAPE GRids about Airfoils using Poisson's Equation (GRAPE) code. Total pressure, total temperature, and flow angle are specified at the inlet. The upstream-running Riemann invariant is extrapolated from the interior. Static pressure is specified at the exit such that mass flow is conserved from blade row to blade row, and the conservative variables are extrapolated from the interior. For viscous flows the noslip condition is imposed at the wall. The normal momentum equation gives the pressure at the wall. The density at the wall is obtained from the wall total temperature.

Jorgenson, P. G. E.↗
NASA NTRS1985

Convergence acceleration of iterative solutions of Euler equations for transonic flow computations

In this paper, two acceleration techniques for Euler calculations are investigated. The first technique is an extrapolation procedure based on the Power Method; it is applicable when the iterative matrix has dominant eigenvalues. Both real and complex conjugate roots are allowed. The second technique is a generalization of the Minimal Residual Method, where the extrapolation step consists of a weighted combination of the corrections at different iteration levels and the weights are chosen to minimize the Lz norm of the residual. Numerical results, using Jameson's Runge-Kutta Multigrid Code, are presented. The extra computational work to apply either technique is negligible and the extra storage is not a problem on current supercomputers.

Hafez, M.↗
NASA NTRS1988

Application of advanced grid generation techniques for flow field computations about complex configurations

In the computation of flowfields about complex configurations, it is very difficult to construct a boundary-fitted coordinate system. An alternative approach is to use several grids at once, each of which is generated independently. This procedure is called the multiple grids or zonal grids approach; its applications are investigated. The method conservative providing conservation of fluxes at grid interfaces. The Euler equations are solved numerically on such grids for various configurations. The numerical scheme used is the finite-volume technique with a three-stage Runge-Kutta time integration. The code is vectorized and programmed to run on the CDC VPS-32 computer. Steady state solutions of the Euler equations are presented and discussed. The solutions include: low speed flow over a sphere, high speed flow over a slender body, supersonic flow through a duct, and supersonic internal/external flow interaction for an aircraft configuration at various angles of attack. The results demonstrate that the multiple grids approach along with the conservative interfacing is capable of computing the flows about the complex configurations where the use of a single grid system is not possible.

Kathong, Monchai↗
DOE OSTI2025

High order interpolation of magnetic fields with vector potential reconstruction for particle simulations

We propose a method for interpolating divergence-free continuous magnetic fields via vector potential reconstruction using Hermite interpolation, which ensures high-order continuity for applications requiring adaptive, high-order ordinary differential equation (ODE) integrators, such as the Dormand-Prince method. The method provides C(m) continuity and achieves high-order accuracy, making it particularly suited for particle trajectory integration and Poincaré section analysis under optimal integration order and timestep adjustments. Through numerical experiments, we demonstrate that the Hermite interpolation method preserves volume and continuity, which are critical for conserving toroidal canonical momentum and magnetic moment in guiding center simulations, especially over long-term trajectory integration. Furthermore, we analyze the impact of insufficient derivative continuity on Runge-Kutta schemes and show how it degrades accuracy at low error tolerances, introducing discontinuity-induced truncation errors. Lastly, we demonstrate performant Poincaré section analysis in two relevant settings of field data collocated from finite element meshes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1988

A conservative approach for flow field calculations on multiple grids

In the computation of flow fields about complex configurations, it is very difficult to construct body-fitted coordinate systems. An alternative approach is to use several grids at once, each of which is generated independently. This procedure is called the 'multiple grids' or 'zonal grids' approach and its applications are investigated in this study. The method follows the conservative approach and provides conservation of fluxes at grid interfaces. The Euler equations are solved numerically on such grids for various configurations. The numerical scheme used is the finite-volume technique with a three-stage Runge-Kutta time integration. The code is vectorized and programmed to run on the CDC VPS-32 computer. Some steady state solutions of the Euler equations are presented and discussed.

Kathong, M.↗
NASA NTRS1988

A conservative approach for flow field calculations on multiple grids

In the computation of flow fields about complex configurations, it is very difficult to construct body-fitted coordinate systems. An alternative approach is to use several grids at once, each of which is generated independently. This procedure is called the multiple grids or zonal grids approach and its applications are investigated in this study. The method follows the conservative approach and provides conservation of fluxes at grid interfaces. The Euler equations are solved numerically on such grids for various configurations. The numerical scheme used is the finite-volume technique with a three-state Runge-Kutta time integration. The code is vectorized and programmed to run on the CDC VPS-32 computer. Some steady state solutions of the Euler equations are presented and discussed.

Kathong, Monchai↗
NASA NTRS1988

Essentially non-oscillatory shock capturing methods applied to turbulence amplification in shock wave calculations

ENO (essentially non-oscillatory) schemes can provide uniformly high order accuracy right up to discontinuities while keeping sharp, essentially non-oscillatory shock transitions. Recently, an efficient implementation of ENO schemes was obtained based on fluxes and TVD Runge-Kutta time discretizations. The resulting code is very simple to program for multi-dimensions. ENO schemes are especially suitable for computing problems with both discontinuities and fine structures in smooth regions, such as shock interaction with turbulence, for which results for 1-D and 2-D Euler equations are presented. Much better resolution is observed by using third order ENO schemes than by using second order TVD schemes for such problems.

Osher, Stanley↗
NASA NTRS1988

Efficient implementation of essentially non-oscillatory shock capturing schemes, 2

Earlier work on the efficient implementation of ENO (essentially non-oscillatory) shock capturing schemes is continued. A new simplified expression is provided for the ENO construction procedure based again on numerical fluxes rather than cell averages. Also considered are two improvements which are labeled ENO-LLF (local Lax-Friedrichs) and ENO-Roe, which yield sharper shock transitions, improved overall efficiency, and lower computational cost than previous implementation of the ENO schemes. Two methods of sharpening contact discontinuities, i.e., the subcell resolution idea of Harten and the artificial compression idea of Yang, which those authors used originally in the cell-average framework, are supplied to the current ENO schemes using numerical fluxes and TVD Runge-Kutta time discretizations. The implementation for nonlinear systems and multi-dimensions is given. Finally, many numerical examples, including a compressible shock turbulence interaction flow calculation, are given.

Shu, Chi-Wang↗
NASA NTRS1986

Euler solutions for high-speed flow about complex three-dimensional configurations

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.

Moitra, A.↗
DOE OSTI2024

A high-order explicit Runge-Kutta approximation technique for the shallow water equations

Here, we introduce a high-order space–time approximation of the Shallow Water Equations with sources that is invariant-domain preserving (IDP), well-balanced with respect to rest states, and employs a novel explicit Runge–Kutta (ERK) introduced in Ern and Guermond (SIAM J. Sci. Comput. 44(5), A3366–A3392, 2022) for systems of non-linear conservation equations. The resulting method is then numerically illustrated through verification and validation.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1986

Development of an explicit multigrid algorithm for quasi-three-dimensional viscous flows in turbomachinery

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.

Chima, R. V.↗
NASA NTRS1985

Development of an explicit multigrid algorithm for quasi-three-dimensional viscous flows in turbo-machinery

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 by 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 scheme. 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.

Chima, R. V.↗
NASA NTRS1989

Non linear evolution of a second mode wave in supersonic boundary layers

Presented here are several direct simulations of one 2-D second mode perturbation wave, superimposed upon a prescribed mean flow. Periodicity is assumed in the streamwise direction (Fourier) and the variables are expanded in Chebyshev series in the direction normal to the flat plate. The code is fully explicit and is time advanced with a 3rd order Runge-Kutta scheme. The second mode wave (R delta prime = 8000), interacts with itself to generate higher streamwise harmonics. Physical parameters are chosen to maximize the linear growth rate at the prescribed Reynolds number. Initial results indicate that the nonlinear processes begin in the critical layer region and are the result of the cubic interactions in the momentum equations, rather than due to the higher streamwise harmonics. Analysis of the various terms in the momentum equations combined with numerical experiments in which various modes are artificially suppressed, lead to the conclusion that asymptotic methods will produce the saturated state in one or two order of magnitude less computer time than that required by the direct numerical simulations.

Erlebacher, Gordon↗
NASA NTRS1986

Particle trajectory computation on a 3-dimensional engine inlet

A 3-dimensional particle trajectory computer code was developed to compute the distribution of water droplet impingement efficiency on a 3-dimensional engine inlet. The computed results provide the essential droplet impingement data required for the engine inlet anti-icing system design and analysis. The droplet trajectories are obtained by solving the trajectory equation using the fourth order Runge-Kutta and Adams predictor-corrector schemes. A compressible 3-D full potential flow code is employed to obtain a cylindrical grid definition of the flowfield on and about the engine inlet. The inlet surface is defined mathematically through a system of bi-cubic parametric patches in order to compute the droplet impingement points accurately. Analysis results of the 3-D trajectory code obtained for an axisymmetric droplet impingement problem are in good agreement with NACA experimental data. Experimental data are not yet available for the engine inlet impingement problem analyzed. Applicability of the method to solid particle impingement problems, such as engine sand ingestion, is also demonstrated.

Kim, J. J.↗