民机翼梢小翼气动特性数值模拟分析研究

针对某民机翼身组合体加装不同翼梢小翼——融合式翼梢小翼和翼尖涡扩散器进行对比分析。用软件ICEM生成原始翼身组合体及加装不同小翼的翼身组合体的点对点对接多块网格技术生成高质量的数值计算网格,运用Roe三阶迎风偏置通量差分裂方法和隐式近似因子分解方法求解雷诺平均N-S方程。对两种不同翼梢小翼在巡航点进行数值模拟,得到合理的翼稍小翼几何参数;对比升阻特性数据,两种翼梢小翼都可以提高升力系数,减少阻力系数。相对原始翼身组合体,融合式翼梢小翼升阻比增加9%,翼尖涡扩散器升阻比增加6%;但翼尖涡扩散器在力矩特性上有较大优势。所得结论对民用飞机翼稍小翼设计工作具有较大的工程应用价值。     

蜜蜂悬停飞行控制的仿生力学研究

研究蜜蜂悬停飞行的控制问题.用数值求解N-S方程的方法计算拍动翅体的控制导数;用特征模态分析方法分析控制特性.文中δu,δω,δq和δθ分别为水平方向速度,垂直方向速度,俯仰角速度和俯仰角度的扰动量;δ(φ)和δ(φ)分别为拍动幅度和拍动平均角增量;δα1表示上、下拍迎角同时增加的增量;δα2表示下拍迎角增加(或减小)而上拍迎角减小(或增加)的增量.获得以下结果:(1)悬停飞行时,改变Φ和α1主要产生垂直力的变化;改变φ主要产生俯仰力矩的变化;改变α2主要产生水平力及俯仰力矩的变化.(2)蜜蜂悬停飞行的纵向扰动运动由3个特征模态构成:不稳定振荡模态,快衰减模态和慢衰减模态;为实现稳定的悬停飞行,不稳定振荡模态和慢衰减模态需要控制.为控制不稳定振荡模态,以δ(φ)(或δα2)反馈δu,δq和δθ这3个量的某种组合便可;为控制慢衰减模态,以δ(φ)(或δα1)反馈δω便可.这就是说,该昆虫只需用控制变量δ(φ)和δ(φ)(或δ(φ)和δα2,δα1和δ(φ)1,δα1和δα2),便可稳定地悬停(当然也可4个控制变量都用).     

应用不同湍流模型计算叶栅流场的结果与比较

采用有限体积TVD-LW格式求解N-S方程,使用了B-L模型、高雷诺数k-ε模型和低雷诺数k-ε模型等多种湍流模型。并且用具体算例说明本文所用程序的可靠性,比较了3种的湍流模型的特点。B-L模型最大的优点是简单、使用方便。k-ε方程适应性强,分为低雷诺数k-ε方程和高雷诺数k-ε方程两种。高雷诺数k-ε方程比低雷诺数k-ε方程精度略低,但所需计算时间和存储量都比后者小,网格也可以粗些。k-ε模型对粘性层、分离区有比较高的分辨能力,高雷诺数k-ε模型比较适于一般数值计算。      

基于 N-S 方程的飞机风挡雨刷气动力计算

风挡雨刷在飞行过程中易出现被气流吹离复位位置、飘起现象,阻碍飞行员视野,严重影响飞行安全。选取某型机风挡雨刷作为优化对象,采用了结构、非结构混合网格的办法,在机身壁面附近生成结构化网格,在雨刷机构表面生成非结构化网格,用 Interface 面将两部分结合生成雨刷最终计算模型。利用 Fluent 软件,基于 N-S 方程和 S-A 湍流模型,采用时间平均法对不可压缩流的参数进行时均化,计算了不同飞行状态下,风挡雨刷在复位位置以不同角度停放及增加扰流板时所受气动力,对风挡雨刷的外形及停放角度进行优化,比较分析了优化前后雨刷所受气动力。仿真计算结果表明,雨刷合理的停放角度、在刷臂上加装扰流板能有效引导风挡雨刷迎风面气流的方向,从而减小雨刷迎风面所受气动力。通过试飞验证了数值计算的可靠性及优化方案的可行性,为风挡雨刷预紧力的设置及复位位置的确定提供依据。     

高阶熵条件格式下Euler方程与 N-S方程的混合算法

针对在工程问题中使用N-S方程计算量大的问题,提出只在粘性作用显著的局部区域求解N-S方程,在流场的其余部分采用Euler方程计算的办法来模拟具有复杂外形的流动问题,并提出了一种基于高阶熵条件格式的算子分裂算法.应用所构造的算法对绕迫击炮弹的亚、跨超声速流动进行了模拟计算.计算结果与试验数据的对比表明,所提出的方法是切实可行的.      

共轭梯度方法在求解不可压N-S方程中的应用

采用非结构网格上的SIMPLE(Semi-implicit Method for Pressure-linked Equations)算法、k-ε湍流模型求解了三维不可压Navier-stokes方程,对低速不可压粘性流场进行数值模拟,在求解压强方程时采用共轭梯度方法,求解其它变量方程时,采用预处理的BICG(biconjugate gradients)算法.对NACA0012翼型绕流流场和飞艇绕流流场进行数值模拟并对结果进行分析,取得较好的结果.     

瞬时起动圆柱的初期流动

 本文将ADI方法和Poisson方程直接法结合起来求解涡量-流函数形式的N-S方程,研究瞬时起动圆柱的初期流动(Re=3000,5000,9500)。成功地算出实验观察到的α流动结构、β流动结构,以及这些结构在流动发展过程中交替出现的情况。本文方法效率高,计算结果与实验结果及他人的数值计算结果相符较好。可期望用该方法计算其它复杂流动。     

物面突变引起的二维分离流动的数值模拟

本文用一个单步二阶精度格式,对由物体截面突变引起的二维超声速分离流动进行了数值模拟。计算了近尾迹流场,并详细研究了来流M数和Re数对绕后台阶分离流场结构的影响,得到了与理论和前人计算一致的结果。另外,对超声速绕竖舨流动进行了数值计算,给出了细致的流场结构图画。     

如何使用CFD来确定风洞实验模型的缩尺比例——风洞试验模型设计准则从感性到理性的飞跃

对上百年传统的风洞试验模型缩尺比例的确定方法,提出了疑问,并且从N-S方程出发,指出了这种“传统”方法不符合相似理论,具有很大的随意性和盲目性。因而按“传统”方法给出的模型尺寸所提供的风洞实验气动力数据是不准确的,不能满足当今世界现代战争精确打击的战术技术要求。因此,必须按文中提出的,满足N-S方程解的相似性要求去确定模型的尺寸,才是今后风洞实验模型设计的方向。从这一点出发,文中还提出了一些新的概念,并且扩展到如何应用CFD来作为风洞实验相似理论辅助研究的设想。     

塞式喷管侧向喷流仿真研究

随着快速、准确等控制要求的提高,越来越多的飞行器采用了直接侧向力控制技术,这对飞行器本体及喷流装置均提出了更高要求,是先进飞行器气动设计中的关键技术之一。参考美国NCADE方案中的塞式喷管控制发动机结构进行了侧向喷流干扰流场的数值仿真研究,并与常规拉瓦尔喷管进行了对比,分析了两种喷管干扰流场结构及性能的异同,为相关研究提供借鉴。