计算海洋声学之波动方程求解思路

摘要

本文按照计算海洋声学中常见计算流程对波动方程的求解过程进行了梳理，给出了求解要点。本文侧重于整体思路，对求解过程细节之处并未详细讨论。
#引言

• 计算海洋声学的研究内容：海洋声学问题转变为特定边界条件下波动方程的求解问题，利用数值计算手段予以解决。
• 波动方程的求解思路：通过层层分解，降低波动方程的维度，得到一维波动方程后予以求解，最后通过反变换还原解。
• 需要关注的要点：方程的维度，边界条件，坐标系，源的情况（点源/线源，有源/无源）。
  #求解流程概述
  
  如上图所示，以各向均匀的理想波导下柱面波传播的声波为例，整个位移势满足的波动方程求解分为以下几步：首先，对波动方程两边取傅里叶变换得到频域的波动方程（亥姆霍兹方程）；其次，通过汉克尔变换得到深度分离的波动方程；然后，求解深度分离的波动方程得到$\psi (k_r,z)$；接下，若采用波束积分方法，则直接离散化计算反汉克尔变换，若采用简正波计算方法，则采用留数定理计算；最后反傅里叶变换得到时域波形。
  #求解流程分析
  首先由声学的基本规律，推导出波动方程，速度势$\psi $满足下述方程
  ${\nabla }^{\text{2}}\psi -\frac{1}{c^2}\frac{{\partial }^2\psi }{\partial t^2}=$