浸入边界法模拟小圆柱对主圆柱涡脱落的抑制

使用浸入边界法研究了小圆柱对主圆柱涡脱落的抑制.方法中使用非贴体笛儿尔网格,易于处理包含复杂边界的流动问题.采用离散附加力直接加入边界条件方法对虚单元进行重构,使边界条件在浸入边界上精确满足.使用隐式分步法解二维非定常不可压Navier-Stokes方程,通过速度和压力解耦提高计算效率.数值模拟单圆柱绕流及不同位置小圆柱和主圆柱的流动干扰,通过分析流场涡结构和升、阻力系数,得到小圆柱对主圆柱涡脱落的延迟和抑制作用.计算结果与已有实验结论和数值结果对比,计算误差不超过5%,说明浸入边界法可以简单有效地处理圆柱涡脱落抑制这类流动干扰问题.     

计入离心力影响的直升机旋翼翼型结冰数值模拟

建立了一套计入离心力影响的直升机旋翼翼型结冰的数值模拟方法.首先生成围绕翼型的贴体正交网格,然后用Navier-Stokes（N-S）方程求解黏性绕流流场.在此基础上,利用拉格朗日法建立水滴运动方程.其中,为提高计算效率,提出了结合位移矢量的水滴所处单元寻觅方法.最后,结合桨叶工作特点,发展了一种计入离心力影响的三维结冰模型.通过与桨叶结冰实验的对比,验证了本文结冰预测方法的可靠性.对比常规结冰模型,桨叶结冰量减少22.3%;若考虑桨叶的挥舞运动影响,桨叶结冰量进一步减少,表明了离心力及桨叶运动在结冰数值模拟中的重要性.通过不同剖面间的结冰量和冰形对比,分析并获得了桨叶结冰特征.结果表明离心力的影响程度随径向位移的增加而增加,下翼面结冰量随挥舞角的增加而减少.      

A spongy icing model for aircraft icing

Researches have indicated that impinging droplets can be entrapped as liquid in the ice matrix and the temperature of accreting ice surface is below the freezing point. When liquid entrapment by ice matrix happens, this kind of ice is called spongy ice. A new spongy icing model for the ice accretion problem on airfoil or aircraft has been developed to account for entrapped liquid within accreted ice and to improve the determination of the surface temperature when enter- ing clouds with supercooled droplets. Different with conventional icing model, this model identifies icing conditions in four regimes： rime, spongy without water film, spongy with water film and glaze. By using the Eulerian method based on two-phase flow theory, the impinging droplet flow was investigated numerically. The accuracy of the Eulerian method for computing the water collection efficiency was assessed, and icing shapes and surface temperature distributions predicted with this spongy icing model agree with experimental results well.     

结冰表面粗糙度预测与分析

通过对结冰表面液态水的受力分析,确定了水膜、水珠和细流特征粗糙度的大小及其之间的界限划分,发展和完善了结冰表面的粗糙度模型,开发了耦合该粗糙度模型的多步结冰计算程序;利用该程序计算结冰表面的粗糙度大小和分布特征,并对粗糙度影响下的流动和换热特性进行了计算分析;研究发现,结冰表面的粗糙度在分布是不均匀的,且随着结冰温度的降低,局部粗糙度的最大值降低,但在整个结冰过程中,局部粗糙度的最大值表现出增加的趋势;粗糙度的存在增强了结冰表面的对流换热。     

涡脱位置及温度对涡声效应压力振荡影响

为探索涡-声效应对固体火箭发动机中压力振荡特性的影响,基于VKI (Von Karman Institute for Fluid Dynamics)发动机,通过改变挡板位置与燃气温度,对旋涡脱落引起的压力振荡进行了大涡模拟数值研究.耦合分析表明:挡板位于速度波腹附近,压力振荡最为严重;旋涡能量在输运过程中易于被湍流耗散,靠近喷管的二阶速度波腹处旋涡脱落压力振幅明显高于其它位置.解耦分析表明:温度对旋涡脱落频率影响不大,当旋涡脱落频率与声振频率分离后,压力振幅显著下降.     

