一种TVD方法在三维非结构网格中的应用

以三维非结构网格的显式有限体积法为基础,采用了一种TVD方法求解三维Euler方程,使用Roe通量来计算网格单元边界处的守恒量通量.为了验证方法的可行性,用该方法模拟三维爆炸问题,得出的结论是我们的方法可行,稳定且有效.     

有二次涡的激波边界层干扰流动的DPNS和NS方程计算

本文针对具有二次涡复杂分离再附现象的激波边界层干扰流动,数值地考察了扩散抛物化Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ（ＤＰＮＳ）方程组的适用情况。壁面摩阻和压力、主涡和二次涡的涡高和涡长、分离再附位置以及流线图等特性的计算表明：ＤＰＮＳ方程组的数值结果均与ＮＳ方程组的数值结果很好相符。     

全流量补燃循环发动机系统的响应特性

以全流量补燃循环氢氧发动机系统为研究对象,对其的动态响应特性进行了研究。建立了描述全流量补燃循环发动机动态特性的非线性数学模型,将免疫策略算法同龙格-库塔法结合起来,提出了求解含有隐式项的常微分方程组的变步长龙格-库塔方法,并应用该方法对全流量补燃循环发动机系统的动态响应特性进行了仿真计算。计算结果表明,当发动机在某一个稳定工况工作时,发动机入口推进剂压强的变化对发动机性能参数的影响不大,发动机的参数都能比较平稳的过渡到一个新的稳定工况;当发动机在短时间内进行小范围的推力调节时,发动机参数的过渡过程的曲线也比较平稳,但是当在短时间内进行大范围推力调节时,参数的过渡过程的曲线振动比较剧烈,因此应当在进行大范围推力调节时,应当对调节时间进行适当延长,或者分级进行推力调节。      

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

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

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

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

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

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

基于成分的燃烧室出口燃气温度计算方法

针对燃烧室不同温度下燃气组成的不同,提出建立不同燃气模型以推算燃气温度的两种计算方法。基于物质守恒、化学平衡和能量守恒原理建立温度计算方程,根据能量守恒方程的不同表达形式给出平均比定压热容法和焓值法,同时给出两种求解非线性方程组的数学解法:三变量迭代法和牛顿法,并采用VB(Visual Basic)编制燃气分析测试程序实现温度算法。列举计算实例比较两种算法的区别,结果显示:随着油气比的增加,两种计算方法得到的燃气温度相差越大,考虑热离解的高温算法结果更小,原因在于燃气中离解组分的摩尔分数随着温度升高而增加,燃烧实际放热量减小,导致燃气温度减小,所以当燃气温度较高时建议采用高温算法更符合实际燃烧情况。      

使用GAO-YONG湍流方程组计算翼型分离流

采用SIMPLE方法求解GAO-YONG不可压湍流方程组, 对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟, 给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数.与实验数据以及大涡模拟结果的比较表明, GAO-YONG不可压湍流方程组能够对翼型绕流的分离点、分离涡形态、表面压力分布、升阻特性做出较好的预测, 能够模拟翼型大攻角分离流动, 计算结果优于FLU-ENT软件中的k-ε、k-ω、k-ωsst模型的计算结果.      

求解非线性方程及方程组的模拟退火算法

探讨了用模拟退火算法求解非线性方程及方程组的方法,该方法把方程组求解问题转化为函数优化问题.计算中不需要使用目标函数的导数信息,在Matlab环境下实现了该算法,数值实验结果表明了该算法的有效性.