胶州湾三维变动边界的潮流数值模拟.docVIP

360 海 洋 与 湖 沼 32 卷 最大流发生时刻 , 随深度增加有所提前, 底层较表层提前 10min 左右。 表1 T ab. 1 站位 C1 ( i, j 40, 20 C1 C6 站 M 2 潮流椭圆要素 的垂直分布 C6 station 椭圆长轴 0. 347 0. 326 0. 302 0. 269 0. 211 0. 292 0. 273 0. 252 0. 224 0. 174 0. 359 0. 341 0. 319 0. 288 0. 230 0. 346 0. 326 0. 303 0. 272 0. 216 0. 803 0. 791 0. 768 0. 724 0. 620 0. 468 0. 449 0. 425 0. 389 0. 319 椭率 0. 141 0. 138 0. 132 0. 126 0. 123 0. 151 0. 147 0. 143 0. 140 0. 136 0. 119 0. 106 0. 095 0. 084 0. 073 0. 135 0. 126 0. 118 0. 111 0. 101 - 0. 156 - 0. 163 - 0. 168 - 0. 171 - 0. 173 - 0. 176 - 0. 200 - 0. 221 - 0. 243 - 0. 266 The vert ical dist ribution of M 2 t idal current ellipses at C1 垂向坐标 - 0. 10 - 0. 30 - 0. 50 - 0. 70 - 0. 90 最大流发生时间 ( h 1. 36 1. 32 1. 28 1. 24 1. 19 1. 47 1. 43 1. 39 1. 35 1. 31 1. 81 1. 76 1. 71 1. 66 1. 60 1. 75 1. 70 1. 65 1. 60 1. 55 1. 37 1. 34 1. 32 1. 29 1. 26 3. 27 3. 21 3. 14 3. 07 3. 00 C2 42, 21 - 0. 10 - 0. 30 - 0. 50 - 0. 70 - 0. 90 C3 37, 22 - 0. 10 - 0. 30 - 0. 50 - 0. 70 - 0. 90 C4 39, 23 - 0. 10 - 0. 30 - 0. 50 - 0. 70 - 0. 90 C5 31, 16 - 0. 10 - 0. 30 - 0. 50 - 0. 70 - 0. 90 C6 30, 3 - 0. 10 - 0. 30 - 0. 50 - 0. 70 - 0. 90 注 : i , j 为网格坐标 ; 值在海面为 0, 在海底为 - 1; 椭率为负表示左旋 , 其余为右旋 湾内潮流椭圆椭率随深度增加减小 , 说明随着深度的增加潮流椭圆变扁。 随深度增加 , 椭圆长轴的偏转方向基本右偏, 偏转角度 3 3. 4. 2 垂直流速 潮流的垂直流 速比水平流速小得多 , 其量级为 10 - 4 m/ s , 底层垂向速度 较大, 可达 30 10- 4 m/ s。表层较小, 最大仅为 7 10- 4 m/ s, 最大值集中在团岛 黄岛间水深梯度 5 , 但也有左偏现象。 4期 孙英兰等 : 胶州湾三维变动 边界的潮流数值模拟 361 较大的水域, 说明垂直速度主要由海底起伏引起。 4 4. 1 4. 2 结论 本模型垂向采用 坐标 , 水平向采用曲线坐标, 因而较好地拟合底形及岸形。 本模型以较高的分辨率揭示了胶州湾潮流的基本特征 , 并得出与理论分析相一致的 结论。模拟结果显示 , 潮流的量值及最大流发生时刻, 随深度增加而减小, 在底层附近, 减 小比率加大, 说明底摩擦作用的结果。 4. 3 底层垂向流速大于表层 , 量级在 10- 4 m/ s 左右 , 最大值为 30 10- 4 m/ s, 发生于胶 州湾口水深水平梯度变化较大的海域。 4. 4 采用二阶紊封闭模式计算各水层的紊动粘滞系数, 克服了垂向紊动系数采用常数 , 或进行各种人为假定的弱点。 4. 5 变动边界模式包括了河口及潮滩水域, 较固定边界模式结果更切合实际。 致谢 美国乔治亚大学陈长胜博士及郑连远、 田晖老师, 在论文完成过程中给予的热情 指导和帮助, 一并致谢。 参 考 文 献 688 100 85 江文胜 , 孙文心 , 2000. 渤海悬浮颗粒物的三维输运模式 江文胜 , 孙文心 , 2001. 渤海悬浮颗粒物的三维输运模式 朱永春 , 蔡启铭 , 1998. 太湖梅梁湾三维水动力学研究 . 模式 . 海洋与湖沼 , 31( 6 : 682 8 . 模拟结果. 海洋与湖沼 , 32( 1 : 94 孙洪