Abstract
It was observed that all the different linear systems arising in an iterative fluid flow simulation algorithm have approximately constant invariant subspaces associated with their smallest eigenvalues. For this reason, we propose to perform one single computation of the eigenspace associated with the smallest eigenvalues, at the beginning of the iterative process, to improve the convergence of the Krylov method used in subsequent iterations of the fluid flow algorithm by means of this pre-computed partial spectral information. The Subspace Inverse Iteration Method with Stabilized Block Conjugate Gradient is our choice for computing the spectral information, which is then used to remove the effect of the smallest eigenvalues in two different ways: either building a spectral preconditioner that shifts these eigenvalues from almost zero close to the unit value, or performing a deflation of the initial residual in order to remove parts of the solution corresponding to the smallest eigenvalues. Under certain conditions, both techniques yield a reduction of the number of iterations in each subsequent runs of the Conjugate Gradient algorithm.
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
Preview
Unable to display preview. Download preview PDF.
References
Patankar, S.V., Spalding, B.B.: A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows. International Journal of Heat and Mass Transfer 15, 1787–1806 (1972)
Golub, G.H., Loan, C.F.V.: Matrix Computation. The Johns Hopkins University Press, Baltimore (1983)
Rhie, C.M., Chow, W.L.: Numerical study of the turbulent flow past an airfoil with trailing edge separation. AIAA Journal 21, 1525–1532 (1983)
Ferziger, J.H., Perić, M.: Computational Methods for Fluid Dynamics, 3rd edn. Springer, Heidelberg (2002)
McGuirk, J.J., Palma, J.M.L.M.: The efficiency of alternative pressure-correction formulations for incompressible turbulent flow problems. Computers & Fluids 22, 77–87 (1993)
Lehoucq, R., Sorensen, D., Yang, C.: Arpack users’ guide: Solution of large scale eigenvalue problems with implicitly restarted arnoldi methods (1997)
Arioli, M., Ruiz, D.: Block conjugate gradient with subspace iteration for solving linear systems. In: Margenov, S., Vassilevski, P. (eds.) Iterative Methods in Linear Algebra, Second IMACS Symposium on Iterative Metohds in Linear Algebra, Blagoevgrad, Bulgaria, June 17-20, pp. 64–79 (1995)
Arioli, M., Duff, I., Ruiz, D., Sadkane, M.: Block lanczos techniques for accelerating the block cimmino method. Journal on Scientific Computing 16, 1478–1511 (1995)
Van der Sluis, A., van der Vorst, H.A.: The rate of convergence of Conjugate Gradients. Numerical Mathematics 48, 543–560 (1986)
Carpentieri, B., Duff, I.S., Giraud., L.: A class of spectral two-level preconditioners. Technical Report TR/PA/02/55, CERFACS, Also Technical Report RAL-TR-2002-020. To appear in SIAM J. Scientific Computing (2002)
Rigal, J.L., Gaches, J.: On the compatibility of a given solution with the data of a linear system. Journal of the ACM 14, 543–548 (1967)
Parlett, B.N.: The Symmetric Eigenvalue Problem. SIAM, Philadelphia (1998)
Arioli, M.: A stopping criterion for conjugate gradient algorithm in a finite element method. Numerical Mathematics, 1–24 (2003) DOI 10.1007/s00211-003-0500-y
Giraud, L., Ruiz, D., Touhami, A.: A comparative study of iterative solvers exploiting spectral information for SPD systems. Technical Report RT/TLSE/04/03, ENSEEIHT - IRIT (2004)
Dongarra, J., Decroz, J., Hammarlung, S., Duff, I.: A set of level 3 basic linear algebra subprograms. ACM Trans. Math. Software 16, 1–28 (1990)
Author information
Authors and Affiliations
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 2005 Springer-Verlag Berlin Heidelberg
About this paper
Cite this paper
Balsa, C., Palma, J.M.L.M., Ruiz, D. (2005). Partial Spectral Information from Linear Systems to Speed-Up Numerical Simulations in Computational Fluid Dynamics. In: Daydé, M., Dongarra, J., Hernández, V., Palma, J.M.L.M. (eds) High Performance Computing for Computational Science - VECPAR 2004. VECPAR 2004. Lecture Notes in Computer Science, vol 3402. Springer, Berlin, Heidelberg. https://doi.org/10.1007/11403937_52
Download citation
DOI: https://doi.org/10.1007/11403937_52
Publisher Name: Springer, Berlin, Heidelberg
Print ISBN: 978-3-540-25424-9
Online ISBN: 978-3-540-31854-5
eBook Packages: Computer ScienceComputer Science (R0)