Hi. As a simple example I model a 2d strip radiating a wavelet pulse into water, with a steel reflector wall at the other end of the grid. I have to set the time scale using the highest velocity in the model, namely steel. This makes things run slowly, and I get numerical dispersion by the time the pulse gets back to the radiator. Is there a way to set different grid spacings or time steps in the two materials, so that the optimal conditions for the modelor exist in each medium?
A MATLAB toolbox for the time-domain
simulation of acoustic wave fields
Variable grid spacing?(3 posts) (3 voices)
If your problem is with numerical dispersion, you might want to consider setting
medium.sound_speed_refto the sound speed in water. While the simulation will no longer be unconditionally stable (you will need to satisfy the CFL condition given in the k-Wave manual), the propagation in the water will be dispersion free.
Brad.Posted 1 day ago #
You must log in to post.