鲁棒的结构网格自动化重叠方法

网格技术是目前数值模拟中的关键技术之一。重叠网格是一种放宽拓扑要求、减小网格生成难度的网格技术。本文以结构重叠网格为基础,分别针对挖洞、寻点以及洞面优化方法进行了研究和改进,同时完成物面网格重叠,形成了一套鲁棒的、自动化的网格重叠系统。在挖洞方面,结合“最小洞映射”方法,提出“复合式挖洞”方法,节省内存开销;在寻点方面,通过构建格心虚网格,保证搜索空间的连续性,同时结合“有效搜索”思想,排除部分对寻点无贡献的网格点,进而减少ADT叉树节点;在洞面优化上,改变填补判别法则并引入两类受保护洞内点,确保两层插值边界建立,提高鲁棒性;在物面网格重叠上,利用物面投影法完成坐标修正,实现物面附近网格流动变量的准确传递。为验证本文方法,分别对定常翼身组合体DLR-F6绕流和非定常机翼挂载分离过程进行了数值模拟,计算结果与实验结果吻合良好,表明该结构重叠网格系统对多物体间定常、非定常扰流具有较好的数值模拟能力和较高的模拟精度,具有较高的工程应用价值。     

基于超网格的重叠网格守恒插值方法

保证重叠网格边界数据插值的守恒性是计算流体力学面临的一大难题。基于超网格方法与格心格式有限体积法,发展了一种新型混合重叠网格守恒插值方法。详细研究了用网格切割技术在重叠网格边界构造局部超网格,用网格求交算法合理地扩大贡献单元模板,建立了一种适用于任意网格类型的隐式并行重叠网格守恒插值方法,以超网格为媒介可实现重叠网格二阶精度的守恒插值。数值结果表明,本文方法对二阶分布的流场变量具有严格的守恒性,相比三线性插值方法和逆向距离权插值方法,本文方法减小了数值误差,提高了重叠网格边界的插值精度,加快了计算收敛速度,改善了重叠区域网格尺度相差较大时流场的光滑性和连续性。     

新型重叠网格洞面优化方法及其应用

高度自动化的洞面优化方法是重叠网格关键技术之一。通过对割补法的分析,针对其洞面切割后洞边界位置不确定性及不可预知性问题,提出将其与物面距离优化准则相结合的新型混合洞面优化方法。该方法针对各物体间体网格,采用物面距离优化准则,将洞边界位置控制在物体中间的位置,使重叠区域更合理、可预测;针对背景网格,将其物面距离设置为较大值,采用物面距离优化准则将其落入物体网格内部的单元一并挖去;针对物面处重叠的网格,采用割补法对其进行切割和填补,保证物面及近壁面网格得以有效的重叠。该方法在保证高适用性的前提下提高了网格重叠质量,自动化程度高,无需人工干预。3个典型复杂流动算例计算结果与实验结果吻合良好,网格重叠区域流场变量传递正确,等值线过渡光滑,流场刻画准确,证明该方法准确可靠。     

基于非结构重叠网格的二维N—S方程求解与应用研究

本文发展了非结构重叠网格技术.通过构造任意单元类型网格重叠的网格间边界确定算法、宿主单元搜索算法以及网格间插值算法,把非结构重叠网格推广应用于求解NS方程,并结合湍流模型模拟二维湍流流场.N-S方程的求解采用有限体积法,通量计算采用Osher格式,湍流模型采用SA模型.为验证本文方法的正确性和有效性,模拟了RAE2822翼型和某三段翼型绕流湍流流场,结果表明采用非结构重叠网格模拟复杂外形湍流流场不但网格生成效率高,而且容易推广应用于具有相对运动的非定常问题,具有较高的工程应用价值.     

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

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

基于重叠网格技术和多重网格算法的悬停旋翼粘性绕流数值模拟

