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

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

高分子溶液壁湍流减阻机理的TRPIV实验研究

利用高时间分辨率粒子图像测速技术(TRPIV)对回流式水槽中低浓度高分子溶液壁湍流的减阻机理进行实验研究。通过对比分析高分子溶液和纯水平板湍流边界层在相同来流速度下的平均速度剖面、湍流强度和雷诺应力,发现高分子溶液的壁面摩擦阻力减小了21.77%,并且其缓冲层增厚,按对数律外移,雷诺应力减小;高分子聚合物主要在近壁区起到抑制湍流脉动的作用,而在主流区的作用不太明显。用流向局部平均多尺度速度结构函数和相干结构条件采样方法,检测并对比了高分子溶液和水的壁湍流相干结构“喷射”和“扫掠”事件中的脉动速度、展向涡量、雷诺应力等物理量的二维拓扑形态,发现高分子溶液近壁区相干结构在猝发时的脉动速度减小,涡量受到抑制,雷诺应力明显减小,说明高分子溶液湍流近壁区相干结构“喷射”和“扫掠”的强度变弱,猝发频率降低,动量和能量的输运减弱,揭示出高分子溶液减阻的重要机理。     

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

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

湍流加速火焰的三维数值模拟

火焰在设有障碍物的管内传播时会自身加速，并可能导致爆炸。本文基于湍流κ-ε模型和改进的EBU—Arrhenius反应模型，对该现象进行了三维空间的数值模拟。计算结果反映了障碍物、湍流和火焰之间相互作用的正反馈机理，描绘了火焰在管内加速传播的三维图像。     

湍流的诱导及其对瓦斯爆炸过程中火焰和爆炸波的作用

在实验的基础上,研究了管内瓦斯爆炸过程中湍流的诱导及其对瓦斯爆炸过程中火焰和爆炸波的影响作用.研究结果表明,管道面积突变对瓦斯爆炸过程中湍流的产生具有重要影响.管道面积突变(变大、变小)时,产生附加湍流,并使下游火焰气流的湍流度增加,瓦斯爆炸过程中火焰的传播速度迅速提高,并可诱导激波的产生.在80×80mm等截面直管中(瓦斯浓度为理论上最猛烈的爆炸浓度9.5%),瓦斯爆炸最大火焰传播速度为40.8m/s,管内各点均为压力波信号,当管道加装一Φ300mm圆管形成面积突扩11倍和突缩11倍两断面后,面积突扩处(L/D=22)火焰速度增大5.05倍,达到64.4m/s,面积突缩处(L/D=28)火焰速度为156.0m/s, 增大4.55倍,并在L/D=48倍处形成激波(超压1.6976atm、波速416.7m/s),在L/D=98倍处,激波强度最大.在面积突变管内加装加速环可使瓦斯爆炸过程中湍流度加剧,火焰的传播速度更高,激波生成的位置(L/D=28)、最强点位置(L/D=70)均前移,激波强度增大.研究结果对指导现场如何防治瓦斯爆炸,减轻瓦斯爆炸的威力具有一定的指导意义.     

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

以正戊烷云雾为研究对象,进行预点火湍流对云雾爆炸参数影响规律的实验研究。首先通过不同气动压力进行喷雾,获得平均特征直径（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时为火焰传播延迟时间的高增长阶段,湍流强度与火焰传播延迟时间在相应的两个湍流强度阶段范围内呈线性增长。     

HyperFLOW软件非结构网格亚跨声速湍流模拟的验证与确认

计算流体力学(computational fluid dynamics,CFD)数值模拟在航空航天等领域发挥越来越重要的作用,然而CFD数值模拟结果的可信度仍然需要通过不断地验证与确认来提高.本文给出了从制造解精度测试、简单到复杂外形湍流模拟网格收敛性研究等三个方面开展CFD软件验证与确认的方法,并对自主研发的CFD软件平台HyperFLOW在非结构网格上模拟亚跨声速湍流问题的能力进行了验证与确认.首先通过基于Euler方程和标量扩散方程的制造解精度测试,分别验证了HyperFLOW在非结构网格上对Euler方程和黏性项的求解精度,结果表明其能够在任意非结构网格上达到设计的二阶精度. 其次,通过NASATurbulence Modeling Resource中的湍流平板、二维翼型近尾迹流动、二维Bump等几个典型的亚声速湍流算例的网格收敛性研究,量化考察了数值结果的观测精度阶和网格收敛性指数,并与国外知名CFD解算器CFL3D,FUN3D的计算结果进行了对比,验证了HyperFLOW对简单湍流问题的模拟能力,且具有良好的网格收敛性和计算精度(阶). 最后,通过NASA CommonResearchModel标模定升力系数的网格收敛性研究和升阻极曲线预测,验证了软件在复杂外形亚跨声速湍流流动数值模拟中也具有良好的可信度.      

大气边界层的一些空气动力学特征

本文在定常和中性温度层结条件下,首先介绍了均匀平坦地面上大气边界层的一些空气动力学特征,并将大气边界层与一般空气动力学边界层作了比较。然后介绍了地面粗糙度分布有阶跃变化以及非平坦地形上的大气边界层特征。文中说明,气流的平均流场结构以及湍流特性与均匀平坦地面的情形相比有显著不同。最后,提出了一些有待今后进一步研究的问题。      

应用谱理论研究压力脉动的湍流模式

本文应用谱分析理论研究了剪切湍流场中的压力脉动,包括功率谱、均方值等.通过对压力脉动Possion方程的Fourier变换,首先将压力脉动谱表示成速度脉动谱的形式.利用Navier-Stokes方程的形式解及准正态分布假设,可以进一步将压力脉动功率谱表达式中所包含的速度脉动的三阶相关与四阶相关表示成速度脉动的二阶相关（功率谱）.最后,引入高雷诺数流的速度脉动功率谱模型,导出了由湍动e₀,耗散ε,雷诺应力-iu_j>及时均速度梯度表示的压力脉动均方值的湍流模式,并同现有数据进行了比较.     

用可实现k-ε模式对细长体大攻角分离流场的数值模拟

利用一种符合可实现性条件的k-ε两方程模式,通过求解可压缩性Navier-Stokes方程,数值模拟了细长旋成体在M=2.5和0.7、中等大攻角下的湍流绕流流场。定性分析了背风面分离涡的发展过程和湍流粘性比等流动特征,对背风面主涡强度和涡核位置、物面压力分布、集中力和力矩与已有的试验结果进行了细致的定量对比,证实了数值模拟的有效性。