A monotone finite element scheme for convection-diffusion equations

A simple technique is given in this paper for the construction and analysis of a class of finite element discretizations for convection-diffusion problems in any spatial dimension by properly averaging the PDE coefficients on element edges. The resulting finite element stiffness matrix is an M-matrix under some mild assumption for the underlying (generally unstructured) finite element grids. As a consequence the proposed edge-averaged finite element scheme is particularly interesting for the discretization of convection dominated problems. This scheme admits a simple variational formulation, it is easy to analyze, and it is also suitable for problems with a relatively smooth flux variable. Some simple numerical examples are given to demonstrate its effectiveness for convection dominated problems.

[1]  Lutz Tobiska,et al.  Numerical Methods for Singularly Perturbed Differential Equations , 1996 .

[2]  S. Selberherr.  Analysis and simulation of semiconductor devices , 1984 .

[3]  Jinchao Xu.  Two-grid Discretization Techniques for Linear and Nonlinear PDEs , 1996 .

[4]  Franco Brezzi,et al.  Numerical simulation of semiconductor devices , 1989 .

[5]  A. H. Schatz,et al.  An observation concerning Ritz-Galerkin methods with indefinite bilinear forms , 1974 .

[6]  Peter A. Markowich,et al.  Inverse-average-type nite element discretizations of self-adjoint second order elliptic problems , 1988 .

[7]  M.-D. D. Huang,et al.  The Constant-Flow Patch Test -- A Unique Guideline for the Evaluation of Discretization Schemes for the Current Continuity Equations , 1985, IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems.

[8]  Claes Johnson.  Numerical solution of partial differential equations by the finite element method , 1988 .

[9]  Paola Pietra,et al.  New mixed finite element schemes for current continuity equations , 1990 .

[10]  Song Wang,et al.  A triangular mixed finite element method for the stationary semiconductor device equations , 1991 .

[11]  P. Bassanini,et al.  Elliptic Partial Differential Equations of Second Order , 1997 .

[12]  Thomas Probst,et al.  Downwind Gauß-Seidel Smoothing for Convection Dominated Problems , 1997, Numer. Linear Algebra Appl..

[13]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[14]  T. Barth.  Aspects of Unstructured Grids and Finite-Volume Solvers for the Euler and Navier-Stokes Equations , 1994 .

[15]  Wolfgang Hackbusch,et al.  Downwind Gauß‐Seidel Smoothing for Convection Dominated Problems , 1997 .

[16]  M. S. Mock,et al.  ANALYSIS OF A DISCRETIZATION ALGORITHM FOR STATIONARY CONTINUITY EQUATIONS IN SEMICONDUCTOR DEVICE MODELS, II , 1983 .

[17]  David Eppstein,et al.  MESH GENERATION AND OPTIMAL TRIANGULATION , 1992 .

[18]  S. Osher,et al.  Triangle based adaptive stencils for the solution of hyperbolic conservation laws , 1992 .

[19]  D. Rose,et al.  Some errors estimates for the box method , 1987 .

[20]  T. Hughes.  Multiscale phenomena: Green's functions, the Dirichlet-to-Neumann formulation, subgrid scale models, bubbles and the origins of stabilized methods , 1995 .

[21]  B. F. Oscillator.  Large-Signal Analysis of a Silicon Read Diode Oscillator , 1969 .

[22]  Gabriel Wittum,et al.  Downwind Numbering : A Robust Multigrid Method for Convection-Diffusion Problems on Unstructured Grids , 1995 .

[23]  I. Babuska,et al.  Generalized Finite Element Methods: Their Performance and Their Relation to Mixed Methods , 1983 .