APP下载

展向周期振动平板的近壁面流场特性实验

2017-10-17邓飞赵建鹏董海军单明明唐勇

哈尔滨工程大学学报 2017年9期
关键词:摩擦阻力边界层湍流

邓飞,赵建鹏,董海军,单明明,唐勇

(1.西北工业大学 航海学院, 陕西 西安 710072; 2.西北工业大学 机电学院, 陕西 西安 710072)

展向周期振动平板的近壁面流场特性实验

邓飞1,赵建鹏1,董海军2,单明明2,唐勇1

(1.西北工业大学 航海学院, 陕西 西安 710072; 2.西北工业大学 机电学院, 陕西 西安 710072)

为了研究壁面展向振动对壁面摩擦阻力及近壁面流场湍流强度的影响,本文开展了在水槽中利用粒子图像测速方法(PIV)分析薄板表面流场信息的实验研究。将一薄板放置于低速水槽中,通过激振器使薄板展向振动,调节激振器的振幅和频率实现不同工况,薄板前缘安置激流丝用以激发壁面湍流。利用粒子图像测速方法连续采集一段时间内的近壁面法向速度剖面,观察近壁面漩涡分布情况,并利用将平均速度剖面与经典湍流对数律拟合的方法估测壁面摩擦应力,为了研究湍流强度影响,对结果进行处理,得到流向速度与法向速度脉动均方根。结果发现:振动后近壁面的漩涡分布明显较少,速度波动明显减弱,流场速度特性趋于层流特性;在振动工况下估测的壁面摩擦阻力均小于非振动工况下,并且在实验工况内振幅和频率增大,减阻效果增强,说明了振动的强度对减阻特性有影响;同时,振动均减小了脉动速度均方根,削弱了近壁面湍流强度,也说明了减阻的机理是降低了近壁面湍流强度,使近壁面流动趋于层流。

展向周期振动; 粒子图像测速; 漩涡分布; 平均速度剖面; 摩擦速度; 摩擦阻力; 近壁面湍流强度

Abstract:In this study, we conducted a particle image velocimeter (PIV) experiment in a water channel to investigate the effects of the spanwise oscillation of a wall on drag reduction and turbulent flow. A plate is oscillated by a vibrator whose amplitude and frequence is adjustable. We recorded turbulent vortex information in a direction normal to the plate, determined the mean velocity profile of the boundary layer, and then fitted the mean velocity profile with the classic logarithmic law to estimate the wall friction drag. To determine the influence of turbulence intensity, we obtained the root- mean- square value of the fluctuation velocity in the flow and normal directions. The results show that the vortices distribution near the wall after vibration becomes less obvious, the speed oscillation decreases sharply, the flow speed tends to be laminar; and the wall friction drag is smaller than that of the non- vibration operation state. In addition, we found that increasing the amplitude and frequency in the experiment led to a better drag reduction effect, which indicates that the vibration intensity influences the drag reduction. Also, we found vibration to reduce the root- mean- velocity value of the fluctuation velocity and to reduce the turbulence intensity near the wall, which illustrates that the drag reduction mechanism lowers the turbulence intensity near the wall and causes the flow near the wall to be laminar.

Keywords:spanwise oscillations of wall; particle image velocity; vortices distribution; mean velocity profile; friction velocity; friction drag; near- wall turbulence intensity

减小水下运动物体的阻力,提高航速、节约能源消耗一直是人类追求的目标。与空气相比,水介质的物理特性(密度、粘性、压缩性、空化特性等)具有巨大差异,因此其减阻技术研究更加复杂。迄今为止,在水下减阻方面的研究成果较之航空与航天领域仍逊色不少,但也日益成为人们关注的热点。

近年来学者提出了许多减小壁面阻力的方法,如改变壁面结构减阻技术、随行波减阻技术、超空泡减阻技术以及壁面展向振动减阻技术等。Jung W J等通过对槽道湍流的直接数值模拟发现壁面展向周期振动能够有效地减小壁面阻力[1]。Laadhari F等则通过热线测速实验研究了壁面展向振动时湍流边界层的流动情况,实验结果同样表明,靠近振动壁面处湍流边界层的平均速度梯度减小,湍流强度降低,从而验证了,壁面展向周期振动可以减小湍流边界层的表面的摩擦力[2]。此外,文献[3-4]通过在风洞中对振动平板表面流场使用热线测速,也发现了湍流对数区上移、湍流强度减弱的规律;发现斯托克斯层依然保持层流,就像由于振动导致湍流边界层消失;振动时的上扫事件信号的时间和强度分别削弱了78%和64%,文献[5-6]也做出了相关研究。

