Space-time formulation, discretization, and computational performance studies for phase-field fracture optimal control problems

TL;DR

This paper develops space-time discretization schemes for phase-field fracture optimal control problems, focusing on minimizing a tracking cost to control crack patterns, and introduces a Newton algorithm with numerical validation.

Contribution

It introduces a novel space-time discretization and a Newton-based solution approach for phase-field fracture optimal control problems.

Findings

Discontinuous Galerkin time discretization effectively handles regularization and crack irreversibility.

The reduced Newton algorithm efficiently computes optimal controls.

Numerical experiments demonstrate the method's accuracy and computational performance.

Abstract

The purpose of this work is the development of space-time discretization schemes for phase-field optimal control problems. Specifically in the optimal control minimization problem, a tracking-type cost functional is minimized to steer the crack via the phase-field variable into a desired pattern. To achieve such optimal solutions, Neumann type boundary conditions need to be determined. First, a time discretization of the forward problem is derived using a discontinuous Galerkin formulation. Here, a challenge is to include regularization terms and the crack irreversibility constraint. The optimal control setting is formulated by means of the Lagrangian approach from which the primal part, adjoint, tangent and adjoint Hessian are derived. Herein the overall Newton algorithm is based on a reduced approach by eliminating the state constraint, namely the displacement and phase-field unknowns, but keeping the control variable as the only unknown. From the low-order discontinuous Galerkin discretization, adjoint time-stepping schemes are finally obtained. Both our formulation and algorithmic developments are substantiated and illustrated with six numerical experiments.