http://www.k-wave.org/forum/topic/simulating-a-plane-wave-with-a-linear-array-tranducer Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 00:06:25 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/simulating-a-plane-wave-with-a-linear-array-tranducer#post-7294 Mon, 16 Mar 2020 12:45:36 +0000 Bradley Treeby 7294@http://www.k-wave.org/forum/ <p>Hi bramth,</p> <p>Some of the wave energy will propagate through the PML. However, with appropriate settings this should be -80 dB less than the incoming wave (see the user manual if you want more details). This means your simulation is only accurate to a few decimal places, but that should be enough for most practical purposes. </p> <p>Brd. </p> http://www.k-wave.org/forum/topic/simulating-a-plane-wave-with-a-linear-array-tranducer#post-7290 Thu, 12 Mar 2020 14:46:25 +0000 bramth 7290@http://www.k-wave.org/forum/ <p>It seems setting the 'PlotScale' to either the max source intensity or to 'auto' solves the problem. Still, I wonder what happens in those low intensity regions. </p> http://www.k-wave.org/forum/topic/simulating-a-plane-wave-with-a-linear-array-tranducer#post-7289 Thu, 12 Mar 2020 14:03:16 +0000 bramth 7289@http://www.k-wave.org/forum/ <p>That is true, thanks. Is it by the way a plotting artefact that the B-mode scan example using linear array, when plotting the simulation, shows black intensities and something that seems PML penetration in the beginning? The "black" wave progression continues afterwards. See for the result after a few time steps <a href="https://imgur.com/a/4ME4RBZ" rel="nofollow">https://imgur.com/a/4ME4RBZ</a> and <a href="https://imgur.com/a/gNVeMEs" rel="nofollow">https://imgur.com/a/gNVeMEs</a>. </p> http://www.k-wave.org/forum/topic/simulating-a-plane-wave-with-a-linear-array-tranducer#post-7288 Wed, 11 Mar 2020 18:53:48 +0000 Bradley Treeby 7288@http://www.k-wave.org/forum/ <p>The particle velocity is typically several orders of magnitude smaller than the pressure. In 1D, they are related by the acoustic impedance (sound speed times density), so to get the same relative amplitude you will need to divide by the amplitude of the velocity source by the impedance. </p> http://www.k-wave.org/forum/topic/simulating-a-plane-wave-with-a-linear-array-tranducer#post-7287 Wed, 11 Mar 2020 18:42:08 +0000 bramth 7287@http://www.k-wave.org/forum/ <p>Dear Bradley,<br /> When using such a dipole transducerin a velocity-source setting (by setting source.p_mask to source.u_mask and source.p to source.ux), why can we not scale the input the same way as with a pressure source? (E.g. the transducer.source_strength, when scaling the tone_burst input with this value (and the impedance).) If I do scale the input, I get 'black' wavefronts (when not using C++ or CUDA implementations). Kind regards, Bram </p> http://www.k-wave.org/forum/topic/simulating-a-plane-wave-with-a-linear-array-tranducer#post-7237 Thu, 20 Feb 2020 15:07:04 +0000 jbishop 7237@http://www.k-wave.org/forum/ <p>Simulation of heating in an air-backed transducer and its substrate was of interest to me, and I too wondered just how well k-Wave could do with it not having access to any FEM solver. Without yet knowing how well the result matches experimental data, I can say so far it is certainly plausible looking. It required a lower CFL number of &lt; 0.1 before the air gap would stabilize. As a result of the low CFL number and fine resolution needed, compute times in 2D matlab were getting too long. I saw the other post on the c++/gpu release for 2D: that will be of immediate use to me. </p> http://www.k-wave.org/forum/topic/simulating-a-plane-wave-with-a-linear-array-tranducer#post-6460 Sun, 13 May 2018 20:10:37 +0000 Bradley Treeby 6460@http://www.k-wave.org/forum/ <p>Hi laughingrice,</p> <p>You could try using a velocity source instead of a pressure source. This will produce a dipole field, which will make the wave-field more directional. If the source is aligned with the grid, there will also be no response detected at the adjacent elements.</p> <p>Alternatively, you could model the source conditions using another model (e.g., FEM) that accounts for the actual source construction, record the acoustic field just in front of the source, and then input this into k-Wave.</p> <p>Accurately modelling the layered structure of the piezo, backing, and air in k-Wave will be pretty difficult due to the small scale structures and the large impedance mismatches.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/simulating-a-plane-wave-with-a-linear-array-tranducer#post-6447 Sun, 06 May 2018 02:55:40 +0000 laughingrice 6447@http://www.k-wave.org/forum/ <p>I'm trying to simulate a linear array clinical ultrasonic transducer emitting a plane wave (64 central elements active on a 128 element transducer). Simulation is running in 2D due to time and memory constraints at the moment.</p> <p>My problem is that just activating the pressure sources along the transducer face in a homogenous medium creates circular waves either at each element, if I make independent elements, or at the end points, if I just make a continuous line source. This waves travel across the transducer face and show up in the output traces.</p> <p>Running the same measurement with an actual probe, those waves traveling across the probe face can barely be seen, if at all, as signal does not propagate well sideways inside the probe. I'm guessing that it's a combination of independent piezo elements which I think are air separated, and attached to a matching layer.</p> <p>There is a little discussion about the physical structure here<br /> <a href="http://www.ndk.com/en/sensor/ultrasonic/basic02.html" rel="nofollow">http://www.ndk.com/en/sensor/ultrasonic/basic02.html</a><br /> and here<br /> <a href="http://www.echocardiographer.org/Echo%20Physics/BasicTransducers.html" rel="nofollow">http://www.echocardiographer.org/Echo%20Physics/BasicTransducers.html</a></p> <p>Anyone have an example of modeling a real probe in a simulation to get a more realistic wave field?<br /> How would I go about distributing the sources, and should I use pressure or velocity sources, if I do a piezo / air / matching layer construction in the domain?</p> <p>Thanks </p>