常曲率弯道二维流动稳定性与非线性演化规律分析1)

2021-03-24李彬徐海珏白玉川冀自青

力学学报 2021年1期

李彬徐海珏白玉川冀自青

∗(天津大学水利工程仿真与安全国家重点实验室,天津 300350)

†(天津大学建筑工程学院,天津 300350)

引言

弯曲河流的稳定性与非线性特征自20 世纪80年代以来,众多学者都未停止对它的探究.Callander[1]在对顺直河道中的水流线性稳定性分析中认为不稳定性是导致河流弯曲或辫状的原因,并指出所有的河流是不稳定的.Ikeda 等[2]基于浅水方程构建了在弯曲河流中的水流积深方程,进行了弯道不稳定性分析时,假定河道宽度为恒定值,弯曲度(河道半宽与曲率半径的比值)远小于1,在坐标转换过程中出现了与河道振幅相关的几何非线性项,但是在稳定性分析中却移除了几何非线性和与动力非线性项,他们在弯曲河流线性稳定性分析中认为最不稳定性波长是弯曲河流自身的有限振幅波长.Parker等[3]对有限振幅弯道的非线性的处理是利用Stokes的水波摄动展开方法的修正进行非线性稳定性分析,忽略了高阶的动力非线性项,仅保留了高阶的几何非线性项,他们分析弯道的形态特征是在正弦派生曲率的基础加上三次谐波的小振幅,其成果是基于水流与河床地形的线性模型,但是河岸侵蚀公式表明几何非线性是受到近岸流速在水流非线性作用下产生高次谐波的驱动而产生,即水流非线性产生了几何非线性.Seminara 和Tubino[4]在对弯曲河流的弱非线性分析中,认为河流具有非线性动力系统的特性,利用弯曲度作为扰动参数,在较小曲率的河流中的共振状态附近进行几何弱非线性分析,揭示了响应的非唯一性以及非线性效应的抑制作用.Imran 等[5]建立了水下与地表的弯曲河道的水流的非线性模型,将曲率振幅的幂次进行展开并代入水深积分的控制方程中,在线性与非线性阶段进行求解,研究了水流与河流几何参数的影响,他指出对于Sun[6]在长期的地质时间尺度上的弯曲河流演变的模拟存在基础性缺陷,原因是忽视了水流动力中的非线性的作用.Seminara[7]指出在对弯曲河流的理论分析与室内试验,都表明弯曲河流的性质与不稳定机制有关,主要为水流自身不稳定性以及水流与可侵蚀边界的流动界面的不稳定性.Pittaluga 等[8]认为弯曲河流的线性理论说明了共振机理,然而仍未建立弯曲河流完全的非线性理论,非线性作用对水流流场的影响不可忽略.Pittaluga 和Seminara[9]提出对于弯曲河流进一步研究的方向应集中在非线性与不稳定性,这与河流自身并不是一个稳定的系统相关,并指出之前对于弯曲河流的研究主要集中于线性稳定性模型的发展上面,非线性与不稳定性的作用被严重的忽略了.徐海珏和白玉川[10]针对顺直河道建立了水流与河床的扰动共振三波与床面泥沙作用的水流非线性理论模型,分析了顺直河道沙纹动力过程与演变特征.Bai 等[11-17]认为非线性流体动力学理论是探究河流在自然或者人为扰动下的演变的重要手段,并在弱非线性理论基础上建立了常曲率与变曲率河湾中的河湾控制方程与稳定性理论,探究了在小扰动情况下的河流响应以及涡量分布,认为在外部扰动影响下,河流类型可以发生短暂的变化,但能够回到原始平衡状态类型.Ail 等[18-21]对弯曲河流的水流不稳定分析中指出,水流不稳定性与水流特性有关,其中流体黏滞系数在弯曲河流的紊流阻力的不稳定性中起到决定性作用,并且在微弯河道的水动力分析中考虑了水流的横向流动的不对称性,此外还指出河流的不稳定性源于水流与松散边界物质的相互作用,并从室内实验、现场观测及理论分析,研究了近床水流动力特征.Tubino 等[22-25]研究了河曲弯道平面发育过程中河道宽度变化和曲率之间的相互作用,不同的河道曲率对水流流场、泥沙输移和平面演变过程具有显著的贡献,并分析了河流分叉时的不稳定性.Nelson 等[26-27]构建了一维与水深积分的正交曲线坐标系下的弯曲河流计算程序,包含水流、泥沙输移和河床演变,但是在其计算程序中缺少对非线性影响的分析与讨论.

