Hello
I'm working on example_tvsp_homogeneous_medium_monopole.m example.
While I increase the mesh grid size ( multiplying nx and ny by 2, so dx lowers to its half) the measured acoustic pressure value also decreases to its half.
Which is unexpected. Can anybody please describe the reason or offer me a way to prevent this situation. How can I compensate this effect?
Thank you
Hi William,
Have you tried setting source.p_mode to 'dirichlet'? That will force the pressure at that point to be what you have input. The default option is 'additive', which adds the source to the current value of the acoustic pressure at the source point. In this latter case, the acoustic pressure amplitude at the source will grow up to a time dx / sound_speed when using a time-varying source.
The reason it grows goes to the fundamentals of k-Wave. When you define a 'point' source, ie. a value at one grid point and zeros at all other grid points, then what k-Wave sees is not a single point source but rather a source function of the form sinc(\pi x / dx) sinc (\pi y / dy) centred at the 'source point'. (It is the band-limited interpolant which matches all the values defined on the grid but it varies in between the grid points. This arises because of the Fourier basis that underlies k-Wave.) It takes a time of dx / sound_speed for the main lobe of this sinc function to propagate away from the source point.
If you are interested in the difference in the source pressure and the pressure at some other point in the field and want to use an additive source type, then one option is to record the pressure at the source point (ie. set one of the sensor points to be the source point).
Hope that helps,
Ben
You must log in to post.