In this contribution we address the issue of efficient finite element treatment for phase-field modeling of brittle fracture. We start by providing an overview of the existing quasi-static and dynamic phase-field fracture formulations from the physics and the mechanics communities. Within the formulations stemming from Griffith's theory, we focus on quasi-static models featuring a tension-compression split, which prevent cracking in compression and interpenetration of the crack faces upon closure, and on the staggered algorithmic implementation due to its proved robustness. In this paper, we establish an appropriate stopping criterion for the staggered scheme. Moreover, we propose and test the so-called hybrid formulation, which leads within a staggered implementation to an incrementally linear problem. This enables a significant reduction of computational cost--about one order of magnitude--with respect to the available (non-linear) models. The conceptual and structural similarities of the hybrid formulation to gradient-enhanced continuum damage mechanics are outlined as well. Several benchmark problems are solved, including one with own experimental verification.

A review on phase-field models of brittle fracture and a new fast hybrid formulation

The phase-field method has become an important and extremely versatile technique for simulating microstructure evolution at the mesoscale. Thanks to the diffuse-interface approach, it allows us to study the evolution of arbitrary complex grain morphologies without any presumption on their shape or mutual distribution. It is also straightforward to account for different thermodynamic driving forces for microstructure evolution, such as bulk and interfacial energy, elastic energy and electric or magnetic energy, and the effect of different transport processes, such as mass diffusion, heat conduction and convection. The purpose of the paper is to give an introduction to the phase-field modeling technique. The concept of diffuse interfaces, the phase-field variables, the thermodynamic driving force for microstructure evolution and the kinetic phase-field equations are introduced. Furthermore, common techniques for parameter determination and numerical solution of the equations are discussed. To show the variety in phase-field models, different model formulations are exploited, depending on which is most common or most illustrative.

/pdf/an-introduction-to-phase-field-modeling-of-microstructure-2cm9coegv9.pdf

An introduction to phase-field modeling of microstructure evolution

Common to liquids, solids and substances in between the former two is that if a stress is applied to them, they will strain. Stress may be visualized by placing a small amount of fluid between two parallel plates. When one plate slides over the other, forces act on the fluid dependent upon the rate of the plate movement. This causes a shear stress on the liquid. Recall laminar flow of fluids through a tubular vessel. Strain is the response to the stress. If solids are elastic, they deform and return to their original shape. Since fluids are not elastic and, hence, viscous, their deformation is irreversible.

Introduction to rheology

Models of cardiac tissue electrophysiology are an important component of the Cardiac Physiome Project, which is an international effort to build biophysically based multi-scale mathematical models of the heart. Models of tissue electrophysiology can provide a bridge between electrophysiological cell models at smaller scales, and tissue mechanics, metabolism and blood flow at larger scales. This paper is a critical review of cardiac tissue electrophysiology models, focussing on the micro-structure of cardiac tissue, generic behaviours of action potential propagation, different models of cardiac tissue electrophysiology, the choice of parameter values and tissue geometry, emergent properties in tissue models, numerical techniques and computational issues. We propose a tentative list of information that could be included in published descriptions of tissue electrophysiology models, and used to support interpretation and evaluation of simulation results. We conclude with a discussion of challenges and open questions.

/pdf/models-of-cardiac-tissue-electrophysiology-progress-3n07jn3xug.pdf

Models of cardiac tissue electrophysiology: progress, challenges and open questions.

Dynamic fracture mechanics

We present a phase-field model of the propagation of fracture under plane strain. This model, based on simple physical considerations, is able to accurately reproduce the different behavior of cracks (the principle of local symmetry, the Griffith and Irwin criteria, and mode-I branching). In addition, we test our model against recent experimental findings showing the presence of oscillating cracks under biaxial load. Our model again reproduces well observed supercritical Hopf bifurcation and is therefore the first simulation which does so.

Dynamic instabilities of fracture under biaxial strain using a phase field model.

Scroll waves are three-dimensional analogs of spiral waves. The linear stability spectrum of untwisted and twisted scroll waves is computed for a two-variable reaction-diffusion model of an excitable medium. Different bands of modes are seen to be unstable in different regions of parameter space. The corresponding bifurcations and bifurcated states are characterized by performing direct numerical simulations. In addition, computations of the adjoint linear stability operator eigenmodes are also performed and serve to obtain a number of matrix elements characterizing the long-wavelength deformations of scroll waves.

Scroll waves in isotropic excitable media: linear instabilities, bifurcations, and restabilized states.

In this study, the phase field model of crack propagation is used to study the dynamic branching instability in the case of inplane loading in two dimensions. Simulation results are in good agreement with theoretical predictions and experimental findings. Namely, the critical speed at which the instability starts is about 0.48cs. They also show that a full 3D approach is needed to fully understand the branching instability. The finite interface effects are found to be neglectable in the large system size limit even though they are stronger than the one expected from a simple one-dimensional calculation.

https://iopscience.iop.org/article/10.1209/0295-5075/83/16004/pdf

Study of the branching instability using a phase field model of inplane crack propagation

We address analytically and numerically the problem of crack path prediction in the model system of a crack propagating under thermal loading. We show that one can explain the instability from a straight to a wavy crack propagation by using only the principle of local symmetry and the Griffith criterion. We then argue that the calculations of the stress intensity factors can be combined with the standard crack propagation criteria to obtain the evolution equation for the crack tip within any loading configuration. The theoretical results of the thermal crack problem agree with the numerical simulations we performed using a phase field model. Moreover, it turns out that the phase-field model allows to clarify the nature of the transition between straight and oscillatory cracks which is shown to be supercritical.

Thermal fracture as a framework for quasi-static crack propagation

/pdf/thermal-fracture-as-a-framework-for-quasi-static-crack-4prukg8dux.pdf

Hervé Henry

Papers

Dynamic instabilities of fracture under biaxial strain using a phase field model.

Scroll waves in isotropic excitable media: linear instabilities, bifurcations, and restabilized states.

Study of the branching instability using a phase field model of inplane crack propagation

Thermal fracture as a framework for quasi-static crack propagation

Thermal fracture as a framework for quasi-static crack propagation