f0 58[kHz] 驱动频率
c0 343[m/s] 声速
lam0 c0/f0 波长
Dt 2*lam0 换能器直径
Dr 3*lam0 反射镜直径
H 5*lam0/2 高度
dvisc 0.22[mm]*sqrt(100[Hz]/f0) 黏滞边界层厚度
d_p lam0/10 粒径
rho_p 500[kg/m^3] 粒子密度
a0 1.5e6[m/s^2] 换能器的法向加速度