Numerical Analysis of Nonlocal Fracture Models in Holder Space

Apr 10, 2018·
Prashant K. Jha
Prashant K. Jha
,
Robert Lipton
· 0 min read
Abstract
In this work, we calculate the convergence rate of the finite difference approximation for a class of nonlocal fracture models. We consider two point force interactions characterized by a double well potential. We show the existence of a evolving displacement field in Holder space with Holder exponent $\gamma \in (0, 1]$. The rate of convergence of the finite difference approximation depends on the factor $C_s h \gamma/\epsilon^2$ where $\epsilon$ gives the length scale of nonlocal interaction, $h$ is the discretization length, and $C_s$ is the maximum of Holder norm of the solution and its second derivatives during the evolution. It is shown that the rate of convergence holds for both the forward Euler scheme as well as general single step implicit schemes. A stability result is established for the semidiscrete approximation. The Holder continuous evolutions are seen to converge to a brittle fracture evolution in the limit of vanishing nonlocality.
Type
Publication
SIAM Journal on Numerical Analysis
Prashant K. Jha
Authors
Assistant Professor of Mechanical Engineering
Our group uses mechanics, applied mathematics, and computational science to understand and represent the complex behavior of materials, e.g., functional soft materials and granular materials.