A finite difference/finite element technique with error estimate for space fractional tempered diffusion-wave equation

An efficient numerical technique is proposed to solve one- and two-dimensional space fractional tempered fractional diffusion-wave equations. The space fractional is based on the Riemann–Liouville fractional derivative. At first, the temporal direction is discretized using a second-order accurate difference scheme. Then a classic Galerkin finite element is employed to obtain a full-discrete scheme. Furthermore, for the time-discrete and the full-discrete schemes error estimate has been presented to show the unconditional stability and convergence of the developed numerical method. Finally, two test problems have been illustrated to verify the efficiency and simplicity of the proposed technique.     

一维二阶椭圆和抛物型微分方程的高精度有限体积元方法

§1.引 言 有限体积元方法作为偏微分方程的求解新技术,日益受到重视,该方法从微分方程的积分守恒形式出发,通过选取试探函数空间为一次有限元空间来导出计算格式[1-4].由于该方法具有非常好的质量守恒性质,在计算流体力学领域得到了广泛的应用[5].就方法而言,有限体积元方法相当于李荣华教授提出的广义差分方法[6-10]的特殊情形,即取试探函数     

Adaptive dimension splitting algorithm for three-dimensional elliptic equations

An adaptive dimension splitting algorithm for three-dimensional (3D) elliptic equations is presented in this paper. We propose residual and recovery-based error estimators with respect to X?Y plane direction and Z direction, respectively, and construct the corresponding adaptive algorithm. Two-sided bounds of the estimators guarantee the efficiency and reliability of such error estimators. Numerical examples verify their efficiency both in estimating the error and in refining the mesh adaptively. This algorithm can be compared with or even better than the 3D adaptive finite element method with tetrahedral elements in some cases. What is more, our new algorithm involves only two-dimensional mesh and one-dimensional mesh in the process of refining mesh adaptively, and it can be implemented in parallel.     

A fully discrete local discontinuous Galerkin method for one-dimensional time-fractional Fisher's equation

In this paper, we consider the local discontinuous Galerkin (LDG) finite element method for one-dimensional time-fractional Fisher's equation, which is obtained from the standard one-dimensional Fisher's equation by replacing the first-order time derivative with a fractional derivative (of order α, with 0<α<1). The proposed LDG is based on the LDG finite element method for space and finite difference method for time. We prove that the method is stable, and the numerical solution converges to the exact one with order O(h^k+1+τ^2?α), where h, τ and k are the space step size, time step size, polynomial degree, respectively. The numerical experiments reveal that the LDG is very effective.     

分数阶对流-弥散方程的移动网格有限元方法

相比经典的对流-弥散方程,分数微分算子的非局部性质导致分数阶对流-弥散方程(FADE)的有限元方法在每个单元上的计算都联系一个带弱奇异核的数值积分.当弥散项分数阶μ接近1时,穿透曲线出现重度拖尾,数值解产生振荡.研究表明:时间半离散后的FADE在特殊的变分形式下,有限元刚度矩阵有直接计算公式;以De Boor算法为基础的移动网格方法能很好地消除数值振荡.     

Numerical analysis of history-dependent hemivariational inequalities and applications to viscoelastic contact problems with normal penetration

In this paper numerical approximation of history-dependent hemivariational inequalities with constraint is considered, and corresponding Céa’s type inequality is derived for error estimate. For a viscoelastic contact problem with normal penetration, an optimal order error estimate is obtained for the linear element method. A numerical experiment for the contact problem is reported which provides numerical evidence of the convergence order predicted by the theoretical analysis.     

A wavelet approach for the multi-term time fractional diffusion-wave equation

In this paper, a Galerkin method based on the second kind Chebyshev wavelets (SKCWs) is established for solving the multi-term time fractional diffusion-wave equation. To do this, a new operational matrix of fractional integration for the SKCWs must be derived and in order to improve the computational efficiency, the hat functions are proposed to create a general procedure for constructing this matrix. Implementation of these wavelet basis functions and their operational matrix of fractional integration simplifies the problem under consideration to a system of linear algebraic equations, which greatly decreases the computational cost for finding an approximate solution. The main privilege of the proposed method is adjusting the initial and boundary conditions in the final system automatically. Theoretical error and convergence analysis of the SKCWs expansion approve the reliability of the approach. Also, numerical investigation reveals the applicability and accuracy of the presented method.     

Quadratic finite volume element method for the air pollution model

In this article, we use the quadratic finite volume element method (FVEM) to solve the problem of the air pollution model, choose the trial function spaces as the Lagrange quadratic element function spaces and the test function spaces as the piecewise constant function spaces, then get the error estimates of L ² and H ¹. Finally, by numerical experiments, we analyse and compare the FVEM with the finite difference method, and the numerical results we get show that the FVEM is much better and more effective, so this article has some practical significance in the improvement and control of air pollution.     

