采用上游扰动体被动控制减阻技术的实验研究

通过风洞测压和水洞流动显示实验,研究了不同直径的扰动体尾迹下圆柱的流动结构和气动特性。结果表明当扰动体和圆柱串列布置时,流动呈现出两种流动模式:空穴流动和尾流撞击模式。当圆柱处于空穴流动模式时,其阻力可以大幅度降低。当扰动体直径为圆柱直径1/2时,圆柱阻力最大降低98%。综合绕流结构和圆柱表面压力分布的变化,指出阻力减小的原因主要是扰动体对圆柱前体的保护作用以及对后体流动的干扰作用,其中后者占主导地位。随着扰动体直径的增加,扰动体对圆柱前体的减阻效果增加。     

合成双射流控制水下圆柱绕流流动分离数值模拟研究

为了探究合成双射流激励器及出口射流参数对圆柱绕流流动分离的控制效果,首次对合成双射流控制水下圆柱绕流流动分离进行了数值模拟。数值计算结果显示:保持激励器出口射流振幅不变的条件下,出口射流频率等于尾迹涡脱落特征频率时,射流控制作用与绕流流场耦合效果最好,控制流动分离效果最佳;保持出口射流频率为尾迹涡脱落特征频率时,在数值计算测试范围内,随着射流振幅的增大,射流对于流场的动量掺混能力增强,控制效果也随之增强。机理分析表明:合成射流位于前驻点的控制,主要通过在圆柱前缘形成虚拟气动外型来达到减阻控制的效果;而合成射流位于后驻点的控制,主要是通过增强回流区的动量掺混来提高回流区抑制分离的能力,从而达到减阻控制的效果。     

亚临界圆柱绕流的DES方法比较

针对地面结冰试验设施喷雾系统尾迹湍流模拟需求,考察了3种脱体涡模拟(DES)方法对三维亚临界圆柱绕流预测的准确性,比较和分析了瞬时流动特征和流场统计量。研究发现:1)从瞬时流动特征来看,k-ω雷诺平均NavierStokes(RANS)控制的区域对剪切应力输运k-ωDES(SSTk-ωDES)求解的准确性有显著影响;2)从流动参数统计量的误差范围来看,回流区长度和流向最小速度的预测质量能在一定程度上反应圆柱绕流数值模拟的准确性;3)综合比较,SSTk-ωDES预测值与试验值和大涡模拟(LES)数据最为吻合,具有应用于地面结冰试验设施喷雾系统尾迹湍流计算的潜力。     

超声速尾迹-剪切流的混合增强

超声速条件下燃料和空气之间的高效混合是超然冲压发动机技术上的主要挑战。基于大涡模拟和流动稳定性分析,针对超声速尾迹-剪切流动开展了混合增强方法研究。尾迹的存在改变了混合层流动的速度剖面,对流动稳定性产生了重要影响,使混合层由三维最不稳定变为二维最不稳定,最不稳定扰动波频率和增长率增大。基于流动稳定性结果引入扰动的混合增强方式依然有效,根据稳定性结果设计了波纹隔板。数值结果表明：二维波纹壁引入的扰动未能增长,不具备混合强化效果,而三维波纹壁引入的扰动能够快速增长,具有混合强化效果,且波纹壁参数越接近最不稳定扰动波参数,混合强化效果越明显。     

折叠V形钝体火焰稳定器原型的稳燃机制研究

为了掌握折叠V形钝体这一新型火焰稳定器的稳燃机制,采用风洞实验与混合雷诺平均/大涡模拟结合的方法对常压条件下10组不同折叠角的折叠V形钝体火焰稳定器原型进行了冷态实验与热态数值模拟研究。其中冷态实验风速10~50m/s,热态数值模拟工况为进气温度700K、来流速度50~150m/s,当量比0.3~1。通过PIV实验测得了折叠V形钝体回流区冷态流动特征,并结合数值模拟方法研究并讨论了折叠V形钝体绕流尾迹区旋涡特征;通过数值模拟研究了燃烧条件下、折叠V形钝体下游槽向驻涡结构发展变化规律。研究初步探明折叠V形钝体总压恢复系数、阻力系数、抗吹熄能力与燃烧效率优于标准钝体的原因,并发现折叠角角度较小的折叠V形钝体,尾迹中心区域无典型回流区,但仍能通过槽向驻涡结构形成驻涡稳燃速度型,该速度型有利于冷态绕流流速的恢复与热态燃气的加速。     

来流扰动对圆盘绕流结构影响的实验研究

主要讨论由水槽流动显示实验得到的圆盘绕流结构在上游小圆柱尾迹干扰下的变化,对圆盘前驻点附近逆压区的流动通过氢气泡显示技术进行观察.由实验结果处理得到圆盘前大涡的结构特征(如周期、最大宽度、在周期内结构随时间的变化等)受扰动几何参数的影响.发现无量纲参数L/d(L为小圆柱中心与圆盘迎风面之间的间距,d为小圆柱直径)对圆盘前驻点区域中大涡的形成和发展有重要的影响.圆柱尾迹中涡通量的增加会在圆盘前产生大尺度的涡对,其周期性演化的频率随L/d增加而增加.     

