http://www.k-wave.org/forum/topic/working-within-memory-limits Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 01:24:22 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/working-within-memory-limits#post-3962 Mon, 14 Oct 2013 16:30:12 +0000 Bradley Treeby 3962@http://www.k-wave.org/forum/ <p>Hi Charles,</p> <p>For a CW source, I think the situation is slightly more complicated. If you have a homogeneous medium, when you enforce the source in your second simulation as a Dirichlet boundary condition, the resulting wavefield downstream from the 2D plane will be the same as if you had run one big simulation. </p> <p>However, in a heterogeneous medium, the waves scattered from your imaging target will then also be reflected from the source plane over which the Dirichlet boundary condition is being applied. This is easiest to imagine in the case of the source pressure being zero. This will result in a pressure of zero being enforced across the 2D plane, which will reflect the sound waves. </p> <p>You could avoid this by using a source signal that has finished being transmitted before any reflections reach the source plane.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/working-within-memory-limits#post-3961 Mon, 14 Oct 2013 16:16:36 +0000 chuck_c 3961@http://www.k-wave.org/forum/ <p>Thank you for the reply. I think this should help me get started. I would possibly like to measure the sound reflected back at some point, but I will try to focus on one thing at a time! It seems the reverse of your suggestion (ie. propagate back via a plane on the left side of the second space into the first space) might be a possible solution. Also, I am using a ~10 msec pulse, so I think it can be viewed as somewhat continuous. </p> http://www.k-wave.org/forum/topic/working-within-memory-limits#post-3960 Mon, 14 Oct 2013 14:55:01 +0000 Bradley Treeby 3960@http://www.k-wave.org/forum/ <p>Hi Charles,</p> <p>One possibility would be to setup a simulation just containing the transducer in water, and record the time varying pressure over a 2D plane on one side of the domain, say the right. In a new simulation containing your imaging target, you could then enforce this pre-recorded pressure over an identical 2D plane as a Dirichlet boundary condition on the left of the domain. This should give you the same acoustic field as if you had performed one larger simulation (keeping in mind the time-staggering between inputs and outputs as described in the k-Wave manual).</p> <p>Do you want to measure the reflected sound back at the transducer surface, or just look at the scattered waves around your imaging target? Also, are you using a CW or pulsed source?</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/working-within-memory-limits#post-3954 Fri, 11 Oct 2013 15:56:17 +0000 chuck_c 3954@http://www.k-wave.org/forum/ <p>Hi,</p> <p>I am trying to simulate a scenario using a GPU but the space is too large for my GPU's memory. The scenario is a spherically focused ultrasound transducer with a large amount of water stand-off, so I think the propagation from the transducer face through the initial water section should be quite similar regardless of my imaging target. I was wondering whether it would be possible to simulate the transducer one time and then use the output from that as input into the smaller target space. I am open to other solutions to this problem but this seemed intuitive.</p> <p>Thanks for the nice simulation package.</p> <p>-Charles </p>