湖泊平面二维水流与水质分布研究.pdfVIP

湖泊平面二维水流与水质分布研究 冯民权1周孝德1赵克玉1 (1．西安理工大学，陕西西安，710048) 摘要：在湖泊大水体研究中引入坐标变换、交替方向、交错网格、线性递推、隐 格式的平面二维水流与水质模型，可完全适应不规则的湖岸边界，实现水流与水 质的快速预报。在新疆博斯腾湖的风生环流与矿化度分布研究中取得良好效果。 关键词：湖泊；坐标变换；水质；交替方向；交错网橇≈隐式。 卜 1 产、 ， 一 卜、 1引言 型也得到研究。这类模型将水库湖泊进行分层分块概化，类似于将一个连续函数 离散表示，因此，模型存在一定误差，而且各层各块之间的边界条件处理比较困 难。Simons等(1977)建立了比较细致地考虑水流运动的三维水温模型，但受实 测资料及计算条件的限制，其实际应用受到限制，而在三维模型基础上的垂向平 均和横向平均二维模型则得到成功的研究和应用。 南京地理所、中国环科院等曾作过平面二维湖泊的水质模型研究。郑邦民做 过湖泊型温度流的数值分析。吴江航、朱咸、刘兰芬做过浅水型冷却池中环流和 温度分布瞬态数学模型的有限元数值分析。大连工学院做过近海浅水环流有限元 数值计算。陈小红(1997)提出水动力与水温水质耦合的水库垂向二维水质分布 预测模型。雒文生等(1997)采用k一占双方程湍流模型计算了水库立面二维流动 与水温水质。江春波(2000)提出了一种预测河道型水库中流速、温度和悬浮污 染物质分布的立面二维数学模型。 本文将坐标变换、交替方向、交错网格、线性递推、隐格式的平面二维水流 与水质模型引入湖泊大水体研究中，并考虑表面风应力的影响，力图完全适应不 规则的湖岸边界，实现水流与水质的快速预报。该模型被应用于新疆博斯腾湖的 风生环流与矿化度分布研究中。 2数学模型 对于湖泊、浅水水库而言，由于水深较小，可以认为沿水体垂直方向掺混比 较均匀，用垂向平均化的二维不可压缩、粘性流体运动方程来描述浅水水域的流 ‘ 场和浓度场分布是较为合理的。 2．1控制方程 【作者简介】冯民权(1964—’，男，山西永济人，高工，博士生，主要从事计算水力学与环境水力学的研究． 183 a)水流连续方程署+昙(肌，)+专(日y)=o．……一乓)a)水流连续方程瓦+石妒“，+荀p’，”………刈， b)水流运动方程 饥………(2) 詈+“罢+v祟+g妾+丘一Qv=刃2一li+饥………(2) 瓦栅瓦w瓦+g瓦+乃一‘2”2内 饥………··(3J 安+甜安+v安+g祟+∥十m：内一2V+饥………．．(3) 百+甜瓦+”面+g瓦+∥十‘Ⅺ。田 式甲，z为水位；丑为水深：“、V为x、Y方向的垂线平均流速；f为阻力 系数，，=z+^，z为河床底部摩阻系数，石=g√矿+哆‘日，^为局部摩 阻系数，五=孝．√Ⅳ2+”／2／Ad，￡为局部阻力系数：o为哥氏力系数；v为紊流涡粘 性系数；x为风应力系数；“。、V。为，、y方向的风速wl，I#2。 c)水质方程 旦字+旦等坠+皇髻堕=丢(舰：塞)+昙(艇j争…………………(4) 2．2控制方程的坐标变换 经过变换后一般曲线坐标下的水动力方程为： 计了1∽(Flu)f一(堋小专【咄(协)，+≮(剐小0．．．．……………(5) 184 在(5)、(6)、(7)式中，左边第一项为时间变化项，它们是非恒定项，如果 没有这一项，流动则是恒定的。在动量方程(6)、(7)中，左边第二项为对流项， 它具有非线性性质。左边第三项为反映重力作用的水位变化项。左边第四项为摩 阻项，它包括底部摩阻和局部摩阻。左边第五项为哥氏力项。右边第一项为涡粘 性应力项，右边第二项为风应力项。 其中“，V仍是原来x，瑚方向的速度分量，而不是孝，玎轴方向的切线分量。Zt，0，z。，‰， V。等分别是z对r，乎，叮；“对孝；v对叩的导数． 经坐标变换后一般曲线坐标下的水质方程： —O(石H_c)+了1抄。