大型飞机后体流动控制及减阻机理研究

通过数值模拟方法研究了大型飞机机身后体的流动特征和机理,并采用涡流发生器对机身后体流动分离进行控制,分析了其控制减阻的机理,讨论了不同参数对减阻效率的影响。计算结果表明:大型飞机机身后体流动分离是导致巡航阻力增加的一个重要因素;在后体底部附近安装涡流发生器,当来流迎角为负时减阻效果较明显,随着迎角的逐渐增大,减阻效果降低;涡流发生器的尺寸、周向距离、安装在机身后体上的前后位置、安装角度等参数对减阻效果具有显著影响。     

利用上游槽道对角区马蹄涡的控制

在平板上放置圆柱形成角区流动,利用布置在圆柱上游平板上的二维和三维槽道来控制或削弱角区马蹄涡,采用风洞试验和数值模拟开展研究。结果表明,二维和三维槽道均能推迟边界层的分离,使圆柱根部马蹄涡的强度减弱、尺度减小;同时槽道上游压力和逆压梯度均有所下降,槽道下游压力显著升高而逆压梯度总体降低。二维槽道对马蹄涡强度的削弱为61.15%~66.51%,而三维槽道对其削弱为66.65%~80.93%。讨论了三维槽道参数(包括槽道宽度、深度以及其中心线与圆柱中心距离)对控制效果的影响。槽道与圆柱的距离在对马蹄涡的控制中起主导作用。槽道控制的机理是,由于槽道的抽吸效应使得其上游靠近壁面的边界层中涡量较高的流体被卷吸入槽道形成槽道涡,槽道涡由三维槽道输运到下游。同时,随着槽道与圆柱的距离减小,更多的边界层流体流入槽道内。正是上述"槽道效应"使得槽道下游的逆压梯度降低,马蹄涡强度减弱,分离区范围减小。     

绕钝体分离流动的非定常数值模拟研究

利用非定常RNGk-ε模型对三种典型的钝体绕流问题进行了研究,这三种流动分别为:方柱绕流、圆柱绕流和绕立方柱流动。研究表明了采用非定常雷诺平均方法可以研究钝体绕流的非定常流动,但所获得的结果各不相同。对于分离角固定的方柱绕流和绕立方柱流动,数值模拟结果和实验以及大涡模拟结果吻合较好,对于分离角不固定的圆柱绕流,随着雷诺数的增加,数值模拟结果和实验结果相差越大,分析认为这是由于壁面模型无法准确预测流动中不固定的分离角所致。改进了绕立方柱流动的模拟方法,采用非定常模拟方法后获得的结果要比定常模拟方法得到的结果有了根本改善。     

等离子体合成射流对钝头激波的控制与减阻

为了验证等离子体合成射流对超声速流动的流动控制和减阻效果,在考虑热完全气体效应的情况下,工程拟合等离子体热力学属性和输运系数,利用能量源项模型对超声速平板和球头等典型流场结构进行了数值模拟,并取得了与实验较为一致的结果。研究结果表明：对于马赫数为2的超声速平板流动,等离子体合成射流能有效干扰边界层的发展,并诱导产生一系列大尺度结构;对于马赫数为3的球头流场,等离子体合成射流能显著改变激波脱体距离与球头的阻力特性;在放电后第1个周期内,合成射流能使球头平均阻力降低6.3%,而在射流峰值情况下使阻力降低32.0%,激波脱体距离增加2倍,实现了激波控制和流动减阻的预期效果。     

趋于临界马赫数的圆柱跨声速绕流特性分析

通过深化认识趋于临界马赫数Ma_cr的圆柱跨声速绕流特性,明确新型飞行器增升减阻设计的空气动力学理论依据。采用大涡模拟方法数值研究了来流马赫数Ma_∞为0.75和0.85、雷诺数Re为2×10⁵的圆柱跨声速绕流。结果表明:当Ma_∞趋于临界马赫数(Ma_cr≈0.9)时,圆柱的阻力下降且升力系数振荡被抑制;通过力的分解,得知圆柱的阻力减小来自旋涡力的影响,而非可压缩性;圆柱的阻力减小与其背压上升有关;剪切层初始阶段的对流马赫数Ma_c随Ma_∞的增加而增大,而增长率相反,这使得剪切层更为稳定、柱体背压更高。此外,由于Ma_∞=0.85时边界层分离点处的激波和尾迹处的激波向下游推移,使得近尾迹处的湍流脉动减弱,进而导致柱体的表面压力振荡和升力系数振荡被抑制。     

Numerical investigation of transitions in flow states and variation in aerodynamic forces for flow around square cylinders arranged inline

