一类求解helmholtz方程高精度紧致差分方法分析-analysis of a high precision compact difference method for solving helmholtz equation.docxVIP

第一章绪论一、 引言实际生活中的很多现象都可归结为界面问题，如不同传导率的两种材料的拼接，相同物质在 不同状态下混存（如水和冰），血液细胞在血管中的流动以及地下水流经不同的物体（如石块、 海绵）等。界面问题涉及环境科学、物理学、生物数学等学科领域，引起了许多科研工作者的浓 厚兴趣，如何正确处理界面问题成为解决实际问题的关键因素之一。许多用以刻画界面问题的模 型方程常常是常微分方程或偏微分方程时，由于模型中的一些物理参数是不连续的，如两种不混 溶流体中的密度、粘性等，源项也是奇异的，因此这类问题的解往往是不光滑的甚至是不连续的。 许多针对具有光滑解的问题所设计的标准有限差分方法在数值求解此类微分方程时计算效果较 差，常常在界面附近达不到预期的精度要求，甚至有时完全不可用。因此基于坐标变换和网格映 射技术的有限差分方法被提出，但这常常需要寻找适合的高精度变换矩阵计算方法进行坐标变换 和网格映射，即便是生成一个高质量正交性好的网格都需要耗费大量的时间进行迭代运算。而随 着处理问题复杂程度的不断提高，生成理想的结构化网格系统变得十分困难。因此发展一类数值 求解界面问题的高效，稳健的差分方法具有一定的实际意义。 二、处理界面问题方法综述1977 年 Peskin 提出了浸入边界方法[5-8] (IB,The Immersed Boundary Method ) 用于模拟血液 在人体心脏中的流动。该方法兼具数学建模和数值离散的特点。它采用欧拉变量描述流体的动态， 利用拉格拉日变量描述结构的运动边界，用光滑 Detla 逼近函数通过分布节点力和插值速度来表 示流场和结构物的交互作用。整个流场计算都使用笛卡尔网格，而不是按照物体形状生成复杂的 贴体网格，无需处理从物理平面到计算平面的坐标和网格转换问题，因而可以大大提高计算效率 而且节省了网格生成所需的时间。尤其对于动态边界问题，它无需在每一时间步长实时更新网格。 由于该方法高效、易于实现引发了大量的后续工作，被广泛应用于解决生物数学及计算流体力学 领域的问题[9,10]，一些学者还将该方法结合涡方法及网格自适应方法来解决一些实际问题[11,12]， 但该方法在界面处通常只能达到一阶精度，想要达到二阶计算精度，Beyer 和 Le Veque 在[13]指 出需谨慎选择所使用的 Delta 函数，且目前还没有完整的关于 IB 方法收敛性的理论证明。 八十年代，Bell 和 Dawson 等人 Samarskii[14]的基础上提出求解椭圆界面问题的调和平均法 (The harmonic averaging for discontinuous coefficients)[15,16],该方法在求解二维问题时不 仅需要花费很大的工作量求解方法中所涉及的积分项，而且在实际计算中并不能达到理论上的二 阶精度[17]。与此同时，对于在不规则区域上求解泊松方程或双调和方程，Mayo 与他的合作者[18-20] 率先提出可结合积分方程，在不规则区域上通过修正已有的差分方法使其在求解该类问题时整体 达到二阶精度。该方法人为将不规则区域嵌入一个更大的矩形区域，所求解的界面问题包括边界在内被重新在该矩形区域上定义。在不规则的网格点上采用泰勒展式，在不规则的边界处则用解 的积分表示，为使用泊松快速求解法，他们从积分方程中推导出跳跃条件，用于整个矩形区域中 的所有网格点计算。该方法还可通过与快速多极法相结合得到加速[21]。 1994 年 Le Veque 和 Li 创造性的发展了浸入界面方法（IM, The Immersed Interface Method） [22,33-35]。用于处理带有不连续系数和奇异源项的椭圆型方程，该方法从因其能够处理和解决流通 量上的不连续性而备受关注。该方法从问题本身的物理背景或受约束的微分方程本身预先给出跳 跃条件，得到更多的界面关系式，在远离界面的区域，通常采用标准的有限差分格式或有限元方 法，界面附近的网格点或元则依据得到的界面关系式对数值方法进行修正，通过比较在无穷范数 意义下的误差量级，浸入界面方法往往具有整体的二阶收敛性。若系数、解及其流通量的不连续 性消失，则浸入界面方法就是标准的有限差分方法或有限元方法。Li 与他的合作者成功将浸入界 面方法用于求解不可压 Stokes 方程组[23-26]及带有奇异源项的 Navier-Stokes 方程组[40]。另外， 对于所求解得界面问题或定义在不规则区域上的问题当流通量或解的跳跃条件是以导数形式给 出时，界面条件关系式的推导会因此变的相对比较复杂，要涉及到高阶导数，大大增加了问题求 解的复杂性。因此 Li 与他的合作者还发展了广义浸入界面方法（AIIM, The Augmented Immersed Interface