In this work we study the numerical solution of nonlinear systems arising from stabilized FEM discretizations of Navier–Stokes equations. This is a very challenging problem and often inexact Newton solvers present severe difficulties to converge. Then, they must be wrapped into a globalization strategy. We consider the classical backtracking procedure, a subspace trust-region strategy and an hybrid approach. This latter strategy is proposed with the aim of improve the robustness of backtracking and it is obtained combining the backtracking procedure and the elliptical subspace trust-region strategy. Under standard assumptions, we prove global and fast convergence of the inexact Newton methods embedded in this new strategy as well as in the subspace trust-region procedure. Computational results on classical CFD benchmarks are performed. Comparisons among the classical backtracking strategy, the elliptical subspace trust-region approach and the hybrid procedure are presented. Our numerical experiments show the effectiveness of the proposed hybrid technique.
Globalization strategies for Newton-Krylov methods for stabilized FEM discretization of Navier-Stokes equations / S. BELLAVIA; S. BERRONE. - In: JOURNAL OF COMPUTATIONAL PHYSICS. - ISSN 0021-9991. - STAMPA. - 226:(2007), pp. 2317-2340. [10.1016/j.jcp.2007.07.021]
Globalization strategies for Newton-Krylov methods for stabilized FEM discretization of Navier-Stokes equations
BELLAVIA, STEFANIA;
2007
Abstract
In this work we study the numerical solution of nonlinear systems arising from stabilized FEM discretizations of Navier–Stokes equations. This is a very challenging problem and often inexact Newton solvers present severe difficulties to converge. Then, they must be wrapped into a globalization strategy. We consider the classical backtracking procedure, a subspace trust-region strategy and an hybrid approach. This latter strategy is proposed with the aim of improve the robustness of backtracking and it is obtained combining the backtracking procedure and the elliptical subspace trust-region strategy. Under standard assumptions, we prove global and fast convergence of the inexact Newton methods embedded in this new strategy as well as in the subspace trust-region procedure. Computational results on classical CFD benchmarks are performed. Comparisons among the classical backtracking strategy, the elliptical subspace trust-region approach and the hybrid procedure are presented. Our numerical experiments show the effectiveness of the proposed hybrid technique.File | Dimensione | Formato | |
---|---|---|---|
YJCPH1763.pdf
Accesso chiuso
Tipologia:
Versione finale referata (Postprint, Accepted manuscript)
Licenza:
Tutti i diritti riservati
Dimensione
528.5 kB
Formato
Adobe PDF
|
528.5 kB | Adobe PDF | Richiedi una copia |
I documenti in FLORE sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.