基于熵限制的Baldwin-Lomax湍流模型

针对Baldwin-Lomax(BL)模型在模拟超声速复杂流动中的不足,提出一种基于熵限制的新型修正方式.分别采用原始的BL模型(BL-origin)、修正后的BL模型(BL-entropy)以及Wilcox k-ω两方程模型对超声速平板、超声速二维压缩拐角以及膨胀-压缩拐角3个典型算例进行计算,结果表明BL-entropy 能够较好地获取长度尺度进而得到均匀合理的涡粘性分布.同时,该修正方式简单高效并且适用于其他与超声速边界层有关的模型,具有较大的发展潜力.     

微重力量值对气液两相流相分布及液相湍流统计量影响的数值研究

气液两相流在空间领域具有广阔的应用前景, 深入理解微重力量值(微重力大小)对相分布和液相湍流的影响十分必要. 采用欧拉elax-elax拉格朗日双向耦合模型深入研究了不同微重力量值对相分布和液相湍流的影响. 液相速度场通过直接数值模拟求解, 气泡的运动轨迹由牛顿运动方程跟踪. 研究表明, 气泡分布和液相湍流与微重力量值均具有直接联系. 在低微重力量值下, 气泡近似均匀分布在槽道内, 且对液相湍流统计量几乎没有影响; 然而当微重力量值较高时, 大量气泡聚集在壁面附近, 液相湍流由于气泡的注入受到极大调制.      

Von Karman模型三维大气紊流仿真理论与方法

基于三维紊流场相关函数矩阵,用蒙特卡罗法仿真生成Von Karman模型三维空间大气紊流场,实现了对复杂的Von Karman模型的精确推导和应用,仿真流场的相关特性在理论上可以任意逼近Von Karman理论模型.通过提高白噪声序列的白化程度及优化蒙特卡罗算法的方法提高了仿真精度和效率.数值仿真结果证明:所生成的空间大气紊流场符合模型特性,具有较好的统计特性和仿真精度,满足实时飞行仿真对紊流场数据的要求.     

激光通信中的分集检测技术

在大气弱湍流模型下,应用Monte Carlo方法,通过产生符合光信号的对数幅度联合概率分布的数据,近似模拟计算了分集接收机的逐符号最大似然检测算法的似然比.研究了减少Monte Carlo方法中仿真点的个数对逐符号最大似然检测的分集接收机检测性能估计的影响.随后提出一种基于训练序列的分集检测方法.该方法在未知衰落分布相关信息的情况下,利用每帧内引入的一定长度的训练序列提供了某一较短时间间隔内的大气湍流信道的信息,基于与似然比表达式结构相似的算法获得改善的检测效果.仿真结果表明在各个分集接收机之间存在较强相关性的情况下,该算法仍可实现有效的检测.      

用GAO-YONG湍流模式对翼型跨音粘流的数值模拟

应用GAO-YONG可压缩湍流模式,数值模拟了NACA0012,RAE2822翼型的定常跨音速粘流算例.对流项采用三阶ROE格式,扩散项采用二阶中心格式,用多步Runge-Kutta显式时间推进法求解空间离散后的控制方程.计算结果很好地预测了翼型表面的压力系数的分布、激波的位置、马赫数等值线的分布等情况,并且对翼型表面激波与边界层相互干扰以及层流向湍流的转捩问题进行了分析计算.计算结果与实验值符合很好,表明GAO-YONG可压缩湍流模式应用合适的计算方法能够高精度模拟翼型跨音粘性流动问题,并且基于GAO-YONG可压缩湍流模式各向异性湍流粘性的机理,提供了一种预测转捩起始位置的判别方法.      

应用修正湍流模型对旋转通道内换热的模拟

针对旋转光滑矩形通道分别应用针对旋转状态修正的 k-ε、标准k-ω以及提出的针对旋转状态修正的k-ω湍流模型进行流动和换热的数值模拟,通过与实验结论的对比,讨论了采用不同湍流模型对计算结果的影响.计算工况为旋转数Ro=0.24,流体进口雷诺数Re=25000.计算结果表明:采用所提出的针对旋转状态修正的k-ω湍流模型的计算结果要比采用针对旋转状态修正的k-ε以及标准k-ω湍流模型的计算结果更接近实验结论.     

标量图像测速原理及数值检验

标量图像测速是一种新型的能够测量微小结构的湍流实验测量系统,其理论依据是Schmidt(Sc)数>1的标量湍流里,标量场脉动的时间尺度和空间耗散尺度要小于对应的速度场脉动的时间尺度和空间耗散尺度,标量脉动比速度脉动具有更强的间歇性,标量场所含信息大于速度场所含信息,从高Sc数标量场提取速度场是可能的.笔者详细介绍了标量图像测速原理,并在积分最小化法思想指导下根据标量输运守恒方程及限定条件,建立了积分最小化标量图像测速的数学表示.通过变分法,获得了积分空间每个点的三个速度分量的欧拉特征方程,对欧拉特征方程离散化即可在积分空间建立一线性稀疏方程组.采用变分迭代法对此线性稀疏方程组进行求解即可获得三维速度场,据此采用积分最小化法建立了标量图像测速图像处理系统.还介绍了四维标量场(三维空间加一维时间)的图像采集原理,并采用DNS数据对标量图像测速处理系统进行了数值检验,检验结果表明标量图像测速技术所求解的速度场和DNS精确解之间在结构上是十分相似的,标量图像测速可以正确地提取速度场.     

空气龄的数值解法

提出了用ＣＦＤ法求解空气龄的具体方法，其步骤为先求解Ｎ－Ｓ方程，后求解空气龄。为求解Ｎ－Ｓ方程，采用了改进型ＨＳＭＡＣ法，并给出了压力和速度同时修正的方法。为求解空气龄，推出了空气龄输送方程式和给出了空气龄和速度成分的边界条件。计算机模拟验证了求解方法的有效性并由空气龄求出了通风效率。     

波瓣喷管混合流动的数值计算

借助Ｎ－Ｓ方程对波瓣喷管混合装置的流场进行了数值分析。文中采用三维贴体曲线坐标,非交错网格法和代数湍流模型,对波瓣喷管内外两侧的气流通道及波瓣喷管下游的混合段波场采用分块计算方法进行了数值计算。     

B/L湍流模型在强压力梯度流场计算中的应用

对在湍流计算中已得到广泛应用的Ｂａｌｄｗｉｎ／Ｌｏｍａｘ（Ｂ／Ｌ）代数湍流模型进行了修正。修正后的Ｂ／Ｌ模型考虑了压力梯度及物面引射速度对内层ＶａｎＤｒｉｅｓｔ衰减因子的影响。与此同时,文中还用修正的Ｂ／Ｌ湍流模型对亚声平板附面层流动、拉伐尔喷管内的正激波／湍流附面层干扰流动进行了计算。计算与实验结果比较表明,修正后的Ｂ／Ｌ模型不仅可以更准确地计算具有强压力梯度的激波／附面层干扰区内流动压力分布,而且还可以准确地预报激波／附面层干扰引起的流动分离,从而为Ｂ／Ｌ模型的工程应用开拓了领域。