On Nonlinear-Evolution of C-type Instability in Nonparallel Boundary Layers

The process of evolution, especially that of nonlinear evolution, of C-type instability of laminar-turbulent flow transition in nonparallel boundary layers are studied by means of a newly developed method called parabolic stability equations （PSE）. Initial conditions, which are very important for the nonlinear problem, are investigated by computing initial solution of the harmonic waves, modifying the mean-flow-distortion, and giving initial value of TS wave and its subharmonic waves at initial station by solving linear PSE. A numerical method with high-order accuracy are developed in the text, the key normalization conditions in the PSE are satisfied, and nonlinear PSE are solved efficiently and implemented stably by the spatial marching. It has been shown that the computed process of nonlinear evolution of C-type instability in Blasius boundary layer is in good agreement with the experimental results.     

多重网格法在混合网格中的应用

设计了一种高效的混合网格聚合方法,在此基础上应用多重网格技术来加速Navier-Stokes方程的收敛过程.方程求解采用Jamson的中心格式有限体积法.对NACA0012翼型和ONERA M6机翼的粘性绕流进行了数值模拟,算例结果与试验结果符合良好,加速收敛效果显著.     

微型旋转翼与拍动翼气动力特性的比较

 用拟压缩性法求解三维非定常不可压N-S方程,研究了旋转翼及拍动翼的气动力和流场。当昆虫正常悬停时,如果翅膀旋转,产生的平均升力以及气动效率并不比拍动时的差,甚至还略好,因此,厘米级微型飞行器可以采用旋转翼,这比拍动模式容易实现。另外,小迎角旋转时升阻比可达3.1（在15°迎角,雷诺数4 000情况下）,气动效率较高而且升力系数不小,多个桨叶可以产生很大的升力。因此,在悬停时,多个桨叶小迎角旋转是厘米级微型飞行器的一个很好的方案。     

跨声速涡轮叶栅三维流场的数值模拟

以具有试验测量数据的一个高压涡轮动叶中截面平面叶栅和一个高压涡轮导向器环形叶栅为对象,采用计算流体动力学(CFD)软件CFX-TASCflow,对两个涡轮叶栅的跨声速三维粘性流场进行了数值计算,验证了CFX-TASCflow软件在航空燃气涡轮设计计算中的可靠性,并对所研究的跨声速涡轮叶栅粘性流场进行了分析.      

共轴刚性旋翼桨毂阻力特性及流动机理

共轴刚性旋翼直升机在高速飞行时,桨毂流动复杂、分离强、阻力大。为明晰其阻力特性和流动机理,采用CFD方法针对已完成风洞试验的共轴桨毂组合模型进行数值模拟研究,获得了桨毂组合模型各单独部件的阻力、表面流动和空间流场特征,阐明了产生阻力最大的部件和影响阻力的主要因素,揭示了中间轴整流罩和塔座设计参数的减阻机制。分析结果表明:上、下旋翼桨毂是产生阻力的主要部件;中间轴和塔座的分离尾流对桨毂表面流动产生较大的干扰作用,使桨毂整流罩表面受干扰区域产生气流分离;具有较缓和逆压梯度的中间轴整流罩和塔座能有效减小分离尾流对桨毂整流罩的干扰,从而降低整个共轴桨毂系统的阻力。     

A transition prediction method for flow over airfoils based on high-order dynamic mode decomposition

