2. 大连理工大学 能源与动力工程学院, 辽宁 大连 116024
2. School of Energy & Power Engineering, Dalian University of Technology, Dalian 116024, China
Taylor-Couette流动是流体力学和磁流体力学领域的经典流动.它是由内、外独立旋转的同心套筒构成, 两筒中间充满不可压缩的流体.
在单内筒旋转、外筒壁固定的圆筒中, 当内筒壁旋转速度较低时, 流动是稳定的, 这种一维流动称为旋转Couette流动; 当增加旋转速度达到某一临界速度后, 施加轴对称扰动后流场失去稳定, 形成Taylor涡流; 继续增加内筒转速后, 经历波状涡等状态后, 流场进入湍流状态[1-2].通过控制套筒的几何尺寸, 以及内、外筒壁面的旋转速度等条件, 可以实现一系列复杂的流动状态.由于该装置配置简单和其中富含的流动现象, Taylor-Couette流动被广泛应用于流动不稳定现象和湍流的研究中.同时, 由于Taylor-Couette流动具有极好的混合特性, 在很多工业中都有着广泛的应用[3-5].
在单内筒旋转的Taylor-Couette流动中, 当旋转速度足够大时, 流动进入湍流状态.此时, 瞬时流场中充满小尺度湍流涡结构, Taylor涡结构被扭曲和破坏, 变得难以识别; 尽管如此, Taylor涡仍然存在于平均流场中.Dong认为湍流Taylor-Couette流可以被看作是大尺度Taylor涡流和湍流脉动的叠加[6].Bilson等的研究也得到相似的结论, 他们通过两种不同的流场平均方法, 将湍流场中脉动来源归为平均Taylor涡的贡献和湍流贡献两部分, 证明了平均Taylor涡对湍流有重要贡献[7].杜珩等也使用类似的流场平均方法, 获得了不同几何尺寸的同心圆筒中的平均Taylor涡对湍流脉动的贡献[8].
导电流体Taylor-Couette流是流体力学一个新的研究热点.Chandrasekhar[9], Donnelly[10]等分别从线性稳定性理论和实验的角度出发, 发现轴向磁场能够增强Taylor-Couette流动稳定性并延迟Taylor涡的出现.目前, 磁场对湍流Taylor-Couette流动的研究还很少, 这是开展本研究的重要原因.
本文采用直接数值模拟方法, 研究无限长圆筒中的湍流Taylor-Couette流动, 通过两种不同的流场平均方法, 研究磁场作用下的平均Taylor涡流对湍流脉动贡献的变化规律.
1 物理问题与数学模型本文所研究的物理模型与坐标系统如图 1所示, 内筒半径为R1、以角速度为ω旋转; 外筒半径为R2, 保持静止; 在轴向(高度方向)上使用周期边界来模拟无限长圆筒, 计算区域高度为H.沿轴向施加强度为B的稳恒磁场.在实验室和一般工业应用中, 导电流体介质满足准静态假设, 即认为磁场是恒定不变的, 磁场会对流场产生作用而流场不会影响磁场[11].基于此, 使用两筒间隙d=R2-R1、内筒壁转速ωR1和磁感应强度B, 作为长度、速度和磁感应强度标尺, 可以得到无量纲控制方程组:
(1) |
(2) |
(3) |
(4) |
其中:U, J和e分别为柱坐标系下的无量纲速度、电流密度和磁感应强度矢量; t, p和φ分别为无量纲时间、压力和电势.描述流场几何形状参数为
无量纲控制参数为
其中:ρ, v和σ分别为导电流体的密度、黏度和电导率.
边界条件采用无滑移速度和等电势边界:
在内筒壁, uθ=1, φ=0, ur=uz=0;
在外筒壁, uθ=0, φ=0, ur=uz=0.
其中, r, θ和z分别为径向、周向和轴向坐标.
2 数值方法在强磁场作用下, 湍流行为将会产生许多独特的性质, 给数值计算带来更多的挑战.本文使用基于电流守恒格式的有限差分方法对控制方程进行离散并求解[12-13].相较于非守恒格式, 本文采用的守恒格式有更好的数值稳定性, 并能够降低非守恒格式引起的误差.
2.1 数值求解过程控制方程的时间离散采用二阶精度、显式的Adams-Bashforth格式和Backward-differentiation格式.这里使用投影算法完成速度压力解耦, 并完成时间步长的推进.在每个时间步中(第n步):
1) 求解电势Poisson方程, 并计算电流密度:
(5) |
(6) |
2) 求解预估速度:
(7) |
3) 求解压力Poisson方程:
(8) |
4) 修正预估速度场, 得到第n+1步的速度,
(9) |
关于数值算法的详细描述参考文献[13].
2.2 算法验证本文研究中, 计算对象为半径比为0.5, 高度为2π的同心圆筒; 雷诺数Re = 4 000, 哈特曼数Ha=0, 20, 40和80.这里选择的计算网格在轴向、周向和径向分别为320, 320和160.
为了验证网格和算法的正确性, 这里首先计算未施加磁场时的湍流Taylor-Couette流动.在Re=8 000的验证计算中, 内、外筒壁的无量纲转矩均为0.013 2, 与文献[6]中的结果0.012 8十分接近.随后, 在图 2中又绘制了Re=4 000和Re=8 000算例中的平均速度分布曲线和文献[6]中的分布曲线.对比之后不难发现, 本文的结果与文献中的结果一致, 说明本文采用的计算方法和网格是可靠的.
本文首先研究平均速度的变化.图 3是施加不同强度磁场之后的平均速度分布.这里最为明显的变化趋势是, 在外筒壁附近的流场中产生了速度为负值的分布; 而且负值的大小随磁场强度的增加而增大.说明在轴向磁场的作用下, 流场中产生了反向的电磁力.由于内筒壁附近速度较大, 惯性力能够克服电磁力的作用而不出现负值, 但仍然可以看到速度降低、梯度增大的趋势; 而外筒壁附近的流动较弱, 反向电磁力抑制正向速度, 并随磁场增大逐渐将流体反向加速.
图 4列出了Ha=0和Ha=80的平均速度矢量图.在每幅图中, 左侧为内筒壁, 右侧为外筒壁.在未施加磁场时, 平均流场中出现三对大小一致的Taylor涡结构, 验证了文献[6-7]中的结论; 而在Ha=80后, 大尺度的涡结构几乎完全消失, 只有内筒壁附近存在少量涡结构.说明磁场能够破坏平均流动; 外筒壁附近区域中, 磁场对流动的抑制作用更强.
下面分析湍流脉动的变化.按照文献[7]中的计算方法, 全部脉动的计算是将流场在周向、轴向和时间方向进行平均, 其均方根脉动的表达式为
(10) |
对于三维流场, 如果过滤掉二维的平均Taylor涡, 就可以获得湍流引起的脉动, 相应的均方根脉动的表达式为
(11) |
将全部脉动减去湍流对脉动的贡献U′-U″, 即为平均Taylor涡对湍流脉动的贡献.
按照这种方法, 类似地, 可以计算出全部湍动能E′tke、湍流对湍动能的贡献E″tke和平均Taylor涡对湍动能的贡献E′tke-E″tke, 如图 5中的分布曲线所示.
在图 5a中, 当Ha=0时, 由于湍流边界层的存在, 两侧边界层中湍动能较大; 特别是在内筒边界层中, 由于内筒的驱动作用, 湍动能更大.增加磁场强度后, 内筒边界层中的湍动能在Ha=20几乎不变化; 而在Ha=40后, 湍动能开始降低.
而在流场中间部位和外筒壁附近, 施加磁场之后, 湍动能大幅降低; 在Ha=40时, 外筒壁附近的湍动能小幅度增加; 最后, 在Ha=80时, 湍动能仍然被抑制.在此过程中, 外筒壁附近湍动能增加的原因, 有可能是反向电磁力在驱动流体形成反向流动的时候形成的.在整个半径方向上, 内筒壁的湍动能始终比外筒壁大.
在图 5b中的湍流对湍动能的贡献中, 当Ha=0时, 内筒壁附近的湍动能比其他区域中的湍动能略大; 而且在Ha=20和Ha=40时, 湍动能大幅度增加, 说明此时该区域中的湍流强度得到局部增强.引起这一现象发生可能的原因是, 反向电磁力增加内壁边界层中的平均剪切率, 湍动能生成量相应增强.在外筒壁附近, 湍动能变化与图 5a中的变化趋势几乎完全一样.
在图 5c中可以看出, 在没有施加磁场时, 平均Taylor涡对湍动能的贡献很大, 约占50%;而在施加磁场之后, 这部分贡献迅速降低为零.说明轴向磁场破坏了湍流中的平均Taylor涡流, 使其对湍动能的贡献大幅降低.这一点也在图 4的矢量图中得到证实.
4 结论1) 在轴向磁场作用下, 等电势的旋转筒壁会产生径向感生电流, 产生电磁力并会改变平均流动的分布.
2) 在本文选择的磁场范围内, 轴向磁场对湍流脉动有明显的抑制作用, 在外筒壁区域抑制作用更加强烈.
3) 随着磁场的增加, 磁场通过破坏平均Taylor涡流, 大幅度抑制平均流动对湍动能贡献; 而对于湍流对湍动能的贡献, 磁场对其抑制作用相对较弱.
致谢: 感谢德国伊尔默瑙工业大学Yuri Kolesnikov教授和Dmitry Krasnov博士的帮助和指导.[1] |
Taylor G I.
Stability of a viscous liquid contained between two rotating cylinders[J]. Philosophical Transactions of the Royal Society A, 1923, 223: 289–343.
DOI:10.1098/rsta.1923.0008 |
[2] |
Swinney H L, Gollub J P.
Hydrodynamics instabilities and transition to turbulence[M]. Berlin: Springer-Verlag, 1981: 139-180.
|
[3] |
Pierson J A.
Photolytic reactors:similitude in Taylor-Couette and channel flows[J]. AIChE Journal, 2003, 49(5): 1285–1292.
DOI:10.1002/(ISSN)1547-5905 |
[4] |
Paule A, Lauga B, Ten-Hage L, et al.
A photosynthetic rotating annular bioreactor(Taylor-Couette type flow)for phototrophic biofilm cultures[J]. Water Research, 2011, 45: 6107–6118.
DOI:10.1016/j.watres.2011.09.007 |
[5] |
Sengupta T K, Kabir M F, Ray A K.
A Taylor vortex photocatalytic reactor for water purification[J]. Industrial & Engineering Chemistry Research, 2001, 40(23): 5268–5281.
|
[6] |
Dong S.
Direct numerical simulation of turbulent Taylor-Couette flow[J]. Journal of Fluid Mechanics, 2007, 587: 373–393.
|
[7] |
Bilson M, Bremhorst K.
Direct numerical simulation of turbulent Taylor-Couette flow[J]. Journal of Fluid Mechanics, 2007, 579: 227–270.
DOI:10.1017/S0022112007004971 |
[8] |
杜珩, 阙夏, 刘难生.
不同半径比Taylor-Couette湍流的直接数值模拟研究[J]. 中国科学技术大学学报, 2014, 44(9): 761–768.
( Du Heng, Que Xia, Liu Nan-sheng. Influence of radial ratio on Taylor-Couette turbulent flow using direct numerical simulation[J]. Journal of University of Science and Technology of China, 2014, 44(9): 761–768. DOI:10.3969/j.issn.0253-2778.2014.09.010 ) |
[9] |
Chandrasekhar S.
Hydrodynamics and hydromagnetic stability[M]. Oxford: Oxford University Press, 1961: 343-426.
|
[10] |
Donnelly R J, Ozima M.
Experiments on the stability of flow between rotating cylinders in the presence of a magnetic field[J]. Proceedings of the Royal Society of London, Series A, 1962, 226: 272–286.
|
[11] |
Davidson P A.
An introduction to magnetohydrodynamics[M]. Cambridge: Cambridge University Press, 2001: 117-119.
|
[12] |
Ni M J, Munipalli R, Morley N B, et al.
A current density conservative scheme for incompressible MHD flows at a low magnetic Reynolds number.part Ⅰ:on a rectangular collocated grid system[J]. Journal of Computational Physics, 2007, 227: 174–204.
DOI:10.1016/j.jcp.2007.07.025 |
[13] |
Krasnov D, Zikanov O, Boeck T.
Comparative study of the finite difference approaches in simulation of magnetohy drodynamic turbulence at low magnetic Reynolds number[J]. Computers & Fluids, 2011, 50: 46–59.
|