http://www.k-wave.org/forum/topic/room-enclouse-acoustic-simulation Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 00:28:39 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/room-enclouse-acoustic-simulation#post-3769 Mon, 08 Jul 2013 14:05:13 +0000 Bradley Treeby 3769@http://www.k-wave.org/forum/ <p>Hi Adrian,</p> <p>Unfortunately, there is no easy way to do this using the functionality of the current release of the code (V1.0). The closest would be to set a region around the outside of your domain with a much higher density than the interior. </p> <p>However, there are some things you could try if you are happy doing some code edits of the simulation functions. For Dirichlet-type boundaries (sound-soft), set the pressure at the edges of the domain to zero, immediately after the pressure is calculated, i.e., in 2D</p> <pre><code>p(1, :) = 0; p(end, :) = 0; p(:, 1) = 0; p(:, end) = 0;</code></pre> <p>Similarly, for Neumann-type boundaries (sound-hard), set the normal velocity at the edges of the domain to zero, immediately after the velocity is updated, i.e., in 2D</p> <pre><code>ux_sgx(1, :) = 0; ux_sgx(end, :) = 0; uy_sgy(:, 1) = 0; uy_sgy(:, end) = 0;</code></pre> <p>Keep in mind I haven't properly tested this, so it's probably best to do some testing to make sure it's behaving as expected. You may need to use a layer more than one-grid point thick, i.e.,</p> <pre><code>ux_sgx(1:2, :) = 0; ux_sgx(end-1:end, :) = 0; uy_sgy(:, 1:2) = 0; uy_sgy(:, end-1:end) = 0;</code></pre> <p>Let me know if you run into any problems.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/room-enclouse-acoustic-simulation#post-3761 Tue, 02 Jul 2013 12:09:33 +0000 ace2112 3761@http://www.k-wave.org/forum/ <p>Hi,<br /> I am trying to figure out how to set and enclosure to simulate wave propagation in a rectangular room.</p> <p>Any one have an idea? </p> <p>Best regards<br /> Adrian </p>