计算水力学06.ppt

第 6 章 二维水流模拟 2 数值方法 第2节 正交边界拟合坐标变换 2 求解方法 第3节 河道二维水流模型 2．边界条件 3 河道二维水流计算模型 动量方程的差分 方程求解 第4节 ADI方法 李光炽 计算水力学 采用线Gauss-Seidel迭代法，求解步骤如下 根据边界条件和工程计算的要求，事先假定区域的初始条件—即给定初始计算网格。 在计算平面上对方程（13）进行离散，由于式(13)为一非线性的偏微分方程组，采用线性化的方法将其离散。其离散过程中，主要是将 用初始计算网格进行计算，而其它项则采用隐式中心差商进行离散。 李光炽 计算水力学 3.对于离散后的差分方程，用线Gauss-Seidel迭代法求解新的节点坐标值。 4.根据正交条件，滑动边界节点坐标值，以便边界处满足正交条件。 5.重复2、3、4的步骤，直到计算的曲线网格全部满足正交条件为止。 李光炽 计算水力学 李光炽 计算水力学 1.基本方程 描述平面二维水流运动的基本方程组为： (1) 李光炽 计算水力学 (4) 李光炽 计算水力学 式中u、v分别为沿ξ和η方向的流速， 分别为曲线网格的长度和宽度， 为曲线网格的面积。 李光炽 计算水力学 与u、v之间的变换关系为： 通过正交变换，把原来在x——y坐标系统中利用方程(1)求解变量Z、u、v变为在ξ——η坐标系统中利用(4)求解Z、u*、v*。下面为了书写方便，略去下标*。 李光炽 计算水力学 (1) 固壁边界 采用不穿透条件。 V·n = 0 （n为固体边界的法向矢量） (2) 自由边界 上边界： Z上为已知 或: u为已知 或: Ｑ为已知及假定水位无横比降 下边界： Z下为已知 或: u为已知 或: Ｑ为已知及假定水位无横比降 李光炽 计算水力学 节点布置 图2 物理平面上的节点示意图 ξ+δξ η+δη ξ η x y ξ η 水位 ξ 流速V η 流速u 李光炽 计算水力学 图3 变换平面上的节点示意图 3 2 1 2M+1 1 2 3 2N+1 ξ η ：流速U节点 ：流速V节点流速V节点 ：水位节点 李光炽 计算水力学 为了便于边界条件的处理，变量采用交错布置，即在网格中心布置水位变量；在网格四边布置相应的流速变量。图(2)、图(3)为物理平面及变换平面上的节点变量示意图。从图(3)中可见，流速u有N·M个节点，流速v有(N+1)·(M+1)个节点，水位Z有(N+1)·M个节点 李光炽 计算水力学 方程离散 连续方程式的差分 在点(2k＋1，2j)对方程(4)的连续方程，采用如下格式差分： 非线性项：gηhu和gξhv的线性化采用如下公式： 李光炽 计算水力学 项 采用中心差分，将上述差分近似和边界条件代入到连续方程可得如下线性差分方程组： 上式写成矢量形式为： A1k·Z2k-1 + B1k·Z2k+1 + C1k·Z2k+3 + D1k·V2k+1 + E1k·U2k + F1k·U2k+2 =H1k (k=1,2,……,N) 式中：A1ｋ、B1k、C1ｋ、E1ｋ、F1ｋ为M·M矩阵，D1k为M·（M－1）矩阵，H1ｋ为M维矢量。 李光炽 计算水力学 李光炽 计算水力学 将上述的差分近似及以边界条件代入到ξ方向的动量方程中得如下线性方程组: 在点(2k+1, 2j+1)对η方向动量方程类似于ξ方向动量方程的处理可得如下线性方程式： 李光炽 计算水力学 A2k·Z2k-1 + B2k·Z2k+1 + C2k·U2k-2 + D2k·U2k + E2k·U2k+2 + F2k·V2k-1 + G2k·V2k+1= H2k (k=1,2,…….,N) (9) A3k·Z2k+1 + B3k·U2k + C3k·U2k+2 + D3k·V2k-1 + E3k·V2k+1 + F3k·V2k+3 = H3k (k=1,2,…….,N) 式中：A2ｋ、B2ｋ、C2ｋ、D2ｋ、E2ｋ为M·M矩阵,F2ｋ、G2ｋ为M·（M－１）矩阵，A3ｋ、B3ｋ、C3ｋ为(M-1)·M矩阵，D3ｋ、E3ｋ、F3ｋ为(M-1)·(M-1)矩阵，H2ｋ为M维矢量,H3ｋ为(M-1)维矢量。 李光炽 计算水力学 矩阵追赶法是一种针对二维河道水流计算特性所提所出的一种求解离散后的线性差分方程组的一种算法 V1=0 上边界水位无横比降无横向流速（11） A21·