http://www.k-wave.org/forum/topic/simulation-of-an-artificial-source Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 12:08:33 +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-of-an-artificial-source#post-5663 Fri, 09 Sep 2016 21:06:39 +0000 Bradley Treeby 5663@http://www.k-wave.org/forum/ <p>Hi Oscar,</p> <p>I've not tried it using a velocity source, but it's certainly possible if you record and re-transmit the pressure using a Dirichlet boundary condition. I created a simple <a href="http://www.k-wave.org/downloads/forum_dirichlet_divided_simulation_1D.m">example</a> in 1D to demonstrate. You can apply the same idea in 2D and 3D.</p> <p>Note, re-transmitting the pressure over a single plane is not sufficient to perfectly recover the spatial gradients, so your output might have a small phase offset compared to a simulation without the split. One of the forum users discovered that this could be overcome if you record and transmit over a plane with a thickness greater than one grid point (see <a href="http://www.k-wave.org/forum/topic/simulation-layer-by-layer">this post</a> and the <a href="dx.doi.org/10.1186/s40349-016-0056-9">related paper</a>).</p> <p>You can test this in the example above by increasing the value of <code>overlap</code>.</p> <p>Hope that helps,</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/simulation-of-an-artificial-source#post-5648 Fri, 26 Aug 2016 08:10:46 +0000 osarbui 5648@http://www.k-wave.org/forum/ <p>Hi Dr. Treeby and Dr. Cox:</p> <p>First, congratulations for your great toolbox, I started using some months ago and it has been really helpful.</p> <p>I am working on a 1D simulation in which I locate the source in one extreme of the grid (outside PML layer) and measure the maximum pressure values along the area, and the time-varying particle velocity in one point in the middle of the grid. After that, in a second simulation I will place the source in that middle point using those time-varying results as inputs for particle velocity and run it (i.e. like creating an artificial source of the first case). The idea is to obtain similar results in the second half of the grid (from the middle point to the end) as those from the first simulation, but unfortunately I haven’t got good results yet.</p> <p>I would like to ask you if such case would be feasible, or if some initial conditions are required to do it. And if it could be reproduced in 2D and 3D cases (with a line and a plane, respectively). Thank you.</p> <p>Best regards,<br /> Oscar </p>