Efficient finite strain elasticity solver for phase-field simulations

Oleg Shchyglo,Muhammad Adil Ali,Hesham Salama

DOI: https://doi.org/10.1038/s41524-024-01235-4

IF: 12.256

2024-03-15

npj Computational Materials

Abstract:We present an effective mechanical equilibrium solution algorithm suitable for finite strain consideration within the phase-field method. The proposed algorithm utilizes a Fourier space solution in its core. The performance of the proposed algorithm is demonstrated using the St. Venant–Kirchhoff hyperelastic model, but the algorithm is also applicable to other hyperelastic models. The use of the fast Fourier transformation routines and fast convergence within several iterations for most common simulation scenarios makes the proposed algorithm suitable for phase-field simulations of rapidly evolving microstructures. Additionally, the proposed algorithm allows using different strain measures depending on the requirements of the underlying problem. The algorithm is implemented in the OpenPhase phase-field simulation library. A set of example simulations ranging from simple geometries to complex microstructures is presented. The effect of different externally applied mechanical boundary conditions and internal forces is also demonstrated. The proposed algorithm can be considered a straightforward update to already existing small strain solvers based on Fourier space solutions.

materials science, multidisciplinary,chemistry, physical

What problem does this paper attempt to address?

The paper aims to address the effective solution methods for finite strain elasticity problems in phase-field simulations. Specifically, the paper proposes a core algorithm based on Fourier space solutions, which is suitable for considering finite strains in phase-field methods. This algorithm has the following features:

1. **Efficiency**: By utilizing fast Fourier transforms and the ability to converge quickly within a few iterations, this algorithm is suitable for phase-field simulations of rapidly evolving microstructures.
2. **Generality**: Although the paper uses the Saint Venant-Kirchhoff hyperelastic model to demonstrate the algorithm's performance, the algorithm is also applicable to other hyperelastic models.
3. **Flexibility**: It allows the use of different strain measures according to the requirements of the underlying problem.
4. **Implementation**: The algorithm has been implemented in the OpenPhase phase-field simulation library and provides a range of example simulations from simple geometries to complex microstructures.

Additionally, the paper demonstrates the effects of different external mechanical boundary conditions and internal forces, and points out that the algorithm can be seen as a direct update to existing small strain solvers based on Fourier space solutions.