Fully scalable solver for frequency-domain visco-elastic wave equations in 3D heterogeneous media: A controllability approach

We develop a controllability strategy for the computation of frequency-domain solutions of the 3D visco-elastic wave equation, in the perspective of seismic imaging applications. We generalize the controllability results for such equations beyond the sound-soft scattering (obstacle) problem. We deta...

Full description

Saved in:
Bibliographic Details
Published in:Journal of computational physics Vol. 468; p. 111514
Main Authors: Tang, Jet Hoe, Brossier, Romain, Métivier, Ludovic
Format: Journal Article
Language:English
Published: Cambridge Elsevier Inc 01.11.2022
Elsevier Science Ltd
Elsevier
Subjects:
Algorithms
Computational physics
Computer Science
Conjugate gradient
Controllability
Discretization
Domain decomposition methods
Elastic wave equations
Elastic waves
Frequency domain analysis
Full waveform inversion
Geophysics
Iterative methods
Modeling and Simulation
Parallelization
Physics
Runge-Kutta method
Spectral element
Viscoelasticity
Wave equations
Waveforms
Controllability
Conjugate gradient
Full waveform inversion
Parallelization
Elastic wave equations
Spectral element
ISSN:0021-9991, 1090-2716
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:We develop a controllability strategy for the computation of frequency-domain solutions of the 3D visco-elastic wave equation, in the perspective of seismic imaging applications. We generalize the controllability results for such equations beyond the sound-soft scattering (obstacle) problem. We detail the conjugate gradient implementation and show how an inner elliptic problem needs to be solved to compute the Riesz representative of the gradient at each iteration. We select a spectral-element spatial discretization and a fourth-order Runge-Kutta time discretization. We implement the controllability method in the framework of the SEM46 full waveform modeling and inversion software, to inherit for its excellent scalability which relies on an efficient domain decomposition algorithm. We perform a series of numerical experiments to validate the approach and illustrate its scalability up to more than fifteen hundred cores. In this case, with an elapsed time of less than 50 minutes, we solve a problem on a cubic domain containing up to 160 wavelengths in each direction, involving more than 1.7 billion unknowns. •Efficient numerical solution of frequency-domain visco-elastic wave equation.•Controllability methods (CM) extended to general visco-elastic problem beyond sound-soft scattering problems.•CM combined with spectral element method, explicit time scheme and parallelization through domain decomposition.•CM offer robust, non-intrusive approach for time-harmonic wave equations.•CM achieve strong scalability on massively parallel architectures.
Bibliography:ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 14
ISSN:0021-9991
1090-2716
DOI:10.1016/j.jcp.2022.111514