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

by   Konstantinos Agathos, et al.

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.



There are no comments yet.


page 17

page 22

page 23

page 26

page 27

page 29


Sparse Approximate Multifrontal Factorization with Butterfly Compression for High Frequency Wave Equations

We present a fast and approximate multifrontal solver for large-scale sp...

Mixed Finite Element Discretization for Maxwell Viscoelastic Model of Wave Propagation

This paper considers semi-discrete and fully discrete mixed finite eleme...

Accelerating linear solvers for Stokes problems with C++ metaprogramming

The efficient solution of large sparse saddle point systems is very impo...

An iterative generalized Golub-Kahan algorithm for problems in structural mechanics

This paper studies the Craig variant of the Golub-Kahan bidiagonalizatio...

Machine learning applied in the multi-scale 3D stress modelling

This paper proposes a methodology to estimate stress in the subsurface b...

Parallel domain decomposition solvers for the time harmonic Maxwell equations

The time harmonic Maxwell equations are of current interest in computati...

Universal Scalable Robust Solvers from Computational Information Games and fast eigenspace adapted Multiresolution Analysis

We show how the discovery of robust scalable numerical solvers for arbit...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.