Hi brad and guys,
Recently, I have been working on the function kspaceFirstOrder3D, which has a calculation method for nonlinear attenuation sound pressure.But I'm a little confused, it doesn't seem to include the Laplace operator.I want to confirm this formula for calculating the sound pressure in the program:
p = c0.^2 .* (...(rhox + rhoy + rhoz) ...+ absorb_tau .* real(ifftn( absorb_nabla1 .* fftn(rho0 .* (duxdx + duydy + duzdz)) ))... absorb_eta .* real(ifftn( absorb_nabla2 .* fftn(rhox + rhoy + rhoz) ))...+ medium.BonA .*(rhox + rhoy + rhoz).^2 ./ (2 .* rho0) ...);
Is the formula the same as eq(2.7) in Kwave_user_manual_1.1?Does it lose the term Lρ?It would be even better if you could explain every item in the program.In fact, I've been wondering how accurately Kwave calculates the sound pressure, and how the sound pressure decays so quickly with distance, even if the attenuation is not set in water,it may be that I don't understand deeply enough. I hope you can help me.Thank you.
kWave
A MATLAB toolbox for the timedomain
simulation of acoustic wave fields
Problems in nonlinear simulation
(4 posts) (2 voices)
Posted 3 years ago #

Posted 3 years ago #

Hi brad,
Thank you for your recommended articles,maybe I still do not fully understand, please forgive my ignorance.In kWave,the coupled governing equations can also be combined to give a modified form of the Westervelt equation valid for heterogeneous media with power law absorption?The penultimate formula in III B?Or is it the first time that I asked the above formula to calculate the sound pressureIn addition, why does the sound pressure decay so quickly in both two  and threedimensional simulations?Even if it decays at 1/ distance.This gives me an incorrect sound pressure gain.Is there any way to prevent the attenuation of my sound pressure source signal from being so serious? Now basically my source signal will attenuate by 10,000 times after passing 2cm (for example, from 1MPa to 100Pa), but I want to simulate the nonlinear sound wave focusing, so I must ensure that the sound pressure signals at the focus position are large enough.Thank you for your patience.
jackYANGPosted 3 years ago # 
Hi Jack,
The waves will decay due to both geometric spreading and acoustic absorption. The former will depend on the shape of your source (e.g., point source vs large piston source), and the latter on the medium properties and acoustic frequency.
One point to check in your code if you think something is wrong: make sure the source and sensor points are not within the perfectly matched layer.
Brad.
Posted 3 years ago #
Reply
You must log in to post.