Raising the order of geometric numerical integrators by composition and extrapolation

We analyse composition and polynomial extrapolation as procedures to raise the order of a geometric integrator for solving numerically differential equations. Methods up to order sixteen are constructed starting with basic symmetric schemes of order six and eight. If these are geometric integrators, then the new methods obtained by extrapolation preserve the geometric properties up to a higher order than the order of the method itself. We show that, for a number of problems, this is a very efficient procedure to obtain high accuracy. The relative performance of the different algorithms is examined on several numerical experiments. AMS subject classification 17B66, 34A50, 65L05     

Symmetric and symplectic exponential integrators for nonlinear Hamiltonian systems

This letter studies symmetric and symplectic exponential integrators when applied to numerically computing nonlinear Hamiltonian systems. We first establish the symmetry and symplecticity conditions of exponential integrators and then show that these conditions are extensions of the symmetry and symplecticity conditions of Runge–Kutta methods. Based on these conditions, some symmetric and symplectic exponential integrators up to order four are derived. Two numerical experiments are carried out and the results demonstrate the remarkable numerical behaviour of the new exponential integrators in comparison with some symmetric and symplectic Runge–Kutta methods in the literature.     

A test problem for molecular dynamics integrators

We derive a test problem for evaluating the ability of time-steppingmethods to preserve the statistical properties of systems inmolecular dynamics. We consider a family of deterministic systemsconsisting of a finite number of particles interacting on acompact interval. The particles are given random initial conditionsand interact through instantaneous energy- and momentum-conservingcollisions. As the number of particles, the particle density,and the mean particle speed go to infinity, the trajectory ofa tracer particle is shown to converge to a stationary Gaussianstochastic process. We approximate this system by one describedby a system of ordinary differential equations and provide numericalevidence that it converges to the same stochastic process. Wesimulate the latter system with a variety of numerical integrators,including the symplectic Euler method, a fourth-order Runge-Kuttamethod, and an energyconserving step-and-project method. Weassess the methods' ability to recapture the system's limitingstatistics and observe that symplectic Euler performs significantlybetter than the others for comparable computational expense.     

A posteriori error estimation for a defect correction method applied to conduction convection problems

We present a posteriori error estimate for a defect correction method for approximating solutions of the stationary conduction convection problems in two dimension. The defect correction method is aiming at small viscosity ν. A reliable a posteriori error estimation is derived for the defect correction method. Finally, two numerical examples validate our theoretical results. The first example is a problem with known solution and the second example is a physical model of square cavity stationary flow. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

Solving parabolic problems using local defect correction in combination with the finite volume method

We present a method for solving partial differential equations characterized by highly localized properties in which the local defect correction (LDC) algorithm for time‐dependent problems is combined with a finite volume discretization. At each time step, LDC computes a numerical solution on a composite grid, a union of a global uniform coarse grid and a local uniform fine grid. The main feature of the method is that the discrete conservation property, typical of the finite volume approach is preserved on the composite grid. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

The application of iterated defect correction to theLod method for parabolic equations

The paper is concerned with the numerical solution of the initial boundary value problem for a class of multi-dimensional parabolic partial differential equations. In particular the time-integration of semi-discrete equations is investigated. An attempt is made to develop integration formulas being computationally attractive and of high accuracy, while possessing unconditional stability properties. To that end iterated defect correction is applied to the LOD method. The convergence properties of this process are investigated. Numerical experiments are reported.     

Conservation properties of numerical integrators for highly oscillatory Hamiltonian systems

^** Email: David.Cohen{at}math.unige.ch. Present address: Mathematisches Institut, Universität Tübingen, D-72076 Tübingen, Germany (cohen{at}na.uni-tuebingen.de) Modulated Fourier expansion is used to show long-time near-conservationof the total and oscillatory energies of numerical methods forHamiltonian systems with highly oscillatory solutions. The numericalmethods considered are an extension of the trigonometric methods.A brief discussion of conservation properties in the continuousproblem and in the multi-frequency case is also given.     

