预点火湍流对正戊烷云雾爆炸参数的影响

以正戊烷云雾为研究对象，进行预点火湍流对云雾爆炸参数影响规律的实验研究。首先通过不同气动压力进行喷雾，获得平均特征直径（SMD）分别为 21.21、14.51 和 8.64 μm 的正戊烷云雾，并得到不同气动压力预点火的湍流均方根速度；随后在 20 L 云雾爆炸参数测量系统中实验获得预点火湍流对正戊烷云雾蒸发速率、爆炸超压峰值、压力上升速率和火焰传播延迟时间的影响。结果表明:(1) 对于圆柱形罐体对称式双喷头分散系统，流场环境可近似认定为零平均速率湍流场；在0.4、0.6和0.8 MPa的气动压力喷雾50 ms的分散作用下，在100～250 ms内，湍流均方根速度在1.0～6.2 m/s范围内，平均湍流积分尺度在40～72 mm范围内，湍流最大湍流尺度的雷诺数在8 000～15 000范围内，柯尔莫哥洛夫微尺度在0.03～0.1 mm范围内；(2) 对于较小的液滴群，随湍流强度的增加，液滴群的蒸发速率有更为明显的提升；(3) 对比云雾三种SMD，粒径8.64 μm的超压峰值与最大压力上升速率随湍流强度增长趋势更显著，并发生爆炸强度显著提升现象，即存在“转变区域”（transition range）现象；(4) 对于SMD在8～22 μm范围内，湍流均方根速度处于1.0～4.0 m/s时为火焰传播延迟时间的低增长阶段，湍流均方根速度处于4.0～6.2 m/s时为火焰传播延迟时间的高增长阶段，湍流强度与火焰传播延迟时间在相应的两个湍流强度阶段范围内呈线性增长。     

非球形颗粒两相流的数值模拟研究进展

非球形颗粒两相流是多相流的重要研究方向之一, 常见于自然界及工业生产过程中. 不同于球形颗粒, 由于非球形颗粒形状的各向异性, 除了颗粒平动行为, 还需要考虑颗粒的转动与取向行为, 颗粒的取向与转动行为会影响颗粒所受的力和力矩. 为了准确模拟非球形颗粒的运动行为, 目前非球形颗粒两相流的数值模拟研究主要基于欧拉?拉格朗日的求解框架展开, 常见的非球形颗粒两相流数值模拟方法主要包括点颗粒法与全分辨颗粒法. 本文将对这两类方法进行介绍, 同时会全面介绍非球形颗粒两相流研究的基础理论模型, 并系统总结非球形颗粒在简单基本流和复杂湍流中的研究进展, 包括对于非球形颗粒在湍流中的取向与转动行为机理, 以及颗粒对湍流减阻调制作用的研究. 最后, 本文提出了非球形颗粒两相流研究存在的问题及未来研究方向.      

相位片湍流仿真装置中的激光漂移效应

 为了在内场中开展激光传输特性实验研究,设计了一种基于Kolmogorov谱的随机相位片的大气湍流物理仿真系统,并用该大气湍流模拟仿真系统对激光漂移效应进行了研究。采用光束质量分析仪对模拟传输的光斑图像进行采集和分析,通过比较实际漂移方差与理论漂移方差,验证了该湍流仿真系统模拟和评价光束在仿真传输中的性能。实验分析表明,基于湍流相位片物理仿真系统模拟的光束漂移方差与理论曲线基本一致。     

Derivation of a second-order model for Reynolds stress using renormalization group analysis and the two-scale expansion technique

With the two-scale expansion technique proposed by Yoshizawa,the turbulent fluctuating field is expanded around the isotropic field.At a low-order two-scale expansion,applying the mode coupling approximation in the Yakhot-Orszag renormalization group method to analyze the fluctuating field,the Reynolds-average terms in the Reynolds stress transport equation,such as the convective term,the pressure-gradient-velocity correlation term and the dissipation term,are modeled.Two numerical examples:turbulent flow past a backward-facing step and the fully developed flow in a rotating channel,are presented for testing the efficiency of the proposed second-order model.For these two numerical examples,the proposed model performs as well as the Gibson-Launder (GL) model,giving better prediction than the standard k-ε model,especially in the abilities to calculate the secondary flow in the backward-facing step flow and to capture the asymmetric turbulent structure caused by frame rotation.     

