• 全国中文核心期刊
  • 中国科技核心期刊
  • 美国工程索引(EI)收录期刊

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

规则波作用下刚性挺水植物波生流试验研究

陈杰 姚震 蒋昌波 曾思益 陈克伦

陈杰, 姚震, 蒋昌波, 曾思益, 陈克伦. 规则波作用下刚性挺水植物波生流试验研究[J]. 水科学进展, 2021, 32(3): 438-448. doi: 10.14042/j.cnki.32.1309.2021.03.012
引用本文: 陈杰, 姚震, 蒋昌波, 曾思益, 陈克伦. 规则波作用下刚性挺水植物波生流试验研究[J]. 水科学进展, 2021, 32(3): 438-448. doi: 10.14042/j.cnki.32.1309.2021.03.012
CHEN Jie, YAO Zhen, JIANG Changbo, ZENG Siyi, CHEN Kelun. Experimental study of the wave-induced current for emergent rigid vegetation under regular wave conditions[J]. Advances in Water Science, 2021, 32(3): 438-448. doi: 10.14042/j.cnki.32.1309.2021.03.012
Citation: CHEN Jie, YAO Zhen, JIANG Changbo, ZENG Siyi, CHEN Kelun. Experimental study of the wave-induced current for emergent rigid vegetation under regular wave conditions[J]. Advances in Water Science, 2021, 32(3): 438-448. doi: 10.14042/j.cnki.32.1309.2021.03.012

规则波作用下刚性挺水植物波生流试验研究

doi: 10.14042/j.cnki.32.1309.2021.03.012
基金项目: 

国家自然科学基金资助项目 51839002

国家自然科学基金资助项目 51979014

详细信息
    作者简介:

    陈杰(1982—), 男, 湖南衡山人, 教授, 博士, 主要从事海岸防灾减灾的研究。E-mail: chenjie166@163.com

  • 中图分类号: TV122

Experimental study of the wave-induced current for emergent rigid vegetation under regular wave conditions

Funds: 

the National Natural Science Foundation of China 51839002

