Mini-Workshop: Mathematical Analysis for Peridynamics

A mathematical analysis for peridynamics, a nonlocal elastic theory, is the subject of the mini-workshop. Peridynamics is a novel multiscale mechanical model where the canonical divergence of the stress tensor is replaced by an integral operator that sums forces at a finite distance. As such, the underlying regularity assumptions are more general, for instance, allowing discontinuous and non-differentiable displacement fields. Although the theoretical mechanical formulation of peridynamics is well understood, the mathematical and numerical analyses are in their early stages. The mini-workshop proved to be a catalyst for the emerging mathematical analyses among an international group of mathematicians

Nonlocal Neumann volume-constrained problems and their application to local-nonlocal coupling

As alternatives to partial differential equations (PDEs), nonlocal continuum models given in integral forms avoid the explicit use of conventional spatial derivatives and allow solutions to exhibit desired singular behavior. As an application, peridynamic models are reformulations of classical continuum mechanics that allow a natural treatment of discontinuities by replacing spatial derivatives of stress tensor with integrals of force density functions. The thesis is concerned about the mathematical perspective of nonlocal modeling and local-nonlocal coupling for fracture mechanics both theoretically and numerically. To this end, the thesis studies nonlocal diffusion models associated with ``Neumann-type'' constraints (or ``traction conditions'' in mechanics), a nonlinear peridynamic model for fracture mechanics with bond-breaking rules, and a multi-scale model with local-nonlocal coupling. In the computational studies, it is of practical interest to develop robust numerical schemes not only for the numerical solution of nonlocal models, but also for the evaluation of suitably defined derivatives of solutions. This leads to a posteriori nonlocal stress analysis for structure mechanical models

A scalable domain decomposition method for FEM discretizations of nonlocal equations of integrable and fractional type

Nonlocal models allow for the description of phenomena which cannot be captured by classical partial differential equations. The availability of efficient solvers is one of the main concerns for the use of nonlocal models in real world engineering applications. We present a domain decomposition solver that is inspired by substructuring methods for classical local equations. In numerical experiments involving finite element discretizations of scalar and vectorial nonlocal equations of integrable and fractional type, we observe improvements in solution time of up to 14.6x compared to commonly used solver strategies

An efficient implementation of an implicit FEM scheme for fractional-in-space reaction-diffusion equations

Fractional differential equations are becoming increasingly used as a modelling tool for processes with anomalous diffusion or spatial heterogeneity. However, the presence of a fractional differential operator causes memory (time fractional) or nonlocality (space fractional) issues, which impose a number of computational constraints. In this paper we develop efficient, scalable techniques for solving fractional-in-space reaction diffusion equations using the finite element method on both structured and unstructured grids, and robust techniques for computing the fractional power of a matrix times a vector. Our approach is show-cased by solving the fractional Fisher and fractional Allen-Cahn reaction-diffusion equations in two and three spatial dimensions, and analysing the speed of the travelling wave and size of the interface in terms of the fractional power of the underlying Laplacian operator