http://www.k-wave.org/forum/topic/validation-of-the-tissue-velocity-recorded-by-sensor Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 03:13:19 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/validation-of-the-tissue-velocity-recorded-by-sensor#post-8125 Mon, 19 Apr 2021 13:46:57 +0000 Bradley Treeby 8125@http://www.k-wave.org/forum/ <p>Hi Ekaterina,</p> <p>I couldn't view the image, however, if you experience an instability you could try reducing the time step size.</p> <p>On the topic of reflections, I'm not sure these will match any physically realistic scenario, as you are imposing a time-varying boundary condition, thus the reflections will depend on the values for the velocity that you are imposing at each time step in the simulation.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/validation-of-the-tissue-velocity-recorded-by-sensor#post-8108 Thu, 25 Mar 2021 19:13:33 +0000 Elseliv 8108@http://www.k-wave.org/forum/ <p>I actually tried using the Dirichlet condition but this brought instability and the simulation got totally crashed. Here is a link with 2 plots for Dirichlet and NoDirichlet cases <a href="https://we.tl/t-Dft5lKj5vd" rel="nofollow">https://we.tl/t-Dft5lKj5vd</a> and it depicts in the top row - source velocities as I imposed them and on the 2-nd and 3-d rows - velocities measure by the sensor exactly next to the source. Could you then assume what might be a reason for such a crash?<br /> Concerning reflections, the source I use is located at the boundary between tissue and water so reflections are expected to happen anyway. </p> http://www.k-wave.org/forum/topic/validation-of-the-tissue-velocity-recorded-by-sensor#post-8106 Thu, 25 Mar 2021 17:52:22 +0000 Bradley Treeby 8106@http://www.k-wave.org/forum/ <p>Hi Ekaterina,</p> <p>By default, a velocity source is a force source. It will inject or "add" to the field that is already there, and thus the actual velocity will be higher than the velocity you specify. If you want the source to act like a boundary condition (exactly impose the values you specify), then you can set <code>source.u_mode = &#39;dirichlet&#39;</code>. Keep in mind that in this case, the imposed boundary condition will also reflect any ways coming back to the source.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/validation-of-the-tissue-velocity-recorded-by-sensor#post-8091 Wed, 24 Mar 2021 01:32:51 +0000 Elseliv 8091@http://www.k-wave.org/forum/ <p>Hello Dr.Treeby,</p> <p>Thank you for the great toolbox!</p> <p>I`m doing shear wave simulation in a thin tissue layer using pstdElastic2D. As a source.ux, I use a line where each element imposes one cycle of a sine wave with a frequency of 50 Hz and max amplitude of 0.03 m/s. Then I record u_split_field with a line of sensors perpendicular to the source line and get values significantly higher (sensor_data.ux_split_s ~ 0.1 m/s sensor_data.uy_split_s ~ 0.5 m/s), than those imposed by source. Even for the sensors placed tight to the source.</p> <p>What could become a reason for differences in the velocities imposed and observed? Is there a way to validate the magnitude of output data? May there be a problem with output units or something else? </p> <p>Best regards,<br /> Ekaterina </p>