浅水问题中方程型式对数值计算的影响.docVIP

浅水问题中方程型式对数值计算的影响 董俊哲 大连理工大学工程力学系，辽宁，大连，116023 julius@dlut.edu.cn 摘要：浅水方程通常假定为无旋流动。基本方程为Boussinesq型方程，待求解的变量除波高外，可以是速度势或者速度，表达形式有多种。本文从零阶近似的Boussinesq方程出发，探讨了不同的方程型式对有限差分格式的影响，提出于方程相称的差分格式计算方法。 关键词：浅水波，Boussinesq方程，稳定性，有限差分。 1. 前言： 零阶近似的Boussinesq浅水方程为双曲型方程[1]，其有限差分显示格式的数值稳定性通常表现为CFL条件[2]。从物理和几何上解释为，差分过程中节点数值对网格的依赖区域应包含方程解本身的依赖区域。零阶近似方程实际上是方程解的主要部分，不仅反映了方程解的主要特征，还在数值模拟过程中对稳定性条件有至关重要的影响。由于在方程的推导过程当中，方程有各种不同的求解变量和不同的表达形式。不同的方程形式导致不同的差分格式，受方程本身的影响，在有些情况下即使满足了CFL条件，差分格式仍然不能够确保满足稳定性条件，使得数值计算不能够收敛。虽然高阶项的添加可以在一定条件下，使差分格式收敛，但由于高阶项的影响在计算过程中只占很小的比例，使得时间步长往往要取得很小才能收敛，严重降低了计算效率。本文对零阶近似方程的不同形式及其相应的差分格式进行了分析，提出了有效的计算方法。 2. 基本方程 Boussinesq方程的零阶近似为： （1） 其中为波高，为速度势，为静水深，为重力加速度。 令，得方程的另一种表形式： （2） 若消去方程（2）中的或，则有： 和 （3） 以上就是浅水波动方程零阶近似的三种表达形式。（1）式和（2）式分别为高阶近似方程的一部分，（3）式为波动方程，波速。 方程性质及差分格式 为了方便考察方程，将方程无量纲化。令：，，，，其中带（‘）为无量纲参数。为方便起见，下面的方程中省略（‘）号。于是，方程（1）、（2）、（3）变为： （4） （5） 和 （6） 方程（6）为典型的波动方程，具有特征线和，显式差分格式的稳定性条件为CFL条件。 对方程（5），令：，则有： （7） 其中：。由于矩阵有实特征值，所以该方程为双曲方程。 如果空间差分取后差，则方程（7）的直接显式差分格式为： （8） 其中为时间步长，为空间步长，为时间步，为空间节点。对该格式进行Fourier分析，令，得增长矩阵为： （9） 由于该矩阵的谱半径，不满足von Neumann条件，故不收敛。如果空间差分采用前差，也同样不收敛。究其根本原因，作为双曲方程，特征线起着重要的作用，因此方程（7）应采用迎风格式。方程（7）的迎风格式与矩阵的特征值有关。特征值的正负号代表着“风向”，令： ，， 则 。为对角阵，对角上的元素即矩阵的特征值。再令： ， 则方程具有迎风差分格式。即：方程 (10) 中的和分别有特征线和，空间差分分别取前差和后差即得到迎风格式，稳定性条件为CFL条件。待求参数和可以由和线形组合而成。实际上，虽然方程（6）从方程（5）推导而来，但方程（5）并不象方程（6）那样，沿未知数和特征线和，而是沿和有特征线和。因此，当求解以波高和速度为未知参数的高阶近似的Boussinesq方程时，可以将关于和的两个方程变换为未知量为和的两个方程，再分别对一阶空间导数取前差和后差，得到迎风格式下的差分方程求解，这样的求解效率要提高很多。对于方程中高阶项和非线性项中难以变换的部分，由于对迎风格式影响较小，仍可使用每一个时间步所求出的波高和速度。 对方程（4），令，则有： ， （11） 其中，。该方程相当于有源的热传导方程，抛物型。方程（11）的古典显示格式为： （12） 对该格式进行Fourier分析，得增长矩阵： （13） 可以得出，谱半径，不满足von Neumann条件，不收敛。 将（12）式中的第二式代入第一式，则有： （14） 很显然，这是一个不收敛的关于波动方程的差分格式。事实上，如果（14）式中的第二项取时间步的物理量，在满足CFL的条件下，该格式就可以收敛了。实现这一步骤的作法是，在计算过程中，先求解差分方程（12）中的第2式，将得到的时刻的速度势代入第一式，这样的作法相当于求解差分方程 （15） 差分方程（15）实际上是一个半隐式的差分格式，所对应的波动方程的差分格式为： （16） 差分格式（15）在满足CFL的条件下是收敛的。与差分格式（10）相比较，格式（15）更容易操作。 4. 结论： 浅水方程