轴流压气机端壁附面层的理论预测

从Ｎ－Ｓ方程出发，推导了轴流压气机端壁边界层动量积分方程，并提出了一种适用范围广泛的叶片力亏损模型，结合边界层的卷吸方程和速度分布表达式，可以方便的预测压气机环壁边界层的发展，与传统的叶片力亏损模型相比，该方法的计算准确性和适用范围均得到提高。     

相对弯度对低雷诺数流动中翼型动态气动力特性的影响

以固定翼微型飞行器为研究背景,研究了相对弯度对低雷诺数流动中翼型动态气动力特性的影响规律。采用Roe迎风差分格式和双时间步迭代方法,数值求解拟压缩性修正不可压Navier-Stokes方程组,给出了数值算法与实验数据的对比验证。以翼型弦长为特征长度,在Re=500~50000情况下,选取不同最大相对弯度和不同最大相对弯度位置的翼型,计算了其等速上仰时的动态气动力,结果表明后者对气动力的影响比较显著,把最大弯度位置布置在翼型弦向40%的地方要比布置在30%和50%两处所获得的动态升阻比大。     

雷诺数对湍流特性的影响

主要对近年来关于雷诺数（Ｒｅ）对湍流特性的影响的研究进行了回顾．研究表明除径向偏斜系数和平坦系数分布外,确实存在一个不受Ｒｅ数影响的区域,但对于不同的物理量这一区域的范围是不同的,即使是同一物理量不同研究者的结果亦不尽相同,显示了近壁区流动的复杂性和测量的困难,激光测速技术（ＬＤＶ）由于其具有测量精度高、不干扰流场、空间分辨率细等优点,近年来广泛地应用于对湍流,特别是对湍流近壁区流动特性及Ｒｅ数对湍流统计特性影响等的精细测量,随着ＬＤＶ自身及其它相关测试技术的进一步改进,ＬＤＶ在了解湍流产生机理、能量转换等方面发挥着越来越重要的作用．     

低雷诺数翼型蒙皮主动振动气动特性及流场结构数值研究

针对低雷诺数(Re)翼型气动性能差的特点,文章通过对翼型柔性蒙皮施加主动振动的方法,提高翼型低Re下的气动特性,改善其流场结构.采用带预处理技术的Roe方法求解非定常可压缩Navier-Stokes方程,对NACA4415翼型低Re流动展开数值模拟.通过时均化和非定常方法对比柔性蒙皮固定和振动两种状态下的升阻力气动特性和层流分离流动结构.初步研究工作表明在低Re下柔性蒙皮采用合适的振幅和频率,时均化升阻力特性显著提高,分离泡结构由后缘层流分离泡转变为近似的经典长层流分离泡,分离点后移,分离区缩小.在此基础上,文章更加细致研究了柔性蒙皮两种状态下单周期内的层流分离结构及壁面压力系数分布非定常特性和演化规律.蒙皮固定状态下分离区前部流场结构和压力分布基本保持稳定,表现为近似定常分离,仅在后缘位置出现类似于卡门涡街的非定常流动现象.柔性蒙皮振动时从分离点附近开始便产生分离涡,并不断向下游移动、脱落,表现为非定常分离并出现大范围的压力脉动.蒙皮振动使流体更加靠近壁面运动,大尺度的层流分离现象得到有效抑制.     

基于低雷诺数串置翼型气动特性研究

基于低雷诺数,应用二维CFD方法,对串置翼型的气动特性进行了数值模拟。比较了串置翼型布局与常规单独翼型的升阻特性。分析了翼差角度对串置翼升阻特性的影响,并且在保持鸭翼±5°偏角不变的同时,研究了两翼之间的距离以及安装的相对高度对整个串置翼型升阻特性的影响。发现上鸭翼与主翼相结合具有相对较高的最大升力系数和临界迎角, 可以显著改善串置翼型的气动特性。     

雷诺数对运输类飞机气动特性影响的试验研究

