In this work, we address the implementation and performance of inexact Newton–Krylov and quasi-Newton algorithms, more specifically the BFGS method, for the solution of the nonlinear elasticity equations, and compare them to a standard Newton–Krylov method. This is done through a systematic analysis of the performance of the solvers with respect to the problem size, the magnitude of the data and the number of processors in both almost incompressible and incompressible mechanics. We consider three test cases: Cook's membrane (static, almost incompressible), a twist test (static, incompressible) and a cardiac model (complex material, time dependent, almost incompressible). Our results suggest that quasi-Newton methods should be preferred for compressible mechanics, whereas inexact Newton–Krylov methods should be preferred for incompressible problems. We show that these claims are also backed up by the convergence analysis of the methods. In any case, all methods present adequate performance, and provide a significant speed-up over the standard Newton–Krylov method, with a CPU time reduction exceeding 50% in the best cases.

Parallel inexact Newton–Krylov and quasi-Newton solvers for nonlinear elasticity / N.A. Barnafi, L.F. Pavarino, S. Scacchi. - In: COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING. - ISSN 0045-7825. - 400:(2022 Oct), pp. 115557.1-115557.22. [10.1016/j.cma.2022.115557]

Parallel inexact Newton–Krylov and quasi-Newton solvers for nonlinear elasticity

N.A. Barnafi
Primo
;
S. Scacchi
Ultimo
2022

Abstract

In this work, we address the implementation and performance of inexact Newton–Krylov and quasi-Newton algorithms, more specifically the BFGS method, for the solution of the nonlinear elasticity equations, and compare them to a standard Newton–Krylov method. This is done through a systematic analysis of the performance of the solvers with respect to the problem size, the magnitude of the data and the number of processors in both almost incompressible and incompressible mechanics. We consider three test cases: Cook's membrane (static, almost incompressible), a twist test (static, incompressible) and a cardiac model (complex material, time dependent, almost incompressible). Our results suggest that quasi-Newton methods should be preferred for compressible mechanics, whereas inexact Newton–Krylov methods should be preferred for incompressible problems. We show that these claims are also backed up by the convergence analysis of the methods. In any case, all methods present adequate performance, and provide a significant speed-up over the standard Newton–Krylov method, with a CPU time reduction exceeding 50% in the best cases.
BFGS; Inexact-Newton; Newton–Krylov; nonlinear elasticity; scalable solvers
Settore MAT/08 - Analisi Numerica
   Modeling the heart across the scales: from cardiac cells to the whole organ
   MINISTERO DELL'ISTRUZIONE E DEL MERITO
   2017AXL54F_003
ott-2022
Article (author)
File in questo prodotto:
File Dimensione Formato  
Superlinear_solvers_in_nonlinear_elasticity.pdf

embargo fino al 01/10/2024

Tipologia: Post-print, accepted manuscript ecc. (versione accettata dall'editore)
Dimensione 725.98 kB
Formato Adobe PDF
725.98 kB Adobe PDF   Visualizza/Apri   Richiedi una copia
1-s2.0-S0045782522005424-main.pdf

accesso riservato

Tipologia: Publisher's version/PDF
Dimensione 2.33 MB
Formato Adobe PDF
2.33 MB Adobe PDF   Visualizza/Apri   Richiedi una copia
2203.05610.pdf

accesso aperto

Tipologia: Pre-print (manoscritto inviato all'editore)
Dimensione 771.72 kB
Formato Adobe PDF
771.72 kB Adobe PDF Visualizza/Apri
Pubblicazioni consigliate

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/2434/944031
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 8
  • ???jsp.display-item.citation.isi??? 5
social impact