在实验研究方面日本横须贺国防研究院Igarashi教授在风洞中对开缝圆柱流场进行了系统的测量，总结了开缝圆柱在不同工况下的流场特征[2-3]；纽卡斯尔大学的Olsen教授和Rajagopalan则重点研究了缝隙大小和圆柱外轮廓形状对斯特劳哈尔数和阻力系数的影响[4]；台湾成功大学Miau教授等[5]在Igarashi[2-3]的研究基础之上，在水洞和风洞里对不同缝隙比的开缝圆柱流场进行了更为细致测量工作。国内北航的王晋军教授等的实验研究也表明开缝圆柱缝隙倾角在一定范围内可以减小流动阻力[6]；中国计量学院的王慧和黄咏梅教授通过实验的方法对带狭缝圆柱钝体流量测量仪器进行了研究，表明带狭缝圆柱的涡街信号强度和信噪比高，斯特劳哈尔数线性关系较好[7]；厦门大学的鲍锋教授等在循环水洞中对开缝圆柱缝隙内部流场进行了细致的测量工作，对开缝圆柱的流态进行了补充[1]。
在数值计算方面，针对圆柱扰流学术界普遍认为采用大涡模拟(large eddy simulation)的方式进行三维数值模拟具有较高的准确性[8-11]，但直接对开缝圆柱绕流特性直接数值计算研究的相对较少[1]。国内厦门大学的鲍锋教授等采用大涡模拟的方式对不同雷诺数下不同缝隙宽的开缝圆柱仿真发现开缝圆柱尾流场具有一定的三维性，其仿真结果与实验数据能较好地吻合[1]；中国计量学院的王慧和黄咏梅也用数值模拟方法发现了带狭缝圆柱产生的涡街信号强度高，旋涡尺寸较大，斯特劳哈尔数线性度好[7]；但值得注意的是也存在一些****从其他领域侧面对开缝圆柱流场特性进行了分析，如对开缝旋转圆筒固定小翼构成的飞机颤振激励器的气动特性[12]和开缝圆柱强化对流换热进[13-15]的仿真结果都从侧面对开缝圆柱绕流流场特性进行了分析。
本文主要工作是基于厦门大学流体&粒子成像测速(PIV)系统实验室的精密循环水洞，采用非介入式流场测试技术，对不同狭缝比的开缝圆柱在不同缝隙倾斜角下的流场特性进行了细致的测量，本文的工作进一步对开缝圆柱绕流流场特性进行了补充。
1 实验模型与研究方法 本文所涉及的实验测量工作均在开放式循环水槽中进行，该型多功能精密循环卧式水槽结构示意如图 1所示，水槽试验段上方设置有光学平板和六分量测力天平，天平下端开设有直径10 mm的模型安装孔，可垂直悬挂模型并实现攻角调节；水槽动力配套有调节转速的数字变频器，频率调节范围为0~70 Hz，对应实验段水位高度为500 mm的水流速度为0~0.3 m/s，水流速度连续可调，在工作范围内湍流度被控制1%以内；实验段水体截面尺寸为500 mm×500 mm，长度3 000 mm，采用高透光率的光学玻璃搭建而成，实现四面良好的光学可达性。
 |
| 图 1 实验水槽 Fig. 1 Experiment water channel |
| 图选项 |
本文需对开缝圆柱绕流的瞬态流场进行分析，因此建立了时序PIV系统，系统参数如表 1所示。该PIV系统包括主要包括了高分辨率相机、激光器、同步器、PIV图像处理分析软件和示踪粒子等。其中PIV相机选用高速CMOS相机，分辨率200万像素，最大触发频率为1 kHz。本文实验图像采集和数据后处理软件采用Dynamic Studio软件，系统误差在1%以内。
表 1 PIV系统参数 Table 1 PIV system parameters
| 组件 | 型号 | 参数及数值 |
| 激光 | 最大输出功率15 W，波长532 nm | |
| 相机 | SpeedSense 9040 | 1 632像素×1 200像素 |
| 镜头 | Micro-Nikkor | 焦距55 mm |
| 处理软件 | Dynamic Studio V3.40 | |
| 示踪粒子 | PSP粒子直径20 μm |
表选项
为保证PIV实验测量的精度，模型外径公差控制在±0.2 μm，并进行黑色阳极氧化处理，减少激光的发射，实验模型如图 2所示。
 |
