On acceleration methods for coupled nonlinear elliptic systems |
| |
Authors: | T Kerkhoven Y Saad |
| |
Institution: | (1) Department of Computer Science, University of Illinois at Urbana-Champaign, 1304 W. Springfield Ave., 61801-2987 Urbana, IL, USA;(2) Research Institute for Advanced Computer Science, NASA Ames Research Center, 94035 Moffett Field, CA, USA |
| |
Abstract: | Summary We compare both numerically and theoretically three techniques for accelerating the convergence of a nonlinear fixed point iterationu T(u), arising from a coupled elliptic system: Chebyshev acceleration, a second order stationary method, and a nonlinear version of the Generalized Minimal Residual Algorithm (GMRES) which we call NLGMR. All three approaches are implemented in Jacobian-free mode, i.e., only a subroutine which returnsT(u) as a function ofu is required.We present a set of numerical comparisons for the drift-diffusion semiconductor model. For the mappingT which corresponds to the nonlinear block Gauß-Seidel algorithm for the solution of this nonlinear elliptic system, NLGMR is found to be superior to the second order stationary method and the Chebychev acceleration. We analyze the local convergence of the nonlinear iterations in terms of the spectrum T
u
(u
(*))] of the derivativeT
u
at the solutionu
(*). The convergence of the original iteration is governed by the spectral radius T
U
(u
(*))]. In contrast, the convergence of the two second order accelerations are related to the convex hull of T
u
(u
(*))], while the convergence of the GMRES-based approach is related to the local clustering in I–T
u
(u
(*))]. The spectrum I–T
u
(u
(*))] clusters only at 1 due to the successive inversions of elliptic partial differential equations inT. We explain the observed superiority of GMRES over the second order acceleration by its ability to take advantage of this clustering feature, which is shared by similar coupled nonlinear elliptic systems. |
| |
Keywords: | 65B 65H10 65N |
本文献已被 SpringerLink 等数据库收录! |
|