Abstract
Standard Galerkin finite element methods or finite difference methods for singular perturbation problems lead to strongly unsymmetric matrices, which furthermore are in general notM-matrices. Accordingly, preconditioned iterative methods such as preconditioned (generalized) conjugate gradient methods, which have turned out to be very successful for symmetric and positive definite problems, can fail to converge or require an excessive number of iterations for singular perturbation problems.
This is not so much due to the asymmetry, as it is to the fact that the spectrum can have both eigenvalues with positive and negative real parts, or eigenvalues with arbitrary small positive real parts and nonnegligible imaginary parts. This will be the case for a standard Galerkin method, unless the meshparameterh is chosen excessively small. There exist other discretization methods, however, for which the corresponding bilinear form is coercive, whence its finite element matrix has only eigenvalues with positive real parts; in fact, the real parts are positive uniformly in the singular perturbation parameter.
In the present paper we examine the streamline diffusion finite element method in this respect. It is found that incomplete block-matrix factorization methods, both on classical form and on an inverse-free (vectorizable) form, coupled with a general least squares conjugate gradient method, can work exceptionally well on this type of problem. The number of iterations is sometimes significantly smaller than for the corresponding almost symmetric problem where the velocity field is close to zero or the singular perturbation parameter ε=1.
Similar content being viewed by others
References
J. Aarden,The numerical solution of non-linear and non-stationary convection-diffusion equations, M.Sc. Thesis, Department of Mathematics, University of Nijmegen, December 1987 (In Dutch).
O. Axelsson,On the numerical solution of convection dominated convection-diffusion problems, Mathematical Methods in Energy Research (K. I. Gross, editor), SIAM, Philadelphia, 1984, pp. 3–21.
O. Axelsson,Incomplete block matrix factorization preconditioning methods, The ultimate answer?, J. Comp. Appl. Math. 12/13 (1985), pp. 3–18.
O. Axelsson,Stability and error estimates of Galerkin finite element approximations for convection-diffusion equations, IMA J. Numerical Analysis 1, 1981, 329–345.
O. Axelsson,A generalized conjugate gradient, least square method, Numer. Math. 51 (1987), pp. 209–227.
O. Axelsson and V. Eijkhout,A note on the vectorization of scalar recursions, Parallel Computing, 3 (1986), pp. 73–83.
O. Axelsson and V. Eijkhout,Robust vectorizable preconditioners for three-dimensional elliptic difference equations with anisotropy, Algorithms and Applications on Vector and Parallel Computers (H. J. J. te Riele, Th. J. Dekker and H. A. van der Vorst, editors), pp. 279–306, North Holland, 1987.
O. Axelsson and I. Gustafsson,A modified upwind scheme for convective transport equations and the use of conjugate gradient method for the solution of non-symmetric systems of equations, J. Inst. Math. Appl. 23 (1979), pp. 321–327.
O. Axelsson and W. Layton,Defect correction methods for convection-diffusion problems, Report 8335, Department of Mathematics, University of Nijmegen.
O. Axelsson and B. Polman,On approximate factorization methods for block-matrices suitable for vector and parallel processors, Lin. Alg. Appl. 77 (1986), 3–26.
O. Axelsson,A generalized conjugate direction method and its application to a singular perturbation problem, Proc. 8th biennial Numerical Analysis Conference (Dundee, Scotland JUne 26–29, 1979). (G. A. Watson, ed.), Lecture Notes in Math. vol. 773, Springer-Verlag, Berlin, 1980, 35–357.
O. Axelsson,A restarted version of a generalised preconditioned conjugate gradient method, Communications in Applied Numerical Methods, 4 (1988), pp. 521–530.
O. Axelsson,Finite element methods for convection diffusion problems, Proceedings, Conference on Numerical Solution of Differential Equations, Halle, DDR, May 25–29, 1987, pp. 171–183.
P. Concus, G. H. Golub and G. Meurant,Block preconditioning for the conjugate gradient method, SIAM, J. Sci. Stat. Comput., 6 (1985) pp. 220–252.
J. E. Dennis and K. Turner,Generalized Conjugate Directions, Lin. Alg. Appl. 88/89 (1987), pp. 187–209.
S. C. Eisenstat, H. C. Elman and M. H. Schultz,Variational iterative methods for nonsymmeric systems of linear equations, SIAM J. Num. Anal. 20 (1983) pp. 345–357.
H. C. Elman,Iterative methods for large sparse nonsymmetric systems of linear equations, Ph.D. Thesis, Technical Report 229, Yale University, New Haven, C.T. 1982.
L. Hageman and D. M. Young,Applied Iterative Methods, Academic press, New York, 1981.
T. J. Hughes and A. Brooks,A multidimensional upwind scheme with no crosswind diffusion, AMD vol. 34, Finite element methods for convection dominated flows, T. J. Hughes (ed.), AMSE, New York, 1979.
T. J. R. Hughes, M. Mallet and Akira Mizukami,A new finite element formulation for computational fluid dynamics: II beyond SUPG, Computer Methods in Applied Mechanics and Engineering 54, 1986, pp. 341–355.
C. Johnson and U. Nävert,An analysis of some finite element methods for advection-diffusion problems, O. Axelsson, L. S. Frank and A. van der Sluis (eds.),Analytical and Numerical Approaches to Asymptotic Problems in Analysis, North Holland, 1981, pp. 99–116.
U. Nävert,A finite element method for convection diffusion problems, Department of Computer Sciences, Chalmers University of Technology, Göteborg, Sweden, (Ph.D. Thesis), 1982.
J. M. Ortega and R. G. Voigt,Solution of Partial Differential Equations on Vector and Parallel Computers, SIAM, Philadelphia, P.A., 1985.
Y. Saad and M. H. Schultz,Conjugate gradient-like algorithms for solving nonsymmetric linear systems, Math. of Comp. 44 (1985) pp. 417–424.
Author information
Authors and Affiliations
Additional information
The 2nd author's research was sponsored by Control Data Corporation through its PACER fellowship program.
The 3rd author's research was supported by the Netherlands organization for scientific research (NWO).
On leave from the Institute of Mathematics, Academy of Science, 1090 Sofia, P.O. Box 373, Bulgaria.
Rights and permissions
About this article
Cite this article
Axelsson, O., Eijkhout, V., Polman, B. et al. Incomplete block-matrix factorization iterative methods for convection-diffusion problems. BIT 29, 867–889 (1989). https://doi.org/10.1007/BF01932750
Received:
Revised:
Issue Date:
DOI: https://doi.org/10.1007/BF01932750