自吸气脉冲爆震发动机超声速进气道气动性能的数值研究

采用有限体积法计算了某轴对称超声速进气道,在不同幅度和频率的出口扰动压强下的进气道内结尾正激波的运动情况,得出了进气道内结尾正激波运动特性和扰动压强的频率和幅度的关系。在计算中,本文采用了多块结构化网格,控制体积的界面无粘通量采用三阶迎风格式插值获得,同时采用了minmod通量限制器,以确保在激波处的解的物理特性;扩散通量采用二阶中心差分格式插值获得。定常计算采用当地时间步法,非定常计算采用双时间步法。离散的代数方程采用交替方向迭代法求解。     

自适应直角切割网格民机增升装置绕流数值模拟

提出并介绍了基于直角切割网格的分区面搭接技术,采用变长宽比网格方法,成功地进行了存在外形间断(剪刀叉)的民机增升装置的网格生成和绕流Euler方程数值模拟.根据外形的特点,运用"根"网格分区算法,降低了整个网格的生成难度;通过基于外形的自适应网格加密技术,详细地描述了外形上的大量缝道和凹槽,提高了网格质量;在分区交界面上,利用重叠面积切割和面搭接算法实现了两侧网格间的流场信息传递,并保证了通量守恒;采用中心有限体积方法,结合双时间推进算法,完成流场的Euler方程数值模拟,计算结果与实验数据吻合良好,说明所述方法的正确性.     

三维直角叉树切割网格Euler方程自适应算法

基于叉树数据结构，实现了一种用于三维直角叉树切割网格的自适应算法，包括对几何外,形以及对流场计算的自适应网格加密技术。在初始网格的生成过程中，根据相邻网格的物面法向向量间的差值，进行针对外形的自适应网格加密；在流场计算中，根据相邻网格间选定物理量梯度的变化，进行针对流场的自适应网格加密。详细地描述了三维直角叉树切割网格的生成过程，以及以任意网格的切割细分算法。在自适应过程中，分别采用了八叉树和全叉树的数据结构，八叉树是基本的数据结构，而全叉树的采用，使网格具有了各向异性的特征，从而大大的减少了自适应网格的数量。采用中心有限体积法，求解Euler方程，并运用上术方法，完成了对外形和流场的逢适应网格加密算法，获得了较好的数值计算结果，证明了自适应算法的正确性，体现了直角叉树切割网格自适应技术的有效性和实用性。     

基于雷达散射截面约束的气动优化设计

采用拉格朗日乘子法优化设计了雷达散射截面约束条件下的锥形融合气动外形.拉格朗日乘子法中的极小化问题采用动态演化的优化设计方法求解.该方法是一种基于非定常演化的优化设计方法,即在求解非定常流动支配方程的时候同时履行优化过程,较其它基于定常解的优化方法具有高得多的计算效率.其中的雷达散射截面通过求解非结构的笛卡儿网格上的时域麦克斯韦方程来得到,而升阻比则通过求解锥形流方程来计算.通过优化设计,得到了M∞=8.0时,升阻比为4.98,雷达散射截面只有1.66m2的锥形融合气动外形.     

带壁面射流燃烧室燃烧数值模拟

利用基于非结构化网格有限体积法对三维有壁面射流的燃烧室内两相流动和燃烧进行了数值研究.对气相流动在Euler坐标系下求解,而对液滴相则利用Lagrange方法进行追踪求解.计算区域采用四面体网格进行划分,气相流场用SIMPLEC计算方法,对液滴相采用了欧拉隐式方法.考虑了液滴相与气相的完全双向耦合作用,分别采用了Spalding液滴蒸发模型和涡破碎(EBU)燃烧模型,数值计算结果与文献中实验数据吻合较好.      

几种典型迎风格式的分析与比较

重点分析了AUSM 格式的构造方法,同时将其与Van Leer, Roe格式等构造方法进行比较和分析,通过求解一维激波管问题,以及二维和三维跨音速流场的计算比较了AUSM 格式,VanLeer格式,Roe格式和中心格式的数值模拟结果和网格适应性,以及在经典四步Runge-Kutta显式时间推进和LU-SGS隐式时间推进下上述四种格式的计算效率.获取了对于基础研究和工程应用很有意义的结论.     

三角翼大迎角绕流特性数值模拟的网格处理技术研究

研究网格对三角翼大迎角绕流特性数值模拟的影响.结果表明,网格的生成需要结合流动现象,对网格拓扑结构和网格点分布进行选择与搭配.C-H型网格适宜模拟尖前缘分离涡流态,法向网格在一定范围内应等距增长,沿流向逆压梯度较大的区域内适当增加网格点,尾迹区网格则应做上翘处理.Euler方程具有模拟三角翼旋涡及预测涡破裂特性的能力,但对二次涡等粘性引起的流动细节把握能力不足.利用层流假设的N-S方程,通过合适的网格,也可得到满意的计算结果,但对涡破裂后的强烈非定常湍流流动模拟能力不足.采用旋涡螺旋度可准确反映主涡与二次涡流动,描述旋涡的破裂现象.用轴向速度迅速减小并小于来流速度的点作为涡破裂判据似应更合理.     

新型碳纤维格栅夹层结构模拟研究

新型飞行器机体构件的设计要求在所允许的有限重量下实现结构超轻型化。格栅夹层结构是当前国际上最有前景的新一代材料结构之一,已开始应用于航天结构。为分析三角形点阵格栅结构的受力性能,本文通过解析分析,用连续体模型模拟平板状格栅结构,得到各种不同尺寸下格栅结构的解析通解。经与ABAQUS有限元软件模拟试验对比,这种方法可将误差控制在较小范围内。与在相同重量条件下传统的蜂窝夹层结构相比较,进行三点弯曲模拟试验,三角形格栅结构具有更高的抗弯性能。     

网格环境下基于移动代理的两级环调度模型

网格环境下基于移动代理的两级环调度模型主要由网格任务管理中心、网格资源节点、网格用户、用户作业组成。网格任务管理中心依据网格用户提交的作业动态生成两级环,即超级环和普通环。两级环中资源节点由超级环管理节点、普通环管理节点和普通环资源节点构成。不同节点充分利用移动代理的自主迁移能力共同维护任务的正常运行。通过对两级环模型的事件响应时间延迟、系统可靠性、负载平衡、作业运行时间的分析和仿真,证实该模型能够有效利用网格计算资源、均衡资源负载。     

悬停共轴双旋翼干扰流动数值模拟

基于非结构动态嵌套网格方法,通过求解三维非定常Euler方程对悬停共轴双旋翼直升机的复杂流场进行了数值模拟研究.双旋翼之间的干扰使得两者的拉力性能下降,上旋翼对下旋翼的干扰程度远大于下旋翼对上旋翼的影响程度.下旋翼处的下洗速度比单旋翼的要大得多,而上旋翼的下洗速度比单旋翼的略大一些,几乎相同.桨叶总星巨角增加,则下洗速...