We study the convergence of an alternate minimization scheme for a Ginzburg–Landau phase-field model of fracture. This algorithm is characterized by the lack of irreversibility constraints in the minimization of the phase-field variable; the advantage of this choice, from a computational stand point, is in the efficiency of the numerical implementation. Irreversibility is then recovered a posteriori by a simple pointwise truncation. We exploit a time discretization procedure, with either a one-step or a multi (or infinite)-step alternate minimization algorithm. We prove that the time-discrete solutions converge to a unilateral $L^2$ -gradient flow with respect to the phase-field variable, satisfying equilibrium of forces and energy identity. Convergence is proved in the continuous (Sobolev space) setting and in a discrete (finite element) setting, with any stopping criterion for the alternate minimization scheme. Numerical results show that the multi-step scheme is both more accurate and faster. It provides indeed good simulations for a large range of time increments, while the one-step scheme gives comparable results only for very small time increments.

Convergence of discrete and continuous unilateral flows for Ambrosio-Tortorelli energies and application to mechanics

Negri Matteo
2019-01-01

Abstract

We study the convergence of an alternate minimization scheme for a Ginzburg–Landau phase-field model of fracture. This algorithm is characterized by the lack of irreversibility constraints in the minimization of the phase-field variable; the advantage of this choice, from a computational stand point, is in the efficiency of the numerical implementation. Irreversibility is then recovered a posteriori by a simple pointwise truncation. We exploit a time discretization procedure, with either a one-step or a multi (or infinite)-step alternate minimization algorithm. We prove that the time-discrete solutions converge to a unilateral $L^2$ -gradient flow with respect to the phase-field variable, satisfying equilibrium of forces and energy identity. Convergence is proved in the continuous (Sobolev space) setting and in a discrete (finite element) setting, with any stopping criterion for the alternate minimization scheme. Numerical results show that the multi-step scheme is both more accurate and faster. It provides indeed good simulations for a large range of time increments, while the one-step scheme gives comparable results only for very small time increments.
File in questo prodotto:
Non ci sono file associati a questo prodotto.

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/11571/1242526
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 11
  • ???jsp.display-item.citation.isi??? 9
social impact