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.

We propose a hybrid spectral element method for fractional two-point boundary value problem (FBVPs) involving both Caputo and Riemann-Liouville (RL) fractional derivatives. We first formulate these FBVPs as a second kind Volterra integral equation (VIEs) with weakly singular kernel, following a similar procedure in [16]. We then design a hybrid spectral element method with generalized Jacobi functions and Legendre polynomials as basis functions. The use of generalized Jacobi functions allow us to deal with the usual singularity of solutions at t = 0. We establish the existence and uniqueness of the numerical solution, and derive a hptype error estimates under L2(I)-norm for the transformed VIEs. Numerical results are provided to show the effectiveness of the proposed methods.

We propose a class of numerical methods for solving nonlinear random differential equations with piecewise constant argument, called gPCRK methods as they combine generalised polynomial chaos with Runge-Kutta methods. An error analysis is presented involving the error arising from a finite-dimensional noise assumption, the projection error, the aliasing error and the discretisation error. A numerical example is given to illustrate the effectiveness of this approach.

We introduce a multiple interval Chebyshev-Gauss-Lobatto spectral collocation method for the initial value problems of the nonlinear ordinary differential equations (ODES). This method is easy to implement and possesses the high order accuracy. In addition, it is very stable and suitable for long time calculations. We also obtain the hp-version bound on the numerical error of the multiple interval collocation method under H 1-norm. Numerical experiments confirm the theoretical expectations.

The h-p version of the continuous Petrov-Galerkin time stepping method is analyzed for nonlinear initial value problems. An L∞-error bound explicit with respect to the local discretization and regularity parameters is derived. Numerical examples are provided to illustrate the theoretical results.

We propose a new Adomian decomposition method (ADM) using an integrating factor for the Emden–Fowler equation. With this method, we are able to solve certain Emden–Fowler equations for which the traditional ADM fails. Numerical results obtained from testing our linear and nonlinear models are far more reliable and efficient than those from existing methods. We also present a complete error analysis and a convergence criterion for this method. One drawback of the traditional ADM is that the interval of convergence of the Adomian truncated series is very small. Some techniques, such as Pade approximants, can enlarge this interval, but they are too complicated. Here, we use a continuation technique to extend our method to a larger interval.

In this paper, we present and analyze a single interval Legendre-Gauss spectral collocation method for solving the second order nonlinear delay differential equations with variable delays. We also propose a novel algorithm for the single interval scheme and apply it to the multiple interval scheme for more efficient implementation. Numerical examples are provided to illustrate the high accuracy of the proposed methods.

First introduced in, the lumped particle framework is a flexible and numerically efficient framework for the modelling of particle transport in fluid flow. In this paper, the framework is expanded to simulate multicomponent particle-laden fluid flow. This is accomplished by introducing simulation protocols to model particles over a wide range of length and time scales. Consequently, we present a time ordering scheme and an approximate approach for accelerating the computation of evolution of different particle constituents with large differences in physical scales. We apply the extended framework on the temporal evolution of three particle constituents in sandladen flow, and horizontal release of spherical particles. Furthermore, we evaluate the numerical error of the lumped particle model. In this context, we discuss the Velocity-Verlet numerical scheme, and show how to apply this to solving Newton’s equations within the framework. We show that the increased accuracy of the Velocity-Verlet scheme is not lost when applied to the lumped particle framework.

We consider an optimal control problem with an 1D singularly perturbed differential state equation. For solving such problems one uses the enhanced system of the state equation and its adjoint form. Thus, we obtain a system of two convection-diffusion equations. Using linear finite elements on adapted grids we treat the effects of two layers arising at different boundaries of the domain. We proof uniform error estimates for this method on meshes of Shishkin type. We present numerical results supporting our analysis.