Robust Multisecant Quasi-Newton Variants for Parallel Fluid-Structure Simulations---and Other Multiphysics Applications


Multisecant quasi-Newton methods have been shown to be particularly suited to solve nonlinear fixed-point equations that arise from partitioned multiphysics simulations where the exact Jacobian is inaccessible. In all these methods, the underdetermined multisecant equation for the approximate (inverse) Jacobian is enhanced by a norm minimization condition. The standard choice is the minimization of the Frobenius norm of the approximate inverse Jacobian. In this setting, it is well known that transient fluid-structure simulations typically require the use of secant information also from previous time steps to achieve a small enough number of iterations per implicit time step. The number of these time steps highly depends on the application, the physical parameters, the used solvers, and the mesh resolution. Using too few leads to a relatively high number of iterations, while using too many leads not only to a computational overhead but also to an increase in the number of iterations as well. Determining the optimal number requires a costly trial-and-error process. In this paper, we present results for two different approaches to overcome this issue: The first approach is to use a modified method (presented in [F. Lindner, M. Mehl, K. Scheufele, and B. Uekermann, “A Comparison of Various Quasi-Newton Schemes for Partitioned Fluid-Structure Interaction,” in Proceedings of ECCOMAS Coupled Problems, Venice, 2015, pp. 1--12]) that minimizes the Frobenius norm of the difference between the current approximate (inverse) Jacobian and that of the previous time step. Thus, previous time step information is taken into account in an implicit and automatized way without magic parameters. The second approach is to use an unrestricted number of previous time steps in combination with a suitable filtering algorithm automatically removing secant information that is outdated (thus, slowing down convergence) or contradicting newer information (deteriorating the condition of the multisecant equation system). We present a novel algorithm to realize the first idea with linear complexity in the number of coupling surface unknowns (note that already storing the approximate inverse Jacobian would induce quadratic complexity) and the efficient parallel implementation for both approaches. This results in highly efficient, parallelizable, and robust iterative solvers applicable for surface coupling in many types of multiphysics simulations using black-box solver software. In addition, our numerical results for the fluid-structure benchmark (FSI3) from Turek and Hron [``Proposal for Numerical Benchmarking of Fluid-Structure Interaction between an Elastic Object and Laminar Incompressible Flow,” in Fluid-Structure Interaction, Springer, Berlin, 2006, pp. 371--385] and for a flow through a flexible tube with a large variety of parameter settings prove the robustness and numerical efficiency of the first approach in particular. The second approach can be shown to be highly sensitive to the choice of the filtering method for the secant information and not always robust for the filters we currently use.


  1. partitioned multiphysics
  2. nonlinear fixed-point solver
  3. quasi-Newton
  4. fluid-structure interactions

MSC codes

  1. 65B99
  2. 65H10
  3. 65M99
  4. 68W10
  5. 74B99
  6. 76D05

Get full access to this article

View all available purchase options and get full access to this article.


