论文部分内容阅读
在声学、电磁学、地震学和一些其他的应用中,有许多在具有分层缓变波导中大尺度的波传播的问题。如果使用直接方法,譬如有限元法、有限差分法等,往往会导致求解大规模的线性方程组,这使得使得算法在内存和计算量的花费非常大,计算效率不高。在本文中,我们针对具有弯曲界面或者弯曲边界的波导,提出了适合于步进计算一般的坐标变换变换的方法,在坐标变换后的坐标系中,使用基于DtN(Dirichlet to Neumann)映射的正向基本解算子步进方法(FFOMM)和逆向基本解算子步进方法(IFOMM)分别求解Helmholtz方程正问题和反问题。在海洋声学中,海洋环境可以看作是一种分层缓变的具有多个弯曲内部界面和多层介质的声波导。在这样的波导中如果要使用步进方法大步长的求解Helmholtz方程,首先需要用适合于步进计算坐标变换将波导中弯曲的界面或者弯曲边界拉直。本文中给出两种坐标变换的方法:第一种是局部的解析正交坐标变换法,局部解析坐标变换用一个解析的积分方程联系新旧坐标,使用牛顿法求解这一积分方程实现新旧坐标之间的转换。坐标变换必然相应的引起方程的变换,我们使用的方程变换将Helmholtz方程转化成为一个可以用步进方法求解的形式。对于具有一个平的顶部、平的底部、n个弯曲的内部界面的二维波导,我们给出了在变换方程系数在各个介质层一般形式的解析表达式,从而完成方程变换。然而,局部解析坐标变换方法需要满足在两个相邻界面之间存在一条水平直线的介质可分性条件。针对不满足可分性条件的情况,我们又提出了一种使用经典求解常微分方程的Runge-Kutta方法计算坐标变换的数值坐标变换方法,这一方法不仅具有更广的适用范围,而且不必引入辅助界面,从而更好的逼近原始问题。对于变换之后的计算坐标系中的正问题,我们使用正向基本解步进方法来求解方程。FFOMM基于DtN算子的单向重建,将边界值问题转化为初始值问题,对于DtN算子和基本解算子满足的方程使用大步长方法离散在传播方向上的变量,用截断特征值展开方法逼近这些算子。FFOMM对于求解大尺度缓变波导中波的传播问题非常有效。在具有分层缓变特性的波导中的波的大尺度传播问题中,反问题有非常重要的实际意义。一般说来,海洋声学中实际的反问题都很复杂。求解海洋中反问题一般使用两类方法:一种是基于正向方法(如FFOMM)的迭代法,另一种为反向方法,反向方法的运算量和计算速度都相比正向方法有很大的优势。基于FFOMM,我们提出了逆向基本解算子步进方法。在IFOMM中使用波导中存在的传播模的个数作为正则化参数,用截断的奇异值分解方法求解步进过程中出现的病态线性方程。数值计算证明,这个方法用于求解Helmholtz方程的反边界值有效,精确并且对于初始波源传播部分稳定,对于海洋声波中实际的反问题研究有着一定的意义。