之前学者们研究主要聚焦于弯曲河流平面形态的几何非线性、床面扰动的几何非线性、水流与床面的不稳定性,但对于水流自身动力的非线性及其演化特征的研究尚有不足,本文通过构建弯道水流控制方程,分析水流自身不稳定引起的非线性作用,在水流的小扰动情况下构建水流扰动的控制方程,构造了在时间模式下的Orr-Sommerfeld 方程,分析弯道水动力的稳定性特征,并利用流动稳定性的弱非线性理论[28-29]构造考虑弯曲度的Landau-Stuart 方程,推导了在水流非线性作用下的常曲率弯道水动力的非线性演化方程,暂不考虑河岸物质组成与床面形态对水流动力的影响.在弯曲河流中,往往存在对于河流的曲率半径沿程不变,宽度受两侧岸壁限制的情况,如黑龙江第一弯、黄河晋陕峡谷段永和乾坤湾等常曲率河湾.在弯曲河流中,水流动力为泥沙的输移提供了动力,进而影响到河流形态与河床形态,因此研究常曲率弯道水动力结构以期为进一步探究河流稳定性与床面形态变化的完全非线性提供水动力非线性基础.

图1 黑龙江第一湾Fig.1 The first bend of Heilongjiang

1 理论模型

1.1 坐标变换

本文研究对象为常曲率河湾,为了反映出河湾的流动情况,首先建立适用于一般蜿蜒河湾的正交曲线坐标系,其中s为河道中轴线,指向下游流动方向,n为垂直于s方向.取B为半河宽;x0和y0分别为河道中心线的坐标值,如图2 所示.

图2 河湾正交曲线坐标系Fig.2 Orthogonal curvilinear coordinate system in meandering river

正交曲线坐标系与直角坐标系的转换关系为

曲率半径为

拉梅系数为

1.2 常曲率河湾控制方程

在二维情况下,对连续性方程与水流运动方程进行无量纲化.在正交曲线坐标下,取为半河宽;为河道最小曲率半径;为基本流的流速峰值,“∗”表示有量纲,因此长度、速度和时间尺度分别为:即有

上式中,p∗为压强,ρ 为流体密度,t∗为时间,ψ 为河流弯曲度,为河流重要的弯曲特征参数,其值越小则弯道越平缓.在常曲率情况下无量纲曲率半径R=1,故hs=1 −ψn,得到的常曲率弯曲河流无量纲控制方程为

式中Re=为雷诺数,υ 为运动黏滞系数,

图3 常曲率河湾Fig.3 Sketch illustrating the constant curve river

1.3 摄动分析

因为弯曲度ψ 为一小参数,按照流体力学中对拟序结构的处理方法,对式(4)～式(6)采用摄动方法求解,将速度与压强分解为

式中,右端第一项为明渠基本流,即ψ=0 时的水流分布,第二项为弯曲修正项,第三部分为水流扰动项,ε 与ψ 皆为小量.Seminara[7]在对接近共振时的弯曲河流的弱非线性分析中,水流速度解的形式采用微小的弯曲度的幂次展开.本研究中的解的形式包含基本流、微小弯曲度的扰动与水流自身扰动,主要考虑了水流自身扰动在弯曲河流中的非线性效应.在基本流与弯曲修正项分析中,依据Lagasse 等[30]的报告,天然河湾的ψ 值大多在0.05～0.20 之内,在微弯下由弯曲引起的影响仅考虑其一阶弯曲度的作用,忽略高阶弯曲度的影响.

将基本流与一阶弯曲修正项作为常曲率弯道基本流,与水流扰动项带到常曲率河湾控制方程中,即

将上式代入控制方程式(4)～式(6)有

1.4 一阶扰动量稳定性的线性方程

本文探究的问题是在时间模式中的水流稳定性与非线性演化过程,即ω=ωr+iωi,其中虚部ωi为与扰动的增长情况相关,即扰动振幅为a=exp(ωit),扰动幅角θ=−ωrt,其满足的关系式为

上式为常曲率弯道水流动力特征的Orr-Sommerfeld方程(简称O-S 方程,具体过程与形式参见附录A),利用O-S 方程可求出在弯曲度影响下的中性曲线分布以及相应的扰动波数与雷诺数的变化情况.

