HR: 0800h
AN: S21C-0717    [Abstracts]
TI: Assessment of a Multidomain Chebyshev Method for seismic wave modeling
AU: * Danecek, P
EM: danecek@bo.ingv.it
AF: Istituto Nazionale di Geofisica e Vulcanologia, Via Donato Creti 12, Bologna, BO 40128, Italy
AU: Seriani, G
EM: gseriani@inogs.it
AF: Istituto Nazionale di Oceanografia e di Geofisica Sperimentale, Borgo Grotta Gigante 42/C, Sgonico, TS 34010, Italy
AB: Simulations of seismic wave propagation for large areas and realistic 3D Earth models require the use of multi- processor computers in order to obtain a relevant range of frequencies. Pseudo-spectral techniques, like other high order methods, provide highly accurate solutions while reducing the grid resolution. Due to the global character of the discrete wave operator, these techniques suffer from very low computational scalability on massive parallel computers and the wave field may suffer from random non-causal errors. The Overlapping Multidomain Chebyshev Method is particularly well suited for developing efficient and scalable parallel algorithms using a domain decomposition approach. The pseudospectral Chebyshev differential operator is applied to a large number of smaller subdomains and is local to each of them. Only a reduced amount of data needs to be exchanged in order to couple adjacent subdomains using a small overlap region. In the present work we use this method in order to implement a velocity-stress formulation of the elasto-dynamic wave equation in cartesian and spherical coordinates and we assess the obtained solutions in dependence of various parameters and implementation alternatives for timestepping scheme and boundaries.
DE: 0545 Modeling (4255)
DE: 0560 Numerical solutions (4255)
DE: 3294 Instruments and techniques
DE: 7290 Computational seismology
SC: Seismology [S]
MN: 2007 Fall Meeting