A Regularized Ensemble Kalman Filter for Stochastic Phase Field Models of Brittle Fracture
This work addresses uncertainty in fracture simulations for structural monitoring applications, but it is incremental as it builds on existing ensemble Kalman filter methods with a regularization step.
The paper tackles the problem of uncertain material parameters in phase-field models of brittle fracture by developing a Bayesian inference procedure that updates the model state (displacements and phase-field) using sensor data, with results showing the updated state matches ground truth reasonably well in 1D and 2D examples.
The phase-field approach to brittle fracture provides a continuum framework for modeling crack initiation and propagation without explicit representation of discrete crack surfaces, provided the spatial discretization is fine enough to resolve the regularization length scale. However, uncertain local material parameters due to material defects can strongly influence simulation results, such as crack paths and remaining structural strength. At the same time, the ability to continuously monitor structures using sensors allows complementing modeling predictions with, e.g., displacement measurements. In this contribution, we connect these two complementary sources of information and present a Bayesian inference procedure that allows updating the current model state with incoming sensor data. We construct a Bayesian prior for the model state (both displacements and phase-field) and employ an ensemble Kalman filter (EnKF) to perform the update. In the EnKF, the update is computed by performing a Kalman shift on each ensemble member. Since the standard EnKF may produce assimilated states that violate common modeling assumptions, we present a phase field-based regularization technique as a proximal step correction toward model-consistent updates. 1D and 2D numerical examples demonstrate the performance and accuracy of the proposed method and show that the updated state matches the ground truth reasonably well. Unlike traditional Bayesian inversion techniques, which have already been applied to brittle fracture, we infer not the model parameters but the model state, i.e., the displacement field and the phase-field. Although only displacements are observed, the strong correlation between both fields also allows inference of the posterior phase-field.