一类二维粘性波动方程的交替方向有限体积元方法

针对二维粘性波动方程模型问题,提出了一类基于双线性插值的交替方向有限体积元方法,并给出了两种具体计算格式,一是基于有限差分方法中Douglas思想的格式,二是一类推广型的局部一维格式.分析证明了该方法按照L~2范数在时间和空间方向均有二阶收敛精度.最后,数值算例验证了算法的有效性和精确性.     

Error Estimates for Semidiscrete Finite Element Approximations of the Stokes Equations Under Minimal Regularity Assumptions

Semidiscrete (spatially discrete) finite element approximations of the Stokes equations are studied in this paper. Properties of L ², H ¹ and H ^–1 projections onto discretely divergence-free spaces are discussed and error estimates are derived under minimal regularity assumptions on the solution.     

二维土壤水中溶质运移问题的全离散混合元法研究

本文建立了二维非粘性土壤水中溶质运移问题的混合元格式,讨论了混合元解的存在唯一性,并分析了误差估计.最后给出数值算例,数值模拟结果表明,用该方法模拟溶质运移问题是合理有效的.     

A Galerkin finite element scheme for time–space fractional diffusion equation

In this paper, a Galerkin finite element scheme to approximate the time–space fractional diffusion equation is studied. Firstly, the fractional diffusion equation is transformed into a fractional Volterra integro-differential equation. And a second-order fractional trapezoidal formula is used to approach the time fractional integral. Then a Galerkin finite element method is introduced in space direction, where the semi-discretization scheme and fully discrete scheme are given separately. The stability analysis of semi-discretization scheme is discussed in detail. Furthermore, convergence analysis of semi-discretization scheme and fully discrete scheme are given in details. Finally, two numerical examples are displayed to demonstrate the effectiveness of the proposed method.     

Adaptive Finite Element Methods for Microstructures? Numerical Experiments for a 2-Well Benchmark

Macroscopic simulations of non-convex minimisation problems with enforced microstructures encounter oscillations on finest length scales – too fine to be fully resolved. The numerical analysis must rely on an essentially equivalent relaxed mathematical model. The paper addresses a prototype example, the scalar 2-well minimisation problem and its convexification and introduces a benchmark problem with a known (generalised) solution. For this benchmark, the stress error is studied empirically to asses the performance of adaptive finite element methods for the relaxed and the original minimisation problem. Despite the theoretical reliability-efficiency gap for the relaxed problem, numerical evidence supports that adaptive mesh-refining algorithms generate efficient triangulations and improve the experimental convergence rates optimally. Moreover, the averaging error estimators perform surprisingly accurate.     

An efficient approximate method for solving linear fractional Klein–Gordon equation based on the generalized Laguerre polynomials

In this paper, a new approximate formula of the fractional derivative is derived. The proposed formula is based on the generalized Laguerre polynomials. Global approximations to functions defined on a semi-infinite interval are constructed. The fractional derivatives are presented in terms of Caputo sense. Special attention is given to study the error and the convergence analysis of the proposed formula. A new spectral Laguerre collocation method is presented for solving linear fractional Klein–Gordon equation (LFKGE). The properties of Laguerre polynomials are utilized to reduce LFKGE to a system of ordinary differential equations, which solved using the finite difference method. Numerical results are provided to confirm the theoretical results and the efficiency of the proposed method.     

Dispersive and Dissipative Properties of Discontinuous Galerkin Finite Element Methods for the Second-Order Wave Equation

Discontinuous Galerkin finite element methods (DGFEM) offer certain advantages over standard continuous finite element methods when applied to the spatial discretisation of the acoustic wave equation. For instance, the mass matrix has a block diagonal structure which, used in conjunction with an explicit time stepping scheme, gives an extremely economical scheme for time domain simulation. This feature is ubiquitous and extends to other time-dependent wave problems such as Maxwell’s equations. An important consideration in computational wave propagation is the dispersive and dissipative properties of the discretisation scheme in comparison with those of the original system. We investigate these properties for two popular DGFEM schemes: the interior penalty discontinuous Galerkin finite element method applied to the second-order wave equation and a more general family of schemes applied to the corresponding first order system. We show how the analysis of the multi-dimensional case may be reduced to consideration of one-dimensional problems. We derive the dispersion error for various schemes and conjecture on the generalisation to higher order approximation in space     

A characteristic centred finite difference method for a 2D air pollution model

In this paper, we propose a characteristic centred finite difference method on non-uniform grids to solve the problem of the air pollution model. Numerical solutions and error estimates of the air pollution concentration and its first-order derivatives for space variables are obtained. The computational cost of the method is the same as that of the characteristic difference method based on a linear interpolation. The error order of the numerical solutions is the same as that of the characteristic difference method based on a quadratic interpolation. At last, we give numerical examples to illustrate feasibility and efficiency of this method.