This paper presents a family of Runge{Kutta type integration schemes of arbitrarily high order for diierential equations evolving on manifolds. We prove that any classical Runge{Kutta method can be turned into an invariant method of the same order on a general homogeneous manifold, and present a family of algorithms that are relatively simple to implement.

In this paper we show that standard preconditioners for parabolic PDEs discretized by implicit Euler or Crank–Nicolson schemes can be reused for higher–order fully implicit Runge–Kutta time discretization schemes. We prove that the suggested block diagonal preconditioners are order–optimal for A–stable and irreducible Runge–Kutta schemes with invertible coefficient matrices. The theoretical inv...

In this paper, we analyze second-order Runge–Kutta approximations to a nonlinear optimal control problem with control constraints. If the optimal control has a derivative of bounded variation and a coercivity condition holds, we show that for a special class of Runge–Kutta schemes, the error in the discrete approximating control is O(h2) where h is the mesh spacing.

The non-linear wave equation is taken as a model problem for the investigation. Different multisymplectic reformulations of the equation are discussed. Multi-symplectic Runge–Kutta methods and multi-symplectic partitioned Runge–Kutta methods are explored based on these different reformulations. Some popular and efficient multi-symplectic schemes are collected and constructed. Stability analyses...

The Maxwell equations are solved by a long-stencil fourth order finite difference method over a Yee grid, in which different physical variables are located at staggered mesh points. A careful treatment of the numerical values near the boundary is introduced, which in turn leads to a “symmetric image” formula at the “ghost” grid points. Such a symmetric formula assures the stability of the bound...

We are concerned with the solution of time-dependent nonlinear hyperbolic partial differential equations. We investigate the combination of residual distribution methods with a consistent mass matrix (discretisation in space) and a Runge-Kutta-type time stepping (discretisation in time). The introduced nonlinear blending procedure allows us to retain the explicit character of the time stepping ...

In this paper, we will present an effective simulation to study the solution behavior of a high dimensional chaos by considering nine-dimensionalLorenz system through Rabotnov fractional-exponential (RFE) kernel fractional derivative. First, derive approximate formula thefractional-order derivative polynomial function $t^{p}$ in terms RFE kernel. work, use spectral collocation method basedon pr...