On the coupling of finite volume and discontinuous Galerkin method for elliptic problems

The coupling of cell-centered finite volume method with primal discontinuous Galerkin method is introduced in this paper for elliptic problems. Convergence of the method with respect to the mesh size is proved. Numerical examples confirm the theoretical rates of convergence. Advantages of the coupled scheme are shown for problems with discontinuous coefficients or anisotropic diffusion matrix.     

Solution of double nonlinear problems in porous media by a combined finite volume–finite element algorithm

The combined finite volume–finite element scheme for a double nonlinear parabolic convection-dominated diffusion equation which models the variably saturated flow and contaminant transport problems in porous media is extended. Whereas the convection is approximated by a finite volume method (Multi-Point Flux Approximation), the diffusion is approximated by a finite element method. The scheme is fully implicit and involves a relaxation-regularized algorithm. Due to monotonicity and conservation properties of the approximated scheme and in view of the compactness theorem we show the convergence of the numerical scheme to the weak solution. Our scheme is applied for computing two dimensional examples with different degrees of complexity. The numerical results demonstrate that the proposed scheme gives good performance in convergence and accuracy.     

On the strong convergence rate for the Euler–Maruyama scheme of one-dimensional SDEs with irregular diffusion coefficient and local time

The strong rate of convergence for the Euler–Maruyama scheme of stochastic differential equations (SDEs) driven by a Brownian motion with Hölder continuous diffusion coefficient or irregular drift coefficient have been widely studied. In the case of irregular diffusion coefficient, however, there are few studies. In this article, under Le Gall's condition on the diffusion coefficient, which leads to conclude the pathwise uniqueness for SDEs, we provide the same result on the strong rate of convergence as in the case of 1/2-Hölder continuous diffusion coefficient. The idea of the proof is to use a version of Avikainen's inequality. As an application, we introduce a numerical scheme for SDEs with local time.     

非光滑时间分数阶齐次扩散方程的修正

当初值不光滑时,时间分数阶齐次扩散方程数值方法的精度会下降.为了得到高阶时间收敛格式,提出加权移位的Grünwald-Letnikov的修正格式,运用Lubich的修正方法,得到非光滑时间分数阶齐次扩散方程的收敛阶仍为O(k²).最后,通过数值算例验证了数值计算结果与理论计算结果一致.     

A finite volume scheme for convection–diffusion equations with nonlinear diffusion derived from the Scharfetter–Gummel scheme

We propose a finite volume scheme for convection–diffusion equations with nonlinear diffusion. Such equations arise in numerous physical contexts. We will particularly focus on the drift-diffusion system for semiconductors and the porous media equation. In these two cases, it is shown that the transient solution converges to a steady-state solution as t tends to infinity. The introduced scheme is an extension of the Scharfetter–Gummel scheme for nonlinear diffusion. It remains valid in the degenerate case and preserves steady-states. We prove the convergence of the scheme in the nondegenerate case. Finally, we present some numerical simulations applied to the two physical models introduced and we underline the efficiency of the scheme to preserve long-time behavior of the solutions.     

三维多面体网格上扩散方程的保正格式

 针对三维任意(星形)多面体网格, 本文构造了扩散方程的一种单元中心型非线性有限体积格式, 证明了该格式具有保正性. 在该格式设计中, 除引入网格中心量外, 还引入网格节点量和网格面中心量作为中间未知量, 它们将用网格中心未知量线性组合表示, 使得格式仅有网格中心未知量作为基本未知量. 在节点量计算中, 利用网格面上的调和平均点, 设计了一种适用于三维多面体网格的局部显式加权方法. 该格式适用于求解非平面的网格表面和间断扩散系数的问题. 数值例子验证了它对光滑解具有二阶精度和保正性.     

Approximations of non-smooth integral type functionals of one dimensional diffusion processes

In this article, we obtain the weak and strong rates of convergence of time integrals of non-smooth functions of a one dimensional diffusion process. We propose the use of the exact simulation scheme to simulate the process at discretization points. In particular, we also present the rates of convergence for the weak and strong errors of approximation for the local time of a one dimensional diffusion process as an application of our method.     

