http://www.k-wave.org/forum/topic/attenuation-in-elastic-wave-simulations Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 01:00:36 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/attenuation-in-elastic-wave-simulations#post-6295 Tue, 20 Feb 2018 14:18:37 +0000 Bradley Treeby 6295@http://www.k-wave.org/forum/ <p>Hi Yishi,</p> <p>The function <code>pstdElastic2D</code> solves a relatively simple visco-elastic wave model as described in <a href="http://bug.medphys.ucl.ac.uk/papers/2014-Treeby-IEEEIUS.pdf">this paper</a>. This doesn't use the fractional Laplacian to implement the absorption, and gives frequency squared absorption. We have experimented with this for the elastic code (as described <a href="http://bug.medphys.ucl.ac.uk/papers/2014-Treeby-JASA.pdf">here</a>), but the computational load is very heavy, so we haven't released this code.</p> <p>If you're working with continuous wave inputs (or at least something with a centre frequency), one option would be to adjust the <code>alpha_coeff</code> values to give the correct absorption at the frequency of interest.</p> <p>Hope that helps,</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/attenuation-in-elastic-wave-simulations#post-6287 Fri, 16 Feb 2018 01:42:14 +0000 stone7121 6287@http://www.k-wave.org/forum/ <p>Hi Brad,</p> <p>First of all, great effort creating such a useful toolbox. The simulation works really well, and I want to thank you for putting the resources developing it. </p> <p>I am simulating acoustic wave field in an elastic medium that is biological in nature. The absorption power does not necessarily follow the frequency square, in fact, it is very close to 1. I know that you used fractional Laplacian to simulate absorption, but is there a way to customize the value of absorption power in the pstdElastic2D solver? I greatly appreciate your response. </p> <p>Yishi </p> http://www.k-wave.org/forum/topic/attenuation-in-elastic-wave-simulations#post-5998 Mon, 19 Jun 2017 20:56:59 +0000 Bradley Treeby 5998@http://www.k-wave.org/forum/ <p>Hi Bill,</p> <p>k-Wave implements the Kelvin-Voigt model of viscoelasticity, which gives absorption with a frequency squared dependence (unless you're operating at very high frequencies / absorptions). You can set these as explained in the <code>Plane Wave Absorption Example</code> under the <code>Elastic Wave Propagation</code> heading in the MATLAB help. As with the fluid model, very high values of absorption will require smaller time steps to remain stable.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/attenuation-in-elastic-wave-simulations#post-5985 Wed, 14 Jun 2017 13:01:10 +0000 bwalker000 5985@http://www.k-wave.org/forum/ <p>MATLAB Version: 9.0.0.341360 (R2016a)<br /> Operating System: Mac OS X Version: 10.12.5 Build: 16F73<br /> Java Version: Java 1.7.0_75-b13 with Oracle Corporation Java HotSpot(TM) 64-Bit Server VM mixed mode<br /> ----------------------------------------------------------------------------------------------------<br /> MATLAB Version 9.0 (R2016a)<br /> Signal Processing Toolbox Version 7.2 (R2016a)<br /> k-Wave Toolbox Version 1.1 </p> http://www.k-wave.org/forum/topic/attenuation-in-elastic-wave-simulations#post-5981 Mon, 12 Jun 2017 20:18:21 +0000 bwalker000 5981@http://www.k-wave.org/forum/ <p>I am unclear from the documentation whether attenuation (shear or longitudinal) is supported in elastic wave simulations. Empirically I see that if I set a fairly high attenuation I get unstable solutions.</p> <p>Any insight about applying attenuation in elastic simulations would be very welcome.</p> <p>Thanks,<br /> Bill Walker </p>