C-D无结构网格上的三维自由水面非静水压力流动模型Ⅰ算法 Three-dimensional non-hydrostatic pressure model for free surface flows on C-D unstructured grid Ⅰ . Scheme.pdfVIP

水 利 学 报 2009年8月 SHUIUXUEBAO 第40卷第8期 文章编号：0559-9350(2009)08—0948—08 C．D无结构网格上的三维自由水面非静水 压力流动模型I：算法 胡德超，张红武，钟德钰 (清华大学水沙科学与水利水电工程国家重点实验室。北京100084) 摘要：为准确模拟垂向运动尺度较大的水流运动，本文基于压力分裂模式和半隐方法提出了C．D无结构、；坐标网 格上的三维非静水压力流动模型的构造方法，并详细论述了模型的求解过程。模型先在忽略动水压强项的隐式部 分的条件下求解控制方程，获得比较接近的初始解；再考虑动水压强项的隐式部分并求解一个关于它的三维 Poisson方程，因此计算容易收敛。本文模型与传统的C网格上的非静水模型相比，增加了对水平切向动量方程的 求解，消除了因C网格的简化处理所带来的精度损失。控制方程用有限体积一有限差分法离散，自由水面用水位 函数法处理。对水位刁和动水压强q分别采用不同的隐式因子以兼顾各自稳定和精度的要求。压力分裂模式、8 半隐方法和欧拉一拉格朗日方法(E删)等的联合使用使得本模型具有简单、稳定、高效的特点。 关键词：三维；数学模型；C．D网格；无结构网格；自由流动；非静水 中图分类号：TVl31．2 文献标识码：A 1研究背景 自20世纪80年代以来，在自由水面水流模拟方面，三维水动力学模型获得了越来越多的应用。当 水流的垂向运动尺度不显著时，常使用静压假定以简化求解，这类模型常称为静水压力模型(后文简称 为静水模型)，其中较著名的有POMⅢ模型等，较新的有ELcirc[21模型。然而，在河流、海洋中由于地形 起伏、短波运动、大密度梯度等原因常使得水流的垂向运动与水平运动相比不可忽略，此时静压假定不 等H1、Casulli等b1较早提出了非静水压力模型(后文简称为非静水模型)框架。随后，出现了一系列基于 显式分步、半隐式、隐式构造思想的三维非静水模型，指导思想不同，它们形成的矩阵形式及求解方法也 不同。前两类模型一般最终得到一个压力Poisson方程，可通过显式或隐式迭代法求解；而隐式模型最 终得到一个以水平流速为变量的块三角方程系统，可直接求解№J。 目前国内使用较多的为源于计算传热学的SLMPLE系列算法，它围绕着“正确的压力场应该使计算 得到的速度场满足连续性方程”¨o的中心思想进行求解，由于求解过程需要显式迭代，计算量较大，因而 它在目前国际水动力学领域的应用并不多。为能快速地求解速度压力耦合方程，国外许多学者采用压 力分裂模式，即将压强分解为静水压和动水压分别计算。在压力分裂模式中，计算分为两步，第一步求 解一个关于自由水面的二维Poisson方程获得I临时水位、流场，第二步求解一个关于动水压强的三维 Poisson方程获得最终解，当不进行第二步计算时模型自动退化静水模型。文献[5，8]为基于压强分裂模 式的两个有代表性的三维非静水模型，尤其是文献[8]对其模型进行了较系统的测试。 在另一方面，Adcroft等归3指出传统的c交错网格由于简化处理会使计算精度降低，并提出增加一组 收稿日期：2008．05．26 基金项目：国家“十一五”科技支撑计划项目(2006SAB06舯2)；国家重点基础研究发展计划(2007CB714100) ·--——948·--—— 万方数据 D网格的变量布置即构成C—D网格可以避免计算噪音，后来Zhang等㈨也报导了C网格的类似问题， 但他们的研究都基于静水压力模型，关于C．D无结构网格上的非静水模型的构建方法，目前还未见有 文献报导。本文在深入研究静水压力模型El_tire‘21的基础上，提出了C—D无结构、z坐标网格上的非静 水压力流动模型的构造方法和求解过程。 2控制方程及边界条件 2．1控制方程借助Boussinesq假定、压力分裂思想¨’81，控制方程的形式可表示如下 连续性方程 ，、 ‘，、 ’ ，、 、‘，(1) 婴+窦+