Journal article

Large scale simulation of pressure induced phase-field fracture propagation using Utopia

  • Zulian, Patrick Euler Institute (EUL), Università della Svizzera italiana, Svizzera
  • Kopaničáková, Alena Euler Institute (EUL), Università della Svizzera italiana, Svizzera
  • Nestola, Maria Giuseppina Chiara Euler Institute (EUL), Università della Svizzera italiana, Svizzera - Institute of Geochemistry and Petrology, ETH Zurich, Zurich, Switzerland
  • Fink, Andreas Swiss National Supercomputing Centre (CSCS), ETH Zurich, Zurich, Switzerland
  • Fadel, Nur Aiman Swiss National Supercomputing Centre (CSCS), ETH Zurich, Zurich, Switzerland
  • VandeVondele, Joost Swiss National Supercomputing Centre (CSCS), ETH Zurich, Zurich, Switzerland
  • Krause, Rolf Euler Institute (EUL), Università della Svizzera italiana, Svizzera
Show more…
    29.06.2021
Published in:
  • CCF transactions on high performance computing. - Springer. - 2021, no. 3, p. 407–426
English Non-linear phase field models are increasingly used for the simulation of fracture propagation problems. The numerical simulation of fracture networks of realistic size requires the efficient parallel solution of large coupled non-linear systems. Although in principle efficient iterative multi-level methods for these types of problems are available, they are not widely used in practice due to the complexity of their parallel implementation. Here, we present Utopia, which is an open-source C++ library for parallel non-linear multilevel solution strategies. Utopia provides the advantages of high-level programming interfaces while at the same time a framework to access low-level data-structures without breaking code encapsulation. Complex numerical procedures can be expressed with few lines of code, and evaluated by different implementations, libraries, or computing hardware. In this paper, we investigate the parallel performance of our implementation of the recursive multilevel trust-region (RMTR) method based on the Utopia library. RMTR is a globally convergent multilevel solution strategy designed to solve non-convex constrained minimization problems. In particular, we solve pressure-induced phasefield fracture propagation in large and complex fracture networks. Solving such problems is deemed challenging even for a few fractures, however, here we are considering networks of realistic size with up to 1000 fractures.
Collections
Language
  • English
Classification
Computer science
License
License undefined
Identifiers
Persistent URL
https://susi.usi.ch/usi/documents/319389
Statistics

Document views: 28 File downloads:
  • Zulian_ccfthpc_2021.pdf: 17