低雷诺数下微型飞行器主动变形机翼非定常气动特性数值分析

以微型飞行器的气动外形设计为研究背景,通过数值求解N-S方程分析了主动变形柔性机翼的气动特性。空间离散格式采用中心有限体积法,时间推进采用双时间推进方法,其中子迭代过程由隐式LU-SSOR方法完成。在动网格技术的基础上,模拟分析了薄翼面上作行波运动的流场。计算结果发现,该行波翼型的气动特性相比于刚性翼型有明显改善,有效地抑制了大攻角下大尺度流动分离,同时升阻比提高了约35%,起到了显著的增升减阻的作用。本文的工作对于进一步的柔性机翼优化具有良好的理论研究与数值模拟的参考价值。     

动能BGK算法在近连续流模拟中的应用

首先比较详细地介绍了动能BGK算法的基本思想,然后应用该算法,数值模拟了克努森数(Knudsen number)为0.0001～1.56,马赫数为1.96绕圆柱的二维粘性流动,研究了滑移边界条件和无滑移边界条件对圆柱阻力的影响,通过与实验结果[5]进行比较,验证了文献[1]关于由连续流到自由分子流的这种流动模型的划分.计算结果表明:采用滑移边界条件,直到克努森数为0.3,计算得到的阻力与实验都符合得很好.因此,采用BGK算法加滑移边界条件,能够模拟由连续流到近连续区的流动.     

虚拟体法数值模拟双圆柱强迫振荡系统

用虚拟体方法研究了均匀流及尾流中圆柱横向振荡的流动特征.单圆柱结果与相关试验吻合,得到了"锁定"区域,发现涡量跳跃现象.对小间距双圆柱的数值模拟发现,在"锁定"区域附近随激励频率的增加流场特征会出现若干显著的变化,包括涡量的跳跃,涡模态的改变,正负涡量带的形成和消失.在大间距情况下发现下游圆柱的振荡对上游尾涡脱落没有明显的影响,下游圆柱尾流场没有发现"锁定"现象.     

适用于激波/边界层相互作用的线性涡粘性湍流模式

本文选择了两个激波 /边界层相互作用诱导分离的跨声速问题 (轴对称圆弧突起和二维管道突起 ) ,采用五个有代表性的湍流模式 (BL模式 ,JL k-ε模式 ,k-ω模式 ,SST模式和双尺度模式 ) ,通过将数值计算结果和实验结果进行比较 ,对有关的湍流模式进行了评估和分析     

飞机机翼表面霜状冰结冰过程的数值模拟

提出了一种处理结冰后结冰表面固壁区域移动的移动边界技术,结合欧拉坐标系下空气-过冷水滴两相流动控制方程的计算,对霜状冰的结冰过程进行了数值模拟,得到了NACA0012机翼在0°和4°攻角下机翼表面结冰后形成的冰形,与文献中的实验数据对比,表明本文的方法是可行和有效的。      

螺旋桨等离子体流动控制的增效实验

实验依托搭建的螺旋桨等离子体流动控制测试平台,基于等离子体附壁射流抑制边界层分离的两种机制,采用等离子体射流与来流方向相同的正向射流方式和与来流方向相反的逆向射流方式,研究了微秒脉冲等离子体射流对螺旋桨三维流动分离的控制效果,对比分析了两种射流方式增效特点.实验结果表明:在螺旋桨转速为300r/min,电压峰值为8.5kV,脉冲频率为10~160Hz范围内,正向射流有利于减小螺旋桨转矩,逆向射流对转矩的效果则相反.两种射流方式均提高了螺旋桨拉力和效率,同时其控制效果受脉冲频率的影响较大;正向和逆向两种射流分别使螺旋桨效率最大提高了11.56%,2.79%.      

平衡空气模型的流动稳定性及转捩预测

研究高超声速平板边界层考虑真实气体效应的流动稳定性问题.采用7组元化学反应平衡模型,黏度和导热系数采用混合律,同时考虑组元浓度扩散引起的能量传递,在马赫数为10~20、壁面温度为500~3500K、飞行高度为20~30.5km等条件下,对平板边界层流动的稳定性进行了分析,给出了扰动演化相对增长的N值.计算结果表明:高马赫数飞行中不稳定扰动的第3模态将与第2模态合并,共同影响转捩;高温真实气体的流动稳定性特征,随着马赫数、壁面温度、飞行高度变化的基本趋势与完全气体的基本一致;与完全气体相比,真实气体的相对增长N值包络线较小,表明高温真实气体将抑制转捩发生.      

S弯进气道内附面层抽吸控制对风扇级性能的影响

为了探讨某S弯进气道出口畸变在其后风扇级内的演变过程以及对风扇级性能的影响,在单独对该进气道进行抽吸控制数值研究并优选出最佳方案后,将最佳方案应用到进气道加风扇级的全流道,进而展开全流道数值研究,着重探讨了抽吸控制前后风扇级性能及内部流场结构变化.结果表明:吸气后风扇级整体性能有较大幅度提升,堵塞流量及最大效率分别增加约0.63%和0.57%;进气道出口低能流体在接触转子之前始终聚集在沿程截面底部,所占区域面积沿流向逐渐减小;接触转子至动叶前缘区间内,畸变流体沿动叶旋转方向的相反方向发生偏移,最终覆盖约3个流道;静叶中畸变流体所处流道内的静叶吸力面发生严重的流动分离,且分离主要发生在50%叶高以下,吸气后略有减弱.      

悬臂边界下纤维增强复合薄板固有频率计算及验证

采用理论与实际相结合的方式,对悬臂状态下纤维增强复合薄板的固有频率进行了计算及验证.首先,针对纤维增强复合薄板的结构特点,考虑了纤维方向的影响,对其进行了理论建模.然后,基于正交多项式法来表示振型函数,并通过Ritz法对该类型复合薄板的固有频率进行求解.最后,搭建了该类型复合薄板结构的固有特性测试系统,并以TC500碳纤维/树脂基复合薄板为例,对其固有频率进行了测试.验证结果表明:基于正交多项式法的纤维增强复合薄板固有频率计算结果与实验结果的相对误差在2.2%~9.7%之间,进而验证了所提出的固有频率计算方法的正确性.      

SA和SST湍流模型对高超声速边界层强制转捩的适应性

凸起物是高超声速流动中常用的一种人工转捩装置。采用高阶精度算法模拟了高超声速进气道压缩面上的强制转捩流动,转捩装置为一排高度为1mm的钻石型凸起物和斜坡型凸起物。考察了Spalart-Allmaras(SA)模型和剪切应力输运(SST)湍流模型对该问题的适应性。在考察过程中通过丰富的算例分析了网格规模、可压缩修正和空间离散格式等对计算结果的影响。在层流区,计算能与试验取得非常一致的结果。但在湍流区,计算得到的热流通常高出试验数据。经分析发现其原因是强制转捩的湍流边界层与自然转捩的湍流边界层在涡结构上存在较大差别,使得湍流模型的效果较差。针对高超强制转捩湍流涡结构丰富的特点,对SST湍流模型进行了修改。计算结果表明,该修改方法对提高热流精度具有一定效果。