没有合适的资源?快使用搜索试试~ 我知道了~
河口海岸近底层潮流速分布模式扩展研究,郝嘉凌,宋志尧,在沿岸水域尤其是在浅海陆架和近岸区,潮流在水流运动中起着主导的作用,在潮流加速进程中,近底流速剖面偏离传统的对数分布。根
资源推荐
资源详情
资源评论
http://www.paper.edu.cn
河口海岸近底层潮流速分布模式扩展研究
1
郝嘉凌,宋志尧,严以新,李海泉
河海大学水文水资源与水利工程科学国家重点实验室、海洋学院 (210098)
E-mail:jialingh@126.com
摘 要:在沿岸水域尤其是在浅海陆架和近岸区,潮流在水流运动中起着主导的作用,在潮
流加速进程中,近底流速剖面偏离传统的对数分布。根据前人的研究以及笔者前一部分的研
究成果,近底流速分布的对数线性模式已经得出。由于先前采用的英国 West Solent 水道
近底 1m 之内有四层测量数据,而这类数据在实际应用时很难得到,因此,将该模式进行扩
展应用于传统的六点法(近底 0.0D,0.2D,0.4D,0.6D,0.8D,表层 1.0D)数据(D 为水
深),通过长江口和王港外西洋海域两个实例计算分析,发现对数线性扩展模式在各层数据
上都显示与实测数据良好的相关性。
关键词:近底 流速剖面 对数线性 对数线性扩展
1.前言
在沿岸水域尤其是在浅海陆架和近岸区,潮流在水流运动中起着主导的作用,近底水流
因底床摩擦而产生的紊动和边界层内水流的滞后对近底泥沙的侵蚀和输送、污染物的迁移扩
散以及底床管线冲刷等产生重要影响。因此正确预测近底水流的垂直结构及其相应的底床剪
切应力对河口海岸工程的研究具有重要意义和应用价值。
Soulsby R.L. & Dyer K.R.指出在加速潮流中近底流速剖面偏离常用的对数流速剖面
[1]
。
通过引入加速长度
**
*
uu
du dt
Λ=
,得到流速分布
0
*
0
ln
zz
uz
u
z
κγ
⎡
⎤
⎛⎞
−
=−
⎢
⎥
⎜⎟
Λ
⎝⎠
⎣
⎦
。Kuo A.Y等通过
理论分析,得出了非恒定潮流模型底部剪切力边界条件的公式。使用常规摄动法求解非恒定
流边界层方程得出邻近边界层的对数线性流速分布
[2]
。但以上二者的研究引入的加速长度和
摄动量缺乏明确的物理意义,所以笔者从流体动力学原理出发,利用Prandtl混合长度和Von
Karman自相似理论,建立河口海岸近底层潮流速分布模式
[3]
0
*
0
ln
4
ZZ
ua
u
Z
κ
Z
⎡
⎤
⎛⎞
+
=−
⎢
⎥
⎜⎟
⎝⎠
⎣
⎦
(1)
其中,
00
,
Z
zDZ z D==,z为距离底床高度,D为水深,z
0
为粗糙长度,u
*
为摩阻流速,
为卡门常数,一般取 0.4。
κ
由于在文献[3]中使用的数据近底一米之内有四层测量数据,这样的资料很难取得,所
以,考虑如何将该模式进行扩展,使之能够应用于容易获得的常规的六点法(近底 0.0D,
0.2D,0.4D,0.6D,0.8D,1.0D)数据。
2.对数线性分布的扩展
考虑到纯潮波情况表面边界条件: 0
zD
u
z
=
∂
=
∂
的特性,对数线性分布扩展后取为
1
本课题得到教育部博士点基金(编号:20020294007),国家自然科学基金(50339010),国家自然科学
基金(40476039)资助。
- 1 -
http://www.paper.edu.cn
2
0
*
0
() ln
zz
uzz
ua
DzD
κ
+
⎛⎞
=++
⎜⎟
⎝⎠
z
b
D
(2)
由
*
112
0
zD
uu
ab
zDDD
κ
=
∂
=++=
∂
,因此得到
*
1
2
u
b
κ
⎛⎞
a
=
−+
⎜
⎝⎠
⎟
。于是(2)式成为
22
0
*
0
1
() ln
22
zz
uzzz
ua
DzDD
κ
⎡⎤
+
⎛⎞ ⎛⎞
=−+−
⎢⎥
⎜⎟ ⎜⎟
⎝⎠ ⎝⎠
⎢⎥
⎣⎦
az
D
(3)
3.求解过程
对于对数线性扩展,采用枚举法结合最小二乘法
[3,4]
进行求解,取目标函数为
()
()
2
22
2
0
*
11
0
1
ln
22
nn
iiii
ii i
ii
zz z z z
ua
f
uz u a u
zDDD
κ
==
⎧⎫
⎡⎤
⎛⎞
+
⎪⎪
⎛⎞ ⎛⎞
=−= −+−−
⎢⎥
⎨⎬
⎜⎟
⎜⎟ ⎜⎟
⎝⎠ ⎝⎠
⎢⎥
⎝⎠
⎪⎪
⎣⎦
⎩⎭
∑∑
(4)
2
22
00
*
11
*0 0
2
0
1
0
11
0ln ln
22
1
ln
2
nn
ii iii
ii
n
ii
i
i
zz z z z zz z
uf
a
uk zD DD zD
zz z
u
zD
==
=
⎡⎤⎡⎤⎡⎤
++
∂
⎛⎞ ⎛⎞ ⎛⎞
=⇒ − + − −
⎢⎥⎢⎥⎢⎥
⎜⎟ ⎜⎟ ⎜⎟
∂
⎝⎠ ⎝⎠ ⎝⎠
⎢⎥⎢⎥⎢⎥
⎣⎦⎣⎦⎣⎦
⎡⎤
+
⎛⎞
=−
⎢⎥
⎜⎟
⎝⎠
⎢⎥
⎣⎦
∑∑
∑
2
1
2
i
(5)
2
22
0
*
11
0
2
1
11 1
0ln
22 2
1
2
nn
i iii ii
ii
n
ii
i
i
zz z z z z z
uf
a
ak z DDD DD
zz
u
DD
==
=
⎡ ⎤⎡⎤⎡⎤
+
∂
⎛⎞ ⎛⎞ ⎛⎞
=⇒ − − + −
⎢ ⎥⎢⎥⎢⎥
⎜⎟ ⎜⎟ ⎜⎟
∂
⎝⎠ ⎝⎠ ⎝⎠
⎢ ⎥⎢⎥⎢⎥
⎣ ⎦⎣⎦⎣
⎡⎤
⎛⎞
=−
⎢⎥
⎜⎟
⎝⎠
⎢⎥
⎣⎦
∑∑
∑
2
⎦
(6)
式中n代表垂线用来拟合的点数。矩阵求解,得出A和B,然后计算目标函数f 。计算 500 组
z
0
对应的 500 个f的最小值,得出根据近底n个点拟合出来的z
0
,u
*
。
4.实例分析
4.1 实例一:长江口 2004 年 5 月水文测验
长江口 2004 年 5 月全潮水文测验定于 2004 年 5 月 3 日-5 月 6 日进行,各垂线同步观
测 32 小时左右,从低潮前开始,保证两涨两落,满足断面潮流闭合要求
[5]
。利用观测南北
港和南北槽分流分沙断面的 6 条固定垂线NG0、BG0、NC1、NC2、NG1 和NG3 测点(其中NG0、
BG0 分别位于南港和北港断面,NC1、NC2 分别位于南槽分流分沙比观测的上下断面,NG1,NG3
位于南港下段~北槽进口段)进行垂线流速分布的分析计算。见图 1。
- 2 -
剩余8页未读,继续阅读
资源评论
weixin_38694355
- 粉丝: 3
- 资源: 964
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功