基于时间谱方法的振荡翼型和机翼非定常黏性绕流数值模拟

 传统的双时间方法在非定常计算中长时间的过渡迭代推进求解导致其计算效率相对较低,针对周期性非定常流动问题的流动特征,发展了一种基于离散傅里叶变换的高效时间谱方法,用于求解振荡翼型和机翼的非定常黏性绕流。在时空耦合的雷诺平均Navier-Stokes(RANS)方程的求解中,对流项的离散应用了Roe的通量差分格式,物理时间项的离散方法为时间谱方法,伪时间推进采用了隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式。考虑到湍流的时空耦合效应,时空耦合的Spalart-Allmaras一方程湍流模型的物理时间项同样采用时间谱方法进行离散。为了进一步提高计算效率,当地时间步长和多重网格技术等加速收敛的措施均被采用。算例对俯仰振荡NACA0012翼型和Lann机翼的周期性非定常流场进行了数值计算。结果表明:对于周期性非定常流场的数值模拟,相比于传统的双时间方法,用时间谱方法近似物理时间项,不仅能够提高流场的计算精度,而且更能够大幅度提高计算效率。     

剪刀式尾桨悬停状态气动力及噪声特性计算研究

建立了一个基于计算流体力学(CFD) /FW-H(Ffowcs Williams-Hawkings)方程的预测剪刀式尾桨悬停状态气动性能和噪声特性的分析方法.该方法首先采用CFD方法对尾桨流场进行求解,并应用嵌套网格技术对流场空间进行离散.控制方程采用非惯性坐标系下的Navier-Stokes方程,空间方向采用二阶迎风格式(Roe格式)进行求解,时间方向采用隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式进行推进.在此基础上,采用FW-H方程将尾桨噪声声压扰动传播至远场,以获得尾桨的噪声特性.应用该方法对两种剪刀式尾桨构型(“L”构型和“U”构型)进行了计算研究,对比分析了剪刀式尾桨在气动力和噪声方面与常规尾桨的差别,以及两个重要构型参数(剪刀角和轴向间距)对剪刀式尾桨气动力和噪声的影响规律.计算结果表明,构型参数对剪刀式尾桨气动力和噪声特性影响很大,合理地选择构型参数可以降低尾桨噪声水平.     

跨声速压气机转子叶尖非定常流场数值研究

采用时间精确求解方法对某高负荷跨声速轴流压气机转子在98%设计转速下的叶尖非定常流场进行了数值研究.结果显示,激波自身振荡不明显,叶尖区域流动的非定常性主要来源于叶尖泄漏涡的破碎及其与激波之间的相互干涉.对比设计状态与失速状态下叶尖泄漏涡的特点发现:在近失速点时,叶尖区域间歇性出现前缘溢流.分析表明:叶尖泄漏涡在激波后破碎是造成堵塞的主要原因,也是造成spike型失速初始扰动的原因.      

微型扑翼飞行器非定常运动对平尾的影响

以西北工业大学自行研制的微型扑翼飞行器ASN211为研究对象,利用其简化的二维扑翼及平尾串列翼模型进行了非定常数值模拟,分析了扑翼俯仰运动及沉浮运动对平尾气动性能的影响。在数值模拟模块中,模型的俯仰运动及沉浮运动由动网格技术实现。通过计算流体力学(CFD)软件Fluent对此非定常流场进行数值计算,重点研究了扑翼非定常运动尾流对平尾气动效率的影响。定常状态与非定常时均条件下平尾升力曲线的对比分析表明,扑翼的非定常运动能够增大平尾的失速迎角及最大升力系数,因而使平尾的失速特性得到改善。     

固体推进剂/衬层界面脱粘裂纹的三区域界面层模型

将固体推进剂/衬层粘结界面层按照剪切模量的变化分成三部分,由此建立了其界面脱粘裂纹的三区域界面层模型。应用Fourier积分变换和传递矩阵方法,推导并数值求解Cauchy奇异积分方程,得到了反平面剪切状态下裂纹问题的半解析解,并讨论了界面层模量和厚度对断裂参量的影响。结果表明,界面层模量随贮存时间的延长而降低时,应力强度因子逐渐减小,但裂纹张开位移变化不大;界面层厚度对应力强度因子的影响则不明显。     

