Hello Dr.Treeby,
Thanks for the useful toolbox. I run a similar simulation medium in this paper (The Contribution of Shear Wave Absorption to Ultrasound Heating in Bones: Coupled Elastic and Thermal Modeling). After I run a 100-microsecond elastic simulation that has bone and water mediums, I get a problem around the boundary. I didn't get the same problem when I run for a 50-microsecond. I am not sure if I get this problem because I put the bone and PML adjacent. Simulation parameters: dx=dy=50e-6, cfl=0.055, same bone properties as you mentioned in your paper. The maximum pressure figure is here: https://imgur.com/a/lMCWxwe
Best regards,
Andrea