大迎角细长体绕流背涡结构与气动特性分析

结合实验成果，对大迎角细长体绕流结构及其气动特征进行了分析。利用势流理论，探讨了大迎角下非对称侧向力沿物体轴线的变化特征与分离背涡结构之间的关系，得出截面侧向力系数沿轴向呈现正弦型曲线的变化特征是由轴向旋涡不断脱落所形成的多涡系统诱导的结果。对于三涡以上的多涡系统中，对称的物面压强分布并非一定对应对称的涡结构。     

不同网格拓扑对细长体绕流计算的影响研究

分别运用扇形(Fan)、阶梯形(Ladder)、交界面形(Interface)网格对细长体小攻角对称、大攻角对称、大攻角非对称绕流流场进行了数值研究．通过涡核位置、涡簇显示、物面压力分布、轴向力分布等的计算结果比较了三种网格的计算精度．数值实验表明：细长体分离涡流场对边界层网格非常敏感,应严格控制边界层网格的正交性;随着攻角增大,流场对网格特性的敏感性有增高的趋势;阶梯形网格可能会对流场带入非物理性扰动,交界面网格对流场捕捉有不连续现象;将三种网格得到的物面压力、侧向力、流动分离位置与实验值进行对比,发现扇形网格误差最小、交界面网格误差最大;大攻角非对称流动时,扇形网格计算的侧向力有整体向细长体头部压缩的趋势,涡脱落位置靠前,第二个及第三个极值更大,说明非对称现象有向尾部发展的趋势．     

Reynolds数对大攻角细长旋成体绕流滚转角效应影响的数值研究

使用低耗散的Roe格式,数值模拟了Reynolds数(Re)对大攻角细长旋成体绕流滚转角效应的影响.模型头部加了几何小扰动块以引发流场的不对称.在较大的Re数(Re=10 5)下,本文的计算结果与实验是相符的,此时细长体的滚转会导致双稳态、双周期现象,即侧向力随滚转角呈现类似方波形式的双周期变化,方波中侧向力基本保持不变的状态对应于流场的正则态,且两个正则态的侧向力方向相反,方波中侧向力基本保持不变的状态对应于流场的正则态,且两个正则态的侧向力方向相反;而在较小的Re数(Re=4 000)下,如果扰动足够大,细长体的滚转将导致不同的双稳态现象,此时两个正则态的侧向力方向相同,而在较小扰动下双稳态现象不再出现;Re数更小时(Re=1 000),即使在较大的扰动下,双稳态现象也不再出现,侧向力随滚动角仍是连续变化的.本文的计算结果表明,Re数越小,流场对头部扰动的感受性越弱.     

矢量喷流对细长体大迎角非对称流动影响研究

采用测压方法研究了矢量喷流对细长旋成体大迎角非对称流动的影响特性.实验结果表明:矢量喷流对细长旋成体大迎角非对称侧向力有明显的抑制作用,该抑制作用是通过喷流诱导作用,改变其空间绕流涡系结构的分布来实现的,但是矢量喷流的存在并不能改变大迎角机身空间绕流涡系的本质结构;随着迎角的增大,矢量喷流对细长旋成体大迎角非对称流动的影响区域不断前移,甚至影响到头部;随着喷流落压比的增加,矢量喷流对细长旋成体大迎角非对称侧向力的抑制作用加强,但当喷流落压比达到临界落压比后(即喷管出口处达到设计马赫数时),喷流影响作用将不会随喷流落压比的增加而改变.     

细长体后部非定常超空泡研究

采用积分方程方法,研究了轴对称细长体后部非定常超空泡问题.应用时间有限差分离散化方法,对积分方程进行了求解.以细长锥体空化器为例,文中分别给出了当锥角和空化数改变(简称扰动)时,空泡长度和形状的变化规律.当流场周期扰动时,分析计算了超空泡的尺度变化.分别采用本文方法和理论公式,对空泡长度与空化数的关系曲线进行了对比.数值结果表明,扰动周期越短,空泡长度的变化越小;在相同的扰动频率下,空泡越长,时间滞后越长;空泡长度相同时,扰动频率越高,时间滞后越长.在高频脉冲扰动下,有脉冲波形沿着空泡表面传播,其传播速度为来流速度.在周期小扰动情况下,扰动波形沿着空泡表面传播,传播速度也是来流速度.本文得到的数值结果为水下航行体空化器的分析和设计提供参考作用.     

