http://www.k-wave.org/forum/topic/different-boundary-conditions-at-the-edges Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 12:13:53 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/different-boundary-conditions-at-the-edges#post-5700 Thu, 06 Oct 2016 22:33:55 +0000 Bradley Treeby 5700@http://www.k-wave.org/forum/ <p>Hi Saptarshi,</p> <p>Unfortunately, there is currently no way to do this, other than by imposing these conditions manually as you have read in other forum posts. I haven't tried this with the elastic code, but you could try increasing the number of grid points used to impose the boundary conditions until the reflection coefficient converges to the correct value. </p> <p>In the fluid code, we have experimented with using discrete cosine and sine transforms instead of the FFT to impose Dirichlet and Neumann boundary conditions, but this has never made it into the release version as MATLAB does not support the discrete trigonometric transform functionality of FFTW (the library it uses for the FFT).</p> <p>Hope that helps,</p> <p>Brad, </p> http://www.k-wave.org/forum/topic/different-boundary-conditions-at-the-edges#post-5695 Mon, 03 Oct 2016 16:03:17 +0000 Saptarshi Das 5695@http://www.k-wave.org/forum/ <p>Hi,<br /> Is there a simple way to impose Dirichlet, Neumann, Robin/mixed and reflecting boundary condition at the edges of the 3D Elastic simulation, similar to the absorbing boundary conditions through PML? I came across the discussions in this thread (Free surface condition) that defining extra sources along the boundaries and setting their amplitude to a fixed value may be a solution. In that case, how many extra layers in each direction would be ideal to impose such boundary conditions accurately, e.g. 10/20 grid points for PML in absorbing case. It would be really useful if these boundary conditions can be implemented in the k-wave elastic simulations.<br /> Thanks a lot,<br /> Saptarshi </p>