Iterative residual-based vector methods to accelerate fixed point iterations. (November 2015)

• Record Type: Journal Article
• Authors: Ramière, Isabelle; Helfer, Thomas
• Is Part Of: Computers & mathematics with applications. Volume 70:issue 9(2015)
• Page Start: 2210
• Page End: 2226
• Publication Date: 2015-11
• Subjects: Electronic data processing -- Periodicals; Mathematics -- Data processing -- Periodicals; 510.28541
• DOI: 10.1016/j.camwa.2015.08.025
• Languages: English
• ISSNs: 0898-1221

Abstract

<abstract xml:lang="en" abstract-type="author" id="a000005"> <title id="st000005">Abstract</title> <sec> <p id="sp000090">Fixed point iterations are still the most common approach to dealing with a variety of numerical problems such as coupled problems (multi-physics, domain decomposition, …) or nonlinear problems (electronic structure, heat transfer, nonlinear mechanics, …). Methods to accelerate fixed point iteration convergence or more generally sequence convergence have been extensively studied since the 1960's. For scalar sequences, the most popular and efficient acceleration method remains the <inline-formula><alternatives><inline-graphic xlink:href="ark:/27927/pgkfv4r8f0" xlink:type="simple" xmlns:xlink="http://www.w3.org/1999/xlink" /><mml:math altimg="si33.gif" display="inline" overflow="scroll" id="d13e830" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msup><mml:mrow><mml:mi>Δ</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:math></alternatives></inline-formula> of Aitken. Various vector acceleration algorithms are available in the literature, which often aim at being multi-dimensional generalizations of the <inline-formula><alternatives><inline-graphic xlink:href="ark:/27927/pgkfv4r8f0" xlink:type="simple" xmlns:xlink="http://www.w3.org/1999/xlink" /><mml:math altimg="si33.gif" display="inline" overflow="scroll" id="d13e839" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msup><mml:mrow><mml:mi>Δ</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:math></alternatives></inline-formula> method.</p> <p id="sp000095">In this paper, we propose and analyse a generic residual-based formulation for accelerating vector sequences. The question of the dynamic use of this residual-based transformation during the fixed point iterations for obtaining a new accelerated fixed point method is then raised. We show that two main classes of such iterative algorithms can be derived and that this approach is generic in that various existing acceleration algorithms for vector sequences are thereby recovered.</p> <p id="sp000100">In order to illustrate the interest of such algorithms, we apply them in the field of nonlinear mechanics on a simplified "point-wise" solver used to perform mechanical behaviour unit testings. The proposed test cases clearly demonstrate that accelerated fixed point iterations based on the elastic operator (quasi-Newton method) are very useful when the mechanical behaviour does not provide the so-called consistent tangent operator. Moreover, such accelerated algorithms also prove to be competitive with respect to the standard Newton–Raphson algorithm when available.</p> </sec> </abstract>