An adapted deflated conjugate gradient solver for robust extended/generalised finite element solutions of large scale, 3D crack propagation problems

An adapted deflation preconditioner is employed to accelerate the solution of linear systems resulting from the discretization of fracture mechanics problems with well-conditioned extended/generalized finite elements. The deflation space typically used for linear elasticity problems is enriched with additional vectors, accounting for the enrichment functions used, thus effectively removing low frequency components of the error. To further improve performance, deflation is combined, in a multiplicative way, with a block-Jacobi preconditioner, which removes high frequency components of the error as well as linear dependencies introduced by enrichment. The resulting scheme is tested on a series of non-planar crack propagation problems and compared to alternative linear solvers in terms of performance.

[1]  I. Babuska,et al.  The design and analysis of the Generalized Finite Element Method , 2000 .

[2]  Stéphane Bordas,et al.  Error-controlled adaptive extended finite element method for 3D linear elastic crack propagation , 2017 .

[3]  V. E. Henson,et al.  BoomerAMG: a parallel algebraic multigrid solver and preconditioner , 2002 .

[4]  Emmanouil Kakouris,et al.  Discrete and Phase Field Methods for Linear Elastic Fracture Mechanics: A Comparative Study and State-of-the-Art Review , 2019, Applied Sciences.

[5]  Andrea Toselli,et al.  Domain decomposition methods : algorithms and theory , 2005 .

[6]  Hans Albert Richard,et al.  Fracture analyses and experimental results of crack growth under general mixed mode loading conditions , 2004 .

[7]  S. Bordas,et al.  A well‐conditioned and optimally convergent XFEM for 3D linear elastic fracture , 2016 .

[8]  Ted Belytschko,et al.  A finite element method for crack growth without remeshing , 1999 .

[9]  Eric de Sturler,et al.  Recycling Krylov Subspaces for Sequences of Linear Systems , 2006, SIAM J. Sci. Comput..

[10]  Eleni Chatzi,et al.  Stable 3D XFEM/vector level sets for non‐planar 3D crack propagation and comparison of enrichment schemes , 2018 .

[11]  T. Belytschko,et al.  A review of extended/generalized finite element methods for material modeling , 2009 .

[12]  Utkarsh Ayachit,et al.  The ParaView Guide: A Parallel Visualization Application , 2015 .

[13]  P Kerfriden,et al.  Bridging Proper Orthogonal Decomposition methods and augmented Newton-Krylov algorithms: an adaptive model order reduction for highly nonlinear mechanical problems. , 2011, Computer methods in applied mechanics and engineering.

[14]  Giulio Ventura,et al.  Stabilized X-FEM for Heaviside and Nonlinear Enrichments , 2016 .

[15]  Michel Salaün,et al.  High‐order extended finite element method for cracked domains , 2005 .

[16]  S. Bordas,et al.  A robust preconditioning technique for the extended finite element method , 2011 .

[17]  I. Babuska,et al.  The partition of unity finite element method: Basic theory and applications , 1996 .

[18]  Barry F. Smith,et al.  Domain Decomposition: Parallel Multilevel Methods for Elliptic Partial Differential Equations , 1996 .

[19]  I. Babuska,et al.  Stable Generalized Finite Element Method (SGFEM) , 2011, 1104.0960.

[20]  Ivo Babuška,et al.  Generalized finite element methods for three-dimensional structural mechanics problems , 2000 .

[21]  Stéphane Bordas,et al.  An extended finite element library , 2007 .

[22]  Robert Scheichl,et al.  High-performance dune modules for solving large-scale, strongly anisotropic elliptic problems with applications to aerospace composites , 2020, Comput. Phys. Commun..

[23]  C. Duarte,et al.  Analysis and applications of a generalized finite element method with global-local enrichment functions , 2008 .

[24]  M. Duflot.  A study of the representation of cracks with level sets , 2007 .

[25]  Olaf Schenk,et al.  Solving unsymmetric sparse systems of linear equations with PARDISO , 2002, Future Gener. Comput. Syst..

[26]  Cornelis Vuik,et al.  Theoretical and numerical comparison of various projection methods derived from deflation, domain decomposition and multigrid methods , 2007 .

[27]  N. Moës,et al.  Improved implementation and robustness study of the X‐FEM for stress analysis around cracks , 2005 .

[28]  Rainald Löhner,et al.  Deflated preconditioned conjugate gradient solvers for linear elasticity , 2011 .

[29]  David E. Keyes,et al.  A Quasi-algebraic Multigrid Approach to Fracture Problems Based on Extended Finite Elements , 2012, SIAM J. Sci. Comput..

[30]  E. Chatzi,et al.  A unified enrichment approach addressing blending and conditioning issues in enriched finite elements , 2019, Computer Methods in Applied Mechanics and Engineering.

[31]  Varun Gupta,et al.  Improved conditioning and accuracy of GFEM/XFEM for three-dimensional fracture mechanics , 2015 .

