基于vof方法的钱塘江三维非稳态数值模拟.docxVIP

基于vof方法的钱塘江三维非稳态数值模拟 1 对该河段流场的影响 池塘与杭州湾相连。下游水域广阔，以喇叭形连接东海，形成典型的洪水河口。钱塘江潮强流急,泥沙易冲易淤,河床冲淤剧烈,主槽摆动频繁。从20世纪50年代开始,钱塘江河口即按照“减少进潮量,增大山水和潮水的比值以及刷低沙坎”的原则进行治理,采取自上而下治江围涂的工程的措施缩窄江道,以控制深槽平面摆动。20世纪70年代之前主要对仓前以上河道进行治理,之后主要对仓前以下河道实施大规模围涂工程。 涌潮是钱塘江河口潮波传播过程中局部变形的结果,是一种特殊的局部水动力学现象。在涌潮过后的数分钟内,涌潮会有复杂的结构。涌潮形成之初,由于潮后的潮位仍在上涨,因此强度逐渐增大,直到涌潮成为整个涨潮面,此时潮后的潮位随即下降。涌潮较强时会分裂成前后两股,并逐渐减弱,直到湮灭。 为了有效地治理涌潮给人们带来的危害,首先需要弄清楚涌潮过程对河道内流场的影响规律。近年来由于计算机技术和计算方法的迅速发展,利用数值模拟和部分观测结果相结合的方法来探索涌潮的发生和发展已经取得了一定的经验。20世纪60年代初金旦华等应用特征线法求解一维圣维南方程组,把向上游传播的同族特征线的第一个交点作为涌潮起点,在涌潮前后采用Rankine-Hugoniot条件(即激波装配法)求解涌潮前后的水流要素。1985年赵雪华在涌潮尚未形成或涌潮过去之后的区域采用特征差分法计算,而在涌潮形成和传播区运用特征线法计算,并采用激波装配法模拟了钱塘江涌潮。1995年谭维炎等采用基于近似Riemann解的Osher格式求解二维浅水方程,同样采用激波装配法模拟了钱塘江涌潮。20世纪末,苏铭德等应用NND格式求解一维圣维南方程和二维浅水方程,并采用激波捕捉法初步模拟了钱塘江涌潮的形成和发展。21世纪初,潘存鸿等[10,11,12,13,14,15]提出了水位床底法,并结合底坡源项离散技术,逐步建立了有结构任意四边形网格、无结构任意四边形网格和三角形网格下的二维涌潮数值模型的Godunov格式及无结构三角形网格下的二维涌潮数值模型的KFVS模式。 钱塘江作为浙江省第一大河,其河道内部天然物质(如泥沙、盐度等)和人为物质(如污染物)的稀释扩散将直接关系到沿岸流域内居民生活及生产用水的安全。因此,研究钱塘江涌潮条件下河道内的流场分布规律显得尤为重要。近半个世纪以来,人们对钱塘江涌潮模拟已经做了很多研究,但大多是以仓前以下至杭州湾河段作为研究对象,对仓前以上至闻堰河段流场特性的研究还较有限。另外从仓前往上游的河道开始更明显地体现出三维特征,其对河道内流场的影响也越来越大,因此二维或准三维模拟已经无法很好地适应对该河段河道的流场特性研究。由于城市饮用水取水口设置在仓前到闻堰的河段,故对城市饮用水安全来说,该河段的流场特性研究则更为重要。因此对闻堰至仓前河段区域的潮流场进行研究,尤其是对河道内三维流场进行研究,就具有重要的实际意义。 本文在利用实测地形数据建立河道模型的基础上,结合VOF方法,采用商用FLUENT软件对钱塘江闻堰至仓前河段涌潮进行了数值模拟,并在取得良好的验证的基础上对河道内流场进行了分析。 2 单元内变量及其属性 VOF方法通过求解单一的流体动力学方程和确定任意时刻流体单元中每一种流体的体积分数来模拟两种或多种不能混合流体的流动。VOF方法中的两种或多种流体之间没有相互穿插,对增加到模型里的每一附加相,就引进一个变量,即计算单元里的相的体积分数。在每个控制体积内,所有相的体积分数的和为1。只要每一相的体积分数在每一位置是可知的,所有变量及其属性的区域被各相共享并且代表了体积平均值。因此任一给定单元内的变量及其属性可能是其中某一相的变量及其属性,也可能是多相混合的变量及其属性,这取决于体积分数值。对于VOF法追踪模拟自由表面的案例,单元内包括空气和水两相,体积分数分别记为和,则在每个单元中有: 对某个计算单元,存在以下三个可能的情况: (1)uf061w(28)0,表示该单元完全被空气充满。 (2)uf061w(28)1,表示该单元完全被水充满。 (3)0(27)uf061w(27)1,表示该单元部分是水,部分是空气,属水气交界面单元。 自由表面的位置由第3种单元确定。 3 控制体积公式 VOF方法的标准k-uf065模型控制方程为: 连续方程 第q相的体积分数方程 其中为第q相流体的体积分数,满足。 动量方程 k方程 方程 其中:和分别为体积分数平均的密度和粘度: P为修正压力;为湍动粘度,可由k和求出:q 为经验常数;和分别为k和的湍流普朗特数;Gk为由平均速度梯度引起的湍动能k的产生项,可由下式计算: 模型常数的取值分别为 在VOF方法中,控制体积公式要求计算穿过控制体表面的对流和扩散通量并与控制体积本身内部的源项平衡。本