This article presents a novel approach for predicting transition locations over airfoils, which are used to activate turbulence model in a Reynolds-averaged Navier-Stokes flow solver. This approach combines Dynamic Mode Decomposition (DMD) with e^N criterion. The core idea is to use a spatial DMD analysis to extract the modes of unstable perturbations from a steady flowfield and substitute the local Linear Stability Theory (LST) analysis to quantify the spatial growth of Tollmien–Schlichting (TS) waves. Transition is assumed to take place at the stream-wise location where the most amplified mode’s N-factor reaches a prescribed threshold and a turbulence model is activated thereafter. To improve robustness, the high-order version of DMD technique (known as HODMD) is employed. A theoretical derivation is conducted to interpret how a spatial high-order DMD analysis can extract the growth rate of the unsteady perturbations. The new method is validated by transition predictions of flows over a low-speed Natural-Laminar-Flow (NLF) airfoil NLF0416 at various angles of attack and a transonic NLF airfoil NPU-LSC-72613. The transition locations predicted by our HODMD/e^N method agree well with experimental data and compare favorably to those obtained by some existing methods (LST/e^N or $\gamma -\stackrel{￣}{R{e}_{\theta \text{t}}}$). It is shown that the proposed method is able to predict transition locations for flows over different types of airfoils and offers the potential for application to 3D wings as well as more complex configurations.     

数值模拟有外挂物的空腔流动

数值模拟了有外挂物存在的开式空腔的超音速流动。数值计算采用全隐式无分裂方法,湍流模型采用修正的B-L代数模型。研究了有外挂物存在的空腔流动的剪切层自持性振荡,并与无外挂物存在时的现象进行类比。压强与时间的变化曲线、空腔底部流向的静压系数分布以及空腔后壁面法向的静压系数分布表明,空腔内的压力脉动与外挂物所处位置是有关的。对于Rossiter振荡模态,数值计算结果能够与实验结果相吻合。     

不规则地形条件下双向DMFT电波传播特性算法研究

 在传统中心差分离散混合傅里叶变换(DMFT)抛物方程模型基础上,针对传统算法计算动态阻抗边界条件时的数值振荡及忽略后向传播与散射造成计算误差增大的问题,提出基于前-后向混合差分方法的递归双向DMFT(TW-DMFT)模型,以一阶前-后向差分方程代替传统的中心差分方程拟合阻抗边界条件,以递归双向传播模型改进传统单向预测模型.仿真分析了双向DMFT模型计算不同边界及媒介条件下单刃峰和实际地形的电波传播分布特性,与传统算法及几何绕射理论(GTD)进行比较,新算法提高稳定性及精度的同时弥补了传统算法的不足,更具有普适性.     

使用Gao-Yong湍流方程计算翼体角偶流动

Gao-Yong湍流方程基于侧偏统计平均方法保留了湍流脉动量的一阶统计信息,并引入加权漂移速度对称性及正交各向异性,导出漂移流连续方程、动量方程和机械能方程,最后依据湍流物理的唯象论,使用动量传输链概念模化封闭了整个方程组.大量算例验证了此方程对广泛范围的复杂湍流问题的适定性.为了进一步的工程应用,将其加入开源计算流体力学软件OpenFOAM(Open Field Operation and Manipulation)中,并就翼体角偶流动进行了数值模拟.翼体角偶流动具有三维分离、马蹄涡等十分复杂的流动特征,通过与实验测量结果的对比研究表明,Gao-Yong湍流方程不论在定性还是定量上,都成功捕捉到这类复杂流动现象的主要特征,得到了令人满意的结果.     

旋翼桨-涡干扰气动特性计算及参数影响研究

建立了一个适用于旋翼桨-涡干扰(BVI)气动特性分析的数值方法。在该方法中,控制方程采用惯性坐标系下的非定常Euler方程,以适合于分析BVI流场的特点;为便于前飞流场分区求解和信息传递,使用了嵌套网格方法;对于流场中涡线的模拟,建立了一种适用于有限体积格式的涡引入方法——广义网格速度法,以简化计算。应用上述方法对旋翼BVI流场进行了计算,并与可得到的试验数据进行对比,验证了方法的有效性。通过对比桨叶弦向位置压强的变化得出,在旋翼BVI过程中,气动载荷主要来自桨叶10%弦长内的前缘部分的压强突变。文中还进一步分析了涡强、干扰距离和干扰角度对BVI气动特性的影响。本文的BVI数值计算表明:当涡接近桨叶前缘时,升力达到最大;而涡位于后缘位置时,诱导速度改变桨叶环量分布,破坏了后缘的库塔条件,但随着涡的远离,桨叶表面环量开始重构。