A backward Euler alternating direction implicit difference scheme for the three‐dimensional fractional evolution equation

A backward Euler alternating direction implicit (ADI) difference scheme is formulated and analyzed for the three‐dimensional fractional evolution equation. In our method, the Riemann‐Liouville fractional integral term is treated by means of first order convolution quadrature suggested by Lubich. Meanwhile, an ADI technique is adopted to reduce the multidimensional problem to a series of one‐dimensional problems. A fully discrete difference scheme is constructed with space discretization by finite difference method. Two new inner products and corresponding norms are defined to analyze the scheme. The verification of stability and convergence is based on the nonnegative character of the real quadratic form associated with the convolution quadrature. Numerical experiments are reported to demonstrate the efficiency of our scheme.     

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.     

Methods for the numerical solution of the Benjamin‐Bona‐Mahony‐Burgers equation

L^∞‐error estimates for finite element for Galerkin solutions for the Benjamin‐Bona‐Mahony‐Burgers (BBMB) equation are considered. A priori bound and the semidiscrete Galerkin scheme are studied using appropriate projections. For fully discrete Galerkin schemes, we consider the backward Euler method and analyze the corresponding error estimates. For a second order accuracy in time, we propose a three‐level backward method. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Finite element approximation of a forward and backward anisotropic diffusion model in image denoising and form generalization

A new forward–backward anisotropic diffusion model is introduced. The two limit cases are the Perona‐Malik equation and the Total Variation flow model. A fully discrete finite element scheme is studied using C⁰‐piecewise linear elements in space and the backward Euler difference scheme in time. A priori estimates are proven. Numerical results in image denoising and form generalization are presented.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

一类带弱奇异核非线性偏积分微分方程的全离散有限元

1引言我们将研究下面一类带弱奇异核非线性偏积分微分方程的数值解:u_t-▽·(a(u)▽u)-integral from n=0 to tβ(t-s)△u(s)ds=f(u),x∈Ω,t∈(?),(1.1) u(·,t)=0,x∈(?)Ω,t∈J,(1.2) u(·,0)=v(x),x∈Ω,(1.3)其中Ω为平面上的凸角域,J=(0,T],α和f为R上的光滑函数,满足0相似文献   

Positive numerical splitting method for the Hull and White 2D Black–Scholes equation

We consider the locally one‐dimensional backward Euler splitting method to solve numerically the Hull and White problem for pricing European options with stochastic volatility in the presence of a mixed derivative term. We prove the first‐order convergence of the time‐splitting. The parabolic equation degenerates on the boundary x = 0 and we apply a fitted finite volume scheme to the equation to resolve the degeneracy and derive the fully discrete problem as we also investigate the discrete maximum principle. Numerical experiments illustrate the efficiency of our difference scheme. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 822–846, 2015     

A finite element method for the Sivashinsky equation

The Sivashinsky equation is a nonlinear evolutionary equation of fourth order in space. In this paper we have analyzed a semidiscrete finite element method and completely discrete scheme based on the backward Euler method and Crank–Nicolson–Galerkin scheme. A linearized backward Euler method have been developed and error bounds are derived for an L² projection.     

Fast evaluation and high accuracy finite element approximation for the time fractional subdiffusion equation

In this article, an efficient algorithm for the evaluation of the Caputo fractional derivative and the superconvergence property of fully discrete finite element approximation for the time fractional subdiffusion equation are considered. First, the space semidiscrete finite element approximation scheme for the constant coefficient problem is derived and supercloseness result is proved. The time discretization is based on the L1‐type formula, whereas the space discretization is done using, the fully discrete scheme is developed. Under some regularity assumptions, the superconvergence estimate is proposed and analyzed. Then, extension to the case of variable coefficients is also discussed. To reduce the computational cost, the fast evaluation scheme of the Caputo fractional derivative to solve the fractional diffusion equations is designed. Finally, numerical experiments are presented to support the theoretical results.     

Superconvergence analysis of finite element method for Poisson–Nernst–Planck equations

