http://www.k-wave.org/forum/topic/instability-modeling-air-and-solids Support for the k-Wave MATLAB toolbox en-US Tue, 12 May 2026 23:37:03 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/instability-modeling-air-and-solids#post-8530 Fri, 13 May 2022 14:51:54 +0000 bencox 8530@http://www.k-wave.org/forum/ <p>Hi bwalker000,</p> <p>Your first plan of making the material properties of air larger, which as you say doesn't change much acoustically, is a trick we've used in the past. That's probably the best option, if you have to model such an interface. If the interface is planar it would be worth thinking about whether you can use an image source to model the reflection. </p> <p>If you want to you can set the time steps yourself using the <code>setTime&#39; method of</code>kWaveGrid', eg.</p> <pre><code>kgrid = kWaveGrid(Nx, dx);&#39;</code></pre> kgrid.setTime(Nt, dt)'</p> <p>Best wishes<br /> Ben </p> http://www.k-wave.org/forum/topic/instability-modeling-air-and-solids#post-8520 Tue, 26 Apr 2022 22:19:43 +0000 bwalker000 8520@http://www.k-wave.org/forum/ <p>I am modeling the reflections from a piece of plastic in air. When I use the following material properties the simulation is unstable and diverges.</p> <pre><code>% source / reflector (PLA) c_plastic = 2730; % [m/s] According to Olympus NDT rho_plastic = 1190; % [kg / m^3] According to ONDA Tables c_air = 345; % [m/s] According to ONDA Tables rho_air = 1.293; % [kg / m^3] According to ONDA Tables</code></pre> <p>However, if I increase the density of the air by a factor of ten, then the simulation is stable:</p> <pre><code>% source / reflector (PLA) c_plastic = 2730; % [m/s] According to Olympus NDT rho_plastic = 1190; % [kg / m^3] According to ONDA Tables c_air = 345; % [m/s] According to ONDA Tables rho_air = 10*1.293; % [kg / m^3] According to ONDA Tables</code></pre> <p>This isn't a bad workaround since the speed of sound stays the same and the reflection coefficient largely stays the same (huge).</p> <p>Alternatively I can reduce the density of the plastic and get a stable solution:</p> <pre><code>c_plastic = 2730; % [m/s] According to Olympus NDT %rho_plastic = 1190; % [kg / m^3] According to ONDA Tables rho_plastic = 1/10*1190; % [kg / m^3] According to ONDA Tables c_air = 345; % [m/s] According to ONDA Tables rho_air = 1.293; % [kg / m^3] According to ONDA Tables</code></pre> <p>My suspicion is that the temporal sampling rate is somehow not properly determined for materials with these extreme mismatches. I would suspect that this can be overridden to resolve the issue. Can anyone comment on this hypothesis or tell me how to do an override of the temporal sampling rate?</p> <p>Thanks. </p>