应用改进的低雷诺数湍流模型预测转捩流动

基于Launder-Sharma(LS)低雷诺数k-ε两方程湍流模型发展了一种改进的具有转捩敏感性的低雷诺数湍流模型.针对LS模型的可实现性(realizability)问题和前缘滞止点湍动能预测过大的不足,模型进行了改进.改进模型只使用当地物理量,不需要求解壁面距离、y+和边界层积分参数.改进模型能够适用于广泛的流动,且容易应用到通用的计算流体动力学(CFD)程序中.对具有详细数据的零压力梯度平板转捩边界层T3A实验的模拟结果显示,改进模型能够预测转捩流动,并能对自由湍流变化给出合理的响应.      

叶栅二次流旋涡结构与损失分析

采用三维粘性程序对某型动力涡轮的第一级进行了数值模拟, 模拟结果捕捉到了该涡轮级叶栅的内部流的流动细节, 展示了涡轮叶栅端壁和型面流动及叶栅通道内的三维流动结构.通过对叶栅中的二次流现象和流动损失机理的分析, 揭示了该涡轮级叶栅通道内二次流旋涡结构(马蹄涡、通道涡、壁角涡、尾迹涡、泄漏涡等)的演变过程, 以及旋涡结构对损失分布的影响.      

液氢空化流动特性研究

对液氢空化流动的实验数据进行了统计分析,数值计算了大范围温度变化下液氢空化流动的变化特征。结果表明采用的数值计算方法计算结果与实验数据吻合良好。随着温度变化,液氢空化流动存在两种空化行为。当温度低于转捩温度时,空化行为的变化受液/汽密度比控制,空穴面积随温度的增加而增加;当温度高于转捩温度时,空化行为的变化热力学效应控制,空穴面积随温度的增加而减小。空化数和来流速度越大,转捩温度越高。     

三维可压和不可压缩流动数值模拟的一种预处理多重网格方法研究(英文)

将矩阵预处理方法与多重网格方法结合,发展了适用于机翼及复杂外形气动数值优化设计的高效数值模拟方法和程序。在不改变定常流动解的前提下,对 Navier-Stokes 方程的时间导数项实施 Choi-Merkle 矩阵预处理,从而改善可压缩控制方程在低速情况下的系统刚性。Choi-Merkle 预处理最初用于二维低马赫数粘性流动数值模拟,本文将其推广应用到计算三维流动,并针对提高跨音速流动计算收敛性进行了相应的改进。预处理后控制方程采用中心格式有限体积进行空间离散,Runge-Kutta 混合多步方法进行显式时间推进求解,并应用 FAS 多重网格方法加速收敛。采用所发展的方法和程序数值模拟了绕 M6 机翼和 M100 机翼的低速和跨音速流动(马赫数从0.01 到 0.839)以及绕 F4 翼身组合体流动。计算结果与实验值吻合良好,验证了所发展的方法的正确性。算例研究表明所发展方法较大幅度地提高了三维流动数值模拟的效率,且同时适用于可压和不可压流动计算,在机翼及复杂外形气动优化设计方法具有应用前景     

钝体俯仰阻尼导数数值计算

分别采用：（1）Euler方程的谐振摄动法；（2）薄层近似的非定常N－S方程模拟有粘扰动流场两种方法计算了钝体外形的高超声速俯仰阻尼导数。通过计算表明：对于复杂外形，方法（1）适用于零攻角附近的流动计算，而方法（2）考虑了流场的粘性作用，能够应用于一般流态（3＜M＜20，0＜a＜20°）的计算。此外，本文发展了文献［6］的动导数转换公式，使之适用于二维转轴情况，并在方法（2）中通过引入“亚迭代”过程提高了非定常流场的计算精度。     

高超声复杂流动中湍流模式应用的评估

本文选择几个高超声速基准流动：二维可压缩拐角、锥柱裙组合体绕流、斜激波与湍流边界层干扰，采用几个常见的湍流模式，BL模式，CHκ-ε模式，κ-ε模式，SST模式，CMOTTκ-ε模式，SHIHκ-ε模式，通过将数值计算结果和实验结果进行比较，对有关的湍流模式地 评估，得到一些有意义的结论。     

基于非线性耦合本构关系的改进边界条件

非线性耦合本构关系（NCCR）模型是在Eu的广义流体动力学方程（GHE）基础上,通过绝热假设、Eu封闭和Myong简化推导出的关于非守恒量（黏性应力与热流）的非线性代数方程,有效拓展了线性的纳维-斯托克斯-傅里叶（NSF）本构模型在非平衡流动中的模拟能力,为快速准确模拟连续与稀薄耦合流动问题提供了强有力的理论工具。针对该模型开展滑移边界条件研究,结合努森层内物理量非线性分布的特点,提出一套在物面处与模型精度相一致的非线性修正滑移边界条件。在有限体积框架下,采用AUSMPW⁺格式和LU-SGS方法以及NCCR的完整耦合求解算法,对不同稀薄程度的高超声速单原子氩气圆柱绕流和平板绕流问题进行数值模拟。研究结果表明,基于NCCR模型的修正边界条件准确刻画出物面努森层内流动的非线性特点,有效提高了固壁滑移边界的精度。采用非线性修正边界的NCCR模型准确预测了连续流、滑移流和过渡流域的物面压力、摩阻与热流系数。     

高速流前端磁场B_z分量的变化

等离子体片高速流在磁层活动中起着重要作用,其形成机制以及与背景等离子体的相互作用日益引起关注.本文利用搭载于Cluster四颗卫星上的磁场和等离子体观测仪器,对2001和2002两年发生在磁尾等离子体片中高速流事件期间的磁场变化进行了统计研究.结果表明,在高速流前端,伴随着等离子体整体速度的增加,绝大多数高速流前端磁场的B_z分量经常出现先短暂减小然后剧烈增大的现象,符合以往利用Geotail卫星观测数据获得的统计特性.然而个例研究发现B_z的下降与上升常常是不对称的,且B_z分量下降的程度并不是总能达到反向的程度,说明这种变化特征并不一定是存在磁结构的表现.我们认为更多时候这种磁场的变化特征是高速流挤压背景等离子体磁场造成的,是相互作用的结果.当偶极化锋面形成后,由类似间断面的磁场界面反射的热离子产生抗磁效应,可能对B_z下降形成部分贡献,而B_z增加则是高速流携带磁通量堆积的效果.     

网络最大流的图单纯形解法

提出网络饱和流的概念，重新定义了网络最大流问题，并提出了求解最大流的图单纯形算法。该方法避免了２Ｆ算法的缺点，其计算复杂性为Ｏ（２ｍｎ）。     

机场终端区跑道配置优化及容流调配模型研究

大型繁忙机场交通需求的持续增长导致的飞行流量与保障能力、机场容量之间的矛盾日益突出。为了充分利用跑道系统资源,合理配置跑道运行容量,优化飞行流,建立非线性0-1整数规划模型解决以下两个问题:确定跑道配置优化序列,匹配进离场飞机流。模型综合考虑机场交通流、跑道配置转换容量折损、跑道容量包络线等约束,以优化区间内航班总延误最小为目标,用LINGO建模求解,使用实际运行数据验证模型的有效性。结果表明,模型实现了跑道资源的优化利用,降低了航班延误。