Couette-Taylor流的力学机理与能量转换
Dynamical Mechanism and Energy Conversion of Couette-Taylor Flow
收稿日期: 2018-03-22
基金资助: |
|
Received: 2018-03-22
Fund supported: |
|
作者简介 About authors
王贺元,E-mail:
同轴圆筒间Couette-Taylor流问题是典型的旋转流动问题,它是层流到湍流过渡的范例,国内外众多学者对其进行了深入的研究.该文探讨Couette-Taylor流问题的力学机理与能量转换,通过将Couette-Taylor流三模混沌系统转换成Kolmogorov形系统,把系统的力矩分为四种类型:惯性力矩,内力矩,耗散力矩和外力矩.通过不同力矩的结合分析和研究了Couette-Taylor流产生混沌的关键因素和物理意义.研究了哈密顿能量,动能和势能之间的相互转换.讨论了能量与雷诺数之间的关系.研究表明四种力矩的耦合是产生混沌的必要条件,而且只有耗散力矩和驱动力矩(外力矩)相匹配时,系统才能产生混沌,其中任何三种力矩耦合均不可能产生混沌.圆筒旋转产生的外力矩供给系统能量,能量增长导致流动失稳,从而产生泰勒漩涡和混沌,进而得出了Couette-Taylor流的能量转换和物理意义.引进Casimir函数分析系统的动力学行为和能量转换,并估计混沌吸引子的界.Casimir函数反映了能量转换和轨道与平衡点间的距离,数值结果仿真出它们之间的关系.
关键词:
There have been a lot of investigations which concern with rotating flow between two concentric cylinders (abbreviate frequently as Couette-Taylor Flow), Couette-Taylor flow is the typical rotation flow problems, It provides a paradigm from laminar to turbulent transition. Dynamical mechanism and energy conversion of Couette-Taylor flow are investigated in this paper, the Couette-Taylor flow chaotic system is transformed into Kolmogorov type system, which is decomposed into four types of torques:inertial torque, internal torque, dissipation and external torque. By the combinations of different torques, key factors of chaos generation and the physical meaning of Couette-Taylor Flow are studied. The conversion among Hamiltonian energy, kinetic energy and potential energy is investigated. The relationship between the energies and the Reynolds number is discussed. It concludes that the combination of the four torques is necessary conditions to produce chaos, and only when the dissipative torques matches the driving (external) torques can the system produce chaos, any combination of three types of torques cannot produce chaos. The external torque, which is provided by the rotation of the cylinder, supply the energy of the system, and that leads to Taylor vortex and chaos, the physical meaning and energy conversion of Couette-Taylor flow system are investigated. The Casimir function is introduced to analyze the system dynamics, and its derivation is chosen to formulate energy conversion. The bound of chaotic attractor is obtained by the Casimir function and Lagrange multiplier. The Casimir function reflects the energy conversion and the distance between the orbit and the equilibria. These relationships are illustrated by numerical simulations.
Keywords:
本文引用格式
王贺元.
Wang Heyuan.
1 引言
同轴圆筒间旋转流动的Couette-Taylor流问题是近一个世纪以来人们普遍关注的热点问题,由于其流动形态的可观测性以及它在湍流研究中的基础性地位及其在流体机械、石油化工等领域的广泛应用,国际上将其列为非线性科学的范例之一.国内外众多学者对其复杂的动力学行为进行了大量深入的研究,相关文献非常丰富(文献[1-12]只是2000多篇文献中较少的一部分).如果内圆筒旋转,外圆筒静止,这种流动历经Couette流、Taylor涡流、Taylor行进波、波状螺旋漩涡、调制波状螺旋漩涡、湍状Taylor涡流等十几种波形最终进入湍流(见图 1),如果内外圆筒同时正向或反向旋转,可以观察到更为丰富的流动,具体见文献[3-12].
图 1
以往的研究工作大都侧重于从流动的稳定性和分岔理论开展研究,主要是利用分岔理论来解释和分析实验中观察到的流动发展到湍流前的各种涡流及其相互演化的过程,以及从层流过渡到湍流的方式及仿真等,而对流动的生成机理,能量转换及物理意义等问题很少有文献涉及.因此探讨Couette-Taylor流的生成机理和物理意义等相关问题是非常有意义和具有挑战性的.
混沌系统按其各项及相互作用是相当复杂的.如果从力学角度对混沌进行研究,可以从中发现更多的混沌基本成因.从力学角度研究混沌系统已经有了良好的开端,阿诺德[13]用Kolmogorov系统来描述具有哈密顿函数的不同的强迫动力系统,流体动力系统等. Pasini和Pelino[14]对Lorenz系统进行了研究,并给出了统一的柯尔莫哥洛夫(Kolmogorov)和洛伦茨(Lorenz)系统.齐和梁[15]引进Casimir函数分析了Chen混沌系统的动力学行为和能量转换,揭示了能量转换和轨道与平衡点间的距离,并估计了混沌吸引子的界.文献[16]研究了Chen混沌系统力学机理与能量转换,将Chen混沌系统转换成Kolmogorov形系统,通过不同力矩的结合分析和研究了Chen混沌系统产生混沌的关键因素和物理意义.研究了哈密顿能量,动能和势能之间的相互转换.讨论了能量与雷诺数之间的关系.文献[17]研究了齐四翼混沌系统的力学机理与能量转换,通过与Kolmogorov系统和欧拉方程的比较,把四翼混沌系统的矢量场分解为惯性力矩、内力矩、耗散和外力矩来探讨产生混沌的基本因素.利用Lie-Poisson括号揭示哈密顿能,动能与势能的相互转换.通过五种情形分析研究了四翼混沌吸引子不同类型力矩的功能和作用以及产生不同类型动力学模式的关键因素.此外,借助扩展的Kolmogorov系统, Pelino等[18]研究了洛伦兹系统的能量转换.上述几篇文献开启了从力学角度研究混沌系统的先河,但文献[15-17]中的混沌系统均为数值混沌系统(通过数值计算发现的混沌系统),没有明确的物理意义.借鉴这三篇文献的思想和方法我们探讨了具有明确物理意义的Couette-Taylor流问题,讨论的Couette-Taylor流三模混沌系统属于物理混沌系统(具有物理背景和物理意义的混沌系统),研究了Couette-Taylor流的动力机制,物理意义和能量转换.把Couette-Taylor流三模态系统(见文献[12])的状态变量视为角动量,将矢量场分解为惯性力矩,内转矩,耗散力矩和外力矩,不同的力矩的耦合视为不同的动力模式,通过分析各种力矩与Couette-Taylor流的不同模态的因果关系,揭示各种Taylor涡和混沌生成的力学机理.从力学角度分析同轴圆筒间旋转流动的各种涡流和混沌的生成机理目前还没有文献涉及.
由于Couette-Taylor流的全局吸引子是非常复杂和难于计算的. Couette-Taylor流系统的轨道从一个不稳定的平衡点移动到另一个不稳定的平衡点,相应地,轨道和平衡点之间的距离随参数和时间的不同而变化.我们发现,卡西米尔(Casimir)函数与距离有密切关系,关于洛伦茨吸引子的可预测性和能量可以通过能量转换来进行研究[18].因此,通过引入卡西米尔函数的导数,我们给出了Couette-Taylor流系统的能量转换,并且分析了它们的动力学和能量转换.
本文的结构如下:第2节将Couette-Taylor系统转换为Kolmogorov系统.第3节分析了不同动力模式的力学机理.在第4节,引入卡西米尔函数来获得混沌吸引子的边界,探讨它的动态行为.第5节按扩展Kolmogorov系统和能量转换,分析了Couette-Taylor流系统的力学机制,势能和动能间的转换和卡西米尔函数的关系.最后,第6部分是总结.
2 Couette-Taylor流三模系统及其Kolmogorov系统
文献[2]利用同轴圆筒间隙区域Stokes算子的特征函数作为基函数对周期性边界条件Navier-Stokes方程进行傅里叶展开,将原偏微分方程组化为常微分方程组,取傅里叶级数的前三项,得到下列方程组
这里
其中
变换之后系统(2.4)平衡点的数量和性质与原系统(2.1)的相同,但平衡点的位置发生了改变.容易获得系统(2.4)有如下三个平衡点:
定义哈密顿能量(the Hamiltonian energy)
其中
注1 变量
3 系统(2.5)的动力学机制及其分析
本节我们研究不同类型的力矩对系统(2.5)的影响,从而发现产生Taylor涡和混沌的关键因素.
情形1 系统只包含惯性力矩(由动能
动能函数的导数
图 2
情形2 系统包含哈密顿函数中的惯性力矩(由动能
哈密顿函数的导数
图 3
情形3 系统包含惯性力矩,内力矩和耗散力矩,即系统包含内能和耗散因素,但不包含驱动因素,相应的方程为
对于
其中V是系统相空间的体积.因此,此情形下系统是耗散的,即相空间的体积在流量下以指数形式收缩[20].哈密顿的能量变化率是
图 4
情形4 系统在惯性力矩,内力矩和外力矩下,即系统包含内能和驱动因素,但不包含耗散因素,相应的方程为
由于有了外力矩
图 5
图 6
为了分析和比较各种力矩对实际Couettt-Taylor流动所起的作用,我们讨论了上面虚拟的四种情形,下面的情形5才具有真实的物理意义.
情形5 系统在全部力矩下,即系统包含内能,同时包含耗散因素和驱动因素,相应的方程为
当
图 7
图 8
4 系统(2.4)的卡西米尔(Casimir)函数及其全局稳定性分析
混沌系统具有一些奇异的性质:如初始灵敏性和解的有界性.对于非混沌系统,正Lyapunov指数意味着系统解无限制地增长.然而,对于一个有界的混沌吸引子来说,其轨线增长的同时在不断的折叠,因此,解的边界性质在混沌系统研究中是至关重要的.通常很难找到混沌吸引子的边界,利用拉格朗日乘数法和卡西米尔函数法,我们得到了Couette-Taylor流系统混沌吸引子的边界.数值模拟表明,这是一个清晰的边界.卡西米尔函数
这意味着在Lie-Poisson括号下卡西米尔函数与每个函数可交换[19].对于Couette-Taylor流系统来说,卡西米尔函数定义为
由(2.3)式和文献[19],我们有
其中
令
显然
定理4.1 卡西米尔函数被限制在如下集合内
证 根据上述分析,
定义拉格朗日函数
其中
通过计算我们得到如下关于变量的导数
因此得到两个解
所以,
证毕.
图 9
这表示系统(2.5)的轨道与平衡点
图 10
图 11
根据以上理论分析和仿真结果,我们得出如下结论,当
表 1 σ=9.35,a=0.758,b=1.45,c=8.4,r取不同值时三模系统(2.1)的动力学行为与能量演化及其相应的Couette-Taylor流实际运动
r值范围 | 0<acr<σ0 r1=1.4685 | acr>σ re=1.7532 rg=27.5452 rh=31.6102 | |||
定点O | 稳定结点 | 鞍结点(一个方向不稳,另两个方向稳定) | |||
定点P+和P- | 不存在 | 稳定结点 | 稳定结点 | 稳定结点 | 鞍点 |
系统(2.4)相空间中的运动情况 | 趋于稳定定态O | 趋于稳定定态P+或P- | 运动最终按螺旋线趋于P+或P- | 同左,但越靠近rh,轨线在P+和P-之间来回跳动的越频繁,出现暂态混沌,最终趋于P+或P- | 不稳定极限环(亚临界霍普夫分岔) |
动能 | 最小值 | 逐渐增大 | 保持增大 | 增长 | |
Casimir函数 | 最小值 | 逐渐增大 | 保持增大 | 增长 | |
D1与D2之和 | 不存在 | 逐渐递增 | 递增 | 逐渐增大 | |
Coutte-Talor流的实际流动 | Couette流 | 形成规则的Talor涡流及Talor行进波等如图 1(a)(b)(c) | 经过波状涡流等达到暂态混沌如图 1(d) | 不规则湍流(混沌)如图 1(d) |
5 能量转换
其中
其中
根据方程(2.3)和(5.2),我们有
那么
其中
这意味着势能的变化率等于它与动能、耗散能量和外部能量的交换.
由于
然后, Couette-Taylor流系统的能量转换可以写成
这考虑了在
在上面的公式中,当
Couette-Taylor流系统有螺旋式轨道,它的轨道从一个不稳定的平衡点移动到另一个不稳定的平衡点.当轨道远离不稳定的平衡点之一时,
图 12
虽然Couette-Taylor流混沌系统与Lorenz系统不同,但能量转换时间演化与文献[18]的结果相似,能量、西米尔函数以及
6 结论
本文研究了Couette-Taylor流类洛伦兹系统的动力学机理和能量转换,通过理论和数值结果揭示了Couette-Taylor流问题的力学机制和物理意义.首先,探讨了Couette-Taylor流类洛伦兹系统作为柯尔莫哥洛夫系统的力学和物理意义.剖析了Couette-Taylor流混沌系统四种不同类型的力矩,研究了这些力矩的耦合的五种情况,探讨了产生泰勒涡和混沌的关键因素.在保守的情况下,哈密顿量是一个常数,相应的方程产生周期解.当耗散力矩或外力矩加入到保守系统中时,哈密顿函数趋近于零或无穷大,对应的系统不会产生混沌.当考虑到所有力矩时, Couette-Taylor流系统才产生混沌,并且动能和势能之间相互转换.对于Couette-Taylor流系统,内能、耗散因素和驱动因素并存是产生混沌的必要条件,而且,只有当耗散与驱动相匹配时,系统才能生成混沌.圆筒旋转的机械能为系统提供能量,增加的动能导致流动不稳定,依次生成Talor涡流、Talor行进波和波状涡流等经过暂态混沌,最终到达混沌.其次,本文将Couette-Taylor流混沌系统作为一个扩展的Kolmogorov系统研究了其力学和物理意义.引进卡西米尔函数来分析系统动力学和能量转换.两种不稳定平衡点的距离之和与卡西米尔函数有密切的关系,它们的时间演化是一致的.卡西米尔函数是内能,因此,它的变化速率是耗散与供给能量(外力矩)之间的交换速率.内能的导数起着能量转换的作用,平衡点距离的时间演化也有类似的动态.
Couette-Taylor流系统是耗散系统,其有界性是不平凡的,通过拉格朗日乘数法和卡西密尔函数法,得到了混沌吸引子的一个清晰的边界,数值仿真显示了方法的有效性.
参考文献
Lorenz systems for the incompressible flow between two concentric rotating cylinders
,DOI:10.4208/jpde.v23.n3.1 [本文引用: 1]
The chaos behavior and simulation of three model systems of Couette-Taylor flow
,
Destabilization of the Couette-Taylor flow by modulation of the inner cylinder rotation
,DOI:10.1016/j.euromechflu.2013.10.006
Optimal Taylor-Couette flow:Direct numerical simulations
,
Dynamical behaviors and numerical simulation of Lorenz systems for the incompressible flow between two concentric rotating cylinders
,
Stability of a viscous liquid contained between two rotating cylinders
,
Nonlinear dynamics of viscoelastic Taylor-Couette flow, effect of elasticity on pattern selection, molecular conformation and drag
,
Flow regimes in a circular Couette system with independent rotating cylinders
,
旋转流动的低模分析及仿真研究
,
Low-dimensional analysis and numerical simulation of rotating flow
A unified view of kolmogorov and lorenz systems
,DOI:10.1016/S0375-9601(00)00620-4 [本文引用: 2]
Mechanical analysis and energy conversion of Chen chaotic system
,
Mechanical analysis of Chen chaotic system
,DOI:10.1016/j.chaos.2017.03.021 [本文引用: 1]
Mechanical analysis of qi four-wing chaotic system
,
Energy cycle for the lorenz attractor
,DOI:10.1016/j.chaos.2013.09.005 [本文引用: 4]
/
〈 | 〉 |