盐水楔垂向结构的二维斜压精细模拟计算.pdfVIP

[9] 不改变磨刀门水道小潮上溯、大潮退出的特殊规律。包芸 等通过对磨刀门水道盐水楔内部垂向盐度的精 细测量，首次发现了磨刀门水道中盐水楔内部存在典型垂向分布结构。 1 河口垂向二维斜压精细模拟 2009 年12 月，课题组对珠江河口的磨刀门水道选取了两个站点进行了8 天的野外精细测量。在后期 数据分析过程中，发现在大潮期间磨刀门水道盐水楔内部涨潮时期具有间断结构，其间断面不到半米，而 盐度跳跃值却超过10‰；落潮时期，盐水楔则是连续结构，如图所示 2 2009.12.16— 17 农历十一月初一至初二 大潮 0 水 -2 深 (m) -4 -6 14514614712 14814915015115215318 1561571581590 1601612 1621634 1641656 1661678 168 10 14 16 22 时间(h) 图1 磨刀门大潮时盐度随时间变化及涨落潮时盐度垂向分布 对于这两种不同情况的盐水楔垂向结构，在流动中会有什么样的斜压力分布，对河道中涨落潮流动将 产生什么样的影响，都值得做进一步的流体力学研究，从而更深入地掌握盐水在河道中的运动规律和动力 特性。目前河口的数值计算模型均以浅水波方程为基础，在模拟计算中垂向运动描述不够精确。为了更好 的研究盐水楔内部不同垂向结构所产生的动力影响，本文从 NS 方程出发，发展新的河口斜压精细模拟方 法。 首先考虑垂向二维流动情况。鉴于河口的流动特征，采用不可压和垂向静压假设。控制方程为： u w  0 x z p g z 盐度守恒方程 2 2 s us ws  s  s   D D t x z x x2 z z2  由于盐水中密度会随盐度变化而有所变化，故不再视为常数。将盐水密度 视为一个密度常数 和  0  局部密度扰动量 之和，即    ， 采用海水状态方程来计算。  0  在计算 p 时，根据垂向静压假设，从 到z 积分得到 p p g (z ) ，其中 为距离水平面的 