BibTex format
@inproceedings{Caunt:2021,
author = {Caunt, E and Nelson, R and Luporini, F and Gorman, G},
pages = {1152--1156},
title = {GENERALISED ALGORITHM AND IMPLEMENTATION OF TOPOGRAPHY WITHIN FINITE DIFFERENCE WAVE SOLVERS},
year = {2021}
}
Publications
Citation
RIS format (EndNote, RefMan)
TY - CPAPER
AB - Partial differential equation solvers based on the finite-difference method have for many years been a keystone of seismic processing and modelling applications, commonly found in practical migration and full-waveform inversion methods. Sharp density contrasts within the computational domain have potential to introduce numerical error: problematic when introducing topography to a seismic model. Including a simple step change in density to approximate an air layer compromises both stability and numerical accuracy, often requiring smoothing of the contrast, and inducing both first and second order errors in space. Topography can instead be implemented via an immersed boundary conforming to the surface. This is achieved by extrapolating the wavefield across the boundary to find solution values at necessary external nodes. As this process is confined to the pre-processing step, it has negligible effect on the computational cost of the simulation. Devitoboundary is a tool in its early stages of development, intended to compliment Devito as a user-friendly means of including immersed boundaries in practical applications. 3D immersed boundaries can be constructed from irregularly sampled topography point clouds, via Delaunay triangulation coupled with a 1D extrapolation scheme. The result is a stable, error-free boundary which can be readily integrated with Devito models.
AU - Caunt,E
AU - Nelson,R
AU - Luporini,F
AU - Gorman,G
EP - 1156
PY - 2021///
SP - 1152
TI - GENERALISED ALGORITHM AND IMPLEMENTATION OF TOPOGRAPHY WITHIN FINITE DIFFERENCE WAVE SOLVERS
ER -