Hello,
How to create a free-surface condition along a model's boudary so that Rayleigh waves can be generated?
Many thanks
Gilles
k-Wave
A MATLAB toolbox for the time-domain
simulation of acoustic wave fields
Free surface condition
(2 posts) (2 voices)-
Posted 10 years ago #
-
Hi Gilles,
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
'dirichlet'
source option, and set your source input to be zero, e.g.,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 = 'dirichlet';
You could also try using very small values for the material parameters, or setting the stress to zero explicitly as explained in this post.
If you come up with a nice solution, I'd be interested to hear about it!
Brad.
Posted 10 years ago #
Reply
You must log in to post.