Aromatic B-series were introduced as an extension of standard Butcher-series for the study of volume-preserving integrators. It was proven with their help that the only volume-preserving B-series method is the exact flow of the differential equation. The question was raised whether there exists a volume-preserving integrator that can be expanded as an aromatic B-series. In this work, we introduce a new algebraic tool, called the aromatic bicomplex, similar to the variational bicomplex in variational calculus. We prove the exactness of this bicomplex and use it to describe explicitly the key object in the study of volume-preserving integrators: the aromatic forms of vanishing divergence. The analysis provides us with a handful of new tools to study aromatic B-series, gives insights on the process of integration by parts of trees, and allows to describe explicitly the aromatic B-series of a volume-preserving integrator. In particular, we conclude that an aromatic Runge–Kutta method cannot preserve volume.

Recent higher-order explicit Runge–Kutta methods are compared with the classic fourth-order (RK4) method in long-term integration of both energy-conserving and lossy systems. By comparing quantity of function evaluations against accuracy for systems with and without known solutions, optimal methods are proposed. For a conservative system, we consider positional accuracy for Newtonian systems of two or three bodies and total angular momentum for a simplified Solar System model, over moderate astronomical timescales (tens of millions of years). For a nonconservative system, we investigate a relativistic two-body problem with gravitational wave emission. We find that methods of tenth and twelfth order consistently outperform lower-order methods for the systems considered here.

In this work, we introduce an IMEX discontinuous Galerkin solver for the weakly compressible isentropic Euler equations. The splitting needed for the IMEX temporal integration is based on the recently introduced reference solution splitting [32, 52], which makes use of the incompressible solution. We show that the overall method is asymptotic preserving. Numerical results show the performance of the algorithm which is stable under a convective CFL condition and does not show any order degradation.

This paper is devoted to an extension of the finite-energy condition for extended Runge-Kutta-Nyström (ERKN) integrators and applications to nonlinear wave equations. We begin with an error analysis for the integrators for multi-frequency highly oscillatory systems , where M is positive semi-definite, . The highly oscillatory system is due to the semi-discretisation of conservative, or dissipative, nonlinear wave equations. The structure of such a matrix M and initial conditions are based on particular spatial discretisations. Similarly to the error analysis for Gaustchi-type methods of order two, where a finite-energy condition bounding amplitudes of high oscillations is satisfied by the solution, a finite-energy condition for the semi-discretisation of nonlinear wave equations is introduced and analysed. These ensure that the error bound of ERKN methods is independent of . Since stepsizes are not restricted by frequencies of M, large stepsizes can be employed by our ERKN integrators of arbitrary high order. Numerical experiments provided in this paper have demonstrated that our results are truly promising, and consistent with our analysis and prediction.

A domain decomposition based spectral collocation method is proposed for numerically solving Lane-Emden equations, which are frequently encountered in mathematical physics and astrophysics. Compared with the existing methods, this method requires less computational cost and is particularly suitable for long-term computation. The related error estimates are also established, indicating the spectral accuracy of the method. The numerical performance and efficiency of the method are illustrated by several numerical experiments.

We study a novel class of numerical integrators, the adapted nested force-gradient schemes, used within the molecular dynamics step of the Hybrid Monte Carlo (HMC) algorithm. We test these methods in the Schwinger model on the lattice, a well known benchmark problem. We derive the analytical basis of nested force-gradient type methods and demonstrate the advantage of the proposed approach, namely reduced computational costs compared with other numerical integration schemes in HMC.

The Jiles-Atherton (J-A) model is a commonly used physics-based model in describing the hysteresis characteristics of ferromagnetic materials. However, citations and interpretation of this model in literature have been non-uniform. Solution methods for solving numerically this model has not been studied adequately. In this paper, through analyzing the mathematical properties of equations and the physical mechanism of energy conservation, we point out some unreasonable descriptions of this model and develop a relatively more accurate, modified J-A model together with its numerical solution method. Our method employs a fixed point method to compute anhysteretic magnetization. We obtain the susceptibility value of the anhysteretic magnetization analytically and apply the 4th order Runge-Kutta method to the solution of total magnetization. Computational errors are estimated and then precisions of the solving method in describing various materials are verified. At last, through analyzing the effects of the accelerating method, iterative error and step size on the computational errors, we optimize the numerical method to achieve the effects of high precision and short computing time. From analysis, we determine the range of best values of some key parameters for fast and accurate computation.

The differential quadrature method has been widely used in scientific and engineering computation. However, for the basic characteristics of time domain differential quadrature method, such as numerical stability and calculation accuracy or order, it is still lack of systematic analysis conclusions. In this paper, according to the principle of differential quadrature method, it has been derived and proved that the weighting coefficients matrix of differential quadrature method meets the important V-transformation feature. Through the equivalence of the differential quadrature method and the implicit Runge-Kutta method, it has been proved that the differential quadrature method is A-stable and s-stage s-order method. On this basis, in order to further improve the accuracy of the time domain differential quadrature method, a class of improved differential quadrature method of s-stage 2s-order has been proposed by using undetermined coefficients method and Padé approximations. The numerical results show that the proposed differential quadrature method is more precise than the traditional differential quadrature method.

In this paper, preconditioned iterative methods for solving two-dimensional space-fractional diffusion equations are considered. The fractional diffusion equation is discretized by a second-order finite difference scheme, namely, the Crank-Nicolson weighted and shifted Grünwald difference (CN-WSGD) scheme proposed in [W. Tian, H. Zhou and W. Deng, A class of second order difference approximation for solving space fractional diffusion equations, Math. Comp., 84 (2015) 1703-1727]. For the discretized linear systems, we first propose preconditioned iterative methods to solve them. Then we apply the D’Yakonov ADI scheme to split the linear systems and solve the obtained splitting systems by iterative methods. Two preconditioned iterative methods, the preconditioned generalized minimal residual (preconditioned GMRES) method and the preconditioned conjugate gradient normal residual (preconditioned CGNR) method, are proposed to solve relevant linear systems. By fully exploiting the structure of the coefficient matrix, we design two special kinds of preconditioners, which are easily constructed and are able to accelerate convergence of iterative solvers. Numerical results show the efficiency of our preconditioners.

This paper is concerned with the numerical stability of implicit Runge-Kutta methods for nonlinear neutral Volterra delay-integro-differential equations with constant delay. Using a Halanay inequality generalized by Liz and Trofimchuk, we give two sufficient conditions for the stability of the true solution to this class of equations. Runge-Kutta methods with compound quadrature rule are considered. Nonlinear stability conditions for the proposed methods are derived. As an illustration of the application of these investigations, the asymptotic stability of the presented methods for Volterra delay-integro-differential equations are proved under some weaker conditions than those in the literature. An extension of the stability results to such equations with weakly singular kernel is also discussed.