壁面湍流流动中高分子减阻等效粘度模型的验证

近来在壁面湍流高分子减阻研究中,一种拉伸的高分子产生自相一致的等效粘度的理论提了出来,这个等效粘度随离开壁面的距离而增长。本文将此线性分布等效粘度置入Navier-Stokes方程,运用雷诺应力模型计算在壁面湍流中的减阻情况,检验这种等效粘度的可行性。可以发现,此模型可以得到湍流减阻的效果,所得到的减阻率随着等效粘度线...     

翼型、机翼非定常运动模式下动态数值模拟

基于RANS方程,通过刚性动网格技术实现对翼型和机翼典型运动模式的描述,采用双时间推进方法和Roe空间离散格式对流场求解,构建了一个非定常气动计算平台;以NACA0012翼型为算倒进行了动态数值模拟可信度验证。数值模拟结果与试验数据吻合较好,升力和俯仰力矩的最大计算误差分别为3％和10％,表明了该平台的可靠性。另外,还数值模拟了M6机翼的动态非定常流场,并分析了两种湍流模型对非定常流场激波的捕捉能力。结果表明非定常流动中S-A湍流模型对激波的捕捉较B-L模型更敏感。文中开发的非定常计算平台对进一步解决三维复杂流场的流动问题有很高的工程应用价值。     

基于混合网格Navier-Stokes方程的并行隐式计算方法研究

针对结构网格很难处理复杂外形和非结构网格无法计算具有边界层的粘性流动的缺点,发展了基于混合网格格点的隐式算法,成功地解决了在工程应用中难于处理的复杂外形粘性流场计算和效率问题。同时针对大规模的工程问题,发展了基于MPI通信技术的染色分层通讯并行计算方法。其中空间离散采用基于Roe格式发展的三阶迎风HLLEW(Harten-Lax-Van Leer-Einfeldt-Wada)或AUSM格式,湍流模型采用k??两方程湍流模型,时间推进考虑到LU-SGS并行等效较困难则采用基于DP-LUR(Data-Parallel Lower-Upper Relaxation)格式的隐式算法,计算CFL数可取到105量级,从2个到128个CPU的并行加速效率都保持在90%以上,大大提高了计算效率。算例对标模M6机翼模型流场进行计算,验证了方法的可靠性;然后对标模DLR-F6翼身组合体进行混合网格粘性与无粘计算结果进行比较,进一步验证混合网格方法;最后计算了DLR-WBNP外挂发动机翼身组合体模型,准确模拟了外挂和超临界机翼的相互干扰流动问题,采用4 CPU 16 CORE到24 CPU 96 CORE,2000步计算时间都不超...     

k-ε双方程湍流模型对制退机内流场计算的适用性分析

为研究不同双方程湍流模型对制退机内复杂流场计算的适用性,以某火炮制退机为研究对象,建立了实际结构下的三维计算模型,利用动网格与滑移网格技术,实现了火炮实际后坐速度下的制退机内部三维运动流场的数值计算。分别采用标准k-ε模型、RNGk-ε模型和Realizable k-ε模型计算制退机内部各腔室压力,与实验曲线对比,结果表明,应用标准k-ε模型对后坐冲击过程的制退机内部压力计算的误差最小,与实验结果吻合最好。     

k-ε双方程湍流模型对制退机内流场计算的适用性分析

为研究不同双方程湍流模型对制退机内复杂流场计算的适用性,以某火炮制堪机为研究对象,建立了实际结构下的三维计算模型,利用动网格与滑移网格技术.实现了火炮实际后坐速度下的制退机内部三维运动流场的数值计算.分别采用标准k-ε模型、RNGk-ε模型和Realizablek-ε模型计算制退机内部各腔室压力,与实验曲线对比,结果表...