the National Natural Science Foundation of China 51979014

  • 摘要: 通过在波浪水槽中进行一系列物理模型试验研究规则波作用下刚性挺水植物波生流特性,测试了2种水深和不同入射波况条件下的波生时均流速并对比了挺水植物有茎与无茎时波生流的情况,基于假设检验分别对挺水植物的拖曳力长度尺度、植物特征参数建立模型。研究结果表明:在茎密度为60株/m2时,波生流仍以根部作用为主要因素,对沉水植物公式进行拓展,据此提出挺水植物时均流速最大值预测公式;茎部对时均流速均方根具有衰减作用,且茎部在根部及根部上方区域产生的相对衰减强度,与厄塞尔数呈抛物线函数关系;在簇状根系分布下,茎部对垂向流速结构起稳定作用,尤其是在时均流速最大值所对应的量纲一垂直位置。
  • 图  1  挺水植物排布方式(非正方形)示意

    Figure  1.  Schematic diagram of arrangement (non-square) of submerged plants

    图  2  试验布置

    Figure  2.  Experimental setup

    图  3  植物模型

    Figure  3.  Plant model

    图  4  希尔伯特法计算均方根流速(M1-工况8)

    Figure  4.  Root mean square velocity by Hilbert(Run M1-8)

    图  5  M1与M2模型Urms垂向变化(M1-工况1、M2-工况12)

    Figure  5.  Table of relevant parameters of M1 and M2 models Urms under the same incoming wave condition

    图  6  γUr变化关系

    Figure  6.  Relationship between the γ and the Ursell number

    图  7  挺水植物根、根茎模型时均流垂向分布

    Figure  7.  Vertical distribution diagram of mean flow in root and stem models of emergent plants

    图  8  量纲一垂直位置随Ur变化

    Figure  8.  Relationship between the normalized vertical position maximum current speed and the Ursell number

    图  9  时均流速最大值与流体质点经历速度差关系

    Figure  9.  Relationship between the maximum of the temporal-averaged flow velocity and the velocity differential experienced by particles that encounter the top of the root

    图  10  簇状挺水植物根、根茎模型的时均流速最大值预测模型参数修正

    Figure  10.  Parameter modification of the Umax in the root and rhizome models of rigid cluster emergent plants

    表  1  试验工况

    Table  1.   Experimental condition

    模型 工况 h/cm H/cm T/s Ld/m hv/m 模型 工况 h/cm H/cm T/s Ld/m hv/m
    根(M1) 1 35 10 1.9 0.122 0.141 根+茎(M2) 12 35 10 1.9 0.124 0.179
    2 35 10 1.6 13 35 10 1.6
    3 35 8 1.9 14 35 8 1.9
    4 35 8 1.6 15 35 8 1.6
    5 35 6 1.6 16 35 6 1.6
    6 30 10 1.9 17 30 10 1.9
    7 30 10 1.6 18 30 10 1.6
    8 30 8 1.9 19 30 8 1.9 0.118 0.17
    9 30 8 1.6 20 30 8 1.6
    10 30 8 1.3 21 30 8 1.3
    11 30 6 1.9 22 30 6 1.9
    下载: 导出CSV

    表  2  相同来波条件下M1、M2模型Urms相关参数表

    Table  2.   Table of relevant parameters of M1 and M2 models′s Urms under the same incoming wave condition

    序号 Ur Urms(M2)r/(cm·s-1) Urms(M1)r/(cm·s-1) Urms(M2)/(cm·s-1) Urms(M1)/(cm·s-1) αr/% α/% γ
    1 25.53 10.14 9.30 12.04 11.43 8.26 5.08 1.63
    2 15.63 10.31 9.66 12.07 10.91 6.25 9.66 0.65
    3 20.40 8.09 7.65 10.07 9.11 5.37 9.50 0.56
    4 12.60 8.58 7.72 10.81 9.88 10.04 8.63 1.16
    5 9.82 6.49 5.93 8.03 7.69 8.54 4.29 1.99
    6 36.46 11.25 9.43 12.53 12.02 16.13 4.05 3.98
    7 26.84 12.80 11.68 15.54 13.69 8.70 11.93 0.73
    8 29.39 9.55 8.45 11.07 9.68 11.49 12.56 0.91
    9 20.69 9.87 8.94 12.35 10.91 9.35 11.68 0.80
    10 10.90 10.41 9.28 12.44 11.90 10.81 4.33 2.50
    11 21.30 6.68 6.36 7.88 7.40 4.85 6.10 0.80
    下载: 导出CSV
  • [1] ABDOLAHPOUR M, GHISALBERTI M, MCMAHON K, et al. Material residence time in marine canopies under wave-driven flows[J]. Frontiers in Marine Science, 2020, 7: 574. doi:  10.3389/fmars.2020.00574
    [2] ZHANG Y H, TANG C H, NEPF H. Turbulent kinetic energy in submerged model canopies under oscillatory flow[J]. Water Resources Research, 2018, 54(3): 1734-1750. doi:  10.1002/2017WR021732
    [3] LUHAR M, COUTU S, INFANTES E, et al. Wave-induced velocities inside a model seagrass bed[J]. Journal of Geophysical Research: Oceans, 2010, 115(C12): C12005. doi:  10.1029/2010JC006345
    [4] ABDOLAHPOUR M, HAMBLETON M, GHISALBERTI M. The wave-driven current in coastal canopies[J]. Journal of Geophysical Research: Oceans, 2017, 122(5): 3660-3674. doi:  10.1002/2016JC012446
    [5] ABDOLAHPOUR M, GHISALBERTI M, LAVERY P, et al. Vertical mixing in coastal canopies[J]. Limnology and Oceanography, 2017, 62(1): 26-42. doi:  10.1002/lno.10368
    [6] ABDOLAHPOUR M, GHISALBERTI M, MCMAHON K, et al. The impact of flexibility on flow, turbulence, and vertical mixing in coastal canopies[J]. Limnology and Oceanography, 2018, 63(6): 2777-2792. doi:  10.1002/lno.11008
    [7] van VEELEN T J, FAIRCHILD T P, REEVE D E, et al. Experimental study on vegetation flexibility as control parameter for wave damping and velocity structure[J]. Coastal Engineering, 2020, 157: 103648. doi:  10.1016/j.coastaleng.2020.103648
    [8] CHEN H F, ZOU Q P. Eulerian-Lagrangian flow-vegetation interaction model using immersed boundary method and OpenFOAM[J]. Advances in Water Resources, 2019, 126: 176-192. doi:  10.1016/j.advwatres.2019.02.006
    [9] 陈杰, 何飞, 蒋昌波, 等. 规则波作用下刚性植物拖曳力系数实验研究[J]. 水利学报, 2017, 48(7): 846-857. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201707012.htm

    CHEN J, HE F, JIANG C B, et al. Experimental investigation on drag coefficient of rigid vegetation influenced by regular waves[J]. Journal of Hydraulic Engineering, 2017, 48(7): 846-857. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201707012.htm
    [10] 何飞, 陈杰, 蒋昌波, 等. 规则波作用下植物带波高衰减特性实验研究[J]. 海洋科学进展, 2018, 36(1): 146-158. doi:  10.3969/j.issn.1671-6647.2018.01.014

    HE F, CHEN J, JIANG C B, et al. Experimental study on wave height attenuation induced by the coastal vegetation under regular waves[J]. Advances in Marine Science, 2018, 36(1): 146-158. (in Chinese) doi:  10.3969/j.issn.1671-6647.2018.01.014
    [11] 何飞, 陈杰, 蒋昌波, 等. 考虑根茎叶影响的刚性植物消浪特性实验研究[J]. 水动力学研究与进展(A辑), 2017, 32(6): 770-778. https://www.cnki.com.cn/Article/CJFDTOTAL-SDLJ201706015.htm

    HE F, CHEN J, JIANG C B, et al. Experimental investigation on wave attenuation under the effects of rigid vegetation with root, stem and crown[J]. Chinese Journal of Hydrodynamics, 2017, 32(6): 770-778. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SDLJ201706015.htm
    [12] 陈杰, 何飞, 蒋昌波, 等. 植物消波机制的实验与理论解析研究进展[J]. 水科学进展, 2018, 29(3): 433-445. doi:  10.14042/j.cnki.32.1309.2018.03.016

    CHEN J, HE F, JIANG C B, et al. Advances in laboratory experiment, theoretical analysis of mechanism of wave attenuation by vegetation[J]. Advances in Water Science, 2018, 29(3): 433-445. (in Chinese) doi:  10.14042/j.cnki.32.1309.2018.03.016
    [13] 高宇, 林光辉. 典型红树林生态系统藻类多样性及其在生态过程中的作用[J]. 生物多样性, 2018, 26(11): 1223-1235. doi:  10.17520/biods.2018080

    GAO Y, LIN G H. Algal diversity and their importance in ecological processes in typical mangrove ecosystems[J]. Biodiversity Science, 2018, 26(11): 1223-1235. (in Chinese) doi:  10.17520/biods.2018080
    [14] 裘晟, 叶勇. 红树林湿地模型研究进展[J]. 生态学杂志, 2020, 39(4): 1330-1337. https://www.cnki.com.cn/Article/CJFDTOTAL-STXZ202004030.htm

    QIU S, YE Y. Research advance of mangrove wetlands models[J]. Chinese Journal of Ecology, 2020, 39(4): 1330-1337. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-STXZ202004030.htm
    [15] LOWE R J, KOSEFF J R, MONISMITH S G. Oscillatory flow through submerged canopies: 1: velocity structure[J]. Journal of Geophysical Research: Oceans, 2005, 110(C10): C10016. doi:  10.1029/2004JC002788
    [16] BELCHER S E, JERRAM N, HUNT J C R. Adjustment of a turbulent boundary layer to a canopy of roughness elements[J]. Journal of Fluid Mechanics, 2003, 488: 369-398. doi:  10.1017/S0022112003005019
    [17] MacDONALD R W. Modelling the mean velocity profile in the urban canopy layer[J]. Boundary-Layer Meteorology, 2000, 97(1): 25-45. doi:  10.1023/A:1002785830512
    [18] COCEAL O, BELCHER S E. A canopy model of mean winds through urban areas[J]. Quarterly Journal of the Royal Meteorological Society, 2004, 130(599): 1349-1372. doi:  10.1256/qj.03.40
    [19] CHEN J, DUAN Z H, JIANG C B, et al. Laboratory study on tsunami erosion and deposition under protection of rigid emergent vegetation[J]. Natural Hazards, 2018, 92(2): 995-1010. doi:  10.1007/s11069-018-3235-x
    [20] HE F, CHEN J, JIANG C B. Surface wave attenuation by vegetation with the stem, root and canopy[J]. Coastal Engineering, 2019, 152: 103509. doi:  10.1016/j.coastaleng.2019.103509
    [21] 李海生, 吴灿雄, 欧阳美霞, 等. 广州市南沙区红树林资源现状与保护[J]. 湿地科学, 2020, 18(2): 158-165. https://www.cnki.com.cn/Article/CJFDTOTAL-KXSD202002004.htm

    LI H S, WU C X, OUYANG M X, et al. The current status and conservation of mangrove resources in Nansha district of Guangzhou[J]. Wetland Science, 2020, 18(2): 158-165. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-KXSD202002004.htm
  • [1] 张英豪, 赖锡军, 张琳, 姚昕, 邓焕广.  风浪作用下水生植物对水流结构的影响——以太湖中两种典型沉水植物为例 . 水科学进展, 2020, 31(3): 441-449. doi: 10.14042/j.cnki.32.1309.2020.03.014
    [2] 陈洪洲, 毕春伟.  近岸斜向波浪不对称特征参数化 . 水科学进展, 2018, 29(3): 424-432. doi: 10.14042/j.cnki.32.1309.2018.03.015
    [3] 陈杰, 何飞, 蒋昌波, 彭浩, 曾思益, 龚尚鹏, 陈克伦.  植物消波机制的实验与理论解析研究进展 . 水科学进展, 2018, 29(3): 433-445. doi: 10.14042/j.cnki.32.1309.2018.03.016
    [4] 聂卫波, 马孝义, 聂坤堃, 费良军.  沟灌土壤湿润体空间矩特征参数估算模型 . 水科学进展, 2017, 28(6): 829-838. doi: 10.14042/j.cnki.32.1309.2017.06.004
    [5] 刘诚, 梁燕, 王其松, 彭石.  磨刀门河口洪季波生流及其泄洪影响 . 水科学进展, 2017, 28(5): 770-779. doi: 10.14042/j.cnki.32.1309.2017.05.015
    [6] 周建银, 邵学军, 假冬冬, 秦翠翠.  弯道出流横向流速的衰减长度 . 水科学进展, 2014, 25(3): 392-400.
    [7] 房克照, 刘忠波, 邹志利, 孙家文, 尹继伟.  波生沿岸流数值模拟 . 水科学进展, 2013, 24(2): 258-265.
    [8] 林声盼, 荆长伟, MOORE Nathan, 吴嘉平.  数字高程模型分辨率对流域地形特征参数的影响 . 水科学进展, 2012, 23(4): 457-463. doi: CNKI:32.1309.P.20120614.2158.004
    [9] 谭冰, 卢鹏, 李志军, HAASChristian, NICOLAUSMarcel.  南极威德尔海冬季冰脊形拖曳力 . 水科学进展, 2012, 23(6): 869-877. doi: CNKI: 32.1309.P.20121101.1755.003
    [10] 解鸣晓, 张玮.  近岸波生流运动三维数值模拟及验证 . 水科学进展, 2011, 22(3): 391-399.
    [11] 周宏飞, 周宝佳, 代琼.  古尔班通古特沙漠植物雾凇凝结特征 . 水科学进展, 2010, 21(1): 56-62.
    [12] 王忖, 王超.  含挺水植物和沉水植物水流紊动特性 . 水科学进展, 2010, 21(6): 816-822.
    [13] 张弛, 王义刚, 郑金海.  波生流垂向结构研究综述 . 水科学进展, 2009, 20(5): 739-746.
    [14] 刘浩, 潘伟然.  营养盐负荷对浮游植物水华影响的模型研究 . 水科学进展, 2008, 19(3): 345-351.
    [15] 李伟, 束龙仓, 李砚阁.  济宁市地下水库特征参数分析 . 水科学进展, 2007, 18(2): 282-285.
    [16] 李志军, 张占海, 卢鹏, 董西路, 程斌, 陈陟.  2003年夏季北冰洋海冰动力学特征参数 . 水科学进展, 2007, 18(2): 193-197.
    [17] 吴福生, 王文野, 姜树海.  含植物河道水动力学研究进展 . 水科学进展, 2007, 18(3): 456-461.
    [18] 周明耀, 余长洪, 钱晓晴.  基于孔隙分形维数的土壤大孔隙流水力特征参数研究 . 水科学进展, 2006, 17(4): 466-470.
    [19] 王文, 许武成.  对水文时间序列混沌特征参数估计问题的讨论 . 水科学进展, 2005, 16(4): 609-616.
    [20] 黄群, 姜加虎.  应用嵌套模式研究梅梁湖区风生流 . 水科学进展, 1998, 9(1): 65-69.
  • 加载中
