k-Wave User Forum » User Favorites: Guillermo Roa

k-Wave User Forum » User Favorites: Guillermo Roa guillermo-roa Support for the k-Wave MATLAB toolbox en-US Tue, 12 May 2026 23:08:59 +0000 http://bbpress.org/?v=1.0.2 <![CDATA[Search]]> q http://www.k-wave.org/forum/search.php biaojiang on "Nondestructive Testing (NDT) using ultrasound propagation" http://www.k-wave.org/forum/topic/nondestructive-testing-ndt-using-ultrasound-propagation#post-6313 Tue, 06 Mar 2018 16:14:20 +0000 biaojiang 6313@http://www.k-wave.org/forum/ Hello, all, I tried to simulate NDT also. When the impedence contrast is too large, like Aluminum or steel, even putting it in water (immersion mode), and set CFL=0.1 (default is 0.3), the k-Wave results are still not as expected. So I can only reduce the material sound speed and density so that they are not too bigger than those of water, then the results are OK. Best Biao RAHUL2011 on "Nondestructive Testing (NDT) using ultrasound propagation" http://www.k-wave.org/forum/topic/nondestructive-testing-ndt-using-ultrasound-propagation#post-125 Thu, 24 Mar 2011 09:37:30 +0000 RAHUL2011 125@http://www.k-wave.org/forum/ hello friends, i am trying to simulate ultrasound in aluminum using focused probe. i wanted to know that is it possible to include elastic properties of material using k-wave toolbox. because if it is not included shear wave will not be visualized .... plzzz help Guillermo Roa on "Nondestructive Testing (NDT) using ultrasound propagation" http://www.k-wave.org/forum/topic/nondestructive-testing-ndt-using-ultrasound-propagation#post-95 Sat, 29 Jan 2011 19:48:50 +0000 Guillermo Roa 95@http://www.k-wave.org/forum/ Hi again; I ran my example successfully, i'm still having problems posting images, so here are the url's: 1) The model: The model was a cross-section of a concrete column measuring 20 cm on each side, 8 reinforcements of steel and 1 inch of diameter were disposed, a PML of 20 pixels were defined all around the column: link: <a href="http://img163.imageshack.us/i/modelod.jpg/" rel="nofollow">http://img163.imageshack.us/i/modelod.jpg/</a> 2) The Source: because in the faces x and z and along them the ultrasound source is going to be fired, i considered two sources and i set them just in the border of the column, one on the x face and the other on the Z face; The ultrasound has a freq of 2 MHz and magnitude 0.4 3) The sensors: Two sensors were disposed of the same size and oposite to the two sources. 4) Simulation results: Those are the images resulting from the reconstruction using time reversal link 2: <a href="http://img218.imageshack.us/i/resultado1d.jpg/" rel="nofollow">http://img218.imageshack.us/i/resultado1d.jpg/</a> link 3: <a href="http://img89.imageshack.us/i/resultado2.jpg/" rel="nofollow">http://img89.imageshack.us/i/resultado2.jpg/</a> 5) Image processing: Finally i filtered the image from time reversal in order to identify the reinforcements link 4: <a href="http://img211.imageshack.us/i/filtroident.jpg/" rel="nofollow">http://img211.imageshack.us/i/filtroident.jpg/</a> so i'll wait further comments about my example; Thanks; Kind regards Guillermo Roa on "Nondestructive Testing (NDT) using ultrasound propagation" http://www.k-wave.org/forum/topic/nondestructive-testing-ndt-using-ultrasound-propagation#post-93 Wed, 05 Jan 2011 01:57:14 +0000 Guillermo Roa 93@http://www.k-wave.org/forum/ Hi; In the last code that i posted i improved a few things: * dx and dz values were defined as (speed_of_sound/freq)/2 * i set the mask sensor just in the border of the column, with 3 sensors as you can see in the step 3 of my preview post * the source of the ultrasound pulse were set just in the border of the column With those changes the simulation ran with better results. Now, i have a few questions: 1) It is possible to make an image reconstruction of the medium with the data acquired by the sensors in the way that those were set in the code? 2) It is possible to use the algorithm of the toolbox (either fft or time-reversal) in order to obtain an approximation image for the medium? I appreciate your help, thanks; Guillermo Bradley Treeby on "Nondestructive Testing (NDT) using ultrasound propagation" http://www.k-wave.org/forum/topic/nondestructive-testing-ndt-using-ultrasound-propagation#post-92 Tue, 04 Jan 2011 03:49:55 +0000 Bradley Treeby 92@http://www.k-wave.org/forum/ Hi Guillermo, did you manage to run your example successfully or are you still having problems? Guillermo Roa on "Nondestructive Testing (NDT) using ultrasound propagation" http://www.k-wave.org/forum/topic/nondestructive-testing-ndt-using-ultrasound-propagation#post-89 Sat, 01 Jan 2011 20:09:02 +0000 Guillermo Roa 89@http://www.k-wave.org/forum/ thanks for fix the images. I apreciate for your comments about the code and your suggestions for my questions. kind regards; Guillermo Guillermo Roa on "Nondestructive Testing (NDT) using ultrasound propagation" http://www.k-wave.org/forum/topic/nondestructive-testing-ndt-using-ultrasound-propagation#post-85 Wed, 08 Dec 2010 01:56:48 +0000 Guillermo Roa 85@http://www.k-wave.org/forum/ Sorry about the images; those are the url's following the order of my post: 1)http://img255.imageshack.us/i/step1q0.jpg/ 2)http://img37.imageshack.us/i/step2u.jpg/ 3)http://img52.imageshack.us/i/step3pn.jpg/ 4)http://img17.imageshack.us/i/step4v0.jpg/ 5)http://img525.imageshack.us/i/step5y.jpg/ by the way, i've tried to found you on facebook but i don't see the picture that you show on your page... [ADMIN EDIT: I've fixed your links above, for imageshack you need to use the direct link to the image to embed them in forum] bencox on "Nondestructive Testing (NDT) using ultrasound propagation" http://www.k-wave.org/forum/topic/nondestructive-testing-ndt-using-ultrasound-propagation#post-84 Wed, 08 Dec 2010 01:45:28 +0000 bencox 84@http://www.k-wave.org/forum/ On your point (1), the dx and dz are the pixel dimensions measured in metres, but they must also be smaller than half the shortest wavelength. Your pictures don't seem to have appeared on the forum. Can you try and put them on again? Thanks, Ben Guillermo Roa on "Nondestructive Testing (NDT) using ultrasound propagation" http://www.k-wave.org/forum/topic/nondestructive-testing-ndt-using-ultrasound-propagation#post-83 Sun, 05 Dec 2010 19:14:01 +0000 Guillermo Roa 83@http://www.k-wave.org/forum/ Thanks for the answer; respect your points: 1) i took note for the data dx and dz; in the help, it says that dx and dz are the dimensions of each pixel of the matrix given in Meters. But now with your help i know that it corresponds to the wavelenght, changing this value the simulation improves its result. 2) You confirm my fear, the contrast between the two mediums is quite large, but i can not change the enviroment medium (air). 3) setting the correct values of dx and dz, the source and waves are now visibles. 4) hehehe excuse me, but english is not my mother language, i used the expression 'crash into the section' trying to explain that it seems that the column resonates. Now, i made a few changes in the model to simulate: a) i don't need that the waves that comes into the column can go out. I can define a square mask sensor just in the border of the structure and make the measurement of the waves at this point. So i made another code lines and the results of the simulation are quite better; I send you the images step by step about what i´m doing (i feels i'm talking with a strong broken english and images are better that 1000 words hehehe) also the new code: step 1: MODEL IDEA: I set only three lines for the square sensor mask, in order to avoid measuring the initial pressure. <img src="http://img255.imageshack.us/img255/8146/step1q0.jpg" /> Step 2: MODEL IN CODE LINES: <img src="http://img37.imageshack.us/img37/2272/step2u.jpg" /> Step 3: SIMULATE <img src="http://img52.imageshack.us/img52/5687/step3pn.jpg" /> Step 4: WHAT I CAN SEE (INTERPRETATE THE WAVES) <img src="http://img525.imageshack.us/img525/9813/step5y.jpg" /> Step 5: (Question 1) MESUREMENT OF THE SIGNAL WITH SENSORS How can i interpretate this graphic? <img src="http://img17.imageshack.us/img17/3840/step4v0.jpg" /> Finally i send you the code, i appreciate further comments and suggestions, thank you. %% hacer la grilla (make the grid) Nx=200; % pixeles en x (columna) Nz=200; % pixeles en z (fila) dx=0.1715e-3; %from 343/2e6 dz=0.1715e-3; kgrid = makeGrid(Nx, dx, Nz, dz); %% modelo a simular (model of mediums) % % % % % % % % % % % % % % % % % % ' - - - - - - - - - - - - - ' % % ' 0 0 0 0 0 0 0 0 0 0 0 0 0 ' % % ' 0 0 0 0 0 0 0 0 0 0 0 0 0 ' % % ' 0 0 + + 0 0 0 0 0 + + 0 0 ' % % ' 0 0 + + 0 0 0 0 0 + + 0 0 ' % % ' 0 0 0 0 0 0 0 0 0 0 0 0 0 ' % % ' 0 0 0 0 0 0 0 0 0 0 0 0 0 ' % % ' 0 0 0 0 0 0 0 0 0 0 0 0 0 ' % % ' 0 0 0 0 0 0 0 0 0 0 0 0 0 ' % % ' 0 0 0 0 0 0 0 0 0 0 0 0 0 ' % % ' 0 0 + + 0 0 0 0 0 + + 0 0 ' % % ' 0 0 + + 0 0 0 0 0 + + 0 0 ' % % ' 0 0 0 0 0 0 0 0 0 0 0 0 0 ' % % ' 0 0 0 0 0 0 0 0 0 0 0 0 0 ' % % % % % % % % % S % % % % % % % % % % % % % % % % % % % % % % % % % % %=air (density: 1.18 Kg/m3) (speed of sound 343 m/s) % 0=concrete (density 2800 kg/m3) (speed of sound 4000 m/s) % +=steel (density 7850 kg/m3) (speed of sound 5100 m/s) % ' sensors % S source of sound %% construct the medium %% speed of sound (air %) medium.sound_speed = 343*ones(Nz, Nx); % [m/s] %% speed of sound (concrete 0) medium.sound_speed(Nz/8:(7*Nz/8),Nz/8:(7*Nz/8)) = 4000; % [m/s] %% speed of sound (4 reinforcements, steel +) medium.sound_speed((2*Nz/8):(3*Nz/8),(2*Nz/8):(3*Nz/8)) = 5100; % [m/s] medium.sound_speed((2*Nz/8):(3*Nz/8),(5*Nz/8):(6*Nz/8)) = 5100; % [m/s] medium.sound_speed((5*Nz/8):(6*Nz/8),(2*Nz/8):(3*Nz/8)) = 5100; % [m/s] medium.sound_speed((5*Nz/8):(6*Nz/8),(5*Nz/8):(6*Nz/8)) = 5100; % [m/s] %% density (air %) medium.density = 1.18*ones(Nz, Nx); % [kg/m^3] %% density (concrete 0) medium.density(Nz/8:(7*Nz/8),Nz/8:(7*Nz/8)) = 2800; % [kg/m^3] %% density (4 reinforcements, steel +) medium.density((2*Nz/8):(3*Nz/8),(2*Nz/8):(3*Nz/8)) = 7850; % [kg/m^3] medium.density((2*Nz/8):(3*Nz/8),(5*Nz/8):(6*Nz/8)) = 7850; % [kg/m^3] medium.density((5*Nz/8):(6*Nz/8),(2*Nz/8):(3*Nz/8)) = 7850; % [kg/m^3] medium.density((5*Nz/8):(6*Nz/8),(5*Nz/8):(6*Nz/8)) = 7850; % [kg/m^3] %% set the sensor mask -' sensor.mask = zeros(kgrid.Nz, kgrid.Nx); sensor.mask(Nz/8, Nz/8:(7*Nz/8)) = 1; sensor.mask(Nz/8:(7*Nz/8),Nz/8) = 1; sensor.mask(Nz/8:(7*Nz/8),(7*Nz/8)) = 1; %% source of ultrasound, S [kgrid.t_array dt] = makeTime(kgrid, medium.sound_speed); % fuente de ultrasonido (source) source.p_mask = zeros(Nz, Nx); source.p_mask(176,100) = 5; % señal (signal) source_freq = 2e6; source_mag = 4; source.p = source_mag*sin(2*pi*source_freq*kgrid.t_array); % suavizar fuente (smooth) source.p = filterTimeSeries(kgrid, medium, source.p); %% simular (run simulation) [sensor_data p_final] = kspaceFirstOrder2D(kgrid, medium, source, sensor,'PlotLayout', true); %% visualizar (visualization) figure; imagesc(sensor_data, [-1, 1]); colormap(getColorMap); ylabel('Sensor Position'); xlabel('Time Step'); colorbar; [ADMIN EDIT: fixed image links] bencox on "Nondestructive Testing (NDT) using ultrasound propagation" http://www.k-wave.org/forum/topic/nondestructive-testing-ndt-using-ultrasound-propagation#post-82 Fri, 03 Dec 2010 11:28:08 +0000 bencox 82@http://www.k-wave.org/forum/ Hi Guillermo, There are a few separate issues here. 1) A frequency of 2MHz has a wavelength in air of 343/2e6 = 0.17 mm so for the grid to support this frequency dx and dz must be half this size at most (there must be at least 2 points per wavelength). 2) There is a very large contrast between air and concrete. Fourier-based methods such as k-Wave are very good when everything is smooth and the contrasts are not too big - the code was originally designed for modelling ultrasound in biological tissue where the contrasts are low - but for a large contrast they can accumulate oscillatory errors due to the fact that steep edges cannot be represented very well by a small number of Fourier components. Smoothing the density and sound_speed using the smooth function (which happens automatically in version 0.2) and keeping the timestep small by setting a small value for CFL in makeTime will help reduce this, but not eliminate when the contrast is this severe. Unfortunately your example is almost the worse case possible for k-Wave. Could you put your column in water? That should work better (although the contrast is still quite large). 3) The default scaling is fixed in version 0.2, so adjust your source magnitude so you can see the wave, or use the 'PlotScale' option in kspaceFirstOrder2D. The amplitude of the wave in the air and in the concrete will be quite different because of their different material properties. The pressure (which is what k-Wave plots) will be higher in the concrete because the density is so much greater. 4) I'm not sure what you mean by 'crash into the section' but any waves that do penetrate the concrete will become trapped in there because the transmission coefficient from concrete to air will be so low (less than 1% transmission I would guess, but you can calculate it simply enough). Thank you for using k-Wave and I'm sorry that it perhaps isn't the best tool for your particular example. Kind regards, Ben Guillermo Roa on "Nondestructive Testing (NDT) using ultrasound propagation" http://www.k-wave.org/forum/topic/nondestructive-testing-ndt-using-ultrasound-propagation#post-81 Thu, 02 Dec 2010 21:51:16 +0000 Guillermo Roa 81@http://www.k-wave.org/forum/ Good day; I'm Ingeneering student and nowadays I'm working on my Bachelor's thesis. I need to do a simulation of the ultrasound propagation across a square section (2D) that represents a column of concrete with its reinforcements in order to generate an image reconstruction of the column based on the ultrasound propagation. I began to use the K-WAVE toolbox, the first step, like the examples, is to "construct" the k-grid and to define the medium properties (speed of sound and density). The medium in my case is composed by: air (that represents the enviroment and where the source of the ultrasound is going to be set), concrete (square section of the column) and steel (square or circular sections that are allowed into the square section of the column representing the reinforcements) giving a model for a Heterogeneous Propagation Medium. As i said, i started defining the medium, i set a simple model avoiding the reinforcements, only the air and in its middle a concrete section were defined with the following data: concrete density: 2730(kg/m3) concrete (speed of sound): 4000(m/s) air density: 1.18(kg/m3) air (speed of sound): 343(m/s) The medium in my case is composed by: 1.air (that represents the enviroment), 2.concrete (square section of the column) and 3.steel (square or circular sections that represents the reinforcements are allowed into the square section of the column). As i said, i started setting the properties of the medium, in this first attemp, i set only air as enviroment and concrete with stone as the column, avoiding the reinforcement. Then i set the ultrasound source and the sensors, then i started the simulation. So, here´s my problem, when the simulation starts, the source of ultrasound doesn´t appear, and in other cases, the waves of sound enter into the square section but they don´t go out, instead, they began to crash into the section. I apreciate that you help me with this problem or telling me what is wrong. I send you the code, thank you!! '%% crear k grid Nx=200; Nz=200; dx=1; dz=1; kgrid = makeGrid(Nx, dx, Nz, dz); %% grid of air (enviroment) (speed of sound) for i=1:Nz for j=1:Nx-150 medium.sound_speed(i,j) = 343; end end for i=1:Nz-150 for j=Nx-150:Nx medium.sound_speed(i,j) = 343; end end for i=Nz-49:Nz for j=Nx-150:Nx medium.sound_speed(i,j) = 343; end end for i=Nz-150:Nz for j=Nx-49:Nx medium.sound_speed(i,j) = 343; end end %% grid of the square section (concrete (4000) and stone (3700)) (speed of sound) for i=Nz-149:Nz-50 for j=Nx-149:Nx-50 if rem((i+j),2)~=0 medium.sound_speed(i,j) = 4000; end end end for i=Nz-149:Nz-50 for j=Nx-149:Nx-50 if rem((i+j),2)==0 medium.sound_speed(i,j) = 3700; end end end %% grid of air (enviroment) (density) for i=1:Nz for j=1:Nx-150 medium.density(i,j) = 1.18; end end for i=1:Nz-150 for j=Nx-150:Nx medium.density(i,j) = 1.18; end end for i=Nz-49:Nz for j=Nx-150:Nx medium.density(i,j) = 1.18; end end for i=Nz-150:Nz for j=Nx-49:Nx medium.density(i,j) = 1.18; end end %% grid of the square section (concrete (2800) and stone (2500)) (density) for i=Nz-149:Nz-50 for j=Nx-149:Nx-50 if rem((i+j),2)~=0 medium.density(i,j) = 2800; end end end for i=Nz-149:Nz-50 for j=Nx-149:Nx-50 if rem((i+j),2)==0 medium.density(i,j) = 2500; end end end %% sensors sensor_radius = 90; num_sensor_points = 500; sensor.mask = makeCartCircle(sensor_radius, num_sensor_points); %% ultrasonido [kgrid.t_array dt] = makeTime(kgrid, medium.sound_speed); % source of ultrasound source.p_mask = zeros(Nz, Nx); source.p_mask(152, 100) = 1; % signal source_freq = 2e6; source_mag = 250; source.p = source_mag*sin(2*pi*source_freq*kgrid.t_array); % smooth the source %source.p = filterTimeSeries(kgrid, medium, source.p); % run the simulation [sensor_data p_final] = kspaceFirstOrder2D(kgrid, medium, source, sensor,'PlotLayout', true); figure; imagesc(sensor_data, [-1, 1]); colormap(getColorMap); ylabel('Sensor Position'); xlabel('Time Step'); colorbar;'