This article concerns with the superconvergence analysis of bilinear finite element method (FEM) for nonlinear Poisson–Nernst–Planck (PNP) equations. By employing high accuracy integral identities together with mean value technique, the superclose estimates in H¹‐norm are derived for the semi‐discrete and the backward Euler fully‐discrete schemes, which improve the suboptimal error estimate in L²‐norm in the previous literature. Furthermore, the global superconvergence results in H¹‐norm are obtained through interpolation postprocessing approach. Finally, a numerical example is provided to confirm the theoretical analysis.     

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.     

A fully discrete spectral method for fractional Cattaneo equation based on Caputo–Fabrizo derivative

Recently Caputo and Fabrizio introduced a new derivative with fractional order without singular kernel. The derivative can be used to describe the material heterogeneities and the fluctuations of different scales. In this article, we derived a new discretization of Caputo–Fabrizio derivative of order α (1 < α < 2) and applied it into the Cattaneo equation. A fully discrete scheme based on finite difference method in time and Legendre spectral approximation in space is proposed. The stability and convergence of the fully discrete scheme are rigorously established. The convergence rate of the fully discrete scheme in H¹ norm is O(τ² + N^1?m), where τ, N and m are the time‐step size, polynomial degree and regularity in the space variable of the exact solution, respectively. Furthermore, the accuracy and applicability of the scheme are confirmed by numerical examples to support the theoretical results.     

Finite difference/$H^1$-Galerkin MFE procedure for a fractional water wave model

In this article, an $H^1$-Galerkin mixed finite element (MFE) method for solving the time fractional water wave model is presented. First-order backward Euler difference method and $L1$ formula are applied to approximate integer derivative and Caputo fractional derivative with order $1/2$, respectively, and $H^1$-Galerkin mixed finite element method is used to approximate the spatial direction. The analysis of stability for fully discrete mixed finite element scheme is made and the optimal space-time orders of convergence for two unknown variables in both $H^1$-norm and $L^2$-norm are derived. Further, some computing results for a priori analysis and numerical figures based on four changed parameters in the studied problem are given to illustrate the effectiveness of the current method     

一类非线性偏积分微分方程二阶差分全离散格式

给出了数值求解一类非线性偏积分微分方程的二阶全离散差分格式.采用了二阶向后差分格式,积分项的离散利用了Lubich的二阶卷积求积公式,给出了稳定性的证明、误差估计及收敛性的结果.     

A high‐order B‐spline collocation scheme for solving a nonhomogeneous time‐fractional diffusion equation

A high‐accuracy numerical approach for a nonhomogeneous time‐fractional diffusion equation with Neumann and Dirichlet boundary conditions is described in this paper. The time‐fractional derivative is described in the sense of Riemann‐Liouville and discretized by the backward Euler scheme. A fourth‐order optimal cubic B‐spline collocation (OCBSC) method is used to discretize the space variable. The stability analysis with respect to time discretization is carried out, and it is shown that the method is unconditionally stable. Convergence analysis of the method is performed. Two numerical examples are considered to demonstrate the performance of the method and validate the theoretical results. It is shown that the proposed method is of order O(Δx⁴ + Δt^{2 ? α}) convergence, where α ∈ (0,1) . Moreover, the impact of fractional‐order derivative on the solution profile is investigated. Numerical results obtained by the present method are compared with those obtained by the method based on standard cubic B‐spline collocation method. The CPU time for present numerical method and the method based on cubic B‐spline collocation method are provided.     

Fully discrete A‐ϕ finite element method for Maxwell's equations with nonlinear conductivity

This article is devoted to the study of a fully discrete A ‐ finite element method to solve nonlinear Maxwell's equations based on backward Euler discretization in time and nodal finite elements in space. The nonlinearity is owing to a field‐dependent conductivity with the power‐law form . We design a nonlinear time‐discrete scheme for approximation in suitable function spaces. We show the well‐posedness of the problem, prove the convergence of the semidiscrete scheme based on the boundedness of the second derivative in the dual space and derive its error estimate. The Minty–Browder technique is introduced to obtain the convergence of the nonlinear term. Finally, we discuss the error estimate for the fully discretized problem and support the theoretical result by two numerical experiments. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 2083–2108, 2014     

