Faculté des sciences

Stability of the standard Crank–Nicolson–Galerkin scheme applied to the diffusion-convection equation: Some new insights

Perrochet, Pierre ; Bérod, Dominique

In: Water Ressources Research, 1993, vol. 29, no. 9, p. 3291–3298

A stability analysis of the classical Crank–Nicolson–Galerkin (CNG) scheme applied to the one-dimensional solute transport equation is proposed on the basis of two fairly different approaches. Using a space-time eigenvalue analysis, it is shown, at least for subsurface hydrology applications, that the CNG scheme is theoretically stable under the condition PeCr≤2, where Pe and... Plus

Ajouter à la liste personnelle
    Summary
    A stability analysis of the classical Crank–Nicolson–Galerkin (CNG) scheme applied to the one-dimensional solute transport equation is proposed on the basis of two fairly different approaches. Using a space-time eigenvalue analysis, it is shown, at least for subsurface hydrology applications, that the CNG scheme is theoretically stable under the condition PeCr≤2, where Pe and Cr are the mesh Péclet and Courant numbers. Then, to assess the computational stability of the scheme, the amplification matrix is constructed, and its norm is displayed in the (Pe, Cr) space. The results indicate that the norm of the amplification matrix is never smaller than unity and exhibits a hyperbolic nature analogous to the above theoretical condition.