1.5 常曲率河湾时间模式弱非线性演化理论

利用Landau-Stuart 方程的弱非线性理论,在时间模式下利用a表示线性理论某一扰动模式的幅值,考虑非线性的影响,扰动幅值与扰动幅角的演化方程为

其中A0=ωia,B0=−ωr,a=θ=−ωrt,根据链式法则有

式(17)中ui为任意阶扰动速度或扰动速度分量(i=1,2,···),将式(17) 代回式(12) 和式(13),则可得到任一阶数下常曲率弯道时间模式弱非线性演化方程.Pittaluga 和Seminara[9]建立的考虑了水流非线性与不稳定性作用影响的弯曲河流地貌动力过程的模型中忽略了侧壁边界层作用,但是他们指出对弯曲河流的分析需要完全的非线性控制方程进行模拟,对于侧壁边界层的作用不可忽略,本文中相应的水流扰动的侧壁上边界条件为

本文将扰动控制方程展开至五阶,即得到了一阶至五阶的常曲率弯道时间模式下弱非线性演化下的扰动控制方程,在附录B 给出了二阶扰动控制方程的展开形式.

2 结果与讨论

2.1 稳定性理论模型验证

利用中心差分与Muller 法,可以计算出在特定的Re,α 下的ωr,ωi变化,其中ωi=0 时,即为中性状态.本次的验证过程采用ψ=0 时计算得到的临界雷诺数与中性曲线与Nishioka 等[31]所做的平面Poiseuille 流实验结果进行验证,临界雷诺数采用Orszag[32]论文中的Recr=5 772.22,α=1.020 59.本文常曲率弯道模型得出对应的临界雷诺数Recr=5772.222 2,α=1.020 59 其结果与国际上公认的理论与实验模型相一致,如图4 所示.

图4 中性曲线结果验证(Re-α)Fig.4 Verification of neutral curve(Re-α)

2.2 河湾弯曲度对中性曲线特征影响

与顺直河道不同,常曲率河湾中弯曲边界会对基本流以及相应的对河湾水流动力稳定性有影响,随着弯曲度逐渐增大,基本流分布向n＜0 方向移动,即水流向凹岸移动,使得微弯弯道逐渐向有限弯曲的弯道演变.

图6 与图7 中黑色点线为不同弯曲度下的临界雷诺数分布曲线,图7 中cr=ωr/α 为扰动波速的实数部分.ψ 对α,cr,Re的影响主要表现在ψ 的增大使中性曲线向Re增大方向迁移,中性曲线外侧区域为稳定区域,内侧为失稳区域,即随着ψ 的增大,常曲率弯道的临界失稳雷诺数增大,弯道趋于稳定.在中性曲线上对于同一个雷诺数(除临界雷诺数Recr外)存在两个中性点,即上支点与下支点,随着ψ 的增大,下支点的α,cr相互趋近,上支点α,cr相差较大,失稳区域缩窄(中性曲线内侧为失稳区、外侧为稳定区),稳定性区域增大,与文献[33]描述的在弯曲度影响下的中性曲线趋势相一致,即弯道的弯曲度的存在使得流动趋于在更大雷诺数下失稳.

图5 弯曲度ψ 对基本流的影响Fig.5 Effect of curvature ratio ψ on basic flow

图6 弯曲度对Re-α 的中性曲线影响Fig.6 Effect of curvature on neutral curve of Re-α

图7 弯曲度对Re-cr的中性曲线影响Fig.7 Effect of curvature on neutral curve of Re-cr

在弯曲度逐渐增大的过程中,临界雷诺数、扰动波数与扰动波速表现出一定的规律,在图8 中,在微弯的限制下,弯曲度从0 增大至0.20 过程中,中性状态下的Recr呈现指数型增长,α 与cr呈现单调下降的趋势.

图8 不同弯曲度下稳定性特征Fig.8 Stability characteristics under different curvature ratio

2.3 非线性幅值与幅角的演化特征分析

在利用可解条件在求得Landau 系数A,B后,设初始扰动幅a0=0.01.取图6 与图7 中的中性曲线附近的值进行分析不同弯曲度对扰动幅值与扰动幅角演化的影响.以下分析中对扰动波速的虚部取ci=−0.005,α 与cr的均取中性曲线上的值.

