太湖水动力学三维数值试验研究.pdfVIP

维普资讯 第 l0卷 第 4期 湖 泊 科 学 VoI．10．No．4 l098年 l2月 JOURNALOFLAKESCIENCES Dec．．1998 《，≥ 太湖水动力学三维数值试验7研孑， — — 2．典型风场风生流的数值计算 胡维平 濮培民 秦伯强 ； 学院南京斑疆写雨鬲 所．南京210o0B) 撩力三者时潮流垂直结构的影响．初步揭示丁太湖潮流垂直结构形成的机理． 舟羹号P343．3 馘 fiV’，J 3 ／乎U ， 流场对湖泊中的物质输移有重大的影响，因而湖泊流场的研究一直受到湖泊研究者高度 重视．由于湖流较小，观测相对 比较困难，因而湖流的数值计算得到了湖泊水动力学研究者的 广泛研究．取得了大量的成果．但是这些模型大多数为整层平均二维模型，少数为两层或为理 想地形多层模型0 ]．这些模型在研究诸如泥沙，藻类等存在明显垂直梯度的物质输移和空间 分布时．就显得过于粗略．胡维平等 建立的 坐标下的太湖风生流三维数值模型因可对垂直 方向进行多层细划分，且各网格点层数一致，比较适于这类研究的湖流数值计算．因而本文用 该模式探讨太湖各典型风场 SE，NW．SW 驱动的湖流，以期为研究太湖诸如泥沙、藻类的输移 和空间分布和机理提供必要的基础． 1 湖水运动控制方程 据胡维平等研究结果 ，在湖底为0，湖面为1的归一化压缩 坐标下，湖水的运动控制方 程为： + + +善dt—o (1) dI 詈+罢“+詈+”差一，口一一嘉+A[喜+翥]+耋[A。妻]+ z 害+“耋+考 耋+^一 善+[+舅]+ ^妻] ㈣ 这里 ，方向分别为在水平面投影指向正东、正北，f为时间变量，为离开湖底的相对高度， 向上为正．f为水面相对平衡位置位移，g为重力加速度，，为科氏力常数，“，，为流速在 ， 方向分量 ；￡^ 为 、v方 向控制方程水平扩散项变换产生的偏差项 ，它们为 ’譬嚣皆 若甜嚣翠 ： ：i兽象学。院五重点项 一一。联台资礓 维普资讯 4期 胡维平等：太潮水动力学三维数值试验研究—— 2．典型风场风生流的数值计算 一 (砉 c一警+窑毒c一警。一z詈再1ic差一-．dH)+ 考击c一 +砉器c一 +~u甭1t~n一等_z詈· H1c詈一等+亲击c券一等) ’ 的表达式和 的表达式相同，只要把 “换成 即可，H为从湖面到湖底的深度 ·边界条件可 简化为： 一 1时 ‘= 0．0 P刍(等，窑)一(，)=cb√“+(，) = 0时 ‘= 0．0 刍(塞，宝)=(，)=c√讲+谢(，) 这里P为水体的密度， 为空气密度， ， 为湖底 ，Y方向的切应力， ，‘为湖面 ，Y方向 的风应力 ·cb，c 为湖面和湖底的拖曳系数 ·方程的侧边界条件为法向速度为零 · 数值计算用的差分方程为交错网格差分方程，变量 “放在 方向和 方向的半整数阿格 点上，变量 放在)r方向和 方向的半整数网格点上， ’，h，H，放在整数网格点上，差分用 中央差分，垂直湍流扩散项用半隐半显式差分．固差分方程十分复杂此处略去 ． 2 计算结果与分析 数值计算的空间步长 、Y方向均为 1kni，方向步长为0．2，时间步长为 120s．水平扩散 系数取 5×10c‘m。／s，垂直方向的湍流扩散系数取4．Ocm。／