I am simulating a 1D problem where the mass density of the fluid changes by several orders of magnitude between the beginning and the end of propagation from 1.2 to 10^-12 over 150 wavelengths. When doing so the code diverges (whether I use linear/nonlinear, attenuating/lossless propagation). When I impose such changes (with several orders of magnitude) to the wave velocity, the code runs without a problem. Did you ever encounter this problem? Any hint on why the code would diverge on big density changes would be useful.
Thanks!
Marcel
k-Wave
A MATLAB toolbox for the time-domain
simulation of acoustic wave fields
code diverges when change of density is too large?
(2 posts) (2 voices)-
Posted 10 years ago #
-
Hi Marcel,
You may find this topic useful:
A Question about The Reflection at The Hetrogeneous Media BoundaryHope it helps :-)
AnthonyPosted 10 years ago #
Reply
You must log in to post.