基于非结构动态网格的翼型积冰过程数值模拟

讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法。外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格 库塔方法。采用四阶龙格 库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率。使用非结构动态网格技术对积冰后的翼型网格进行快速调整。在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性。     

基于非结构网格的非定常流数值模拟方法研究

本文给出了两种基于非结构网格的非定常流动问题的数值模拟方法.一种是基于非结构网格的重叠网格方法,该方法结合了非结构网格方法和重叠网格法的优点,适合于处理复杂外形飞行器的非定常运动问题.另一种是基于Denauley图映射的动弹网格方法,该方法具有较强处理粘性网格的变形运动的能力.本文通过多个数值算例(三维外挂物投放、直升机机身+旋翼前飞、后缘襟翼大偏角运动以及机翼弹性变形等)对这些方法的正确性和有效性进行了验证,结果表明这两类非结构网格的非定常流数值模拟方法是精确高效的.     

聚合多重网格法在旋翼前飞流场计算中的应用

基于非结构动态嵌套网格, 将聚合多重网格方法引入到旋翼的前飞流场计算中.从初始的整体非结构网格出发, 快速建立相互重叠的两个子域网格, 同时无需建立背景网格便可以确定初始洞边界和初始插值贡献单元, 且动态嵌套时无需更新洞边界, 并采用了一种无需背景网格便可快速查找新的贡献单元的方法.分区利用聚合法生成多重网格, 利用双时间方法分区运用三层网格V循环求解非定常Euler方程, 两区的流场参数在重叠区实现耦合, 对旋翼的前飞流场进行了计算.结果表明, 多重网格方法能够有效地加快旋翼前飞流场收敛.      

一种预测平板尾迹噪声的时域无网格方法

基于离散涡方法和涡声理论建立了一种预测二维平板尾迹发声的时域无网格方法。该方法应用解耦方式完成声场计算,首先使用离散涡方法计算了均匀来流作用下的平板尾迹流场,得到了流场中点涡的涡量、位置和速度等关键参数,然后基于涡声理论建立了自由空间中点涡发声模型,并引入了时域边界元方法来模拟平板表面对声场的散射作用,计算得到了平板尾迹涡发声的偶极子声场分布和指向性等关键特征。通过对上下表面涡以及平板散射对声场贡献的深入分析表明,从平板尾缘上下角点脱落并卷起的涡团均为偶极子源,平板的散射作用使得声场在一定程度上得到加强,并且使声场具有极大值方向垂直于平板表面的偶极子指向性特征。所建立的无网格方法计算快速,能同时获得流场和声场分布的关键特征,可提升对气动噪声产生机理的基本认识,同时还为尾迹噪声的理论研究提供了一种具有工程应用价值的可靠计算方法。     

运输机空投舱内压力动态变化数值模拟研究

运输机重装空投时,机舱横截面方向堵塞度较大,快速移动的空投物会造成舱内压力剧烈变化,甚至损害舱内乘员的身体健康,因此,有必要对运输机空投舱内压力的变化展开数值模拟研究。以运输机机舱内大型货箱在牵引伞拉力、重力、摩擦力以及气动力共同作用下 x 方向的变加速直线运动这一流固耦合问题为研究对象,基于运动嵌套重叠网格,对运输机重装空投前舱动态压力变化进行数值模拟,分析模型尺寸、来流速度、机舱堵塞度、投放位置等参数对舱内动态压力的影响规律,进而提出运输机重装空投风洞试验时与舱内压力相关的相似准则。结果表明：运输机执行重装空投任务时,空投物距离前舱位置越近,来流速度越大,堵塞度越大,空投时前舱动态压力变化越大,反之越小;风洞试验时,用空投前前舱压力值为参考量的舱内无量纲数是最准确、合适的舱内压力相似准则。     

球坐标系六片网格下三维定态行星际太阳风模拟

采用二阶MacCormack差分格式, 利用稳态的磁流体(MHD)方程组在球坐标系六片网格下模拟研究了行星际太阳风. 六片网格系统能有效避免极区奇性和网格收敛性. 迭代按径向方向推进求解, 很大程度上减少了计算量, 节约了计算时间. 内边界条件根据太阳与行星际观测确定, 比较测试了5种内边界条件, 模拟给出了1922卡林顿周的背景太阳风结构. 几种内边界条件所得模拟结果与行星际观测基本吻合. 太阳风速度采用McGregor 等的经验公式给出, 磁场由水平电流片(HCCS)模型得到, 密度和温度分别根据动量守恒和气压守恒得到, 研究表明采用这样的边界条件模拟结果最佳.      

围绕旋翼飞行器的三维结构化运动嵌套网格生成方法

考虑到旋翼飞行器的旋翼各种运动和机身复杂外形特点,建立了一种鲁棒性高的三维结构化嵌套网格生成方法。采用三维Poisson方程对初始网格进行迭代,获得了高质量的旋翼桨叶三维结构化网格。以此方法为基础,以复杂外形几何体表面网格为初始网格,建立了一种抛物型法向外推方法。推进过程中,使用经抛物化处理的新Poisson方程进行光顺迭代,源项由Hilgenstock法确定,有效克服了椭圆型方程生成网格时边界点较难调整的困难。通过网格最小正交度的检测,表明了该抛物型方法的先进性。进一步将高效的挖洞法("Top map"法)和贡献单元搜寻法("Inverse map"法)引入上述三维结构化网格生成中,建立了围绕旋翼飞行器的运动嵌套网格生成方法,采用RANS方程分别对Caradonna-Tung旋翼和V-22倾转旋翼机流场进行计算流体力学(Computational fluid dynamics,CFD)计算,表明该网格生成方法能够满足CFD计算的需求。     

气动弹性计算中网格变形方法研究进展

网格变形是气动弹性计算中实现计算网格随运动边界变形的主要方法。在总结目前网格变形方法发展现状的基础上,对近几年常用的网格变形方法,即弹簧法、弹性体法、超限插值法、Delaunay背景网格法、径向基函数插值法和温度体法等做了简要的总结。根据各方法构建模型的不同,将它们分成物理模型法、数学插值法和混合方法3类,简要介绍了各方法的基本思想和研究进展,重点比较了各方法的网格变形特性（变形能力、变形质量和变形效率）和优缺点。总结了当前流场边界与结构边界的数据传递方法,对当前气动弹性计算中遇到的网格变形的难点问题作了简要评述并对未来网格变形方法的发展方向进行了探讨。     

一种改进的ORB特征匹配算法

针对现有的ORB特征匹配算法在图像模糊、光照变化、图像压缩、噪声条件下,匹配准确率下降问题,提出了一种改进的ORB特征匹配算法。首先,在提取特征点过程中,对图像进行网格化处理,并引入四叉树结构,使提取的特征点在图像中均匀分布,解决传统的特征提取方法遇到的特征点集中问题。然后,利用暴力匹配进行初步匹配,并采用交叉验证的方式,剔除部分误匹配,改善暴力匹配的结果。最后,利用高斯核对网格运动统计的结果做加权处理,优化统计结果,进一步剔除误匹配,得到准确率更高的匹配集合。实验结果表明:改进后的算法在图像模糊、光照变化、图像压缩和噪声条件下,平均准确率分别提高了3.5%、4.2%、2.2%和6%。      

分区多重网格法在复杂区域中的数值研究

突破多重网格方法在多连通域中应用的局限性,将分区法与多重网格法结合,形成分区多重网格法,应用在复杂几何域及复杂多连通域中,结果表明,它不仅能提高精度,而且能加快收敛。