Laminar flow control research at TsAGI: Past and present

This paper presents a brief review of activities in laminar flow control being performed at the Central Aerohydrodynamic Institute named after Prof. N.E. Zhukovsky (TsAGI). These efforts are focused on the improvement of the existing laminar flow control methods and on the development of new ones. The investigations have demonstrated the effectiveness of aircraft surface laminarization applications with the aim of friction drag reduction. The opportunity of considerable delaying of laminar-turbulent transition due to special wing profile geometry and using boundary layer suction and surface cooling has been verified at sub- and supersonic speeds through various wind tunnel testing at TsAGI and during flying laboratory experiments at the Flight Research Institute (LII). The investigations on using hybrid laminar flow control systems for friction drag reduction were also carried out. New techniques of laminar flow control were proposed, in particular, the method of local heating of the wing leading edge, boundary layer laminarization by means of receptivity control, and electrohydrodynamic methods of boundary layer stability control.     

中国当务之急应转变负利率

我国目前利率水平低于CPI涨幅,处于负利率状态,负利率正在使居民的消费能力下降,财富缩水,扭曲资本的市场化配置功能,加剧通货膨胀预期,更造成社会资金一定程度的错配,助长资产泡沫。通过分析负利率的成因和负面影响,提出可提高定期存款利率、稳定物价和转变经济增长模式来消除负利率。     

自适应遥测信道速率的分包遥测调度方案

在卫星分包遥测静态分配适应单一遥测信道的方案基础上,提出了一种改进型动态自适应遥测信道速率的遥测分包方法,将包传输的控制参数由绝对参数优化为相对参数,从而满足不同信道速率、不同包源速率下的包传输适应性.该方法通过将遥测包周期计数与遥测包静态分配的传输周期的比值作为遥测包传输等待迫切比例系数,用相对的传输迫切比例系数值对全部遥测包进行排序,从大到小组入遥测信道下传.通过这种遥测包调度策略解决了遥测信道速率变换对包传输控制的影响.此外通过调节遥测包传输比例系数也能很好地解决遥测包数据率变化对信道传输的适应能力的影响,如数据率增加时整体降低包传输频率、数据率减少时整体加快包传输频率.     

导引头耦合度对控制系统影响机理分析

以经典三环控制结构自动驾驶仪、纯比例导引律为例,从自动驾驶仪的角度分析了导引头耦合度(隔离度)对控制系统的影响.耦合度通过改变自动驾驶仪内部的结构参数,极大地影响了控制系统性能,从而引发一系列其他相关问题.     

带自由表面水流的分离涡模拟

目前水动力学问题的数值求解仍然广泛采用RANS模型,特别对于大尺度的地表水流的数值模拟,该类模型计算效率较高,但其难以给出湍流场的小尺度涡结构.深入研究复杂地形条件下的水流运动,需要发展诸如LES等的高级数值模型.天然地表水流,地形、边界复杂,雷诺数通常较高,高级数值模拟受限于计算能力的限制,还难以工程应用.RANS和LES混合模型是目前具有工程应用的一种混合数值模型,其已经在CFD领域获得了长足发展.然而对于地表水流运动的数值模拟,相关工作还较少.本文将分离涡模型(DES),即一种RANS和LES的混合模型,应用于带自由表面的地表水流运动,建立了一套数值仿真模型.模型基于有限体积法,水平面内采用非结构计算网格,垂向为结构化网格,对流项离散格式采用二阶TVD格式,并行基于OpenMP语言库.算例表明DES模型有助于揭示复杂地形条件下带自由表面水流的大涡拟序结构.     

基于非结构网格的不可压N-S方程多矩有限体积法

