Numerical Simulation of the Influence of Riparian Vegetation on Hydrodynamic Characteristics of a Compound River
-
摘要: 广泛分布于天然河道岸坡的植被能够改变河道水流特性，进而影响坡脚淘刷、河岸稳定及河道形态演变趋势。为精细揭示岸坡植被对复式河道水动力特性的影响，将植被作为固壁边界，利用雷诺应力模型封闭紊流控制方程构建了复式河道3维水动力数学模型，并对复式河道水流特性进行了模型验证，模拟结果与实测结果吻合较好。将该模型应用于岸坡有植被的复式河道开展水流特性研究，对比种植非淹没刚性植被、种植淹没刚性植被、淹没与非淹没刚性植被交替种植和无植被情况下流速、流量、河床剪应力、雷诺应力及紊动能分布的差异。结果表明：岸坡植被加剧了滩槽间横向动量交换，引起二次流涡旋数量、强度和范围的增加，同时使得横断面流速、流量分配、河床剪应力、雷诺应力及紊动能在滩槽交界区显著减小，紊动交换在临近植被的周边区域明显增加。其中，岸坡种植非淹没刚性植被时对水流特性影响最为显著，相比无植被情况，种植非淹没刚性植被后，二次流结果显示增加了一对涡旋；岸坡植被区流速、流量分配、河床剪应力、雷诺应力及紊动能大幅减小，甚至趋近于0；而临近植被的周围区域紊动交换剧烈，无量纲化的紊动能值可达到20。此外，随着河岸坡比减小，无植被河道滩槽间动量交换逐渐衰减，而坡面种植非淹没刚性植被后，则会减弱动量交换的衰减程度。Abstract: The riparian vegetation can change the characteristics of the river flow, thereby affecting the bank erosion, bank stability and the river evolution trends. In order to clarify the influence of riparian vegetation on the hydrodynamic characteristics, the method treating the vegetation as a solid wall boundary was proposed. Using the Reynolds stress model to close the turbulence equation, a three-dimensional hydrodynamic numerical model for the compound channel flows was constructed. The flow characteristics model of the compound channel under the condition of no vegetation was verified first. The model was then applied to the study of the water flow characteristics of a compound river with vegetation on the bank slope. The differences in the distribution of velocity, flow, river bed shear stress, Reynolds stress, and turbulence kinetic energy under emerged rigid vegetation, submerged rigid vegetation, alternating emerged and submerged rigid vegetation, and no vegetation were compared. The results show that vegetation intensifies the momentum exchange between the floodplain and the main channel, causing the increase in the number, the intensity and the range of secondary current vortexes. At the same time, the cross-sectional flow velocity, flow, river bed shear stress, Reynolds stress and turbulence kinetic energy are significantly reduced in the beach-trough junction area, and the turbulence exchange was significantly increased in the area around the vegetation. The effect is most significant when emerging rigid vegetation is planted on the bank slope. Compared with the case of no vegetation, the number of secondary flow vortices increases by one pair, and the velocity, flow distribution, river bed shear stress, Reynolds stress and turbulence kinetic energy of the bank slope vegetation area are greatly reduced and even close to zero. The turbulence exchange in the surrounding area near the vegetation is intense, and the dimensionless turbulence kinetic energy value can be increased to 20. In addition, as the riparian slope ratio decreases, the momentum exchange between the floodplain and the main channel without vegetation gradually decreases. However, when non-submerged rigid vegetation is planted on slope, the attenuation degree of momentum exchange will be weakened.
-
天然河流大多是由主槽和漫滩组成的复式河道，枯水期主槽承担泄流功能，洪水期水流溢出主槽，由主槽和漫滩共同泄流[1]。而在自然河流近岸临水区，即主槽岸坡，经常会观察到乔木和灌木等植被。传统观点认为，岸坡植被不仅可以减缓水流速度，减弱河岸冲刷侵蚀，避免崩岸频发，从而稳定河势，还能起到净化水体、美化环境的作用，是一种既经济又有效的生态治河措施[2]。 然而在洪水期，岸坡植被在一定程度上会改变滩槽内部原有水流结构，同时植被和水流间的相互作用会使河道中部分水流的能量转化为植被附近产生的紊动能，从而降低河道行洪能力。因此，了解复式河道岸坡植被对水流的扰动影响，进而分析岸坡植被的利害关系，对于河道防洪、护岸工程设计及水生态保护与修复均具有重要的研究意义[3]。
关于植被水流研究，在水槽试验方面，Nepf[4]利用示踪剂显示了圆形非淹没植被阵列内的紊流形态变化。López[5]、Nezu[6]和闫静[7]等以木棒或条状薄片模拟淹没刚性植被，研究了水流的紊动特性。赵芳[8]和杨克君[9]等研究了不同植被类型条件下河道的水流特性。随着数值计算理论和计算机技术的发展，解析方法和数值模拟技术得到成功的应用，且与试验相比，易得到整个计算域的水流特性。Huai等[10]提出了预测平均流速垂向分布的分析模型，很好地模拟了含双层刚性植被河道水流的垂向分布。Wu[11]和Zhang[12]等在动量方程和k–ε模型中考虑水流对植被的阻力影响，建立了植被水流2维模型。但是，含植被的河道水流具有复杂的3维各向异性，因此，3维模型的建立非常必要。吴梦瑶等[13]将多孔介质理论应用于自由表面水流3维模型，较好地反映了有限植被群的尾流区流动特征。刘诚[14]将Wu[11]等将模型扩展到非正交曲线坐标系下建立了3维模型。以上3维模拟方法可实现植被对水流的宏观效应模拟，但对植被群内紊流特性的捕捉仍存在一定困难。槐文信[15]和Nadaoka[16]等采用大涡模拟技术建立了植被水流模型，较好地模拟了植被尾流和植被群内细部水流结构和大涡的发展过程，但当计算域较大，而细部涡流复杂区域占比相对较小时，采用大涡模拟运算量较大。
复式河道受自身断面形态特点影响，水流特性复杂，再加上河道内植被与水体间相互作用，水体流动更加复杂。然而，当前研究大多以河床或漫滩植被为研究对象，如：米云彤[17]和任俊韬[18]等考虑了河床植被对水流特性的影响；蒋北寒[19]和袁素勤[20]等探究了漫滩植被对流速分布的影响。而实际河流中，植被在复式河道岸坡广泛分布，其对河道行洪的影响尚不明确。以往研究也多以单个植被或植被斑块为研究对象，而对顺水流方向植被沿程覆盖的研究较少，且现有研究结果主要集中于植被对河道流速分布和紊流特性的影响，对二次流结构、河道流量分配及河床剪应力分布影响的研究尚不多见，而同时考虑这些影响能够更加全面认识植被对水流的扰动。因此，岸坡植被对复式河道水动力特性的影响还需进一步研究。
本文将植被作为固壁边界，并采用雷诺应力模型封闭紊流控制方程，建立3维植被河道水动力数学模型，相比于2维模型，可以获取河道横断面的水流特性；同时，与传统研究在动量方程中添加植被阻力项相比，本文方法可有效捕捉植被区域内部的紊流结构，同时又避免了使用大涡模拟时的高运算量。为了验证模型的准确性，对复式河道的水流特性进行了模拟，并将结果与Tominaga和Nezu[21]的试验数据进行比较。随后，利用该模型模拟了复式河道岸坡种植不同淹没状态的植被和不同河岸坡比条件下的河道水流结构，以揭示岸坡植被对复式河道水流扰动的影响。
1. 方 法
1.1 数学模型构建
1.1.1 控制方程
利用CFD计算流体力学软件建立基于有限体积法的3维植被河道水动力数学模型，其中紊流控制方程采用连续性方程和雷诺时均纳维–斯托克斯方程[8]，并通过雷诺应力模型对紊流模型进行封闭。
对于不可压缩均匀流体，连续性方程与雷诺时均纳维–斯托克斯方程形式如下：
$$ \frac{\partial {U}_{{i}}}{\partial {x}_{i}}=0 $$ (1) $$\begin{gathered} \quad\quad \rho {U_{{j}}}\frac{{\partial {U_i}}}{{\partial {x_j}}} = - \frac{{\partial p}}{{\partial {x_i}}} + \frac{\partial }{{\partial {x_j}}}\left( {\mu \frac{{\partial {U_i}}}{{\partial {x_j}}} - \rho \overline {{u_i}^\prime {u_j}^\prime } } \right) + {g_i} \end{gathered} $$ (2) 式（1）～（2）中：i,j = 1, 2, 3为空间指数，x1、x2、x3分别代表x、y、z方向，即横向、纵向、垂向；ρ为水的密度，取998.2 kg/m3；μ为水的动力黏滞系数，取1.002×10–3
${\text{Pa}} \cdot {\rm{s}}$ ；p为时均压力；Ui（Uj）和ui（uj）分别为在xi（xj）方向上的时均流速和瞬时流速；$ \overline {{u_i}^\prime {u_j}^\prime } $ 为雷诺应力项；gi为重力加速项。紊流封闭采用雷诺应力模型，该模型通过求解雷诺应力输运方程对方程组进行封闭[22]。雷诺应力输运方程具体形式如下：
$$ \quad\quad\frac{{\partial \left( {\overline {{u_{{i}}}^\prime {u_{{j}}}^\prime } } \right)}}{{\partial t}} + {C_{C,ij}} = {D_{T,ij}} + {D_{L,ij}} + {\tau _{ij}} + {\varphi _{ij}} - {\varepsilon _{ij}} $$ (3) 其中：
对流项为：
$$ {C_{C,ij}} = \frac{{\partial \left( {{u_k}\overline {{u_i}^\prime {u_j}^\prime } } \right)}}{{\partial {x_k}}} $$ (4) 式中：xk为垂直于xi和xj的方向；uk为xk方向的瞬时流速。
紊动扩散项可通过Daly等[23]给出的广义梯度扩散模型进行计算：
$$ {D_{T,ij}} = \frac{\partial }{{\partial {x_k}}}\left( {\frac{{{\mu _{\text{t}}}}}{{{\partial _k}}}\frac{{\partial \overline {{u_i}^\prime {u_j}^\prime } }}{{\partial {x_k}}}} \right) $$ (5) 式中：
$\; {\mu _{\text{t}}} $ 为紊动黏度，采用公式${\mu _{\text{t}}} = \rho {C_\mu }\dfrac{{{k^2}}}{\varepsilon }$ 计算，$ {C_\mu } $ 取为0.09，k为紊动能，ε为紊动耗散率；$ {\sigma _k} $ 为紊动扩散系数，取为0.82。分子黏性扩散项为：
$$ {D_{L,ij}} = \frac{\partial }{{\partial {x_k}}}\left[ {\mu \frac{\partial }{{\partial {x_k}}}\left( {\overline {{u_i}^\prime {u_j}^\prime } } \right)} \right] $$ (6) 剪应力产生项为：
$$ \begin{gathered}\quad\quad{\tau _{ij}} = - \rho \left( {\overline {{u_i}^\prime {u_k}^\prime } \frac{{\partial {u_j}}}{{\partial {x_k}}} + \overline {{u_j}^\prime {u_k}^\prime } \frac{{\partial {u_i}}}{{\partial {x_k}}}} \right) \end{gathered}$$ (7) 压力应变项为：
$$ {\varphi _{ij}} = \overline {p'\left( {\frac{{\partial {u_i}^\prime }}{{\partial {u_j}^\prime }} + \frac{{\partial {u_j}^\prime }}{{\partial {u_i}^\prime }}} \right)} $$ (8) 式中，p′为单位时间内的压力。
黏性耗散项为：
$$ {\varepsilon _{ij}} = 2\mu \overline {\frac{{\partial {u_i}^\prime }}{{\partial {x_k}}}\frac{{\partial {u_j}^\prime }}{{\partial {x_k}}}} $$ (9) 雷诺应力输运方程中，i、j 取1、2、3时，得到3个正应力方程，将3个方程相加得到紊动能的准确方程；而方程式中除了剪应项外，均含有2阶或3阶相关矩，因此，需给定边界条件使雷诺应力输运方程封闭，才能对紊动能、紊动耗散率等参数进行求解。
1.1.2 边界条件
模型计算域及边界条件如图1所示，边界条件类型主要包含速度进口边界、自由流出口边界、对称边界和壁面边界。复式河道中，水面线几乎与底坡平行，自由水面高程的起伏变化很小，在水深较大的情况下可以用对称边界条件近似代替自由水面，因此，自由水面采用刚盖假定法做近似处理。模型进口处设定为水流速度进口边界条件，并给定平均流速Uave。出口处设为自由出流边界条件，即出口处沿水流流向的速度梯度和紊动强度梯度为零。由于河道形态沿主槽中心对称分布，为了减少计算工作量，计算区域选取整个断面的一半，将主槽中心纵剖面设为对称边界。河床底部、滩地边界及刚性植被均设为壁面边界，黏性流动中壁面默认为静止无滑移边界条件，即Ui = 0。由于Tominaga等[21]研究发现相对于断面形态与水力条件，壁面糙率对水流结构的影响十分微弱[24]，因此，该模型中未考虑壁面糙率的影响。
图 1 模型计算域及边界条件示意图Fig. 1 Schematic diagram of model calculation domain and boundary conditions
下载:
全尺寸图片
1.1.3 植被固壁边界处理方法及网格划分
网格划分方法和网格质量对模拟结果有很大的影响，因此合理的网格划分可以提高模型的计算精度和计算效率。针对本文模型的特点，要保证在靠近河床和植被边壁的流速大梯度变化区域网格有较细的间距，而在主槽和漫滩中心的低速梯度区域网格有更宽的间距。此外，主槽和漫滩区水平面上划分为均匀的笛卡尔网格块，而岸坡斜面上划分为向植被圆柱体拉伸的“O”型网格块。此划分方法能够在保证捕捉植被区细部水流结构的同时减小网格数量，节省计算资源。模拟时采用了不同分辨率的网格验证了模型的收敛性和模拟精度，进而选取适当的计算网格。以岸坡种植非淹没刚性植被为例，其3维几何模型不同剖面的网格划分如图2所示。
图 2 几何模型在不同剖面的网格划分Fig. 2 Grid division of the geometric model in different sections下载:
全尺寸图片
图2中：近床面和岸坡植被区网格细密，主槽和漫滩中心网格较粗；同时，在x–z剖面和植被区水平剖面图中，x1∶x2、z1∶z2及y1∶y2应小于网格尺度比1∶5。
1.1.4 求解方法
采用基于有限体积法的计算流体力学软件建立含植被河道3维紊流数值模型，方程离散时，压力速度耦合方式选择SIMPLE算法，压力的插值格式采用体积力加权格式，动量等其余量的离散格式采用2阶迎风插值格式。当各计算参量残差小于0.000 01且监测某点的参量不随迭代次数变化时，可认为计算收敛，迭代终止。
1.2 模型验证
采用Tominaga和Nezu[21]矩形复式河道紊流试验中工况 S-2的结果对所建数学模型进行验证。该试验中，水槽长12.5 m，宽0.4 m，高0.4 m；水槽底板为铁板，边壁为有机玻璃。河道主槽水深hm为0.08 m，滩地水深hf为0.04 m，主槽宽度b为0.2 m，总宽度B为0.4 m，水力参数见表1。
表 1 验证工况相关水力参数Table 1 Related hydraulic parameters of verification conditions摩阻流速U*/(m·s–1) 最大流速Umax/(m·s–1) 平均流速Uave/(m·s–1) 0.016 4 0.389 0.349 选择河道中心横断面作为代表断面，绘制了二次流、主流速及河床剪切应力分布，如图3～5所示。图3～5中，横、纵坐标轴均使用主槽水深hm进行无量纲处理，主流速使用最大流速Umax进行无量纲处理，河床剪切应力
$ \tau $ 采用平均剪切应力$ \overline \tau $ 进行无量纲处理。
图 3 二次流分布验证Fig. 3 Secondary flow distribution verification下载:
全尺寸图片
图 4 主流速分布验证Fig. 4 Main flow velocity distribution verification下载:
全尺寸图片
图 5 河床剪切应力分布验证Fig. 5 River bed shear stress distribution verification下载:
全尺寸图片
由图3可知，试验和模拟结果在复式河道滩槽交接处均存在一对方向相反的二次流漩涡，二次流漩涡的影响范围从滩槽交界处开始延伸至自由表面，水平方向的影响范围也一致且漩涡的位置、大小和强度十分接近，故拟合效果较好。
由图4可知，试验与模拟结果规律相近，只在自由水面处略有误差，这是由于数值计算时将自由水面作为对称面边界处理，无法较真实地模拟自由表面涡，导致近水面流速误差相对较大。但数值模拟与试验结果均显示，复式河道滩槽交界处的主流速等值线向主槽一侧的自由水面凸起，滩槽交界区主流速有明显减小趋势，说明数值模拟结果仍能较好地反映断面主流速分布的主要特点。
由图5可知，复式河道河床剪切应力分布的数值模拟与试验结果基本一致，均在滩槽交界处剪应力迅速减小并达到最小值，并在交界处向主河道区和漫滩区两个方向过渡时剪切应力值迅速增加至最大值，然后再分别减小，直至边界处达到最小值。
1.3 工况设置
模型设置为复式梯形河道，如图2所示，模拟类似洪水期漫滩水流河道，设置主槽宽度为1.5 m，漫滩宽度为2.0 m，岸坡高度为0.4 m，水深0.8 m，流向长度30.0 m，河道坡降为0.1%。岸坡顺水流方向沿程布设刚性植被，植被以平行方式排布，间距为0.14 m，植被直径为20 mm。模拟工况设置了不同的植被淹没方式及不同的河岸坡比（河岸横向宽度与高度之比）进行对比，所有模拟工况中的进口平均流速Uave均取0.5 m/s，具体工况设置见表2。
表 2 数值模拟工况Table 2 Numerical simulation conditions工况编号 坡比 植被淹没方式 岸坡断面植被数量 1 1∶1.0 无 0 2 1∶1.0 非淹没 4 3 1∶1.0 淹没与非淹没交替 4 4 1∶1.0 淹没 4 5 1∶1.5 无 0 7 1∶1.5 非淹没 6 6 1∶2.0 无 0 8 1∶2.0 非淹没 8 2. 结果与讨论
2.1 流速分布与流量分配
图6为岸坡种植刚性非淹没植被时（z=0.6 m）处的河道水平面流速分布。由图6可知，左侧进口处流场为非充分发展阶段，该区域流速分布与单个圆柱绕流相似[25]，也说明了模拟方法可以很好地捕捉植被区细部水流特点，具有一定的精度。由于本文考虑植被在岸坡沿程覆盖的情况，因此重点关注水流充分发展阶段的水流结构，故下文河道断面水流特性数据也均取自充分发展阶段的断面进行分析。
图 6 水平面流速分布Fig. 6 Velocity distribution on the horizontal surface下载:
全尺寸图片
图7为8种工况下的复式河道横断面流速分布。
图 7 断面主流速分布Fig. 7 Distribution of main flow velocity in section下载:
全尺寸图片
对比图7（a）～7（d）发现：岸坡植被使得岸坡区流速显著减小，且流速减小区域与植被分布相关；其中，刚性非淹没植被所占岸坡范围最大，因此对滩槽交界区流速的减小作用最强。
图7（e）～7（h）分别为1∶1.5和1∶2.0坡比条件下，岸坡有、无植被时的断面主流速分布，与图7（a）和7（b）对比可发现：无植被情况下，主流速在滩槽交界处凸起，这与Xiao等[26]对不同坡比下复式河道水流研究得出的结果相似，这是由于此处边界凸起，引起了二次流的产生；但当岸坡植被存在时，滩槽交互区流速仍显著减小。
岸坡植被改变了河道断面流速分布，进而影响流量在主槽、河岸和漫滩3个区域的分配，图8为工况1～4的断面流量分配占比。由图8可知，在来流量相同条件下，岸坡植被的存在大大削弱了岸坡区过流量，但主槽和漫滩区过流量增加，这在一定程度上增加了河道的过流负荷。其中，河岸种植非淹没刚性植被影响最显著，与岸坡无植被相比，岸坡区流量分配由9.21%缩减到2.07%，减小了80%，使得主槽区流量分配增加了6%，漫滩区流量分配增加了13%，这些结果可为确定河道行洪负荷提供依据。
图 8 断面流量分配Fig. 8 Cross-section flow distribution下载:
全尺寸图片
2.2 二次流分布
图9给出8种工况下复式河道的二次流分布，横、纵坐标轴均使用主槽水深hm进行无量纲处理。
图 9 断面二次流分布Fig. 9 Secondary flow distribution in section下载:
全尺寸图片
图9（a）～（d）显示：同一坡比条件下，岸坡无植被时，滩槽交界区存在一对方向相反的二次流涡旋。当岸坡种植非淹没刚性植被时，二次流涡旋数量增加，滩槽交界区有两对方向相反的二次流涡旋；当岸坡交替种植淹没与非淹没刚性植被时，滩槽交界区有3个相邻且方向相反的二次流涡旋；当岸坡种植淹没刚性植被时，滩槽交界区有一对方向相反的二次流涡旋，但相比无植被时，二次流涡旋水平影响范围增加。此外，根据二次流流速矢量图例可发现，岸坡有植被时的二次流约为无植被情况的3倍，说明岸坡植被可以增加二次流涡旋的数量、强度及影响范围，其中非淹没刚性植被影响最显著。
由图9（a）、（e）及（g）可见，当河岸坡比逐渐减小时，涡旋水平方向的影响范围进一步减小，二次流的形状也逐渐模糊，这是由于坡角变缓后，主槽与滩地之间存在一个流速过渡区，滩槽间横向动量交换减弱，二次流强度也随之减弱。与图9（b）、（f）及（h）对比同样可发现，岸坡植被能够促使二次流涡旋的产生，并扩大二次流的影响范围，增加滩槽间横向动量交换，但随着坡比的减小，植被的影响会随之减弱。
2.3 河床剪应力分布
图10给出了8种工况下，对平均剪应力
$ \overline \tau $ 进行无量纲化后的河床剪应力分布。图10中，河床剪切应力$ \tau $ 可由CFD计算后直接得出，再将$ \tau $ 沿床面平均即可求得平均剪应力$ \overline \tau $ 。
图 10 河床剪应力分布Fig. 10 River bed shear stress distribution下载:
全尺寸图片
由图10可见：各工况下河床剪应力分布趋势基本一致，均在主槽和滩地中心接近均匀分布，主槽和漫滩向岸坡过渡时河床剪应力逐渐增大，进入河岸边坡区又发生突变减小；但不同工况下，主槽和滩地中心河床剪应力值及滩槽交界区河床剪应力的变化幅度有所不同。
对比工况1～4可以发现：同一坡比条件下，岸坡无植被时，主槽和漫滩河床剪应力最小，无量纲化剪应力值约为1；岸坡种植植被后，主槽和漫滩河床剪应力均有所增大，可增大至1.1～1.3。而对于岸坡区无植被时，河床剪应力最大，无量纲化剪应力值约为0.7～1.0，该结果与Xiao等[26]对坡比1∶1复式河道河床剪应力的模拟结果基本吻合；有植被时，河床剪应力大幅减小，甚至趋近于零。在主槽和漫滩向岸坡过渡区，岸坡种植非淹没刚性植被时，河床剪应力增加幅度最大，交替种植淹没与非淹没刚性植被时次之，种植淹没刚性植被时最小。这说明非淹没刚性植被能够引起强烈的滩槽动量交换，从而使得该处河床剪应力急剧增加。
对比工况1～2和5～8可以发现，在岸坡有、无植被情况下，随着坡比的减小，主槽、漫滩及滩槽交界区的河床剪应力均逐渐增大，说明随着坡比减小，二次流引起的动量交换强度变小，对边界剪切应力分布的影响减弱，河床剪应力分布向均匀化趋近。
2.4 雷诺应力分布
为了解岸坡植被对河道横断面雷诺应力的影响，重点研究了坡比为1∶1条件下4种工况的复式河道横断面雷诺应力分布，如图11所示。图11中，雷诺应力（
$ {{ - }}\overline {{u{{'}}}{v{{'}}}} $ ）由平均摩阻流速$\overline {{U}}_*$ 进行无量纲化，摩阻流速${{U}}_*$ 根据$ {U_{\text{*}}} = \sqrt {\dfrac{\tau }{\rho }} $ 计算，其中，ρ为水的密度，$ \tau $ 为边界剪切应力；再将$ {U_*} $ 沿床面平均即可求得平均摩阻流速$ \overline {{U}}_* $ 。
图 11 断面雷诺应力分布Fig. 11 Reynolds stress distribution of section下载:
全尺寸图片
由图11（a）可知，岸坡无植被时，雷诺应力在岸坡区呈凸起分布。与图11（b）～（d）对比可以发现：岸坡种植植被后，滩槽交界区雷诺应力分布与无植被时有很大差异，雷诺应力分布不再呈凸起状态，而是与岸坡植被分布相关；岸坡主槽侧存在一个紧邻植被的雷诺应力值增大区，且为正值；岸坡漫滩侧存在一个紧邻植被的雷诺应力值减小区，且为负值。这说明岸坡植被的存在显著影响了雷诺应力分布的变化，使得临近岸坡的主槽和漫滩与岸坡之间的横向动量交换更加剧烈。
2.5 紊动能分布
图12为1∶1坡比条件下4种工况的复式河道横断面紊动能分布，其中，紊动能k由平均摩阻流速
$ \overline {{U}}_* $ 进行无量纲化处理。
图 12 断面紊动能分布Fig. 12 Distribution of section turbulent kinetic energy下载:
全尺寸图片
由图12可知：无植被时，滩槽交界区紊动能呈凸起分布；而岸坡种植植被后，同雷诺应力分布相似，该区域紊动能不再呈凸起状态，而是与岸坡植被分布相关。岸坡种植非淹没刚性植被时，紊动能在岸坡植被区两侧显著增加，无量纲化的紊动能值可达到20；岸坡交替种植淹没与非淹没刚性植被时，紊动能也在岸坡植被区两侧显著增加，但幅度比前者略小；岸坡种植淹没刚性植被时，紊动能在岸坡植被区两侧及植被顶部显著增加，增加幅度也小于前者。这说明岸坡植被的存在能够增加紧邻植被周围区域的紊动交换，而植被区域内部紊动能会显著减小。米云彤等[17]对植被群周围水流特性进行测定，其植被群处断面紊动能结果与本文结果趋势相同。
3. 结 论
本文采用雷诺应力模型封闭紊流模型，建立了3维植被河道水动力数学模型，在模型验证相似的基础上，利用所建模型研究了岸坡刚性植被对水流的影响，主要研究结论为：
1）岸坡植被增加滩槽间动量交换，使得岸坡区的流速、流量及河床剪应力减小，而二次流涡旋数量、强度和范围均有所增加；此外，随着河岸坡比减小，无植被河道滩槽间动量交换逐渐衰减，但在坡面种植非淹没刚性植被后，滩槽间动量交换的衰减程度有所减小。植被区域内紊动有所减弱，而临近植被的周边区域紊动交换剧烈，且岸坡种植非淹没刚性植被时影响最显著。
2）根据断面流量分配结果，可确定岸坡植被对行洪负荷的影响，而河床剪应力结果又可与河岸土体起动建立联系，进一步为坡脚冲刷及河岸稳定提供判据，因此该模拟方法和结果有助于评估和权衡植被的利弊影响，进而制定适宜的河岸带生态护岸植被配置方案，使得植被发挥到最大的生态效益和工程安全效益。
-
图 1 模型计算域及边界条件示意图
Fig. 1 Schematic diagram of model calculation domain and boundary conditions
下载:
全尺寸图片
图 2 几何模型在不同剖面的网格划分
Fig. 2 Grid division of the geometric model in different sections
下载:
全尺寸图片
图 3 二次流分布验证
Fig. 3 Secondary flow distribution verification
下载:
全尺寸图片
图 4 主流速分布验证
Fig. 4 Main flow velocity distribution verification
下载:
全尺寸图片
图 5 河床剪切应力分布验证
Fig. 5 River bed shear stress distribution verification
下载:
全尺寸图片
图 6 水平面流速分布
Fig. 6 Velocity distribution on the horizontal surface
下载:
全尺寸图片
图 7 断面主流速分布
Fig. 7 Distribution of main flow velocity in section
下载:
全尺寸图片
图 8 断面流量分配
Fig. 8 Cross-section flow distribution
下载:
全尺寸图片
图 9 断面二次流分布
Fig. 9 Secondary flow distribution in section
下载:
全尺寸图片
图 10 河床剪应力分布
Fig. 10 River bed shear stress distribution
下载:
全尺寸图片
图 11 断面雷诺应力分布
Fig. 11 Reynolds stress distribution of section
下载:
全尺寸图片
图 12 断面紊动能分布
Fig. 12 Distribution of section turbulent kinetic energy
下载:
全尺寸图片
表 1 验证工况相关水力参数
Table 1 Related hydraulic parameters of verification conditions
摩阻流速U*/(m·s–1) 最大流速Umax/(m·s–1) 平均流速Uave/(m·s–1) 0.016 4 0.389 0.349 表 2 数值模拟工况
Table 2 Numerical simulation conditions
工况编号 坡比 植被淹没方式 岸坡断面植被数量 1 1∶1.0 无 0 2 1∶1.0 非淹没 4 3 1∶1.0 淹没与非淹没交替 4 4 1∶1.0 淹没 4 5 1∶1.5 无 0 7 1∶1.5 非淹没 6 6 1∶2.0 无 0 8 1∶2.0 非淹没 8 -
[1] 杨颖宜,李东子,马永顺,等.基于Fluent仿真软件的滩地植被作用下复式河道水动力特性的精细化数值模拟[J].大连海洋大学学报,2021,36(2):317–324. doi: 10.16535/j.cnki.dlhyxb.2020-076 Yang Yingyi,Li Dongzi,Ma Yongshun,et al.Refinement three-dimensional simulation of effect of floodplain vegetation on hydrodynamic characteristics in compound open channel based on Fluent software[J].Journal of Dalian Ocean University,2021,36(2):317–324 doi: 10.16535/j.cnki.dlhyxb.2020-076 [2] Liu Chao,Nepf H.Sediment deposition within and around a finite patch of model vegetation over a range of channel velocity[J].Water Resources Research,2016,52:600–612. doi: 10.1002/2015WR018249 [3] 黄雄合.滨河植被影响下河湾演变机理研究[D].天津:天津大学,2018. Huang Xionghe.A study on the mechanism of meandering rivers under the influence of riparian vegetation[D].Tianjin:Tianjin University,2018. [4] Nepf H M.Drag,turbulence,and diffusion in flow through emergent vegetation[J].Water Resources Research,1999,35(2):479–489. doi: 10.1029/1998wr900069 [5] López F,García M.Open-channel flow through simulated vegetation:Suspended sediment transport modeling[J].Water Resources Research,1998,34(9):2341–2352. doi: 10.1029/98wr01922 [6] Nezu I,Onitsuka K.Turbulent structures in partly vegetated open-channel flows with LDA and PIV measurements[J].Journal of Hydraulic Research,2001,39(6):629–642. doi: 10.1080/00221686.2001.9628292 [7] 闫静,陈扬,唐洪武,等.含植物明渠流动分区和特征尺度研究进展[J].中国水利水电科学研究院学报,2015,13(6):428–434. doi: 10.13244/j.cnki.jiwhr.2015.06.005 Yan Jing,Chen Yang,Tang Hongwu,et al.Advances in research on sub-division and characteristic scales of open channel flows with vegetation[J].Journal of China Institute of Water Resources and Hydropower Research,2015,13(6):428–434 doi: 10.13244/j.cnki.jiwhr.2015.06.005 [8] 赵芳.刚性植被作用下明渠水流的水动力特性研究[D].武汉:武汉大学,2017. Zhao Fang.Research on hydraulic characteristics of open channel flows with rigid vegetation[D].Wuhan:Wuhan University,2017. [9] 杨克君,刘兴年,曹叔尤,等.植被作用下的复式河槽漫滩水流紊动特性[J].水利学报,2005,36(10):1263–1268. doi: 10.3321/j.issn:0559-9350.2005.10.022 Yang Kejun,Liu Xingnian,Cao Shuyou,et al.Turbulence characteristics of overbank flow in compound river channel with vegetated floodplain[J].Journal of Hydraulic Engineering,2005,36(10):1263–1268 doi: 10.3321/j.issn:0559-9350.2005.10.022 [10] Huai Wenxin,Wang Weijie,Hu Yang,et al.Analytical model of the mean velocity distribution in an open channel with double-layered rigid vegetation[J].Advances in Water Resources,2014,69(7):106–113. doi: 10.1016/j.advwatres.2014.04.001 [11] Wu W,Wang S.Application of a depth-averaged 2-D model in river restoration[C]//World Environmental & Water Resources Congress.Reston:ASCE,2014. [12] Zhang M,Hao Z,Zhang Y,et al.Numerical simulation of solitary and random wave propagation through vegetation based on VOF method[J].Acta Oceanologica Sinica,2013,37(2):38–46. doi: 10.1007/s13131-013-0330-4 [13] 吴梦瑶,张景新.基于多孔介质模型的有限柱群绕流模拟[J].水动力学研究与进展(A辑),2019,34(4):467–474. doi: 10.16076/j.cnki.cjhd.2019.04.006 Wu Mengyao,Zhang Jingxin.Numerical simulation of flow around a finite circular array of cylinders using porous media model[J].Chinese Journal of Hydrodynamics(Series A),2019,34(4):467–474 doi: 10.16076/j.cnki.cjhd.2019.04.006 [14] 刘诚,沈永明.水生植物对水沙运动影响的三维湍流模型[J].水科学进展,2008,19(6):851–856. doi: 10.3321/j.issn:1001-6791.2008.06.014 Liu Cheng,Shen Yongming.3D turbulence model for the flow and sediment transport with aquatic vegetation[J].Advances in Water Science,2008,19(6):851–856 doi: 10.3321/j.issn:1001-6791.2008.06.014 [15] 槐文信,柳梦阳,杨中华.基于LES的圆柱尾流漩涡特性分析[J].华中科技大学学报(自然科学版),2018,46(5):95–99. doi: 10.13245/j.hust.180517 Huai Wenxin,Liu Mengyang,Yang Zhonghua.Analysis of vortex characteristics in wake of circular cylinder based on LES[J].Journal of Huazhong University of Science and Technology(Natural Science Edition),2018,46(5):95–99 doi: 10.13245/j.hust.180517 [16] Nadaoka K,Yagi H.Shallow-water turbulence modeling and horizontal large-eddy computation of river flow[J].Journal of Hydraulic Engineering,1998,124(5):493–500. doi: 10.1061/(asce)0733-9429(1998)124:5(493 [17] 米云彤,王丹,蔡暾,等.植被群河道水流紊动特性研究[J].中国农村水利水电,2021(7):107–112. doi: 10.3969/j.issn.1007-2284.2021.07.017 Mi Yuntong,Wang Dan,Cai Tun,et al.Research on the flow turbulence characteristics of vegetation channels[J].China Rural Water and Hydropower,2021(7):107–112 doi: 10.3969/j.issn.1007-2284.2021.07.017 [18] 任俊韬.刚性淹没植被对河道生态环境水流特性的影响研究[D].绵阳:西南科技大学,2021. Ren Juntao.Study on the hydraulic characteristics of the ecological environment in rivers with submerged rigid vegetation[D].Mianyang:Southwest University of Science and Technology,2021. [19] 蒋北寒,杨克君,曹叔尤,等.植被化复式河道流速分布模型的参数敏感性研究[J].四川大学学报(工程科学版),2013,45(增刊2):1–7. Jiang Beihan,Yang Kejun,Cao Shuyou,et al.Parameter sensitivity of the velocity distribution model in compound channels with vegetated floodplains[J].Journal of Sichuan University(Engineering Science Edition),2013,45(Supp2):1–7 [20] 袁素勤,雷明慧,杨克君,等.滩地一排植被作用下的复式河槽水流特性研究[J].四川大学学报(工程科学版),2015,47(增刊2):24–28. doi: 10.15961/j.jsuese.2015.s2.004 Yuan Suqin,Lei Minghui,Yang Kejun,et al.Research on characteristics of the flow in a compound channel with one-line vegetation along the floodplain[J].Journal of Sichuan University(Engineering Science Edition),2015,47(Supp2):24–28 doi: 10.15961/j.jsuese.2015.s2.004 [21] Tominaga A,Nezu I.Turbulent structure in compound open-channel flows[J].Journal of Hydraulic Engineering,1991,117(1):21–41. doi: 10.1061/(asce)0733-9429(1991)117:1(21 [22] Jing H,Guo Y,Li C,et al.Three-dimensional numerical simulation of compound meandering open channel flow by the Reynolds stress model[J].International Journal for Numerical Methods in Fluids,2009,59(8):927–943. doi: 10.1002/fld.1855 [23] Daly B J,Harlow F H.Transport equations in turbulence[J].The Physics of Fluids,1970,13(11):2634–2649. doi: 10.1063/1.1692845 [24] Janssen F,Cardenas M B,Sawyer A H,et al.A comparative experimental and multiphysics computational fluid dynamics study of coupled surface–subsurface flow in bed forms[J].Water Resources Research,2012,48(8):W08514. doi: 10.1029/2012WR011982 [25] Kasiteropoulou D,Liakopoulos A,Michalolias N,et al.Numerical modelling and analysis of turbulent flow in an open channel with submerged vegetation[J].Environmental Processes,2017,4(1):47–61. doi: 10.1007/s40710-017-0235-x [26] Xiao Yang,Wang Nairu,Liang Dongfang,et al.Flow structures in trapezoidal compound channels with different side slopes of main channel[J].International Journal of Civil Engineering,2018,16(7):823–835. doi: 10.1007/s40999-017-0212-9
