曲线坐标下三维水动力学生态综合模型建立.pdfVIP

水 利 学 报 2005 年 8 月 SHUILI XUEBAO 第 36 卷 第 8 期 文章编号：0559-9350(2005)08-0891-09 曲线坐标下三维水动力学生态综合模型的建立 倪浩清，李福田 (中国水利水电科学研究院，北京 100044) 摘要：为了获得高分辨率、高精度及快速的数值模拟计算成果，以便解决复杂边界的工程实际问题，本文研究了通用 形式表示的湍流代数应力模型流动基本方程组在拟合坐标系中的转换。通过该模型方程组的转换，并结合具有生物化 学作用过程的生态子模型，建立了曲线坐标下三维水动力学的生态综合模型，便于复杂边界工程问题的方程组求解及 对未来情况的预测和评价。 关键词：湍流基本方组转换；贴体曲线坐标；水动力生态综合模型 中图分类号：X171 文献标识码：A 本文中所提出的水动力学生态综合模型实质上是污染物质在水环境中因物理、化学、生物作用及影响因素 之间相互关系的数学描述。这方面比较有代表性的模型有Grenney[1]在QUAL-I综合模型基础上发展的QUAL-Ⅱ综 合水质模型；倪浩清等[2] [3] 提出的“三维湍流代数应力/热流模型(简化的统一二阶矩湍流模型)”；沈永明 发表 的“近海环境准三维k-ε湍流流动和水质的精细模拟”；雒文生等[4]发表的“垂向二维湍流与水温水质耦合模型”； 李清雪[5]提出的“海湾浮游生物及氮营养盐生态水动力学模拟”等。这些模型对水环境数学模型发展均取得积 极的推进作用，但在这些模型中没有综合考虑海湾的弯曲边界、各向异性和生态动力作用对湍流数值模拟计算 的影响。因此，考虑到湍流、温度、生态及流动边界相互耦合影响的关系而建立一个完整的曲线坐标下三维水 动力学生态综合模型是十分必要的。 1 控制方程在拟合坐标系中的守恒方程 1.1 模型的选择与拟合坐标系 由于江河湖海的岸边界弯弯曲曲极不规则，其中的流体质点受到离心力及 各向异性湍流雷诺应力场的作用而导致产生二次流，它与弯道轴线方向的主流合成各向异性的螺旋流，这种二 次流能导致主流流场、温度场及污染物浓度场再分布，采用湍流代数应力/热流模型，则能模拟出离心力和各 向异性湍流应力场所引起的二次流。因此，选用具有考虑生化反应影响的各向异性湍流代数应力/热流模型作 为水环境数学模型是符合天然水域生态水动力学特性的。同时，针对天然河流、海湾边界弯曲不规则性，采用 拟合坐标系变换办法，则是解决数值模拟计算问题中重要途径之一。 1.2 变换坐标系中的守恒方程式 时均化的N-S方程及质量输运方程可以写成如下的笛卡儿坐标向量微分 方程的守恒形式： ∂M ∂F G H ∂ ∂ + + + Q (1) ∂t ∂x ∂y ∂z 收稿日期：2004-04-02 作者简介：倪浩清(1931-)，男，江苏无锡人，教授级高级工程师，主要从事水利水电工程、海洋工程和火、核电冷却水工 程中的水动力学、波动力学及水环境动力学研究。 1 水 利 学 报 2005 年 8 月 SHUILI XUEBAO 第 36 卷 第 8 期 u