A mortar mimetic finite difference method on non-matching grids

We consider mimetic finite difference approximations to second order elliptic problems on non-matching multiblock grids. Mortar finite elements are employed on the non-matching interfaces to impose weak flux continuity. Optimal convergence and, in certain cases, superconvergence is established for both the scalar variable and its flux. The theory is confirmed by computational results. Supported by the US Department of Energy, under contractW-7405-ENG-36. LA-UR-04-4740. Partially supported by NSF grants EIA 0121523 and DMS 0411413, by NPACI grant UCSD 10181410, and by DOE grant DE-FGO2-04ER25617. Partially supported by NSF grants DMS 0107389, DMS 0112239 and DMS 0411694 and by DOE grant DE-FG02-04ER25618.     

Uniform approximation of singularly perturbed reaction‐diffusion problems by the finite element method on a Shishkin mesh

We consider the numerical approximation of singularly perturbed reaction‐diffusion problems over two‐dimensional domains with smooth boundary. Using the h version of the finite element method over appropriately designed piecewise uniform (Shishkin) meshes, we are able to uniformly approximate the solution at a quasi‐optimal rate. The results of numerical computations showing agreement with the analysis are also presented. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 89–111, 2003     

Pointwise error estimates of the bilinear SDFEM on Shishkin meshes

A model singularly perturbed convection–diffusion problem in two space dimensions is considered. The problem is solved by a streamline diffusion finite element method (SDFEM) that uses piecewise bilinear finite elements on a Shishkin mesh. We prove that the method is convergent, independently of the diffusion parameter ε, with a pointwise accuracy of almost order 11/8 outside and inside the boundary layers. Numerical experiments support these theoretical results. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

The characteristic finite difference streamline diffusion method for convection-dominated diffusion problems

In this paper, we consider the characteristic finite difference streamline diffusion method for two-dimensional convection-dominated diffusion problems. The scheme is combined the method of characteristics with the finite difference streamline diffusion (FDSD) method to create the characteristic FDSD (C-FDSD) procedures. Stability analysis and error estimate of the C-FDSD method are deduced. The scheme not only realizes the purpose of lowering the time-truncation error, using larger time step for solving the convection-dominated diffusion problems, but also keeps the favorable stability and high precision of the FDSD method. Finally, numerical experiments are presented to illustrate the availability of the scheme.     

High order methods on Shishkin meshes for singular perturbation problems of convection–diffusion type

In this paper we construct and analyze two compact monotone finite difference methods to solve singularly perturbed problems of convection–diffusion type. They are defined as HODIE methods of order two and three, i.e., the coefficients are determined by imposing that the local error be null on a polynomial space. For arbitrary meshes, these methods are not adequate for singularly perturbed problems, but using a Shishkin mesh we can prove that the methods are uniformly convergent of order two and three except for a logarithmic factor. Numerical examples support the theoretical results. This revised version was published online in June 2006 with corrections to the Cover Date.     

Parallel Galerkin domain decomposition procedures based on the streamline diffusion method for convection-diffusion problems

Based upon the streamline diffusion method, parallel Galerkin domain decomposition procedures for convection-diffusion problems are given. These procedures use implicit method in the sub-domains and simple explicit flux calculations on the inter-boundaries of sub-domains by integral mean method or extrapolation method to predict the inner-boundary conditions. Thus, the parallelism can be achieved by these procedures. The explicit nature of the flux calculations induces a time step limitation that is necessary to preserve stability. Artificial diffusion parameters δ are given. By analysis, optimal order error estimate is derived in a norm which is stronger than L²-norm for these procedures. This error estimate not only includes the optimal H¹-norm error estimate, but also includes the error estimate along the streamline direction ‖β⋅∇(u−U)‖, which cannot be achieved by standard finite element method. Experimental results are presented to confirm theoretical results.     

A new fitted operator finite difference method to solve systems of evolutionary reaction-diffusion equations

Abstract

In recent years, fitted operator finite difference methods (FOFDMs) have been developed for numerous types of singularly perturbed ordinary differential equations. The construction of most of these methods differed though the final outcome remained similar. The most crucial aspect was how the difference operator was designed to approximate the differential operator in question. Very often the approaches for constructing these operators had limited scope in the sense that it was difficult to extend them to solve even simple one-dimensional singularly perturbed partial differential equations. However, in some of our most recent work, we have successfully designed a class of FOFDMs and extended them to solve singularly perturbed time-dependent partial differential equations. In this paper, we design and analyze a robust FOFDM to solve a system of coupled singularly perturbed parabolic reaction-diffusion equations. We use the backward Euler method for the semi-discretization in time. An FOFDM is then developed to solve the resulting set of boundary value problems. The proposed method is analyzed for convergence. Our method is uniformly convergent with order one and two, respectively, in time and space, with respect to the perturbation parameters. Some numerical experiments supporting the theoretical investigations are also presented.     

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.     

A supercloseness result for the discontinuous Galerkin stabilization of convection–diffusion problems on Shishkin meshes

We consider a convection–diffusion problem with Dirichlet boundary conditions posed on a unit square. The problem is discretized using a combination of the standard Galerkin FEM and an h–version of the nonsymmetric discontinuous Galerkin FEM with interior penalties on a layer–adapted mesh with linear/bilinear elements. With specially chosen penalty parameters for edges from the coarse part of the mesh, we prove uniform convergence (in the perturbation parameter) in an associated norm. In the same norm we also establish a supercloseness result. Numerical tests support our theoretical estimates.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Uniform superconvergence of a Galerkin finite element method on Shishkin‐type meshes

