We study the performance of the multigrid method applied to spectral element (SE) discretizations of the Poisson and Helmholtz equations. Smoothers based on finite element (FE) discretizations, overlapping Schwarz methods, and point-Jacobi are considered in conjunction with conjugate gradient and GMRES acceleration techniques. It is found that Schwarz methods based on restrictions of the originating SE matrices converge faster than FE-based methods and that weighting the Schwarz matrices by the inverse of the diagonal counting matrix is essential to effective Schwarz smoothing. Several of the methods considered achieve convergence rates comparable to those attained by classic multigrid on regular grids.
We describe recent algorithmic developments that have enabled large eddy simulations of reactor flows on up to P = 65, 000 processors on the IBM BG/P at the Argonne Leadership Computing Facility.
An overlapping mesh methodology that is spectrally accurate in space and up to third-order accurate in time is developed for solution of unsteady incompressible flow equations in three-dimensional domains. The ability to decompose a global domain into separate, but overlapping, subdomains eases mesh generation procedures and increases flexibility of modeling flows with complex geometries. The methodology employs implicit spectral element discretization of equations in each subdomain and explicit treatment of subdomain interfaces with spectrallyaccurate spatial interpolation and high-order accurate temporal extrapolation, and requires few, if any, iterations, yet maintains the global accuracy and stability of the underlying flow solver. The overlapping mesh methodology is thoroughly validated using two-dimensional and three-dimensional benchmark problems in laminar and turbulent flows. The spatial and temporal convergence is documented and is in agreement with the nominal order of accuracy of the solver. The influence of long integration times, as well as inflow-outflow global boundary conditions on the performance of the overlapping grid solver is assessed. In a turbulent benchmark of fully-developed turbulent pipe flow, the turbulent statistics with the overlapping grids is validated against published available experimental and other computation data. Scaling tests are presented that show near linear strong scaling, even for moderately large processor counts.
Summary. The performance of multigrid methods for the standard Poisson problem and for the consistent Poisson problem arising in spectral element discretizations of the Navier-Stokes equations is investigated. It is demonstrated that overlapping additive Schwarz methods are effective smoothers, provided that the solution in the overlap region is weighted by the inverse counting matrix. It is also shown that spectral element based smoothers are superior to those based upon finite element discretizations. Results for several large 3D Navier-Stokes applications are presented.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.