非定常的热传导-对流问题的非线性Galerkin混合元法(Ⅱ):向后一步的Euler全离散化格式

In this paper, a fully discrete format of nonlinear Galerkin mixed element method with backward one-step Euler discretization of time for the non stationary conduction-convection problems is presented. The scheme is based on two finite element spaces XH and Xh for the approximation of the velocity, defined respectively on a coarse grid with grids size H and another fine grid with grid size h<< H, a finite element space Mh for the approximation of the pressure and two finite element spaces AH and Wh, for the approximation of the temperature,also defined respectivply on the coarse grid with grid size H and another fine grid with grid size h. The existence and the convergence of the fully discrete mixed element solution are shown. The scheme consists in using standard backward one step Euler-Galerkin fully discrete format at first L0 steps (L0 2) on fine grid with grid size h, but using nonlinear Galerkin mixed element method of backward one step Euler-Galerkin fully discrete format through L0 + 1 step to end step. We have proved that the fully discrete nonlinear Galerkin mixed element procedure with respect to the coarse grid spaces with grid size H holds superconvergence.     

One‐level and two‐level finite element methods for the Boussinesq problem

Previous works on the convergence of numerical methods for the Boussinesq problem were conducted, while the optimal L²‐norm error estimates for the velocity and temperature are still lacked. In this paper, the backward Euler scheme is used to discrete the time terms, standard Galerkin finite element method is adopted to approximate the variables. The MINI element is used to approximate the velocity and pressure, the temperature field is simulated by the linear polynomial. Under some restriction on the time step, we firstly present the optimal L² error estimates of approximate solutions. Secondly, two‐level method based on Stokes iteration for the Boussinesq problem is developed and the corresponding convergence results are presented. By this method, the original problem is decoupled into two small linear subproblems. Compared with the standard Galerkin method, the two‐level method not only keeps good accuracy but also saves a lot of computational cost. Finally, some numerical examples are provided to support the established theoretical analysis.     

非定常的热传导-对流问题的非线性Galerkin混合元法(Ⅱ):向后一步的Euler全离散化格式

1．引言本文的工作主要是讨论非定常的热传导一对流问题的向后一步的Euler全离散化的非线性Galerkin混合元解的存在性及其误差估计．该工作是对山中的同一问题研究的第二部分．在第一部分［1］，我们已经讨论了此问题的半离散化的情形．由于所研究的目标都是非定常的热传导一对流问题，其背景是相同的，在此将不重复了，请参考［1］．本文的安排如下，52先回顾非定常的热传导一对流问题的混合元解的经典性质．53回顾半离散化的非线性Galerkin混合元解的性质，并导出后续讨论需要的一些关于时间导数的估计．54讨论向后一步的Euler全离散化…     

非定常Stokes方程的全离散稳定化有限元格式

本文研究二维非定常Stokes方程全离散稳定化有限元方法.首先给出关于时间向后一步Euler半离散格式,然后直接从该时间半离散格式出发,构造基于两局部高斯积分的稳定化全离散有限元格式,其中空间用P_1—P_1元逼近,证明有限元解的误差估计.本文的研究方法使得理论证明变得更加简便,也是处理非定常Stokes方程的一种新的途径.     

计算高维带弱奇异核发展型方程的交替方向隐式欧拉方法

本文主要研究高维带弱奇异核的发展型方程的交替方向隐式(ADI)差分方法.向后欧拉(Euler)方法联立一阶卷积求积公式处理时间方向的离散,有限差分方法处理空间方向的离散,并进一步构造了ADI全离散差分格式.然后将二维问题延伸到三维问题,构造三维空间问题的ADI差分格式.基于离散能量法,详细证明了全离散格式的稳定性和误差分析.随后给出了2个数值算例,数值结果进一步验证了时间方向的收敛阶为一阶,空间方向的收敛阶为二阶,和理论分析结果一致.