二维结构/非结构混合网格的生成及其流场模拟

采用分区求解方法在二维结构/非结构混合网格中对粘性绕流作了数值模拟。利用层面推进法生成绕物体的O形结构网格,利用推进阵面法在结构网格外生成非结构网格。在结构网格区内,采用B-L湍流模型,求解雷诺平均N-S方程。在非结构网格区内求解Euler方程。两区的解通过在交界面上满足通量守恒条件进行耦合。数值算例表明了方法的有效性和实用性。     

飞机外形网格拓扑和多块网格拼接技术

多块网格方法成功的关键在于网格拓扑的选择和块与决之间的拼接方法。针对真实飞机外形多块网格生成,本文分析研究了翼类和体类部件的各种网格拓扑;论述了翼身组合体的多块网格非结构连续拼接技术;指出连续拼接方法的缺点,以及相应的解决方法;同时,还讨论了分块策略和确保块与块之间网格线完全连续拼接的方法;通过实例进一步阐述了本文工作的另外一些独到之处。     

基于非结构重叠网格的螺旋桨滑流非定常数值模拟

基于各向异性非结构混合网格及多套网格重叠技术,通过求解非定常雷诺平均Navier-Stokes(URANS)方程,分别研究了螺旋桨动力特性和螺旋桨滑流对涡桨飞机气动特性影响。运用物面相交准则实现多套网格间的挖洞,采用距离权和三线性插值技术传递重叠区变量信息,并针对定轴旋转优化了网格重叠边界。采用双时间步求解控制方程,内迭代计算采用LU-SGS(lower-upper symmetric Gauss-Seidel)方法。模拟了单独螺旋桨旋转的非定常流动,计算的拉力系数和扭矩系数与试验结果一致,表明采用的数值方法和网格技术能有效模拟螺旋桨滑流效应。数值模拟了某涡桨飞机,对比了有无滑流下飞机的气动特性,给出了滑流对飞机影响。结果表明:滑流导致其后部气流加速和旋转,改变飞机气动特性,增大飞机的阻力和俯仰力矩,导致气流下洗,影响机翼及平尾压力分布。      

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

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

一种新的基于位置信息的路由算法研究

 为航空移动Ad hoc网络(MANET)提出一种基于位置信息的路由算法,即位置网格路由（LBGR）算法,以此来解决节点高速移动引起的路径重建问题。算法中,数据包沿路由发现时得到的网格轨迹进行传输,每个收到数据包的中转节点根据它与相邻节点、目的节点的位置关系及网络轨迹下游的路由网格位置决定下一跳节点。由于利用了网格位置信息,LBGR算法在航空节点高速移动和拓扑变化无法事先预测的情况下,使用较少的开销就可获得稳定的路由。仿真结果表明,与GRID路由算法相比,LBGR算法在节点中高速移动时,能够有效地减少路由开销,提高网络吞吐量。     

大规模混合网格的分区策略研究

在飞行器阻力数值模拟中常常需要数百万网格点的巨型混合网格,甚至需要网格点数达到数千万的超巨型混合网格.在并行计算前往往需要对这些大规模混合网格进行分区操作.针对Metis库函数所需要输入的大规模边表的生成问题,提出了一种生成边表的并行算法,并运用此算法,在微机机群上对具有4,787,893个网格点的DLR-F6翼型(翼身组合体 挂架 发动机短舱)的巨型混合网格成功地生成21,556,110条边的大规模边表,进而采用Metis提供的多级循环二分法顺利地对该巨型混合网格实现均匀分区,从而验证了该并行算法的有效性.     

PARAMETERIZATION OF TRIANGULAR MESHES

A new algorithm called the weighted least square discrete parameterization (WLSDP) is presented for the parameterization of triangular meshes over a convex planar region. This atgorithm is the linear combination of the discrete Conformal mapping (DCM) and the discrete Authalic mapping (DAM). It provides the good properties of both DCM and DAM, such as robustness and low distortion. By adjusting the scaling factor q embedded in the WLSDP, satisfactory parameterizations in different special applications can be achieved.     

对三向测力仪性能评价的探讨

三向干扰系数通常作为评价三向测力仪静态性能的一个重要参数。分析与实验证明这种评价方法不能真实反映测力仪在使用中的干扰情况,因而提出用平面干扰系数代替三向干扰系数来评价测力仪性能,并对测力仪在设计、制造、使用中的一些问题提出了看法。     

应用于DSMC方法的二维Voronoi网格技术

研究了Voronoi网格技术并将其应用于直接模拟蒙特卡洛DSMC(Direct Simulation Monte Carlo)计算.基于Dirichlet镶嵌与Voronoi图理论,Voronoi网格利用特征点表征网格单元,具有建立粒子与网格单元之间映射关系的独特算法,适合于DSMC方法的统计特点.在剔除过于靠近边界的特征点以及必要情况下边界细化的基础上,通过区分由边界节点表征的非完整Voronoi网格单元以及由计算区域内镶嵌点表征的完整Voronoi网格单元,解决了Voronoi网格的二维边界匹配问题.Voronoi网格技术支持自适应DSMC计算.映射效率对比表明,Voronoi网格的DSMC计算效率高于三角形网格,低于多级直角网格.通过MEMS微喷管流动数值模拟,验证了Voronoi网格技术在DSMC方法中的有效性.      

迎风格式在二维非结构网格中的应用

将原来在结构网格中提出的AUSM⁺(Advection Upstream Splitting Method)格式,推广应用于二维非结构网格中.利用格林-高斯公式对控制体内的变量进行线性重构,获得空间二阶精度,并采用Barth型限制器以抑制数值解的振荡.为提高计算效率,采用了自适应多重网格法.最后给出了NACA0012翼型和平板激波反射的几个Euler方程的算例,并进行了分析比较.