滚波演化中聚合过程的数值模拟研究1)
NUMERICAL STUDY OF COALESCENCE IN ROLL WAVE DEVELOPMENT1)
通讯作者: 2)王晓亮, 特别副研究员, 主要研究方向: 环境流体力学和颗粒物质力学. E-mail:wangxiaoliang36@bit.edu.cn
收稿日期: 2020-12-31 接受日期: 2021-01-29 网络出版日期: 2021-05-19
| 基金资助: |
|
Received: 2020-12-31 Accepted: 2021-01-29 Online: 2021-05-19
作者简介 About authors
滚波是一种重力作用下自由液面失稳诱发的水面波动现象, 通常可分为具有相对稳定波形和波速的周期性滚波与波形和波速不断变化的不规则滚波(自然滚波). 不规则滚波的相互作用和发展演化过程十分复杂, 至今对其认识尚不成熟. 本文采用基于雷诺平均Navier-Stokes方程的立面二维数值模型, 对不规则滚波发展过程中的吸收聚合和追赶聚合现象进行了数值模拟研究. 分析了两种聚合模式的演化过程, 给出了滚波聚合过程中完整的波形、波速、速度剖面以及湍流黏性等重要信息. 结果表明滚波的聚合过程是不规则滚波演化和增长的重要机制, 在特定条件下滚波增长由自然增长模式转变为以吸收聚合和追赶聚合为主的增长模式. 滚波聚合过程中, 依次经历后波追赶、爬升、与前波合并、内部流场调制等多个步骤, 最终形成一个具有更大波长和波高的滚波. 本文发现了在3个滚波间距较近的情况下, 会发生二重聚合现象, 即后两个滚波首先聚合, 然后与前波进一步聚合形成一个新的滚波.
关键词:
Roll wave is a kind of unstable gravity driven free surface fluctuation in inclined open channel, which could be classified into two types: periodic roll wave with relatively stable waveform and wave speeds, and irregular roll wave with constantly changing waveform and wave speeds (natural roll wave). The development of irregular roll wave is a complicated process, which is studied however far from enough as compared with the regular roll waves. In this study, a numerical model based on the two-dimensional Reynolds averaged Navier-Stokes equation and the renormalization group $k$-$\varepsilon$ turbulence model was adopted. A large number of numerical simulations and statistical analyses were conducted to provide more dynamical and turbulence information about coalescence processes in irregular roll waves. The evolution processes of both absorption and overtake modes of coalescence were systematically studied to obtain the waveform, wave speeds, velocity profile and turbulence viscosity etc. Results imply that coalescence plays an important role on natural roll wave development, where under certain conditions, roll wave growth changes from regular natural growth mode to irregular mode dominated by coalescence, during which a following wave chases, climbs, and merges with the leading wave successively, along with adjustments of internal flow fields, to finally produce a new combined roll wave with larger wave length and wave height. Furthermore, it is found for the first time that when three roll waves were located closely enough, double coalescences would take place. The two tailing roll waves first coalesce to form a new following wave, which then interacts with the leading wave to produce another stronger roll wave propagating downstream.
Keywords:
本文引用格式
陈歆怡, 王晓亮, 刘青泉, 张静.
Chen Xinyi, Wang Xiaoliang, Liu Qingquan, Zhang Jing.
引言
Brock[22]最早通过实验观测发现, 当滚波的波高大小不一时, 传播较快的滚波会追赶传播较慢的滚波, 在其背面爬升, 最终聚合形成一个新的滚波, 并记录了滚波聚合前后几个位置的波高时程曲线. Kranenburg[24]忽略水流黏性和湍流作用, 采用圣维南方程模拟了滚波的聚合过程, 给出了滚波聚合前后的波形图, 并指出尽管以往的理论研究表明存在永久周期性滚波, 然而实验的结果则暗示滚波的周期性充其量是准周期. Huang[23]进一步使用带黏性项的圣维南方程对滚波的聚合现象进行了探讨, 认为滚波的聚合存在吸收聚合与追赶聚合的不同形式, 将小波赶超大波并被大波吞并的现象称为吸收聚合(absorption), 而大波追赶吞并小波的聚合现象被称为追赶聚合(overtake). 按照Huang的定义, Brock实验中展示的三组滚波聚合中[22], 有两次吸收聚合现象和一次追赶聚合现象.
Cao等[15]首先将$k$-$\varepsilon $湍流模型引入圣维南方程, 模拟了滚波的形成和演化过程. 发现即便入口处施加周期性扰动, 在一定的条件下, 也可能会形成非周期性滚波. 他们认为浅水流动具有稳定周期性发展的固有波, 其频谱由流动的水深、速度、坡度和河床粗糙度确定. 当周期足够小的初始扰动向下游传播时, 固有波可以适应这些扰动并与之发生了共振, 最终形成了周期性滚波; 而当初始扰动的周期过大时, 固有波的周期无法容纳下这些扰动, 则会发展为不规则滚波, 最终形成自然滚波.
Cao等[15]和Kranenburg[24]都发现了滚波中一个非常重要的现象, 即便是周期性扰动, 发展过程中也可能会发生聚合现象. 事实上早在Brock的实验中就发现, 滚波波高的增长主要有两种方式, 分别是自然增长和聚合增长. 自然增长往往发生在滚波发展的初期阶段; 聚合增长发生在滚波的后期演化阶段, 不同大小滚波通过聚合而形成新的滚波, 其波高和波长往往比聚合前的单个滚波更大. Brock的实验观测[22]、Kranenburg[24]对滚波聚合现象的分析, 以及Cao等[15]对自然滚波形成过程的研究均表明, 滚波的发展是一个非常复杂的过程, 聚合增长是滚波发展演化的重要形式. 首先, 滚波由初始的小扰动开始增长、加强, 经历自然增长阶段; 当自然增长发展到一定阶段后, 扰动停止增长, 形成一个较为平稳的平台期. 这一阶段往往被称为``周期性滚波''. 然而, 周期性滚波是否可以长期保持, 取决于滚波的频率与浅水流动的固有波频率是否一致. 如果一致, 那么滚波的状态将维持很长时间, 发展为所谓的永久周期性滚波; 如果不一致, 在发展过程中单个滚波可以发展出微小的波速差, 导致滚波以不同的波速传播, 进一步诱发滚波聚合, 并最终形成自然滚波.
1 数值模型
沿斜槽水流方向建立求解域, 数值模型采用雷诺平均不可压缩Navier-Stokes方程模型, 为更好地模拟包含充分发展湍流和低强度湍流并存的滚波发展和演化过程[34], 采用重整化群湍流模式, 界面捕捉采用TruVOF方法. 模型的建立和验证详见文献[34]. 本文采用入口出口边界条件, 初始小扰动以入口处水深的变化给出, 如图1所示. 入口处水深采用周期扰动, $h=h_{\rm n} +\alpha \sin [(2\pi/T)t]$, 其中$h_{\rm n}$为计算域中的初始水深, $\alpha $为入口扰动的幅值, $T$为入口扰动的周期, $t$为时间. 图1中的$v_{\rm n} $为计算域中的入口速度, $g$为重力加速度, $\theta$为斜坡与水平面的夹角. 该模型成功模拟了多组周期性滚波的发生和演化过程[34].
图1
图1
求解域及边界条件示意图
Fig.1
Schematic diagram of computational domain and boundary conditions
图2
图2
(a) 当扰动周期$T= 1.12$ s时滚波运行至$t= 50$ s时的波形图; (b) 当扰动周期$T= 0.449$ s时滚波运行至$t= 50$ s时的波形图
Fig.2
(a) Roll wave form at $t= 50$ s when the perturbation period $T= 1.12$ s; (b) roll wave form at $t = 50$ s when the perturbation period $T = 0.449$ s
图3
图3
吸收聚合滚波波形发展演化
Fig.3
Waveform evolution during the absorption process of roll wave
表3 Brock实验测量滚波聚合前后波速[22]
Table 3
![]() |
2 吸收聚合
2.1 波形演化
为研究滚波的吸收聚合现象, 本文计算了大量的算例, 并进行了分析, 表4给出了两组典型的滚波吸收聚合算例的参数.
图4给出了组次1的波形演变图, 发生了2次吸收聚合事件, 称为吸收聚合事件A和吸收聚合事件B. 当$t=26.4$ s时, 后波的波高均略小于前方滚波的波高. 在追赶过程中, 后波以略高于前方滚波的波速逐渐拉近与前方滚波的距离, 并经历追赶、爬升、吞并和融合的过程, 最终形成一个波长和波高更大的滚波向下游传播.
图4
图4
滚波吸收聚合过程中波形演化图, 左侧框图表征吸收聚合事件A, 右侧框图表征吸收聚合事件B
Fig.4
Waveform evolution during the absorption process of roll wave. The left box depicts absorption A, and the right box shows absorption B
滚波聚合前后最大波高的演化如图5所示, 聚合前, 后波的波高小于前波, 在聚合过程中, 后波逐渐爬升, 波高逐渐变大, 而前波的波高在组次1的两个例子中基本保持不变, 在组次2中略微增大, 但是幅度都低于后波, 最终后波与前波合并形成一个新的滚波继续向下游传播.
图5
图5
滚波吸收聚合过程中波高演化图
Fig.5
Maximum wave height evolution during the absorption process of roll wave
2.2 滚波波速
图6
图6
滚波吸收聚合过程中波峰位置演变
Fig.6
Evolution of wave crest position during the absorption process of roll wave
滚波在聚合过程中内部流场结构复杂, 单一的波形图并不能完整表征滚波聚合作用, 本文的立面二维模型可以给出滚波吸收聚合过程中波峰处的速度剖面分布与演变过程. 图7给出了组次2滚波吸收聚合过程中两个波峰处的速度剖面演变过程, 其中黑色点代表前波, 红色点代表后波, 蓝色点为聚合后形成的滚波.
图7
图7
组次2滚波聚合过程中滚波峰值处速度剖面随时间的变化
Fig.7
Evolution of velocity profile at the roll wave crest along with time during the absorption process of roll wave from case 2
图7表明后波追赶前波过程中速度剖面基本维持不变, 在滚波发生作用时, 后波的波高由于爬升逐步增大, 直至追上前波, 逐步融合形成一个新滚波的速度剖面. 滚波聚合初始时刻速度剖面和单个滚波速度剖面相比有一些差别, 随时间推进会稍作调整演变成单个滚波的速度剖面, 这主要是由于湍流作用稍微落后于波形的演变, 湍流对聚合滚波具有一定的调制所致, 这一点将在第2.3节给出.
2.3 组次2湍流黏性
图8给出了组次2滚波吸收聚合过程中湍流黏性的分布和演化过程. 在两个滚波聚合时, 滚波波高到达最大值, 而滚波聚合后湍流黏性分布并没有直接转变为类似单个滚波的湍流黏性分布. 刚聚合形成的滚波中, 有两个湍流黏性高能区, 其形状类似聚合前的湍流黏性分布, 随着时间的推移, 这两个湍流黏性区逐步合并调整变成单个滚波的湍流黏性分布. 表明湍流黏性的演变稍微落后于滚波波形的发展, 对滚波聚合过程起着一定的调制作用, 基于浅水波方程的滚波模型无法表征这个机制.
图8
图8
组次2滚波吸收聚合过程中湍流黏性分布随时间的演化
Fig.8
Evolution of turbulent viscosity distribution with time during the absorption process of roll wave in case 2
3 追赶聚合
3.1 波形演化
表5中给出了3组典型滚波追赶聚合算例的参数.
图9
图9
组次3滚波追赶聚合过程中波形演化图
Fig.9
Waveform evolution during the overtake process of roll wave in case 3
滚波追赶聚合前后最大波高的演化如图10所示, 与吸收聚合不同的是, 追赶聚合前后波的波高大于前波; 在后波爬升过程中, 后波的波高变化较为平缓, 在短时间内最大波高的增大较有限, 不像吸收聚合那样会经历最大波高短时间内的激增. 这是由于追赶聚合中滚波虽然也需要在追赶时爬升前方滚波, 但由于后波本身波高相较于被爬升的滚波就较高, 因此较小的变化对其影响不明显. 在被吞并的过程中, 前波的形貌没有明显的变化, 波高基本维持不变. 与吸收聚合过程相同的是, 最终所形成新滚波的波高和波长相比于聚合前的两个滚波均更大.
图10
图10
滚波追赶聚合过程中波高演化图
Fig.10
Maximum wave height evolution during the overtake process of roll wave
3.2 滚波波速
图11展示了滚波追赶聚合过程中波峰位置的位置变化. 后波和前波在运动过程中的波速基本维持不变, 由其斜率不同可以看出, 后波的波速略高于前波, 追赶聚合后形成的全新滚波也基本维持等速传播.
图11
图11
滚波追赶聚合过程中波峰位置演变
Fig.11
Evolution of wave crest position during the overtake process of roll wave
表 6 组次3中追赶聚合过程中滚波波速变化
Table 6
![]() |
图12
图12
组次3滚波追赶聚合过程中滚波峰值处速度剖面随之间的变化
Fig.12
Evolution of velocity profile at the roll wave crest along with time during the overtake process of roll wave from case 3
3.3 组次4湍流黏性
图13给出了组次4追赶聚合过程中湍流黏性的分布和演化. 滚波追赶聚合过程中, 在两个波刚刚发生融合的时刻, 滚波的波高到达最大值, 有两个湍流黏性高能区, 其形状类似聚合前的湍流黏性分布. 随着时间的推移, 来自后波的能量逐渐转移进了前方滚波中, 原先保留的两个湍流黏性高能区逐渐融合成一个高湍流黏性区. 这也说明追赶聚合中湍流作用稍微落后于波高演化, 对追赶聚合形成一定的调制作用.
图13
图13
组次4滚波追赶聚合过程中湍流黏性分布随时间的演化
Fig.13
Evolution of turbulent viscosity distribution with time during the overtake process of roll wave in case 4
3.4 二重追赶聚合
图14
图14
二重聚合滚波波峰位置
Fig.14
Evolution of wave crest positions during double coalescences
图15
图15
滚波二重聚合过程中滚波峰值处速度剖面随时间的变化
Fig.15
Evolution of velocity profile at the roll wave crest along with time during the double coalescences
4 结论
本文采用基于雷诺平均Navier-Stokes方程和重组化群$k\mbox{-}\varepsilon $湍流模式的立面二维数值模型, 系统研究了不规则滚波发展和演化过程中的聚合现象, 给出了滚波聚合过程中的波速、波形、速度剖面以及湍流黏性分布等关键信息, 结论总结如下.
(1) 滚波在特定条件下, 在足够长的空间自然增长后, 会演变成以吸收聚合和追赶聚合为主要增长模式的不规则滚波, 并进一步向自然滚波发展.
(2) 在滚波的吸收聚合和追赶聚合中, 后波经历追赶、爬升、合并以及调制等多个步骤, 最终和前波聚合成一个波长和波高更大的滚波向下游传播. 滚波聚合过程中湍流黏性演化比波形发展稍慢, 湍流作用对吸收聚合和追赶聚合起一定的调制作用.
(3) 本文还首次发现, 在多个滚波位置相互接近的情况下, 将发生多重聚合现象, 后两个滚波首先聚合形成一个后波, 该后波追赶前波并聚合形成一个新的滚波向下游传播.
本文研究, 初步分析了两类典型滚波聚合中的相互作用特征, 这一领域还有许多工作可以开展, 如滚波从规则滚波转变成不规则滚波增长的条件, 湍流在滚波聚合过程中所扮演的作用等.
参考文献
Mathematical solution of the problem of roll-waves in inclined open channels
Wave formation in laminar flow down an inclined plane
On roll waves down an open inclined channel
A new model of roll waves: comparison with Brock's experiments
Influence of relative roughness and Reynolds number on the roll-waves spatial evolution
明渠水流自由面失稳与滚波发育研究
Instabilities of free surface in open-channel flows and development of roll waves
泄水陡槽中滚波的危害及预防
The hazards and prevention of roll waves in the drain chute
Roll waves in overland flow
可冲刷坡面滚波影响因素的数值模拟研究
Factors influenced on the morphology and hydrodynamics of roll waves in overland flow over mobile bed
人工加糙床面薄层滚波流水力学特性试验
Experiments on hydraulic characteristics of roll wave for sheet flow with artificial rough bed
用一种新TVD有限体积格式模拟间断行进滚波的演化
Modeling evolution of discontinuous progressing roll waves using a novel TVD finite volume scheme
Modelling roll waves with shallow water equations and turbulent closure
坡面薄层水流滚波演变规律试验研究. [博士论文]
Experimental study on evolution law of thin layer rolling wave flow on slope surface. [PhD Thesis]
坡面薄层水流滚波水力参数影响因素
Influence factors of hydraulic parameters of roll waves in overland flow
波浪破碎的一种混合湍流模拟模式
A hybrid turbulence model for wave breaking simulation
波浪扰动下河口幂律异重流的动力场分布特性
Dynamic distribution of power-law density current in estuary under wave
非线性海洋内波的理论、模型与计算
Theory, modelling and computation of nonlinear ocean internal waves
波浪作用下直立结构物附近强湍动掺气流体运动的数值模拟
Numerical simulation of the hydrodynamic characteristics of violent aerated flows near vertical structure under wave action
Development of roll waves in open channels. [PhD Thesis]
Open channel flow instabilities: Modeling the spatial evolution of roll waves. [PhD Thesis]
On the evolution of roll waves
The flow of water in an inclined channel of rectangular section
Arrested coarsening of granular roll waves
Nonlinear evolution of small disturbances into roll waves in an inclined open channel
Roll waves on a shallow layer of mud modelled as a power-law fluid
Numerical solutions for unsteady gravity-driven flows in collapsible tubes: Evolution and roll-wave instability of a steady state
Coherent structures, self-similarity, and universal roll wave coarsening dynamics
Roll waves simulation using shallow water equations and Weighted Average Flux method
Numerical proof of stability of roll waves in the small-amplitude limit for inclined thin film flow
梯级溃坝洪水洪峰增强机制
Mechanism of peak discharge enhancement of cascade dam break floods
Numerical study of roll wave development based on Navier-Stokes equations: A two-dimensional simulation
/
| 〈 |
|
〉 |