[32]  Christophe Geuzaine,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[33]  Marc Duflot,et al.  Derivative recovery and a posteriori error estimate for extended finite elements , 2007 .

[34]  Chupeng Ma,et al.  Novel design and analysis of generalized FE methods based on locally optimal spectral approximations , 2021, ArXiv.

[35]  N. Chevaugeon,et al.  Improved crack tip enrichment functions and integration for crack modeling using the extended finite element method , 2013 .

[36]  Eleni Chatzi,et al.  Multiple crack detection in 3D using a stable XFEM and global optimization , 2018, Computational Mechanics.

[37]  Xiang Ren,et al.  A novel hybrid approach for level set characterization and tracking of non-planar 3D cracks in the extended finite element method , 2016 .

[38]  R. Nicolaides.  Deflation of conjugate gradients with applications to boundary value problems , 1987 .

[39]  S. Bordas,et al.  A posteriori error estimation for extended finite elements by an extended global recovery , 2008 .

[40]  G. Ventura,et al.  Well Conditioned Extended Finite Elements and Vector Level Sets for Three-Dimensional Crack Propagation , 2017 .

[41]  Miss A.O. Penney.  (b) , 1974, The New Yale Book of Quotations.

[42]  Olaf Schenk,et al.  Two-level dynamic scheduling in PARDISO: Improved scalability on shared memory multiprocessing systems , 2002, Parallel Comput..

[43]  Thomas de Quincey.  [C] , 2000, The Works of Thomas De Quincey, Vol. 1: Writings, 1799–1820.

[44]  E. Chatzi,et al.  Improving the conditioning of XFEM/GFEM for fracture mechanics problems through enrichment quasi-orthogonalization , 2019, Computer Methods in Applied Mechanics and Engineering.

[45]  T. Belytschko,et al.  Vector level sets for description of propagating cracks in finite elements , 2003 .

[46]  Cornelis Vuik,et al.  On POD-based Deflation Vectors for DPCG applied to porous media problems , 2018, J. Comput. Appl. Math..

[47]  Haim Waisman,et al.  Updating strategy of a domain decomposition preconditioner for parallel solution of dynamic fracture problems , 2020, J. Comput. Phys..

[48]  A. G. Sanchez-Rivadeneira,et al.  A stable generalized/eXtended FEM with discontinuous interpolants for fracture mechanics , 2019, Computer Methods in Applied Mechanics and Engineering.

[49]  T. Belytschko,et al.  The extended/generalized finite element method: An overview of the method and its applications , 2010 .

[50]  Stéphane Bordas,et al.  Enriched finite elements and level sets for damage tolerance assessment of complex structures , 2006 .

[51]  C. C. Law,et al.  ParaView: An End-User Tool for Large-Data Visualization , 2005, The Visualization Handbook.

[52]  Frédéric Guyomarc'h,et al.  A Deflated Version of the Conjugate Gradient Algorithm , 1999, SIAM J. Sci. Comput..

[53]  Kurt Maute,et al.  A simple and efficient preconditioning scheme for heaviside enriched XFEM , 2013, 1312.6092.

[54]  Stéphane Bordas,et al.  Stable 3D extended finite elements with higher order enrichment for accurate non planar fracture , 2016 .

[55]  Frédéric Nataf,et al.  Spillane, N. and Dolean Maini, Victorita and Hauret, P. and Nataf, F. and Pechstein, C. and Scheichl, R. (2013) Abstract robust coarse spaces for systems of PDEs via generalized eigenproblems in the overlaps , 2018 .

[56]  P Kerfriden,et al.  A partitioned model order reduction approach to rationalise computational expenses in nonlinear fracture mechanics. , 2012, Computer methods in applied mechanics and engineering.

[57]  Wolfgang Fichtner,et al.  Efficient Sparse LU Factorization with Left-Right Looking Strategy on Shared Memory Multiprocessors , 2000 .

[58]  Richard Butler,et al.  dune-composites – A New Framework for High-Performance Finite Element Modelling of Laminates , 2017, 1707.04228.

[59]  X. Han,et al.  A novel multi-grid based reanalysis approach for efficient prediction of fatigue crack propagation , 2019, Computer Methods in Applied Mechanics and Engineering.

[60]  H. Waisman,et al.  AN ADAPTIVE DOMAIN DECOMPOSITION PRECONDITIONER FOR CRACK PROPAGATION PROBLEMS MODELED BY XFEM , 2013 .

[61]  Nicolas Moës,et al.  Extended finite element method in computational fracture mechanics: a retrospective examination , 2015, International Journal of Fracture.

[62]  Travis B. Fillmore,et al.  Preconditioned Conjugate Gradient Solvers for the Generalized Finite Element Method , 2017, Meshfree Methods for Partial Differential Equations IX.

[63]  Raymond S. Tuminaro,et al.  An algebraic multigrid approach to solve extended finite element method based fracture problems , 2013 .

[64]  Stefan Loehnert,et al.  A stabilization technique for the regularization of nearly singular extended finite elements , 2014 .