提出了一种基于三角形及四面体非结构网格的有限体积法(FVM),用以鲁棒且精确地求解不可压粘性流动问题.与传统的FVM方法仅将体积分平均值(VIA)作为计算变量的做法不同,本文提出的方法将VIA及点值(PV)同时作为计算变量并在每个迭代步进行计算更新.VIA以通量形式进行计算以确保数值守恒,PV可以通过控制方程的不同形式进行求解更新,无需守恒,因此可以采用非常高效的方法进行求解.将PV作为增加的变量使得紧致网格模板得以实现更高阶精度的重构,而且由此获得的数值模型对于非结构网格变得更鲁棒.本文针对二维/三维的三角形/四面体非结构网格提出了数值格式,给出了几个基准测试算例,验证了本文提出的数值方法在采用非结构网格求解不可压粘性流动问题时的精确性和鲁棒性.     

跨声速串列转子失速机制的数值研究

为了解决跨声速串列转子的低裕度问题,就必须了解跨声速串列转子的流场结构与失速机制.设计了叶尖切线速度为450m/s,负荷系数为0.56的高负荷跨声速串列转子.基于数值模拟的结果,分析了该串列转子在0.5mm叶尖间隙下的叶尖流场结构与失速机制,并在此基础上分别探讨了叶尖间隙和前、后排叶片周向位置对串列转子特性的影响和失速机制的变化.结果表明:前排叶片的叶尖区域是 影响串列转子稳定性的关键;随着叶尖间隙的增加,串列转子的失速机制也发生变化,从前排叶片叶尖区域的尾迹与径向潜流堵塞后排叶片通道转变为前排叶片叶尖泄漏流堵塞;在较大周向相对位置(后排叶片压力面周向远离前排吸力面)的情况下,串列转子获得最好的效果,随着周向相对位置(PP)的增加,失速部位从后排叶片转移至前排叶片.      

蜂窝式轴心通风器油气分离性能计算

为对航空发动机蜂窝式轴心通风器油气分离效率进行研究,建立考虑油气双向耦合的流场计算方法及油滴/壁面相互作用模型,在验证通风阻力及油气分离效率可靠性的基础上,对不同转速、通风流量和环境温度下蜂窝式轴心通风器的油气分离效果进行计算和分析.结果表明:转速的增加会使油气分离效率得到提升,而通风流量和环境温度的增加则导致油气分离效率的降低.蜂窝孔结构的加入对通风阻力影响不大,却对通风器的滑油分离过程起主要作用,计算表明其对滑油分离贡献率在80%以上.      

翼吊发动机安装结构等效建模及其隔振设计

为了研究翼吊发动机安装结构隔振特性并优化其隔振器设计,建立了发动机安装节-吊架-机翼结构理论分析及有限元模型.利用有限元方法进行了模态验证并分析了安装结构的隔振特性.进行发动机3种典型工况下的结构动响应分析确定了振动传递的主路径.基于振动传递路径法研究了隔振器参数和安装位置对安装系统隔振性能的影响规律.结果表明:振动载荷经安装结构后低压转子转频和高压转子转频峰值响应分别降低22.03%和14.65%.低压转子转频振动传递主路径为发动机-前安装节-吊架-机翼,高压转子转频为发动机-后安装节-上连杆-机翼.通过合理设置隔振器位置可以使安装系统隔振率达到50.41%,隔振器的频率比为5和阻尼比为0.25时安装系统隔振率可达70.67%.为了优化整个发动机安装系统的隔振效果,设计隔振器时必须选取合适的安装位置和参数.      

New constitutive relation for eddy viscosity models

A new turbulent constitutive relation was directly derived from Boussinesq's hypothesis and mixing length theory, and then implemented in the standard k-ε model. The performance of this constitutive relation was validated in zero pressure gradient flat-plate boundary layer flow, fully-developed turbulent channel flow and separated flow in a plane asymmetric diffuser. The investigation demonstrated that, this new constitutive relation gave very accurate results in the former two basic cases and provided significant improvement in prediction of separated and reattachment points in the plane asymmetric diffuser. Separation and reattachment points at x/H=7.5 and 29 were calculated accurately in comparison to experimental results, and the static pressure coefficient of 0.82 was very close to large eddy simulation calculation. These results are very encouraging but further verification and extensive application of the new constitutive relation to other two-equation eddy viscosity model are needed.