LARGE EDDY SIMULATION OF FLOW PAST TWO CONICAL CYLINDERS IN TANDEM ARRANGEMENT
-
摘要: 利用大涡模拟研究了雷诺数Re = 3900下串列双锥柱在间距比L/Dm = 2 ~ 10下的升阻力特性及三维流动结构. 研究发现: 上游锥柱在后方形成的两个展向不对称回流区, 使其后方压力分布不对称. 上游锥柱发展的上洗、下洗和侧面剪切层作用在下游锥柱的附着点位置不同是上游和下游锥柱时均阻力系数和脉动升力系数变化的主要原因, 串列双锥柱间流动结构随间距比变化可分为三种状态: 剪切层包裹状态, 过渡状态及尾流撞击状态. 剪切层包裹状态. 上游锥柱的自由端主导来流在下游锥柱迎风面影响范围广, 上游锥柱剪切层完全包裹住下游锥柱, 从而抑制下游锥柱后方回流区形成, 导致下游锥柱时均阻力系数降低; 尾流撞击状态; 上游锥柱尾流得到充分发展, 其回流区大小随间距比增大不再发生变化, 上游锥柱尾流出现周期性脱落, 撞击在下游锥柱表面, 从而使脉动升力系数大幅增加, 最大脉动升力系数较单直圆柱提升约20.7倍; 过渡状态, 此时时均阻力系数和脉动升力系数均会较剪切层包裹状态增加. 该研究可以为风力俘能结构群列阵布局提供理论支持.Abstract: In order to explore the time-averaged drag coefficient, fluctuating lift coefficient characteristics and flow field mechanism of two conical cylinders in tandem arrangement, large eddy simulation is used to simulate two conical cylinders in tandem arrangement with a spacing ratio of 2−10 at Re = 3900. The two spanwise asymmetric reflux zones formed behind the upstream conical cylinder make the en dash pressure distribution behind it asymmetric. The upwash, downwash and side shear layers developed by the upstream conical cylinder are the main reasons for the variation of the time-averaged drag coefficient and the fluctuating lift coefficient of the upstream and downstream conical cylinders. The flow structure between tandem two conical cylinders can be divided into three states with the change of spacing ratio: in the shear layer wrapping state, transition state and wake impact state. Shear layer wrapping state, the dominant incoming flow at the free ends of the upstream conical cylinder has a wide range of influence on the windward side of the downstream conical cylinder. The shear layer of the upstream conical cylinder completely wraps the downstream conical cylinder, inhibiting the formation of the backflow zone behind the downstream conical cylinder, causing a decrease in the time-averaged drag coefficient of the downstream conical cylinder. In the wake impact state, the wake of the upstream conical cylinder is fully developed, and the size of the recirculation zone does not change with the spacing ratio. The wake of the upstream conical cylinder periodically falls off and hits the surface of the downstream conical cylinder, which greatly increases the pulsating lift coefficient. The maximum fluctuating lift coefficient is about 20.7 times higher than that of a single straight cylinder. In the transition state, the time-averaged drag coefficient and the fluctuating lift coefficient will both increase compared with the shear layer wrapped state. This research can provide theoretical support for the layout of wind energy harvesting structures.
-
引 言
风能作为一种重要的可再生清洁能源, 对其开发利用受到国内外学者的广泛关注[1-2]. 对于风能的获取, 如涡激振动发电[3]、风力压电效应发电[4], 增振就尤其重要.
圆柱作为绕流的经典研究对象, 通常通过改变圆柱几何形状来改变周围流动状态以实现流动控制, Zhang等[5]通过改变波浪外形圆柱的波长比以及波幅比, 发现波幅比为0.2时波浪圆柱受到的时均阻力系数显著降低. Zhang等[6]通过对比研究无限长和一端自由的圆柱, 发现有限长圆柱的时均阻力系数和脉动升力系数远小于无限长圆柱. 赵桂欣[7]通过对比四种不同端面以及侧面的两端自由的有限长圆柱, 发现改变外形对柱体阻力系数均有降低作用, 在特定雷诺数下, 脉动升力系数较单直圆柱有所提升. 赵萌等[8]、乔永亮[9]的研究发现有限长圆柱长径比会明显改变圆柱的流场结构, 存在一个临界长径比使得有限长圆柱的流动模式由自由端来流主导的无序涡脱模式转向侧面来流主导的卡门涡街涡脱模式. 杨耀宗[10]的研究发现, 有限长锥柱在斜率k = 0.05时锥柱振动性能有显著提升. 由此可见, 对于有限长圆柱, 通过改变其表面形状对脉动升力系数有较大提升.
双圆柱较单圆柱的流动模式更为复杂. Alam等[11-12]通过实验, 研究不同布置形式的双圆柱, 对两个圆柱之间的尾流形态进行分类, 研究发现在间距比为2.4 ~ 3.0, 交错角为10°时的错列双圆锥脉动阻力出现最大值. Zdravkovich[13]研究两个串列圆柱间距比小于5的情况, 发现间距比小于3时, 上、下游圆柱之间不会产生明显涡脱, 在间距比大于3时, 上下游圆柱之间逐渐产生涡脱. 杜晓庆等[14-16]对不同布置状态下的双圆柱进行了数值仿真, 按照上、下游圆柱相互作用的不同形式将流态进行分类, 解释了下游圆柱在上游圆柱尾流的不同状态下的升阻力特性, 发现下游圆柱被上游包裹状态下, 下游圆柱的阻力系数显著降低, 脉动升力系数反而受到一定抑制, 但是相对于单圆柱, 下游圆柱脉动升力系数均有不同程度增幅. 文献[17-19]的研究发现, 串列双圆柱随间距比增加, 脉动升力系数和时均阻力系数先是缓慢降低, 但是在间距比为3.0 ~ 4.0附近发生突增跳跃到一个较高的水平, 上游圆柱对下游圆柱的影响由剪切层包裹的抑制状态逐渐变为涡脱撞击状态, 使得下游圆柱在超过临界间距比后脉动升力相对于单圆柱大幅提升. 对于改变圆柱外形的串列双圆柱, 唐涛等[20]数值模拟了间距比为3的不同长度比的梯形截面串列柱, 发现两圆柱之间的间隙流会随着梯形长度比的改变而呈现不同的耦合模式, 导致升力发生改变. 因此上游圆柱的尾流结构对下游圆柱脉动升力的提升有着决定性的作用, 但是对于有限长串列双锥柱, 其流动模式和升阻力特性缺乏研究.
综上所述, 下游柱体脉动升力系数有可能相对于单柱体有大幅提升, 同时表面的形状改变对于柱体脉动升力也有很大影响. 根据课题组前期研究, 发现锥形圆柱在斜率k = 0.05时相对于直圆柱增振最明显[10], 因此, 本文将以提升脉动升力系数为目的, 利用商用软件Fluent中的大涡模拟, 在亚临界雷诺数(Re = 3900)下, 对斜率k = 0.05的串列双锥柱在间距比L/Dm = 2 ~ 10下进行数值模拟, 通过分析上、下游锥柱之间的流动机理, 探索下游锥柱时均阻力系数和脉动升力系数变化规律, 为风力俘能结构群的列阵布局提供理论支持.
1. 数值方法
1.1 控制方程
大涡模拟(large eddy simulation, LES)[21-22]采用空间滤波技术, 通过截止尺度和滤波函数将涡识别为大涡和小涡, 大涡直接解析, 小涡则通过亚格子尺度应力模型(sub-grid-scale stress, SGS)封闭求解, 过滤后的不可压缩黏性流体Navier-Stokes方程为
$$ \qquad \frac{{\partial {{\bar u}_i}}}{{\partial {{\bar x}_i}}} = 0 $$ (1) $$ \qquad \frac{{\partial \overline {{u_i}} }}{{\partial t}} + \frac{{\partial \overline {{u_i}{u_j}} }}{{\partial {x_j}}} = - \frac{1}{\rho }\frac{{\partial \overline p }}{{\partial {x_i}}} + v\frac{{{\partial ^2}\overline {{u_i}} }}{{\partial {x_j}{x_j}}} - \frac{{\partial {\tau _i}_j}}{{\partial {x_j}}} $$ (2) 式中, xi, yi为笛卡儿坐标, i, j = x, y, z, ui, uj为滤波速度矢量,
$ \overline p $ 为流体所受压力, ρ为流体密度, ν为流体运动黏度, τij为亚格力应力.根据Smagorinsky[23]的基础SGS模型, 设应力表达式
$$ {\tau }_{ij}-\frac{{\tau }_{kk}{\delta }_{ij}}{3}=-2{u}_{t}\overline{{S}_{ij}} $$ (3) 式中,
$ {\tau _{kk}} $ 是SGS应力各向同性部分,$ {\delta _{ij}} $ 为Kronecker符号,$ {u_t}{\text{ = }}{\left( {{C_S}\varDelta } \right)^2}\left| {\overline S } \right| $ 为亚格子尺度湍流黏度$$ \overline {{S_{ij}}} {\text{ = }}\frac{1}{2}\left( {\frac{{\partial {{\overline u }_i}}}{{\partial {x_j}}} + \frac{{\partial {{\overline u }_j}}}{{\partial {x_i}}}} \right) $$ (4) $$ \left|\overline{S}\right|=\sqrt{2\overline{{S}_{ij}}\overline{{S}_{ij}}} $$ (5) $$ \varDelta {\text{ = }}{\left( {\varDelta _x\varDelta _y\varDelta _z} \right)^{1/3}} $$ (6) 式中,
${\varDelta _i}$ 为网格在i方向的尺度, Cs为Smagorinsky尺度, 本文取0.2[24].1.2 控制方程离散
本文计算控制方程离散均采用有限体积法, 压力速度耦合方式为PISO, 离散化均采用二阶格式[6].
2. 计算域及网格
2.1 模型简介
基于多数人研究的直圆柱, 本文加入斜率k, 将直圆柱转化为锥柱, 锥柱如图1所示, 几何表达式如下
$$ {D}_{z}={D}_{m} + 2kz $$ (7) $$ {D_m} = \frac{{{D_{\max }} + {D_{\min }}}}{2} $$ (8) 式中, Dz为锥柱在高度为z处的横截面直径, Dm = 0.01 m为锥柱的平均直径, k为锥柱的斜率, 圆柱高度比H/Dm = 7, Dmax和Dmin分别为锥柱最大直径和最小直径.
本文中出现的相关参数定义如表1.
表 1 相关参数定义Table 1. Parameter definitionParameter Remark Parameter Remark Dm 0.01 m, average diameter of conical cylinder $\nu $ 1.48 × 10−5 m2/s, kinematic viscosity H 7Dm, the length of the conical cylinder ρ 1.225 kg/m3, density of air L center distance between two conical cylinders U, U∞ 5.772 m/s, velocity of inlet Δy height of the first level grid Re 3900, Reynolds number Cd 2Fd/(ρU2HDm), drag coefficient, Fd is the total
drag forceCp (p−p0)/(0.5ρU2), pressure coefficient, p, p0 are the static
pressure and reference pressure, respectivelyCl 2Fl/(ρU2HDm), lift coefficient, Fl is the total lift force Nfe number of free ends Cdmean time-averaged drag coefficient Clrms RMS value of lift coefficient 2.2 计算域及网格划分
本文研究的串列双锥柱计算域大小为(30Dm + L) × 20Dm × 10Dm, X轴正方向为顺流向, Y轴为横流向, Z轴为圆柱展向, 笛卡尔坐标系圆心位于上游锥柱中心处. 上游锥柱中心距离计算域入口10Dm, 下游锥柱中心距离计算域出口20Dm, 锥柱的上下端面距离计算域顶部和底部均为1.5Dm, 如图2所示. 入口采用速度入口(velocity-inlet), U∞ = 5.772 m/s, 出口采用压力出口(pressure-outlet), 背压设置为0 Pa, 计算域两个侧面以及顶面采用对称边界(symmetry), 计算域底面和圆柱体表面为无滑移壁面(no-slip-wall). 流体介质为空气, 密度ρ为1.225 kg/m3, 运动黏性系数ν为1.48 × 10−5 m2/s.
本文使用ICEM对计算域进行结构化网格划分, 利用O-block技术对圆柱周围以及顶部网格进行加密, 同时避免了片状低质量网格的生成. 大涡模拟对于边界层的要求控制y+ ≤1, 每个圆柱周围有3Dm × 3Dm指数形式增长的加密区域, 同时也对两个圆柱之间以及下游锥柱尾流区的网格进行了加密, 如图2所示.
2.3 计算模型验证
本文首先验证使用单直圆柱计算模型, 探究时间步长Δt、网格参数(圆周节点数, 第一层边界层高度)对圆柱的时均阻力系数(Cdmean), 脉动升力系数(Clrms)和Strouhal数(St)的影响, 计算结果如表2. 综合考虑计算成本和计算精度的影响, 本文选用Case2的网格精度和时间步长进行下面仿真.
表 2 时间步长、圆周节点数和第一层高度对时均阻力系数, 脉动升力系数和Strouhal数的影响Table 2. The influence of time step and grid parameters on the Cdmean , Clrms and StCase Gird number Δy/D Nfe Δt Re H/Dm Cdmean Clrms St Case1 1 035 424 0.001 2 0.0002 3 900 7.0 0.749 0.0151 0.184 Case2 2 862 840 0.001 2 0.0002 3 900 7.0 0.756 0.0123 0.167 Case3 5 407 524 0.001 2 0.0002 3 900 7.0 0.758 0.0121 0.164 Case4 2 956 920 0.0005 2 0.0002 3 900 7.0 0.759 0.0119 0.161 Case5 2 709 560 0.002 2 0.0002 3 900 7.0 0.741 0.0112 0.154 Case6 2 862 840 0.001 2 0.00002 3 900 7.0 0.760 0.0126 0.164 Case7 2 862 840 0.001 2 0.002 3 900 7.0 0.743 0.0110 0.115 num.[25] — — 1 — 3 900 6.0 0.932 0.0190 — exp.[26] — — 2 — 88 000 5.0 0.742 — — 由于自由端数目以及长径比不同, 本文单圆柱数值结果与仿真数据(num.)[25]以及实验数据(exp.)[26]有一定误差, 为了进一步验证本次仿真所用的算法, 网格精度以及时间步长能够适用于串列双锥柱, 对串列双锥柱的流动结构进行烟线实验, 其实验装置如图3所示. 风洞为直流式风洞(CSU-WT2450-40), 烟线控制器(SW-2) 通过充电放电加热涂有石蜡油的铜丝, 石蜡油在瞬间高温下产生白色烟雾, 风洞试验段均匀来流经过铜丝时, 白烟在风的作用下绕过固定在风洞低端串列布置的锥柱, 通过相机(Canon EOS 850 D)捕捉串列双波浪锥柱之间白烟形成的流场. 对串列双波浪锥柱四个截面Y = 0, S1(z/H = 3/4), S2(z = 0), S3(z/H = −3/4)与同截面下仿真结果进行对比.
实验结果如图4所示, 串列双锥柱实现显示, Y = 0截面可以看到, 间距比L/Dm = 5的串列双锥柱两个自由端来流撞击在下游波浪锥柱迎风面, 这与仿真的结论一致, 观察S1, S2, S3截面发现, 实验显示锥柱后方回流区发展完成, 剪切层已经收缩, 尾流撞击在下游锥柱表面, 仿真结果与实验结论基本一致. 因此可以认为本文验证中Case2的算法, 网格精度及时间步长可用于串列双锥柱仿真研究.
3. 结果分析
3.1 时均阻力系数和脉动升力系数
时均阻力系数和脉动升力系数是风力俘能结构获取风能的两个重要参数. 本文与文献[18]上游柱和下游柱及Case2中单圆柱数据对比结果, 如图5所示, 其中“UC”表示上游柱, “DC”表示下游柱.
本文上、下游锥柱与文献中上、下游圆柱的Cdmean和Clrms随间距比变化规律基本一致, 但是数值大小存在很大差异, 可能是因为本文锥柱自由端的影响, 下文将通过流场结构分析其原因. 本研究发现, Cdmean和Clrms在L/Dm = 4 ~ 5时存在突变, 在这个区间可能发生流动状态的改变; 在临界间距比之前上、下游锥柱Clrms呈下降趋势, 而Cdmean则随间距比先增加后减小; 在间距比超过临界值之后, 上游锥柱Clrms在临界值之后基本保持不变, 且与单直圆柱基本一致. 下游锥柱Clrms在间距比L/Dm = 5 ~ 6有一个略微降低然后突增的过程, 在L/Dm = 7的时候出现峰值, 此时较单直圆柱增长约20.7倍,此阶段可能存在流动结构改变. 上游锥柱Cdmean在低于单直圆柱一定值处波动, 下游锥柱Cdmean由一个较小值逐渐增大, 趋近于上游锥柱, 下游锥柱相对于单直圆柱Cdmean最大减小约90.7%, 在L/Dm = 7时, 减小约19.8%.
3.2 表面压力系数
表面压力系数Cp能够直接反映出锥柱表面受力情况, 为了方便观察表面Cp分布情况, 如图6(a)所示, 将锥柱的表面沿母线剪开, 并以角度“θ”和高度“Z/H”作为横、纵坐标铺平, 表面时均Cp分布如图6(b)所示.
从整体来看, 无论上游锥柱还是下游锥柱在不同间距比下的时均Cp分布基本关于180°对称, 说明来流对上、下游锥柱两侧的影响在时间平均上基本相同. 根据云图分布的不同, 可将上游圆柱的Cp分为两类, 第一类为低压区(压力系数小于−0.6的区域)集中在锥柱顶部(DIS1), 如L/Dm = 2, 3, 5, 第二类分为高压区(压力系数大于0.6的区域)主要集中在锥柱底部(DIS2), 如L/Dm = 4, 6, 7, 8, 9. 可以看到上游锥柱除L/Dm = 2, 3, 4, 低压区强度相对较大以外, 其他间距比分布基本相近, 这是Cdmean在临界区域前端逐渐降低的主要原因; 由于下游锥柱受到上游锥柱尾流的影响, 在小间距比下压力分布强度较弱, 随着间距比增加, 高压区和低压区范围和强度都有所增加, 并且分布形式逐渐接近上游圆锥, 这也是下游锥柱阻力系数随间距比增加而增大的主要原因(临界区域除外). 锥柱的表面压力分布主要是由其流场决定, 下面从流场的本质对压力分布进行分析.
3.3 时均流线图
三维流场较为复杂, 但串列双锥柱的流场分布具有一定对称性, 下面取Y/Dm = 0, S1, S2, S3截面(见图3(b))对串列双锥柱的时均流线进行分析, 如图7所示.
由于自由端的存在, 锥柱上、下端部流体分别存在上洗(流经下自由端的空气向上运动)和下洗(流经上自由端的空气向下运动)现象, 上洗、下洗作用会抑制侧面来流在锥柱后方形成周期性涡脱落, 从而导致脉动升力系数降低, 这也是串列双锥柱脉动升力远小于串列无限长圆柱的主要原因[27-28]. 空气流经上游锥柱在其后方形成一对Y截面的展向回流区(arch-vortex), 同时空气流过柱体侧面,与arch-vortex相互作用, 在柱体后方形成Z截面的横向回流区(side-vortex, 以下简称回流区), 除此之外, 在柱体端面处能观察到顶涡(tip-vortex)[29-30], 如图7(g). 课题组前期研究发现[10], 由于锥柱斜率的存在, 导致自由端上洗作用强于下洗, 下端形成的arch-vortex挤压上端形成的arch-vortex使其变为椭圆形, 对于串列双锥柱, 在L/Dm ≥ 7时, 上游锥柱后方能够明显观察到一对arch-vortex, 但在间距比L/Dm = 2 ~ 6时, 受到下游锥柱影响, 上游锥柱自由端形成的涡结构发生明显变化, 根据上洗和下洗在两个圆柱之间形成的arch-vortex不同, 将两锥之间的流场形式分为两类, 第一类是由下洗作用占主导, 流体经上自由端向下运动一部分锥柱上端形成一个涡, 另一部分与下自由端的上洗流相作用在锥柱下端形成一个涡, 如L/Dm = 2, 3, 5, 第二类则与第一类相反, 如L/Dm = 4, 6, 7, 8, 9. 这两种分类方式正好也解释上一节中DIS1和DIS2的分布规律, 占主导作用的一端, 空气流速更快, 上游锥柱背风面对应位置压力也就越低. 而无限长串列圆柱上游柱背风面不存在自由端影响, 背风面形成的回流区主要来自侧面, 压力系数展向对称分布[6]. 串列双锥柱间流动结构随间距比变化可分为三种状态: 剪切层包裹状态, 过渡状态, 尾流撞击状态. 剪切层包裹状态, 上游锥柱回流区受下游锥柱影响并未完全发展, 剪切层完全包裹住下游锥柱,如图7中L/Dm = 2, 3; 尾流撞击状态, 上游锥柱当前位置回流区已经完全发展并且大小随间距比基本不发生改变, 且下游锥柱已位于回流区之外, 不再被上游剪切层包裹, 如图7中间距比L/Dm = 7, 8, 9; 过渡状态, 即上游锥柱自由端来流作用到下游锥柱表面, 而侧面剪切层部分包裹下游锥柱, 如图7中间距比L/Dm = 4, 5, 6. 随着三种状态变化, 上游锥柱的上洗和下洗直接作用在下游锥柱表面的范围逐渐减小, 同时下游锥柱后方的回流区随着间距比增加而逐渐扩大, 这也是下游锥柱随间距比增加, 其高压区和低压区范围逐渐扩大并增强的主要原因.
3.4 涡结构图
涡的识别方法可以分为三代[31], 其中第二代中的Q方法运用较为广泛, 表达式如下
$$ Q = \frac{1}{2}\left( {{{\left\| {\boldsymbol{B}} \right\|}_{{\text{Fro}}}}^2 - {{\left\|{\boldsymbol{ A}} \right\|}_{{\text{Fro}}}}^2} \right) $$ (9) 式中, A为对称张量, B为反对称张量,
${\left\| {\;\;} \right\|_{{\rm{Fro}}}}^2$ 为矩阵的Frobenius范数. 下面对串列双锥柱的涡采用Q方法识别, 如图8所示, 每个间距比提供两个视角view1, view2.从涡量图中可以看到, 前面提到的两种分类, DIS1在这里显示的是大量涡“下泄”集中在两个锥柱之间的中下部, 如图8(a)(b)(d), DIS2在这里反应的是大量涡“上泻”集中在两个锥柱之间的中上部, 如图8(c)~图8(h). 从涡的结构上来看, 间距比较小时, 上游锥柱脱落的涡碎且无序, 与下游圆柱产生的涡相互作用, 在下游圆柱后方形成“肋状涡”(rib vortex), 如图8(a)中所示, 当间距比逐渐增大时, 上游锥柱脱落的涡逐渐变大并且有周期特性, 下游锥柱后方“肋状涡”也逐渐变大变长, 如图8(e)~图8(h)所示. 而对于无限长圆柱, 不同间距比下均能观察到卡门涡街现象[19], 所以其脉动升力系数能明显高于本文有限长双串列锥柱. 上游锥柱脱落的尾涡作用在下游锥柱表面, 与来自下游锥柱端部以及侧面来流相互作用, 产生较大的升力, 相较于低间距比串列双锥柱, 当间距比为L/Dm = 7 ~ 9时, 下游锥柱处于上游锥柱回流区临界位置或回流区下游, 上游锥柱尾流能够充分发展, 形成类似卡门涡街特征, 周期性的漩涡作用在下游锥柱表面, 使下游圆锥表面更容易产生大的升力. 间距比L/Dm = 7时, 上游锥柱尾流发展最充分, 其尾流对下游锥柱的作用更为强烈, 这可能是下游锥柱Clrms相对于间距比L/Dm = 8, 9更大的原因. 下面通过分析几个截面的涡量, 来具体分析上下游锥柱之间的流动情况.
3.5 Z截面涡量图
为了方便观察串列双锥柱之间的流动形态, 下面取S1, S2, S3(见图3(b))三个截面的涡量图, 如图9所示. 尾流之间的干扰可分为很多不同形式, 根据锥柱截面直径以及流场特性, 参照文献[12-13]对流场的分析, 将Z截面涡量分为两类, V1: 上游锥柱尾流包裹下游锥柱, 如图9中双点划红线外框中涡量图, 上、下游锥柱共同形成的剪切层极不稳定, 在下游锥柱后方形成的回流区相较于单锥柱[8]很小, 在后方形成非周期性的涡脱; V2: 下游锥柱位于上游锥柱回流区外端或恰好位于其回流区边界处, 如图7中虚线黑线外框中涡量图, 上游锥柱尾涡得到充分发展, 尾流撞击在下游锥柱表面, 结合下游锥柱产生的涡在下游锥柱后方产生大量不规则碎涡, 相较于V1为下游锥柱提供更大的升力. 结合图8中涡量图显示, 上游锥柱的上洗和下洗作用强度不同, 使得作用在下游锥柱前端的涡集中位置不同, 如图8(f), 上游锥柱尾涡集中在锥柱中下部, 导致图9中S1截面两个圆柱相互作用较弱,这也是图6中间距比L/Dm = 4的下游圆柱顶部时均压力系数高的主要原因. 图9中, 当间距比L/Dm = 2, 3时, 三个截面流动模式均为V1, 即4.3节提到的剪切层包裹状态, 间距比L/Dm = 4, 5, 6中, 三个截面中V1, V2模式共同存在, 即过渡状态, 而当间距比L/Dm = 7, 8, 9时, 三个截面流动模式全部变为V2, 即尾流撞击状态, 这是间距比L/Dm = 7, 8, 9有更高的脉动升力值主要原因. 从图9中可以观察到, 图9(f)相对于9(g)和9(h)每个截面下游锥柱距离上游锥柱回流区更近, 上游锥柱发展的周期性尾涡作用在下游锥柱表面更剧烈, 这是下游锥柱在L/Dm = 7取得Clrms幅值的主要原因.
4. 结 论
本文在Re = 3900下利用大涡模拟对间距比为L/Dm = 2 ~ 10的串列双锥柱进行了固定绕流仿真. 首先, 验证了Case2时间步长以及网格达到收敛, 通过烟线实验证明本次网格精度, 时间步长和算法可用于串列双锥柱仿真; 其次, 分析了其升阻力系数随间距比变化的特性; 再次, 分析了其时均流场信息分析其阻力变化的主要原因; 最后, 通过涡量图分析了脉动升力系数变化的主要原因, 得出以下结论.
(1)随间距比增加, 上游锥柱Clrms先减小逐渐趋近于单直圆柱, 在L/Dm = 4之后基本保持不变; 上游锥柱Cdmean先减小, 在L/Dm = 4之后保持低于单直圆柱的波动状态.
(2)下游锥柱Clrms随间距比基本保持增加后减小趋势(除临界区间), 在L/Dm = 7时, 取得峰值, 约为单直圆柱的21.7倍; 下游锥柱Cdmean由一个较小值随间距比增加逐渐接近上游锥柱的Cdmean, 下游锥柱相对于单直圆柱Cdmean最大减小约90.7%, 在L/Dm = 7时, 减小约19.8%.
(3)相较于无限长串列双圆柱, 串列双锥柱的上游锥柱自由端的上洗和下洗作用能显著影响侧面来流在锥柱后方形成的涡脱的发展, 抑制其产生规则卡门涡街现象; 上洗和下洗作用强的一端会使上游锥柱背风面形成大且强的低压区, 同时使同端下游锥柱迎风面的高压区减小.
(4)随间距比增加, 下游锥柱与上游锥柱回流区的关系有剪切层包裹状态、过渡状态、尾流撞击状态, 其中尾流撞击状态时, 上游锥柱尾流能够充分发展, 形成交替脱落的肋状涡, 撞击在下游锥柱表面, 使下游锥柱产生较大的脉动升力系数, 间距继续增加, 这种相互作用的强度会随之减弱, 脉动升力系数随之减小. 本文的研究能为风力俘能结构的列阵布局提供理论支持.
-
表 1 相关参数定义
Table 1 Parameter definition
Parameter Remark Parameter Remark Dm 0.01 m, average diameter of conical cylinder $\nu $ 1.48 × 10−5 m2/s, kinematic viscosity H 7Dm, the length of the conical cylinder ρ 1.225 kg/m3, density of air L center distance between two conical cylinders U, U∞ 5.772 m/s, velocity of inlet Δy height of the first level grid Re 3900, Reynolds number Cd 2Fd/(ρU2HDm), drag coefficient, Fd is the total
drag forceCp (p−p0)/(0.5ρU2), pressure coefficient, p, p0 are the static
pressure and reference pressure, respectivelyCl 2Fl/(ρU2HDm), lift coefficient, Fl is the total lift force Nfe number of free ends Cdmean time-averaged drag coefficient Clrms RMS value of lift coefficient 表 2 时间步长、圆周节点数和第一层高度对时均阻力系数, 脉动升力系数和Strouhal数的影响
Table 2 The influence of time step and grid parameters on the Cdmean , Clrms and St
Case Gird number Δy/D Nfe Δt Re H/Dm Cdmean Clrms St Case1 1 035 424 0.001 2 0.0002 3 900 7.0 0.749 0.0151 0.184 Case2 2 862 840 0.001 2 0.0002 3 900 7.0 0.756 0.0123 0.167 Case3 5 407 524 0.001 2 0.0002 3 900 7.0 0.758 0.0121 0.164 Case4 2 956 920 0.0005 2 0.0002 3 900 7.0 0.759 0.0119 0.161 Case5 2 709 560 0.002 2 0.0002 3 900 7.0 0.741 0.0112 0.154 Case6 2 862 840 0.001 2 0.00002 3 900 7.0 0.760 0.0126 0.164 Case7 2 862 840 0.001 2 0.002 3 900 7.0 0.743 0.0110 0.115 num.[25] — — 1 — 3 900 6.0 0.932 0.0190 — exp.[26] — — 2 — 88 000 5.0 0.742 — — -
[1] May R, Jackson CR, Middel H, et al. Life-cycle impacts of wind energy development on bird diversity in Norway. Environmental Impact Assessment Review, 2021, 90: 106635 doi: 10.1016/j.eiar.2021.106635
[2] 王同光, 田琳琳, 钟伟等. 风能利用中的空气动力学研究进展Ⅰ: 风力机气动特性. 空气动力学学报, 2022, 出版中 Wang Tongguang, Tian Linlin, Zhong Wei, et al. A review of aerodynamic research progress in wind energy Ⅰ: wind turbine aerodynamic characteristics. Acta Aerodynamica Sinica, 2022, in press(in Chinese)
[3] 白旭, 陈云. 串列不等直径双圆柱海流能发电振子涡激振动数值模拟. 可再生能源, 2020, 38(12): 1698-1704 (Bai Xu, Chen Yun. Numerical simulation of vortex induced vibration of a series of unequal diameter double cylindrical current energy generator. Renewable Energy Resources, 2020, 38(12): 1698-1704 (in Chinese) doi: 10.3969/j.issn.1671-5292.2020.12.021 [4] Bilgen O, Kochersberger K, Diggs E, et al. Morphing wing aerodynamic control via macro-fiber-composite actuators in an unmanned aircraft//AIAA Infotech@ Aerospace 2007 Conference and Exhibit, 2007: 2741
[5] Zhang ZH, Tu JH, Zhang K, et al. Vortex characteristics and flow-induced forces of the wavy cylinder at a subcritical Reynolds number. Ocean Engineering, 2021, 222(2): 108593
[6] Zhang H, Yang JM, Xiao LF, et al. Large-eddy simulation of the flow past both finite and infinite circular cylinders at Re = 3900. Journal of Hydrodynamics, 2015, 27(2): 195-203 doi: 10.1016/S1001-6058(15)60472-3
[7] 赵桂欣, 桂洪斌. 亚临界状态下两端自由圆柱绕流减阻数值模拟分析. 中国舰船研究, 2020, 15(S1): 12-19 (Zhao Guixin, Gui Hongbin. Numerical simulation analysis on drag reduction of circular cylinder with two free ends in subcritical state. Chinese Journal of Ship Research, 2020, 15(S1): 12-19 (in Chinese) [8] 赵萌, 毛军, 郗艳红. 高雷诺数下有限长圆柱绕流阻力特性研究. 机械工程学报, 2015, 51(22): 176-182 (Zhao Meng, Mao Jun, Xi Yanhong. Research on drag characteristic of flow around finite circular cylinder at High Reynolds Numbers. Journal of Mechanical Engineering, 2015, 51(22): 176-182 (in Chinese) doi: 10.3901/JME.2015.22.176 [9] 乔永亮. 有限长圆柱绕流的流动特性及其机理研究. [硕士论文]. 哈尔滨: 哈尔滨工业大学, 2016 Qiao Yongliang. Study on the flow characteristics and me chanism of flow around a finite-length circular cylinder. [Master Thesis], Harbin: Harbin Institute of Technology, 2016 (in Chinese))
[10] 杨耀宗. 波浪锥柱流固耦合振动机理的研究. [硕士论文]. 武汉: 武汉理工大学, 2020 Yang Yaozong. Study on the fluid-solid coupling vibration mechanism of wavy cone cylinder. [Master Thesis]. Wuhan: Wuhan University of Technology, 2020 (in Chinese))
[11] Alam MM, Sakamoto H, Zhou Y. Determination of flow configurations and fluid forces acting on two staggered circular cylinders of equal diameter in cross-flow. Journal of Fluids and Structures, 2005, 21(4): 363-394 doi: 10.1016/j.jfluidstructs.2005.07.009
[12] Alam MM, Meyer JP. Two interacting cylinders in cross flow. Physical Review E, 2011, 84: 056304
[13] Zdravkovich MM. Review of flow interference between two circular cylinder s in various arrangements. Asme Journal of Fluids Engineering, 1977, 99: 618-633 doi: 10.1115/1.3448871
[14] 杜晓庆, 顾李敏, 吴葛菲等. 错列双圆柱的脉动升力特性及其流场机理. 哈尔滨工业大学学报, 2020, 52(8): 62-71 (Du Xiaoqing, Gu Limin, Wu Gefei. Fluctuating lift and flow field mechanisms of two staggered circular cylinders. Journal of Harbin Institute of Technology, 2020, 52(8): 62-71 (in Chinese) [15] 杜晓庆, 林伟群, 施春林等. 高雷诺数下并列双圆柱绕流的大涡模拟. 哈尔滨工业大学学报, 2019, 51(6): 193-200 (Du Xiaoqing, Lin Weiqun, Shi Chunlin, et al. Large eddy simulation of flow around two side - by-side circular cylinders at a high reynolds number. Journal of Harbin Institute of Technology, 2019, 51(6): 193-200 (in Chinese) [16] 杜晓庆, 张盛华, 王玉梁等. 错列双圆柱下游圆柱的升力机理. 同济大学学报(自然科学版), 2019, 47(4): 491-498 (Du Xiaoqing, Zhang Shenghua, Wang Yuliang, et al. Lift mechanism of downstream cylinder on two staggered circular cylinders. Journal of Tongji University (Natural Science) , 2019, 47(4): 491-498 (in Chinese) [17] Papaioannou GV, Yue D, Triantafyllou MS, et al. Three-dimensionality effects in flow around two tandem cylinders. Journal of Fluid Mechanics, 2006, 558: 387-413 doi: 10.1017/S0022112006000139
[18] 刘小兵, 姜会民, 王世博等. 串列三圆柱的脉动气动力特性试验研究. 振动与冲击, 2021, 40(13): 96-103 (Liu Xiaobing, Jiang Huimin, Wang Shibo, et al. Tests for fluctuating aerodynamic force characteristics of three cylinders in series. Journal of Vibration and Shock, 2021, 40(13): 96-103 (in Chinese) [19] 涂佳黄, 曹波, 谭潇玲等. 串列双圆柱体绕流特性与互扰效应研究. 应用力学学报, 2019, 36(4): 869-875, 998 (Tu Jiahuang, Cao Bo, Tan Xiaolin, et al. Study on flow characteristics and mutual interference effect around tandem double cylinders. Chinese Journal of Applied Mechanics, 2019, 36(4): 869-875, 998 (in Chinese) [20] 唐涛, 朱红钧, 周同明等. 串列梯形柱绕流水动力特性数值分析. 水动力学研究与进展, 2021, 36(4): 471-480 (Tang Tao, Zhu Hongjun, Zhou Tongming, et al. Numerical investigation of flow past two tandem trapezoidal cylinders and their hydrodynamic characteristics. Chinese Journal of Hydrodynamics, 2021, 36(4): 471-480 (in Chinese) [21] 陈林烽. 基于Navier-Stokes方程残差的隐式大涡模拟有限元模型. 力学学报, 2020, 52(5): 1314-1322 (Chen Linfeng. A residual-based unresolved-scale finite element modelling for implict large eddy simulation. Chinese Journal of Theoretical and Applied Mechanics, 2020, 52(5): 1314-1322 (in Chinese) [22] 吴霆, 时北极, 王士召等. 大涡模拟的壁模型及其应用. 力学学报, 2018, 50(3): 453-466 (Wu ting, Shi Beiji, Wang Shizhao, et al. Wall-model for large-eddy simulation and its applications. Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(3): 453-466 (in Chinese) [23] Smagorinsky J. General circulation experiments with the primitive equations. Monthly Weather Review, 1963, 91: 99-165 doi: 10.1175/1520-0493(1963)091<0099:GCEWTP>2.3.CO;2
[24] Ma HT, He L, et al. Large eddy simulation of natural convection heat transfer and fluid flow around a horizontal cylinder. International Journal of Thermal Sciences, 2021, 162: 106789 doi: 10.1016/j.ijthermalsci.2020.106789
[25] 赵桂欣, 桂洪斌, 王晓聪. 有限长波浪形圆柱绕流数值模拟. 哈尔滨工业大学学报, 2021, 53(6): 163-170 (Zhao Guixin, Gui Hongbin, Wang Xiaocong. Numerical simulation of flow around finite ength wavy cylinders. Journal of Harbin Institute of Technology, 2021, 53(6): 163-170 (in Chinese) [26] Zdravkovich M, Brand V, Mathew G, et al. Flow past short circular cylinders with two free ends. Journal of Fluid Mechanics, 2006, 203(203): 557-575
[27] Essel EE, Tachie MF, Balachandar R. Time-resolved wake dynamics of finite wall-mounted circular cylinders submerged in a turbulent boundary layer. Journal of Fluid Mechanics, 2021, 917: A8
[28] 梁圣召, 周军伟, 梅蕾. 不同工况下有限长圆柱绕流实验测量及分析. 舰船科学技术, 2021, 43(11): 19-24 (Liang Shengzhao, Zhou Junwei, Mei Lei. Experiments and numerical simulations of the three-dimensional flow around a finite-length circular cylinder. Ship Science and Technology, 2021, 43(11): 19-24 (in Chinese) [29] Gao W, Nelias D, Lű Y, et al. Numerical investigations on drag coefficient of circular cylinder with two free ends in roller bearings. Tribology International, 2018, 123: 43-49 doi: 10.1016/j.triboint.2018.02.044
[30] Frederich O, Scouten J, Luchtenburg DM, et al. Large-scale dynamics in the flow around a finite cylinder with a ground plate. Fluid Dynamics Research, 2011, 43(1): 015504 doi: 10.1088/0169-5983/43/1/015504
[31] 王义乾, 桂南. 第三代涡识别方法及其应用综述. 水动力学研究与进展(A辑), 2019, 34(4): 413-429 (Wang Yiian, Gui Nan. A review of the third-generation vortex identification method and its applications. Chinese Journal of Hydrodynamics, 2019, 34(4): 413-429 (in Chinese) -
期刊类型引用(3)
1. 邹琳,吴伟男,王家辉,郑云龙,徐汉斌. 三角形排列波浪锥柱绕流流动特性数值分析. 东南大学学报(自然科学版). 2024(01): 224-231 . 百度学术
2. 邓斌,尹龙斌,黄姣凤,熊凯,蒋昌波. 波浪与新型双排开孔圆筒防波堤相互作用三维数值模拟. 力学学报. 2023(04): 845-857 . 本站查看
3. 朱健申,王小静,王宗龙,刘腾岩. 串列水翼绕流非定常特征实验研究. 实验室研究与探索. 2023(08): 30-34 . 百度学术
其他类型引用(1)