Analysis of local defect correction and high‐order compact finite differences

We study the possibility of combining the LDC technique with high‐order compact schemes. An algorithm is shown first for the 1D stationary convection‐diffusion equation, and then it is extended to 2D. The results of testing show that we get the same accuracy of the solution as on the reference fine grid with much less points in the domain (up to 50% fewer points for the examples presented here). © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2006     

Error estimate of a fully discrete defect correction finite element method for unsteady incompressible Magnetohydrodynamics equations

In this study, a fully discrete defect correction finite element method for the unsteady incompressible Magnetohydrodynamics (MHD) equations, which is leaded by combining the Back Euler time discretization with the two-step defect correction in space, is presented. It is a continuous work of our formal paper [Math Method Appl Sci. 2017. DOI:10.1002/mma.4296]. The defect correction method is an iterative improvement technique for increasing the accuracy of a numerical solution without applying a grid refinement. Firstly, the nonlinear MHD equation is solved with an artificial viscosity term. Then, the numerical solutions are improved on the same grid by a linearized defect-correction technique. Then, we introduce the numerical analysis including stability analysis and error analysis. The numerical analysis proves that our method is stable and has an optimal convergence rate. Some numerical results [see Math Method Appl Sci. 2017. DOI:10.1002/mma.4296] show that this method is highly efficient for the unsteady incompressible MHD problems.     

Modified defect correction algorithms for ODEs. Part II: Stiff initial value problems

As shown in part I of this paper and references therein, the classical method of Iterated Defect Correction (IDeC) can be modified in several nontrivial ways, extending the flexibility and range of applications of this approach. The essential point is an adequate definition of the defect, resulting in a significantly more robust convergence behavior of the IDeC iteration, in particular, for nonequidistant grids. The present part II is devoted to the efficient high-order integration of stiff initial value problems. By means of model problem investigation and systematic numerical experiments with a set of stiff test problems, our new versions of defect correction are systematically evaluated, and further algorithmic measures are proposed for the stiff case. The performance of the different variants under consideration is compared, and it is shown how strong coupling between non-stiff and stiff components can be successfully handled. AMS subject classification 65L05 Supported by the Austrian Research Fund (FWF) grant P-15030.     

A defect correction method for the time‐dependent Navier‐Stokes equations

A method for solving the time dependent Navier‐Stokes equations, aiming at higher Reynolds' number, is presented. The direct numerical simulation of flows with high Reynolds' number is computationally expensive. The method presented is unconditionally stable, computationally cheap, and gives an accurate approximation to the quantities sought. In the defect step, the artificial viscosity parameter is added to the inverse Reynolds number as a stability factor, and the system is antidiffused in the correction step. Stability of the method is proven, and the error estimations for velocity and pressure are derived for the one‐ and two‐step defect‐correction methods. The spacial error is O(h) for the one‐step defect‐correction method, and O(h²) for the two‐step method, where h is the diameter of the mesh. The method is compared to an alternative approach, and both methods are applied to a singularly perturbed convection–diffusion problem. The numerical results are given, which demonstrate the advantage (stability, no oscillations) of the method presented. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

SBA splitting for approximation of invariants in time-dependent Hamiltonian systems

Most physical phenomena are described by time-dependent Hamiltonian systems with qualitative features that should be preserved by numerical integrators used for approximating their dynamics. The initial energy of the system together with the energy added or subtracted by the outside forces, represent a conserved quantity of the motion. For a class of time-dependent Hamiltonian systems [8] this invariant can be defined by means of an auxiliary function whose dynamics has to be integrated simultaneously with the system’s equations. We propose splitting procedures featured by a SB³A property that allows to construct composition methods with a reduced number of determining order equations and to provide the same high accuracy for both the dynamics and the preservation of the invariant quantity.     

A local defect correction technique for time‐dependent problems