国内学者在振动减阻方法的数值模拟方面也做了大量的研究,黄乐萍等利用Fourier- Chebyshev 谱方法,通过直接数值模拟,对槽道湍流壁面展向周期振动抑制壁湍流、实现减阻的内在机理进行了研究,认为振动壁面的减阻主要有两个机理,一是振动引起的涡与条带的倾斜,这导致负展向涡的产生;二是振动引起的涡与条带间的滑移, 这导致条带的加宽以及流向涡的减弱[7-9]。同时,存在一个最优的振动波数,在其调制下,固有流场对诱导流场的影响最弱,而诱导流场对固有流场的影响显著,减阻效果最好。

文献[10-12]分别采用传统Smagorinsky模式、动力Smagorinsky模式和Cui模式基于Kolmogorov方程对壁面在展向作周期振动的槽道湍流进行了大涡模拟,通过对雷诺应力输运方程的分析研究了壁面展向周期振动的减阻机理,并从DNS结果的统计平均角度入手,指出Reynolds应力输运方程中压力变形项的变化在展向壁面振动减阻过程中起到了关键作用。

董海军等应用分子动力学方法研究壁面展向振动条件下流体的流动特性[13]。通过对固体平板壁面施加不同振动参数的展向振动,结果表明对壁面施加展向振动可以有效地减小固液界面间的摩擦力,减阻的幅度与振动参数关系很大。另外,文献[14-19]在其他相关方面也进行了研究。

本文为了实验观测振动减阻现象及机理,实验研究振动对湍流摩擦阻力的影响,探究减阻的机理,开展了在低速水槽中使用流场粒子图像分析方法研究振动平板流场的实验。

1 实验装置和过程

实验原理如图1所示,实验平台为低速层流水槽,平板通过激振器驱动进行展向振动,激光器产生面激光,激光平面为迎来流方向的铅垂面。如图2所示,通过高速摄像机采集平板表面流场的速度信息,通过处理平均速度剖面,分析壁面摩擦阻力所受的影响;通过不同高度速度脉动均方根曲线,比较振动对流场湍流强度的影响;通过对平板法向流场平面速度矢量信息观察,发现振动对流场造成的影响。

图1 粒子图像测速方法Fig.1 Experimental principle diagram

图2 平板放置位置Fig.2 Plate position

激振器在频率0.1~100 Hz、振幅0~5 mm范围内可调,安装在可调节高度的云台上。振动平板尺寸160×80×1.5 mm,水槽工作段截面边长100 mm,示踪粒子采集平面为20×20 mm,水流速度0.1 m/s,基于流速与截面边长的雷诺数Re=10 000。振动平板前端粘贴了一根激流丝,用于激发湍流,流场法向湍流场脉动情况如图3所示。高速摄像机设置为曝光时间2 000 ns、帧频率300 Hz。

由于光源照射在介质表面时会产生反射,反射光强度较高,会干扰平板表面附近处的示踪粒子采集,因此,将平板表面喷一层黑漆,来削弱壁面处反光的干扰。由于平板为薄板,振动频率较低,平板展向振动不会对平板表面的流场产生太大的影响。

平板展向正弦振动方式为

y=Asin(2πft),

(1)

式中:A为激振器振动的振幅,f为激振器的频率。本实验根据振幅和频率的不同,共有6个振动工况,分别为

1)A=0;

2)A=4 mm,f=1.2 Hz;

3)A=4 mm,f=1.8 Hz;

4)A=4 mm,f=2.5 Hz;

5)A=2 mm,f=1.2 Hz;

6)A=6 mm,f=1.2 Hz。

2 摩擦阻力及湍流强度结果分析

2.1 平板近壁面涡量分布

