http://www.k-wave.org/forum/topic/free-surface-condition Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 06:46:29 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/free-surface-condition#post-4836 Tue, 11 Nov 2014 18:37:50 +0000 Bradley Treeby 4836@http://www.k-wave.org/forum/ <p>Hi Gilles,</p> <p>There is currently no general way in k-Wave to enforce a Dirichlet boundary condition (i.e., a pressure-release or solid-vacuum type interface), and thus generate true Rayleigh waves. However, if you're not using a stress source, one option to try is to use the <code>&#39;dirichlet&#39;</code> source option, and set your source input to be zero, e.g.,</p> <pre><code>source.s_mask = zeros(Nx, Ny); source.s_mask(1:Nx/4, :) = 1; source.sxx = zeros(1, kgrid.Nt); source.syy = source.sxx; source.sxy = source.syy; source.s_mode = &#39;dirichlet&#39;;</code></pre> <p>You could also try using very small values for the material parameters, or setting the stress to zero explicitly as explained in <a href="http://www.k-wave.org/forum/topic/room-enclouse-acoustic-simulation">this post</a>.</p> <p>If you come up with a nice solution, I'd be interested to hear about it!</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/free-surface-condition#post-4830 Thu, 06 Nov 2014 16:24:14 +0000 grandjean 4830@http://www.k-wave.org/forum/ <p>Hello,<br /> How to create a free-surface condition along a model's boudary so that Rayleigh waves can be generated?<br /> Many thanks<br /> Gilles </p>