Fast second-order accurate difference schemes for time distributed-order and Riesz space fractional diffusion equations

The aim of this paper is to develop fast second-order accurate difference schemes for solving one- and two-dimensional time distributed-order and Riesz space fractional diffusion equations. We adopt the same measures for one- and two-dimensional problems as follows: we first transform the time distributed-order fractional diffusion problem into the multi-term time-space fractional diffusion problem with the composite trapezoid formula. Then, we propose a second-order accurate difference scheme based on the interpolation approximation on a special point to solve the resultant problem. Meanwhile, the unconditional stability and convergence of the new difference scheme in $L_2$-norm are proved. Furthermore, we find that the discretizations lead to a series of Toeplitz systems which can be efficiently solved by Krylov subspace methods with suitable circulant preconditioners. Finally, numerical results are presented to show the effectiveness of the proposed difference methods and demonstrate the fast convergence of our preconditioned Krylov subspace methods.     

The Maximum Principle for Time-Fractional Diffusion Equations and Its Application

Using an equivalent reformulation of the definition of the Caputo fractional derivative, we present the stability and convergence analysis for a finite difference scheme that we use to solve a time-fractional diffusion equation. The analysis is based on a maximum principle for such equations.     

A new reliable algorithm based on the sinc function for the time fractional diffusion equation

This paper deals with the numerical solution of time fractional diffusion equation. In this work, we consider the fractional derivative in the sense of Riemann-Liouville. At first, the time fractional derivative is discretized by integrating both sides of the equation with respect to the time variable and we arrive at a semi–discrete scheme. The stability and convergence of time discretized scheme are proven by using the energy method. Also we show that the convergence order of this scheme is O(τ2?α). Then we use the sinc collocation method to approximate the solution of semi–discrete scheme and show that the problem is reduced to a Sylvester matrix equation. Besides by performing some theorems, the exponential convergence rate of sinc method is illustrated. The numerical experiments are presented to show the excellent behavior and high accuracy of the proposed hybrid method in comparison with some other well known methods.     

A mixed-step algorithm for the approximation of the stationary regime of a diffusion

In some recent papers, some procedures based on some weighted empirical measures related to decreasing-step Euler schemes have been investigated to approximate the stationary regime of a diffusion (possibly with jumps) for a class of functionals of the process. This method is efficient but needs the computation of the function at each step. To reduce the complexity of the procedure (especially for functionals), we propose in this paper to study a new scheme, called the mixed-step scheme, where we only keep some regularly time-spaced values of the Euler scheme. Our main result is that, when the coefficients of the diffusion are smooth enough, this alternative does not change the order of the rate of convergence of the procedure. We also investigate a Richardson–Romberg method to speed up the convergence and show that the variance of the original algorithm can be preserved under a uniqueness assumption for the invariant distribution of the “duplicated” diffusion, condition which is extensively discussed in the paper. Finally, we conclude by giving sufficient “asymptotic confluence” conditions for the existence of a smooth solution to a discrete version of the associated Poisson equation, condition which is required to ensure the rate of convergence results.     

Convergence of an implicit Voronoi finite volume method for reaction–diffusion problems

We investigate the convergence of an implicit Voronoi finite volume method for reaction–diffusion problems including nonlinear diffusion in two space dimensions. The model allows to handle heterogeneous materials and uses the chemical activities of the involved species as primary variables. The numerical scheme works with boundary conforming Delaunay meshes and preserves positivity and the dissipative property of the continuous system. Starting from a result on the global stability of the scheme (uniform, mesh‐independent global upper, and lower bounds), we prove strong convergence of the chemical activities and their gradients to a weak solution of the continuous problem. To illustrate the preservation of qualitative properties by the numerical scheme, we present a long‐term simulation of the Michaelis–Menten–Henri system. Especially, we investigate the decay properties of the relative free energy over several magnitudes of time, and obtain experimental orders of convergence for this quantity. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 141–174, 2016     

Convergence of the Grünwald–Letnikov scheme for time-fractional diffusion

