Elastic Mode Solver

2+1D eigenmode solver for linear elastic waves in isotropic media.

This Matlab program uses a finite-difference, split-step grid approach to solve for the eigenmodes of elastic wave propagation in a two-dimensional cross-section of a z-invariant medium. The mode solver assumes linear, elastic, isotropic media without viscous damping. Several boundary conditions have been implemented: free, fixed, and symmetry/anti-symmetry; a perfectly matched layer (absorbing/radiation boundary condition) is also available. The solver and all boundary conditions have been tested against COMSOL and demonstrated to return eigenmodes/eigenfrequencies within ~5% of the COMSOL solutions. See the "Scripts" folder for example simulations, including Rayleigh waves, suspended beam modes, and leaky waveguide modes.