Fast finite element solver for incompressible Navier-Stokes equation by parallel Gram-Schmidt process based GMRES and HSS

Di Zhao

Abstract


After finite element discretization of incompressible Navier-Stokes equation, a sparse linear system is obtained in every iteration, and GMRES provides an efficient approach to solve this system. However, if the size of the original PDE model is large, the solution speed of the incompressible Navier-Stokes equation is still slow. Since the linear system from the incompressible Navier-Stokes equation is saddle point problem, we find that large portion of computational efforts for solving the linear system is occupied by the vector projection in GMRES. In this paper, by our parallel Gram-Schmidt process based GMRES and newly developed preconditioner Hermitian/Skew-Hermitian Separation (HSS), we develop a fast solver HSS-pGMRES for the saddle point problem from incompressible Navier-Stokes Equation. Theoretical analysis shows that, HSS-pGMRES decreases the computational complexity of finite element solver for incompressible Navier-Stokes equation from O(m2n) to O(mn), where m is the grid size. Computational experiments show that, the fast solver HSS-pGMRES significantly increases the solution speed for the saddle point problem of incompressible Navier-Stokes equation than the conventional solvers.

Full Text: PDF

How to Cite this Article:

Di Zhao, Fast finite element solver for incompressible Navier-Stokes equation by parallel Gram-Schmidt process based GMRES and HSS, Journal of Mathematical and Computational Science, Vol 5, No 3 (2015), 280-296

Copyright © 2015 Di Zhao. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

J. Math. Comput. Sci.

ISSN: 1927-5307

Editorial Office: jmcs@scik.org

 

Copyright ©2017 JMCS