Open AccessA conservative, unconditionally stable, second‐order three‐point differencing scheme for the diffusion–convection equation
Author(s) -
Clauser Christoph,
Kiesner Stefan
Publication year - 1987
Publication title -
geophysical journal of the royal astronomical society
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.302
H-Index - 168
eISSN - 1365-246X
pISSN - 0016-8009
DOI - 10.1111/j.1365-246x.1987.tb01658.x
Subject(s) - central differencing scheme , convection–diffusion equation , scheme (mathematics) , numerical diffusion , numerical solution of the convection–diffusion equation , diffusion , mathematics , diffusion equation , convection , point (geometry) , upwind scheme , numerical analysis , mathematical analysis , computer science , physics , mechanics , geometry , discretization , finite element method , finite difference coefficient , economy , mixed finite element method , economics , thermodynamics , service (business)
A three‐point differencing scheme for the diffusion–convection equation is presented that offers all the advantages of both the central and the one‐sided ('upwind') differencing scheme without suffering from their drawbacks. Specifically, the scheme is conservative, unconditionally stable, and second‐order‐accurate in space. It is free of oscillations and over‐ or undershoots, simple to code, and requires essentially no more computing time than the one‐sided scheme. Although known for a relatively long time in numerical mathematics, the scheme apparently has not received sufficient attention from modellers of hydrothermal systems or contaminant transport in the geosciences. In order to fill this gap a comparison is made between this scheme and the widely used one‐sided scheme for the transient diffusion–convection equation in different time discretizations. The results are discussed taking into account other approaches towards minimizing numerical diffusion.