发展了一种基于重叠网格的多重网格算法,并应用于悬停旋翼粘性绕流的Navier-Stokes方程数值模拟.多重网格算法在重叠网格上实施主要面临如下两个问题:粗网格间几何关系的建立,即洞单元、边界单元的确定;洞边界存在的情况下如何进行粗细网格间流场信息的传递.本文针对以上问题给出了相应的解决方案.采用两层网格的V循环,应用有限体积法对一跨声速悬停旋翼粘性流场进行了数值计算.计算结果表明该多重网格方法可以显著加快悬停旋翼粘性流场数值解收敛速度,加速比约为3.     

用动态非结构重叠网格法模拟三维多体相对运动绕流

 发展了动态非结构重叠网格技术。本文提出了基于阵面推进的相邻单元搜索算法和适合于高阶空间离散格式的网格间插值关系定义方法，前者提高了非结构重叠网格在模拟多体间具有大幅相对运动的非定常问题时网格处理的效率，而后者保证了网格间边界附近流场求解精度；本文非定常Euler方程的求解采用有限体积法，通量计算采用Osher格式，时间积分采用双时间步长推进。为验证本文提出的方法，分别对定常M6机翼绕流和机翼下外挂物投放过程的非定常绕流进行了模拟,结果表明本文提出的方法提高了非结构重叠网格在模拟多体间具有大幅相对运动问题时的效率，保证了网格边界处的精度，具有较高的工程应用价值。     

旋翼流场计算嵌套网格并行装配方法改进研究

针对嵌套网格系统下进行旋翼流场计算时的网格装配问题,给出了一种改进的并行化网格装配模型。在该模型中,使用压缩存储的辅助结构网格进行并行"洞切割",生成插值边界;应用基于插值点坐标范围的局部可自适应辅助结构网格定位,缩小了处理器之间进行并行贡献单元搜索时的交互数据范围及生成辅助结构网格的内存消耗,结合局部范围的Neighbor to Neighbor(N2N)搜索重启和穷举方法,消除了网格分区、物面、外场等边界面对贡献单元搜索过程的影响,且方法对网格的拓扑没有限制,适合于结构、非结构、直角网格等网格块的装配过程。在流场计算方面,采用并行化的隐式双时间格式数值求解Navier-Stokes方程,用于旋翼数值模拟问题。理论分析和计算测试表明,在使用粗化的辅助网格时,改进的计算方法能够达到传统方法贡献单元搜索速度近似10倍,在低内存开销情况下降低网格装配过程对流场计算时间占比。在内迭代大于7次时,网格装配时间小于流场计算的1%。该方法是对旋翼多周期计算时反复网格装配过程的很好的效率改进。     

带制退器复杂流场并行算法研究

基于集群系统(CLUSTER),在Linux和MPI并行环境下,研究了带制退器的复杂膛口流场三维数值模拟并行算法。计算网格是三维多区结构化网格,数值格式采用Jameson格式及LU隐式分解方法。针对流场的复杂拓扑结构,建立分区信息文件记录各分区的相关信息,并提出了通过合并相邻网格分区的方法均衡各个进程的负载,提高程序的并行效率。在处理需要通讯的信息时,设计了自定义数据类型和"设置中间参数"相结合的方法,并通过动态分配内存减少了程序对计算机系统资源的占用。数值实验表明,本文设计的算法可以适应拓扑结构较复杂的流场,并能取得理想的并行效率,将有助于制退器性能的研究。     

带控制律导弹投放数值模拟

基于非结构重叠网格技术,耦合舵面控制律模块,采用非定常雷诺平均NS方程,建立了模拟控制律作用下导弹投放分离动态特性的CFD数值方法.非结构重叠网格生成技术以“物面相交”方法确定洞边界,在多套网格重叠部分采用算术平均进行物理量的计算传递.采用机翼/挂架/带舵外挂物模型,以及不同的舵面控制规律展开了对比计算,分析了结果.并比较了某导弹有无控制律作用时从某战斗机分离的动态特性,研究了增加控制律的必要性.     

二维DRAGON网格自动生成技术的研究