| 图 2 实验模型 Fig. 2 Experimental model |
| 图选项 |
实验时，模型垂直吊装在拖曳台车上方的天平上，通过控制天平内部伺服电机旋转角度达到改变开缝圆柱与来流方向夹角的目的。测量时激光切面从水槽侧面玻璃垂直入射，切面与水流流向平行，并照射在模型中部。实验模型安装与激光切面位置如图 3所示，实验堵塞比低于8%。
 |
| 图 3 实验模型及激光切面安装位置 Fig. 3 Mounting position of experimental model andlaser section |
| 图选项 |
实验模型选用直径d=40 mm，缝隙宽度s=4，6 mm轴向有效长度l=440 mm的开缝圆柱，β为开缝圆柱缝隙中心线与来流方向的夹角，具体参数见表 2。
表 2 开缝圆柱尺寸参数 Table 2 Geometrical parameters of slotted cylinder
| 参数 | 数值 |
| d/mm | 40 |
| l/mm | 440 |
| d/l | 11 |
| s/mm | 4，6 |
| s/d | 0.10，0.15 |
| α/(°) | 8.63 |
| 注：d/l为长宽比；α＝arcsin(s/d)为开缝角度。 | |
表选项
2 开缝圆柱流场特性的流动显示 首先进行了开缝圆柱近场区流场特性流动显示实验。为了更加细腻地显示流场信息，选用激光粒子，采用波长为532 nm的绿色激光片光照明，摄像机拍摄图像。图 4为来流速度u=0.15 m/s，β=90°(基于圆柱直径的雷诺数Red=5 837，下同)工况下开缝圆柱缝隙内部及近壁面的流体流动一个完整周期T的演化过程。
 |
| 图 4 Red=5 837时缝隙内部流动 Fig. 4 Fluid flow in slit-vent at Red=5 837 |
| 图选项 |
t0时刻，开缝圆柱上侧边界层被“吸进”缝隙，缝隙下方的流体被“吹出”，此时在圆柱下表面附着逆时针旋转的小旋涡。
t1时刻，缝隙上端生成2个旋转方向相反的协动小旋涡，t0时刻在圆柱下表面形成的旋涡涡核发展变大，逐渐从圆柱表面开始脱落。
t2时刻，缝隙内的流体短暂静止，该时刻缝隙内部无旋涡存在。
t3时刻，缝隙内部经历短暂静止后，缝隙上下两端流向进入逆转，缝隙下端口吸入边界层，并在缝隙内部形成2个旋向相反的协动小漩涡，而缝隙下端口吹出边界层。
总结边界层吹吸与旋涡演化过程之间的演化关系可概述为当缝隙下端口发生边界层吸入效应，则缝隙入口处依上而下形成顺时针和逆时针方向旋涡，同时在圆柱的上表面后方生成顺时针旋涡；当下方由边界层吸入效应转变为边界层吹出效应时，圆柱后方旋涡逐渐发展至脱落；而当缝隙下方处于边界层吹出的状态时，流体沿缝隙右表面底部向外“溢出”，同时圆柱下表面后方生成逆时针旋涡，当下方由边界层吹出向边界层吸入过渡时，旋涡逐渐发展至脱落。在整个周期内缝隙内部形成近似于“活塞”的流动现象，即流体上下周期性地振荡，并会形成一系列的协动小漩涡。
3 PIV实验结果及分析 3.1 不同来流速度对其脱落涡的影响 在PIV实验过程中，激光切面从水槽侧面玻璃垂直入射，切面与水流流向平行，相机在水槽底部拍摄，记录开缝圆柱尾缘的流场信息，根据实验设备的性能和流场测量需求，在实验过程中设置高频相机的工作频率为100 Hz，在一定的相机的存储空间下尽可能捕捉较长时域的流场状态，相机捕捉的图幅大小和相应的坐标说明如图 5所示。
 |