图(10) / 表 (2)
计量
  • 文章访问数:  79
  • HTML全文浏览量:  19
  • PDF下载量:  20
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-08-10
  • 网络出版日期:  2021-03-11
  • 刊出日期:  2021-05-30

规则波作用下刚性挺水植物波生流试验研究

doi: 10.14042/j.cnki.32.1309.2021.03.012
    基金项目:

    国家自然科学基金资助项目 51839002

    国家自然科学基金资助项目 51979014

    作者简介:

    陈杰(1982—), 男, 湖南衡山人, 教授, 博士, 主要从事海岸防灾减灾的研究。E-mail: chenjie166@163.com

  • 中图分类号: TV122

摘要: 通过在波浪水槽中进行一系列物理模型试验研究规则波作用下刚性挺水植物波生流特性,测试了2种水深和不同入射波况条件下的波生时均流速并对比了挺水植物有茎与无茎时波生流的情况,基于假设检验分别对挺水植物的拖曳力长度尺度、植物特征参数建立模型。研究结果表明:在茎密度为60株/m2时,波生流仍以根部作用为主要因素,对沉水植物公式进行拓展,据此提出挺水植物时均流速最大值预测公式;茎部对时均流速均方根具有衰减作用,且茎部在根部及根部上方区域产生的相对衰减强度,与厄塞尔数呈抛物线函数关系;在簇状根系分布下,茎部对垂向流速结构起稳定作用,尤其是在时均流速最大值所对应的量纲一垂直位置。

