http://www.k-wave.org/forum/topic/simulation-result-of-cylindrical-transducer Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 01:28:01 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/simulation-result-of-cylindrical-transducer#post-4642 Tue, 15 Jul 2014 11:35:05 +0000 Bradley Treeby 4642@http://www.k-wave.org/forum/ <p>Hi cge,</p> <p>Perhaps you could post your code and a sketch of the experiment, and I can take a look.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/simulation-result-of-cylindrical-transducer#post-4632 Tue, 08 Jul 2014 23:21:53 +0000 cge 4632@http://www.k-wave.org/forum/ <p>Moreover, for the simulation of flat disk and flat square transducer, I am still using the settings identical to the cylindrical one. i.e:<br /> attenuation coefficient=0.002 exponent power=2<br /> propagation speed=1540m/s; frequency=1.5Mhz<br /> dx=dy=dz=0.25mm</p> <p>But for these simulation, the result is quite normal </p> http://www.k-wave.org/forum/topic/simulation-result-of-cylindrical-transducer#post-4631 Tue, 08 Jul 2014 23:18:36 +0000 cge 4631@http://www.k-wave.org/forum/ <p>Another thing is that for the flat square and flat disk simulation, especially the flat disk one, the simulation results matches with measurement and theories of beamforming for flat disk simulation quite well.</p> <p>I think the cylindrical one(generating waves on its side wall as shown in previous figure) can be described by a square disk equiavelent model since if you cut the cylindrical one along one of its ridge and put the side wall on a plane, it is a square.</p> <p>My measurement support this. But still I am wondering why I can't directly get the simulation result for the cylindrical one </p> http://www.k-wave.org/forum/topic/simulation-result-of-cylindrical-transducer#post-4630 Tue, 08 Jul 2014 20:45:14 +0000 cge 4630@http://www.k-wave.org/forum/ <p>The frequency is 1.5MHz, with propagation speed to be 1540 m/s</p> <p>I set the dx=dy=dz=0.25mm, the wavelength is 1.02 mm, will this affect? </p> http://www.k-wave.org/forum/topic/simulation-result-of-cylindrical-transducer#post-4629 Tue, 08 Jul 2014 20:43:39 +0000 cge 4629@http://www.k-wave.org/forum/ <p>Hi Brad</p> <p>Yes, I am trying to replicate the experiment.</p> <p> Here I am using a cylindrical source and a square sensor mask. So that I can know the distribution along the circle at certain height.</p> <p>For the pressure plot, it is just the same.</p> <p>I don't understand what do you mean by the source frequency close to two points per wavelength. Does that mean there are two grid contained per wavelength?</p> <p>Best Regards </p> http://www.k-wave.org/forum/topic/simulation-result-of-cylindrical-transducer#post-4626 Tue, 08 Jul 2014 14:51:59 +0000 Bradley Treeby 4626@http://www.k-wave.org/forum/ <p>Hi cge,</p> <p>I'm not exactly sure what you are trying to simulate. Are you trying to replicate your experiment? From your figures above, it seems you have a square source and a cylindrical sensor mask. </p> <p>A couple of things to check:<br /> 1) Is your source frequency close to two points per wavelength?<br /> 2) What do you see if you plot the pressure rather than intensity?</p> <p>The current release of k-Wave only simulates longitudinal waves. This is valid if you are trying to model ultrasound waves in a fluid.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/simulation-result-of-cylindrical-transducer#post-4621 Tue, 08 Jul 2014 06:00:53 +0000 cge 4621@http://www.k-wave.org/forum/ <p>High Brad</p> <p>Sorry for the late reply.</p> <p>I am using a pressure source, a sine wave with an amplitude of 0.6MPa</p> <p>The source frequency is 1.5 MHz</p> <p>Media alpha power is set to 2</p> <p>Media attenuation coefficient is set to 0.0022</p> <p>Media propagation speed is set to 1540 m/s, density is the 1000 kg/m^3</p> <p>I didn't swiched off the PML at Z direction, I set the PMLAlpha to 10 with thickness at 3 direction being 10 respectively.</p> <p>The last firgure is for RMS pressure and average intensity has the similar profile.</p> <p>Recently, I carried out a measurement of the cylindrical transducer, the result of pressure distribution is quite similar to the simulation result of a flat square transducer. The measurement was carried out using hydrophone. </p> <p>For the distribution in the last figure, I remove the data of internal distribution in the figure. </p> <p>According to my literature study and measurement, it seems that the simulation above is for a shear wave , which experience huge attenuation in water, rather than longitudinal wave. However I didn't set anything about the type of wave generated. Also, the waves generated in simulation for flat square and flat disk transdcuers(with their surface perpendicular to Z direction), were longitudinal waves. In my measurement for real transducer, it seems that the wave generated is albo longitudinal waves. </p> <p>Does this related to the k-wave simulation tool? </p> http://www.k-wave.org/forum/topic/simulation-result-of-cylindrical-transducer#post-4586 Fri, 20 Jun 2014 08:42:15 +0000 Bradley Treeby 4586@http://www.k-wave.org/forum/ <p>Hi cge,</p> <p>A couple of things to check:<br /> 1) Are you using a pressure or velocity source?<br /> 2) What is your source frequency, and how close is it to two points per wavelength? Do you see the same behaviour if you use a lower frequency?<br /> 3) Have you switched off the PML in the z-direction, i.e., by setting the optional input <code>&#39;PMLAlpha&#39;, [2, 2, 0]</code>?</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/simulation-result-of-cylindrical-transducer#post-4585 Fri, 20 Jun 2014 01:41:42 +0000 cge 4585@http://www.k-wave.org/forum/ <p>Hello all</p> <p>I set all the parameters, like attenuation coefficient, as close as to the real water. Then, when I increase the excitation amplitude(in MPa), the uniformity of the intensity distribution becomes better.</p> <p>Is there any explanation for this? </p> http://www.k-wave.org/forum/topic/simulation-result-of-cylindrical-transducer#post-4584 Thu, 19 Jun 2014 05:58:13 +0000 cge 4584@http://www.k-wave.org/forum/ <p>Sorry for the fail of uploading. Here is another try from imageshack</p> <p><img src="http://imageshack.com/a/img843/3716/e2bw9.jpg" /><br /> <img src="http://imageshack.com/a/img822/3453/9yfmd.jpg" /><br /> <img src="http://imageshack.com/a/img819/4809/fylr.jpg" /> </p> http://www.k-wave.org/forum/topic/simulation-result-of-cylindrical-transducer#post-4583 Thu, 19 Jun 2014 05:55:43 +0000 cge 4583@http://www.k-wave.org/forum/ <p>Sorry for the fail of uploading. Here is another try from imageshack</p> <p><img src="https://imageshack.com/i/nf1n3wyj" /><br /> <img src="https://imageshack.com/i/npnqutj" /><br /> <img src="https://imageshack.com/i/np9dsdj" /> </p> http://www.k-wave.org/forum/topic/simulation-result-of-cylindrical-transducer#post-4582 Thu, 19 Jun 2014 05:47:07 +0000 cge 4582@http://www.k-wave.org/forum/ <p>Hi everyone</p> <p>I simulate the ultrasound pressure and intensity distribution for a cylindrical transducer, the transducer, sensor mask, and the result is like below.</p> <p>For the profile of distribution, it is really interesting. Does anyone know is there any principle of ultrasound physics that explain this? </p> <p><img src="http://xiangce.baidu.com/picture/detail/830b2e08a41f08319b6e8c9010c232e7bcbc64d9?from=dialog" /></p> <p><img src="http://xiangce.baidu.com/picture/detail/edc8a72b826c7ade630b73d85d13e02d81df52d9?from=dialogg" /></p> <p><img src="http://xiangce.baidu.com/picture/detail/247283539dd3d0ed1f0755126823826c7bde5d82?from=dialog" /> </p>