| 图 5 PIV数据坐标系 Fig. 5 Coordinate system of PIV data |
| 图选项 |
本文对开缝圆柱脱落频率f的统计采用环量监控的方式，在原始采集图片划定区域作为脱落涡周期监控窗口(图 5红色框内)，对于速度场解算，采用自适应相关(adaptive correlation)算法，查问区域(interrogation areas)32像素×32像素，水平和垂直方向的重叠度(overlap)为50%，窗口监控指标取其顺时针旋涡环量Γ变化，采用这种方式可以在一定程度上避免由于PIV分辨率问题带来的相关统计信息的误差。环量是流场中速度在某一封闭周线切线上的分量V→沿封闭曲线S→的积分，环量Γ的计算可以根据斯托克斯公式表述为式(1)，即可以根据对涡量的面积分A求得环量，二维情况下，涡量ω的计算的如式(2)所示，u与v分别为x与y方向的速度分量，L为封闭曲线。
 | (1) |
 | (2) |
下面以狭缝比s/d=0.15，β=90°的开缝圆柱周期统计进行分析，实验来流速度u设置在0.035、0.070、0.100、0.145和0.180 m/s时，截取监控窗口环量Γ在12 s内的变化，如图 6所示。
 |
| 图 6 β=90°时不同Red下环量发展变化曲线 Fig. 6 Circulation development curves at differentRed when β=90° |
| 图选项 |
在本文设置的雷诺数范围内，随着雷诺数的提高，监控区内脱落涡整体环量增强，每个周期内的环量Γ变化存在差异。
在环量统计的基础上运用FFT变换，得到不同雷诺数下环量变化频谱图，见图 7。图中：P为能量。
 |
| 图 7 β=90°时不同Red下的频谱图 Fig. 7 Spectrogram at different Red when β=90° |
| 图选项 |
从不同来流速度下的频谱图可以看出，随着流场速度的提高，开缝圆柱的脱落涡脱落频率增大，特征频率信号也逐渐增强。
斯特劳哈尔数St是当地惯性力与迁移力的比值，反映了流动非定常运动的相似性，对于周期性的非定常运动就反映周期性相似，所以和周期性相关的非定常流动由St决定[1]，在频率统计f的基础上，这里采用式(3)计算斯特劳哈尔数。
 | (3) |
不同u下的Red、 f和St如表 3所示，通过统计数据，可以发现本文测量计算的结果与文献[5]结论吻合，开缝圆柱的斯特劳哈尔数St具有很好的线性关系，其实验值在0.22~0.23之间。
表 3 不同u下的Red、 f和St Table 3 Red，f and St at different u
| u/(m·s-1) | 0.035 | 0.070 | 0.100 | 0.145 | 0.180 |
| Red | 1 393 | 2 786 | 3 980 | 5 771 | 7 164 |
| f/Hz | 0.196 | 0.393 | 0.566 | 0.822 | 1.020 |
| St | 0.224 | 0.225 | 0.226 | 0.227 | 0.227 |
表选项
本节采用了本征正交分解(POD)方法对开缝圆柱流场做做更为细致的分析。POD方法是基于主元分析思想，依据Karhumen-Loève展开，将流场进行分解，并按照能量比计算出各阶模态，识别出流动中大尺度的拟序结构[16]。本节中采用的POD方法是Dynamic Studio软件自带的快照POD(Snapshot POD)算法。该算法将瞬时流场分解为时均速度场和脉动速度场，然后对脉动速度场进行POD分解，这样流场的全部信息可分解为时均流场和脉动流场，任意瞬时流场可视为各阶模态一定的线形组合。在实际运用中POD方法对于流场大尺度结构解算具有明显的优势[16-17]。
本节取能量占比90%以上(0~5阶)的模态，对来流速度u=0.10 m/s(Red=3 980)下基准圆柱、狭缝比s/d=0.10和s/d=0.15的开缝圆柱流场进行流场重构，得到一个周期内的速度矢量图(取图 5蓝色实线框区域)，如图 8~图 10所示。
 |
| 图 8 基准圆柱一个周期内的流态 Fig. 8 Fluid flow of a base cylinder in a period |
| 图选项 |
 |
| 图 9 s/d=0.10的开缝圆柱一个周期内的流态 Fig. 9 Fluid flow of a slotted cylinder with s/d=0.10 in a period |
| 图选项 |
 |
| 图 10 s/d=0.15的开缝圆柱一个周期内的流态 Fig. 10 Fluid flow of a slotted cylinder with s/d=0.15 in a period |
| 图选项 |
通过对比以上不同实验模型在一个周期内的流态结构，在相同雷诺数下，缝隙的存在对圆柱尾流结构产生明显的影响，其卡门涡街特征更为明显，漩涡脱落周期缩短；狭缝比s/d=0.10和s/d=0.15的开缝圆柱在脱落涡周期的差异不大。
3.2 缝隙倾斜角对脱落涡的影响 为进一步探究开缝圆柱的流场特性，本节在u=0.10 m/s(Red =3 980)下通过改变开缝圆柱缝隙与来流方向的夹角β，探究开缝圆柱St对缝隙倾斜角β的敏感性，如图 11所示。
 |