D. G. Anderson, Iterative procedures for nonlinear integral equations, J. Assoc. Comput. Mach., 12 (1965), pp. 547--560.
K.-J. Bathe and G. A. Ledezma, Benchmark problems for incompressible fluid flows with structural interactions, Comput. & Structures, 85 (2007), pp. 628--644, https://doi.org/10.1016/j.compstruc.2007.01.025.
D. S. Blom, A. H. van Zuijlen, and H. Bijl, Acceleration of strongly coupled fluid-structure interaction with manifold mapping, in Proceedings of the 11th World Congress on Computational Mechanics, 5th European Congress on Computational Mechanics, 6th European Congress on Computational Fluid Dynamics, E. On͂ate, X. Oliver, and A. Huerta, eds., 2014, pp. 4484--4495.
A. E. J. Bogaers, S. Kok, B. D. Reddy, and T. Franz, Quasi-Newton methods for implicit black-box FSI coupling, Comput. Methods Appl. Mech. Engrg., 279 (2014), pp. 113--132.
M. Brand, Fast low-rank modifications of the thin singular value decomposition, Linear Algebra Appl., 415 (2006), pp. 20--30, https://doi.org/10.1016/j.laa.2005.07.021.
C. G. Broyden, A class of methods for solving nonlinear simultaneous equations, Math. Comp., 19 (1965), pp. 577--593.
H.-J. Bungartz, F. Lindner, B. Gatzhammer, M. Mehl, K. Scheufele, A. Shukaev, and B. Uekermann, preCICE---a fully parallel library for multi-physics surface coupling, Comput. & Fluids, 141 (2016), pp. 250--258.
H.-j. Bungartz, F. Lindner, M. Mehl, K. Scheufele, A. Shukaev, and B. Uekermann, Partitioned Fluid-Structure-Acoustics Interaction on Distributed Data: Coupling via preCICE, Lect. Notes Comput. Sci. Eng., Springer, Berlin, Heidelberg, New York, submitted.
J. W. Daniel, W. B. Gragg, L. Kaufman, and G. Stewart, Reorthogonalization and stable algorithms for updating the Gram-Schmidt $QR$ factorization, Math. Comp., 30 (1976), pp. 772--795.
M. Darwish, I. Sraj, and F. Moukalled, A coupled finite volume solver for the solution of incompressible flows on unstructured grids, J. Comput. Phys., 228 (2009), pp. 180--201, https://doi.org/10.1016/j.jcp.2008.08.027.
A. de Boer, A. H. vanZuijlen, and H. Bijl, Radial Basis Functions for Interface Interpolation and Mesh deformation, in Advanced Computational Methods in Science and Engineering, Lect. Notes Comput. Sci. Eng. 71, Springer, Berlin, Heidelberg, 2010, pp. 143--178.
J. Degroote, K. J. Bathe, and J. Vierendeels, Performance of a new partitioned procedure versus a monolithic procedure in fluid-structure interaction, Comput. & Structures, 87 (2009), pp. 793--801, https://doi.org/10.1016/j.compstruc.2008.11.013.
J. Degroote, P. Bruggeman, R. Haelterman, and J. Vierendeels, Stability of a coupling technique for partitioned solvers in FSI applications, Comput. & Structures, 86 (2008), pp. 2224--2234, https://doi.org/10.1016/j.compstruc.2008.05.005.
H.-R. Fang and Y. Saad, Two classes of multisecant methods for nonlinear acceleration, Numer. Linear Algebra Appl., 16 (2008), pp. 197--221, https://doi.org/10.1002/nla.617.
M. A. Fernández and M. Moubachir, A Newton method using exact Jacobians for solving fluid-structure coupling, Comput. & Structures, 83 (2005), pp. 127--142, https://doi.org/10.1016/j.compstruc.2004.04.021.
B. Gatzhammer, Efficient and Flexible Partitioned Simulation of Fluid-Structure Interactions, Ph.D. thesis, Technische Universität München, 2014.
R. Haelterman, Analytical Study of the Least Squares Quasi-Newton Method for Interaction Problems, Ph.D. thesis, Ghent University, 2009, http://lib.ugent.be/fulltxt/RUG01/001/333/190/RUG01-001333190_2010_0001_AC.pdf.
R. Haelterman, A. Bogaers, B. Uekermann, K. Scheufele, and M. Mehl, Improving the performance of the partitioned QN-ILS procedure for fluid-structure interaction problems: Filtering, Comput. & Structures, 171 (2016), pp. 9--17.
R. Haelterman, J. Degroote, D. Van Heule, and J. Vierendeels, The quasi-Newton least squares method: A new and fast secant method analyzed for linear systems, SIAM J. Numer. Anal., 47 (2009), pp. 2347--2368, https://doi.org/10.1137/070710469.
F. Lindner, M. Mehl, K. Scheufele, and B. Uekermann, A comparison of various quasi-Newton schemes for partitioned fluid-structure interaction, in Proceedings of ECCOMAS Coupled Problems, Venice, 2015, pp. 1--12.
J. Loffeld and C. S. Woodward, Considerations on the Implementation and Use of Anderson Acceleration on Distributed Memory and GPU-Based Parallel Computers, Technical Report LLNLPROC-675918, Lawrence Livermore National Laboratory, 2015.
P. Lott, H. Walker, C. Woodward, and U. Meier-Yang, An accelerated Picard method for nonlinear systems related to variably saturated flow, Adv. Water Resources, 38 (2012), pp. 92--101, https://doi.org/10.1016/j.advwatres.2011.12.013.
M. Mehl, B. Uekermann, H. Bijl, D. Blom, B. Gatzhammer, and A. van Zuijlen, Parallel coupling numerics for partitioned fluid--structure interaction simulations, Comput. Math. Appl., 71 (2016), pp. 869--891, https://doi.org/10.1016/j.camwa.2015.12.025.
C. Michler, An interface Newton--Krylov solver for fluid--structure interaction, Internat. J. Numer. Methods Fluids, 47 (2004), pp. 1189--1195, https://doi.org/10.1002/fld.850.
C. Michler, H. Van Brummelen, and R. De Borst, An investigation of interface-GMRES(R) for fluid-structure interaction problems with flutter and divergence, Comput. Mech., 47 (2011), pp. 17--29, https://doi.org/10.1007/s00466-010-0519-8.
K. Scheufele, Robust Quasi-Newton Methods for Partitioned Fluid-Structure Simulations, Master thesis, University of Stuttgart, 2015.
A. Toth and C. T. Kelley, Convergence analysis for Anderson acceleration, SIAM J. Numer. Anal., 53 (2015), pp. 805--819, https://doi.org/10.1137/130919398.
S. Turek and J. Hron, Proposal for numerical benchmarking of fluid-structure interaction between an elastic object and laminar incompressible flow, in Fluid-Structure Interaction, Lect. Notes Comput. Sci. Eng. 53, Springer, Berlin, 2006, pp. 371--385.
B. Uekermann, Partitioned Fluid-Structure Interaction/Multi-Physics on Massively Parallel Systems, Ph.D. thesis, Technische Universität München, 2016.
B. Uekermann, H.-J. Bungartz, B. Gatzhammer, and M. Mehl, A parallel, black-box coupling for fluid-structure interaction, in Computational Methods for Coupled Problems in Science and Engineering, COUPLED PROBLEMS 2013, S. Idelsohn, M. Papadrakakis, and B. Schrefler, eds., Stanta Eulalia, Ibiza, Spain, 2013.
J. Vierendeels, J. Degroote, S. Annerel, and R. Haelterman, Stability issues in partitioned FSI calculations, in Fluid Structure Interaction II, H.-J. Bungartz, M. Mehl, and M. Schäfer, eds., Lect. Notes Comput. Sci. Eng., Springer, Heidelberg, Berlin, 2010, pp. 83--102, http://link.springer.com/chapter/10.1007/978-3-642-14206-2_4.
H. F. Walker and P. Ni, Anderson acceleration for fixed-point iterations, SIAM J. Numer. Anal., 49 (2011), pp. 1715--1735, https://doi.org/10.1137/10078356X.

Information & Authors


Published In

cover image SIAM Journal on Scientific Computing
SIAM Journal on Scientific Computing
Pages: S404 - S433
ISSN (online): 1095-7197


Submitted: 30 June 2016
Accepted: 22 February 2017
Published online: 26 October 2017


  1. partitioned multiphysics
  2. nonlinear fixed-point solver
  3. quasi-Newton
  4. fluid-structure interactions

MSC codes

  1. 65B99
  2. 65H10
  3. 65M99
  4. 68W10
  5. 74B99
  6. 76D05



Metrics & Citations



If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.

Cited By

View Options

View options


View PDF







Copy the content Link

Share with email

Email a colleague

Share on social media

The SIAM Publications Library now uses SIAM Single Sign-On for individuals. If you do not have existing SIAM credentials, create your SIAM account https://my.siam.org.