English Abstract

陈杰, 姚震, 蒋昌波, 曾思益, 陈克伦. 规则波作用下刚性挺水植物波生流试验研究[J]. 水科学进展, 2021, 32(3): 438-448. doi: 10.14042/j.cnki.32.1309.2021.03.012
引用本文: 陈杰, 姚震, 蒋昌波, 曾思益, 陈克伦. 规则波作用下刚性挺水植物波生流试验研究[J]. 水科学进展, 2021, 32(3): 438-448. doi: 10.14042/j.cnki.32.1309.2021.03.012
CHEN Jie, YAO Zhen, JIANG Changbo, ZENG Siyi, CHEN Kelun. Experimental study of the wave-induced current for emergent rigid vegetation under regular wave conditions[J]. Advances in Water Science, 2021, 32(3): 438-448. doi: 10.14042/j.cnki.32.1309.2021.03.012
Citation: CHEN Jie, YAO Zhen, JIANG Changbo, ZENG Siyi, CHEN Kelun. Experimental study of the wave-induced current for emergent rigid vegetation under regular wave conditions[J]. Advances in Water Science, 2021, 32(3): 438-448. doi: 10.14042/j.cnki.32.1309.2021.03.012
  • 在水生群落中, 水动力过程会对动植物的物质循环、幼虫分散、固碳过程等产生影响, 从而影响群落的生物生产。对沿海水生植物而言, 由波生流引起的水生植物流是影响植物周围水体中溶解物质和颗粒物质在生态系统中的交换强度、维持距离以及停留时间的直接动因[1-2]。研究由波生流引起的水生植物流对于分析水生生态系统循环机制具有重要意义。

    近年来, 众多学者针对由波生流引起的水生植物流特性开展了大量的试验及理论推导工作。Luhar等[3]基于边界层流动理论提出了沉水植物冠顶波生流理论模型;Abdolahpour等[4-6]对沉水植物刚柔性以及植物密度进行试验研究, 提出了含植物特征参数的沉水植物时均流速最大值预测公式, 认为沿海沉水植物冠层与周围环境之间的物质交换可以由植物带顶部的垂直混合速率来控制, 以及研究了柔性模型与刚性模型的垂直扩散系数;van Veelen等[7]通过采用不同刚度的植被模型进行物理模型试验, 研究了柔性植被摇摆作用对流速结构的影响;Chen和Zou[8]采用欧拉-拉格朗日混合网格的方法对柔性植物进行数值模拟, 研究了植被弯曲刚度对经验阻力系数和惯性系数的影响。从目前国内外研究现状可知, 以往由波生流引起的水生植物流的研究对象大多为沉水植物, 而对挺水植物的研究更关注于消波特性方面[9-12], 对水生植物流的研究较少。

    中国沿海水生生态群落中, 红树林生态系统具有较高的生物多样性[13], 其生态系统具有强大的汇碳能力, 对缓解全球气候变暖有重大意义[14]。本文以红树林为原型, 设计具有单茎与簇状根特征的刚性挺水植物波生流试验, 并通过与Abdolahpour等[4]的沉水植物试验进行对比分析, 对刚性挺水植物的植物特征参数及拖曳力长度尺度进行假设并验证, 分析挺水植物茎部对根部波生流的影响作用, 并将沉水植物时均流速最大值预测公式拓展至挺水植物适用范围。

    • 沉水植物几何形状及分布条件决定了植物带内部的水流结构, 假设沉水植物带中每株植物对水流结构的作用效应相同, 且不考虑边壁效应, 并将这种作用进行参数化, 则可用植物特征参数lflp表示[15]:

      $$ \lambda_{\mathrm{f}}=A_{\mathrm{f}} / A_{\mathrm{T}} $$ (1)
      $$ \lambda_{\mathrm{p}}=A_{\mathrm{p}} / A_{\mathrm{T}} $$ (2)

      式中: λf表示植被群迎水阻力参数;λp表示植被群单元阻力参数;Af为单株沉水植物迎流投影面积;Ap为单株沉水植物垂向投影面积;AT为单株沉水植物下底自由表面面积(植物带所占总面积除以植物带内总株数)。若植物为矩形(正方形)排布, 且植株间距为S, 每柱沉水植物模型直径为d, 高为hr, 则该模型中植物特征参数为λf =hrd/(S+d)2λp =(πd2/4)/ (S+d)2

      对于簇状挺水植物, 本研究作以下假设: ①单株挺水植物簇状根系中的每个单根对水流结构的作用效应相同, 且不考虑边壁效应;②挺水植物在规则波作用下的时均淹没部分等于静水面时的淹没部分, 对挺水植物计算植物特征参数时只考虑静水面情况。

      基于上述假设, 并将下标r表示为根部, s表示为茎部, r+s表示为根茎整体。如图 1所示, 均匀分布(非正方形)沉水植物中茎部直径为Ds, 高为hs。单株植株具有单根数量为n, 其中单根直径为Dr, 高为hr, 单植株间距分别为LW, 单植株根部所占最大直径为D(图 1中虚线圆部分), 则该挺水植物根部植物特征参数为:

      $$ {\lambda _{{\rm{f}},{\rm{r}}}} = \frac{{n{D_{\rm{r}}}{h_{\rm{r}}}}}{{(W + D)(L + D)}} $$ (3)
      $$ {\lambda _{{\rm{P}},{\rm{r}}}} = \frac{{n{\rm{ \mathsf{ π} }}D_{\rm{r}}^2}}{{4(W + D)(L + D)}} $$ (4)

      图  1  挺水植物排布方式(非正方形)示意

      Figure 1.  Schematic diagram of arrangement (non-square) of submerged plants

      因此挺水植物根部与茎部总植物特征参数为:

      $$ {\lambda _{{\rm{f}},{\rm{r}} + {\rm{s}}}} = \frac{{n{D_{\rm{r}}}{h_{\rm{r}}} + {D_{\rm{s}}}{h_{\rm{s}}}}}{{(W + D)(L + D)}} $$ (5)
      $$ {\lambda _{{\rm{P}},{\rm{r}} + {\rm{s}}}} = \frac{{\left( {nD_{\rm{r}}^2 + D_{\rm{s}}^2} \right){\rm{ \mathsf{ π} }}}}{{4(W + D)(L + D)}} $$ (6)
    • 在沉水植物冠顶拖曳力研究中, Belcher等[16]把淹没粗糙体(沉水植物)的冠顶拖曳力与电动力学中的测试电荷产生的作用力进行类比, 得出沉水植物的冠顶拖曳力与淹没粗糙体(沉水植物)尾流区域中的流体动量衰减程度有关, 并提出了淹没粗糙体(沉水植物)的冠顶拖曳力物理模型:

      $$ f_{\mathrm{id}}(x)=\frac{\left|U_{\mathrm{i}}\right| U_{\mathrm{i}}}{L_{\mathrm{d}}} $$ (7)

      式中: fid为冠层单位体积拖曳力, fid=(fix, fiy, fiz);Ui为时空平均流速, Ui=(Ux, Uy, Uz);Ld为冠顶拖曳力长度尺度, 与Ui不相关;当在粗糙体外时fid=0.

      Belcher等[16]根据MacDonald[17]提出的平均截面拖曳力系数(Cd), 推导出了Ld模型如下:

      $$ \frac{1}{L_{\mathrm{d}}}=\frac{N_{\mathrm{e}} \frac{1}{2} \bar{C}_{\mathrm{d}} A_{\mathrm{f}}}{(1-\beta) h_{\mathrm{c}} A_{\mathrm{T}}} $$ (8)

      式中: Ne为淹没粗糙体(沉水植物)区域中粗糙个体的总株数目;Cd为拖曳力系数;β为粗糙体区域中的固体体积分数;hc为粗糙体个体的高度。采用上述植物特征参数λfλp表示, 可得:

      $$ L_{\mathrm{d}}=\frac{2 h_{\mathrm{c}}\left(1-\lambda_{\mathrm{p}}\right)}{C_{\mathrm{d}} \lambda_{\mathrm{f}}} $$ (9)

      对于簇状挺水植物, 本研究作以下假设: 当茎密度未达到一定阈值时, 挺水植物波生流机制仍以根部作用为主要因素, 且茎部作用可以类比为同一植物参数的根部作用。

      基于上述假设, 并参考Coceal和Belcher[18]在研究平均空气流对城市区域内建筑物顶部作用中, 采用体积平均法计算平均建筑物高度。对挺水植物根部与茎部淹没部分采用体积法, 建立挺水植物计算模型, 可得计算模型高度(hv):

      $$ h_{\mathrm{v}}=\frac{n h_{\mathrm{r}} D_{\mathrm{r}}^{2}+h_{\mathrm{s}} D_{\mathrm{s}}^{2}}{n D_{\mathrm{r}}^{2}+D_{\mathrm{s}}^{2}} $$ (10)

      将式(5)、式(6)和式(9)代入式(10)得挺水植物根部顶端拖曳力长度尺度:

      $$ L_{\mathrm{d}, \mathrm{r}+\mathrm{s}}=\frac{2 h_{\mathrm{v}}\left(1-\lambda_{\mathrm{p}, \mathrm{r}+\mathrm{s}}\right)}{C_{\mathrm{d}} \lambda_{\mathrm{f}, \mathrm{r}+\mathrm{s}}} $$ (11)
    • Luhar等[3]认为当波浪经过沉水植物带时产生的拖曳力导致了波浪雷诺应力的产生, 且该波浪雷诺应力以及产生的波生时均流速可以用边界层流动理论进行求解, 基于线性波理论推导出了沉水植物冠顶波生时均流模型:

      $$ {U_{{\rm{c}},{\rm{m}}}} = \sqrt {\frac{4}{{{\rm{3 \mathsf{ π} }}}}\frac{{{C_{{\rm{DW}}}}}}{{{C_{{\rm{DC}}}}}}\frac{k}{\omega }u_{{\rm{w}},{\rm{m}}}^3} $$ (12)

      式中: Uc, m为沉水植物冠层内垂向时均流的流速;CDWCDC分别为拖曳力系数的时变分量与稳态分量;k为波数;ω为角频率;uw, m为沉水植物带内部振荡流速的幅值, 通常该值可假定为无植物带时边界层流速, 且由线性波理论可解得uw, m=$\frac{{H\omega }}{{2\sinh \left( {kh} \right)}}$, 其中, H为波高;h为水深。

      但上述理论公式未考虑沉水植物带密度以及沉水植物植株排布方式等因素, 存在一定的局限性。由此Abdolahpour等[4]针对植物带密度以及植物刚柔性开展试验研究, 认为沉水植物冠层内任意高度的平均波生流强度与该高度(静水面状态下)的流体质点在振荡运动时的最大速度差(δU)成正比。根据该假设理论, 结合试验数据拟合, 提出了沉水植物时均流速最大值预测模型:

      $$ \bar{U}_{\max }=0.5 U_{\mathrm{rms}}^{\infty}\left(\frac{\delta_{\mathrm{T}}}{L_{\mathrm{d}}}\right)^{0.3} $$ (13)

      式中: umax为沉水植物沿垂向时均流速最大值;Urms为沉水植物带以上时均流速均方根流速稳定值;δT为在沉水植物冠顶处水质点作振荡运动时的垂直轨道偏移距离, 且由线性波理论可以解得:

      $$ \delta_{\mathrm{T}}=\frac{H}{2} \frac{\sinh \left(k h_{\mathrm{r}}\right)}{\sinh (k h)} $$ (14)
    • 试验在长沙理工大学波浪水槽内进行, 水槽长40.0 m, 宽0.5 m, 高0.8 m。试验布置如图 2所示, 水槽左端为一台推板式造波机, 另一端设有由多孔吸波材料构成的消能网, 以减小波浪反射。

      图  2  试验布置

      Figure 2.  Experimental setup

      在Chen等[19]、He等[20]的试验基础之上, 根据学者对海岸红树林的研究成果[21], 同时考虑试验条件, 本文采用考虑根和茎的挺水植物模型。如图 3(a)所示, 该模型以新加坡Rhizophora为植物原型, 由于该种红树林其树根根须数量、根须直径、对底面投影面积等几何因素非常复杂, 本试验仅对植物带根茎部分做概化, 并采用变态模型, 横向比尺为8∶39、竖向比尺为1 ∶ 310。刚性挺水植物茎部、根部采用有机玻璃圆管进行模拟, 有机玻璃圆管动态弹性模量(Edy)为4.35×109 N/m2、静态弹性模量(Est)为3.15×109 N/m2。如图 3(b)所示, 单株挺水植物的茎部由直径为2.0 cm、高为45.0 cm的圆管组成;根部采用簇状结构分布, 由28根圆管组成, 每根圆管直径为0.8 cm、高度为14.1 cm, 两者在本文试验波浪作用下均不会产生变形, 可以较好地模拟刚性挺水植物的根茎部分。为使植物模型固定在水槽底部, 安置有机玻璃底板进行辅助固定, 底板尺寸长×宽×高为200 cm×49.5 cm×0.8 cm, 并按图 3(c)所示平面布置方式在底板上进行穿孔。

      图  3  植物模型

      Figure 3.  Plant model

      本试验共设计2种模型, 如图 3(b)所示, M1为挺水植物根部模型, M2为挺水植物根部与茎部模型, 如图 3(c)所示, 模型各植株间均为矩形分布。考虑实际情况中的波况及试验条件限制, 试验水深(h)分别采用30.0 cm和35.0 cm。采用规则波, 入射波高(H)分别为6.0 cm、8.0 cm和10.0 cm, 入射波周期(T)分别为1.3 s、1.6 s和1.9 s。具体植物模型特征参数及试验工况如表 1所示。

      表 1  试验工况

      Table 1.  Experimental condition

      模型 工况 h/cm H/cm T/s Ld/m hv/m 模型 工况 h/cm H/cm T/s Ld/m hv/m
      根(M1) 1 35 10 1.9 0.122 0.141 根+茎(M2) 12 35 10 1.9 0.124 0.179
      2 35 10 1.6 13 35 10 1.6
      3 35 8 1.9 14 35 8 1.9
      4 35 8 1.6 15 35 8 1.6
      5 35 6 1.6 16 35 6 1.6
      6 30 10 1.9 17 30 10 1.9
      7 30 10 1.6 18 30 10 1.6
      8 30 8 1.9 19 30 8 1.9 0.118 0.17
      9 30 8 1.6 20 30 8 1.6
      10 30 8 1.3 21 30 8 1.3
      11 30 6 1.9 22 30 6 1.9

      试验开始前, 进行了大量预备试验, 确保造波机的稳定性和重复可靠性。试验时, 首先放置植物模型, 再摆放仪器设备位置, 并标定流速仪测量位置。随后对造波机进行预热, 以保证造波机性能稳定。同时开启数据采集系统, 并确保浪高仪与流速仪性能良好。调整水位达到试验水深, 待水面平静后开始造波, 测量波高与流速数据。完成试验数据测量后, 停止造波, 调整流速仪垂向位置, 待水面平静后继续该工况造波。如此直至该工况下整个植物断面流速测量完毕, 则为1个工况试验完成。完成1个工况试验后, 需待水面平静后, 再进行下1个工况造波, 重复上述步骤完成下1个工况试验。

      图 2所示, 在植物带前设置2台ULS 80D超声波浪高仪(G1、G2)测量入射波数据, 可通过式(15)计算得到入射波的厄塞尔数(Ur)。

      $$ U r=\frac{H L_{\mathrm{w}}^{2}}{h^{3}} $$ (15)

      式中: Lw为波长。

      在植物带中, 使用三维剖面流速仪(S1)沿垂直剖面测量植物带内流速数据, 通过计算均方根流速来探究在植物带内的振荡流强度变化。本试验在计算均方根流速时采用Abdolahpour等[4]提出的通过希尔伯特变换求解均方根流速的方法, 如图(4)所示, 以剖面流速仪中的1处测量点数据为例, 对波传播方向上的流速时间序列数据取10个整周期数据(如图 4(a)所示), 并对这10个整周期数据进行希尔伯特变换, 得到10个整周期内的所有相位速度点, 其中包括了平均、湍流、振荡3个分量速度(如图 4(b)所示), 再通过平均法得到相位平均振荡流速(如图 4(c)所示)。因而在任意高度上的波传播方向上均方根振荡流速Urms

      $$ {U_{{\rm{rms}}}} = \sqrt {\frac{1}{{2{\rm{ \mathsf{ π} }}}}\int_{ - {\rm{ \mathsf{ π} }}}^{\rm{ \mathsf{ π} }} {(\langle } u(\phi )\rangle - \bar u{)^2}{\rm{d}}\phi } $$ (16)

      图  4  希尔伯特法计算均方根流速(M1-工况8)

      Figure 4.  Root mean square velocity by Hilbert(Run M1-8)

      式中: u表示与时间无关的定常时均流速;〈·〉表示对多个波周期平均得到的平均相位速度。

    • 图 5给出了相同来波条件下的挺水植物根模型(M1)与根茎模型(M2)2种模型的Urms沿垂向变化情况。本文参考Abdolahpour等[4]与Lowe等[15]Urms的定义方法, 将Urms定义为h>hr时的Urms稳定段平均值;将Urmsr定义为h < hr时的Urms稳定段平均值。为使上述Urms便于量化分析, 定义Urms衰减度(α)为

      $$ \alpha=\frac{U_{\mathrm{rms}(M 1)}-U_{\mathrm{rms}(\mathrm{M} 2)}}{U_{\mathrm{rms}(\mathrm{M} 1)}} \times 100 \% $$ (17)

      图  5  M1与M2模型Urms垂向变化(M1-工况1、M2-工况12)

      Figure 5.  Table of relevant parameters of M1 and M2 models Urms under the same incoming wave condition

      αrUrmsr衰减度, αUrms衰减度。如表 2所示各参数随Ur变化情况, 结果表明当Ur变化时, αrα始终大于零, 且α平均值为7.984, 方差为11.013;αr平均值为8.679, 方差为14.917。说明挺水植物茎部对UrmsrUrms都具有衰减作用, 但在两区域的衰减作用强度不同。为对比UrmsrUrms的衰减强度, 定义相对衰减强度(γ)为

      $$ \gamma=\frac{\alpha_{\mathrm{r}}}{\alpha_{\infty}} $$ (18)

      表 2  相同来波条件下M1、M2模型Urms相关参数表

      Table 2.  Table of relevant parameters of M1 and M2 models′s Urms under the same incoming wave condition

      序号 Ur Urms(M2)r/(cm·s-1) Urms(M1)r/(cm·s-1) Urms(M2)/(cm·s-1) Urms(M1)/(cm·s-1) αr/% α/% γ
      1 25.53 10.14 9.30 12.04 11.43 8.26 5.08 1.63
      2 15.63 10.31 9.66 12.07 10.91 6.25 9.66 0.65
      3 20.40 8.09 7.65 10.07 9.11 5.37 9.50 0.56
      4 12.60 8.58 7.72 10.81 9.88 10.04 8.63 1.16
      5 9.82 6.49 5.93 8.03 7.69 8.54 4.29 1.99
      6 36.46 11.25 9.43 12.53 12.02 16.13 4.05 3.98
      7 26.84 12.80 11.68 15.54 13.69 8.70 11.93 0.73
      8 29.39 9.55 8.45 11.07 9.68 11.49 12.56 0.91
      9 20.69 9.87 8.94 12.35 10.91 9.35 11.68 0.80
      10 10.90 10.41 9.28 12.44 11.90 10.81 4.33 2.50
      11 21.30 6.68 6.36 7.88 7.40 4.85 6.10 0.80

      γ>1时, 表明挺水植物茎部对根部区域的Urmsr相对衰减度大于根部上方区域, 根部区域为主要受影响区域;当γ < 1时, 表明挺水植物茎部对根部上方的Urms相对衰减度大于根部区域, 根部上方为主要受影响区域。

      图 6给出了根部相对衰减强度(γ)随厄塞尔数(Ur)的变化情况。结果显示γ总体上随Ur增大呈现先减小后增大的趋势。综上所述, 挺水植物茎部对Urms起到了衰减作用, 且随Ur由小至大, 挺水植物茎部主要作用区域从根部区域逐渐转移至根部上方区域, 当Ur到超过某一阈值时再增大, 主要作用区域从部上方区域再次逐渐转移至根部区域。

      图  6  γUr变化关系

      Figure 6.  Relationship between the γ and the Ursell number

      γUr进行非线性回归拟合, 可得:

      $$ \gamma=0.013\ 7 U r^{2}+0.575\ 2 U r+6.540\ 1 $$ (19)

      式(19)的拟合度R2为0.832 8, 表明γUr函数关系较好, 说明挺水植物茎部对Urms的主要作用区域随着来波的非线性程度变化而变化。

    • 图 7(a)图 7(k)分别给出了在不同Ur值(如表 2中序号1至11所示Ur值)的情况下, 挺水植物根模型(M1)与根茎模型(M2)2种模型时均流垂向分布情况。可以看出当Ur在发生变化时, M1与M2模型时均流速最大值所对应的量纲一垂直位置(z-hr)/δT均在[-1, 1]的区间内, 这与Abdolahpour等[4]在对淹没刚性植物冠层进行试验时结果一致。但Abdolahpour等[4]并未考虑植株排布方式, 与本试验采用的簇状排布方式不同, 造成本试验中M1模型在不同Ur值的情况下, 出现了时均流沿垂向先增大后减小的分布情况, 但并未全部出现如Abdolahpour等[4]所提出的平均波生流坍塌现象。

      图  7  挺水植物根、根茎模型时均流垂向分布

      Figure 7.  Vertical distribution diagram of mean flow in root and stem models of emergent plants

      此外, 对比M1与M2模型时均流速最大值所对应的量纲一垂直位置随Ur的变化关系可以发现, 当Ur值逐渐增大时, M1与M2模型时均流速最大值所对应的量纲一垂直位置沿着某一平均值上下波动。如图 8所示: M1模型平均值为15.027 2, 方差为0.832 9; M2模型平均值为15.018 1, 方差为0.396 0。这表明波浪非线性程度变化会引起时均流速最大值所对应的量纲一垂直位置的波动, 且M1模型的波动要明显大于M2模型。综上所述, 在来波不同波浪非线性程度即不同水动力条件下, 挺水植物茎部对时均波生流垂向分布具有稳定作用, 其主要表现在时均流速最大值所对应的量纲一垂直位置。

      图  8  量纲一垂直位置随Ur变化

      Figure 8.  Relationship between the normalized vertical position maximum current speed and the Ursell number

    • 图 9分别给出了本试验结果与Abdolahpour等[4]试验结果的时均流速最大值与流体质点经过淹没植物层顶部所经历的流速差关系图。可以看出, 本文M1模型与M2模型试验结果与Abdolahpour等[4]试验结果一致。表明当淹没植株(根部)排布为簇状、挺水植物茎部密度为60株/m2时, 由植物淹没部分引起的波生时均流强度符合Abdolahpour等[4]提出的理论, 即任意高度的时均流强度与以围绕该高度做非规则圆周运动的粒子所经历的速度差成正比, 且以时均流速最大值为主要研究对象。因此, 当茎密度为60株/m2, 挺水植物与具有同一植物特征参数的沉水植物产生的时均流速最大值相同, 故与前文假设一致, 当挺水茎密度未达到一定阈值时, 挺水植物波生流机制仍以根部作用为主要因素, 且茎部作用可以类比为同一植物特征参数的根部作用。

      图  9  时均流速最大值与流体质点经历速度差关系

      Figure 9.  Relationship between the maximum of the temporal-averaged flow velocity and the velocity differential experienced by particles that encounter the top of the root

      将挺水植物计算模型相关参数代入Abdolahpour等[4]提出的预测模型式(13)计算预测值, 结果显示本试验结果与预测值的拟合度R2为0.766 9。考虑到本文的植株根系分布方式为簇状分布与Abdolahpour等[4]不同, 因此基于式(13)进行非线性回归拟合, 对式中的系数项进行修正, 可由式(20)表示:

      $$ \bar{U}_{\max }=0.434 U_{\mathrm{rms}}^{\infty}\left(\frac{\xi_{\mathrm{T}}}{L_{\mathrm{d}}}\right)^{0.174} $$ (20)

      图 10所示, 修正后式(20)的拟合度较高且拟合度R2为0.881。值得指出的是式(20)是本文通过对变态模型试验分析得出。所以该处理方式可能会导致式(20)用于实际植物时, 存在一定的不确定因素。

      图  10  簇状挺水植物根、根茎模型的时均流速最大值预测模型参数修正

      Figure 10.  Parameter modification of the Umax in the root and rhizome models of rigid cluster emergent plants

    • 基于物理模型试验, 探究了规则波作用下挺水植物波生流特性。试验采用根茎形式与根形式2种概化模型进行对比, 主要结论如下:

      (1) 挺水植物茎部对均方根振荡流速具有衰减作用。主要作用区域随着来波的非线性程度增大而上下移动, 先从根部区域逐渐移动至根部上方区域, 再从根部上方区域逐渐返回至根部区域。并结合本试验数据, 提出了挺水植物茎部的主要作用区域预测公式。

      (2) 挺水植物茎部对时均流垂向分布具有稳定作用。挺水植物在不同水动力环境下, 茎部可以稳定根系内的时均流速分布, 起到了保护根系及其附属生物的作用。

      (3) 与前人研究结果进行对比, 发现挺水植物产生的波生流机理与沉水植物类似, 当挺水茎密度未达到一定阈值时, 挺水植物波生流仍以根部作用为主要因素, 茎部作用可以类比为同一植物参数的根部作用。通过体积平均法建立挺水植物计算模型, 对沉水植物预测公式进行系数项修正, 提出当植株根系分布为簇状且茎密度为60株/m2时, 挺水植物时均流速最大值预测模型。

参考文献 (21)

目录

    /

    返回文章
    返回