Error estimates to smooth solutions of semi‐discrete discontinuous Galerkin methods with quadrature rules for scalar conservation laws

In this article, we focus on error estimates to smooth solutions of semi‐discrete discontinuous Galerkin (DG) methods with quadrature rules for scalar conservation laws. The main techniques we use are energy estimate and Taylor expansion first introduced by Zhang and Shu in (Zhang and Shu, SIAM J Num Anal 42 (2004), 641–666). We show that, with (piecewise polynomials of degree k) finite elements in 1D problems, if the quadrature over elements is exact for polynomials of degree , error estimates of are obtained for general monotone fluxes, and optimal estimates of are obtained for upwind fluxes. For multidimensional problems, if in addition quadrature over edges is exact for polynomials of degree , error estimates of are obtained for general monotone fluxes, and are obtained for monotone and sufficiently smooth numerical fluxes. Numerical results validate our analysis. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 467–488, 2017     

A fast discontinuous finite element discretization for the space‐time fractional diffusion‐wave equation

In this article, we study fast discontinuous Galerkin finite element methods to solve a space‐time fractional diffusion‐wave equation. We introduce a piecewise‐constant discontinuous finite element method for solving this problem and derive optimal error estimates. Importantly, a fast solution technique to accelerate Toeplitz matrix‐vector multiplications which arise from discontinuous Galerkin finite element discretization is developed. This fast solution technique is based on fast Fourier transform and it depends on the special structure of coefficient matrices. In each temporal step, it helps to reduce the computational work from required by the traditional methods to log , where is the size of the coefficient matrices (number of spatial grid points). Moreover, the applicability and accuracy of the method are verified by numerical experiments including both continuous and discontinuous examples to support our theoretical analysis.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 2043–2061, 2017     

Decoupled schemes for unsteady MHD equations. I. time discretization

In this article, a time discretization decoupled scheme for two‐dimensional magnetohydrodynamics equations is proposed. The almost unconditional stability and convergence of this scheme are provided. The optimal error estimates for velocity and magnet are provided, and the optimal error estimate for pressure are deduced as well. Finite element spatial discretization and numerical implementation are considered in our article (Zhang and He, Comput Math Appl 69 (2015), 1390–1406). © 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 956–973, 2017     

Unconditional convergence and error estimates for bounded numerical solutions of the barotropic Navier–Stokes system

We consider a mixed finite‐volume finite‐element method applied to the Navier–Stokes system of equations describing the motion of a compressible, barotropic, viscous fluid. We show convergence as well as error estimates for the family of numerical solutions on condition that: (a) the underlying physical domain as well as the data are smooth; (b) the time step and the parameter of the spatial discretization are proportional, ; and (c) the family of numerical densities remains bounded for . No a priori smoothness is required for the limit (exact) solution. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1208–1223, 2017     

Fourth‐order compact and energy conservative scheme for solving nonlinear Klein‐Gordon equation

In this article, a fourth‐order compact and conservative scheme is proposed for solving the nonlinear Klein‐Gordon equation. The equation is discretized using the integral method with variational limit in space and the multidimensional extended Runge‐Kutta‐Nyström (ERKN) method in time. The conservation law of the space semidiscrete energy is proved. The proposed scheme is stable in the discrete maximum norm with respect to the initial value. The optimal convergent rate is obtained at the order of in the discrete ‐norm. Numerical results show that the integral method with variational limit gives an efficient fourth‐order compact scheme and has smaller error, higher convergence order and better energy conservation for solving the nonlinear Klein‐Gordon equation compared with other methods under the same condition. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1283–1304, 2017     

A fourth‐order Hermitian box‐scheme with fast solver for the Poisson problem in a cube

We present a fourth‐order Hermitian box‐scheme (HB‐scheme) for the Poisson problem in a cube. A single‐nonstaggered regular grid is used supporting the discrete unknowns u and . The scheme is fourth‐order accurate for u and in norm. The fast numerical resolution uses a matrix capacitance method, resulting in a computational complexity of . Numerical results are reported on several examples including nonseparable problems. The present scheme is the extension to the three‐dimensional case of the HB‐scheme presented in Abbas and Croisille [J Sci Comp 49 (2011), 239–267]. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 609–629, 2015     

Two unconditionally stable and convergent difference schemes with the extrapolation method for the one‐dimensional distributed‐order differential equations

The Grünwald formula is used to solve the one‐dimensional distributed‐order differential equations. Two difference schemes are derived. It is proved that the schemes are unconditionally stable and convergent with the convergence orders and in maximum norm, respectively, where and are step sizes in time, space and distributed order. The extrapolation method is applied to improve the approximate accuracy to the orders and respectively. An illustrative numerical example is given to confirm the theoretical results. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 591–615, 2016     

Error analysis of mixed finite element method for Poisson‐Nernst‐Planck system

To improve the convergence rate in L² norm from suboptimal to optimal for both electrostatic potential and ionic concentrations in Poisson‐Nernst‐Planck (PNP) system, we propose the mixed finite element method in this article to discretize the electrostatic potential equation, and still use the standard finite element method to discretize the time‐dependent ionic concentrations equations. Optimal error estimates in norm for the electrostatic potential, and in and norms for the ionic concentrations are attained. As a by‐product, the electric field can also achieve a higher approximation order in contrast with the standard finite element method for PNP system. Numerical experiments are performed to validate the theoretical results.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1924–1948, 2017     

Error and stability analysis of numerical solution for the time fractional nonlinear Schrödinger equation on scattered data of general‐shaped domains

In present work, a kind of spectral meshless radial point interpolation (SMRPI) technique is applied to the time fractional nonlinear Schrödinger equation in regular and irregular domains. The applied approach is based on erudite combination of meshless methods and spectral collocation techniques. The point interpolation method with the help of radial basis functions is used to construct shape functions which play as basis functions in the frame of SMRPI. It is proved the scheme is unconditionally stable with respect to the time variable in and also convergent by the order of convergence , . In the current work, the thin plate spline are used as the basis functions and to eliminate the nonlinearity, a simple predictor‐corrector (P‐C) scheme is performed. It is shown that the SMRPI solution, as a complex function, is suitable one for the time fractional nonlinear Schrödinger equation. The results of numerical experiments are compared to analytical solutions to confirm the reliable treatment of these stable solutions. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1043–1069, 2017     

On a nonoverlapping additive Schwarz method for h‐p discontinuous Galerkin discretization of elliptic problems

The condition number of a discontinuous Galerkin finite element discretization preconditioned with a nonoverlapping additive Schwarz method is analyzed. We improve the result of Antonietti and Houston (J Sci Comput 46 (2011), 124–149), where a bound has been proved for a two‐level nonoverlapping additive Schwarz method with coarse problem using polynomials of degree on a coarse mesh size . In a more general framework, where the concurrency of the algorithm is increased by applying solvers on subdomains smaller than the coarse grid cells, we prove that the condition number of the preconditioned system is where is the coarse space element degree polynomial and is the size of subdomain where local problems are solved in parallel. Our result also extends to the case of discontinuous coefficient, piecewise constant on the coarse grid, for a composite continuous–discontinuous Galerkin discretization. © 2016Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1572–1590, 2016