In this article a local defect correction technique for time‐dependent problems is presented. The method is suitable for solving partial differential equations characterized by a high activity, which is mainly located, at each time, in a small part of the physical domain. The problem is solved at each time step by means of a global uniform coarse grid and a local uniform fine grid. Local and global approximation are improved iteratively. Results of numerical experiments illustrate the accuracy, the efficiency, and the robustness of the method. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Defect correction method for viscoelastic fluid flows at high Weissenberg number

We study a defect correction method for the approximation of viscoelastic fluid flow. In the defect step, the constitutive equation is computed with an artificially reduced Weissenberg parameter for stability, and the resulting residual is corrected in the correction step. We prove the convergence of the defect correction method and derive an error estimate for the Oseen‐viscoelastic model problem. The derived theoretical results are supported by numerical tests for both the Oseen‐viscoelastic problem and the Johnson‐Segalman model problem. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Boosting the accuracy of finite difference schemes via optimal time step selection and non-iterative defect correction

In this article, we present a unified analysis of the simple technique for boosting the order of accuracy of finite difference schemes for time dependent partial differential equations (PDEs) by optimally selecting the time step used to advance the numerical solution and adding defect correction terms in a non-iterative manner. The power of the technique, which is applicable to time dependent, semilinear, scalar PDEs where the leading-order spatial derivative has a constant coefficient, is its ability to increase the accuracy of formally low-order finite difference schemes without major modification to the basic numerical algorithm. Through straightforward numerical analysis arguments, we explain the origin of the boost in accuracy and estimate the computational cost of the resulting numerical method. We demonstrate the utility of optimal time step (OTS) selection combined with non-iterative defect correction (NIDC) on several different types of finite difference schemes for a wide array of classical linear and semilinear PDEs in one and more space dimensions on both regular and irregular domains.     

Convergence Analysis for the Iterated Defect Correction Scheme of Finite Element Methods on Rectangle Grids

This paper develops a new method to analyze convergence of the iterated defect correction scheme of finite element methods on rectangular grids in both two and three dimensions. The main idea is to formulate energy inner products and energy (semi)norms into matrix forms. Then, two constants of two key inequalities involved are min and max eigenvalues of two associated generalized eigenvalue problems, respectively. Local versions on the element level of these two generalized eigenvalue problems are exactly solved to obtain sharp (lower) upper bounds of these two constants. This and some essential observations for iterated solutions establish convergence in 2D and the monotone decreasing property in 3D. For two dimensions the results herein improve those in literature; for three dimensions the results herein are new. Numerical results are presented to examine theoretical results.     

Defect correction and a posteriori error estimation of Petrov-Galerkin methods for nonlinear Volterra integro-differential equations

We present two defect correction schemes to accelerate the Petrov-Galerkin finite element methods [19] for nonlinear Volterra integro-differential equations. Using asymptotic expansions of the errors, we show that the defect correction schemes can yield higher order approximations to either the exact solution or its derivative. One of these schemes even does not impose any extra regularity requirement on the exact solution. As by-products, all of these higher order numerical methods can also be used to form a posteriori error estimators for accessing actual errors of the Petrov-Galerkin finite element solutions. Numerical examples are also provided to illustrate the theoretical results obtained in this paper.     

A new defect correction method for the Navier-Stokes equations at high Reynolds numbers

In this paper, a new defect correction method for the Navier-Stokes equations is presented. With solving an artificial viscosity stabilized nonlinear problem in the defect step, and correcting the residual by linearized equations in the correction step for a few steps, this combination is particularly efficient for the Navier-Stokes equations at high Reynolds numbers. In both the defect and correction steps, we use the Oseen iterative scheme to solve the discrete nonlinear equations. Furthermore, the stability and convergence of this new method are deduced, which are better than that of the classical ones. Finally, some numerical experiments are performed to verify the theoretical predictions and show the efficiency of the new combination.