High accuracy error estimates of a Galerkin finite element method for nonlinear time fractional diffusion equation

In this work, an effective and fast finite element numerical method with high-order accuracy is discussed for solving a nonlinear time fractional diffusion equation. A two-level linearized finite element scheme is constructed and a temporal–spatial error splitting argument is established to split the error into two parts, that is, the temporal error and the spatial error. Based on the regularity of the time discrete system, the temporal error estimate is derived. Using the property of the Ritz projection operator, the spatial error is deduced. Unconditional superclose result in H¹-norm is obtained, with no additional regularity assumption about the exact solution of the problem considered. Then the global superconvergence error estimate is obtained through the interpolated postprocessing technique. In order to reduce storage and computation time, a fast finite element method evaluation scheme for solving the nonlinear time fractional diffusion equation is developed. To confirm the theoretical error analysis, some numerical results are provided.     

Convergence analysis for second‐order accurate schemes for the periodic nonlocal Allen‐Cahn and Cahn‐Hilliard equations

In this paper, we provide a detailed convergence analysis for fully discrete second‐order (in both time and space) numerical schemes for nonlocal Allen‐Cahn and nonlocal Cahn‐Hilliard equations. The unconditional unique solvability and energy stability ensures ? ⁴ stability. The convergence analysis for the nonlocal Allen‐Cahn equation follows the standard procedure of consistency and stability estimate for the numerical error function. For the nonlocal Cahn‐Hilliard equation, because of the complicated form of the nonlinear term, a careful expansion of its discrete gradient is undertaken, and an H ^?1 inner‐product estimate of this nonlinear numerical error is derived to establish convergence. In addition, an a priori bound of the numerical solution at the discrete level is needed in the error estimate. Such a bound can be obtained by performing a higher order consistency analysis by using asymptotic expansions for the numerical solution. Following the technique originally proposed by Strang (eg, 1964), instead of the standard comparison between the exact and numerical solutions, an error estimate between the numerical solution and the constructed approximate solution yields an O (s ³+h ⁴) convergence in norm, in which s and h denote the time step and spatial mesh sizes, respectively. This in turn leads to the necessary bound under a standard constraint s ≤C h . Here, we also prove convergence of the scheme in the maximum norm under the same constraint.     

Finite difference discretization of the Rosenau-RLW equation

In this paper, numerical solutions of the Rosenau-RLW equation are considered using Crank–Nicolson type finite difference method. Existence of the numerical solutions is derived by Brouwer fixed point theorem. A priori bound and the error estimates as well as conservation of discrete mass and discrete energy for the finite difference solutions are discussed. Unconditional stability, second-order convergence and uniqueness of the scheme are proved using discrete energy method. Some numerical experiments have been conducted in order to verify the theoretical results.     

A fully Galerkin method for the damped generalized regularized long‐wave (DGRLW) equation

In this article, a fully discrete Galerkin scheme based on a nonlinear Crank–Nicolson method to approximate the solution of the DGRLW equation is constructed. Some a priori bounds are proved as well as error estimates. Then, a linearized modification scheme by an extrapolation method is discussed. The two schemes are time second order convergence. The last part is devoted to some numerical results. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Cahn-Hilliard方程非条件能量稳定二阶方法的稳定性和收敛性分析

在这项工作中,我们构建了一种非条件能量稳定的高效不变能量四分法(IEQ)来求解Cahn-Hilliard方程.所构建的数值格式是线性的、具有二阶时间精度和非条件能量稳定性.我们仔细分析了数值格式的唯一可解性、稳定性和误差估计.结果表明,所构建的格式满足唯一可解性、非条件能量稳定性和时间方向的二阶收敛性.通过大量的二维和三维数值实验,我们进一步验证了所提出格式的收敛阶、非条件能量稳定性和有效性.     

On the convergence of a conservative numerical scheme for the usual Rosenau-RLW equation

In this paper, we study the initial-boundary value problem of the usual Rosenau-RLW equation by finite difference method. We design a conservative numerical scheme which preserves the original conservative properties for the equation. The scheme is three-level and linear-implicit. The unique solvability of numerical solutions has been shown. Priori estimate and second order convergence of the finite difference approximate solutions are discussed by discrete energy method. Numerical results demonstrate that the scheme is efficient and accurate.     