由弯曲度影响下的扰动幅值演化曲线可知,在同一雷诺数下,随着弯曲的增大,其a随时间的衰减速度增大,a随时间的变化率A在ψ 为0.15 时,衰减率大于ψ 为0.10 与0.00 的变化率,但当时间增大到一定程度时,其变化率A为0,相应的扰动幅值也降为0,即在初始扰动在经历了一定的时间后扰动会逐渐消失,这种效应可以被称为“非线性衰减”.扰动幅角随时间呈现单调下降的趋势,扰动幅角的衰减率B基本保持不变,在雷诺数相同时,弯曲度越大,扰动幅角的衰减率B更大,即扰动幅角随时间下降的更快,如图9 与图10 所示.

2.4 扰动速度分布与变化

2.4.1 扰动幅角对扰动速度变化分析

图9 不同弯曲度下扰动幅值与扰动幅角演化图(Re=10 000,ci=−0.005)Fig.9 Evolution diagram of disturbance amplitude and angle under different curvature ratio Re(Re=10 000,ci=−0.005)

图10 不同弯曲度下扰动幅值与扰动幅角变化率曲线图(Re=10 000,ci=−0.005)Fig.10 Curve chart of amplitude and angle change rate of disturbance under different curvature ratio(Re=10 000,ci=−0.005)

上式中,Re 为取实部,Im 为取虚部,usr,usi分别为一阶至五阶扰动流速(i=1,2,3,4,5) 累和的实部与虚部分量,unr,uni分别为一阶至五阶扰动流速(i=1,2,3,4,5)累和的实部与虚部分量.考虑到ε 为小量,故以扰动振幅a取代,令k=αs+θ,将欧拉变化式(eik=cosk+i sink)代入式(19),在此将对基本流的修正项略去,仅仅分析扰动速度的分布.

在Seminara 和Tubino[4]的分析中,弯曲河道频率高阶谐波的产生是与水流动力非线性相关,即空间上的扰动与水流的弱非线性相互作用产生高次谐波过程,这一点与文献[2-3]研究的几何非线性相区别,并且认为高阶谐波不可忽略.本文研究中各阶扰动量的高次谐波的产生是源于弯道中的基本流与水流弱非线性相互作用,这一点与以上学者的研究相区别.为了探讨弯曲度在不同幅角下对扰动速度的影响,选取ψ=0.03,0.05,0.10,0.15,Re=10 000,在初始位置s=0,分析在θ=−4π,−4π −π/2,−4π −2π/3,−5π,−6π 下各个弯曲度下的扰动速度变化特征.

在图11 中,受弯曲度与扰动幅角影响下的扰动流速中,沿水流流向的usr与usi的在横断面上都呈现“S”形分布,这与白玉川和冀自青[11]所得到的不同弯曲程度河流中扰动流速的形状函数分布相似.当弯曲度较大时,在n=−1 端附近(凹岸)usr的扰动速度,如ψ=0.15 时存在一个速度转折点,使得在凹岸近壁处的速度与在n=1 端(凸岸)的速度方向相同,但远小于凸岸的速度,在此转折点之前,即在稍微远离n=−1 处,存在一个速度的极大值,该速度与在n=1 端的速度极大值相对应,但在弯曲度逐渐增大的过程中,该速度极大值大小逐渐减小,如在θ=−4π时,ψ=0.03 的usr两个速度极大值的大小几乎相等,ψ=0.15 时,在n=−1 端的usr速度极大值约为n=1端的一半(图11(a)).随着ψ 增大,沿着河道轴线方向的usr与usi对称性逐渐减弱,即在ψ=0 时,usr与usi的两个极大值点呈现关于n=0 对称的特点,在ψ 逐渐增加到0.15 时,在n=1 端的极大值位置点不发生改变,但在n=−1 端对应的极大值点向n=0 移动(图11(a)).此外在n=1 端的usr与usi的极大值随着θ 的减小(负向增大),如当θ=−4π 与θ=−6π,该点处的usr与usi在ψ 较大时衰减较慢,即弯曲度的增大使得扰动流速对扰动幅角的反应减弱(图11(a),图11(i)).

在横断面方向的扰动速度分布中,unr与uni主要呈现“抛物线”形分布,但unr在扰动幅角在θ=−4π,−5π 时(图11(b),图11(h)),弯曲度较大,即ψ=0.10,0.15,unr的分布呈现“S”形分布,即两侧速度方向相反,但在n＜0 一侧的速度明显小于n>0 一侧,这种现象在uni的速度分布中同样出现,如θ=−4π −π/2(图11(d)).