| 图 11 开缝圆柱示意图 Fig. 11 Sketch of circular cylinder with a slit |
| 图选项 |
根据文献[2-3]，开缝圆柱缝隙倾斜角在60°≤β≤90°，开缝圆柱缝隙对边界层的抽吸现象主导尾流流动，故本节设置了从70°~90°共计21组不同倾斜角β工况，统计得到在不同倾斜角下，开缝圆柱的斯特劳哈尔数St的变化见图 12，可以看出随着缝隙倾斜角β的变大，其St逐渐变大，相比狭缝比s/d=0.15和s/d=0.10的斯特劳哈尔数，s/d=0.15的开缝圆柱St稳定性更强。
 |
| 图 12 St变化曲线在β=70°~90° Fig. 12 St development curves at β=70°-90° |
| 图选项 |
本节还选取了狭缝比s/d=0.15开封圆柱在β为90°、85°和80°条件下一个周期内的速度云图变化，如图 13所示。
 |
| 图 13 β=90°，85°，80°一个周期内的流态 Fig. 13 Fluid flow of a slotted cylinder in a period at β=90°,85°，80° |
| 图选项 |
对比可以发现，在实验参数范围内，随着缝隙倾斜β的减小，整体流态变化不大，但其脱落涡尾迹会出现轻微的上扬，对缝隙倾斜角存在一定的敏感性。
4 结 论 本文采用流动显示技术，观察了开缝圆柱的近场区流场特性，探究缝隙对其流场的影响，并采用PIV测量技术测量了开缝圆柱的缝隙倾斜角对其脱落涡的影响，得到以下结论：
1) 缝隙垂直于来流的开缝圆柱在实验雷诺数范围内，当流体流经开缝圆柱时，在圆柱缝隙内部形成了类似于“活塞”的流动现象，并导致了近壁面区边界层出现的“吹吸”效应，缝隙内部会形成一定数量的协动小漩涡。
2) 当缝隙倾斜角β=90°时，实验雷诺数范围内，随着雷诺数的增大，其脱落涡频率也明显增大，涡街特征更加显著，开缝圆柱的St与雷诺数具有很强的线性关系。
3) 开缝圆柱缝隙在不同倾斜角下，随着缝隙倾斜角β的变大，其脱落涡的频率对开缝圆柱缝隙倾斜角存在一定的敏感性，且缝隙倾斜角对尾流形态也产生了一定的影响，相比狭缝比s/d=0.15与s/d=0.10的开缝圆柱，s/d=0.15开缝圆柱的St稳定性更佳。
参考文献
| [1] | 杨锦文, 何意, 鲍锋. 开缝圆柱近场尾流的大涡模拟和实验测量[J].空气动力学报, 2014, 32(3): 308–314.YANG J W, HE Y, BAO F. Large eddy stimulation and experimental measurement of near wake structures of a slotted circular cylinder[J].Acta Aerodynamica Sinica, 2014, 32(3): 308–314.(in Chinese) |
| [2] | IGARASHI T. Flow characteristics around a circular cylinder with a slit:1st report, flow control and flow patterns[J].Bulletin of the JSME, 1978, 21(154): 656–664.DOI:10.1299/jsme1958.21.656 |
| [3] | IGARASHI T. Flow characteristics around a circular cylinder with a slit:2nd report, effect of boundary layer flow[J].Bulletin of the JSME, 1982, 25(207): 1389–1397.DOI:10.1299/jsme1958.25.1389 |
| [4] | OLSEN J F, RAJAGOPALAN S. Vortex shedding behind modified circular cylinders[J].Journal of Wind Engineering and Industrial Aerodynamics, 2000, 86(1): 55–63.DOI:10.1016/S0167-6105(00)00003-9 |
| [5] | PENG B H, MIAU J J, BAO F, et al. Performance of vortex shedding from a circular cylinder with a slit normal to the steam[J].Flow Measurement and Instumentation, 2012, 25: 54–62.DOI:10.1016/j.flowmeasinst.2011.07.003 |
| [6] | 谭广琨, 王晋军, 李秋胜. 圆柱体减阻技术及其机理初步研究[J].北京航空航天大学学报, 2001, 27(6): 658–661.TAN G K, WANG J J, LI Q S. Drag reduction technique of cylinder and mechanism research[J].Journal of Beijing University of Aeronautics and Astronautics, 2001, 27(6): 658–661.(in Chinese) |
| [7] | 王慧, 黄咏梅. 带狭缝圆柱钝体流量测量特性研究[J].传感技术学报, 2013, 26(12): 1661–1665.WANG H, HUANG Y M. Study of flowrate measurement characteristics of a circular cylinder with a slit[J].Chinese Journal of Sensors and Actuators, 2013, 26(12): 1661–1665.(in Chinese) |
| [8] | 郝鹏, 李国栋, 杨兰, 等. 圆柱绕流流场结构的大涡模拟研究[J].应用力学学报, 2012, 29(4): 437–443.HAO P, LI G D, YANG L, et al. Large eddy simulation of the circular cylinder flow in different regimes[J].Chinese Journal of Applied Mechanics, 2012, 29(4): 437–443.(in Chinese) |
| [9] | BREUER M. A challenging test case for large eddy simulation:High Reynolds number circular cylinder flow[J].International Journal of Heat and Fluid Flow, 2000, 21(5): 648–654.DOI:10.1016/S0142-727X(00)00056-4 |
| [10] | SARKAR S, SARKAR S. Large-eddy simulation of wake and boundary layer interactions behind a circular cylinder[J].Journal of Fluid Engineering, 2009, 131(9): 091201–14.DOI:10.1115/1.3176982 |
| [11] | SARKAR S, SARKAR S. Vortex dynamics of a cylinder wake in proximity to a wall[J].Journal of Fluids and Structures, 2010, 26(1): 19–40.DOI:10.1016/j.jfluidstructs.2009.08.003 |
| [12] | 杨智春, 党会学, 孙勇军, 等. 带开缝旋转圆筒固定小翼气动特性的数值模拟[J].工程力学, 2009, 29(11): 222–227.YANG Z C, DANG H X, SUN Y J, et al. Numerical simulation of the aerodynamic characteristics of a fixed wing with trailing edge rotating slotted cylinder[J].Engineering Mechanics, 2009, 29(11): 222–227.(in Chinese) |
| [13] | 申春赟, 杨茉, 王津, 等. 圆内开缝圆不同开缝方向自然对流换热[J].上海理工大学学报, 2013, 35(5): 425–429.SHEN C Y, YANG M, WANG J, et al. Natural convection heat transfer characteristics of a circle with an internal concentric round slotted in different direction[J].Journal of University of Shanghai for Science and Technology, 2013, 35(5): 425–429.(in Chinese) |
| [14] | 汪健生, 刘志毅, 张金凤, 等. 斜截椭圆柱式涡流发生器强化传热的大涡模拟[J].机械工程学报, 2007, 43(10): 55–61.WANG J S, LIU Z Y, ZHANG J F, et al. Large eddy simulation on heat transfer enhancement of inclined-cut ellipsoidal vortex generator[J].Journal of Mechanical Engineering, 2007, 43(10): 55–61.DOI:10.3901/JME.2007.10.055(in Chinese) |
| [15] | 汪健生, 王晓, 朱强, 等. 湍流边界层内钝体扰流的流动与传热特性[J].机械工程学报, 2015, 51(24): 168–176.WANG J S, WANG X, ZHU Q, et al. Flow and heat transfer characteristics of flow past the bluff body in turbulent boundary layer[J].Journal of Mechanical Engineering, 2015, 51(24): 168–176.DOI:10.3901/JME.2015.24.168(in Chinese) |
| [16] | MA L Q, FENG L H, PAN C, et al. Fourier mode decomposition of PIV data[J].Science China Technological Sciences, 2015, 58(11): 1935–1948.DOI:10.1007/s11431-015-5908-y |
| [17] | 王汉封, 徐胜金. 用POD方法研究有限长方形棱柱尾流的双稳态现象[J].空气动力学学报, 2014, 32(6): 827–833.WANG H F, XU S J. On the bistable phenomenon of a finite length square cylinder wake with POD method[J].Acta Aerodynamica Sinica, 2014, 32(6): 827–833.(in Chinese) |
