65°三角翼亚音速复杂流场计算和数据可视化

提出一种基于非结构混合网格和有限体积法的有效计算策略,对第二期国际涡流试验项目(second international vortex flow experiment,VFE-2)的尖前缘65°三角翼在马赫数0.4,迎角20.3°,雷诺数2×10~6条件下的亚音速复杂流场结构进行数值模拟,重点探讨了基于计算数据进行该类型复杂涡系干扰表面和空间流场关键特征提取和数据可视化问题.通过与相关试验类比,建立了与先进试验流动显示技术相比拟的定性和定量分析方法,为三角翼这类复杂流场结构的精细分析奠定了技术基础.采用上述方法,细致分析了亚音速三角翼的大迎角复杂旋涡流场结构,得到了与试验一致的结论.研究证实:在大迎角条件下,三角翼流动物理复杂,黏性效应耦合严重,只有通过N-S方程计算才能准确地捕捉主涡和二次涡的发展.     

基于雨燕翅膀的仿生三角翼气动特性计算研究

针对低雷诺数微型飞行器的气动布局, 设计出类似雨燕翅膀的一组具有不同前缘钝度的中等后掠($\varLambda =50^{\circ}$)仿生三角翼. 为了定量对比研究三角翼后缘收缩产生的气动效应, 设计了一组具有同等后掠的普通三角翼. 为了深入研究仿生三角翼布局的前缘涡演化特性以及总体气动特性, 采用数值模拟方法详细地探索了低雷诺数($Re=1.58\times 10^{4})$流动条件下前缘涡涡流结构和气动力随迎角的变化规律. 分析结果表明, 前缘钝度和后缘收缩对仿生三角翼前缘涡的涡流强度和涡破裂位置有显著影响. 相对于钝前缘来说, 尖前缘使仿生三角翼上下表面的压力差增大, 涡流强度也更大, 增升作用也更显著. 相对于普通三角翼构型, 仿生三角翼的前缘斜切使其阻力更大, 但后缘的收缩使涡破裂位置固定在此位置, 因此整个上翼面保持低压, 总的升力更大. 由于小迎角时升力增大更明显, 因此仿生三角翼的气动效率在小迎角时明显大于普通三角翼. 这些结论对于揭示鸟类的飞行机理以及未来微型仿生飞行器的气动布局设计具有重要的研究价值.      

多种因素对三角翼俯仰/滚转运动特性影响的数值研究

耦合求解NS方程和刚体动力学方程数值模拟80°后掠三角翼强迫俯仰、自由滚转双自由度耦合运动特性,研究了转动惯量、轴承机构阻尼、翼面流态以及俯仰运动频率、振幅、平均俯仰角等因素对三角翼俯仰、滚转双自由度耦合运动特性的影响。结果表明:机翼的转动惯量和轴承的机械阻尼显著影响自由滚转的频率和振幅;在转动惯量、轴承摩擦和湍流等多种因素的共同影响下,三角翼的双自由度运动可能会形成台阶形式的振荡曲线;俯仰运动的振幅、频率以及平均俯仰角对强迫俯仰、自由滚转双自由度耦合运动特性存在不同程度的影响。     

细长机翼摇滚的数值模拟及物理特性分析

耦合三维非定常Navier-Stokes方程与Euler刚体运动方程数值研究80°尖前缘后掠三角翼的机翼摇滚问题. 采用高精度的WNND(weighted non-oscillatory, containing no free parameters and dissipative)格式离散流动控制方程、采用时间二阶精度单边差分离散刚体运动方程数值模拟了马赫数为0.35, 攻角为10°, 22°, 30°下三角翼受扰后的自由滚转运 动. 结果表明：22°攻角附近为所给三角翼出现横向不稳定的摇滚运动的临界攻角；当攻角小于临界值时，受扰后的机翼滚转运动收敛，而当攻角大于临界值时，受扰后的机翼滚转运动发散并形成极限环形式的机翼摇滚.     

双三角翼拉升流场特性数值模拟研究

发展了适用于双三角翼大攻角非定常分离流场模拟的数值方法,开展双三角翼拉升运动的动态流场特性研究。通过减缩频率、转轴位置和起始攻角的变化,仔细分析了这些运动参数对动态流场施加影响的物理机制,有利于提高对双三角翼在拉升条件下的非定常特性和流场滞后效应等非线性现象的认识。     

起始攻角附近三角翼涡破裂的演化特性研究

