Premium
A semi‐local spectral/hp element solver for linear elasticity problems
Author(s) -
Yu Yue,
Bittencourt Marco L.,
Karniadakis George Em
Publication year - 2014
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.4739
Subject(s) - preconditioner , conjugate gradient method , solver , linear elasticity , discretization , finite element method , elasticity (physics) , linear system , augmented lagrangian method , mathematical optimization , computer science , mathematics , algorithm , mathematical analysis , structural engineering , physics , thermodynamics , engineering
SUMMARY We develop an efficient semi‐local method for speeding up the solution of linear systems arising in spectral/hp element discretization of the linear elasticity equations. The main idea is to approximate the element‐wise residual distribution with a localization operator we introduce in this paper, and subsequently solve the local linear system. Additionally, we decouple the three directions of displacement in the localization operator, hence enabling the use of an efficient low energy preconditioner for the conjugate gradient solver. This approach is effective for both nodal and modal bases in the spectral/hp element method, but here, we focus on the modal hierarchical basis. In numerical tests, we verify that there is no loss of accuracy in the semi‐local method, and we obtain good parallel scalability and substantial speed‐up compared to the original formulation. In particular, our tests include both structure‐only and fluid‐structure interaction problems, with the latter modeling a 3D patient‐specific brain aneurysm. Copyright © 2014 John Wiley & Sons, Ltd.