三维打印结合化学气相渗透制备Si3N4-SiC复相陶瓷

采用三维打印(3DP)技术成型Si3N4多孔陶瓷并结合化学气相渗透(CVI)SiC制备了Si3N4-SiC复相陶瓷。研究了烧结工艺对3DPSi3N4陶瓷线收缩率和孔隙率的影响。结果表明,3DPSi3N4坯体经热解除碳后再烧结,可以获得较小的线收缩率(<6%)及较大的气孔率(77.5%)。对其进行CVISiC近尺寸强化,研究了Si3N4-SiC复相陶瓷的抗弯强度随SiC体积分数的变化规律。     

用转捩模型预测转捩及确定最佳粗糙带高度

采用SST k-ω二方程湍流模型,附加γ-Reθ转捩模型对不同迎角下的NACA0012翼型绕流进行了转捩预测,计算预测的转捩位置与实验结果符合较好,说明γ-Reθ转捩模型对转捩位置具有较好的预测能力。进一步用模型对带有粗糙带的翼型进行转捩预测,获得了不同雷诺数下刚好激发转捩又不引起过大局部压力扰动和附加阻力的最低粗糙带高度,对固定转捩风洞实验有指导意义。     

某中空长航时无人机机载浮标分离安全性研究

研究机载浮标投放过程中气动干扰特性及浮标运动规律,对机载浮标安全投放方案的设计具有重要意义。采用计算流体力学方法对机载浮标投放过程进行数值模拟。通过动态嵌套网格技术模拟浮标运动,结合刚体六自由度运动方程求解三维非定常 N-S ( Navier-Stokes )方程,模拟基于察/打一体无人机平台的浮标自由投放过程,获得详细的包括浮标运动姿态、运动轨迹和受力情况等在内的数据信息;通过算例,验证本文数值模拟程序的准确性。结果表明:浮标投放过程中,轻质浮标口盖自由分离会对螺旋桨后推式无人机的安全性造成重大影响;考虑浮标口盖固连在浮标降落伞上的状态下,初速度为 0m / s 时投放浮标整体(包含浮标、降落伞、口盖)离螺旋桨最近距离为 0.151m ,安全余量不够;初速度为 1.8m / s 时投放浮标整体离螺旋桨最近距离为 0.671m ,满足投放安全性要求。     

亚格子模型对泵喷推进器宽带非定常力预测的影响研究

为了研究亚格子模型对泵喷推进器非定常流动与宽带非定常力预报结果的影响,采用分块结构化网格建立了模型尺度下艇后泵喷推进器的计算模型,并进行了大涡模拟数值计算。从艇尾非定常流场特征量和推进器转子脉动载荷两个方面对比了三种不同亚格子模型计算结果的差异,并分析了泵喷内部流动与转子非定常力间的内在联系。研究结果表明：三种亚格子模型均能得到含有叶频宽带峰的轴向推力谱,且整体趋势相近。泵喷转子上游的湍流强度和尺度的分布对亚格子模型较敏感,其中Smargrinsky-Lilly模型得到的湍流强度较强,尺度较大,该模型下的湍流谱在非平衡区量级较大,但由低频向高频的衰减较快,并且预测到的分离区范围大,导叶尾缘脱落涡分散,导致转子上游来流空间分布不均程度较强。对于转子叶片上的载荷脉动,Smargrinsky-Lilly模型预测的推力谱中线谱成分明显,并且叶频处宽带谱峰“陡峭”,而WALE模型和KET模型的结果宽带特性较强,对比非定常推力测试结果可知,WALE模型和KET模型更适于该问题宽带非定常力预测。     

轨道转移计算中的逐次消元迭代法

本文针对导弹飞行轨道转移的一个算例,叙述应用消元法思想构造数值求解非线性方程组的迭代算法的一般原理,具体算法步骤和一些计算结果。给出算法成功的理论依据,最后给出实现算法时的一些变形。