对重叠网格技术中新出现的Cut-Paste方法进行了研究和改进,结合网格洞边界标识的缝合与重组织技术,首次将其应用在DRAGON网格生成中,提出了二维DRAGON网格的自动生成方法,最后编制出通用程序,给出了多体网格的复杂示例.实际计算表明,本文方法具有自动化程度高、可靠性好、DRAGON网格质量高等特点.     

判断力权值的确定方法及AHP中群组判断矩阵的形成

在ＡＨＰ模型中，请一组专家对某一属性进行评价，根据专家的判断能力，给出判断力权值，再通过适当的数学处理方法，可得最优判断矩阵，这是ＡＨＰ中构造判断矩阵的一种新方法。     

AHP中判断矩阵元素最优估计值的误差分析

在AHP模型中，请一组专家对某一属性进行评价，通过适当的数学处理方法，可得判断矩阵元素的最优估计值，并根据残差来确定所求估计值是否达到精度要求，若不符合精度要求则需返还给专家进行重新评估。这是AHP中构造判断矩阵的一种新方法。     

基于遗传算法的多机器人系统集中协调式路径规划

根据多机器人系统无碰撞运动的需要,对其工作空间进行了分解,确定了机器人运行路线上的各个可能路径点,从而得到了规划空间的多路径点链接图描述。基于这种对规划空间的链接图建模描述,开发了一种混合遗传算法用于寻找多个机器人的无碰撞协调运动路线。仿真结果表明,这种方法可有效地解决复杂规划空间下的多机器人路径规划问题。     

基于聚类分析和投影距离融合的航空发动机

针对类别判定时对类中个体信息挖掘不足的情况,提出1种新的类别判别方法,充分利用了某1类中各成员的信息,构造待测样本的投影距离,并以此作为分类标准.考虑单一分类标准的片面性,将滑油中金属变化量及变化率进行融合.取距离倒数作为融合参数,得到最终分类结果,从而判断航空发动机所处的性能状况.实例表明:对于测量参数存在小幅值随机噪声的情况,也能得到准确的分类结果,为发动机视情维修提供参考依据.     

判断矩阵的最少元素构造法

在层次分析法中构造判断矩阵并使判断矩阵具有一致性是非常关键的。为减少正互反矩阵两两比较的复杂性,降低计算难度,从判断矩阵的构造入手,首先确定相对独立的主对角线上(下)侧元素,而非独立元素就按照一致性原理计算。利用最少的元素构造判断矩阵,使构造矩阵过程简洁,效率提高,满足一致性的要求。     

一种新的判断二维直角坐标网格角点在任意外形物面内外的算法

 根据计算机图形学中的视线遮挡原理,提出了二维平面上判断直角坐标网格角点在任意外形物面内外的一种新算法,并计算了一个 27× 15的直角坐标点阵元素与 3个等腰三角形物面的位置关系,作为该方法的一个算例。计算结果表明,该方法精度高,计算量小,通用性好,并且易于推广至三维复杂情况。     

航空发动机性能趋势监控定量化的方法研究

通过对某型飞机飞行参数记录系统记录信息的破译、分析和处理,对得到的与发动机性能相关的数据在航空发动机性能趋势监控中的应用进行了较为深入的研究。提出了如何利用这些数据得到综合参数,并采用综合参数作为评判发动机性能状况的量化指标以及在此基础上建立发动机预警信号系统的实现方法。最后,利用实际的发动机数据进行了验证。      

计算机辅助FGF判断法解直线度误差

介绍一种按FGF判断法求符合最小条件直线度误差的快速算法，该算法很容易在计算机上实现，具有简单、明了的特点。     

图像分析法测量液滴接触角

接触角是表征液体对固体润湿程度的重要参数,本文采用基于Matlab计算软件的图像分析法测量液滴接触角。首先使用图像增强、去噪、锐化手法对CCD采集到的液滴图像进行预处理;再运用Canny多阶段优化算子提取液滴轮廓线;最后以多项式为拟合框架,获得了接触角稳定状态下液滴轮廓区间和接触角测量结果。本文为接触角测量技术提供了一种有效实用的方法。