On acceleration methods for coupled nonlinear elliptic systems (Q758154)

From MaRDI portal
scientific article
Language Label Description Also known as
English
On acceleration methods for coupled nonlinear elliptic systems
scientific article

    Statements

    On acceleration methods for coupled nonlinear elliptic systems (English)
    0 references
    0 references
    0 references
    0 references
    1992
    0 references
    We compare both numerically and theoretically three techniques for accelerating the convergence of a nonlinear fixed point iteration \(u\to 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 returns T(u) as a function of u is required. We present a set of numerical comparisons for the drift-diffusion semiconductor model. For the mapping T which corresponds to the nonlinear block Gauss-Seidel algorithm for the solution of this nonlinear elliptic system, NLGMR is found to be superior to the second order stationary method and the Chebyshev acceleration. We analyze the local convergence of the nonlinear iterations in terms of the spectrum \(\sigma [T_ u(u^{(*)})]\) of the derivative \(T_ u\) at the solution \(u^{(*)}\). The convergence of the original iteration is governed by the spectral radius \(\rho [T_ u(u^{(*)})].\) In contrast, the convergence of the two second order accelerations is related to the convex hull of \(\sigma [T_ u(u^{(*)})]\), while the convergence of the GMRES-based approach is related to the local clustering in \(\sigma [I-T_ u(u^{(*)})]\). The spectrum \(\sigma [I- T_ u(u^{(*)})]\) clusters only at 1 due to the successive inversions of elliptic partial differential equations in T. 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.
    0 references
    0 references
    0 references
    0 references
    0 references
    0 references
    0 references
    0 references
    0 references
    convergence acceleration
    0 references
    compactness
    0 references
    matrix-free methods
    0 references
    nonlinear fixed point iteration
    0 references
    Chebyshev acceleration
    0 references
    Generalized Minimal Residual Algorithm
    0 references
    numerical comparisons
    0 references
    drift-diffusion semiconductor model
    0 references
    nonlinear block Gauss-Seidel algorithm
    0 references
    nonlinear elliptic system
    0 references
    second order stationary method
    0 references
    local convergence
    0 references
    0 references
    0 references
    0 references
    0 references
    0 references
    0 references
    0 references
    0 references
    0 references