k-Wave

1. anil
   

   Member
   Joined: Mar '16
   Posts: 15

   Hi Bradley,

   I am working on transcranial ultrasound simulations with the "pstdElastic2D" function at 0.5, 0.75, 1, and 1.25 MHz.
   I use 19.5 for shear wave attenuation and 8.83 for compressional wave attenuation.
   I define my source waveform and surface intensity with following part of the code:

   source_freq = 1e6; % [Hz]
   source_surface_intensity = 5; % [W/cm^2]
   source_strength = sqrt(2*source_surface_intensity*1e4*rho0*c0); % [Pa]
   source.sxx = -createCWSignals(kgrid.t_array, source_freq, source_strength, 0);
   source.syy = source.sxx;

   I was wondering if the definition above and the attenuation constants are valid for k-Wave.

   Best,
   Prof. AK

   Posted 3 years ago #

2. anil
   

   Member
   Joined: Mar '16
   Posts: 15

   * Since the source is in the water.
   rho0 = water's density
   c0 = water's compressional wave speed

   Posted 3 years ago #

3. Bradley Treeby
   

   Developer
   Joined: Mar '10
   Posts: 1,640

   Hi Prof. AK,

   There are a variety of different coefficients for the attenuation of bone. This paper uses values quite a bit lower if you are only interested in the absorption, and are modelling the scattering by using a heterogeneous skull.

   The other equations look correct at first glance.

   Hope that helps,

   Brad.

   Posted 3 years ago #

4. anil
   

   Member
   Joined: Mar '16
   Posts: 15

   Bradley,

   Thank you for your detailed explanation. It is really interesting to see this variety of numbers.

   Best,
   Prof. AK

   Posted 3 years ago #

RSS feed for this topic

Reply

You must log in to post.