Summation-by-parts--simultaneous approximation term finite differences for anisotropic elasticity on curvilinear multiblock grids

Publications that use this code

1. Almquist, M., Dunham, E. M, Elastic wave propagation in anisotropic solids using energy-stable finite differences with weakly enforced boundary and interface conditions. J. Comput. Phys., 2021. doi:10.1016/j.jcp.2020.109842

Introduction

This MATLAB code demonstrates the advantages of a newly developed summation-by-parts--simultaneous approximation term (SBP-SAT) finite difference discretization for general anisotropic elasticity on curvilinear multiblock grids. The spatial discretization is energy stable and self-adjoint. The code reproduces figures 1-5 in Almquist and Dunham, 2021. The code is based on the sbplib repository.

Running the code

1. Clone the elastic-curvilinear repository: hg clone http://hg.code.sf.net/p/elastic-curvilinear/code elastic-curvilinear-code
2. Start MATLAB and navigate to the repository.
3. To run all computations, run the script runExperiments. Results from the computations will be stored in .mat-files.
4. To generate figures 1-5 in Almquist and Dunham, 2021, run the script generateFigures (after running runAllSimulations). The figures will appear in a folder named paperFigures.