用数值方法求解N-S方程,并用"业迭代"方法保证时间方向二阶精度,数值模拟了75°后掠三角翼非定常绕流,给出了详细的涡破裂起始攻角附近涡破裂的形成、发展和演化过程.涡破裂从开始的螺旋破裂变为向泡状破裂演化的"过渡态",然后发展到泡状破裂,再由泡状破裂变为向螺旋破裂演化的"过渡态",最终发展到典型的螺旋破裂形态.从计算结果中还发现,当泡状破裂形成后破裂点前移速度明显减慢,同时当从泡状破裂向螺旋破裂演变时,破裂点出现后移现象.可从最近相关的实验和数值模拟结果中发现这些现象的存在.     

基于输运方程类空化模型的通气空泡流数值模拟

超空泡流动涉及多相、相变、湍流、可压缩性和非定常特性等复杂情况,亟待发展相应的精细数值模拟方法和计算软件.本文基于均质平衡流理论和输运方程类空化模型,提出了复杂气、汽、水多相的通气空泡流数学模型,采用有限体积方法和SIMPLE算法,建立了求解混合介质的雷诺平均N-S方程、湍流方程、各组分质量分数输运方程的数值模拟方法,自主开发了相应的三维计算程序.模拟了绕圆盘的轴对称通气超空泡流,计算结果与经典的近似公式、半经验公式及试验数据吻合良好,验证了数学模型和数值计算方法能够正确预报通气空泡流场的基本空泡形态和流场动力特性.     

三角翼的双襟翼控涡作用的数值模拟研究

对装有“前端襟翼”和“前缘襟翼”的７４°后掠三角翼的不可压缩流场作了数值模拟,以研究襟翼的旋涡控制作用．数值模拟是用拟压缩性方法求解一般曲线坐标系下的三维不可压缩Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,时间离散用向后Ｅｕｌｅｒ差分,空间无粘项的离散用二阶迎风ＴＶＤ格式,所得的离散方程用对角化形式的近似隐式因子分解格式求解．湍流模型用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数模式．计算了三种平面形状的机翼在迎角范围为１０°～５０°的绕流和气动特性．计算和实验的比较表明,襟翼向下偏转可以推迟旋涡破裂,且对提高机翼的减阻能力、升阻比和改善失速前后的气动特性有明显效果,双襟翼具有更佳的控涡效果．     

三角翼大迎角不可压粘流的数值模拟

研究了人工压缩法拟压缩性系数β的选取，采用函数形式的β有效地加速了收敛过程．采用求解不可压Ｎ－Ｓ方程，对三角翼大迎角绕流进行了数值模拟，得到了与实验吻合很好的结果．分析和讨论了大迎角旋涡流动的复杂物理现象     

三角翼动态大迎角气动力特性数值分析研究

采用数值计算方法,对三角翼从 0°上仰至 90°的动态流场和气动力特性进行了计算,并对俯仰角速度对三角翼流场和气动力特性的影响进行了计算分析。给出了三角翼纵向动态情况下的气动力系数变化,特别是大迎角横侧力矩系数的变化特征。结果表明,随着机翼俯仰角速度的提高,前缘分离涡破裂位置相对滞后,机翼升力和阻力系数明显增加,机翼抵抗旋涡非对称破裂的能力明显增强。     

Advances in the study of separated flows

This paper consists of three parts. The first deals with the separation conditions for three dimensional steady viscous separated flows, in which the behaviour of separated flow described by Navier-Stokes equations and boundary layer equations is studied. The second part involves an application of differential topology to qualitative analysis of flow fields. Here the distribution rule of singular points on the separation line is studied. The last part discusses the numerical method solving Navier-Stokes equations for separated flows. The obtained computational results are analysed by the above mentioned theories and methods.     

三角翼大迎角流场结构和气动力特性的数值分析研究

采用数值计算方法对亚音速三角翼纵向及带有小侧滑情况下的流场结构和气动力特性进行了计算。文中给出了三角翼大迎角纵向情况下气动力、机翼前缘分离涡轴线位置和旋涡破裂位置随迎角的变化规律，以及带有横侧小扰动和小侧滑情况下流场结构的非对称性对气动力的影响。计算结果表明与实验结果符合较好。     

Looking for Turbulence Structures: A Numerical Exploration

The purpose of this paper, dedicated to the memory of our friend Hieu Ha Minh, is to wander through the turbulence universe using deterministic computational tools based on large-eddy simulations (LES) in physical and spectral space. We first briefly recall the subgrid models used, then apply them to mixing layers, jets, separated flows (the backstep in particular, on which Hieu worked a lot) and boundary layers. The influence of compressibility will be also considered. These fine-grain LES allow us to decipher very nicely the intimate vortical structure of turbulence, and predict their statistics. Finally, we discuss of the applicability of these methods to industrial flows. This revised version was published online in August 2006 with corrections to the Cover Date.     