This study focuses on the transitions in flow states around two-, three- and four-inline square cylinders under the effect of Reynolds numbers at two different gap spacing values using the lattice Boltzmann method. For this purpose, Reynolds number is varied in the range 1–130 while two different values of spacing taken into account are gap spacing?=?2 and 5. Before going to actual problem, the code is tested for flow around a single square cylinder by comparing the results with experimental and numerical results of other researchers, and good agreement is found. The current numerical computations yield that for both spacing values and all combinations of cylinders there exist three different sates of flow depending on Reynolds numbers: steady state, transitional state and unsteady state. It is found that the range of Reynolds numbers for these flow states is different for both spacing values. At gap spacing?=?2 the range of Reynolds numbers for each flow state decreases by increasing the number of cylinders while at gap spacing?=?5 opposite trend is observed. The results also show that at gap spacing?=?2 the reduction in drag force is greater than the corresponding reduction at gap spacing?=?5. The maximum reduction in drag force is observed at Reynolds numbers?=?1 at both spacing values. Similarly, at both spacing values and all Reynolds numbers, the maximum reduction in drag force is observed for the case of four-inline square cylinders.     

浸入边界法模拟小圆柱对主圆柱涡脱落的抑制

使用浸入边界法研究了小圆柱对主圆柱涡脱落的抑制.方法中使用非贴体笛儿尔网格,易于处理包含复杂边界的流动问题.采用离散附加力直接加入边界条件方法对虚单元进行重构,使边界条件在浸入边界上精确满足.使用隐式分步法解二维非定常不可压Navier-Stokes方程,通过速度和压力解耦提高计算效率.数值模拟单圆柱绕流及不同位置小圆柱和主圆柱的流动干扰,通过分析流场涡结构和升、阻力系数,得到小圆柱对主圆柱涡脱落的延迟和抑制作用.计算结果与已有实验结论和数值结果对比,计算误差不超过5%,说明浸入边界法可以简单有效地处理圆柱涡脱落抑制这类流动干扰问题.     

ALE方法求解圆柱的涡致振动

数值求解基于拉格朗日-欧拉(ALE)描述的不可压缩流体的N-S方程,计算了较小雷诺数下圆柱的涡致振动现象.N-S方程的对流项和扩散项分别采用三阶迎风紧致格式和四阶中心紧致格式离散,计算网格采用H-O非交错网格系统,并结合分块耦合方法.柱体运动简化为弹簧-阻尼-质量系统,柱体运动方程采用经典龙格-库塔方法求解.通过模拟柱体和流体之间的非线性耦合作用,成功地捕捉到了"锁定"、"拍"和"相位开关"等现象,并与试验数据相吻合.另外,本文详细分析了Re＝200时圆柱尾涡形态、升力、阻力以及圆柱位移等随圆柱自振频率变化的过程,捕捉到尾涡结构变化的频率转折点.     

斜、直圆柱绕流的CFD模拟

采用CFD软件CFX5.5对直圆柱和斜圆柱绕流进行数值模拟,计算了直、斜置圆柱的阻力系数、升力系数、周向平均压力系数及其脉动压力系数的均方根值,并与试验结果进行了比较.计算表明,斜圆柱绕流与直圆柱绕流有很大的差异:斜圆柱绕流的平均阻力系数减小,平均升力系数不等于零,升力系数的频域除包含有传统的卡门涡脱频率外,还包含有许多低频成份.     

吸气条件对圆柱非定常分离流影响的数值研究

本文根据二维可压非定常Ｎ－Ｓ方程，基于Ｂｅａｍ ａｎｄ Ｗａｒｍｉｇｎ及Ｓｔｅｇｅｒ的方法，在圆柱绕流出现周期分离后，施加吸气条件，研究气对于圆柱非定常绕流的影响，计算结果表明，在恰当地区域施加恰当的吸气条件，可以希望通过减小压差阻力来减小总阻力，吸气还有抑制尾迹中非对称流动的作用，从而使圆柱所受的垂直于来流方向的升力的变化幅度减小，变化频率减小，但是吸气会增加磨擦阻力，过份的吸气会导致摩阻增大太     

Study of riblet drag reduction for an infinite span wing with different sweep angles

This paper studies the riblet drag reduction effect for an infinite swept wing under a low Reynolds number using a large-eddy simulation. The results show that the drag reduction ratio is not linear under different sweep angles. The maximum drag reduction ratio in this study is 9.5% for a wing with a 45° sweep angle. The local surface streamline angle and turbulence quantities are calculated to analyze the drag reduction mechanism. The results demonstrate that the riblets considerably suppress the Reynolds stresses above the wing upper surface, while the turbulence kinetic energy in the near wake is increased. A possible relaminarization phenomenon is observed at the middle part of the wing. Quasi-two-dimensional flow structures are observed near the wall, and a peak frequency is considered as the dominant frequency of the region.