Hi,
Could you help me figuring out why k-wave crashes when defining high values (e.g. 1e6 Pa) of the source pressure for arbitrary sources such as in the example "example_tvsp_transducer_field_patterns.m"?
Thanks,
Hermes
k-Wave
A MATLAB toolbox for the time-domain
simulation of acoustic wave fields
Error using example_tvsp_transducer_field_patterns.m
(5 posts) (2 voices)-
Posted 9 years ago #
-
Hi Hermes,
What do you mean by crashes? Is it just a plot scaling issue? What happens if you add:
'PlotScale', [-1, 1]*source_magto the list of input options?
Brad.
Posted 9 years ago # -
Hi Dr. Treeby,
Thanks for the prompt answer.
I am simulating the standing wave formation in structures defined by CT images. I was wondering if the artifact seen in the image in the link below is only a scale problem or if it's something else. In the image you can see the formation of waves in the bottom before the wavefront reaches it. My source is similar to the referred example.https://drive.google.com/file/d/0B2glsz_fmxT0Rk81SUZ5ODNMc2M/view?usp=sharing
Thanks,
HermesPosted 9 years ago # -
Without rescaling the example it's also possible to see the wave formation from the sides. Is it a PML setting that need to be changed depending on the source magnitude?
https://drive.google.com/file/d/0B2glsz_fmxT0MV9vendnM0dCM1k/view?usp=sharing
Thanks,
HermesPosted 9 years ago # -
Hi Hermes,
The performance of the PML at normal incidence is a reduction in the pressure amplitude of somewhere between 60 and 80 dB (independent of the initial pressure). This corresponds to three or four decimal places of accuracy. If your source magnitude is of the order 10^6, and your plot scale is set to the default (-1 to 1), then you will undoubtedly see different behaviour, including wave wrapping and the underlying behaviour of the band-limited interpolant.
Brad.
Posted 9 years ago #
Reply
You must log in to post.