Stochastic phase-field modeling of brittle fracture: computing multiple crack patterns and their probabilities

by   Tymofiy Gerasimov, et al.

In variational phase-field modeling of brittle fracture, the functional to be minimized is not convex, so that the necessary stationarity conditions of the functional may admit multiple solutions. The solution obtained in an actual computation is typically one out of several local minimizers. Evidence of multiple solutions induced by small perturbations of numerical or physical parameters was occasionally recorded but not explicitly investigated in the literature. In this work, we focus on this issue and advocate a paradigm shift, away from the search for one particular solution towards the simultaneous description of all possible solutions (local minimizers), along with the probabilities of their occurrence. Inspired by recent approaches advocating measure-valued solutions (Young measures as well as their generalization to statistical solutions) and their numerical approximations in fluid mechanics, we propose the stochastic relaxation of the variational brittle fracture problem through random perturbations of the functional. We introduce the concept of stochastic solution, with the main advantage that point-to-point correlations of the crack phase fields in the underlying domain can be captured. These stochastic solutions are represented by random fields or random variables with values in the classical deterministic solution spaces. In the numerical experiments, we use a simple Monte Carlo approach to compute approximations to such stochastic solutions. The final result of the computation is not a single crack pattern, but rather several possible crack patterns and their probabilities. The stochastic solution framework using evolving random fields allows additionally the interesting possibility of conditioning the probabilities of further crack paths on intermediate crack patterns.



There are no comments yet.


page 4

page 8

page 13

page 14

page 16

page 32

page 33

page 35


On the Numerical Approximation of the Karhunen-Loève Expansion for Lognormal Random Fields

The Karhunen-Loève (KL) expansion is a popular method for approximating ...

Efficient Parallel Estimation for Markov Random Fields

We present a new, deterministic, distributed MAP estimation algorithm fo...

Analysis of Spherical Monofractal and Multifractal Random Fields with Cosmological Applications

The Rényi function plays an important role in the analysis of multifract...

Numerical approximation of statistical solutions of the incompressible Navier-Stokes Equations

Statistical solutions, which are time-parameterized probability measures...

Variational phase recovering without phase unwrapping in phase-shifting interferometry

We present a variational method for recovering the phase term from the i...

Enhanced Multi-Index Monte Carlo by means of Multiple Semi-Coarsened Multigrid for Anisotropic Diffusion Problems

In many models used in engineering and science, material properties are ...

Solution of the 3D density-driven groundwater flow problem with uncertain porosity and permeability

As groundwater is an essential nutrition and irrigation resource, its po...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.