飞机热气防冰系统研究

介绍了飞机热气防冰系统的应用背景及基本结构,回顾了其发展历程及国内外研究现状.从防冰表面水滴撞击特性计算、蒙皮外部结冰/溢流水模型建立、热气防冰腔内部结构参数研究以及多物理场防冰表面内外耦合仿真这4个关键问题角度,分析了热气防冰系统的主要研究内容.分析表明:准确预测三维溢流冰形成过程、优化热气防冰腔内部结构参数、确定防冰安全边界和防冰裕度以及冰风洞防冰系统实验的参数缩比等关键研究点是热气防冰系统的发展趋势.      

涡扇发动机短舱结冰试验相似方法

通过分析短舱表面结冰缩比试验相似要求,在保证绕流流场、水滴撞击特性和结冰冰型相似的基础上,分别采用基于气流雷诺数,气流韦伯数和水滴韦伯数3种速度选取方法建立了试验相似准则。利用该准则获得了明冰和霜冰条件下涡扇发动机短舱1/2缩比模型的结冰试验参数,开展了短舱进气道表面结冰冰型预测,对参考模型和缩比模型表面过冷水滴撞击特性、结冰特性、溢流水分布特性进行了分析。结果表明:采用基于气流韦伯数和水滴韦伯数速度选取方法的相似准则,在明冰和霜冰条件下均能保证缩比模型与参考模型表面水滴撞击特性相似,且相似性不受结冰温度影响;溢流水变化趋势与参考模型相似,流动极限相近;缩比模型表面冰型均与参考模型相似,该方法能够为短舱结冰试验参数选取提供依据。      

翼型表面明冰的数值模拟

提出一种翼型表面明冰的数值模拟方法.基于欧拉两相流理论,建立空气、水滴两相流控制方程;提出一种可穿透型壁面边界模拟水滴对翼型表面的撞击;采用考虑了粗糙度影响的边界层积分法求解结冰表面的对流换热特性,基于经典Messinger模型提出明冰数值模拟的热力学模型;利用冰层时间推进法模拟积冰过程,采用积冰法向生长假设生成积冰外形.对NACA 0012翼型在不同环境温度下的明冰冰形进行了预测,并与实验结果进行了对比.      

翼型积冰的数值模拟

本文采用计算流体力学的方法对NACA0012翼型前沿的Rimeice积冰进行了预测。在流场计算中,用有限体积(FV)法对二维定常不可压粘流的时均N-S方程进行离散求解;采用四阶龙戈-库塔(RK)法求解水滴运动方程;假定水滴在与翼型相碰撞的点处完全凝结,并且冰沿着与当地翼型表面法向一致的方向增长,以此来预测积冰的形状。文中同时对迎角为4°时结冰翼型的气动特性进行了分析。     

冰灾对车八岭山地常绿阔叶林区系和群落的影响

研究了2008年年初南方冰雪冻雨灾害对车八岭常绿阔叶林不同区系的冠层植被所造成的影响.结果显示,研究样地内温带分布的乔木受损比例高于热带分布的乔木受损比例.X~2检验说明不同区系间的冠层乔木个体受损比例差异显著(P<0.01);东亚分布及其变型受损个体比例最小,热带亚洲至热带非洲分布及变型、东亚和北美间断分布及其变型、热带亚洲分布及其变型受损个体比例较大.冠层乔木都以重度受损为主,中度受损和轻度受损个体数都远小于前者.CCA多元分析显示,株数密度、平均胸高断面积、坡度、海拔各自对不同物种的受损乔木的分布有显著影响(Monte Carlo test,P=0.010 0).优势种米槠Castanopsis carlesii受损植株的分布与群落的平均胸高断面积显著正相关;亚优势种栲Castanopsis fargesii、鸭公青Neolitsea chuii、尾尖叶柃Eurya acuminata受损植株的分布与坡度显著负相关;亚优势种荷木Schima superb与海拔显著正相关.总之,在相同的生境中能够占据林冠优势层的区系类型往往受灾更严重.