http://www.k-wave.org/forum/topic/transition-from-sound-in-liquids-to-solids Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 09:43:06 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/transition-from-sound-in-liquids-to-solids#post-5350 Fri, 18 Dec 2015 01:02:20 +0000 mavheart 5350@http://www.k-wave.org/forum/ <p>Hi Brad,</p> <p>thank you for your help :)</p> <p>Marvin </p> http://www.k-wave.org/forum/topic/transition-from-sound-in-liquids-to-solids#post-5339 Tue, 01 Dec 2015 11:31:52 +0000 Bradley Treeby 5339@http://www.k-wave.org/forum/ <p>Hi Marvin,</p> <p>It seems to be a plot scale issue. With auto plot scaling, the limits are constantly adjusted, so it looks like there are waves remnant in your grid after the pulse reaches the PML, even though these are very small.</p> <p>If you change your code to <code>&#39;PlotScale, [-1, 1]</code> it should work as expected. Also, I would advise using a slightly fatter PML (perhaps 20 grid points).</p> <p>Hope that helps,</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/transition-from-sound-in-liquids-to-solids#post-5326 Tue, 24 Nov 2015 22:19:57 +0000 mavheart 5326@http://www.k-wave.org/forum/ <p>Hi Brad,</p> <p>the new version does not seem to solve the issue. I would be glad if you took a look at my code :)</p> <p>Best regards,<br /> Marvin</p> <p>clear all;</p> <p>DATA_CAST = 'double';<br /> MASK_PLANE = 'xy';<br /> USE_STATISTICS = true;</p> <p>PML_X_SIZE = 10;<br /> PML_Y_SIZE = 10;</p> <p>Nx = 100;<br /> Ny = 100;</p> <p>c_water = 1500; % [m/s]<br /> c_plastic = 2230; % [m/s]<br /> c_poly = 1560; % [m/s]<br /> c_air = 343; % [m/s] @ 20°C</p> <p>rho_water = 998; % [kg/m^3]<br /> rho_plastic = 1330; % [kg/m^3]<br /> rho_poly = 1100; % [kg/m^3]<br /> rho_air = 1.22; % [kg/m^3]</p> <p>freq = 15e3; % [Hz]</p> <p>lambda = c_water/freq;</p> <p>x = 1;</p> <p>dx = x/Nx; % [m]<br /> dy = dx; % [m]</p> <p>kgrid = makeGrid(Nx, dx, Ny, dy);</p> <p>tr_pos = [round(0.5*Nx), round(0.5*Ny)];<br /> tr_rad = 10.8e-2/(2);</p> <p>tr_pos_abs = [tr_pos(1)*dx, tr_pos(2)*dy]</p> <p>ITC = makeDisc(Nx, Ny, tr_pos(1), tr_pos(2), round(tr_rad/dx));</p> <p>medium.sound_speed_compression = c_water*ones(Nx, Ny);<br /> medium.sound_speed_shear = zeros(Nx, Ny);<br /> medium.density = rho_water*ones(Nx, Ny);</p> <p>t_end = 10e-4; % [s]</p> <p>kgrid.t_array = makeTime(kgrid, c_water*1.5);</p> <p>delta_t = kgrid.t_array(2) - kgrid.t_array(1);</p> <p>ncycles = 2;</p> <p>input_signal = toneBurst(1/delta_t, freq, ncycles);</p> <p>source.sxx = input_signal;<br /> source.syy = input_signal;</p> <p>source.s_mask = ITC;</p> <p>sensor.mask = zeros(Nx, Ny);<br /> sensor.mask(50, 100) = 1;</p> <p>input_args = {'PlotScale', 'auto','PMLSize', [PML_X_SIZE, PML_Y_SIZE],...<br /> 'DataCast', DATA_CAST, 'DataRecast', true,'PlotPML', true};</p> <p>sensor_data = pstdElastic2D(kgrid, medium, source, sensor, input_args{:}); </p> http://www.k-wave.org/forum/topic/transition-from-sound-in-liquids-to-solids#post-5321 Tue, 24 Nov 2015 15:20:50 +0000 Bradley Treeby 5321@http://www.k-wave.org/forum/ <p>Hi Marvin,</p> <p>That certainly sounds suspicious. If you post your code, I can take a look.</p> <p>Also, have you updated to k-Wave version 1.1.1? This corrected a bug in the <code>pstdElastic</code>codes that might affect stability. </p> <p>Brad. </p> http://www.k-wave.org/forum/topic/transition-from-sound-in-liquids-to-solids#post-5314 Tue, 17 Nov 2015 23:26:33 +0000 mavheart 5314@http://www.k-wave.org/forum/ <p>Hi Brad!</p> <p>Thank you for your reply. I used a test simulation in 2D to have a look at the behaviour of pstdElastic for liquids before applying it to more complex volumes and I encountered a problem:</p> <p>I put a stress source (disc) into a grid of water. Then I created a burst signal and applied it to the source like this:</p> <p>source.sxx = input_signal;<br /> source.syy = input_signal;</p> <p>It looked like it would work fine - I was trying to achieve homogenous radiation - but as the pressure waves reached the corners of my grid, everything inside the medium but the PML seemed to explode, reaching high pressure values. Is there a way to mess up your PML so that this happens? I tried turning it off and it was a completely different result.</p> <p>Thanks for your help!</p> <p>Marvin </p> http://www.k-wave.org/forum/topic/transition-from-sound-in-liquids-to-solids#post-5289 Fri, 23 Oct 2015 10:04:11 +0000 Bradley Treeby 5289@http://www.k-wave.org/forum/ <p>Hi Marvin,</p> <p>To do this, you can use the <code>pstdElastic</code> functions and set the shear properties in the fluid part of the domain to zero. There is an example of this included in the toolbox (see <a href="http://www.k-wave.org/documentation/example_ewp_layered_medium.php">here</a>).</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/transition-from-sound-in-liquids-to-solids#post-5283 Thu, 22 Oct 2015 20:43:55 +0000 mavheart 5283@http://www.k-wave.org/forum/ <p>Hello!</p> <p>I'm interested in sound propagation from liquid to solid and vice versa.<br /> I understand the new functions pstdElastic are meant for solids while kSpaceFirstOrder is meant for liquids. Is there a way built in to have those simulations interact in a way that allows for pressure or stress waves to travel from one to the other?</p> <p>Thank you very much for this powerful toolbox!</p> <p>Marvin Willam </p>