定义流向、法向、展向为x、y、z,流向、法向、展向速度为u、v、w。利用PIV程序处理连续拍摄到的粒子分布图像,得到如图3所示的速度矢量图。图3(a)为非振动工况下的平板法向瞬时速度矢量,水平方向为流向,竖直方向为法向,灰度代表展向涡量,处理得到的边界层厚度约为10 mm,如图3(a)所示,边界层底部速度很小,随着高度的上升,速度逐渐增大,当达到边界层上边界时,与外部势流速度趋于一致;边界层底部的漩涡包含上升段和下降段,上升段中低速流体向上层冲击,下降段中边界层上层的流体向下掠扫,下降段是造成湍流摩擦阻力增大的主要原因。

振动工况A=4 mm,f=2.5 Hz下的平板法向瞬时速度矢量如图3(b)所示,通过对比发现,振动减少了边界层底部漩涡的数量,漩涡位置相对上移,底部流动更加平缓,较好地改善了底部的湍流流场。

图3 平板法向速度湍动信息Fig.3 Fluctuation velocity information in plate normal direction

2.2 不同频率和振幅对壁面摩擦阻力的影响规律

湍流中边界层的对数区速度分布规律符合y+=2.5lnu++5的对数规律,本实验中对连续采集的间隔0.008 s的1 000张速度矢量图做平均处理,得到平均速度剖面,将其中对数区部分与上述对数律拟合,无量纲过程如下

(2)

(3)

(4)

式中:u+为无量纲速度,y+为无量纲高度,u*为摩擦速度,先拟合u*,再根据u*求出壁面摩擦应力τω。振动工况A=4 mm,f=1.8 Hz时拟合的平均速度剖面如图4所示。

图4 f=1.8 Hz振动下的平均速度剖面图Fig.4 The mean velocity profile in f=1.8 Hz oscillation case

图中虚线为粘性底层和对数区的理论速度剖面曲线。图中可以看出,粒子图像测速方法处理得到的曲线与粘性底层衔接良好,与对数规律充分符合,说明测量结果基本可靠。

定义减阻率为D:

(5)

式中:u*为振动情况下的摩擦速度,u0为非振动平板的拟合摩擦速度,拟合出的u0为0.006 506 m/s。不同频率和振幅下的摩擦速度以及减阻率情况分别见表1、2。

表1不同频率下的摩擦速度以及减阻率(A=4mm)

Table1Frictionvelocityanddragreductionindifferentfrequency(A=4mm)

f/Hz1.21.82.5u*/(m·s-1)0.006190.006040.005805D/%9.513.820.4

表2不同振幅下的摩擦速度以及减阻率(f=1.2Hz)

Table2Frictionvelocityanddragreductionindifferentamplitude(f=1.2Hz)

A/mm246u*/(m·s-1)0.006230.006190.006002D/%8.39.514.9

表中可以看出,在各个振动情况内阻力均有降低,并且减阻幅度最大达到20%,说明了振动起到了一定的减阻效果,也验证了Jung W J关于振动降低摩擦阻力的观点[1]。相同振幅下,随着频率增大,减阻效果越好,说明在一定范围内,振动减阻效果随频率呈正相关。相同频率下,随着振幅增大,减阻效果也越好,说明振动减阻效果随振幅也呈正相关。以上说明了在本实验范围内,增加振动的频率和振幅,能够增大振动减阻的效果,同时有利于进一步理解湍流形成及减阻的机理。

2.3 振动对壁面湍流强度的影响

图5、6为不同频率与振幅下的流向脉动速度均方根分布,其中横坐标为测量点距壁面高度,纵坐标为流向脉动均方根值。

图5 不同频率的流向脉动速度均方根Fig.5 RMS value of u’ in different frequency

图6 不同振幅的流向脉动速度均方根Fig.6 RMS value of u’ in different amplitude

图5可以看出,流向脉动速度均方根随着频率的增大而降低;图6可以看出,流向脉动速度均方根随着振幅的增大而降低,湍流强度可降低60%之多。说明了振动使近壁面湍流强度降低,并且湍流强度降低的程度与摩擦阻力降低的程度正相关。

以上结论说明了振动通过减弱近壁面上抛和下扫等湍流活动对壁面的影响来降低摩擦阻力,验证了文献[3-4]提出的振动减阻的机理是由于振动减弱近壁面湍流强度的结论。

3 结论

1)振动使摩擦阻力降低,并且在一定范围内,减阻效果随着频率和振幅的增大而增大,预示着通过改变振动特性可以达到更理想的减阻效果。

2)通过对流向脉动速度均方根的对比,证明了振动能减弱近壁区的湍流强度,说明了振动减阻的机理是振动削弱了平板近壁面的湍流强度。

