This paper is concerned with time-stepping numerical methods for computing sti semi-discrete systems of ordinary di erential equations for transient hypersonic ows with thermo-chemical nonequilibrium. The sti ness of the equations is mainly caused by the viscous ux terms across the boundary layers and by the source terms modeling nite-rate thermo-chemical processes. Implicit methods are needed ...

A popular approach for the numerical solution of boundary value ODE problems involves the use of collocation methods. Such methods can be naturally implemented so as to provide a continuous approximation to the solution over the entire problem interval. On the other hand, several authors have suggested as an alternative, certain subclasses of the implicit Runge-Kutta formulas, known as mono-imp...

A one-step 4-stage Hermite-Birkhoff-Taylor method of order 12, denoted by HBT(12)4, is constructed for solving nonstiff systems of first-order differential equations of the form y′ = f(x, y), y(x0) = y0. The method uses derivatives y′ to y(9) as in Taylor methods combined with a 4-stage Runge-Kutta method. Forcing an expansion of the numerical solution to agree with a Taylor expansion of the tr...

Practical, structure-preserving methods for integrating classical Heisenberg spin systems are discussed. Two new integrators are derived and compared, including (1) a symmetric energy and spin-length preserving integrator based on a Red-Black splitting of the spin sites combined with a staggered timestepping scheme and (2) a (Lie-Poisson) symplectic integrator based on Hamiltonian splitting. Th...

A general class of one-step methods for index 2 differential-algebraic systems in Hessenberg form is studied. This family of methods, which we call partitioned Runge-Kutta methods, includes all one-step methods of Runge-Kutta type proposed in the literature for integrating such DAE systems, including the more recently proposed classes of half-explicit methods. A new family of super-convergent p...

We apply a Runge-Kutta-based waveform relaxation method to initial-value problems for implicit differential equations. In the implementation of such methods, a sequence of nonlinear systems has to be solved iteratively in each step of the integration process. The size of these systems increases linearly with the number of stages of the underlying Runge-Kutta method, resulting in high linear alg...

Hamiltonian systems arise in many areas of physics, mechanics, and engineering sciences as well as in pure and applied mathematics. To their symplectic integration certain Runge–Kutta– type methods are profitably applied (see Sanz–Serna and Calvo [10]). In this paper Runge–Kutta and partitioned Runge–Kutta methods are considered. Different features of symmetry are distinguished using reflected ...

We present a new class of adaptivity algorithms for time-dependent partial differential equations (PDE) that combines adaptive higher-order finite elements (hp-FEM) in space with arbitrary (embedded, higher-order, implicit) Runge-Kutta methods in time. Weak formulation is only created for the stationary residual of the equation, and the Runge-Kutta method is supplied via its Butcher’s table. Ar...

Runge-Kutta methods are the classic family of solvers for ordinary differential equations (ODEs), and the basis for the state of the art. Like most numerical methods, they return point estimates. We construct a family of probabilistic numerical methods that instead return a Gauss-Markov process defining a probability distribution over the ODE solution. In contrast to prior work, we construct th...