Fine numerical analysis of the crack-tip position for a Mumford-Shah minimizer

DOI10.4171/IFB/357MaRDI QIDQ2630802zbMATH OpenOpenAlexWikidataFDO

Authors H. Mikayelyan, Zhilin Li

Publication date 22 July 2016

Published in Interfaces and Free Boundaries (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1511.07733

zbMATH Keywords

augmented immersed interface methodcrack-tipfast/Poisson solverfree boundaryfree discontinuityirregular domainMumford-Shah energy

Mathematics Subject Classification ID

Free boundary problems for PDEs (35R35) Optimality conditions for problems involving partial differential equations (49K20) Discrete approximations in optimal control (49M25) Finite difference methods for initial value and initial-boundary value problems involving PDEs (65M06) Fictitious domain methods for initial value and initial-boundary value problems involving PDEs (65M85) Theories of fracture and damage (74A45) Finite difference methods applied to problems in solid mechanics (74S20)

Abstract: A new algorithm to determine the position of the crack (discontinuity set) of certain minimizers of Mumford-Shah functional in situations when a crack-tip occurs is introduced. The conformal mapping $w=sqrtz$ in the complex plane is used to transform the free discontinuity problem to a new type of free boundary problem, where the symmetry of the free boundary is an additional constraint of a non-local nature. Instead of traditional Jacobi or Newton iterative methods, we propose a simple iteration method which does not need the Jacobian but is way fast than the Jacobi iteration. In each iteration, a Laplace equation needs to be solved on an irregular domain with a Dirichlet boundary condition on the fixed part of the boundary; and a Neumann type boundary condition along the free boundary. The augmented immersed interface method is employed to solve the potential problem. The numerical results agree with the analytic analysis and provide insight into some open questions in free discontinuity problems.

Review text

Summary: A new algorithm to determine the position of the crack (discontinuity set) of certain minimizers of Mumford-Shah functional in situations when a crack-tip occurs is introduced. The conformal mapping \(\tilde{z}=\sqrt{z}\) in the complex plane is used to transform the free discontinuity problem to a new type of free boundary problem, where the symmetry of the free boundary is an additional constraint of a non-local nature. Instead of traditional Jacobi or Newton iterative methods, we propose a simple iteration method which does not need the Jacobian but is way fast than the Jacobi iteration. In each iteration, a Laplace equation needs to be solved on an irregular domain with a Dirichlet boundary condition on the fixed part of the boundary; and a Neumann type boundary condition along the free boundary. The augmented immersed interface method is employed to solve the potential problem. The numerical results agree with the analytic analysis and provide insight into some open questions in free discontinuity problems.

Recommendations

• Topological indices of intersection ideal graph of \(\mathbb{Z}_n\)
• Cracktip is a global Mumford-Shah minimizer
• Implementation of an adaptive finite-element approximation of the Mumford-Shah functional
• Stationarity of the crack-front for the Mumford-Shah problem in 3D
• Some recent results on regularity of the crack-tip/crack-front of Mumford-Shah minimizers

Cites work

• A model for the quasi-static growth of brittle fractures: Existence and approximation results
• Cracktip is a global Mumford-Shah minimizer
• Higher integrability for minimizers of the Mumford-Shah functional
• scientific article; zbMATH DE number 1448982 (Why is no real title available?)
• Revisiting brittle fracture as an energy minimization problem
• Stable regular critical points of the Mumford-Shah functional are local minimizers

Cited in

(2)

• Stationarity of the crack-front for the Mumford-Shah problem in 3D
• Effective matrix-free preconditioning for the augmented immersed interface method

Describes a project that uses

Uses Software

• IIMPACK

This page was built for publication: Fine numerical analysis of the crack-tip position for a Mumford-Shah minimizer

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2630802)