细长体水下运动空化流场及弹道特性实验

为了研究细长体水下高速运动时空泡的产生、闭合及脱落特性,及影响细长体空泡形态及弹道特性的复杂因素等,初步开展了细长体模型水下高速运动的实验研究,分析了不同初始空化数下细长体模型在水中高速运动的一系列流动现象,重点研究了空泡的发展、闭合、尾部回射流和尾部脱落特性,以及轴对称细长体模型弹道特性与空泡形态变化之间的关系和转动特性随时间的变化历程等。结果表明:细长体水下高速运动时形成超空泡,空泡头部光滑透明,尾部凝结有汽水混合物且有交替脱落的含气漩涡;初始空化数对细长体的速度衰减有所影响;受初始扰动影响,细长体水下运动伴随有绕头部的转动且初始扰动影响细长体俯仰角随时间的变化历程。     

不同剪切率来流作用下柔性圆柱涡激振动数值模拟

采用浸入边界法对细长柔性圆柱在线性剪切流条件下的涡激振动进行三维数值模拟.细长柔性圆柱振动采用三维索模型模拟,其两端铰接,质量比为6,长细比为50,无量纲顶张力为496.来流为线性剪切流,剪切率从0到0.024变化,最大雷诺数为250.研究发现:剪切流作用下柔性立管横流向振动表现为驻波模式,而顺流向振动表现为行波与驻波混合模式.随着剪切率增大,振动频谱呈现多频响应,振动能量逐渐向低频转移.阻力系数平均值随着展向变化,脉动阻力系数和升力系数的均方根值均表现为"双峰"模式.流固能量传递系数沿立管轴向的分布表明,振动激励区集中于高流速区,而振动阻尼区多位于低流速区.剪切率较小时,圆柱的泻涡为平行交叉模式;剪切率较大时,圆柱的泻涡为倾斜泻涡模式,且由于泻涡频率沿立管轴向变化,尾流发生涡裂现象,形成泻涡频率不同的胞格结构.     

微扰动对液体螺旋槽机械密封性能的影响

针对非接触式液体螺旋槽机械密封,分析了周期性轴向微扰动和角微扰动对密封环端面间液膜厚度的影响规律,建立了单周期螺旋槽液膜模型,采用数值法求解时变雷诺方程,研究了微扰动对液膜承载力和泄漏率等密封特性参数的影响.结果表明,膜厚对液膜承载力和泄漏率的影响显著,膜厚增加,液膜承载力减小,泄漏率增大.液膜承载力、泄漏率的变化幅度和频率主要受到轴向扰动的影响.角扰动造成沿圆周方向分布的液膜承载力不均匀而使液膜稳定性变差.     

并列圆柱绕流的格子Boltzmann数值模拟

采用非均匀不可压格子Boltzmann模型对低雷诺数下并列圆柱绕流进行了数值模拟,给出了数值计算结果,分析了间距g对圆柱尾流及升力、阻力的影响,并在此基础上得到了4种尾迹模式.此外,研究了流场的初始扰动对流动分岔现象的影响,发现在适当的扰动下可以很快得到同步同相的尾流.对Re=160和200下圆柱的升、阻力进行了对比,结果表明升力和阻力受间距g的影响大于雷诺数.     

不同剪切率来流作用下柔性圆柱涡激振动数值模拟

采用浸入边界法对细长柔性圆柱在线性剪切流条件下的涡激振动进行三维数值模拟。细长柔性圆柱振动采用三维索模型模拟,其两端铰接,质量比为6,长细比为50,无量纲顶张力为496。来流为线性剪切流,剪切率从0到0.024变化,最大雷诺数为250。研究发现:剪切流作用下柔性立管横流向振动表现为驻波模式,而顺流向振动表现为行波-驻波混合模式。随着剪切率增大,振动频谱呈现多频响应,振动能量逐渐向低频转移。阻力系数平均值随着展向变化,脉动阻力系数和升力系数的均方根值均表现为“双峰”模式。流固能量传递系数沿立管轴向的分布表明,振动激励区集中于高流速区,而振动阻尼区多位于低流速区。剪切率较小时,圆柱的泻涡为平行交叉模式;剪切率较大时,圆柱的泻涡为倾斜泻涡模式,且由于泻涡频率沿立管轴向变化,尾流发生涡裂现象,形成泻涡频率不同的胞格结构。      