开展了某背景飞机的高低速变雷诺数风洞试验,并对试验结果进行了分析;结合国外的一些变雷诺数风洞试验结论,给出了大展弦比运输类飞机雷诺数对升力、阻力以及俯仰力矩特性的影响规律.试验研究表明:升力特性方面,在中小迎角下,随着雷诺数增大,迎角相同时升力增加,大迎角时最大升力系数增大,失速迎角增大;阻力特性方面,随着雷诺数的增大阻力减小,并且最小阻力系数随着雷诺数对数的增加基本呈线性减小;俯仰力矩特性方面,随着雷诺数增大机翼后部载荷增大,低头力矩增大,稳定性增强.除此之外,文中还从雷诺数对气动特性影响机理角度出发,分析了这些规律形成的原因.     

无壁面参数低雷诺数非线性涡黏性模式研究

建立了一个低雷诺数的非线性涡黏性湍流模式,该模式的一个显著特性是它不包含壁面参数（如y^ ,n等）,因而特别适用于复杂几何流场的计算,本模式在几种包括回流、分离、激波等典型流动中进行了验证,结果令人满意。     

随机湍流工况低雷诺数风力机翼型优化研究

为改善小型风力机随机湍流工况适应性,以NACA0012翼型为研究对象,采用非嵌入式概率配置点法,获得随机湍流工况下小型风力机叶片翼型运行攻角分布规律;在气动优化中耦合层流分离预测,基于Transition SST模型、拉丁超立方试验设计、Kriging模型和带精英策略非支配排序遗传算法NSGA-II进行高湍流低雷诺数风力机翼型气动优化。结果表明,优化翼型叶片平均风能捕获效率分别提高3.01%和4.76%,标准差分别降低4.76%和14.93%,优化翼型湍流适应性增强。该方法将翼型设计与湍流风况相匹配,为湍流工况低雷诺数翼型及小型风力机设计提供参考。     

利用修正低雷诺数k—ε模型对管道分离流动的数值研究

采用适合于分离流的修正低雷诺数ｋ－ε模型对突扩园管中的流动进行了数值研究，重点是对近壁流动和再附着点附近流动的研究。从理论上证明了所采用的湍流模型既满足近壁渐进行为，又避免了在通过再附着点截面上μｔ恒为０的困难。和已有的实验数据比较表明，所采用的模型对再附点位置及其附近的流动预测能力优于普通低雷诺数ｋ－ε模型。     

超高雷诺数下矩形截面气动特性研究

为研究工程常见的超高雷诺数下矩形柱流场特性和气动特性,以常见的起重机单箱梁截面(梁高h=1m,雷诺数Re=6.85×105)为例,采用雷诺平均RANS Realizable k-ε方法,对中等湍流强度来流下矩形柱绕流进行数值模拟。首先进行网格无关性验证,然后选取宽高比为0.6、0.8、1.0、1.5、2.0共5组不同截面的矩形柱结构进行流场分析与研究。通过对比不同宽高比柱体截面周围的气动特性参数,以及尾流区的尾涡、回流参数等特征量表明:阻力系数、平均背压系数、升力系数均方根、斯托罗哈数等均随着宽高比增大而降低,尾涡宽度和距离随着宽高比增大而减小,回流长度和最大回流速度随宽高比增大而增大。同时,通过与相关雷诺数(Re=104～105)的风洞实验和仿真结果进行对比,发现其在数值上和参考值相差很小,说明高雷诺数对绕流参数的影响很小,可为相关工程设计提供参考。     

The effect of Reynolds number on boundary layer turbulence

A new facility for studying high Reynolds number incompressible turbulent boundary layer flows has been constructed. It consists of a moderately sized wind tunnel, completely enclosed by a pressure vessel, which can raise the ambient air pressure in and around the wind tunnel to 8 atmospheres. This results in a Reynolds number range of about 20:1, while maintaining incompressible flow. Results are presented for the zero pressure gradient flat plate boundary layer over a momentum thickness Reynolds number range 1500–15?000. Scaling issues for high Reynolds number non-equilibrium boundary layers are discussed, with data comparing the three-dimensional turbulent boundary layer flow over a swept bump at Reynolds numbers of 3800 and 8600. It is found that successful prediction of these types of flows must include length scales which do not scale on Reynolds number, but are inherent to the geometry of the flow.     