We consider a Galerkin finite element method that uses piecewise bilinears on a class of Shishkin‐type meshes for a model singularly perturbed convection‐diffusion problem on the unit square. The method is shown to be convergent, uniformly in the diffusion parameter ϵ, of almost second order in a discrete weighted energy norm. As a corollary, we derive global L₂‐norm error estimates and local L_∞‐norm estimates. Numerical experiments support our theoretical results. © 2000 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16:426–440, 2000     

Superconvergence of conforming finite element for fourth‐order singularly perturbed problems of reaction diffusion type in 1D

We consider conforming finite element approximation of fourth‐order singularly perturbed problems of reaction diffusion type. We prove superconvergence of standard C¹ finite element method of degree p on a modified Shishkin mesh. In particular, a superconvergence error bound of in a discrete energy norm is established. The error bound is uniformly valid with respect to the singular perturbation parameter ?. Numerical tests indicate that the error estimate is sharp. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 550–566, 2014     

Error estimation of a quadratic finite volume method on right quadrangular prism grids

In this paper, we develop a finite volume element method with affine quadratic bases on right quadrangular prism meshes for three-dimensional elliptic boundary value problems. The optimal H¹$H^1$-norm error estimate of second order accuracy is proved under certain assumptions about the meshes. Numerical results are presented to illustrate the theoretical analysis.     

A defect–correction parameter-uniform numerical method for a singularly perturbed convection–diffusion problem in one dimension

Monotone second order parameter-robust numerical methods for singularly perturbed differential equations can be designed using the principles of defect-correction. However, the proofs of second order parameter-uniform convergence can be difficult, even in one dimension. In this paper, we examine a variant of the standard defect-correction scheme. This variant is proposed in order to simplify the analysis of a defect-correction method in the case of a one dimensional convection–diffusion problem. Numerical results are presented to validate the theoretical results.*This research was partially supported by the project MEC/FEDER MTM 2004-019051 and the grant EUROPA XXI of the Caja de Ahorros de la Inmaculada.     

Parallel characteristic finite difference method for convection–diffusion equations

Based on the overlapping domain decomposition, an efficient parallel characteristic finite difference scheme is proposed for solving convection‐diffusion equations numerically. We give the optimal convergence order in error estimate analysis, which shows that we just need to iterate once or twice at each time level to reach the optimal convergence order. Numerical experiments also confirm the theoretical analysis. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 854–866, 2011     

A robust and accurate finite difference method for a generalized Black-Scholes equation

In this paper we present a numerical method for a generalized Black-Scholes equation, which is used for option pricing. The method is based on a central difference spatial discretization on a piecewise uniform mesh and an implicit time stepping technique. Our scheme is stable for arbitrary volatility and arbitrary interest rate, and is second-order convergent with respect to the spatial variable. Furthermore, the present paper efficiently treats the singularities of the non-smooth payoff function. Numerical results support the theoretical results.     

A second-order hybrid finite difference scheme for a system of singularly perturbed initial value problems

A system of coupled singularly perturbed initial value problems with two small parameters is considered. The leading term of each equation is multiplied by a small positive parameter, but these parameters may have different magnitudes. The solution of the system has boundary layers that overlap and interact. The structure of these layers is analyzed, and this leads to the construction of a piecewise-uniform mesh that is a variant of the usual Shishkin mesh. On this mesh a hybrid finite difference scheme is proved to be almost second-order accurate, uniformly in both small parameters. Numerical results supporting the theory are presented.     

Uniformly convergent novel finite difference methods for singularly perturbed reaction–diffusion equations

In this paper, we construct a kind of novel finite difference (NFD) method for solving singularly perturbed reaction–diffusion problems. Different from directly truncating the high‐order derivative terms of the Taylor's series in the traditional finite difference method, we rearrange the Taylor's expansion in a more elaborate way based on the original equation to develop the NFD scheme for 1D problems. It is proved that this approach not only can highly improve the calculation accuracy but also is uniformly convergent. Then, applying alternating direction implicit technique, the newly deduced schemes are extended to 2D equations, and the uniform error estimation based on Shishkin mesh is derived, too. Finally, numerical experiments are presented to verify the high computational accuracy and theoretical prediction.     

A finite element method for the numerical solution of convection-dominated anisotropic diffusion equations

Summary. The proposed method is based on an additive decomposition of the differential operator and the subsequent fitted discretization of the resulting components. For standard situations, the derived stability and error estimates in the energy norm qualitatively coincide with well-known estimates. In the case of small diffusion, a uniform error estimate with reduced order is obtained. Received August 7, 1997 / Revised version received July 15, 1998 / Published online December 6, 1999     

The role of coefficients of a general SPDE on the stability and convergence of a finite difference method

In this paper for the approximate solution of stochastic partial differential equations (SPDEs) of Itô-type, the stability and application of a class of finite difference method with regard to the coefficients in the equations is analyzed. The finite difference methods discussed here will be either explicit or implicit and a comparison between them will be reported. We prove the consistency and stability of these methods and investigate the influence of the multiplier (particularly multiplier of the random noise) in mean square stability. From stochastic version of Lax-Richtmyer the convergence of these methods under some conditions are established. Numerical experiments are included to show the efficiency of the methods.