Prediction of asymmetric vortical flows around slender bodies using Navier-Stokes equations

Steady and unsteady asymmetric vortical flows around slender bodies at high angles of attack are solved using the unsteady, compressible, this-layer Navier-Stokes equations. An implicit, upwind-biased, flux-difference splitting, finite-volume scheme is used for the numerical computations. For supersonic flows past point cones, the locally conical flow assumption has been used for efficient computational studies of this phenomenon. Asymmetric flows past a 5° semiapex-angle circular cone at different angles of attack, free-stream Mach numbers, and Reynolds numbers has been studied in responses to different sources of disturbances. The effects of grid fineness and computational domain size have also been investigated. Next, the responses of three-dimensional supersonic asymmetric flow around a 5° circular cone at different angles of attack and Reynolds numbers to short-duration sideslip disturbances are presented. The results show that flow asymmetry becomes stronger as the Reynolds number and angles of attack are increased. The asymmetric solutions show spatial vortex shedding which is qualitatively similar to the temporal vortex shedding of the unsteady locally conical flow. A cylindrical afterbody is also added to the same cone to study the effect of a cylindrical part on the flow asymmetry. One of the cases of flow over a cone-cylinder configuration is validated fairly well by experimental data.     

Some Observations of Vortex Breakdown in a Confined Flow with Solid Body Rotation

The occurrence of breakdown in slender vortex flows as a ``bubble'' or ``spiral'' pattern depends on the degree of radial deflection of the vortex core from its original axis as shown in [1]. A smooth transition from a bubble to a spiral-type ``mode'' can be forced by inducing a small asymmetric disturbance which led to the conclusion, that the patterns do not represent different fundamental modes of breakdown. The subject presented herein addresses the following question: how does breakdown evolve in a swirling flow in which the vortex core is forced on a straight axis? In addition, what is the effect of turbulent inflow conditions? This type of vortex conditions is achieved in a spinning tube flow. The swirl is introduced at the entrance of the rotating tube with a honeycomb package and maintained by the viscous action in the boundary layer of the spinning tube. A diffuser at the end induces an adverse pressure gradient to force the breakdown. Flow visualization experiments are carried out to characterize the nature of breakdown over a range of different flow conditions. For some selected characteristic stages, detailed velocity fields were obtained using the method of Digital Particle-Image-Velocimetry (DPIV). The results show, that for the range of parameters investigated, breakdown is initiated at Rossby-numbers below a critical value of Ro ≈ 0.6 similar to those observed in other experiments. The bursted part of the vortex has a near axi-symmetric slender conical shape containing approximately stagnant flow. Its downstream end is characterized by a jump-like contraction where the flow evolves into a jet with enhanced swirl on the axis. It is only in this region downstream of the jump-like contraction that asymmetric instabilities and wavy flow patterns could be observed. Perturbations caused by them travel upstream but do not change the near-axisymmetric shape of the bursted part of the vortex.     

The combined effects of non‐planarity and asymmetry on primary and secondary flows in the small bronchial tubes

The laminar flow in the small bronchial tubes is quite complex due to the presence of vortex‐dominated, secondary flows. In this paper, we report the results of a numerical investigation of the simultaneous effects of asymmetric and non‐planar branching on the primary and secondary flows in the small bronchial tubes, i.e. generations 6–12. We simulate steady‐state inspiratory flow at a Reynolds number of 1000 in three‐generation, asymmetric planar and non‐planar bronchial tube models. The non‐planar model was defined by applying a 90° out‐of‐plane rotation to the third‐generation branches. A detailed mesh refinement study was performed in order to demonstrate mesh independence. Significant differences were observed between flows in the planar and non‐planar models. An uneven mass flow distribution was observed in the non‐planar model in contrast to the evenly distributed mass flow in the planar model. The secondary flows created symmetric vortex patterns in the planar model, whereas vortex symmetry was lost in the non‐planar model. These results illustrate the importance of incorporating asymmetry in addition to non‐planarity in the geometric models. Copyright © 2008 John Wiley & Sons, Ltd.     

