A Crosswind Block Iterative Method for Convection-Dominated Problems
Abstract
This paper is on a block iterative algorithm for convection-dominated equations. The algorithm uses crosswind thin blocks in a block Gauss--Seidel method. The relaxation sweep is carried out successively along the downwind direction and exact solvers are used for the block systems. This method is efficient for convection-dominated problems discretized by monotone finite element/finite difference schemes, such as the edge-average finite element method. An optimal partitioning and ordering algorithm, Tarjan's algorithm, is used to partition the nodes into crosswind blocks and order the blocks in the downwind direction. The convergence of this block iterative method is analyzed and exponential convergence rates are proved for both one- and two-dimensional cases on both structured and unstructured meshes. Our empirical and analytical studies indicate that crosswind grouping is essential for the rapid convergence of the method and merely ordering the nodes along the downwind direction is not good enough. Some numerical examples are given to illustrate the effectiveness of the proposed algorithm for convection-dominated problems and to compare it with other popular algorithms.
MSC codes
Keywords
Get full access to this article
View all available purchase options and get full access to this article.
References
1.
Jürgen Bey, Gabriel Wittum, Downwind numbering: robust multigrid for convection‐diffusion problems, Appl. Numer. Math., 23 (1997), 177–192, Multilevel methods (Oberwolfach, 1995)
2.
A. Brandt, I. Yavneh, On multigrid solution of high‐Reynolds incompressible entering flows, J. Comput. Phys., 101 (1992), 151–164
3.
P. G. Ciarlet, The Finite Element Method for Elliptic Problems, North–Holland, Amsterdam, New York, 1978.
4.
Udi Manber, Introduction to algorithms, Addison‐Wesley Publishing Company Advanced Book Program, 1989xiv+478, A creative approach
5.
I. Duff, A. Erisman, J. Reid, Direct methods for sparse matrices, Monographs on Numerical Analysis, The Clarendon Press Oxford University Press, 1989xiv+341, Oxford Science Publications
6.
Howard Elman, Michael Chernesky, Ordering effects on relaxation methods applied to the discrete one‐dimensional convection‐diffusion equation, SIAM J. Numer. Anal., 30 (1993), 1268–1290
7.
Howard Elman, Michael Chernesky, Ordering effects on relaxation methods applied to the discrete convection‐diffusion equation, IMA Vol. Math. Appl., Vol. 60, Springer, New York, 1994, 45–57
8.
Paul Farrell, Flow conforming iterative methods for convection dominated flows, IMACS Ann. Comput. Appl. Math., Vol. 1, Baltzer, Basel, 1989, 681–686
9.
Alan George, Joseph Liu, Computer solution of large sparse positive definite systems, Prentice‐Hall Inc., 1981xii+324, Prentice‐Hall Series in Computational Mathematics
10.
Wolfgang Hackbusch, Thomas Probst, Downwind Gauss‐Seidel smoothing for convection dominated problems, Numer. Linear Algebra Appl., 4 (1997), 85–102
11.
Hou Han, V. Il’in, R. Kellogg, Wei Yuan, Analysis of flow directed iterations, J. Comput. Math., 10 (1992), 57–76
12.
I. Munro, Efficient determination of the transitive closure of a directed graph, Inform. Process. Lett., 1 (1971), pp. 55–58.
13.
I. Munro, Some Results in the Study of Algorithms, Ph.D. thesis, University of Toronto, Toronto, Ontario, Canada, 1971.
14.
J. Oliger, W. Skamarock, and W.‐P. Tang, Schwarz Alternating Methodand Its SOR Accelerations, Tech. Report, Computer Science Dept., Stanford University, Stanford, CA, 1986.
15.
James O’Neil, Daniel Szyld, A block ordering method for sparse matrices, SIAM J. Sci. Statist. Comput., 11 (1990), 811–823
16.
R. Read, Graph Theory and Computing, Academic Press, New York and London, 1972.
17.
R. W. H. Sargent and A. W. Westerberg, Speed‐up in chemical engineering design, Trans. Inst. Chem. Engrg., 42 (1964), pp. 190–197.
18.
D. Scharffetter and H. Gummel, Large‐signal analysis of a silicon read diode oscillator, IEEE Trans. Electron Devices, ED‐16 (1969), pp. 64–77.
19.
S. Selberherr, Analysis and Simulation of Semiconductor Devices, Springer‐Verlag, New York, 1984.
20.
W.‐P. Tang, Schwarz Splitting and Template Operators, Ph.D. thesis, Stanford University, Stanford, CA, 1987.
21.
R. E. Tarjan, Depth‐first search and linear graph algorithms, SIAM J. Comput., 1 (1972), pp. 146–160.
22.
Robert Tarjan, Efficiency of a good but not linear set union algorithm, J. Assoc. Comput. Mach., 22 (1975), 215–225
23.
Richard Varga, Matrix iterative analysis, Prentice‐Hall Inc., 1962xiii+322
24.
J. Xu, The EAFE scheme and CWS method for convection dominated problems, in Proceedings 9th International Symposium on Domain Decomposition Methods in Science and Engineering, P. Bjorstad, M. Espedal, and D. Keyes, eds., John Wiley & Sons Ltd., Chichester, UK, 1997.
25.
Jinchao Xu, Ludmil Zikatanov, A monotone finite element scheme for convection‐diffusion equations, Math. Comp., 68 (1999), 1429–1446
Information & Authors
Information
Published In
SIAM Journal on Scientific Computing
Pages: 620 - 645
ISSN (online): 1095-7197
Copyright
Copyright © 1999 Society for Industrial and Applied Mathematics.
History
Published online: 25 July 2006
MSC codes
Keywords
Authors
Metrics & Citations
Metrics
Citations
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.
Cited By
- Stabilized finite element discretizations of general convection-diffusion problemsSCIENTIA SINICA Mathematica, Vol. 54, No. 1 | 3 November 2023
- A multilevel approach for trace system in HDG discretizationsJournal of Computational Physics, Vol. 407 | 1 Apr 2020
- High-order implicit palindromic discontinuous Galerkin method for kinetic-relaxation approximationComputers & Fluids, Vol. 190 | 1 Aug 2019
- iHDG: An Iterative HDG Framework for Partial Differential EquationsSIAM Journal on Scientific Computing, Vol. 39, No. 5 | 26 October 2017
- Numerical studies of a class of linear solvers for fine-scale petroleum reservoir simulationComputing and Visualization in Science, Vol. 18, No. 2-3 | 21 December 2016
- Cost-Effective Parallel Reservoir Simulation on Shared MemoryAll Days | 25 October 2016
- A multilevel preconditioner and its shared memory implementation for a new generation reservoir simulatorPetroleum Science, Vol. 11, No. 4 | 4 October 2014
- Multilevel Preconditioners for a New Generation Reservoir SimulatorAll Days | 16 September 2013
- The single-grid multilevel method and its applicationsInverse Problems & Imaging, Vol. 7, No. 3 | 1 Jan 2013
- Parallel Implementation of the Gauss-Seidel Algorithm on k-Ary n-Cube MachineApplied Mathematics, Vol. 04, No. 01 | 1 Jan 2013
- A Block Solver for the Exponentially Fitted IIPG-0 MethodDomain Decomposition Methods in Science and Engineering XX | 9 May 2013
- PARALLEL GAUSS-SEIDEL ON A TORUS NETWORK-ON-CHIP ARCHITECTUREJournal of Interconnection Networks, Vol. 13, No. 01n02 | 10 March 2013
- Towards a New Generation of Physics-Driven Solvers for Black-Oil and Compositional Flow SimulationAll Days | 2 February 2009
- Scalable Heuristic Algorithms for the Parallel Execution of Data Flow Acyclic DigraphsSIAM Journal on Scientific Computing, Vol. 31, No. 5 | 25 September 2009
- Fast computation of multiphase flow in porous media by implicit discontinuous Galerkin schemes with optimal ordering of elementsJournal of Computational Physics, Vol. 227, No. 24 | 1 Dec 2008
- Kernel Preserving Multigrid Methods for Convection-Diffusion EquationsSIAM Journal on Matrix Analysis and Applications, Vol. 27, No. 4 | 31 July 2006
- A Phase Error Analysis of Multigrid Methods for Hyperbolic EquationsSIAM Journal on Scientific Computing, Vol. 25, No. 3 | 25 July 2006
- Singularly perturbed problems in partial differential equations: a surveyApplied Mathematics and Computation, Vol. 134, No. 2-3 | 1 Jan 2003
View Options
Get Access
- Access via your Institution
- Questions about how to access this content? Contact SIAM at [email protected].