Hi
Thank you for your great toolbox.
I was reading about the implementation of absorption modell which is implemented as a integrodifferential operator in the equation of state, and I think there was a very little error. As stated in the manual "To avoid needing to explicitly calculate the time derivative of the acoustic density (which would require storing a copy of at least ρ n and ρ n−1 in memory), the temporal derivative
of the acoustic density is replaced using the linearized mass conservation equation", but if we derive the equation from Navier equation for viscous fluid (also considering thermal conductivity which lead to Kuznetsov equation) the density should be outside of the Fourier transform.
i think this error happened because, equation of state supposed to relate density and pressure and we tried to replace time derivative of density with divergence of flow.
kWave
A MATLAB toolbox for the timedomain
simulation of acoustic wave fields
Acoustic absorption modell for inhomogeneous medium
(3 posts) (2 voices)
Posted 2 months ago #

Hi 1369mohammadi,
Thanks for your comment. I'm not sure I follow what you mean by outside the Fourier transform? If you add a classical viscous absorption term to the momentum equation, I agree that there will be no time derivative. However, if you add a classical thermal loss term to the equation of state, you will still have a time derivative of the density. When you solve this numerically, this will need to be discretised somehow (e.g., as the manual mentions). In some sense, the fractional loss term in kWave replaces the thermal loss term.
Brad.
Posted 1 month ago # 
Thank you for your response.
I'm sorry for my bad explanations. I was concerned about inhomogeneous medium (not the fractional loss). I think rho_0 in the first term in RHS of equation 2.24 (in the manual) should be outside of the spatial Fourier transform. If the medium were homogeneous it would be a constant and won't matter if we place it inside or outside of spatial Fourier transform, but for a medium with inhomogeneous density it would be different. A simply viscous medium was a supporting example for my logic (a counterexample for current formula).Posted 3 weeks ago #
Reply
You must log in to post.