本实验尚有待改善,首先流场信息不够精细,要想获得更精细的流场信息,必须采用体积更小、荧光效果更强的示踪粒子、分辨率更高的高速相机,其次本实验中近壁面湍流边界层发展不够充分,需改善雷诺数、平板长度等条件获得充分发展的湍流边界层。

[1] JUNG W J, MANGIAVACCHI N, AKHAVAN R. Suppression of turbulence in wall- bounded flows by high- frequency spanwise oscillations[J]. Physics of fluids a: fluid dynamics,1992, 4(8): 1605-1607

[2] LAADHARI F, SKANDAJI L, MOREL R. Turbulence reduction in a boundary layer by a local spanwise oscillating surface[J]. Physics of fluids, 1994, 6(10): 3218-3220.

[3] CHOI K S, DEBISSCHOP J R, CLAYTON B R. Turbulent boundary- layer control by means of spanwise- wall oscillation[J]. AIAA journal, 1998,36(7): 1157-1163.

[4] CHOI K S, CLAYTON B R. The mechanism of turbulent drag reduction with wall oscillation[J]. International journal of heat and fluid flow, 2001,22(1): 1-9.

[5] BARON A, QUADRIO M. Turbulent drag reduction by spanwise wall oscillations[J]. Applied scientific research, 1995,55(4): 311-326.

[6] DHANAK M R, SI C. On reduction of turbulent wall friction through spanwise wall oscillations[J]. Journal of fluid mechanics, 1999,383: 175-195.

[7] 黄乐萍, 范宝春, 董刚. 槽道湍流壁面展向周期振动减阻机理研究[J]. 南京理工大学学报:自然科学版, 2010, 34(3):361-366. HUANG Leping,FAN Baochun,DONG Gang.Mechanism of drag reduction due to spanwise wall oscillation in turbulent channel flow[J]. Journal of Nanjing University of Science and Technology, 2010, 34(3): 361-366.

[8] 黄乐萍, 范宝春, 董刚.壁面展向周期振动减阻机理及其二维波谱分析[J].船舶力学, 2010,4(14): 325-332. HUANG Leping,FAN Baochun,DONG Gang.Mechanism of turbulent drag reduction with spanwise wall oscillation and spectral analysis of two dimensions[J]. Journal of ship mechanics, 2010,4(14): 325-332.

[9] 许春晓, 吴超, 崔桂香. 壁面在展向作周期运动的槽道湍流的大涡模拟[J]. 计算物理, 2006, 23(5): 537-544. XU Chunxiao,WU Chao,CUI Guixiang.Large eddy simulation of turbulent channel flows withan oscillating wall in spanwise direction[J]. Chinese journal of computational physics, 2006, 23(5): 537-544.

[10] 黄伟希,许春晓,崔桂香,等. 壁面展向周期振动对雷诺应力输运的影响[J]. 全国流体力学青年研讨会论文集(C), 2003: 205-211.

[11] 黄伟希,许春晓,崔桂香,等.壁面展向周期振动的槽道湍流减阻机理的研究[J]. 力学学报, 2004, 36(1): 24-30. HUANG Weixi,XU Chunxiao,CUI Guixiang,et al.Mechanism of drag reduction by spanwise walloscillation in turbulent channel flow[J].Chinese journal of theoretical and applied mechanics, 2004, 36(1): 24-30.

[12] 董海军, 卢晓青, 张浩鹏,等. 壁面展向振动流体特性的分子动力学[J].系统仿真学报, 2015, 27(9): 1680-1686. DONG Haijun,LU Xiaoqing,ZHANG Haopeng,et al.Molecular dynamics of fluid properties by normal oscillations of solid wall[J]. Journal of system simulation, 2015, 27(9): 1680-1686.

[13] 吴文堂, 洪延姬, 范宝春,等. 槽道湍流在确定分布的展向电磁力调制下的壁面减阻机理[J]. 推进技术, 2015, 36(8): 1135-1142. WU Wentang,HONG Yanji,FAN Baochun,et al.Turbulent drag reduction in channel flow modulated by spanwise Lorentz force with steady distriction[J]. Journal of propulsion technology, 2015, 36(8): 1135-1142.

