| Texto completo | |
| Autor(es): |
Ravelo, Fernando V.
;
Schreiber, Martin
;
Peixoto, Pedro S.
Número total de Autores: 3
|
| Tipo de documento: | Artigo Científico |
| Fonte: | COMPUTATIONAL GEOSCIENCES; v. N/A, p. 21-pg., 2024-09-04. |
| Resumo | |
Seismic imaging is a major challenge in geophysics with broad applications. It involves solving wave propagation equations with absorbing boundary conditions (ABC) multiple times. This drives the need for accurate and efficient numerical methods. This study examines a collection of exponential integration methods, known for their good numerical properties on wave representation, to investigate their efficacy in solving the wave equation with ABC. The purpose of this research is to assess the performance of these methods. We compare a recently proposed Exponential Integration based on Faber polynomials with well-established Krylov exponential methods alongside a high-order Runge-Kutta scheme and low-order classical methods. Through our analysis, we found that the exponential integrator based on the Krylov subspace exhibits the best convergence results among the high-order methods. We also discovered that high-order methods can achieve computational efficiency similar to low-order methods while allowing for considerably larger time steps. Most importantly, the possibility of undertaking large time steps could be used for important memory savings in full waveform inversion imaging problems. (AU) | |
| Processo FAPESP: | 21/06176-0 - Métodos numéricos para uma nova geração de modelos de previsão de tempo e clima |
| Beneficiário: | Pedro da Silva Peixoto |
| Modalidade de apoio: | Auxílio à Pesquisa - Programa de Pesquisa sobre Mudanças Climáticas Globais - Jovens Pesquisadores - Fase 2 |