Computation of field structure and aerodynamic characteristics of delta wings at high angles of attack

A numerical investigation of the structure of the vortical flowfield over delta wings at high angles of attack in longitudinal and with small sideslip angle is presented. Three-dimensional Navier-Stokes numerical simulations were carried out to predict the complex leeward-side flowfield characteristics that are dominated by the effect of the breakdown of the leading-edge vortices. The methods that analyze the flowfield structure quantitatively were given by using flowfield data from the computational results. In the region before the vortex breakdown, the vortex axes are approximated as being straight line. As the angle of attack increases, the vortex axes are closer to the root chord, and farther away from the wing surface. Along the vortex axes, as the adverse pressure gradients occur, the axial velocity decreases, that is, A is negativee, so the vortex is unstable, and it is possible to breakdown. The occurrence of the breakdown results in the instability of lateral motion for a delta wing, and the lateral moment diverges after a small perturbation occurs at high angles of attack. However, after a critical angle of attack is reached the vortices breakdown completely at the wing apex, and the instability resulting from the vortex breakdown disappears.     

Properties and Averaged Equations for Flows of Bubbly Liquids

Averaged properties of bubbly liquids in the limit of large Reynolds and small Weber numbers are determined as functions of the volume fraction, mean relative velocity, and velocity variance of the bubbles using numerical simulations and a pair interaction theory. The results of simulations are combined with those obtained recently for sheared bubbly liquids [19] and the mixture momentum and continuity equations to propose a complete set of averaged equations and closure relations for the flows of bubbly liquids at large Reynolds and small Weber numbers.     

Possibilities and Limitations of Computer Simulations of Industrial Turbulent Dispersed Multiphase Flows

We give an overview on the usage of computer simulations in industrial turbulent dispersed multiphase flows. We present a few examples of industrial flows: bubble columns and bubbly pipe flows, stirred tanks, cyclones, and a fluid catalytic cracking unit. The fluid catalytic cracking unit is used to illustrate the complexity of the physical phenomena involved, and the possibilities and limitations of the different approaches used: Eulerian–Lagrangian (particle-tracking) and Eulerian–Eulerian (two-fluid). In the first approach, the continuous phase is solved using either RANS simulations (Reynolds-Averaged Navier–Stokes simulations) or DNS/LES (Direct Numerical Simulations/Large-Eddy Simulations), and the individual particles are tracked. In the second approach, the dispersed phase is averaged, leading to two sets equations, which are quite similar to the RANS equations of single-phase flows. The Eulerian–Eulerian approach is the most commonly used in industrial applications, however, it requires a significant amount of modelling. Eulerian–Lagrangian RANS can be simpler to use; in particular in situations involving complex boundary conditions, polydisperse flows and agglomeration/breakup. The key issue for the success of the simulations is to have good models for the complex physics involved. A major weakness is the lack of good models for: the turbulence modification promoted by the particles, the inter-particle interactions, and the near-wall effects. Eulerian–Lagrangian DNS/LES can play an important role as a research tool, in order to get a better physical understanding, and to improve the models used in the RANS simulations (either Eulerian–Eulerian or Eulerian–Lagrangian).     

Numerical simulation of vortical and coherent structures in compressible jet flows

The compressible flows of plane free jets and jets of the intake-stroke of a rectangular piston-engine model are investigated by numerical simulations. The observed vortical structures appear to be the well-known coherent structures of turbulent shear layers. The simulated structures are compared to experimental data by means of density fields and turbulent statistics taken from different authors. The computed flow depends on physical as well as on numerical parameters. The good agreement with the experimental data is obtained by direct simulation without any turbulence model.     

On Plane-Parallel Separated Flows of an Incompressible Fluid near Flat Boundaries

On the basis of Stokes separated flows, examples of separated flows described by the Navier-Stokes equations of a viscous incompressible fluid are constructed. These flows are represented by series convergent in a certain non-zero neighborhood of a flat contour immersed in the flow. In this neighborhood, the series have the same structure as those for the basic Stokes flows. Examples of the regions in which the series segments chosen give only a slight deviation from the numerical solutions of the Navier-Stokes equations are presented. The comparison between inviscid separated flows (without the no-slip condition on the contour) and viscous flows of the same structure (with the no-slip condition) shows that the viscosity does not play a decisive role in the formation of separation or the type of streamline approach to or departure from the contour.