Recovering nonlinear terms in an inverse boundary value problem for Laplace's equation: A stability estimate

Stationary thermography can be used for investigating the functional form of a nonlinear cooling law that describes heat exchanges through an inaccessible part of the boundary of a conductor. In this paper, we obtain a logarithmic stability estimate for the associated nonlinear inverse problem. This stability estimate is obtained from the convergence and sensitivity analysis of a finite difference method for the numerical solution of the Cauchy problem for Laplace's equation, based on the Störmer–Verlet scheme.     

抛物方程时间连续有限元全离散格式的超收敛性

利用张量积分解和时间方向单元正交分解,证明了线性抛物型方程的时间连续全离散有限元在单元节点和内部的特征点的超收敛性.并用连续有限元计算了非线性Schrodinger方程,验证了能量的守恒性.计算结果与理论相吻合.     

二维半线性反应扩散方程的交替方向隐格式

本文研究一类二维半线性反应扩散方程的差分方法.构造了一个二层线性化交替方向隐格式.利用离散能量估计方法证明了差分格式解的存在唯一性、差分格式在离散H~1模下的二阶收敛性和稳定性.最后给出两个数值例子验证了理论分析结果.     

A fast numerical solution of the 3D nonlinear tempered fractional integrodifferential equation

In this paper, we investigate the numerical solution of the three-dimensional (3D) nonlinear tempered fractional integrodifferential equation which is subject to the initial and boundary conditions. The backward Euler (BE) method in association with the first-order convolution quadrature rule is employed to discretize this equation for time, and the Galerkin finite element method is applied for space, which is combined with an alternating direction implicit (ADI) algorithm, in order to reduce the computational cost for solving the three-dimensional nonlocal problem. Then a fully discrete BE ADI Galerkin finite element scheme can be obtained by linearizing the non-linear term. Thereafter we prove a positive-type lemma, from which the stability and convergence of the proposed numerical scheme are derived based on the energy method. Numerical experiments are performed to verify the effectiveness of the proposed approach.     

Convergence of a finite volume scheme for coagulation-fragmentation equations

This paper is devoted to the analysis of a numerical scheme for the coagulation and fragmentation equation. A time explicit finite volume scheme is developed, based on a conservative formulation of the equation. It is shown to converge under a stability condition on the time step, while a first order rate of convergence is established and an explicit error estimate is given. Finally, several numerical simulations are performed to investigate the gelation phenomenon and the long time behavior of the solution.

     

On the convergence of operator splitting for the Rosenau–Burgers equation

We present convergence analysis of operator splitting methods applied to the nonlinear Rosenau–Burgers equation. The equation is first splitted into an unbounded linear part and a bounded nonlinear part and then operator splitting methods of Lie‐Trotter and Strang type are applied to the equation. The local error bounds are obtained by using an approach based on the differential theory of operators in Banach space and error terms of one and two‐dimensional numerical quadratures via Lie commutator bounds. The global error estimates are obtained via a Lady Windermere's fan argument. Lastly, a numerical example is studied to confirm the expected convergence order.     

Numerical and theoretical discussions for solving nonlinear generalized Benjamin–Bona–Mahony–Burgers equation based on the Legendre spectral element method

This article is devoted to solving numerically the nonlinear generalized Benjamin–Bona–Mahony–Burgers (GBBMB) equation that has several applications in physics and applied sciences. First, the time derivative is approximated by using a finite difference formula. Afterward, the stability and convergence analyses of the obtained time semi‐discrete are proven by applying the energy method. Also, it has been demonstrated that the convergence order in the temporal direction is O(dt) . Second, a fully discrete formula is acquired by approximating the spatial derivatives via Legendre spectral element method. This method uses Lagrange polynomial based on Gauss–Legendre–Lobatto points. An error estimation is also given in detail for full discretization scheme. Ultimately, the GBBMB equation in the one‐ and two‐dimension is solved by using the proposed method. Also, the calculated solutions are compared with theoretical solutions and results obtained from other techniques in the literature. The accuracy and efficiency of the mentioned procedure are revealed by numerical samples.     

