A note on high order schemes for the one dimensional wave equation

We construct in two ways stable schemes of arbitrarily high order for one dimensional wave equation. In both cases the stability limit is independent of the order. The other construction yields also automatically stable schemes for the heterogeneous case. In this case the estimated stability limit decreases slightly with order.     

Compact finite difference schemes of sixth order for the Helmholtz equation

New compact finite difference schemes of sixth order are derived for the three dimensional Helmholtz equation, Δu-κ²u=-f$\mathrm{\Delta }u-{\kappa }^{2}u=-f$. Convergence characteristics and accuracy are compared and a truncation error analysis is presented for a broad range of κ$\kappa$-values.     

Parameterized families of finite difference schemes for the wave equation

This article is devoted to an analysis of simple families of finite difference schemes for the wave equation. These families are dependent on several free parameters, and methods for obtaining stability bounds as a function of these parameters are discussed in detail. Access to explicit stability bounds such as those derived here may, it is hoped, lead to optimization techniques for so‐called spectral‐like methods, which are difference schemes dependent on many free parameters (and for which maximizing the order of accuracy may not be the defining criterion). Though the focus is on schemes for the wave equation in one dimension, the analysis techniques are extended to two dimensions; implicit schemes such as ADI methods are examined in detail. Numerical results are presented. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 20: 463–480, 2004.     

Geometric finite difference schemes for the generalized hyperelastic-rod wave equation

Geometric integrators are presented for a class of nonlinear dispersive equations which includes the Camassa-Holm equation, the BBM equation and the hyperelastic-rod wave equation. One group of schemes is designed to preserve a global property of the equations: the conservation of energy; while the other one preserves a more local feature of the equations: the multi-symplecticity.     

Finite difference schemes for the two‐dimensional semilinear wave equation

In this article, two finite difference schemes for solving the semilinear wave equation are proposed. The unique solvability and the stability are discussed. The second‐order accuracy convergence in both time and space in the discrete H¹‐norm for the two proposed difference schemes is proved. Numerical experiments are performed to support our theoretical results.     

New estimates of mixed finite element method for fourth‐order wave equation

The conforming bilinear mixed finite element method is established for a fourth‐order wave equation. By applying the interpolation and projection simultaneously, the superclose and superconvergence results of order O (h ²) for the original variable u and intermediate variable p  =? a ₁(X )Δu in H¹‐norm are achieved for the semi‐discrete and fully discrete schemes, respectively (where a ₁(X ) is the coefficient appeared in the equation). The distinct advantage of this method is that it can reduce the smoothness of solutions u ,u _t , and p compared with the existing literature for getting optimal estimates alone. At last, some numerical results are carried out to validate the theoretical analysis. Copyright © 2017 John Wiley & Sons, Ltd.     

Locally one-dimensional difference schemes for the fractional order diffusion equation

Locally-one-dimensional difference schemes for the fractional diffusion equation in multidimensional domains are considered. Stability and convergence of locally one-dimensional schemes for this equation are proved.     

The maximal accuracy of stable difference schemes for the wave equation

We consider three time-level difference schemes, symmetric in time and space, for the solution of the wave equation,u _tt =c ² u _xx , given by

High‐order finite difference schemes for numerical solutions of the generalized Burgers–Huxley equation

In this article, up to tenth‐order finite difference schemes are proposed to solve the generalized Burgers–Huxley equation. The schemes based on high‐order differences are presented using Taylor series expansion. To establish the numerical solutions of the corresponding equation, the high‐order schemes in space and a fourth‐order Runge‐Kutta scheme in time have been combined. Numerical experiments have been conducted to demonstrate the high‐order accuracy of the current algorithms with relatively minimal computational effort. The results showed that use of the present approaches in the simulation is very applicable for the solution of the generalized Burgers–Huxley equation. The current results are also seen to be more accurate than some results given in the literature. The proposed algorithms are seen to be very good alternatives to existing approaches for such physical applications. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1313‐1326, 2011     

Asymptotically exact a posteriori local discontinuous Galerkin error estimates for the one‐dimensional second‐order wave equation

In this article, we analyze a residual‐based a posteriori error estimates of the spatial errors for the semidiscrete local discontinuous Galerkin (LDG) method applied to the one‐dimensional second‐order wave equation. These error estimates are computationally simple and are obtained by solving a local steady problem with no boundary condition on each element. We apply the optimal L² error estimates and the superconvergence results of Part I of this work [Baccouch, Numer Methods Partial Differential Equations 30 (2014), 862–901] to prove that, for smooth solutions, these a posteriori LDG error estimates for the solution and its spatial derivative, at a fixed time, converge to the true spatial errors in the L²‐norm under mesh refinement. The order of convergence is proved to be , when p‐degree piecewise polynomials with are used. As a consequence, we prove that the LDG method combined with the a posteriori error estimation procedure yields both accurate error estimates and superconvergent solutions. Our computational results show higher convergence rate. We further prove that the global effectivity indices, for both the solution and its derivative, in the L²‐norm converge to unity at rate while numerically they exhibit and rates, respectively. Numerical experiments are shown to validate the theoretical results. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1461–1491, 2015     

Conservative arbitrary order finite difference schemes for shallow-water flows