Experimental investigation of side force control on cone-cylinder slender bodies with flexible micro balloon actuators

Side forces on slender bodies of revolution at medium to high angles of attack (AOA > 30°) has been known from a large number of investigations. Asymmetric vortex pairs over a slender body are believed to be the principle cause of the side forces. Under some flight conditions, this side force may be as large as the normal force acting on the slender body. In this paper, experimental results are presented for side force control on a cone-cylinder slender body by using microfabricated balloon actuators. The micro balloon actuators are made of polydimethylsiloxane (PDMS) elastomer by using micro molding techniques. They can be packaged on curve surfaces of a cone-cylinder slender body. As actuator is actuated, the micro balloon actuator inflates about 1.2 mm vertically, which is about 2.4% of the cylinder diameter D (=50 mm) of the cone-cylinder slender body. Micro balloon actuators are actuated at different roll angles of a cone-cylinder slender body. Aerodynamic force measurement results indicate the effects of micro balloon actuators vary at different actuation locations on the cone-cylinder slender body. The side forces can be significantly reduced if the actuators are actuated in the weak vortex side (the side corresponding to the asymmetric vortex which is far from the surface) and actuation angles are located at about 50–60° (the actuation angle here is measured from stagnation line of the incidence plane toward weak vortex side direction). Significant changes are noticed from the surface pressure, as well as leeside vortex flow field, measurement. Micro balloon actuators change nose shapes of the slender body which decide adverse-pressure-gradient values and directly influence the origin of the separation lines and characteristics of the separated vortices over the leeside surface.     

Numerical investigation of the mechanism of asymmetric vortex flows

A robust iterative method suitable for the numerical simulation of high angle of-attack vortex flows is established based upon the multiple line vortex model (MLVM). With symmetric or asymmetric positions of separation lines given, the first converged solution at an angle of attack as high as 60 degree is obtained by means of the present method. Numerical experiments for a tangent-ogive forebody indicate the viscous onset mechanism of asymmetric vortex flows over a body of revolution at high angles of attack and zero sideslip.     

Supersonic flow past a slender cone with asymmetric nosetip

Supersonic flow of an inviscid gas that does not conduct heat past a slender cone with asymmetric nosetip has been investigated at zero and nonzero angle of attack. The flow is calculated by the Babenko-Voskresenskii numerical method. It is shown that the asymmetry of the nosetip has a strong influence on the distribution of the parameters of the gas over the surface of the slender cone even at a large distance from the nosetip. The results of the calculations are compared with experimental data.     

非对称涡流机理的数值研究

基于多重线涡模型(MLVM),本文建立了一种计算机时少、收敛性强的适用于大迎角涡流势流数值模拟的迭代算法,在给定物面上对称或不对称分离线位置条件下,首次得到了迎角大到60°的收敛解,用本文算法对一切拱头体进行的数值实验表明,旋转体在大迎角零侧滑时产生菲对称涡流的机理本质上是粘性的。     

Experimental study on the wake behind tapered circular cylinders

The flow around tapered cylinders can act as basic models for numerous bluff body flows with a spanwise variation of either the body shape or the inflow conditions. The well-known vortex street is influenced by strong three-dimensional effects from the spanwise variation of the shedding frequency, namely oblique vortex shedding and vortex dislocations. Stereo-PIV was chosen to study these phenomena, since it allows analyzing planes with the full three-component, instantaneous velocity fields and local, time-dependent variations in the same setting. Hence, detailed aspects of the vortex dislocation phenomenon are presented. Single vortex dislocation events are presented through the local variation of the three measured velocity components u, v and w. Longer time-series reveal both period and location of these dislocation events, as well as quantity and sizes of the cells of constant shedding velocity in between them. The influence of the Reynolds number and the cylinder aspect ratio on the vortex cells could be shown. The analysis of the vortex shedding behavior shows good agreement with previously published results. At the same time, the applied PIV technique provides more spatial information than point-based measurements and offers insight into a Reynolds number range that is currently out of reach of Direct Numerical Simulations.