[14] 郭春风, 范宝春. 流向传播的展向电磁力壁面减阻[J]. 船舶力学, 2014,18(7): 754-759. GUO Chunfeng,FAN Baochun,Turbulent drag reduction along streamwise direction via a transverse wave travelling induced by spanwise Lorentz force[J]. Journal of ship mechanics, 2014, 18(7): 754-759.

[15] 刘璐璐, 张军,翟树成,等. 湍流边界层流向-展向平面中发卡涡的涡迹特征[J]. 水动力学研究与进展, 2015,30(4): 376-381. LIU Lulu,ZHANG Jun,ZHAI Shucheng,et al.Hairpin vortex signature in the streamwise- spanwise plane in a turbulent boundary layer[J].Journal of hydrodynamics, 2015(4): 376-381.

[16] 徐驰, 周本谋,黄亚冬,等. 展向电磁力控制鱼雷模型湍流边界层减阻实验研究[J]. 中国科技论文, 2015,10(17): 2033-2035. XU Chi,ZHOU Benmou, HUANG Yadong,et al.Experimental study on the drag reduction of torpedo model controlled by spanwise electromagnetic force[J].China sciencepaper, 2015, 10(17): 2033-2035.

[17] 李山, 杨绍琼, 姜楠,等. 沟槽面湍流边界层减阻的TRPIV测量[J]. 力学学报, 2013,45(2): 183-192. LI Shan, YANG Shaoqiong,JIANG Nan, et al. TRPIV measurement of drag- reduction in the turbulent boundary layer over riblets plate[J]. Chinese journal of theoretical and applied mechanics, 2013,45(2): 183-192.

[18]全鹏程, 易仕和, 武宇,等.激波与层流/湍流边界层相互作用实验研究[J]. 物理学报, 2014, 63(8): 243-247. QUAN Pengcheng, YI Shihe, WU Yu,et al. Experimental investigation of interactions between laminar or turbulent boundary layer and shock wave[J]. Acta phys.sin., 2014, 63(8): 243-247.

[19]王洪平, 高琪, 王晋军,等. 基于层析PIV的湍流边界层涡结构统计研究[J]. 中国科学:物理学 力学 天文学, 2015, 45(12): 73-86. WANG Hongping, GAO Qi,WANG Jinjun,et al. The statistical study of vortex structure in turbulent boundary layer flow based on Tomographic PIV[J]. Scientia sinica:physica,mechanica & astronomica, 2015, 45(12): 73-86.

Experimentalstudyofeffectsofspanwiseoscillationofwallondragreductionandturbulentflows

DENG Fei1, ZHAO Jianpeng1, DONG Haijun2, SHAN Mingming2, TANG Yong1

(1.School of Marine Science and Technology, Northwestern Polytechnical University, Xi′an 710072, China; 2.School of Mechanical Engineering, Northwestern Polytechnical University, Xi′an 710072, China)

10.11990/jheu.201605044

http://www.cnki.net/kcms/detail/23.1390.u.20170427.1523.116.html

O357.5

A

1006- 7043(2017)09- 1380- 05

2016-05-14. < class="emphasis_bold">网络出版日期

日期:2017-04-27.

国家自然科学基金项目(51275412).

邓飞(1973-), 男, 副教授; 赵建鹏(1989-),男,硕士研究生.

赵建鹏,E- mail:zhaojianpeng6992@163.com.

本文引用格式:邓飞,赵建鹏,董海军,等. 展向周期振动平板的近壁面流场特性实验[J]. 哈尔滨工程大学学报, 2017, 38(9): 1380-1384.

DENG Fei, ZHAO Jianpeng, DONG Haijun, et al. Experimental study of effects of spanwise oscillation of wall on drag reduction and turbulent flows[J]. Journal of Harbin Engineering University, 2017, 38(9): 1380-1384.

猜你喜欢

摩擦阻力边界层湍流
一维摄动边界层在优化网格的一致收敛多尺度有限元计算
空间机构用推力滚针轴承摩擦阻力矩分析
航空发动机起动过程摩擦阻力矩计算分析
Bakhvalov-Shishkin网格上求解边界层问题的差分进化算法
“湍流结构研究”专栏简介
基于HIFiRE-2超燃发动机内流道的激波边界层干扰分析
重气瞬时泄漏扩散的湍流模型验证
超大型集装箱船靠泊分析
一类具有边界层性质的二次奇摄动边值问题
湍流十章