Treadmilling swimmers near a no-slip wall at low Reynolds number

The motion of a circular treadmilling low Reynolds number swimmer near a no-slip wall is studied analytically. First, the exact solution of Jeffrey and Onishi [Q. J. Mech. Appl. Math., 34 (1981)] for a translating and rotating solid cylinder near a no-slip wall is rederived using a novel conformal mapping approach that differs from the original derivation which employed bipolar coordinates. Then it is shown that this solution can be combined with the reciprocal theorem, and the calculus of residues, to produce an explicit non-linear dynamical system for the treadmilling swimmer's velocity and angular velocity. The resulting non-linear dynamical system governing the swimmer motion is used to corroborate the qualitative results obtained by an approximate model of the same swimmer recently presented in Crowdy and Or [Phys. Rev. E., 81 (2010)].     

Fluid forces on a very low Reynolds number airfoil and their prediction

This paper presents the measurements of mean and fluctuating forces on an NACA0012 airfoil over a large range of angle (α) of attack (0-90°) and low to small chord Reynolds numbers (Re_c), 5.3 × 10³-5.1 × 10⁴, which is of both fundamental and practical importance. The forces, measured using a load cell, display good agreement with the estimate from the LDA-measured cross-flow distributions of velocities in the wake based on the momentum conservation. The dependence of the forces on both α and Re_c is determined and discussed in detail. It has been found that the stall of an airfoil, characterized by a drop in the lift force and a jump in the drag force, occurs at Re_c ? 1.05 × 10⁴ but is absent at Re_c = 5.3 × 10³. A theoretical analysis is developed to predict and explain the observed dependence of the mean lift and drag on α.     

The effect of the Reynolds number on the Reynolds stresses and vorticity in a turbulent far-wake

Using two orthogonal arrays of 16 X-wires, eight in the (x,y)-plane and eight in the (x,z)-plane, the effect of the Reynolds number in a turbulent plane far-wake has been investigated for two values of Re_θ (based on the free stream velocity and the momentum thickness), i.e. 1350 and 4600. It is observed that as the Reynolds number increases the magnitudes of the measured Reynolds stresses increase, as does the size of two-point vorticity correlation iso-contours. Discernible differences are also observed in probability density function, spectra and three-dimensional topologies. The Reynolds number dependence seems to vanish when Re_θ5000.     

A low Reynolds number variant of partially-averaged Navier-Stokes model for turbulence

A low Reynolds number (LRN) formulation based on the Partially Averaged Navier-Stokes (PANS) modelling method is presented, which incorporates improved asymptotic representation in near-wall turbulence modelling. The effect of near-wall viscous damping can thus be better accounted for in simulations of wall-bounded turbulent flows. The proposed LRN PANS model uses an LRN k-ε model as the base model and introduces directly its model functions into the PANS formulation. As a result, the inappropriate wall-limiting behavior inherent in the original PANS model is corrected. An interesting feature of the PANS model is that the turbulent Prandtl numbers in the k and ε equations are modified compared to the base model. It is found that this modification has a significant effect on the modelled turbulence. The proposed LRN PANS model is scrutinized in computations of decaying grid turbulence, turbulent channel flow and periodic hill flow, of which the latter has been computed at two different Reynolds numbers of Re = 10,600 and 37,000. In comparison with available DNS, LES or experimental data, the LRN PANS model produces improved predictions over the standard PANS model, particularly in the near-wall region and for resolved turbulence statistics. Furthermore, the LRN PANS model gives similar or better results - at a reduced CPU time - as compared to the Dynamic Smagorinsky model.     

Lift force on rotating sphere at low Reynolds numbers and high rotational speeds

