三维椭圆型方程网格生成中的新源项修正法

构造一种新的源项修正方法用于三维椭圆型方程网格生成，可以实现对网格与边界正交及网格与边界间距的双重直接控制．     

分区对接网格下多段翼型N-S方程计算

对多段翼型的复杂粘性流场，各翼段物面及尾迹区生成C型网格．外场采用O型网格，中间用H型附加区过渡．自动生成高效的分区对接C-O-H型网格。采用时间精确有限体积高精度格式解可压缩N—S方程。层流、紊流的N—S方程求解结果，说明了本方法的有效性。     

用k-ω紊流模型和Osher格式在混合网格上计算紊流流场

用有限体积法求解了雷诺平均 Navier-Stokes方程。空间离散采用了 Osher逆风格式 ,紊流模型采用了两方程模型 ,并对绕 NACA0 0 1 2翼型及双椭球外型的紊流流动进行了数值模拟。结果表明用上述方法求解紊流流动非常有效     

调和微分求积法权系数矩阵的一种显式计算式

简要介绍了调和微分求积法，导出了求一阶导数权系数矩阵的显式计算公式。利用该公式和其中反心对称的性能，可进一步提高计算效率。由于均匀网点有时不能给出可靠的解，本文导出了几种能出可靠结果的不等距网点公式，其中一种公式虽然用不同的方法导出，但结果与Ｇａｕｓｓ－Ｌｏｂａｔｔｏ方法等价，本文还证明了调和微分求积法权系数矩阵具有中心对称或中心反对称的性质（取决于导数的阶数），利用这些性质可以进一步减少计算工程     

用非交错网格求解曲线坐标系下的流动控制方程

本文采用非交错网格法,求解任意曲线坐标系下的流动控制方程。文中阐述了压力振荡的根本原因和相应的消除措施,导出了压力修正方程中控制体边界上曲线速度分量的修正表达式。在这些表达式中,出现了一个表征一个网格间距和两个网格间距压力梯度差分值的附加差值项。当出现压力振荡时,该差值很大,可以有效地消除振荡,而当压力场没有振荡时,这一差值又很小。作为计算方法可靠性的检验,本文分别计算了二维直通道和渐缩形混合管内的气流流动问题,结果是令人满意的。     

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

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

三维动态非结构重叠网格Navier-Stokes方程并行算法

 提出了一种三维动态非结构重叠网格Navier-Stokes(N-S)方程的并行计算方法。N-S方程的空间离散采用格点有限体积方法,时间离散采用隐式的双时间步长方法。应用一方程Spalart-Allmaras(S-A)模型来计算湍流黏性。并行计算采用动态的区域分裂方法,在每一物理时间步利用METIS网格分区系统对网格进行分区。为了实现各CPU之间的负载均衡,每块网格都按CPU个数进行分区并对活动节点和非活动节点进行了加权处理。最后,通过对外挂物投放无黏流动的数值模拟和内埋武器弹舱开启黏性流动的数值模拟,验证了该并行程序的准确性、高性能并行计算以及处理复杂几何外形的能力。     

机群环境中基于信度模型的作业调度算法

提出了一种新的机群环境下的作业调度算法,该算法基于信度模型.在信度模型 中,每一个计算节点是一个信度实体,系统根据信度实体的物理能力、当前状态、历史作业完成情况,对节点进行信度评价,赋予其相应的信度值,调度算法根据信度实体信度值为作业分配合适的资源运行作业.同时,调度系统还对机群系统的整体信度进行评价,并根据估计的作业运行时间和完成期限对作业做信度要求折算,调度系统根据机群系统的整体信度情况和作业信度要求来决定是否接受该作业请求.在OpenPBS机群调度系统上实现了该算法,并通过GridSim工具包与先入先出FIFO(First Input First Output) 算法和Libra调度算法进行了模拟对比试验.试验结果表明,该算法是有效的.      

复杂外形网格生成技术及应用

以某型双立尾战斗机为例,给出了实际复杂型号飞行器结构网格生成的一种分块方案及实现过程,用一种矢性三次多项式插值的方法生成相邻块公共交界面的网格,用Hilgenstock方法修正源项,求解椭圆型方程生成块内空间网格,使网格线与边界面夹角、与边界面间距控制为期望的值.在所生成的双立尾战斗机分块结构网格中求解Euler方程,获得了合理的全机气动力系数、翼面压力分布和绕流图谱结果.     

空间在轨环境下的30 cm离子推力器三栅极组件间距变化仿真分析

采用有限元仿真(FEM)与地面热平衡试验验证相结合的方法，计算并模拟了30 cm离子推力器处于在轨环境时，有、无主动热控对三栅极相对位移变化造成的影响，并对目前离子推力器设置的工作启动流程可能造成的打火风险进行了预估。结果显示：三栅极组件的热形变方向均为法向方向，且栅极中心区域的间距最小；在 -269 ℃ 在轨极限环境温度下，推力器在5 kW工作模式下温度平衡后的屏栅与加速栅最大热态间距为0.14 mm，加速栅和减速栅则已发生贴合；在受太阳辐照以及卫星帆板恒温边界的影响下，栅面最低初始温度为-102 ℃；当推力器主动热控保证温控点为20 ℃时，栅面最低启动温度为-25 ℃，且推力器工作8000 s后，屏栅与加速栅、加速栅与减速栅的最小间距分别稳定在0.25 mm和0.20 mm；当推力器主动热控保证温控点为50 ℃时，推力器工作9000 s后，屏栅与加速栅、加速栅和减速栅最小间距分别稳定在0.31 mm和0.30 mm，能够满足0.25 mm的栅极安全打火间距要求。