Hi everyone,
I am trying to simulate the behaviour of a focused spherical ultrasound transducer (64 mm diameter, 100 mm focal length, operating at 1.1 MHz) in water. I am particularly interested in obtaining an axial plot of the pressure and I was using the example_at_focused_bowl_3D.m file to do so, however, it seems that the example does not complete simulating my transducer, as the pressure field map obtained is constant (zero) before reaching the focal position (the simulation is completed, taking almost 2h to complete). Then the axial plot is obtained in which the exact O’Neil solution is perfectly plotted, but the k-wave numerical computation goes to be zero before reaching the focal position. I tried to compare the results with other software’s (like HIFU Beam) but since the k-wave computation is null for most of the region, I couldn’t.
The issue didn’t occur with the original example parameters and I am using these for my simulation:
% medium parameters
c0 = 1500; % sound speed [m/s]
rho0 = 1000; % density [kg/m^3]
% source parameters
source_f0 = 1.1e6; % source frequency [Hz]
source_roc = 100e-3; % bowl radius of curvature [m]
source_diameter = 64e-3; % bowl aperture diameter [m]
source_amp = 0.879e6; % source pressure [Pa]
% grid parameters
axial_size = 120e-3; % total grid size in the axial dimension [m]
lateral_size = 120e-3; % total grid size in the lateral dimension [m]
Does anybody have come across something like this before?
As an alternative I am running the example_diff_focused_ultrasound_heating.m and I was trying to extract the 2D plot of the axial pressure. Any ideas to do so?
Thank you in advance.
All the best,
Tony