k-Wave

1. laughingrice
   

   Member
   Joined: May '13
   Posts: 5

   I've been trying to match the k-wave simulation to a standard finite difference acoustic simulation with a Ricker wavelet forcing function.

   i.e solve u_tt - c^2 u_xx = f

   As Ricker is the second order time derivative of a Gaussian, the source pressure is a Gaussian. I tried using that for source.p which did not produce the same results. Looking at the governing equations defined in eq 2.8 in the user manual, what is defined as the pressure source term matches the mass source in that equation (d\rho/dt), and using that (the normalized time derivative of the Gaussian pressure source) does match the result of the finite difference simulations.

   I wanted to verify what does source.p actually define, a pressure source term as in the function help file and section 3.4 of the manual or a mass source term (flow of mass) as defined in section 2.2 of the manual.

   Thanks

   Posted 5 years ago #

2. Bradley Treeby
   

   Developer
   Joined: Mar '10
   Posts: 1,643

   That's correct - the mass source term in k-Wave is equivalent to a source term with a time derivative in the second order wave equation (see Eq. 2.9 in the manual V1.1).

   The source term in the code appears as in Eq. 2.17 (as a mass source). The user inputs are given in units of pressure, so it's scaled appropriately.

   Posted 5 years ago #

RSS feed for this topic

Reply

You must log in to post.