Hybrid Fast Sweeping Methods for Anisotropic Eikonal Equation in Two-Dimensional Tilted Transversely Isotropic Media
Autor: | Guangnan Huang, Songting Luo |
---|---|
Rok vydání: | 2020 |
Předmět: |
Numerical Analysis
Point source Eikonal equation Applied Mathematics Mathematical analysis General Engineering Degrees of freedom (statistics) 010103 numerical & computational mathematics 010502 geochemistry & geophysics 01 natural sciences Theoretical Computer Science Computational Mathematics Singularity Computational Theory and Mathematics Factorization Discontinuous Galerkin method Transverse isotropy Convergence (routing) 0101 mathematics Software 0105 earth and related environmental sciences Mathematics |
Zdroj: | Journal of Scientific Computing. 84 |
ISSN: | 1573-7691 0885-7474 |
DOI: | 10.1007/s10915-020-01280-3 |
Popis: | We present hybrid fast sweeping methods for computing first-arrival traveltime of the qP, qSV and qSH waves in two-dimensional tilted transversely isotropic media, based on solving the anisotropic eikonal equation. A factorization approach is applied to resolve the source singularity near the point source, which leads to a factored anisotropic eikonal equation whose solutions can be computed with high accuracy. The proposed methods solve the factored equation in a neighborhood of the point source with the size of the neighborhood independent of the mesh, and solve the original equation outside the neighborhood. The methods enjoy all the appealing features, such as efficiency, accuracy and convergence, of the usual fast sweeping method. Furthermore, the “super-convergence” property of the first-order fast sweeping method, i.e., both its numerical solution and gradient are first-order accurate, allows us to design a second-order fast sweeping method based on a linear discontinuous Galerkin formulation. As a post-processing procedure of the first-order method, the second-order method reduces the local degrees of freedom from three to one in the linear discontinuous Galerkin formulation, which implies a simple local updating formula, hence an efficient second-order scheme. Numerical experiments are presented to demonstrate the proposed methods. |
Databáze: | OpenAIRE |
Externí odkaz: |