【声辐射】——不同坐标系下的格林函数

概述

本文整理了不同坐标系下Helmholtz方程对应的格林函数
$$({\nabla }^2+k^2)G(x)=\delta (x)$$

直角坐标系(描述一维无限区域)

物理背景：在$x=x_0$位置处放置一个平行于YOZ平面的无限大面声源，$G(x,x_0)$为其产生的声场，声压不随$y$和$z$方向的坐标变化，波阵面垂直于$x$轴，即沿$x$方向传播的平面波。
$$\frac{d^{2}G(x,x_0)}{d{x}^2}+k^{2}G(x,x_0)=\delta (x-x_0)$$
$$G(x,x_0)=\frac{j}{2k}{e}^{\pm jk|x-x_0|}$$
注：正负号代表传播方向，具体方向要看表达式中$\omega t$前面的正负号。

• 平面波按柱面波展开
  $$e^{jkx}=J_0(k\rho )+2\sum _{m=0}^{\infty }{j}^m{J}_m(k\rho )cos(m\phi )$$
  其中，$J_m(k\rho )$是柱贝塞尔函数，物理意义是柱面波的驻波声场(可以展开成柱汉克尔函数，物理意义为柱面波的行波声场)；而$x=\rho cos(\phi )$，$\rho$对应径长，$\phi$对应极角，将直角坐标系转换成了柱坐标系。

• 平面波按球面波展开
  $$e^{ikz}=e^{ikrcos\theta }=\sum _{l=0}^{\infty }{i}^l(2l+1)j_l(kr)P_l(cos\theta )$$
  其中，$j_l(kr)=\sqrt{\frac{\pi }{2kr}}{J}_{l+\frac{1}{2}}(kr)$是球贝塞尔函数，物理意义是球面波的驻波声场。

柱(极)坐标系(描述二维无限区域)

物理背景：在原点放置一个与$z$轴平行的无限长的线声源，声场与极角$\phi$无关，$G(k\rho )$为其产生的柱面波。
$$\frac{1}{\rho }\frac{d}{d\rho }(\rho \frac{dG(k\rho )}{d\rho })+k^{2}G(k\rho )=-\frac{\delta (\rho )}{2\pi \rho }$$
$$\frac{d^{2}G}{d{\rho }^2}+\frac{1}{\rho }\frac{dG}{d\rho }+k^{2}G=0,\rho >0$$
上述方程的解是0阶贝塞尔函数，为了描述行波，故用0阶汉克尔函数去描述解的通式。根据信号与系统柜求解冲激响应的方法，可以知道通过$G^{{}^{\prime }}(0_{+})=1$可以求出通式前面的的系数，即：
$$G(k\rho )=\frac{i}{4}{H}_0^{(1)}(k\rho )$$

• 柱面波按平面波进行展开
  $$H_0^{(1)}(k\rho )=\frac{1}{\pi }{\int }_{-\infty }^{\infty }\frac{e^{i(k_{x}x+k_{y}y)}}{k_y}\mathrm{d}{\mathrm{k}}_{\mathrm{y}}$$
  即角谱法，可以通过傅里叶变换证明，也可以通过求得二维无限区域的Helmholtz方程在直角坐标系下的解等价之来证明。

球坐标系(描述三维无限区域)

物理背景：在原点放置一个点声源，声场与$\phi ,\theta$无关，$G(kr)$为球面波。
$$\frac{1}{r^2}\frac{d}{dr}(r^2\frac{dG(kr)}{dr})+k^{2}G(k\rho )=-\frac{\delta (r)}{4\pi r^2}$$
$$G(kr)=\frac{e^{jkr}}{4\pi r}$$
实际上，三维无限空间下的格林函数$G(kr)$也可以用球汉克尔函数$h_0^{(1)}(kr)$去表示，只是可以化简：
$$h_0^{(1)}(kr)=\frac{e^{jkr}}{jr}$$

• 球面波按平面波展开
  $$\frac{e^{jkr}}{r}=\iint \frac{j}{2\pi k_z}{e}^{j(k_{x}x+k_{y}y\pm k_{z}z)}\mathrm{d}{\mathrm{k}}_{\mathrm{x}}\mathrm{d}{\mathrm{k}}_{\mathrm{y}}$$
  角谱法，没什么好解释的。
• 球面波按柱面波展开
  $$\frac{e^{jkr}}{r}=i{\int }_0^{\infty }{J}_0(k_{\rho }\rho )\frac{e^{i{k}_z|z|}}{k_z}{k}_{\rho }\mathrm{d}{\mathrm{k}}_{\rho }=\frac{i}{2}{\int }_{-\infty }^{\infty }\frac{1}{k_z}{H}_0^{(1)}(k_{\rho }\rho )e^{j(k_{z}z)}{k}_{\rho }\mathrm{d}{\mathrm{k}}_{\rho }$$
  在球面波反射那篇博文中有这个展开式的应用。

Reference：
一本不知来源的数学物理方程书中的格林函数篇