Using bivariate generating functions, we prove convergence of the Grünwald–Letnikov difference scheme for the fractional diffusion equation (in one space dimension) with and without central linear drift in the Fourier–Laplace domain as the space and time steps tend to zero in a well-scaled way. This implies convergence in distribution (weak convergence) of the discrete solution towards the probability of sojourn of a diffusing particle. The difference schemes allow also interpretation as discrete random walks. For fractional diffusion with central linear drift we show that in the Fourier–Laplace domain the limiting ordinary differential equation coincides with that for the solution of the corresponding diffusion equation.     

求解奇异摄动问题的一个耦合差分格式

本文考察奇异摄动问题(1.1).在一特殊的非均匀网格上,将不稳定、二阶精度的中心差格式和稳定、一阶精度的Abrahamsson-Keller-Kreiss箱子格式相耦合,得到了一个二阶一致收敛的差分格式.最后给出了数值结果.     

A new approximation for the generalized fractional derivative and its application to generalized fractional diffusion equation

In this paper, we derive a fourth order approximation for the generalized fractional derivative that is characterized by a scale function z(t) and a weight function w(t) . Combining the new approximation with compact finite difference method, we develop a numerical scheme for a generalized fractional diffusion problem. The stability and convergence of the numerical scheme are proved by the energy method, and it is shown that the temporal and spatial convergence orders are both 4. Several numerical experiments are provided to illustrate the efficiency of our scheme.     

Numerical simulation for the variable-order Galilei invariant advection diffusion equation with a nonlinear source term

In this paper, we consider the variable-order Galilei advection diffusion equation with a nonlinear source term. A numerical scheme with first order temporal accuracy and second order spatial accuracy is developed to simulate the equation. The stability and convergence of the numerical scheme are analyzed. Besides, another numerical scheme for improving temporal accuracy is also developed. Finally, some numerical examples are given and the results demonstrate the effectiveness of theoretical analysis.     

一阶双曲问题的间断流线扩散法

1．引言众所周知，求解一阶双曲问题的Galerkin有限元法，仅具有次最优LZ一收敛阶估计，难于建立H＼误差估计【‘1，且Gderkill有限元解常呈现伪数值振荡．为改善计算精度与稳定性，诸多非标准有限元解法相继提出，其中，间断（Discontinuous）Galerkin有限元法（以下简称DG方法）与流线扩散（StreamlineDiffusion）有限元法（以下简称SD法）是两种具有鲜明特点，较为成功的算法．具体地，DG方法是一种迎风型显式算法，它从入流边界开始，沿流场方向，自上游往下游，逐个单元进行解算，计算十分简便且可局部并行化．SD方法则是一种P…     

A uniformly convergent scheme for a system of reaction–diffusion equations

In this work a system of two parabolic singularly perturbed equations of reaction–diffusion type is considered. The asymptotic behaviour of the solution and its partial derivatives is given. A decomposition of the solution in its regular and singular parts has been used for the asymptotic analysis of the spatial derivatives. To approximate the solution we consider the implicit Euler method for time stepping and the central difference scheme for spatial discretization on a special piecewise uniform Shishkin mesh. We prove that this scheme is uniformly convergent, with respect to the diffusion parameters, having first-order convergence in time and almost second-order convergence in space, in the discrete maximum norm. Numerical experiments illustrate the order of convergence proved theoretically.     

Numerical solutions of diffusion mathematical models with delay

In this paper an explicit numerical difference scheme for mixed problems for the delay diffusion equation is proposed, as a generalization of the classic difference scheme for the diffusion problem. A sufficient condition for the asymptotic stability of the new scheme is proved. Consistence, convergence and some properties of stability for this scheme are studied. Illustrative examples of numerical results are also included.     

Strong convergence rates for backward Euler–Maruyama method for non-linear dissipative-type stochastic differential equations with super-linear diffusion coefficients

In this work, we generalize the current theory of strong convergence rates for the backward Euler–Maruyama scheme for highly non-linear stochastic differential equations, which appear in both mathematical finance and bio-mathematics. More precisely, we show that under a dissipative condition on the drift coefficient and super-linear growth condition on the diffusion coefficient the BEM scheme converges with strong order of a half. This type of convergence gives theoretical foundations for efficient variance reduction techniques for Monte Carlo simulations. We support our theoretical results with relevant examples, such as stochastic population models and stochastic volatility models.