复杂边界下三维水流数学模型的建立和验证.pdfVIP

第4期 黄国鲜等：复杂边界下三维水流数学模型的建立和验证 67 l平面正交坐标、立面盯坐标系统下三维流动控制方程 连续方程掣+专掣+专掣+警=o ㈤ e方向动量方程型号≠堕+号【型掣+型掣】+堕等字堕+ 学(‰碧一n，等)～咖≯=一观+矾 ． (z) _方向动量方程堕号≠堕+专[型掣+型掣】+型号字堕+ 竽(％誓一‰誓)+曲加一卯：…： (3) 一方向动量方程鱼孕+专[掣+掣】+攀立=一观+。B (4) 式中I为时间，P为水流密度，D为水深，J=^，也，^．和^：分别为平面正交网格在#，口方向上的拉梅系数， ‰，％，矽分别为模型坐标系数下水流在{，q，a方向上的流速分量，；为水位，口坐标下垂向伪流速妒和实际垂 向流速”的相互转化关系为： Ⅳ=D·”=”一‰(篓+一警)一u，(葛+一鬻)一(薯+一·警)cs) 别为水流紊动动能和紊动动能耗散率，采用标准^．e紊流模型求解水流紊动粘性系数，模型坐标系下＆，e方程 具体形式见文献[4]。 2模型的边界条件 2．1 自由水面的求解及其边界条件 模型计算自由水面的方法有：(1)钢盖假定方法；(2)使用集成在三维模型中的平面二维模型来计算水位； (3)采用wu㈨提出的基于静压假定的二维泊松方程求解自由水面；(4)在方法3的基础上考虑动水压力和二次环 流对自由水面的影响的求解法“1。数值试验表明：钢盖假定法较为简单且使得模型的计算容易稳定和收敛；在较 长河段的非恒定三维水流模拟过程中采用一、二维模型计算自由水面就能满足要求；在空间尺度较小的恒定流计 算中采用动水压力和二次环流修正的泊松方程求解法能够使得自由水面的计算合理，采用交替隐式法求解自由 水面的离散方程。自由面上变量u}，‰，，p一，％裙垂向的一阶导致为O，￡在自由水面上的边界条件由经验公式 直接给出，即：塞=o，关=o．赛=o，尝=o，蒜=o川一=胖，(o．43口， 2．2床面和岸壁边界条件 为了简化计算和提高计算效率，本文采用标准壁函数来计算第一个内网格节点P处紊动牯性系数‰紊动 动能＆和紊动动能耗散率e。 2．3进口边界条件 在进口断面各计算节点上给定变量(u。，”。，”，t，￡，s)的具体数值，即： Ⅱ{2“缸，u{=“枷，坤=坤“，知=k，￡=e_ 2．4出口边界条件 出口处流速、动水压力、紊动能量和紊动动能耗散率沿纵向的一阶导数为0，即： 警=o；警=o；警=o；等=o；篓=o；篓=o，下游给定水位控制边界条件：f=k。 水力发电学报 2啷年 3模型的变量布设、离散和求解 为j求解方便，3}哥控制方程写戚统一形式： 警+专【掣+掣】+学= 。 其中：输运方程的源项函包括对流源项s¨、扩散源项s蜥、和其它源项s№(如压强梯度源项、科氏力源项 及其物质的产生和降解源项等)： S=s。。-+s^甜+s‘。h=s讲+sH声 (7) 为了使得计算易于稳定和收敛，系数sⅣ≤O，根据模型控制 方程的形式不难写出s．，s_，sH的表达式。采用控制体积法对 方程(8)在控制体积(如图1)进行离散得到线性方程如下： 口P拳，=口E≯f+口r乒r+口一庐Ⅳ+口5庐5+ (81 … 口r乒r+口口≠口+6+口二≠； o。，口，，口。，Ⅱ。，钾，d。为不同方向的邻节点变量对计算节 点P的影响系数．其取值与流速分布、网格尺寸和离散格式有 关，6为源项，n：为前一步计算结果的影响系数。本文