k-Wave

1. osarbui
   

   Member
   Joined: Apr '16
   Posts: 1

   Hi Dr. Treeby and Dr. Cox:

   First, congratulations for your great toolbox, I started using some months ago and it has been really helpful.

   I am working on a 1D simulation in which I locate the source in one extreme of the grid (outside PML layer) and measure the maximum pressure values along the area, and the time-varying particle velocity in one point in the middle of the grid. After that, in a second simulation I will place the source in that middle point using those time-varying results as inputs for particle velocity and run it (i.e. like creating an artificial source of the first case). The idea is to obtain similar results in the second half of the grid (from the middle point to the end) as those from the first simulation, but unfortunately I haven’t got good results yet.

   I would like to ask you if such case would be feasible, or if some initial conditions are required to do it. And if it could be reproduced in 2D and 3D cases (with a line and a plane, respectively). Thank you.

   Best regards,
   Oscar

   Posted 8 years ago #

2. Bradley Treeby
   

   Developer
   Joined: Mar '10
   Posts: 1,643

   Hi Oscar,

   I've not tried it using a velocity source, but it's certainly possible if you record and re-transmit the pressure using a Dirichlet boundary condition. I created a simple example in 1D to demonstrate. You can apply the same idea in 2D and 3D.

   Note, re-transmitting the pressure over a single plane is not sufficient to perfectly recover the spatial gradients, so your output might have a small phase offset compared to a simulation without the split. One of the forum users discovered that this could be overcome if you record and transmit over a plane with a thickness greater than one grid point (see this post and the related paper).

   You can test this in the example above by increasing the value of overlap.

   Hope that helps,

   Brad.

   Posted 8 years ago #

RSS feed for this topic

Reply

You must log in to post.