The lift force on an isolated rotating sphere in a uniform flow was investigated by means of a three-dimensional numerical simulation for low Reynolds numbers (based on the sphere diameter) (Re<68.4) and high dimensionless rotational speeds (Г5). The Navier-Stokes equations in Cartesian coordinate system were solved using a finite volume formulation based on SIMPLE procedure. The accuracy of the numerical simulation was tested through a comparison with available theoretical, numerical and experimental results at low Reynolds numbers, and it was found that they were in close agreement under the above mentioned ranges of the Reynolds number and rotational speed. From a detailed computation of the flow field around a rotational sphere in extended ranges of the Reynolds number and rotational speed, the results show that, with increasing the rotational speed or decreasing the Reynolds number, the lift coefficient increases. An empirical equation more accurate than those obtained by previous studies was obtained to describe both effects of the rotational speed and Reynolds number on the lift force on a sphere. It was found in calcttlations that the drag coefficient is not significantly affected by the rotation of the sphere. The ratio of the lift force to the drag force, both of which act on a sphere in a uniform flow at the same time, was investigated. For a small spherical particle such as one of about 100μm in diameter, even if the rotational speed reaches about 10^6 revolutions per minute, the lift force can be neglected as compared with the drag force.     

A study on the mechanism of high-lift generation by an airfoil in unsteady motion at low reynolds number

The aerodynamic force and flow structure of NACA 0012 airfoil performing an unsteady motion at low Reynolds number (Re=100) are calculated by solving Navier-Stokes equations. The motion consists of three parts: the first translation, rotation and the second translation in the direction opposite to the first. The rotation and the second translation in this motion are expected to represent the rotation and translation of the wing-section of a hovering insect. The flow structure is used in combination with the theory of vorticity dynamics to explain the generation of unsteady aerodynamic force in the motion. During the rotation, due to the creation of strong vortices in short time, large aerodynamic force is produced and the force is almost normal to the airfoil chord. During the second translation, large lift coefficient can be maintained for certain time period and , the lift coefficient averaged over four chord lengths of travel, is larger than 2 (the corresponding steady-state lift coefficient is only 0.9). The large lift coefficient is due to two effects. The first is the delayed shedding of the stall vortex. The second is that the vortices created during the airfoil rotation and in the near wake left by previous translation form a short “vortex street” in front of the airfoil and the “vortex street” induces a “wind”; against this “wind” the airfoil translates, increasing its relative speed. The above results provide insights to the understanding of the mechanism of high-lift generation by a hovering insect. The project supported by the National Natural Science Foundation of China (19725210)     

A three-dimensional simulation of a steady approach flow past a circular cylinder at low Reynolds number

The three-dimensional (3D) unsteady viscous wake of a circular cylinder exposed to a steady approach flow is calculated using a fractional-step finite-difference/spectral-element method. The calculated flow fields at Reynolds numbers of 100 (2D) and 200 (3D) are examined in detail. The flow field at Re = 100 is 2D as expected, while the flow field at Re = 200 has distinct 3D features, with spanwise wavelengths of about 3.75 cylinder diameters. The calculated results produce drag and lift coefficients and Strouhal numbers that agree extremely well with the experimental values. These 3D values at Re = 200 are in better agreement with experimental values than the results of a 2D calculation at Re = 200, which is expected. © 1998 John Wiley & Sons, Ltd.     

Effect of an adverse pressure gradient on the streamwise Reynolds stress profile maxima in a turbulent boundary layer

It is widely accepted that in a turbulent boundary layer （TBL） with adverse pressure gradient （APG） an outer peak usually appears in the profile of streamwise Reynolds stress. However, the effect of APG on this outer peak is not clearly understood. In this paper, the effect of APG is analysed using the numerical and experimental results in the literature. Because the effect of upstream flow is inherent in the TBL, we first analyse this effect in TBLs with zero pressure gradient on flat plates. Under the individual effect of upstream flow, an outer peak already appears in the profile of streamwise Reynolds stress when the TBL continues developing in the streamwise direction. The APG accelerates the appearance of the outer peak, instead of being a trigger.     

Enhancement of finite Reynolds number effects: Inner-outer sublayer interaction in the turbulent boundary layer

Residual Reynolds number effects in the established data for the velocity profile in turbulent boundary layers (and in pipe or channel flows) are found to be remarkably large. We combine two eddy-viscosity models (with overlapping validity in the inertial sublayer) and show (both analytically and numerically) that this enhancement (which involves a viscous correlation length) arises from inner-outer sublayer interaction.