HYDRODYNAMIC INSTABILITY CHARACTERISTICS OF LAMINAR FLOW IN A MEANDERING CHANNEL WITH MOVING BOUNDARY
-
摘要: 河流形态与水动力结构息息相关,形态约束水动力结构,水动力结构则通过泥沙运动进一步塑造形态,在自然界河流中形成一对辩证互馈关系.天然河流形态形式多样,大致可以分为顺直、微弯、分叉和散乱游荡几种类型,其中微弯及多个弯曲构成的河型为河流动力演化中最重要的一环.多个弯曲构成的河型可用正弦派生曲线来描述,它也是天然河流主槽与水动力结构复杂相互作用的结果.作为探讨这一过程的力学作用机理,构建摆动槽道并研究槽道摆动与其内部流动结构的互馈关系,既是流体力学研究的热点内容,也是目前河流动力过程研究的基础内容.在此重点讨论这一互馈关系前一部分,即水流对摆动边界的响应.文中建立了随体坐标系下摆动河槽与内部水流动力响应理论模型,通过给定摆动弯曲槽道的不同特征参数,研究讨论了正弦派生型摆动边界下的槽道水流动力稳定性特征,明确了弯曲槽道摆动对其内部主流及扰动水流结构的影响,确定弯曲槽道摆动波数、摆动频率对扰动流发展影响的相应参数定量关系,得到了槽道弯曲度和摆动特征对其内部水流不同尺度扰动影响的阈值选择性范围.Abstract: Configuration of river is closely related with hydrodynamic structures of flows, for the shape of a channel influences the flow structures in it, and the flow structures also affect the developing trend of the channel through the movement of sediment, forming a pair of dialectical interactions in the river system. The natural rivers are different in configurations, which can generally be divided into such types as straight, bending, branching and wandering. Among them, the bending river or the river composed of several curved channels, the result configuration of interaction between the natural river and the complex hydrodynamic flow structure in it, become one of the most important types in the study of river dynamics. As the basis of theoretical research, the establishment of model and the study on flows within a moving channel has become the focus not only from researchers of fluid mechanics, but also from investigators of river dynamics. Therefore, this study first established a theoretic model on the flow in a meandering channel with a moving boundary by using a streamwise-transverse coordinate system. It next discussed the hydro-dynamic instability characteristics of the laminar flow within the sine-generated moving boundaries. Then it quantitatively analyzed the influences of various character parameters to the velocity distributions of main flow. Finally, it obtained the selective influences from the curvature and meandering properties to the flow structure.
-
引言
复杂边界与内部水流的相互作用既是研究河流摆动机理的基础,也是近来流体力学研究的热门课题.
平行槽道内的流动在小雷诺数时为层流,其流速分布为二次分布;然而,在复杂可动边界内流速将产生变化,此时其水动力稳定性特征也会相应变化.
目前对于复杂可动边界内流体流动特征的研究主要沿着两个方向进行:一个是弹性固壁在水流作用下的自由振动.这方面主要代表学者有Davies和Carpenter[1-2], Hell和Waters[3], Guaus等[4], Pitman和Lucey[5],其核心思想是将边壁构造成无数弹性支承的薄壁梁,对弹性边界和内部水流分别建立控制方程,研究弹性边壁对平面Poiseuille流、Taylor--Couette流水动力稳定性的影响.另一个是以水流减阻为主要目标,研究静态或动态复杂边界内层流运动的水动力稳定性和湍流拟序结构的发展.例如,Hall[6]、Thomas等[7]探讨了边壁纵向振动对平面Poiseuille流水动力稳定性特征的影响;Kuhlmann等[8]针对两壁逆向旋转的空腔,研究了其内部水流的稳定性;Hell和Waters[3]则针对弹性圆柱壳中的水流结构进行研究;Ren等[9]分析并得到了边壁纵向匀速平移对内部水流水动力稳定性的影响;Sen等[10]研究了固定边壁和可动边壁交替出现时T-S波在边界上的传播.
由于研究对象的不同,以上这些成果大多研究边壁纵向运动对内部水流结构和水动力稳定性的影响,很少有模型考虑边壁横向摆动对内部水流的影响.然而对于弯曲河流的演化机理[11],槽道横向摆动边界内水流运动规律的研究更具有实际意义[12-13].
因此,本文首先假设河槽的摆动趋势,建立了随体坐标系下摆动性槽道内水流流动的理论模型,讨论了正弦派生型摆动边界下的槽道水流动力稳定性特征,其次确定了摆动波数与摆动频率对主流流速和内部扰动发展影响的定量关系,最后分析并得到了槽道弯曲度和摆动特征对其内部水流不同尺度扰动的选择性影响.
1. 基本方程
天然河道,尤其是下游河段一般是宽浅型的 (河宽与水深比值不小于1).因此,河道在下游弯曲河段会呈现出典型的三维流动特征.然而,河道的中上游河段,若经过大峡谷,则一般会形成窄深型的河道 (河宽与水深比值小于1).此时,河道水深方向上变量的变化就可忽略不计,表现出典型的平面二维水流流动特征,这种河槽中的二维流动在Parker等[14],Ikeda等[15]、Bai和Xu[16]以及Xu和Bai [17]的文中都有专门讨论和分析.因此,本文也针对窄深型河槽建立平面二维坐标系统.
边界采用正弦派生曲线形式,其中心线表述方程为$\theta (s^{\ast })=\theta _{\rm m}\sin (\alpha_{\rm c}^{\ast}s^{\ast }-\omega_{\rm c}^{\ast }t^{\ast })$,其中,$s^{\ast }$为纵向坐标,$\theta $为轴线与$x^{\ast}$轴偏角,$\theta_{\rm m}$为曲线与$x^{\ast }$轴的最大夹角,$\alpha_{\rm c}^{\ast }$和$\omega_{\rm c}^{\ast}$分别为摆动波数和摆动角频率,$\alpha_{\rm c}^{\ast }=2\pi/M^{\ast }$,$\vert \omega_{\rm c}^{\ast }\vert=2\pi /T_{\rm c}^{\ast }$,$M^{\ast }$为曲线波长,$T_{\rm c}^{\ast }$为摆动周期.变量上标*表示有量纲的物理量.
采用如图 1中所示的坐标系,河道中心线是沿着流向且处于河道中心位置的曲线,$s^{*}$方向的其他曲线都是平行于河道中心线$y^{\ast }_{0}=y^{\ast }(x^{\ast}_{0})$,而$n^{*}$方向是处处垂直于河道中心线并在水平面内的方向,因此$s^{*}$和$n^{*}$线是处处正交的关系.这种坐标系在Parker等[14]以及Ikeda等[15]的文章中首次使用,并在Bai等[18]的文章中多次讨论.
长度尺度、速度尺度、时间尺度分别用半宽$B^{\ast }$, 零阶基本流流速峰值$U_{\rm m}^{\ast }$, $B^{\ast }/U_{\rm m}^{\ast }$无量纲化;河道曲率用河道最大曲率$c^{\ast }_{\rm m}$无量纲化,对于正弦派生曲线$c^{\ast }_{\rm m}=2\pi / M^{\ast }$.
正弦派生曲线相关物理量的无量纲参数:偏角幅值$\theta_{\rm m}$,摆动波数$\alpha _{\rm c}$,摆动角频率$\omega_{\rm c}$,其与实际物理量直接的关系式为
$$\left\{ \begin{array}{l} {\alpha _{\rm {c}}} = 2\pi {B^ * }/{M^ * }\\ {\left| \omega \right|_{\rm {c}}} = 2\pi {B^ * }/U_{\rm {m}}^ * T_{\rm {c}}^ * \\ {u_{\rm {c}}} = {M^ * }/U_{\rm {m}}^ * T_{\rm {c}}^ * \end{array} \right.$$ 式中,$u_{\rm c}$为槽道摆动的无量纲相速度,亦可写为$u_{\rm c}=\omega _{\rm c}/\alpha _{\rm c}$.
物理量无量纲化关系式为
$$\left. \begin{array}{l} ({s^ * }, {n^ * }) = {B^ * }(s, n)\\ {c^ * } = c_{\rm {m}}^ * \\ c({u^ * }, {v^ * }) = U_{\rm {m}}^ * ({u_s}, {u_n})\\ {p^ * } = {\rho ^ * }U_{\rm {m}}^{ * 2}\\ p{t^ * } = t{B^ * }/U_{\rm {m}}^ * \end{array} \right\}$$ (1) 摆动相位:$\phi_{\rm c}=\alpha_{\rm c}s-\omega_{\rm c}t$;偏角写为复数形式:$\theta=\theta_{\rm m}\exp ({\rm i}\phi_{\rm c})$;无量纲曲率函数的复数形式:$c=\theta_{, s}={\rm i}\alpha_{\rm c}\theta$;无量纲曲率幅值$\psi=B^{\ast }c^{\ast }_{\rm m }=\alpha_{\rm c}\theta_{\rm m}$.
无量纲形式的控制方程为
$$\dfrac 1{h_s} \dfrac{\partial u_s}{\partial s}+\dfrac{\partial u_n}{\partial n}-\psi \dfrac{cu_n}{h_s}=0$$ (2) $$\dfrac{\partial u_s}{\partial t}+\dfrac{u_s}{h_s}\dfrac{\partial u_s}{\partial s}+u_n \dfrac{\partial u_s}{\partial n}-\psi \dfrac{cu_su_n}{h_s}= \\ \qquad -\dfrac 1{h_s} \dfrac{\partial p}{\partial s} +\dfrac 1R \bar \Delta u_s - \\ \qquad \psi R^{-1} \Big ( 2\dfrac c{h^2_s}\dfrac{\partial u_n}{\partial s}+\dfrac{\partial u_s}{\partial n}+\psi \dfrac{c^2u_s}{h^2_s} \Big ) \\ \qquad \dfrac{\psi}{R} \dfrac{\partial c}{\partial s} \Big ( \dfrac n{h^3_s} \dfrac{\partial u_s}{\partial s}-\dfrac{u_n}{h^2_s} \Big ) $$ (3) $$\frac{{\partial {u_n}}}{{\partial t}} + \frac{{{u_s}}}{{{h_s}}}\frac{{\partial {u_n}}}{{\partial s}} + {u_n}\frac{{\partial {u_n}}}{{\partial n}} + \psi \frac{{cu_s^2}}{{{h_s}}} = \\ \qquad -\dfrac{\partial p}{\partial n}+\dfrac 1R \bar \Delta u_n + \\ \qquad \dfrac \psi R \Big ( 2 \dfrac c{h^2_s} \dfrac{\partial u_s}{\partial s}- \dfrac c{h_s}\dfrac{\partial u_n}{\partial n} -\psi \dfrac{c^2u_n}{h^2_s}\Big )+ \\ \qquad \dfrac\psi R \dfrac{\partial c}{\partial s} \Big ( \dfrac n{h^3_s} \dfrac{\partial u_n}{\partial s}-\dfrac{u_s}{h^2_s} \Big)$$ (4) 无滑移边界条件
$$n=\pm 1, \ \ u_s=u_{s0}=-n\theta_{, t}, \ \ u_n=u_{n0}$$ 式中,$\bar\Delta=\dfrac 1{h^2_s} \dfrac{\partial ^2}{\partial s^2}+ \dfrac{\partial ^2}{\partial n^2}$为拉普拉斯算子; $h_s=1-\psi nc$为Lame系数; $u_s$和$u_n$分别为s方向和n方向无量纲流速;槽道中心线在曲线坐标系中的运动速度 ($u_{s0}$, $u_{n0})$,可表述为以下形式
$$n=0, \ \ \dfrac{\partial u_{s0}}{\partial s}=-u_{n0}\dfrac{\partial \theta}{\partial s}, \ \ \dfrac{\partial u_{n0}}{\partial s}=\dfrac{\partial \theta}{\partial t} \dfrac{\partial S}{\partial s}$$ 式中,$S (s, t)$为槽道中心线在拉格朗日坐标.槽道中心线伸缩比:$S_{, s}=1+\int^t_0 u_{s0, s} {\rm{d}} t$;当槽道中心线无伸缩时,$ \partial S/\partial s=1$,亦即$u_{s0, s}=0$或$u_{s0}=u_{s0}(t)$.
值得注意的是,河道形态参数很多具有复数形式,其中包括:河道波数$\alpha_{\rm c}$、河道频率$\omega_{\rm c}$、摆动相位$\phi_{\rm c}$、无量纲曲率函数c、Lame系数$h_{s}$等.基本所有参数的实部具有物理含义,而虚部没有具体物理含义,因此参数的取值一般为参数的实部.然而,在计算中需要用到参数的实部和虚部.如:偏角$\theta $可以写作$\theta=\theta_{\rm m }\exp ({\rm i}\phi_{\rm c})=\theta_{\rm m }\exp[{\rm i} (\phi_{\rm cr}+{\rm i} \phi_{\rm ci})]=\theta_{\rm m }\exp (-\phi_{\rm ci}) \exp ({\rm i}\phi_{\rm cr})$.其中,$\phi_{\rm cr}$和$\phi_{\rm ci}$分别是摆动相位$\phi_{\rm c}$的实部和虚部.因此,摆动相位实部$\phi_{\rm cr}$体现了河道周期性变化对偏角的影响,而其虚部$\phi_{\rm ci}$则体现了河道非周期摆动幅度改变对偏角的影响.
2. 基本流摄动解
摄动法又称小参数展开法,就是将系统视为理想模型的参数或结构作了微小扰动的结果来研究其运动过程的数学方法.这种方法最早应用于天体力学,用来计算小天体对大天体运动的影响,后来广泛应用于物理学和力学的理论研究.利用摄动法求解方程,也称摄动分解,通常需要在无量纲方程中选择一个能反映物理特征的无量纲小参数作为摄动量,然后假设解可以按小参数展成幂级数,将这一形式级数代入无量纲方程后,可得各级近似方程,依据这些方程可确定幂级数的系数,对级数进行截断,便得到原方程的渐近解.
无量纲曲率幅值$\psi=B^{\ast }/r^{\ast } $,表示半河宽与最小曲率半径之比,它代表了槽道的空间弯曲程度,是重要的弯曲特征参数,这个幅值越大就证明弯道越尖锐.根据实际测量资料的结果显示,天然形成河流槽道的$\psi $值大多在0.05$\sim$0.1的范围内[19-21],几乎没有超过0.3的.因此,在对方程的处理中,一般认为$\psi $是可用于方程摄动展开的小参数量.
摆动角频率$\omega_{\rm c}$反映摆动槽道的时间特征,中心线横向速度$u_{n0}$的量级为$u_{\rm c}\theta_{m}$,槽道缓慢摆动时可认为其与$\psi $同量级.
2.1 基本流
以无量纲曲率幅值$\psi $为小参数,对${\pmb X}_{\psi }$进行摄动分解,只保留一阶摄动一次谐波量,可得
$${\pmb X}_{\psi}={\pmb X}_{\psi 0}+\psi {\pmb X}_{\psi 1}+O(\psi^2) $$ (5) 式中,第1项${\pmb X}_{\psi 0}$为顺直项部分 ($\psi^{0})$,对应$\psi=0$时的直槽流,为实向量函数;第2项为一阶弯曲修正项,与槽道中心线摆动同频同波数,表示为复数形式,${\pmb X}_{\psi 1}=\hat {\pmb X}_{\psi 1}\exp ({\rm i} \phi_{\rm c})+{\rm c.c}$,其幅值:$\hat {\pmb X}_{\psi 1}=[\hat u_{s\psi 1}, \hat u_{n\psi 1}, \hat p_{\psi 1}]^{\rm T}$,顶标$\hat { \ \ }$表示波动量的幅值,为复向量函数,${\pmb X}_{\psi 1}=\hat {\pmb X}_{\psi 1}(n)$.
$\psi^{0}$阶对应固定直槽道流动,其对应基本流解为Poiseuille流流速分布
$$u_{s\psi 0}=1-n^2 $$ (6) $\psi^{1}$阶反映了槽道弯曲对水流流动造成的影响,其解对应的方程为
$$\left. \begin{array}{l} {\rm{i}}{\alpha _{\rm{c}}}{{\hat u}_{s\psi 1}} + {\rm{D}}{{\hat u}_{n\psi 1}} = 0\\ \left( {{{\rm{D}}^2} - \alpha _{\rm{c}}^2} \right){{\hat u}_{s\psi 1}} - {\rm{i}}Re{\alpha _c}{{\hat p}_{\psi 1}} + {\rm{i}}Re\left( {{\omega _{\rm{c}}} - {\alpha _{\rm{c}}}{u_{s\psi 0}}} \right){{\hat u}_{s\psi 1}} - \\ Re{\rm{D}}{u_{s\psi 0}}{{\hat u}_{n\psi 1}} = nRe{p_{\psi 0, s}} + {\rm{D}}{u_{s\psi 0}}\\ \left( {{{\rm{D}}^2} - \alpha _{\rm{c}}^2} \right){{\hat u}_{n\psi 1}} - ReD{{\hat p}_{\psi 1}} + {\rm{i}}Re{\omega _{\rm{c}}}{{\hat u}_{n\psi 1}} - \\ {\rm{i}}Re{\alpha _{\rm{c}}}{u_{s\psi 0}}{{\hat u}_{n\psi 1}} = Reu_{s\psi 0}^2 + {\rm{i}}{\alpha _{\rm{c}}}{u_{s\psi 0}} \end{array} \right\}$$ (7) 边界条件:$n=\pm 1, \ \hat u_{s\psi 1}=nu_{\rm c}, \ \hat u_{n\psi 1}=\dfrac{{\rm i}u_{\rm c}}{\alpha_{\rm c}}$.
方程 (7) 中,$\hat u_{s\psi 1}, \hat u_{n\psi 1}, \hat p_{ \psi 1}$为未知量,导数符号:${\rm D}={\rm d}/{\rm d}n$,${\rm D}^{2}={\rm d}^{2} / {\rm d}n^{2}$. s方向驱动力$n Re \hat p_{ \psi 0, s}$和${\rm D}\hat u_{s\psi 0}$分别由零阶压力梯度和曲面切应力提供,压力梯度$\hat p_{ \psi 0, s}=Re^{-1} {\rm D}^2 \hat u_{s\psi1}$;n方向驱动力$Re \hat u_{ s\psi 0 }^2$和${\rm i}\alpha_{\rm c} \hat u_{s\psi 0}$分别由离心力和曲面切应力提供, Re为雷诺数.
对$\psi^{1}$阶方程 (7) 简化,可得一个O-S方程形式的常微分方程
$$ {\rm L}\hat u_{n\psi 1}=-Re \alpha^2_{\rm c} u_{s\psi 0}^2-{\rm i} \alpha_{\rm c}({\rm D}^3+{\rm D}^2)u_{s\psi 0}-\\ \qquad {\rm i} \alpha_{\rm c}({\rm D}^2+\alpha^2_{\rm c}) u_{s\psi 0} $$ (8) 边界条件:$n=\pm 1$, $\hat u_{n\psi 1}=\dfrac{{\rm i}u_{\rm c}}{\alpha_{\rm c}}$, ${\rm D}\hat u_{n\psi1}=\dfrac{{\rm i}nu_{\rm c}}{\alpha_{\rm c}}$.
式中,算子${\rm L}=({\rm D}^2-\alpha^2_{\rm c})^2-{\rm i} Re [(\alpha_{\rm c} u_{s\psi 0}-\omega_{\rm c})({\rm D}^2-\alpha^2_{\rm c})-\alpha_{\rm c}{\rm D}^2 u_{s\psi 0}]$, 等式右侧3项分别由离心力、压力梯度和曲面应力引起.
结合边界条件,差分法求解方程,采用横向网格数$N=1 000$,由式 (8) 可解得横向流速幅值$\hat u_{n\psi1}$,将其代回方程 (7) 的连续方程和s向动量方程,可解得一阶纵向流速和压力的幅值$\hat u_{s\psi 1}$和$ \hat p_{\psi 1}$.
值得注意的是,为满足摄动方法对式 (8) 中边界条件$\hat u_{s\psi 1} $和$\hat u_{n\psi 1}$的要求,摆动相速度$u_{\rm c}$和摆动频率$\omega_{\rm c}$与摆动波数$\alpha_{\rm c}$要求一定的量级关系:$u_{\rm c}\sim O (\alpha_{\rm c})$; $\omega_{\rm c }\sim O (\alpha_{\rm c}^{2})$.这是因为,摆动频率$\omega_{\rm c}$较小时,Stokes厚度$d^*_{\rm s}=\sqrt{2v^*/\omega^*_{\rm c}}$足够大,其值远大于前述半宽$B^{\ast}$,水流仍以压力梯度引起的流动为主.而当摆动频率$\omega_{\rm c}$增大时,Stokes厚度$d_{\rm s}^{\ast}$逐渐接近特征长度$B^{\ast }$,摆动引起的振荡流开始与压力流等量齐观,甚至占据主导地位.因此,摄动方法不再适应.若要进行高摆动频率影响下水流流动特征的研究可具体参考Thomas等[7]以及Luo和Wu[22]的成果.
2.2 摆动波数$\alpha_{\rm c}$对一阶基本流解的影响
对于固定的弯曲槽道,摆动角频率$\omega_{\rm c}=0$.雷诺数$Re=2 000$时,求解控制方程 (7),可得一阶流速分量$\hat {u}_{s\psi 1}$和$\hat {u}_{n\psi 1}$沿n方向分布图,见图 2和图 3.
图 2为一阶流速分布在复数空间的分解,函数Re ( ), Im ( ) 分布表示复数$\hat {u}_{s\psi 1}$和$\hat {u}_{n\psi1} $的实部和虚部.实部${\rm Re}(\hat {u}_{i\psi 1})$为相位$\phi=0$的正弯顶处的流速分布,虚部${\rm Im}(\hat {u}_{i\psi 1})$则为相位$\phi=-\pi /2$的零曲率的直槽处流速分布,提前${\rm Re} (\hat {u}_{i\psi 1})\pi /2$个相位.
图 2显示,一阶s方向流速$\hat {u}_{s\psi 1}$呈反对称分布,而n方向流速$\hat {u}_{n\psi 1}$则呈对称分布;一阶s方向流速$\hat {u}_{s\psi 1}$在近两壁处有最大值和最小值,极值大小和所处位置随摆动波数$\alpha_{\rm c}$增变化而变化.
随摆动波数$\alpha_{c}$增大,一阶n方向流速$\hat {u}_{n\psi 1}$随摆动波数$\alpha _{\rm c}$增大明显增大;一阶s方向流速$\hat {u}_{s\psi 1}$极值所处位置随摆动波数$\alpha_{\rm c}$增大而逐渐靠近两壁,在$\alpha_{\rm c}=(0, 0.5)$的范围内实部${\rm Re}(\hat {u}_{s\psi 1})$逐渐增大,而虚部${\rm Im}(\hat {u}_{s\psi1})$则先增大后减小.
这是因为,无量纲曲率幅值$\psi $是摆动波数$\alpha_{\rm c}$和偏角幅值$\theta _{\rm m}$的乘积,当偏角幅值保持不变时,$\psi $会随着$\alpha_{\rm c}$的增大而增大.因此,$\alpha_{\rm c}$增大就意味着河槽的摆动幅度增加,此时,s方向的最大水流流速会由于惯性的作用偏向凹岸一侧,而且摆动幅度越大,最大流速的偏离就越严重,甚至开始慢慢贴近凹岸边界.
图 3则为一阶流速分布在幅值-相位空间的分解,Am ( ) 和Ph ( ) 分别表示复数$\hat u_{s\psi1}$和$\hat u_{n\psi 1}$的幅值和相位 (弧度制).由图 3和图 2可知,纵向流速$\hat u_{s\psi 1}$和横向流速$\hat u_{n\psi 1}$的幅值呈对称分布,$\hat u_{n\psi 1}$的相位呈对称分布,但$\hat u_{s\psi1}$的相位在槽道两侧有180°的相位差,亦即
$$\left\{ \begin{array}{l} {\rm{Am}}({{\hat u}_{s\psi 1}}) = {\rm{Am}}( - {{\hat u}_{s\psi 1}})\\ {\rm{Ph}}({{\hat u}_{s\psi 1}}) = {\rm{Ph}}( - {{\hat u}_{s\psi 1}}) + \pi \\ {\rm{Am}}({{\hat u}_{n\psi 1}}) = {\rm{Am}}( - {{\hat u}_{n\psi 1}})\\ {\rm{Ph}}({{\hat u}_{n\psi 1}}) = {\rm{Ph}}( - {{\hat u}_{n\psi 1}}) \end{array} \right.$$ 一阶纵向流速的幅值${\rm Am}(\hat u_{s\psi 1})$在近壁处有最大值,槽道中心流速相位比边壁处略小.随摆动波数$\alpha_{\rm c}$增大,纵向流速的幅值${\rm Am}(\hat {u}_{s\psi 1})$最大值所处位置更加靠近边壁.一阶横向流速的幅值${\rm Am}(\hat {u}_{n\psi 1})$在随摆动波数$\alpha_{\rm c}$增大而增大.
一阶纵向流速$\hat u_{s\psi 1}$和横向流速$\hat u_{n\psi 1}$随摆动波数$\alpha_{\rm c}$增大而减小.当摆动波数$\alpha _{\rm c}$较小时,相位的横向分布变化较小,其分布类似抛物线,如图 3中$\alpha_{\rm c}=0.01$的曲线所示,边壁处的相位差大约为0.8$\pi $.当摆动波数$\alpha_{\rm c}$较大时,相位的横向分布变化较大,其分布类似圆管湍流的分布情况,边界相位梯度很大,中心区域很平坦,相位变化很小,如图 3中$\alpha_{\rm c}=0.50$的曲线所示,这种趋势随摆动波数$\alpha_{\rm c}$增大越加明显.
2.3 摆动角频率$\omega_{\rm c}$对一阶基本流解的影响
对于摆动的弯曲槽道,摆动角频率$\omega_{\rm c} \ne 0$.摆动波数$\alpha _{\rm c}=0.3$, 雷诺数$Re=2 000$时,求解控制方程 (1),所得一阶流速分量$\hat u_{s\psi1}$和$\hat {u}_{n\psi 1}$的横向分布绘图,见图 4和图 5.
图 4为一阶流速分布在复数空间的分解.一阶纵向流速$\hat {u}_{s\psi 1}$呈反对称分布,横向流速$\hat {u}_{n\psi 1}$则呈对称分布.同时,由控制方程 (7) 的边界条件可知:$n=\pm 1$时,${\rm Re}(\hat {u}_{s\psi 1})$和${\rm Im}(\hat {u}_{n\psi 1})$均不为零,尤其${\rm Im}(\hat {u}_{n\psi 1})$在图 4(b)中极为明显.这是因为在所设的条件下,固体边界是以波的形式传播的,相位$\phi_{\rm c}=\pi / 2$处的边界会以相对较大的速度沿n方向运动,从而造成了此断面边界处 ($n=\pm 1$) 的速度不为零.而在槽道中心线 ($n=0$) 附近,零阶s方向的流速较大,若河槽摆动 ($\omega_{\rm c} \ne0$),则会大大影响到原有流速分布,具体表现为一阶流速${\rm Re}(\hat u_{n\psi 1})$和${\rm Im}(\hat {u}_{n\psi1})$对零阶流速的大幅度修正,即${\rm Re}(\hat {u}_{n\psi 1})$和${\rm Im}(\hat {u}_{n\psi 1})$在$n=0$处获得极值.
摆动槽道的一阶流速分布取决于两方面的影响,一方面源自于槽道弯曲, 即摆动波数$\alpha _{\rm c}$的影响, 一方面源自于摆动频率$\omega_{\rm c}$的影响. 图 4中黑色实线为$\omega _{\rm c}=0$的固定弯曲槽道中的一阶流速分布,取$\omega_{\rm c}=-0.10$, $-0.05$, 0, 0.05, 0.10五种不同的摆动情况,并以实线和虚线区分摆动传播方向,实线对应$\omega _{\rm c}>0$,虚线对应$\omega_{\rm c} < 0$.由图 4可见,摆动向下游传播时,摆动频率$\omega _{\rm c}$对一阶流速影响比向上游传播大,而且,整体来说,摆动频率$\omega_{\rm c}$对一阶流速的影响与摆动波数对它的影响是相反的 (图 2).这是因为河槽中的水流流速方向是向下游的,因此,水流的信息主要是向下游传播的,摆动向下游传播间接推动了水流信息的传播效果.然而,又由于水流和河槽边界同时向下游运动,它们之间的相对速度就减小了,就相当于间接减小了波数的影响.因此,摆动频率对一阶流速的影响与摆动波数对它的影响效果是相反的.
随着摆动频率$\omega_{\rm c}$由零减小,一阶纵向流速和横向流速分布迅速反向,且$\omega _{\rm c}=0.05$和$\omega_{\rm c}=0.10 $两种摆动频率下相差不大.
随着摆动频率$\omega_{\rm c}$由零增大,$\omega_{\rm c}=0.05$时,弯顶处两壁附近出现一阶纵向流速${\rm Re}(\hat u_{s\psi 1})$的反向分布 (与$\omega_{\rm c}=0$时的黑线分布对比),这时摆动波数$\alpha _{\rm c}$和摆动频率$\omega_{c}$的影响相当;摆动频率继续增大$\omega_{\rm c}=0.10$时,反向流速分布扩展到全断面,这时摆动频率$\omega_{\rm c}$的影响占主导地位.直槽处的一阶横向流速${\rm Im}(\hat {u}_{n\psi 1})$分布随摆动波数的变化与弯顶处${\rm Re}(\hat {u}_{s\psi 1})$分布相似.二者随$\omega_{\rm c}$的复杂变化规律是由于边界流速非零,因此其近壁流速才呈现局部反向分布的现象.
弯顶处一阶横向流速${\rm Re}(\hat {u}_{n\psi 1})$和直槽处一阶纵向流速${\rm Im}(\hat {u}_{s\psi 1})$则在$\omega_{\rm c}=0.05$时远大于$\omega_{\rm c}=0.10$的情况.
图 5为一阶流速分布在幅值-相位空间的分解,一阶纵向流速$\hat {u}_{s\psi1}$和横向流速$\hat {u}_{n\psi 1}$的幅值和相位均呈对称分布.
图 5(a)显示,一阶纵向流速幅值${\rm Am}(\hat {u}_{s\psi 1})$在近壁处有最大值,该值随摆动频率$\omega_{\rm c}$由零向正负两个方向增大时均增大,且摆动向下游传播时明显大于向上游传播.随着摆动频率$\omega_{\rm c}$由零减小,一阶纵向流速相位${\rm Ph} (\hat {u}_{s\psi 1})$迅速反向分布,且$\omega_{\rm c}=0.05$和$\omega_{\rm c}=0.10$两种摆动频率下相差不大;随着摆动频率$\omega_{\rm c}$由零增大,相位分布也迅速反向,中心处相位大于边壁,$\omega_{\rm c}=0.10$时的相位分布比$\omega_{\rm c}=0.05$大.
图 5(b)显示,一阶横向流速$\hat {u}_{n\psi 1}$的幅值${\rm Am} (\hat {u}_{n\psi 1})$对称分布.随着摆动频率$\omega_{\rm c}$由零减小,幅值${\rm Am} (\hat {u}_{n\psi 1})$增大,$\omega_{\rm c}=0.05$和$\omega_{\rm c}=0.10$两种摆动频率下相差较大;随着摆动频率$\omega_{\rm c}$由零减小,幅值${\rm Am}(\hat {u}_{n\psi 1})$在近壁处非零,且局部呈反向分布,$\omega_{\rm c}=0.05$和$\omega_{\rm c}=0.10$时相差不大.一阶横向流速相位${\rm Ph} (\hat {u}_{n\psi 1})$分布变化规律与前述${\rm Ph} (\hat {u}_{s\psi 1})$相同.
同时,从图 5(b)中可以观察到$\omega_{\rm c}=0.10$对应的蓝色实线从1.0突然跳转到$-1.0$,这是因为图 5(b)中第2个图代表的是$u_{n\psi 1}$的相位,而通过反正切计算相位的值域在$(-\pi, \pi)$之间.这种跳转体现了相位的值域范围的约束,其本身只有数学上的意义,在物理上${\rm Ph} (u_{n\psi 1})/\pi$仍然是一个连续函数.
3. 水动力稳定性及拟序扰动发展的边界效应
按照描述湍流拟序结构理论方法[23-25],方程的解${\pmb X}=[u_{s}, u_{n}, p]^{\rm T}$分解为基本流解${\pmb X}_{\psi }$与扰动解${\pmb X}_{T}$之和的形式,${\pmb X}={\pmb X}_{\psi } + {\pmb X}_{T}$.
与平面Poiseuille流相比,摆动槽道基本流中的波动项部分${\pmb X}_{\psi1}$的使内部水流既有流动的自然失稳又有由于摆动弯曲诱发的失稳,扰动量${\pmb X}_{T}=[\hat {u}_{sT}, \hat {u}_{nT}, \hat p_T]^{\rm T}$可表述为以下形式
$${\pmb X}_T=\varepsilon_T\sum^1_{m=-1} \hat {\pmb X}_{Tm}\exp [{\rm i}(\phi_{ T}+m\phi_{\rm c})]+o(\varepsilon^2_T) $$ (9) 式中,扰动相位$\phi_T=\alpha_Ts-\omega_Tt$, 由于基本量${\pmb X}_{\psi }$对$\psi $取一阶近似,扰动解中包含$-1$, 0, 1次谐波量,故式 (9) 在$m=\pm 1$处截断,${\pmb X}_{T}$的幅值$\hat {\pmb X}_{Tm}$为n的复数函数,$\hat {\pmb X}_{Tm}=\hat {\pmb X}_{Tm}(n)$.
将式 (9) 中的拟序扰动${\pmb X}_{T}$、基本解${\pmb X}_{\psi }$代入原始控制方程 (2) $\sim $方程 (4),得到拟序扰动的控制方程为
$$\dfrac 1{h_s} \dfrac{\partial u_{sT} }{\partial s}+\dfrac{\partial u_{nT} }{\partial n}-\psi \dfrac{ncu_{nT}}{h_s}=0 $$ (10) $$\dfrac{\partial u_{sT} }{\partial t} + \Big(\dfrac{u_{s\psi 0} }{h_s}\dfrac{\partial u_{sT} }{\partial s} + u_{nT} \dfrac{\partial u_{s\psi 0}}{\partial n} - \psi \dfrac{cu_{s\psi 0} u_{nT}}{h_s } \Big )+ \\ \qquad \dfrac{\psi }{h_s }\left( {u_{s\psi 1} \dfrac{\partial u_{sT} }{\partial s} + u_{sT}\dfrac{\partial u_{s\psi 1} }{\partial s}} \right) + \psi \left( {u_{n\psi 1} \dfrac{\partial u_{sT}}{\partial n} + }\right. \\ \qquad \left.{u_{nT} \dfrac{\partial u_{s\psi 1} }{\partial n}} \right)-\psi ^2c\dfrac{u_{s\psi 1} u_{nT} + u_{n\psi 1} u_{sT} }{h_s } = -\dfrac{1}{h_s }\dfrac{\partial p_T }{\partial s} + \\ \qquad \dfrac{1}{Re}\bar {\Delta }u_{sT} -\dfrac{\psi }{R}\left( {\dfrac{2c}{h_s^2}\dfrac{\partial u_{nT} }{\partial s} + \dfrac{c}{h_s }\dfrac{\partial u_{sT} }{\partial n} + \psi\dfrac{c^2u_{sT} }{h_s^2 }} \right) + \\ \qquad \dfrac{\psi }{Re}\dfrac{\partial c}{\partial s}\left( {\dfrac{n}{h_s^3 }\dfrac{\partial u_{sT} }{\partial s} - \dfrac{u_{nT} }{h_s^2 }} \right) $$ (11) $$ \dfrac{\partial u_{nT} }{\partial t} +\Bigg ( h^{-1}_s \dfrac{u_{s\psi 0} }{h_s}\dfrac{\partial u_{nT} }{\partial s} + u_{nT} \dfrac{\partial u_{s\psi 0} }{\partial n} + \psi c\dfrac{2u_{s\psi 0} u_{sT} }{h_s }\Bigg) + \\ \qquad \dfrac{\psi }{h_s }\left( {u_{s\psi 1} \dfrac{\partial u_{nT} }{\partial s} + u_{sT}\dfrac{\partial u_{n\psi 1} }{\partial s}} \right) + \psi u_{n\psi1} \dfrac{\partial u_{nT} }{\partial n}+ \\ \qquad \psi u_{nT} \dfrac{\partial u_{n\psi 1} }{\partial n} + \psi ^2c\dfrac{2u_{s\psi 1}u_{sT} }{h_s } = - \dfrac{\partial p_T }{\partial n} + \dfrac{1}{Re}\bar {\Delta }u_{nT}+ \\ \qquad \dfrac{\psi }{Re}\left( {\dfrac{2c}{h_s^2 }\dfrac{\partial u_{sT} }{\partial s} -\dfrac{c}{h_s }\dfrac{\partial u_{nT} }{\partial s} - \psi\dfrac{c^2u_{nT} }{h_s^2 }} \right)+ \\ \qquad \dfrac{\psi }{Re}\dfrac{\partial c}{\partial s}\left( {\dfrac{n}{h_s^3}\dfrac{\partial u_{nT} }{\partial s} - \dfrac{u_{sT} }{h_s^2 }} \right) $$ (12) 边界条件$n=\pm 1$, $u_{sT}=u_{nT}=0$.
式 (10) $\sim $式 (12) 为拟序结构扰动项控制方程.方程中既有边界波动 (拉梅系数$h_{s}$) 对扰动解的直接影响,也有边界波动引起的一阶基本流${\pmb X}_{\psi 1}$对扰动解的间接影响.扰动方程的非线性使得各扰动谐波量相互耦合,这种耦合有两方面,一种是由高阶基本流 ($u_{s\psi i}$, $u_{n\psi i}$) 的色散作用引起,一种是由拉梅系数$h_{s}$的Taylor高阶展开项的色散作用引起.由于本文只保留一阶基本流和三种扰动谐波量,前者只对相邻扰动谐波起作用,而拉梅系数$h_{s}$则在各扰动分量之间都能起到影响, 故需保留$h_{s}$的$\psi $的二阶展开项,即
$$h_s=1+\psi nc+\dfrac 12 \psi^2n^2c^2+{\rm c.c}+o(\psi^3) $$ (13) 式中,c.c表示复共轭.将式 (13) 代入拟序扰动方程 (10) $\sim $方程 (12),可得扰动幅值方程
$$\left[{ {\rm i}\left( {\alpha _T + m\alpha _{\rm c} } \right)\hat {u}_{sTm} + {\rm D}\hat {u}_{nTm} } \right] + \\ \qquad \dfrac{1}{2}\psi \varPhi _{1\psi 1s} \hat {u}_{sTm} + \dfrac{1}{2}\psi \varPhi _{1\psi1n} \hat {u}_{nTm} + \\ \qquad \dfrac{1}{4}\psi ^2\left( {\varPhi _{1\psi 2s} \hat {u}_{sTm} + \varPhi _{1\psi 2n}\hat {u}_{nTm} } \right) = 0 $$ (14) $$\left\{ {{\rm D}^2 - \left( {\alpha _T + m\alpha _{\rm c} } \right)^2 - {\rm i}Re \left[{\omega _T + m\omega _{\rm c}-}\right. } \right. \\ \qquad \left.{\left.{\left( {\alpha _T + m\alpha _{\rm c} } \right)u_{s\psi 0} } \right]}\right\}\hat {u}_{sTm} + \\ \qquad Re {\rm D}u_{s\psi 0} \hat {u}_{nTm} + {\rm i}Re\left( {\alpha _T + m\alpha _{\rm c}} \right)\hat {p}_{Tm} + \\ \qquad \dfrac{1}{2}\psi \left( {\varPhi _{2\psi 1s} \hat {u}_{sTm} + \varPhi _{2\psi 1sd}{\rm D}\hat {u}_{sTm} + \varPhi _{2\psi 1n} \hat {u}_{nTm} } \right) + \\ \qquad \dfrac{1}{4}\psi ^2\left( {\varPhi _{2\psi 2s} \hat {u}_{sTm} + \varPhi _{2\psi 2sd}{\rm D}\hat {u}_{sTm} + \varPhi _{2\psi 2n} \hat {u}_{nTm} } \right) = 0$$ (115) $$ \left\{ {{\rm D}^2 - \left( {\alpha _T + m\alpha _{\rm c} } \right)^2 - {\rm i}Re\left[{\omega _T + m\omega _{\rm c}-}\right.}\right. \\ \qquad \left.{\left.{\left( {\alpha _T + m\alpha _{\rm c} } \right)u_{s\psi 0} } \right]}\right\}\hat {u}_{nTm} + Re {\rm D}u_{s\psi 0} \hat {u}_{nTm} + Re {\rm D}\hat {p}_{Tm} + \\ \qquad \dfrac{1}{2}\psi \left( {\varPhi _{3\psi 1s} \hat {u}_{sTm} + \varPhi _{3\psi 1n} \hat {u}_{nTm} + \varPhi _{3\psi 1nd} {\rm D}\hat {u}_{nTm} } \right) + \\ \qquad \dfrac{1}{4}\psi ^2\left( {\varPhi _{3\psi 2s} \hat {u}_{sTm} + \varPhi _{3\psi 2n}\hat {u}_{nTm} + \varPhi _{3\psi 2nd} {\rm D}\hat {u}_{nTm} } \right) = 0 $$ (16) 式中,$m=0$, $\pm 1$,式 (14) $\sim$式 (16) 构成9元耦合方程组. $\alpha_{T }$和$\omega_{T}$分别为无量纲形式的扰动波数和扰动角频率;$\alpha_{{\rm c}}$和$\omega_{\rm c}$分别为无量纲边界摆动波数和摆动角频率.
式 (14) $\sim$式 (16) 代表弯曲河道水动力稳定性特征.方程组中未知量包含9个线性无关的特征变量:$\hat u_{sT0}, \hat u_{nT0}, \hat p_{T0}$, $\hat u_{sT+1}, \hat u_{nT+1}, \hat p_{T+1}$, $\hat u_{sT-1}, \hat u_{nT-1}, \hat p_{T-1}$.这些特征变量均为自变量n的特征函数,且相互耦合,以矩阵形式可表述为
$${\pmb A}\hat {\pmb X}_T+{\pmb B}\hat {\pmb X}_{T, n}+{\pmb C}\dot{\pmb X}_{T, nn}=\boldsymbol{0}$$ (17) 式中,$\hat {\pmb X}_T$为扰动幅值方程的基础解系,${\pmb A}, {\pmb B}, {\pmb C}$为复数域上的9$\times$9阶系数矩阵,系数矩阵各元素均由 (Re, $\alpha_{T }$, $\omega_{T})$,($\theta_{\rm m}$, $\alpha_{\rm c}$, $\omega_{\rm c})$,($u_{s\psi 0}$, $\hat {u}_{s\psi 1}$, $\hat {u}_{n\psi 1})$决定.
基础解系$\hat {\pmb X}_T$表述为向量形式,如下
$$\begin{array}{c} {{\boldsymbol{\hat X}}_T} = \left[{{{\hat u}_{sT-1}}\left( n \right)\;\;{{\hat u}_{sT0}}\left( n \right)\;{{\hat u}_{sT + 1}}\left( n \right)\;\;{{\hat u}_{nT-1}}\left( n \right)} \right.\\ {\left. {{{\hat u}_{nT0}}\left( n \right)\;\;{{\hat u}_{nT + 1}}\left( n \right)\;\;{{\hat p}_{T-1}}\left( n \right)\;\;{{\hat p}_{T0}}\left( n \right)\;\;{{\hat p}_{T + 1}}\left( n \right)} \right]^{\rm{T}}} \end{array}$$ 系数矩阵${\pmb A}, {\pmb B}, {\pmb C}$各元素的显示表述可由方程 (14)$\sim$方程 (16) 及附录得到,限于篇幅不再赘述.系数矩阵中共有两组参数,一组为与直槽一样的拟序扰动参数组 ($\alpha_{T}, \omega_{T}, Re)$,还有一组为边界弯曲参数组$ (\theta_{\rm m}$, $\alpha_{\rm c}$, $\omega_{\rm c})$.
特征值方程 (17) 采用时间模式,$\alpha_{T }$为实数,$\omega_{T}$为复数,用muller法求解扰动幅值方程的特征值.扰动角频率$\omega_{T}$为 ($\alpha_{T }$, $Re)$和 ($\theta _{\rm m}$, $\alpha_{\rm c}$, $\omega_{\rm c})$的函数,亦即:$\omega_{T}=\omega _{T}(\alpha _{T}$, Re, $\theta_{\rm m}$, $\alpha_{\rm c}$, $\omega_{\rm c})$,其虚部即为扰动增长率:$\omega_{Ti }=\omega _{Ti }(\alpha_{T}$, Re, $\theta_{\rm m}$, $\alpha_{\rm c}$, $\omega_{\rm c})$.
正弦派生曲线河槽内的扰动波,不仅与扰动波数$\alpha_{T}$和雷诺数R的有关,还受到边界摆动弯曲 ($\theta_{\rm m}$, $\alpha_{\rm c}$, $\omega_{\rm c})$的影响.直槽水流雷诺数Recr=5 772.222(Orszag[26]),而弯曲槽道流的临界雷诺数$Re_{\rm cr}$受到边界摆动弯曲参数 ($\theta_{\rm m}$, $\alpha_{\rm c}$, $\omega_{\rm c})$的影响,亦即:$ Re_{\rm cr}=Re_{\rm cr}$($\theta _{\rm m}$, $\alpha_{\rm c}$, $\omega_{\rm c})$,与此同时,临界雷诺数对应的临界扰动参数 ($\alpha _{T\rm cr}$, $\omega_{T\rm cr})$均为 ($\theta_{\rm m}$, $\alpha_{\rm c}$, $\omega_{\rm c})$的函数,亦即,$\alpha_{T\rm cr}=\alpha_{T\rm cr}$ ($\theta_{\rm m}$, $\alpha_{\rm c}$, $\omega_{\rm c})$,$\omega_{T\rm cr}=\omega_{T\rm cr}$ ($\theta _{\rm m}$, $\alpha_{\rm c}$, $\omega_{\rm c})$.
由于$\psi=0$时矩阵奇异,本文以$\theta_{\rm m}=\alpha_{\rm c}=1\times 10^{-6}$代替顺直河道 ($\theta_{\rm m}=\alpha_{\rm c }=0$) 来验证临界雷诺数及中性曲线,采用目前国际最为公认的理论和实验结果验证模型的准确性,并与Reynold和Potter[27]的数值结果和Nishioka和Ichikawa[28]的实验结果进行了对比,见图 6所示.可以看出理论计算结果与已发表的最为公认的中性曲线结果吻合较好,说明程序计算順直河道计算结果的正确性.
3.1 边界摆动波数的影响
图 7为扰动频率$\omega_{Tr}$及扰动增长率$\omega_{Ti}$随摆动波数$\alpha_{\rm c}$的变化关系,对应平面状态参数:$\theta _{\rm m}=0.1$, $\omega_{\rm c}=0$, 拟序扰动参数:$\alpha_{T }=1$, $ Re=6 000$;图 8临界雷诺数随摆动波数$\alpha_{\rm c}$的变化,对应平面状态参数:$\theta_{\rm m}=0.1$, $\omega_{\rm c}=0$.
图 7 扰动频率$\omega_{Tr}$及扰动增长率$\omega_{Ti}$随摆动波数$\alpha_{\rm c}$的变化(平面状态参数:$\theta_{\rm m}=0.1$, $\omega_{\rm c}=0$; 拟序扰动参数:$\alpha _{T }=1.02$, $ Re=5 772.222$)Figure 7. Variation of disturbance frequency $\omega_{Tr}$ and growth rate $\omega_{Ti}$ with swinging wavenumber $\alpha_{\rm c}$(for: $\theta_{\rm m}=0.1$, $\omega_{\rm c}=0$; $\alpha_{T }=1.02$, $ Re=5 772.222$)从图 7(a)中可以看出,随着摆动波数的增加,扰动频率有持续减少的趋势.造成这种变化趋势的原因可能是扰动频率与扰动波数、摆动频率和摆动波数的一种适应过程.即,摆动波数在逐渐接近于扰动波数过程中,扰动频率也逐渐接近于摆动频率 ($\omega_{\rm c}=0$).这种结果虽然与Bai等[18]以及Xu等[17]的结果由于采用不同坐标系而产生一定的差别,但总体趋势是一致的,体现出模型计算的合理性.
从图 7(b)中可以看出,随着摆动波数的增加,扰动增长率出现先增大后减小的趋势.槽道弯曲对其内部水流各种尺度的扰动都有选择性影响.摆动波数$\alpha_{\rm cm}$在 (0, 0.48) 区间上,槽道摆动加速扰动增长,水流比直槽道水流更不稳定.当摆动波数$\alpha _{\rm cm}=0.3$,扰动增长最快,水流最不稳定,此时临界雷诺数不足直槽流的一半.
直槽水流扰动临界参数$Re_{\rm cr}=5 772.222$,$\alpha_{T\rm cr}=1.020$,$\omega_{T\rm cr}=0.27$,摆动槽道中心曲线和临界扰动显然受摆动参数 ($\theta_{\rm m}$, $\alpha_{\rm c}$, $\omega_{\rm c}$) 的影响.
图 9绘制了随摆动波数$\alpha_{\rm c}$增大水流中心曲线的移动,图中给出了$\alpha _{\rm c}$=0,0.1,0.2,0.3,0.4,0.5时的中性曲线.由图 9可以看出,在摆动波数$\alpha_{\rm c}$较小时,中性曲线变化缓慢,摆动波数$\alpha_{\rm c}$较大时,中性曲线变化迅速,且趋于平坦化.随摆动波数$\alpha_{\rm c}$增大,中性曲线向左移动,对应临界扰动也随之左移,临界雷诺数减小,至$\alpha_{\rm c}=0.3$时达到最小 (亦可见图 8),然后中心曲线右移,临界雷诺数增大.
图 10只显示了5条中心曲线的详细信息,计算$\alpha_{\rm c}$在 (0, 0.5) 区间上共50条中性曲线,将各中性曲线临界扰动对应的点列 ($Re_{\rm cr}$, $\alpha _{T\rm cr})$和 ($Re_{\rm cr}$, $\omega_{T\rm cr})$分别连接形成临界扰动曲线,绘制在图 8.由图 8可见,随着摆动波数$\alpha _{\rm c}$增大,临界扰动对应的临界扰动波数$\alpha_{T\rm cr}$先增大后减小然后再次增大,并在临近$\alpha_{\rm c}=0.3$附近形成一个"绳套";而临界扰动频率$\omega_{T\rm cr}$.则先增大后减小,形成类似于直槽中性曲线的"拇指"形状.
3.2 边界摆动频率的影响
图 11为扰动频率$\omega_{Tr}$及扰动增长率$\omega_{Ti}$随摆动频率$\omega_{\rm c}$的变化关系,对应平面状态参数:$\theta_{\rm m}=0.1$, $\alpha_{c}=0.3$, 拟序扰动参数:$\alpha_{T }=1.016$, $ Re=2 307$,该参数组合为图 8中最小临界雷诺数对应的扰动情况.由图可见,槽道摆动频率对其内部水流各种尺度的扰动有选择性影响;拟序扰动对摆动频率$\omega_{\rm c}$非常敏感,这是因为扰动频率的微小变化都能极大改变一阶流速分布,如图 4和图 5所示,流速分布的改变进一步影响拟序扰动的发展.
图 11 扰动频率$\omega_{Tr}$和扰动增长率$\omega_{Ti}$随摆动频率$\omega_{\rm c}$的变化(平面状态参数:$\theta_{\rm m}=0.1$, $-0.1 \leq \omega_{\rm c} \leq 0.1$; 拟序扰动参数:$\alpha_{T}=1.016$, $ Re=2 307$)Figure 11. Variation of disturbance frequency $\omega_{Tr}$ and growth rate $\omega _{Ti}$ with swinging frequency $\omega_{\rm c}$(for: $\theta_{\rm m}=0.1$, $-0.1 \leq \omega_{\rm c} \leq 0.1$; $\alpha _{T }=1.016$, $ Re=2 307$)由图 11扰动特征值随摆动频率$\omega_{\rm c}$的变化曲线可以看出,不同模态之间的交叉使得曲线偶有不可导的尖点和急剧变化的区段.这是因为,任意一种流动中都包含有各种尺度不同的扰动波,这种扰动波体现在特征方程的计算中则是可以计算出一系列的特征值.而湍流流动通常是由其中最不稳定/最少稳定的特征值模态所控制,当流动条件发生改变时有时会出现第二甚至第三模态增长率急剧增加,超过了第一模态的增长率,从而发生模态交叉现象.这种现象在波纹状底边界[16]水流结构的讨论中首次被发现,而在波状侧边界[18]的流动中被详细讨论,是目前弯曲型河槽水流结构与順直型河槽水流结构的重要区分之一.
扰动增长率$\omega_{Ti}$在$\omega_{\rm c}=0.085$和$\omega_{\rm c}=-0.039$处有极大值,在$\omega_{\rm c}=-0.006$附近有极小值,该极小值处模态交叉;$\omega_{\rm c}=-0.009$和$\omega_{\rm c}=0.064$附近,扰动频率$\omega_{Tr}$分别有极大值和极小值.
层流转捩为湍流,能量消耗会由于分子运动的加剧而急剧增加.然而根据最小耗能原理,自然的壁面边界和水流结构的相互适应过程中会自然遵循使水流的能耗降低,即从自然现象上则表现为推迟层流向湍流的转捩过程.Crosato[29]在论文中与VanBalen讨论了蜿蜒型河道的摆动规律,发现河道的发展不仅会向下游移动,在一些特定的条件下会向上游移动,其方向主要取决于冲刷池相对于河道顶点的位置[30].在相当陡峻的 (窄深型) 河岸,最大近岸流速将出现在河道顶点的上游,上游的河岸在水流流速的作用下变弱,从而造成河道向上游方向迁移.在计算中就体现两个重要的趋势:(1) 模态出现交叉,控制湍流结构的不稳定模态发生交换,从而最大的近岸流速位置出现在河道顶点的上游而非通常的下游;(2) 窄深型河道在水流结构的影响下,开始向上游迁移.表现在方程特征值上则是向上游迁移的摆动频率在一定范围内减小了扰动增长率 (如图 11(b)),即河道向上游迁移时,能量耗散更小,水流结构更加稳定,如图 12中所显示的趋势.
图 12为临界雷诺数$Re_{\rm cr}$随摆动频率$\omega_{\rm c}$的变化关系,对应平面状态参数:$\theta_{\rm m}=0.1$, $\alpha_{\rm c}=0.3$.随着摆动频率增大,临界雷诺数减小,内部水流更加容易失稳;与槽道摆动向下游传播 ($\omega_{\rm c}>0$) 相比,槽道摆动向上游传播 ($\omega_{\rm c} < 0$) 时水流更加稳定.
图 13为中性曲线随摆动频率$\omega_{\rm c}$的变化关系,对应平面状态参数:$\theta _{\rm m}=0.1$, $\alpha_{\rm c}=0.3$, $-0.02 \leq \omega_{\rm c} \leq 0.02$).图中绘出了$\omega_{\rm c}=-0.015$,0,0.015时的中性曲线.由图 13可以看出,随摆动频率$\omega_{\rm c}$增大,中性曲线左移,临界点 ($Re_{\rm cr}$, $\alpha_{T\rm cr})$向左上方移动,即临界雷诺数$Re_{\rm cr}$减小,临界扰动波数$\alpha _{T\rm cr}$增大,扰动频率$\omega_{T\rm cr}$减小,该变化趋势与图 11和图 12对应一致.槽道摆动向上游传播时,中性曲线随摆动频率增大的速度大于向下游传播时减小的速度.
4. 结论
本文建立了随体坐标系下槽道内水流流动的理论模型,讨论了正弦派生型摆动边界下的槽道水流动力稳定性特征,确定了摆动波数与摆动频率对主流流速和内部扰动发展影响的定量关系,得到了槽道弯曲度和摆动特征对其内部水流不同尺度扰动的选择性影响.具体如下:
(1) 一阶纵向流速$\hat {u}_{s\psi 1}$呈反对称分布,而横向流速$\hat {u}_{n\psi1}$则呈对称分布.
(2) 一阶纵向流速$\hat {u}_{s\psi 1}$在近两壁处有最大值和最小值,极值大小和所处位置随摆动波数$\alpha_{\rm c}$增大而逐渐靠近两壁.槽道中心处流速相位小于近壁处的,且随着摆动波数$\alpha_{\rm c}$增大而减小.
一阶纵向流速$\hat {u}_{s\psi 1}$和横向流速$\hat {u}_{n\psi 1}$随摆动波数$\alpha _{\rm c}$增大而减小.当摆动波数$\alpha_{\rm c}$较小时,相位的横向分布变化较小,其分布类似抛物线,边壁处的相位差大约为0.8$\pi$;当摆动波数$\alpha_{\rm c}$较大时,相位的横向分布变化较大,其分布类似圆管湍流的分布情况,边界相位梯度很大,中心区域趋于平坦,相位变化很小,且这种趋势随摆动波数$\alpha _{\rm c}$增大越加明显.
(3) 一阶流速分布对摆动频率$\omega_{\rm c}$的影响相当敏感,由$\omega_{\rm c}=0$向正负两个方向增大时流速分布迅速变化,且开始出现与$\omega_{\rm c}=0$时相反的分布规律;当$\omega_{\rm c} < -0.05$或$\omega_{\rm c} >0.1$时反向分布占据主导地位.
随着摆动频率$\omega_{\rm c}$由零减小,一阶流速、相位分布均迅速反向,且在$\omega _{\rm c} < -0.05$以后不再有明显变化;随着摆动频率$\omega_{\rm c}$由零增大,一阶流速首先在两壁附近开始反向,$\omega_{\rm c}>0.1$后反向分布占据主导地位,而流速相位则迅速反向,槽道中心处相位随摆动频率$\omega _{\rm c}$变化及其明显.
(4) 槽道弯曲对其内部水流各种尺度的扰动有选择性影响.摆动波数$\alpha_{\rm cm}$在 (0, 0.48) 区间上,槽道摆动加速扰动增长,水流比直槽道水流更不稳定.当摆动波数$\alpha _{\rm cm}=0.3$,扰动增长最快,水流最不稳定,此时临界雷诺数不足直槽流的一半.大摆动波数下中性曲线更加平坦.各中性曲线临界扰动对应的点列 ($Re_{\rm cr}$, $\alpha_{T\rm cr})$和 ($Re_{\rm cr}$, $\omega_{T\rm cr})$分别呈现"绳套"和"拇指"形状.
(5) 槽道摆动对水流内部各种尺度的扰动有较大影响.拟序扰动发展对摆动频率$\omega _{\rm c}$变化十分敏感.随着摆动频率$\omega_{\rm c}$增大,拟序扰动的增长率增大,临界雷诺数减小,水流更加容易失稳.槽道摆动向上游传播时水流更加稳定.
附录
式 (14) $\sim$式 (16) 中各参数为:
$$\varPhi_{1\psi 1s}=n({\rm e}^{{\rm i}\phi_{\rm c}}+{\rm c.c} )\cdot{\rm i}(\alpha_T+m \alpha_{\rm c}) \\ \varPhi_{1\psi 1n}=-({\rm e}^{{\rm i}\phi _{\rm c}}+{\rm c.c}) \\ \varPhi_{1\psi 2s}=n^2({\rm e}^{2{\rm i}\phi _{\rm c}}+{\rm c.c}+2)\cdot {\rm i}(\alpha_T+m\alpha_{\rm c}) \\ \varPhi_{1\psi 2n}=-n ({\rm e}^{2{\rm i}\phi _{\rm c}}+{\rm c.c}+2) $$ $$\varPhi _{2\psi 1s} = Ren\left( {{\rm e}^{{\rm i}\phi _{\rm c}} + c.c} \right) \cdot {\rm i}\left({\alpha _T + m\alpha _c } \right)u_{s\psi 0} + \\ \qquad Re\left( {\hat {u}_{s\psi 1} {\rm e}^{{\rm i}\phi _{\rm c}} + c.c} \right) \cdot {\rm i}\left( {\alpha _T + m\alpha_c } \right) + \left( {{\rm e}^{{\rm i}\phi _{\rm c}} + c.c} \right) + \\ \qquad {\rm i}\alpha _c Re\left({\hat {u}_{s\psi 1} {\rm e}^{{\rm i}\phi _{\rm c}}- c.c} \right) + 2\left( {\alpha _T +m\alpha _c } \right)^2n\left( {{\rm e}^{{\rm i}\phi _{\rm c}} + c.c} \right) + \\ \qquad {\rm i}\alpha_c n\left( {{\rm e}^{{\rm i}\phi _{\rm c}} - c.c} \right) \cdot {\rm i}\left( {\alpha _T + m\alpha _c } \right) \\ \varPhi _{2\psi 1sd} = Re\left( {\hat {u}_{n\psi 1}{\rm e}^{{\rm i}\phi _{\rm c}} + c.c} \right) - \left( {{\rm e}^{{\rm i}\phi _{\rm c}} + c.c} \right) \\ \varPhi _{2\psi 1n} = - Re\left( {{\rm e}^{{\rm i}\phi _{\rm c}} + c.c} \right) \cdot u_{s\psi0} + Re\left( {{\rm D}\hat {u}_{s\psi 1} {\rm e}^{{\rm i}\phi _{\rm c}} + c.c} \right) - \\ \qquad 2\left( {{\rm e}^{{\rm i}\phi _{\rm c}}+ c.c} \right) \cdot {\rm i}\left( {\alpha _T +m\alpha _c } \right) - {\rm i}\alpha _c \left( {{\rm e}^{{\rm i}\varphi _c } - c.c}\right)$$ $$\varPhi _{2\psi 2s} = Ren^2\left( {{\rm e}^{2{\rm i}\phi _{\rm c}} + c.c + 2} \right) \cdot{\rm i}\left( {\alpha _T + m\alpha _c } \right)u_{s\psi 0} + \\ \qquad Re n\left({{\rm e}^{{\rm i}\phi _{\rm c}} + c.c} \right)\left( {\hat {u}_{s\psi 1} {\rm e}^{{\rm i}\phi _{\rm c}}+ c.c} \right) \cdot {\rm i}\left( {\alpha _T + m\alpha _c } \right) + \\ \qquad Ren\left({{\rm e}^{{\rm i}\phi _{\rm c}} + c.c} \right) \cdot i\alpha _c \left( {\hat {u}_{s\psi 1}{\rm e}^{{\rm i}\phi _{\rm c}} - c.c} \right) - \\ \qquad Re\left( {{\rm e}^{{\rm i}\phi _{\rm c}} + c.c}\right)\left( {\hat {u}_{n\psi 1} {\rm e}^{{\rm i}\phi _{\rm c}} + c.c} \right) - \\ \qquad 3\left( {\alpha _T + m\alpha _c } \right)^2n^2\left( {{\rm e}^{2{\rm i}\phi _{\rm c}}+c.c + 2} \right) - \left( {{\rm e}^{2{\rm i}\phi _{\rm c}} + c.c + 2} \right) + \\ \qquad 3{\rm i}\alpha_c n^2\left( {{\rm e}^{2{\rm i}\phi _{\rm c}} - c.c} \right) \cdot {\rm i}\left( {\alpha _T + m\alpha _c } \right) \\ \varPhi _{2\psi 2sd} = - n\left( {{\rm e}^{2{\rm i}\phi _{\rm c}} + c.c+ 2} \right)$$ $$\varPhi _{2\psi 2n} = - Ren\left( {{\rm e}^{2{\rm i}\phi _{\rm c}} + c.c + 2}\right) \cdot u_{s\psi 0} - \\ \qquad Re\left( {\phi _c + c.c} \right)\left( {\hat {u}_{s\psi 1} {\rm e}^{{\rm i}\phi _{\rm c}} + c.c} \right) - 4n\left( {{\rm e}^{2{\rm i}\phi _{\rm c}} + c.c + 2} \right) \cdot \\ \qquad {\rm i}\left( {\alpha _T + m\alpha _c} \right) - 2{\rm i}\alpha _c n\left( {{\rm e}^{{\rm i}\phi _{\rm c}} - c.c} \right)$$ $$\varPhi _{3\psi 1s} = Re\left( {{\rm e}^{{\rm i}\phi _{\rm c}} + c.c} \right) \cdot 2u_{s\psi0} + Re {\rm i}\alpha _c \left( {\hat {u}_{n\psi 1} {\rm e}^{{\rm i}\phi _{\rm c}} - c.c} \right) \\ \varPhi _{3\psi 1n} = Re n\left( {{\rm e}^{{\rm i}\phi _{\rm c}} + c.c} \right) \cdot {\rm i}\left( {\alpha _T + m\alpha _c } \right)u_{s\psi 0} + \\ \qquad Re\left( {\hat {u}_{s\psi 1} {\rm e}^{{\rm i}\phi _{\rm c}} + c.c} \right) \cdot {\rm i}\left( {\alpha _T + m\alpha _c } \right) + \\ \qquad Re\left( {D\hat {u}_{n\psi 1} {\rm e}^{{\rm i}\phi _{\rm c}} +c.c} \right) + 2\left( {\alpha _T + m\alpha _c } \right)^2\cdot \\ \qquad n\left( {{\rm e}^{{\rm i}\phi _{\rm c}} + c.c} \right)-{\rm i}\alpha _c n\left( {{\rm e}^{{\rm i}\phi _{\rm c}}- c.c} \right) \cdot {\rm i}\left( {\alpha _T + m\alpha _c } \right) \\ \varPhi _{3\psi 1nd} = Re \left( {\hat {u}_{n\psi 1} {\rm e}^{{\rm i}\phi _{\rm c}} + c.c}\right) + {\rm e}^{{\rm i}\phi _{\rm c}} + c.c $$ -
图 7 扰动频率$\omega_{Tr}$及扰动增长率$\omega_{Ti}$随摆动波数$\alpha_{\rm c}$的变化
(平面状态参数:$\theta_{\rm m}=0.1$, $\omega_{\rm c}=0$; 拟序扰动参数:$\alpha _{T }=1.02$, $ Re=5 772.222$)
Figure 7. Variation of disturbance frequency $\omega_{Tr}$ and growth rate $\omega_{Ti}$ with swinging wavenumber $\alpha_{\rm c}$
(for: $\theta_{\rm m}=0.1$, $\omega_{\rm c}=0$; $\alpha_{T }=1.02$, $ Re=5 772.222$)
图 11 扰动频率$\omega_{Tr}$和扰动增长率$\omega_{Ti}$随摆动频率$\omega_{\rm c}$的变化
(平面状态参数:$\theta_{\rm m}=0.1$, $-0.1 \leq \omega_{\rm c} \leq 0.1$; 拟序扰动参数:$\alpha_{T}=1.016$, $ Re=2 307$)
Figure 11. Variation of disturbance frequency $\omega_{Tr}$ and growth rate $\omega _{Ti}$ with swinging frequency $\omega_{\rm c}$
(for: $\theta_{\rm m}=0.1$, $-0.1 \leq \omega_{\rm c} \leq 0.1$; $\alpha _{T }=1.016$, $ Re=2 307$)
-
[1] Davies C, Carpenter PW. Instabilities in a plane channel flow be-tween compliant walls. Journal of Fluid Mechanics, 1997, 352:205-243 doi: 10.1017/S0022112097007313
[2] Davies C, Carpenter PW. Numerical simulation of the evolution of Tollmien-Schlichting waves over finite compliant panels. Journal of Fluid Mechanics, 1997, 335:361-392 doi: 10.1017/S0022112096004636
[3] Hell M, Waters SL. Transverse flows in rapidly oscillating elastic cylindrical shells. Journal of Fluid Mechanics, 2006, 547:185-214 doi: 10.1017/S0022112005007214
[4] Guaus A, Airiau C, Bottaro A, et al. Effects of wall compliance on the linear stability of Taylor-Couette flow. Journal of Fluid Mechanics, 2009, 630:331-365 doi: 10.1017/S0022112009006648
[5] Pitman MW, Lucey AD. Stability of plane-Poiseuille flow interacting with a finite compliant panel. 17th Australasian Fluid Mechanics Conference, Auckland, 2010-12-5-9. New Zealand, 2010
[6] Hall P. The stability of the Poiseuille flow modulated at high frequencies//Proceedings of the Royal Society of London, London, 1975. England, 1975. 453-464
[7] Thomas C, Bassom AP, Blennerhassett PJ, et al. The linear stability of oscillatory Poiseuille flow in channels and pipes//Proceedings of the Royal Society A, 2011. Australia, 2011. 467:2643-2662
[8] Kuhlmann C, Wanschura M, Rath HJ. Flow in two-sided lid-driven cavities:non-uniqueness, instabilities, and cellular structures. Journal of Fluid Mechanics, 1997, 336:267-299 doi: 10.1017/S0022112096004727
[9] Ren L, Chen JG, Zhu KQ. Dual role of wall slip on linear stability of plane Poiseuille flow. Chinese Physics Letters, 2008, 25(2):601-603 doi: 10.1088/0256-307X/25/2/067
[10] Sen PK, Carpenter PW, Hegde S, et al. A wave driver theory for vortical waves propagating across junctions with application to those between rigid and compliant walls. Journal of Fluid Mechanics, 2009, 625:1-46 doi: 10.1017/S0022112008005545
[11] Friedkin JF. A laboratory study of the meandering of alluvial rivers. U. S. Army Engineer Waterways Experiment Station, Vicksburg, Mississippi, 1945. USA
[12] Brice JC. Platform properties of meandering processes, In:River Meandering//Elliott CM. ed. Proceedings of the conference Rivers, Louisiana, 1983-10-24-26. New Orleans:ASCE, 1984. 1-15
[13] Hooke JM, Redmond CE. Use of cartographic sources for analyzing river channel change with examples from Britain//Petts GE, ed. Historical of Large Alluvial Rivers:Western Europe, Wiley, 1989
[14] Parker G, Sawai K, Ikeda S. Bend theory of river meanders. Part 2. Nonlinear deformation of finite-amplitude bends. Journal of Fluid Mechanics, 1982, 115:303-314 doi: 10.1017/S0022112082000767
[15] Ikeda S, Parker G, Sawai K. Bend theory of river meanders. Part 1. Linear development. Journal of Fluid Mechanics, 1981, 112:363-377 doi: 10.1017/S0022112081000451
[16] Bai YC, Xu HJ. A study on the stability of laminar open-channel flow over a sandy rippled bed. Science in China, Series E:Engineering & Material Science, 2005, 35:53-73 http://mall.cnki.net/magazine/Article/JEXG200501007.htm
[17] Xu HJ, Bai YH. Theoretical analyses on hydrodynamic instability in narrow deep river with variable curvature. Applied Mathematics and Mechanics, 2015, 36(9):1147-1168 doi: 10.1007/s10483-015-1971-6
[18] Bai YC, Ji ZQ, Xu HJ. Instability and self-adaption character of turbulence coherent structure in narrow-deep river bend. Science China:Technology Science, 2012, 55:2990-2999 doi: 10.1007/s11431-012-4977-4
[19] Hooke JM. Magnitude and distribution of rates of river bank erosion. Earth Surface Processes and Landforms, 1980, 5:143-157 doi: 10.1002/(ISSN)1931-8065
[20] Hickin, E J. Lateral migration rates of rivers bends//Cheremisinoff PN, Chereminisoff NP, Cheng SL eds. Handbook of Civil Engineering. Lancaster, Pennsylvinia:Technomic Publishing, 1988
[21] Pyle CJ, Richards KS, Chandler JH. Digital photogrammetric monitoring of river bank erosion. Photogrammetric Record, 1997, 15(89):753-764 doi: 10.1111/phor.1997.15.issue-89
[22] Luo JS, Wu XS. On the linear instability of a finite stokes layer:Instantaneous versus Floquet modes. Physics of Fluids, 2010, 22:054106-1-054106-13 doi: 10.1063/1.3422004
[23] Herbert T. Secondary of boundary layers. Annual Review of Fluid Mechanics, 1988, 20:487-526 doi: 10.1146/annurev.fl.20.010188.002415
[24] Zhang ZS, Lilley GM. A theoretical model of coherent structure in a plate turbulent boundary layer//Turbulent Shear Flow Ⅲ. Berlin:Springer-Verlag, 1981
[25] 张兆顺.湍流.北京:国防工业出版社, 2002 Zhang Zhaoshun. Turbulence. Beijing:National Defence Industry Press, 2002(in Chinese)
[26] Orszag SA. Accurate solution of the Orr-Sommerfeld stability equation. Journal Fluid Mechanics, 1971, 50(4):689-703 doi: 10.1017/S0022112071002842
[27] Reynolds WC, Potter MC. Finite-amplitude instability of parallel shear flows. Journal of Fluid Mechanics, 1967, 27:465-492 doi: 10.1017/S0022112067000485
[28] Nishioka M, Ichikawa Y. An experimental investigation of the stability of plane Poiseuille flow. Journal of Fluid Mechanics, 1975, 72:731-751 doi: 10.1017/S0022112075003254
[29] Crosato A. Analysis and Modelling of River Meandering.[PhD Thesis]. Italia:University of Padua, 2008
[30] Seminara G, Zolezzi G, Tubino M, et al. Downstream and upstream influence in river meandering. Part 2. Planimetric development. Journal of Fluid Mechanics, 2001, 438:213-230 https://www.researchgate.net/publication/231859426_Downstream_and_upstream_influence_in_river_meandering_Part_2_Planimetric_development
-
期刊类型引用(5)
1. 冀自青 ,白玉川 ,徐海珏 . 蜿蜒边界下平面流的线性流动稳定性. 力学学报. 2023(05): 1075-1086 . 本站查看
2. 李彬,徐海珏,白玉川,冀自青. 常曲率弯道二维流动稳定性与非线性演化规律分析. 力学学报. 2021(01): 168-183 . 本站查看
3. 徐海珏,吴金森,白玉川. 波浪扰动下河口幂律异重流的动力场分布特性. 力学学报. 2019(06): 1699-1711 . 本站查看
4. 胡璐,闫寒,张文明,彭志科,孟光. 黏性流体环境下Ⅴ型悬臂梁结构流固耦合振动特性研究. 力学学报. 2018(03): 643-653 . 本站查看
5. 覃金帛,曾志强,梁藉,杨明祥,张健. GPU并行优化技术在水利计算中的应用综述. 计算机工程与应用. 2018(03): 23-29+63 . 百度学术
其他类型引用(2)