Multilevel augmentation methods for solving the sine-Gordon equation

In this paper we develop the multilevel augmentation method for solving nonlinear operator equations of the second kind and apply it to solving the one-dimensional sine-Gordon equation. We first give a general setting of the multilevel augmentation method for solving the second kind nonlinear operator equations and prove that the multilevel augmentation method preserves the optimal convergence order of the projection method while reducing computational cost significantly. Then we describe the semi-discrete scheme and the fully-discrete scheme based on multiscale methods for solving the sine-Gordon equation, and apply the multilevel augmentation method to solving the discrete equation. A complete analysis for convergence order is proposed. Finally numerical experiments are presented to confirm the theoretical results and illustrate the efficiency of the method.     

Fourth‐order compact scheme for the one‐dimensional sine‐Gordon equation

Finite difference scheme to the generalized one‐dimensional sine‐Gordon equation is considered in this paper. After approximating the second order derivative in the space variable by the compact finite difference, we transform the sine‐Gordon equation into an initial‐value problem of a second‐order ordinary differential equation. Then Padé approximant is used to approximate the time derivatives. The resulting fully discrete nonlinear finite‐difference equation is solved by a predictor‐corrector scheme. Both Dirichlet and Neumann boundary conditions are considered in our proposed algorithm. Stability analysis and error estimate are given for homogeneous Dirichlet boundary value problems using energy method. Numerical results are given to verify the condition for stability and convergence and to examine the accuracy and efficiency of the proposed algorithm. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

带强耗散项的拟线性波方程的数值解

研究了一类拟线性波方程的数值解．构造了带强耗散项的拟线性波方程的三级差分格式,并证明其收敛性,估计了差分解的误差．最后给出数值例子．     

一维非线性Schrödinger 方程的两个无条件收敛的守恒紧致差分格式

本文对一维非线性 Schrödinger 方程给出两个紧致差分格式, 运用能量方法和两个新的分析技 巧证明格式关于离散质量和离散能量守恒, 而且在最大模意义下无条件收敛. 对非线性紧格式构造了 一个新的迭代算法, 证明了算法的收敛性, 并在此基础上给出一个新的线性化紧格式. 数值算例验证 了理论分析的正确性, 并通过外推进一步提高了数值解的精度.     

Conservative compact finite difference scheme for the N‐coupled nonlinear Klein–Gordon equations

In this article, a compact finite difference method is developed for the periodic initial value problem of the N‐coupled nonlinear Klein–Gordon equations. The present scheme is proved to preserve the total energy in the discrete sense. Due to the difficulty in obtaining the priori estimate from the discrete energy conservation law, the cut‐off function technique is employed to prove the convergence, which shows the new scheme possesses second order accuracy in time and fourth order accuracy in space, respectively. Additionally, several numerical results are reported to confirm our theoretical analysis. Lastly, we apply the reliable method to simulate and study the collisions of solitary waves numerically.     

Fully discrete linear approximation scheme for electric field diffusion in type-II superconductors

Superconductors are attracting physicists thanks to their ability to conduct electric current with virtually zero resistance. Their nonlinear behaviour opens, on the other hand, challenging problems for mathematicians. Our model of the diffusion of electric field in superconductors is based on three pillars: the eddy-current version of Maxwell’s equations, power law model of type-II superconductivity and linear dependence of magnetic induction on magnetic field. This leads to a time-dependent nonlinear degenerate partial differential equation. We propose a linear fully discrete approximation scheme to solve it. We have proven the convergence of the method and derived the error estimates describing the dependence of the error on the discretization parameters. These theoretical results were successfully confronted with numerical experiments.     

High-order linear compact conservative method for the nonlinear Schrödinger equation coupled with the nonlinear Klein–Gordon equation

In this paper, we design a linear-compact conservative numerical scheme which preserves the original conservative properties to solve the Klein–Gordon–Schrödinger equation. The proposed scheme is based on using the finite difference method. The scheme is three-level and linear-implicit. Priori estimate and the convergence of the finite difference approximate solutions are discussed by the discrete energy method. Numerical results demonstrate that the present scheme is conservative, efficient and of high accuracy.