论文信息 - A Parallel Algorithm for Systems of Convection-Diffusion Equations

A Parallel Algorithm for Systems of Convection-Diffusion Equations

The numerical solution of systems of convection-diffusion equations is considered. The problem is described by a system of second order partial differential equations (PDEs). This system is discretized by Courant-elements. The preconditioned conjugate gradient method is used for the iterative solution of the large-scale linear algebraic systems arising after the finite element discretization of the problem. Discrete Helmholtz preconditioners are applied to obtain a mesh independent superlinear convergence of the iterative method. A parallel algorithm is derived for the proposed preconditioner. A portable parallel code using Message Passing Interface (MPI) is developed. Numerical tests well illustrate the performance of the proposed method on a parallel computer architecture.

[1] Thomas A. Manteuffel,et al. Optimal equivalent preconditioners , 1993 .

[2] O. Axelsson. Iterative solution methods , 1995 .

[3] Tamás Kurics,et al. Superlinearly convergent PCG algorithms for some nonsymmetric elliptic systems , 2008 .

[4] J. Karátson,et al. Symmetric Part Preconditioning for the Conjugate Gradient Method in Hilbert Space , 2003 .

[5] Thomas A. Manteuffel,et al. On the theory of equivalent operators and application to the numerical solution of uniformly elliptic partial differential equations , 1990 .

[6] Jack Dongarra,et al. LAPACK Users' Guide, 3rd ed. , 1999 .

[7] Jack Dongarra,et al. MPI: The Complete Reference , 1996 .

[8] Zahari Zlatev,et al. Computer Treatment of Large Air Pollution Models , 1995 .

[9] James Demmel,et al. LAPACK Users' Guide, Third Edition , 1999, Software, Environments and Tools.

[10] O. Axelsson. A generalized conjugate gradient, least square method , 1987 .

[11] Owe Axelsson,et al. Superlinearly convergent CG methods via equivalent preconditioning for nonsymmetric elliptic operators , 2004, Numerische Mathematik.