A Rotated Staggered Grid (RSG) finite-difference elastic wave solver based on Devito

<p>Built upon the open-source Python finite-difference framework, Devito, this Rotated Staggered Grid elastic wave solver is designed to model seismic wave propagation in intricate anisotropic media.</p>

<p>Key Components:</p>

<ul>
<li>
<p><code>model.py</code>: A refined file derived from Devito's original <code>model.py</code>. It introduces a new object, <code>AnisoSeismicModel</code>, tailored to support the input of 21 distinct elastic parameters for the grid.</p>
</li>
<li>
<p><code>wavesolver.py</code>: Solver object provides the operators for forward modeling computations. It encapsulates both temporal and spatial discretizations pertinent to a specified problem setup.</p>
</li>
<li>
<p><code>operators.py</code>: Construct Method for the forward modeling operator handling source excitation and monitoring component recorded by receivers.</p>
</li>
<li>
<p><code>stiffnessoperator.py</code>: Serving as an operator for the stiffness tensor, this object accommodates isotropy, Vertical Transverse Isotropy (VTI), and 3D rotations. Users can either opt for manual input of the 21 elastic parameters or utilize the automated setup. When paired with a mask, it aids in crafting the <code>AnisoSeismicModel</code>.</p>
</li>
</ul>

<p>Getting Started:</p>

<p>For a practical introduction and hands-on experience, refer to the tutorial Jupyter notebook <code>Tutorial.ipynb</code>. This guide walks users through the application of the package, showcasing wave propagation modeling within a two-layered structure (comprising an upper water layer and a lower TTI media) in 2D.</p>

<p>Open Research Archive (<code>OpenResearchArchive.zip</code>) contains the result and Python code to conduct the numerical tests in terms of stability, numerical dispersion and performance benchmark described in the manuscript: Zhang, O. and Schmitt, D. R. (2023) An optimized 2D/3D Finite-difference Seismic Wave Propagator Using Rotated Staggered Grid for Complex Elastic Anisotropic Structures. <i>ESS Open Archive.</i> November, 2023. File <code>Readme.md</code> details the content in each subdirectory.</p>