Hello
I m' developing a simulation of ultrasound for a HIFU application in order to measure the acoustic pressure. However, I have a problem for order of magnitude of found pressure when I do a temporal shoot with a spherical transducer with 256 elements (256 source k wave) .
the grid is :
Nx = 300; % number of grid points in the x direction
Ny = 300; % number of grid points in the y direction
Nz = 400; % number of grid points in the z direction
dx = 0.5e-3; % grid point spacing in the x direction [m]
dy = 0.5e-3; % grid point spacing in the y direction [m]
dz = 0.5e-3; % grid point spacing in the z direction [m]
dt=1e-7;%s
Nt=1200;%s
kgrid.dt=dt;
kgrid.Nt=Nt;
t_end = kgrid.Nt*kgrid.dt;
kgrid.t_array=[0:kgrid.dt:t_end];
the paramater of pressures sources:
power_source=1e6;%Pa
freq=1000e3;%Hz
t_shoot=50e-6;%s
source.p = source_mag * sin(2 * pi * source_freq * t_shoot);
the focusing works but I find pressures order to 10^-23( Pa??)when I use p_max_all for the sensors. I don't understand why?
Math