k-Wave

1. babakkh
   

   Member
   Joined: Feb '17
   Posts: 10

   Hi,

   I need to define absorption for part of my simulation domain.
   I am using medium.alpha_power the same for all domain and the medium.alpha_coeff is non-zero for the part of domain that must be absorbing.
   It works for high frequency sinusoidal sources but when the frequency of the source is low or it is broadband non-sinusoidal the absorption doesn't work. it seems there is leakage of acoustic energy to the absorbing domain. Would you please help me if k-wave can absorb non-harmonic or low-frequency excitation?

   Thanks,

   Babak

   Posted 6 years ago #

2. Bradley Treeby
   

   Developer
   Joined: Mar '10
   Posts: 1,643

   Hi Babak,

   Yes, the absorption also works for broadband sources. Foe example, see the "Modelling Power Law Absorption Example". I'm not sure I follow what you mean by "leakage of acoustic energy to the absorbing domain".

   Brad.

   Posted 6 years ago #

3. babakkh
   

   Member
   Joined: Feb '17
   Posts: 10

   Hi Brad,

   Thanks for your response. I couldn't find the example Power Law Absorption. Would you please tell me where that file is?

   Regards,

   Babak

   Posted 6 years ago #

4. Bradley Treeby
   

   Developer
   Joined: Mar '10
   Posts: 1,643

   Hi Babak,

   All the k-Wave examples are accessible via the MATLAB help browser. Open help, then select k-Wave. You should find the example under the heading "Numerical Analysis".

   Brad.

   Posted 6 years ago #

5. babakkh
   

   Member
   Joined: Feb '17
   Posts: 10

   Hi again,

   I have tried several times to model absorption in the part of the domain for the case of transient pressure source. But it becomes unstable. The same model works for a harmonic source.

   I would highly appreciate it if you could have an example of this scenario to see if it works.

   Regards,

   Babak

   Posted 6 years ago #

6. Bradley Treeby
   

   Developer
   Joined: Mar '10
   Posts: 1,643

   Hi Babak,

   Could you post a simple example of your code working for one case, but not the other? I'll take a look.

   Brad.

   Posted 6 years ago #

RSS feed for this topic

Reply

You must log in to post.