图11 扰动幅角影响下多弯曲度的扰动速度变化特征(左侧: usr与usi;右侧: unr与uni.实线为实部,虚线为虚部)Fig.11 Characteristics of disturbance velocity variation of multi curvature ratio under the influence of disturbance amplitude(left side: usrand usi;right side: unrand uni.Solid line is real part,dotted line is imaginary part)

图11 扰动幅角影响下多弯曲度的扰动速度变化特征(左侧: usr与usi;右侧: unr与uni.实线为实部,虚线为虚部)(续)Fig.11 Characteristics of disturbance velocity variation of multi curvature ratio under the influence of disturbance amplitude(left side: usrand usi;right side: unrand uni.Solid line is real part,dotted line is imaginary part)(continued)

此外,在不同的弯曲度下,unr与uni的分布都会随弯曲度的增大而向n=1 端(凸岸)偏移,扰动流速的对称性减弱.因此弯曲度的增大,对于usr与usi来说其扰动速度分布的对称性则越弱,即在n>0 侧速度加强,在n＜0 侧速度减弱,对于unr与uni的扰动速度分布逐渐向n=1 端偏移,并且弯曲度的存在使得扰动速度对扰动幅角的变化反应减弱.

2.4.2 扰动速度空间分布变化特征

对于弯曲度对扰动速度的沿程分布的影响,主要分析θ=−4π,Re=10 000 时,ψ=0.05,0.07,0.10,在s=0,π/3,π/2,2π/3,5π/6,π 位置上usr与unr变化特征.

在一个弯道内,在入口端与出口端的usr呈现对称分布,即在入口端(s=0)正向流速出现在凸岸区域(n>0),负向流速出现在凹岸区域(n＜0),在出口端(s=π)则相反,出口端的流速大小在靠近凹岸区域略小于入口端,而在靠近凸岸区域则相反.从入口到出口,扰动流速在空间分布上具有周期性,本次分析中(Re=10 000,θ=−4π),在s=π/3 处靠近两岸usr出现双峰结构(图12(b)),在s=π/2 处usr的大小达到最小值(图12(c)),随着s的增大,usr在n>0 侧流速向负向增大,n＜0 侧向正向增大,呈现出扰动水流分离的情况,分离点在s=π/2 附近.在弯曲度较大时,usr的扰动速度曲线向n>0 方向偏斜,在n＜0侧的最大值逐渐减小,在n>0 侧的最大值则稍微加强.

unr在弯曲度逐渐增大的过程中,流速向n>0方向偏斜,其速度抛物线顶点也随弯曲度的增大而向n>0 方向移动.可以看出在不同弯曲度下unr并不发生较大的变化,但在s=5π/6 至s=π 的区域内(图12(e),图12(f)),随着弯曲度的增大,unr逐渐减小,其扰动速度曲线不断向内移动.从s=0 到s=π,unr经历了从增长至峰值再衰减,其峰值的位置为s=π/2附近,这与usr的分离点一致.

图12 不同弯曲度下的扰动速度在空间位置上的变化特征(实线:usr;虚线:unr)Fig.12 Change characteristics of the disturbance velocity in spatial position under different curvature(Solid line: usr,dotted line: unr)

图12 不同弯曲度下的扰动速度在空间位置上的变化特征(实线:usr;虚线:unr)(续)Fig.12 Change characteristics of the disturbance velocity in spatial position under different curvature(Solid line: usr,dotted line: unr)(continued)

2.4.3 扰动速度空间全域分布特征

