二维翼型结冰数值计算

翼型结冰计算是随着结冰时间步长变化的准定常过程,因此在每个时间步长后都要更新各个模块的计算工作,如此反复直到达总结冰计算时间.在各模块计算中,分别采用了泊松方程法生成网格、求解雷诺时均方程计算流场、拉格朗日追踪法求解水滴轨迹、基于Messinger模型计算冰形、采用变距滤点法重构边界.引入动态链接库技术将各模块集成为结冰计算平台,并将三种冰形(楔形冰、槽状冰及混合冰)计算结果与国外软件及试验数据进行了对比.      

叶片前缘双出口孔射流冷却效率数值模拟

利用数值模拟的方法研究了叶片前缘双出口孔射流冷却效率.该孔型是一种新概念气膜孔.吹风比变化范围为1.0～3.5,主流通道入口湍流度30%.结果表明,圆柱孔射流冷却效率和实验数据吻合较好,随着吹风比增加,单孔射流冷却效率减小.双出口孔有效地改善了叶片前缘冷却效率.吹风比从1.0增加到2.5时,冷却效率显著增加,吹风比为2.5和3.0时的冷却效率没有显著差别,吹风比为3.5时的冷却效率较低.双出口孔射流冷却效率径向分布比较均匀.      

基于两级直角网格结构的三维DSMC算法研究

结合直角网格计算的高效率以及表面非结构网格对飞行器几何外形的精细描述,本文建立了基于两级直角坐标网格结构的DSMC数值模拟算法。一级网格为计算区域离散后的背景网格,二级网格为独立的碰撞网格和取样网格,是在背景网格的基础上根据流场性质的变化经过自适应后获得的亚网格,碰撞分子对在亚网格内选取,以保证选取的分子是"最相邻的"。计算的Orion飞船外形从自由分子流到稀薄过渡流的气动力系数,与DS3V和DAC软件计算结果的一致性非常好,验证了本文算法的可靠性和有效性。同时对过渡区的流动特征和表面气动参数的分布规律进行了分析。     

一种三维结构网格并行算法

基于三维结构网格,以MPICH2作为并行开发环境,研究了对等模式下流场的并行计算问题。在三维结构网格已经生成的前提下,针对计算区域的划分,网格输出格式和数据通信,提出了一种并行算法。算法的核心思想是将需要通讯的其他区域的单元按区号从1到n排列,并获得本区域通信单元的I-J-K编号与其他区域通讯单元的排列序号的映射关系。最后,在IBM BladeCenter HS22刀片机平台上进行了算例验证,结果表明算法具有很好的可靠性与通用性。     

C/E复合材料网格缠绕结构三维建模与模具设计

通过对圆柱螺旋线和圆锥螺旋线参数方程的推导,实现了网格结构及模具的三维建模。三维模型的建立,可实现结构空间协调,缩短设计周期,节省研制经费,提高生产效率。在此基础上提出了可行的金属模具设计方案,找到了切实可行的模具加工方法,并在工程中加以应用,取得了较好效果。     

F-35战斗机气动及隐身特性分析

以F-35战斗机为研究目标,对目标体进行三维外形重建,并对重建后的模型进行气动及隐身特性的计算和分析,首先采用了基于飞机三视图进行轮廓线提取来重建F-35全机理论外形的方法,引入了基于草图跟踪的CATIA三视图的标定,大大地提高了模型重建的精度。其次采用非结构网格对F-35的计算区域进行网格划分,采用EULER方程完成了F-35在亚声速、跨声速及超声速等飞行条件下的流场计算,分析了不同状态下的升力、阻力和大迎角气动特性。最后利用曲面像素法对F-35全机高频雷达目标特性进行了计算,提出了一种基于IGS数据转换格式的隐身计算网格生成方法,比较了不同俯仰角及方位角下的RCS特性曲线,分析了对RCS影响比较大的部件。     

脉冲等离子体推力器放电波形设计评估仿真研究

运用磁流体动力学的方法,对脉冲等离子体推力器推力室工作过程开展三维双温MHD数值仿真,并开展对推力室的工作机理的分析研究。通过该模型计算得到的推进剂烧蚀质量和元冲量与实验结果相符,同时针对不同初始电压、电容以及两种特征波形的情况进行评估,为推力器的放电波形设计提供指导依据。结果显示在推力器波形设计的初期,是采用电流高变化率低能量水平,还是采用电流低变化率高能量水平,需要结合实际,进行优化评估。     

带旋转修正的弹簧-TFI混合动网格方法

网格存在大变形时,传统的超限插值(TFI)动网格方法易于引起网格正交性问题.依据几何关系并考虑其插值特点,提出了一种带旋转修正的TFI动网格方法.对结构网格块进行分块,在子块角点间建立弹簧元,并基于弹簧类比法计算角点位移;子网格内部结点位移采用修正后的TFI动网格方法进行计算,从而形成具有弹簧-TFI混合特征的动网格新...     

应用多块对接结构网格方法的直升机涵道尾桨气动特性分析

基于湍流Navier-Stokes方程,建立了一套直升机涵道尾桨流场及气动特性的CFD分析方法.在该方法中,针对直升机涵道尾桨流场的特点以及动量源方法对网格系统的要求,提出了一种多块对接网格生成方法,分块贴体网格采用求解Poisson方程的方法生成.通过结合涵道尾桨的运动方式、几何特征及气动特征,建立一个包含动量源项的N-S方程的涵道尾桨流场计算方法和迭代流程;为较好地捕捉涵道壁附近存在的旋涡及分离现象,采用了一方程S-A湍流模型.通过对ONERA M6机翼的绕流特性的计算分析,以及对涵道尾桨在悬停状态下的拉力、诱导速度、载荷分布等的计算分析,验证了该CFD方法的有效性.在此基础上,采用该方法进行直升机悬停、侧飞和前飞状态下涵道尾桨流场与气动特性的数值分析,得到了关于涵道尾桨流场和气动特性的一些有意义的结论.     

NACA4412翼型分离涡计算的混合网格策略研究

以NACA4412翼型为例,针对低速大攻角含有分离涡的复杂流动,提出了一种采用推高局部附面层网格的方法来改善分离涡模拟的新网格生成策略。数值算例表明,这种推高各向异性单元区域总层高的混合网格能够较合理地捕捉到分离流动涡区的位置和大小,数值实验和分析表明,这种网格生成策略是合理有效的。