改进σ坐标下的三维水流传质数学模型 3-D flow and matter transport model for improved σ-coordinate system.pdfVIP

第18卷第5期 水科学进展 v01．18，No．5 20昕年9月 ADVANCESINWATERSCIENCE sep．，2007 改进d坐标下的三维水流传质数学模型 董 壮1，龚 政2，韩 青2 (1．河海大学水利水电工程学院．江苏南京21∞9B；2．水利部太潮流域管理局，上棒卸0434) 摘要：针对J坐标变换的不足，引A垂向坐标变挟系数．通过严密的数学推导，建立了改进的三维本流传质数学 模型。鲭出了改进后的控制方程及相应的边界条件、初始条件。采用有限差分穗对模型进行了离散，针对瓯江河口 局部海区进行了潮流、盐度过程的验证计算，取得了较好的成果。 美■词：三维水流；口坐标；散学模型；FDM 中围分类号：Tvl33．2 文献标识码：^ 文章编号：J∞l—679l(20町105明17—a7 在三维水流数学模型中，目前应用较多的是a坐标变换[1_3]。口坐标变换的主要优势在于自由水面的处理 相对简单，但也存在一定不足，如变换后的控制方程形式更为复杂，计算结果需加以转换，使计算难度增大， 可能带来新的误差；某些情况下容易发散；地形变化复杂处误差较大等。除口坐标变换外，绝对分层的方法 也有一定的应用，与一坐标变换相比。这种方法控制方程形式相对简单，不必进行坐标变换，同时与实际地 形更为吻合。但绝对分层需对自由水面单独处理，也给计算网格的划分带来一定困难【4】。 本文从实际工程的需要出发。吸收上述两种方法的特色，将a坐标变换加以改进，既保留原有a坐标变 换便于自由水面的处理的优势。又结合绝对分层方法与与实际地形吻合更好的特点，建立新的数学模型。 1 盯坐标系中的三维水流数学模型 一坐标变换其形式一般为 一=宁 (1) 式中z为笛卡儿坐标系中的垂向坐标；≈为自由水面相对于基准水位的波动高程；D为全水深(D=^+々，^ 为平均水位到水底的距离)．自由水位处==o。变换后，自由水面口=0，水底d=一1．0。 采用静压假定，经过上述变换后，雷诺方程组的形式为【2】： 娑+孕+攀+翌：o (2) oI 04 oo o， c，)”， 鲁+警+等+等=．s磅+咿雾+们雾+丢【告笔】百+百+百+百一gu蠢+叫孬刈^Ⅳ孬+万【万万j 等+警+等+智=一s咯ⅥD骞¨D骞+易【告舅】百+百+百+百一∥苟”∥孬”∥孬+-I万苏j c。，”7 竽+孪+掣：o (5) o‘ 口* oy 沿垂向积分得到。在以上各式中，l为时间；￡，y，d分别为水平方向及垂向坐标；u，，m分别为#，y，口 pd口；‰。q分别为水平方向及垂 方向流速分量；i，i分别为*，，方向垂线平均流速，i=Ind口，i=I 向的紊动粘性系数，一般认为z，y方向的水平紊动枯性系数相等并取为常数，这一做法往往用于大体积水域 牧稿日期：加06．O车．∞；修订日期：2006-11—2l 作者简介：t壮(1973一)，男，天津人，讲师。主要从事计算流体力学方面的研究。Bm矗丑：叠d呻萨0晒@y且h∞帅，皿 万方数据 718 水科学进展 第18卷