The classical nonlinear shallow-water model (SWM) of an ideal fluid is considered. For the model, a new method for the construction of mass and total energy conserving finite difference schemes is suggested. In fact, it produces an infinite family of finite difference schemes, which are either linear or nonlinear depending on the choice of certain parameters. The developed schemes can be applied in a variety of domains on the plane and on the sphere. The method essentially involves splitting of the model operator by geometric coordinates and by physical processes, which provides substantial benefits in the computational cost of solution. Besides, in case of the whole sphere it allows applying the same algorithms as in a doubly periodic domain on the plane and constructing finite difference schemes of arbitrary approximation order in space. Results of numerical experiments illustrate the skillfulness of the schemes in describing the shallow-water dynamics.     

Application of transparent boundary conditions to high-order finite-difference schemes for the wave equation in waveguides

We propose a method for generating finite-difference approximations of transparent boundary conditions (TBCs) with the fourth and sixth order in space. It is based on the wave equation solution continuation extra two or three layers of grid points outside the computational domain to use them in central-difference operators on approaching the boundary. We present the theoretical background of the method, give estimates of computational resources, and discuss accuracy and stability results of numerical tests in 1D and 2D cases.     

Superconvergence of the local discontinuous galerkin method applied to the one‐dimensional second‐order wave equation

We analyze the superconvergence properties of the local discontinuous Galerkin (LDG) method applied to the second‐order wave equation in one space dimension. With a suitable projection of the initial conditions for the LDG scheme, we prove that the LDG solution and its spatial derivative are super close to particular projections of the exact solutions for pth‐degree polynomial spaces. We use these results to show that the significant parts of the discretization errors for the LDG solution and its derivative are proportional to ‐degree right and left Radau polynomials, respectively. These results allow us to prove that the p‐degree LDG solution and its derivative are superconvergent at the roots of ‐degree right and left Radau polynomials, respectively, while computational results show higher convergence rate. Superconvergence results can be used to construct asymptotically correct a posteriori error estimates by solving a local steady problem on each element. This will be discussed further in Part II of this work, where we will prove that the a posteriori LDG error estimates for the solution and its derivative converge to the true spatial errors in the L ²‐norm under mesh refinement. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 862–901, 2014     

Some temporal second order difference schemes for fractional wave equations

In this article, motivated by Alikhanov's new work (Alikhanov, J Comput Phys 280 (2015), 424–438), some difference schemes are proposed for both one‐dimensional and two‐dimensional time‐fractional wave equations. The obtained schemes can achieve second‐order numerical accuracy both in time and in space. The unconditional convergence and stability of these schemes in the discrete H¹‐norm are proved by the discrete energy method. The spatial compact difference schemes with the results on the convergence and stability are also presented. In addition, the three‐dimensional problem is briefly mentioned. Numerical examples illustrate the efficiency of the proposed schemes. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 970–1001, 2016     

The combination of collocation,finite difference,and multigrid methods for solution of the two‐dimensional wave equation

In this article, we apply compact finite difference approximations of orders two and four for discretizing spatial derivatives of wave equation and collocation method for the time component. The resulting method is unconditionally stable and solves the wave equation with high accuracy. The solution is approximated by a polynomial at each grid point that its coefficients are determined by solving a linear system of equations. We employ the multigrid method for solving the resulted linear system. Multigrid method is an iterative method which has grid independently convergence and solves the linear system of equations in small amount of computer time. Numerical results show that the compact finite difference approximation of fourth order, collocation and multigrid methods produce a very efficient method for solving the wave equation. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

High‐order compact difference schemes for the modified anomalous subdiffusion equation

In this article, two kinds of high‐order compact finite difference schemes for second‐order derivative are developed. Then a second‐order numerical scheme for a Riemann–Liouvile derivative is established based on a fractional centered difference operator. We apply these methods to a fractional anomalous subdiffusion equation to construct two kinds of novel numerical schemes. The solvability, stability, and convergence analysis of these difference schemes are studied by using Fourier method. The convergence orders of these numerical schemes are and , respectively. Finally, numerical experiments are displayed which are in line with the theoretical analysis. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 213–242, 2016     

Higher order finite difference schemes for the magnetic induction equations

We describe high order accurate and stable finite difference schemes for the initial-boundary value problem associated with the magnetic induction equations. These equations model the evolution of a magnetic field due to a given velocity field. The finite difference schemes are based on Summation by Parts (SBP) operators for spatial derivatives and a Simultaneous Approximation Term (SAT) technique for imposing boundary conditions. We present various numerical experiments that demonstrate both the stability as well as high order of accuracy of the schemes.      

On L1 decay problem for the dissipative wave equation

We study the decay estimates of solutions to the Cauchy problem for the dissipative wave equation in one, two, and three dimensions. The representation formulas of the solutions provide the sharp decay rates on L¹ norms and also L^p norms. Copyright © 2003 John Wiley & Sons, Ltd.     

Homotopy perturbation method for two dimensional time-fractional wave equation

The aim of this paper is to present an efficient and reliable treatment of the homotopy perturbation method (HPM) for two dimensional time-fractional wave equation (TFWE) with the boundary conditions. The fractional derivative is described in the Caputo sense. The initial approximation can be determined by imposing the boundary conditions. The method provides approximate solutions in the form of convergent series with easily computable components. The obtained results shown that the technique introduced here is efficient and easy to implement.     

Modification of basis functions in high order discontinuous Galerkin schemes for advection equation

High order discontinuous Galerkin (DG) discretization schemes are considered for an advection boundary-value problem on 2-D unstructured grids with arbitrary geometry of grid cells. A number of test cases are developed to study the sensitivity of a high order DG scheme to local grid distortion. It will be demonstrated how to modify the formulation of a DG discretization for the advection equation. Our approach allows one to maintain the required accuracy on distorted grids while using a fewer number of basis functions for the solution approximation in order to save computational resources.