Stability of the standard Crank‐Nicolson‐Galerkin Scheme applied to the diffusion‐convection equation: Some new insights

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 Peclet 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.

[1]  W. G. Gray,et al.  Finite Element Simulation in Surface and Subsurface Hydrology , 1977 .

[2]  Richard S. Varga,et al.  Application of Oscillation Matrices to Diffusion-Convection Equations , 1966 .

[3]  Joannes J. Westerink,et al.  Consistent higher degree Petrov–Galerkin methods for the solution of the transient convection–diffusion equation , 1989 .

[4]  P. S. Huyakorn,et al.  Solution of transient transport equation using an upstream finite element scheme , 1979 .

[5]  L. Quartapelle,et al.  An analysis of time discretization in the finite element solution of hyperbolic problems , 1987 .

[6]  Juan C. Heinrich,et al.  Petrov-Galerkin method for multidimensional time-dependent convective-diffusion equations , 1987 .

[7]  Chin-Fu Tsang,et al.  A perspective on the numerical solution of convection‐dominated transport problems: A price to pay for the easy way out , 1992 .

[8]  O. C. Zienkiewicz,et al.  A note on upwinding and anisotropic balancing dissipation in finite element approximations to convective diffusion problems , 1980 .

[9]  W. Gray,et al.  An analysis of the numerical solution of the transport equation , 1976 .

[10]  Emil O. Frind,et al.  Comparative error analysis in finite element formulations of the advection-dispersion equation , 1985 .

[11]  P. Perrochet.  A streamline-upwind-full-Galerkin method for space-time convection dominated transport problems , 1993 .

[12]  G. Pinder,et al.  Computational Methods in Subsurface Flow , 1983 .

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