The classical overlapping Schwarz algorithm is here extended to the spectral element discretization of linear elastic problems, for both homogeneous and heterogeneous compressible materials. The algorithm solves iteratively the resulting preconditioned system of linear equations by the conjugate gradient or GMRES methods. The overlapping Schwarz preconditioned technique is then applied to the numerical approximation of elastic waves with spectral elements methods in space and implicit Newmark time advancing schemes. The results of several numerical experiments, for both elastostatic and elastodynamic problems, show that the convergence rate of the proposed preconditioning algorithm is independent of the number of spectral elements (scalability), is independent of the spectral degree in case of generous overlap, otherwise it depends inversely on the overlap size. Some results on the convergence properties of the spectral element approximation combined with Newmark schemes for elastic waves are also presented.
Overlapping Schwarz and spectral element methods for linear elasticity and elastic waves / L.F. PAVARINO, E. ZAMPIERI. - In: JOURNAL OF SCIENTIFIC COMPUTING. - ISSN 0885-7474. - 27:1-3(2006), pp. 51-73. ((Intervento presentato al 6. convegno International Conference on Spectral and High Order Methods (ICOSAHOM) tenutosi a Providence nel 2004.
Overlapping Schwarz and spectral element methods for linear elasticity and elastic waves
L.F. PAVARINOPrimo
;E. ZAMPIERIUltimo
2006
Abstract
The classical overlapping Schwarz algorithm is here extended to the spectral element discretization of linear elastic problems, for both homogeneous and heterogeneous compressible materials. The algorithm solves iteratively the resulting preconditioned system of linear equations by the conjugate gradient or GMRES methods. The overlapping Schwarz preconditioned technique is then applied to the numerical approximation of elastic waves with spectral elements methods in space and implicit Newmark time advancing schemes. The results of several numerical experiments, for both elastostatic and elastodynamic problems, show that the convergence rate of the proposed preconditioning algorithm is independent of the number of spectral elements (scalability), is independent of the spectral degree in case of generous overlap, otherwise it depends inversely on the overlap size. Some results on the convergence properties of the spectral element approximation combined with Newmark schemes for elastic waves are also presented.File | Dimensione | Formato | |
---|---|---|---|
JSC2006.pdf
accesso riservato
Tipologia:
Publisher's version/PDF
Dimensione
698.79 kB
Formato
Adobe PDF
|
698.79 kB | Adobe PDF | Visualizza/Apri Richiedi una copia |
Pubblicazioni consigliate
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.