This merge request provides a restart simulation capability for the phase field process, allowing the simulation to be resumed from the end of a previous run by reading the stress at Gauss points.
A ctest is provided. Basically, the existing surfing benchmark is divided into two parts:
The first part, initial simulation, runs from time 0 to 0.5, during which the displacement, phase field at mesh nodes, and stress at Gauss points are stored in the mesh.
The second part, restart simulation, continues the simulation from time 0.5 to 1, using the stored data as initial conditions.
In the project file, a new computational mesh is utilized (vtu file of last time step of initial simulation), storing the stress, displacement, and phase field from the initial simulation. These values are then used as the initial conditions for the restart simulation.
The restart simulation runs from time 0.5 ( t_initial) to 1 (t_end). We expect the phase field and displacement at the end of this restart simulation to match those in the original surfing benchmark.
The following plot shows the initial stress at Gauss points at the start of the restart simulation, generated using the ipDataToPointCloud utility tool in OGS. The plot highlights the stress singularity at the fracture tip.
