An element-based spectrally-optimized approximate inverse preconditioner for the Euler equations
Abstract
We introduce a method for constructing an element-by-element sparse approximate inverse (SAI) preconditioner designed to be effective in a massively-parallel spectral element modeling environment involving non- symmetric systems. This new preconditioning approach is based on a spectral optimization of a low-resolution pre- conditioned system matrix (PSM). We show that the local preconditioning matrices obtained via this element-based, spectrum-optimized (EBSO) approach may be applied to arbitrarily high-resolution versions of the same system matrix without appreciable loss of preconditioner performance. We demonstrate the performance of the EBSO precondition- ing approach using 2-D spectral element method (SEM) formulations for a simple linear conservation law and for the fully-compressible 2-D Euler equations with various boundary conditions. For the latter model running at suffi- ciently large Courant number, the EBSO preconditioner significantly reduces both iteration count and wall-clock time regardless of whether a generalized minimum residual (GMRES) or a stabilized biconjugate gradient (BICGSTAB) iterative scheme is employed. To assess the value added by this new preconditioning approach, we compare its perfor- mance against two other equally-parallel SAI preconditioning methods: low-order Chebyshev generalized least-squares polynomials and an element-based variant of the well-known Frobenius norm optimization preconditioner which we also develop herein. The EBSO preconditioner significantly out-performs both the Chebyshev polynomials and the element-based Frobenius-norm-optimized (EBFO) preconditioner regardless of whether the GMRES or BICGSTAB iterative scheme is employed. Moreover, when the EBSO preconditioner is combined with the Chebyshev polynomial method dramatic reductions in iterations per time-step can be achieved while still achieving a significant reduction in wall-clock time.
Description
Non-hydrostatic Unified Model of the Atmosphere (NUMA)
The first NUMA papers appeared in 2008. From 2008 through 2010, all the NUMA papers appearing involved the 2D (x-z slice) Euler equations. All the theory and numerical implementations were first developed in 2D.
Rights
This publication is a work of the U.S. Government as defined in Title 17, United States Code, Section 101. Copyright protection is not available for this work in the United States.Collections
Related items
Showing items related by title, author, creator and subject.
-
A polynomial-based nonlinear least squares optimized preconditioner for continuous and discontinuous element-based discretizations of the Euler equations
Carr, L.E., III; Borges, C.F.; Giraldo, F.X. (SIAM, 2014);We introduce a method for constructing a polynomial preconditioner using a nonlinear least squares (NLLS) algorithm. We show that this polynomial-based NLLS-optimized (PBNO) preconditioner significantly im- proves the ... -
Matrix-Free Polynomial-Based Nonlinear Least Squares Optimized Preconditioning and its Application to Discontinuous Galerkin Discretizations of the Euler Equations
Carr, L.E.; Borges, C.F.; Giraldo, F.X. (2015-02-12);We introduce a preconditioner that can be both constructed and applied using only the ability to apply the underlying operator. Such a preconditioner can be very attractive in scenarios where one has a highly efficient ... -
Preconditioner circuit analysis
Nye, Matthew J. (Monterey, California. Naval Postgraduate School, 2011-09);Voltages up to 10,000 volts or higher must be attenuated and measured to provide control feedback for many applications like medium voltage generators or pulsed power systems. How these medium voltage signals can be ...