k-Wave

1. Saptarshi Das
   

   Member
   Joined: Nov '15
   Posts: 1

   Hi,
   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.
   Thanks a lot,
   Saptarshi

   Posted 8 years ago #

2. Bradley Treeby
   

   Developer
   Joined: Mar '10
   Posts: 1,643

   Hi Saptarshi,

   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.

   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).

   Hope that helps,

   Brad,

   Posted 8 years ago #

RSS feed for this topic

Reply

You must log in to post.