Thank you for providing such a convenient and free simulation tool.
Currently, I am using pstdElastic3D to simulate shear wave propagation in thin-layered phantoms and to calculate shear wave speed. However, it is known that when the thickness of the layer is smaller than the wavelength of the shear wave, guided waves can occur. In such cases, the shear modulus cannot be derived using the simplified formula: shear modulus = density*(shear wave speed)^2.
In pstdElastic3D.m, we input medium.sound_speed_shear to calculate lame parameter using
mu = medium.sound_speed_shear.^2 .* medium.density;.
I would like to ask if it is possible to:
1. Modify the input to directly provide the shear modulus (mu) instead of medium.sound_speed_shear while maintaining simulation accuracy.
2. Set flags.kelvin_voigt_model to false to simulate an anisotropic, incompressible, and nonlinear thin layer.
Best ragards,
Ting-Wei