三维Euler方程组隐式LU分解的高性能并行计算

许多非定常无粘流体力学问题的数值模拟都需要利用Euler方程组来进行计算,而由于在隐格式下,所选取的时间步长可以比在显格式下时大得多,所以隐格式越来越受到重视,其中隐式LU分解是最常用的方法之一。对三维Euler方程组,采用隐式LU分解进行计算时,网格点所在的各个对角阵面之间存在数据依赖关系,本文分析了采用区域分解且边界上用显格式代替隐格式进行计算的高效性,在长方体建筑物内的爆炸模拟表明,在有112个CPU的某MPP巨型机上,并行计算效率超过60%。本文还分析了计算结果与串行计算时的差异,以及利用区域重叠减小这种差异的方法,同时考虑了对处理器进行合理的逻辑组织,将计算网格映射到处理器网格,以最大限度减少通信开销的方法。文中最后以一个爆炸毁伤的例子实际说明了所述方法的可行性与高效性。     

一种适用于压气机非定常流动的相位滞后算法

应用相位滞后边界条件,以1∶1的转、静叶计算域,在不改变转子尾迹及静子势干扰的频率的条件下,计算叶轮机非定常流场.以NASAstage₆₇为算例,计算结果与实验数据相吻合,表明采用的方法具有较大的工程实用性.      

多段翼型绕流的数值模拟研究

采用C—H型混合的结构化网格,选用Wilcox提出的k-ω模型求解可压的N-S方程,模拟了多段翼型的流场。结果表明：在4°和6°迎角下的非定常计算结果与试验结果的对比,计算结果与试验值吻合的较好。本文应用了亚迭代时间法,正确地描述了30P30N多段翼型前缘结构区域的流动特征和流动机理。     

基于Hamilton-Jacobi方程的编队飞行控制

 针对椭圆参考轨道的编队飞行控制系统,基于哈密尔顿 雅克比（Hamilton-Jacobi）方程,给出一种新的生成函数半解析近似解的计算方法,同时推导了编队飞行控制系统相对运动的线性状态转移矩阵,并且验证了模型的精度。在此基础上,结合一个具体的编队飞行任务,采用线性规划算法得到了队形保持的闭环控制器,该算法得到的是脉冲形式的控制序列,因此更具有实际意义。同时还提出了一种简单的避免生成函数奇异点的方法,从而减少了生成函数间相互转换所带来的计算量。最后的非线性仿真结果验证了控制方法的有效性。     

战损飞机动力学特性分析及杀伤研究

现代战争中,防空武器性能不断提升,对飞机造成的威胁越来越大,确定飞机受到打击时的杀伤状态及杀伤等级成为飞机设计时需要考虑的重要问题。飞机被防空武器击中时,可能发生机翼脱离损伤,其质量、面积产生明显变化,将导致飞机的重心、气动特性的变化,严重危害飞机的稳定性。针对上述情况,以某型作战飞机为例,建立损伤情况下的飞机模型,综合分析机翼脱离损伤对飞机重量、重心和机翼面积带来的影响;建立损伤飞机的气动模型,计算损伤飞机气动特性的变化情况。基于这些变化,分析损伤对飞机动力学特性的影响,建立损伤状态下的动力学方程,计算损伤飞机的失控时间,依此确定飞机的杀伤等级。本文所做研究可用于在飞机设计过程中分析其杀伤状态及杀伤等级,为飞机设计时飞机杀伤的分析评估提供一种有效的解法。     

旋转机翼悬停气动特性研究

鸭式旋转机翼(CRW)是一种先进的高速直升机方案,旋翼同时也是机翼,是其最关键气动部件之一。采用结构分区拼接网格技术进行空间离散,分区建立参考系,通过求解多重参考系下的N-S方程来计算旋翼流场,首先以传统的Caradonna-Tung实验旋翼的亚、跨声速悬停流场分析为例验证该方法的可靠性,进而采用该方法对旋转机翼悬停流场进行了数值计算,旋翼拉力计算值和地面实验值吻合较好,结果分析表明旋转机翼的悬停流场有着不同于传统旋翼的流场特性。     

Friction moment analysis of space gyroscope bearing with ribbon cage under ultra-low oscillatory motion

This paper presents the model of calculating the total friction moment of space gyroscope ball bearings which usually work under ultra-low oscillatory motion and are very sensitive to the friction moment. The aim is to know the proportion of the friction moment caused by each frictional source in the bearing＇s total friction moment, which is helpful to optimize the bearing design to deduce the friction moment. In the model, the cage dynamic equations considering six degree-of-freedom and the balls dynamic equations considering two degree-of-freedom were solved.The good trends with different loads between the measured friction moments and computational results prove that the model under constant rate was validated. The computational results show that when the speed was set at 5 r/min, the bearing＇s maximum total friction moment when oscillation occurred was obviously larger than that occurred at a constant rate. At the onset of each oscillatory motion, the proportion of the friction moment caused by cage in the bearing＇s total friction moment was very high, and it increased with the increasing speed. The analyses of different cage thicknesses and different clearances between cage pocket and ball show that smaller thickness and clearance were preferred.     

曲线物面边界条件在DGM中的应用研究

基于二维非结构网格,采用高精度间断Galerkin有限元方法数值求解定常Euler方程。为了有效克服传统的反射物面边界条件在物面处易于生成伪熵降低计算精度的缺陷,提出一种曲线物面边界条件,从而实现对物理模型而不是数值模型进行数值模拟。对经典圆柱亚音速无粘绕流进行数值模拟,结果表明:采用曲线物面边界条件之后,流场解的精度得到很好的提高;此外,在非常稀疏的网格上,通过提高基函数的阶次仍可以得到高精度的数值解。     

考虑背景孔隙的单开孔两空间结构的风致内压响应研究

风致内压对建筑结构的安全性有着显著的影响,在结构抗风设计时不容忽视。本文基于合理的假定,利用非定常伯努利方程、质量守恒定律和绝热气体状态方程推导了有背景孔隙的单开孔两空间结构的内压响应非线性控制方程组,并通过数值算例分析了背景孔隙对内压响应的影响。结果表明:本文导出的内压控制方程组可很好地分析有背景空隙的单开孔两空间结构的风致内压响应;无背景孔隙时,内侧房间的内压响应要高于外侧房间;背景孔隙所引入的附加阻尼使得两房间的内压响应均受到抑制。这些研究成果对结构风致内压有理论与实际应用的价值。     

基于自适应差分进化算法的变循环发动机模型求解方法研究

为了降低传统迭代算法在求解变循环发动机非线性模型时对初值的依赖性,将模型的求解问题转换为求最小值的优化问题,引入差分进化算法进行模型的求解,并提出一种自适应差分进化算法(ADE)。ADE借助轮盘赌选择法,利用种群的进化经验可以自适应的选择最适合当前种群的差分策略与算法控制参数。针对变循环发动机四个典型工作点的模型求解问题,研究了标准差分进化算法(SDE)的控制参数对其性能的影响,获取了SDE在求解四个典型工作点时的最优控制参数组合,对比分析了ADE与SDE的性能差异,最后研究了种群规模对ADE性能的影响。结果表明：SDE在求解发动机模型时具有较好的鲁棒性,在求解不同工作点时算法的最优控制参数并不完全相同;相比于使用最优控制参数的SDE,ADE可以在不影响算法鲁棒性的情况下提升效率50%以上;减少ADE的种群规模会在提升算法效率的同时破坏鲁棒性。