Frequency-domain Acoustic Wave Modeling Using a Hybrid Direct-iterative Solver Based on a Domain Decomposition Method
- Others:
- Géoazur (GEOAZUR 6526) ; Institut de Recherche pour le Développement (IRD)-Université Pierre et Marie Curie - Paris 6 (UPMC)-Université Nice Sophia Antipolis (1965 - 2019) (UNS) ; COMUE Université Côte d'Azur (2015-2019) (COMUE UCA)-COMUE Université Côte d'Azur (2015-2019) (COMUE UCA)-Institut national des sciences de l'Univers (INSU - CNRS)-Observatoire de la Côte d'Azur ; COMUE Université Côte d'Azur (2015-2019) (COMUE UCA)-Université Côte d'Azur (UCA)-Université Côte d'Azur (UCA)-Centre National de la Recherche Scientifique (CNRS)
- Laboratoire de Géophysique Interne et Tectonophysique (LGIT) ; Observatoire des Sciences de l'Univers de Grenoble (OSUG) ; Université Joseph Fourier - Grenoble 1 (UJF)-Institut polytechnique de Grenoble - Grenoble Institute of Technology (Grenoble INP )-Institut national des sciences de l'Univers (INSU - CNRS)-Institut national de recherche en sciences et technologies pour l'environnement et l'agriculture (IRSTEA)-Université Savoie Mont Blanc (USMB [Université de Savoie] [Université de Chambéry])-Centre National de la Recherche Scientifique (CNRS)-Université Joseph Fourier - Grenoble 1 (UJF)-Institut polytechnique de Grenoble - Grenoble Institute of Technology (Grenoble INP )-Institut national des sciences de l'Univers (INSU - CNRS)-Institut national de recherche en sciences et technologies pour l'environnement et l'agriculture (IRSTEA)-Université Savoie Mont Blanc (USMB [Université de Savoie] [Université de Chambéry])-Centre National de la Recherche Scientifique (CNRS)-Laboratoire Central des Ponts et Chaussées (LCPC)-Centre National de la Recherche Scientifique (CNRS)
- SEISCOPE
Description
Designing an efficient modeling tool is a key point for large 3D frequency-domain full-waveform inversion problems. We present a frequency-domain acoustic wave modeling using a hybrid direct-iterative solver based on a parallel domain decomposition method and Schur complement approach. The main interest of mixing solvers is to overcome the huge memory complexity of direct solvers while partially preserving efficient multi-RHS simulations and mitigating the iteration count in iterative solvers. To improve the convergence rate of the iterative solver, a preconditioning based on an additive Schwartz approach is used. Discretization of the Helmholtz equation is based on a parsimonious finite-difference method but the domain decomposition method could apply to any numerical scheme such as finite-element or finite-volume methods and to any media such as elastic, anisotropic ones ... To asses the efficiency of the hybrid approach, we computed simulations in the 2D Marmousi II and 3D SEG/EAGE Overthrust model, and compared results with that of a direct solver. ...
Abstract
International audience
Additional details
- URL
- https://hal.archives-ouvertes.fr/hal-00408386
- URN
- urn:oai:HAL:hal-00408386v1
- Origin repository
- UNICA