模型的建立-Read.DOC

一类全场数学模型的建立 王泽良 中国水利水电科学研究院水环境所 摘要 对流项在水动力学模型中起着重要作用，它的非线性特征，不仅使得求解这类问题变得复杂，而且也正是由于它的非线性，使得垂向积分的数学模型的表达形式各异。本文研究了不同流速分布水深积分水动力学模型对流项的表达方式，在此基础上结合湍流模式，建立了一类水深积分的全场数学模型（包括风生流的数学模型）。并采用全场模型对香港拟建中的南丫新电厂对下尾湾流态的影响进行了数值模拟研究, 把数值计算结果与物理模型试验结果进行比较。 关键词 全场模型，流速分布，离散作用。 1 前言 深度平均的二维水动力学模型在水环境研究领域应用十分广泛。由于模型中存在一些综合系数，使得模型应用中的参数确定是一直成为一个难题。为了研究热水或污染物在大水体中各向异性的输移扩散特性，1991年提出了深度平均的全场模型。这些模型的建立为更加精确地模拟大水体中热水或污染物的输移扩散提供了一类湍流封闭模式。但是在深度平均的二维水动力学模型中，由于积分作用引起的离散作用在模式的建立过程中影响颇大。但在这类全场模型中，没有细致考虑由于积分引起的离散作用。使得模型在应用效果仍缺乏应有的说服力。水深积分引起的离散作用的产生项，主要是非线性对流项产生的，对对流项进行深入的研究，建立相应的模式，对于确立完备的全场模型显得尤为重要。本文在假定了多类型速度分布的基础之上，推导了各类速度分布条件下，离散作用项的模拟形式，建立了基于修正对流项的一类全新的模型封闭的二维水深积分水动力学模式。采用物理模型试验结果对本文提出的模型进行了验证。 2 模型的建立 2.1 深度平均的水动力学模型 连续方程 （1） 运动方程 （2） 2.2 深度平均的湍流模型 湍动能方程 （3） （4） 湍动扩散量方程 （5） （6） 2.3 方程的模式化 2.3.1 项的模拟 根据推导条件可知，也既是。不同流速剖面的积分形式为， (1) 对数流速分布 假定，为无量纲垂向坐标（），为水深，为摩阻流速，为常数,为卡门常数。 （7） （2） 抛物线型流速分布 ①风应力处于次要因素 假定，为无量纲垂向坐标（），为水深。（需要说明得是，这里假定自由表面处垂向速度梯度为零，不考虑风应力的作用）。 （8） ②风应力处于主要因素（风生流） 假定为无量纲垂向坐标(），为水深。（需要说明得是，这里考虑风应力在自由表面处的影响）。自由表面应力条件： （9） （10） （3）指数速度分布 ①1/7次率 假定，为无量纲垂向坐标（），为水深。 （11） ②1/8次率 假定，为无量纲垂向坐标（），为水深。 （12） 表 1 修正系数的选择 （这里不含考虑风生流的情况） 对数流速分布 抛物型速度分布 （风应力为次要因素） 指数型(1/7) 指数型(1/8) 6/5 64/63 81/80 风应力处于主要因素（风生流），需要说明得是，这里考虑风应力在自由表面处的影响。自由表面应力条件： （13） （14） 式中： 2.3.2 项的模拟： 采用雷诺比拟，可得：, （15） 2.4 深度平均的湍动全场数学模型 连续方程 （16） 运动方程 （17） 湍动能方程 （18） 深度平均湍动扩散量方程 （19） 式中： （20） , （21） （22） （23） 表2 通用常数 1.44 1.92 0.09 1.0 1.3 风生流的二维水深积分平均的水动力学模型: （24） 3 模型的验证 在南丫岛的西海岸，建有香港电灯公司电厂。由于发展的需要，香港电灯公司拟在南丫岛现南丫电厂南侧海域填海造地，修一新发电厂。这些新的填海造地计划，必将在一定程度上影响海湾的流态，改变海流的输移规律，为了解填海工程的影响，进行了物理模型试验研究。本文采用此次物理模型试验结果对本文提出的模式进行验证[4]。 3.1 物理模型试验 采用多通道粒子图像测速系统，该系统包括示踪粒子、图像接收与记录系统（12台CCD摄像机），信息处理与存储系统（图像卡、计算机及配套软件）。 测量范围及精度 使用于平面二维流速场测量 最短采集间隔时间：1/25秒 图像分辨率：768X576 测量误差５％ 使用该测量系统可实时监测水面流况，并获取不同潮时模型水面流场