在阐明了在不同的扰动幅角以及在不同空间位置上的扰动速度随弯曲度变化特征后,对常曲率弯道中的全域扰动流速的分布特征进行分析,取ψ=0.10,θ=−4π 至−5π −π/3 中的多个特征扰动幅角进行分析,θ=−5π−π/3 作为−4π−π/3 相差半个周期进行对比分析,Re=10 000,结果如图13,图中UXRAO 为扰动速度矢量在x方向的分量.在全域中的常曲率弯道中的扰动速度分布具有明显的特点,半个周期内,θ 在减小(负向增大)过程中,在常曲率弯道内的扰动漩涡的位置也逐渐改变,主要是在θ 从−4π 至−5π中,扰动漩涡开始出现并逐渐从s=0 的位置逐渐向s=π 移动并消失的过程,计算各个漩涡中心所在的位置,发现在一个周期内扰动幅角的大小与扰动漩涡的位置有一定的相关性,如在=θ −4π −π/3 时,相应的扰动漩涡中心所在位置为s=0.338π(图13(b)).在θ 相差半个周期时,如θ 为−4π −π/3 与−5π −π/3时(图13(b)与图13(c)),扰动速度和扰动漩涡的方向相反,即θ=−4π −π/3 时,在凸岸附近的流速方向指向下游,在凹岸附近的流速则指向上游,扰动漩涡呈现顺时针旋转,当θ=−5π −π/3 时,在凸岸附近的流速方向指向上游,在凹岸附近的流速则指向下游,扰动漩涡呈现逆时针旋转.扰动漩涡随着扰动幅角周期性变化而变化,表现在:(1)在扰动幅角减小(负向增大)过程中,即θ=−4π 至−5π −π/3,扰动漩涡的位置和方向也会出现周期性变化,扰动漩涡总是从上游产生并逐渐向下游移动;(2)随着扰动幅角的减小(负向增大),扰动流场的水流大小逐渐减小,相应的扰动漩涡的强度也会逐渐衰减并最终消亡.Seminara[7]在弯曲河流失稳的性质分析中指出,如果一个初始小扰动产生时,当时间趋于无穷大整个流场不受干扰,则这种不稳定是对流不稳定;反之,当初始小扰动同时向上游与下游两个方向传播,随着时间的增长,最终影响个水流区域,则这种不稳定是绝对不稳定.在上述分析中,usr与unr随着时间延长而衰减,扰动漩涡逐渐消亡,因此本次进行的分析属于扰动传播的对流不稳定性.

为分析ψ=0.05,0.07,0.10 时的扰动速度的全域分布特征,取θ=−4π −π/3,Re=10 000 进行分析,结果如图14 所示.在同一扰动幅角下对于不同的弯曲度,扰动漩涡中心的位置大致相同,扰动漩涡边缘的速度在弯曲度较小时大于弯曲度为较大值时,相应大的速度梯度更大.在相同弯曲度下的不同扰动幅角以及相同幅角下的不同弯曲度下的扰动速度空间分布可知,在全域中扰动流速最大值出现在扰动漩涡所在的侧壁处,并且随着扰动漩涡的移动而发生迁移.

3 结论

(1)弯曲度的增大使中性曲线向雷诺数增大方向迁移,即在微弯情况下,随着弯曲度的增大,常曲率河湾的临界失稳雷诺数增大,河湾趋于稳定.

图13 不同扰动幅角下的弯道全域扰动流速分布Fig.13 Disturbance velocity distribution in the whole bend under different disturbance angles

(2)在相同雷诺数下,随着弯曲度的增大,扰动幅值随时间的衰减速度增大,即在初始扰动在经历了一定的时间后会逐渐消失,扰动流场具有非线性衰减效应.扰动幅角随时间呈现单调下降的趋势,弯曲度越大,扰动幅角随时间下降的更快.

(3)在不同的弯曲度下,沿流向的扰动速度呈现“S”形分布.随弯曲度增大,沿着流动方向的扰动速度对称性逐渐减弱,向凸岸偏斜,弯曲度越大,扰动流速对扰动幅角的反应减弱.在垂直于流动方向的扰动速度呈现“抛物线”形分布,扰动流速unr会随弯曲度的增大而向凸岸偏移.

(4)在全域中常曲率弯道中的扰动速度分布具有明显的特点,(1)在扰动幅角减小(负向增大)过程中,扰动漩涡的位置和方向呈出现周期性变化,在一个周期内,扰动漩涡总是从上游产生并逐渐向下游移动;(2)随着扰动幅角的减小(负向增大),扰动流场的水流大小逐渐衰减,相应的扰动漩涡的强度也会逐渐衰减并最终消亡,这种衰减效应属于对流不稳定性.在相同扰动幅角下的不同弯曲度扰动漩涡中心位置大致相同,但侧壁速度在弯曲度较小的时候要大于弯曲度较大的情况.在全域中扰动流速最大值出现在扰动漩涡所在的侧壁处,并且随着扰动漩涡的移动而发生迁移.

(5)对于未来的探究应该构建三维完全非线性的河湾控制方程,应考虑泥沙输移和床面形态特征的影响,构建水流动力非线性与床面形态几何非线性相互作用模型.

附录A

将式(14) 代入式(11)～式(13) 中即可得到一阶扰动的控制方程