http://www.k-wave.org/forum/topic/long-run-time-problem Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 02:28:28 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/long-run-time-problem#post-7722 Wed, 05 Aug 2020 11:24:50 +0000 Bradley Treeby 7722@http://www.k-wave.org/forum/ <p>Hi Andrea,</p> <p>Can you put a colorbar on the images, and compare the size of the artefact to the size of the waves when they are still within the solid? I'm wondering whether this is some kind of instability, or just very small waves that haven't been absorbed by the PML. You could also try extending the bone region completely into the PML at both ends (so it is periodic without any corners).</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/long-run-time-problem#post-7678 Mon, 29 Jun 2020 15:34:11 +0000 Andrea Addison 7678@http://www.k-wave.org/forum/ <p>Dr.Treeby,</p> <p>Yes, it is after the waves have left the domain. I am not sure how I can compare them because the initial wave starts in the liquid and artefacts' region is in the solid medium. I am sure bone is not within the PML.</p> <p>Best regards,<br /> Andrea </p> http://www.k-wave.org/forum/topic/long-run-time-problem#post-7649 Sat, 27 Jun 2020 13:04:45 +0000 Bradley Treeby 7649@http://www.k-wave.org/forum/ <p>Hi Andrea,</p> <p>Is this after all the waves have left the domain? How big are the artefacts compare to the initial amplitude of your waves? If you make sure the bone is not within the PML, do you see the same effect?</p> <p>Brad </p> http://www.k-wave.org/forum/topic/long-run-time-problem#post-7618 Wed, 17 Jun 2020 06:36:50 +0000 Andrea Addison 7618@http://www.k-wave.org/forum/ <p>Hello Dr.Treeby,</p> <p>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: <a href="https://imgur.com/a/lMCWxwe" rel="nofollow">https://imgur.com/a/lMCWxwe</a></p> <p>Best regards,<br /> Andrea </p>