气固两相圆湍射流中湍流调制的数值模拟

对空间发展的气固两相圆湍射流进行了双向耦合的大涡模拟,系统地给出了不同颗粒尺寸和不同气载比时颗粒相对气相圆湍射流的调制规律.结果表明: 在下游区域,所有尺寸的颗粒均使气相湍流削弱,而在上游区域,颗粒对气相的调制存在临界尺寸.该临界尺寸随颗粒材料密度的增加而减少,而气载比的增加会加剧颗粒对气相的调制程度.     

入口扰动对圆湍射流流动大涡模拟预报的影响

运用大涡模拟方法分别对不加扰动和加扰动的圆湍射流流动进行了数值模拟,比较了两者在大尺度涡结构的演化过程以及时均速度和湍流强度等统计结果方面的预报差异,并与实验结果进行了对照。模拟结果表明,施加了入口扰动的模拟结果在物理真实性方面和统计结果的准确性方面均要明显优于未加扰动的模拟结果。分析了不加扰动和加扰动模拟的本质差异,指出准确预报实际圆湍射流流动应采用施加入口扰动的三维模拟。     

矩形射流流动的大涡模拟

为探讨矩形射流中的拟序运动特征,用大涡模拟方法(LES)对三维矩形射流流动过程进行了数值模拟。统计平均的流向速度横向分布曲线在发展区某些位置上呈现"马鞍型"分布,速度峰值偏离中心轴线,与圆射流显著不同,这与已有的实验和直接数值模拟(DNS)结果吻合良好。计算结果表明,在发展区湍流脉动最强,且应力的横向分布不均匀。随着流向距离的增大,速度功率谱(PSD)曲线呈现峰值向低频移动的现象。速度的PSD分布表明:矩形射流中拟序运动主要位于射流的核心区末端和发展区,射流发展区中的拟序运动具有明显的时间准周期性。     

平面自由湍射流拟序结构的大涡模拟研究

对空间发展的平面不可压缩湍射流拟序结构的非定常演化过程进行了大涡模拟。采用标志物浓度等值线分布示踪平面射流气相流动 ,Reynolds数为 1130 0 ,模拟结果再现了平面射流中 Kelvin- Helm holtz不稳定性的发生以及展向大尺度涡的卷起、合并、破碎过程。捕获到了射流拟序结构剧烈相互作用的“偶极子”和“三极子”现象。在拟序结构的配对过程中 ,大尺度旋涡的尾迹混合形成了流向“发卡”型拟序涡结构。平面射流拟序结构的时空演化过程所呈现出的卡门涡街特征主要集中于初始段的后部和过渡段的前端。数值流场显示结果和染线法的实验结果吻合得很好     

离心泵叶轮内部固液两相流动的大涡模拟

为提高离心泵的抗磨蚀能力 ,对泵内固液两相流的流态进行深入研究。在多相紊流运动双流体模型的基础之上 ,建立了两相紊流运动的大涡模拟。利用大涡模拟 ,对离心泵叶轮内部的固液两相流动进行数值计算。计算结果得出液相压力分布和相对速度场 ,及固相的颗粒数分布和相对速度场。其中压力和颗粒的分布与同一离心泵在日本所做合作研究量测的实验结果进行比较 ,两者相互吻合     

气固流动大涡模拟和两相湍流模型的评价

用基于气体Smagorinsky亚网格应力模型和颗粒动力论模型的双流体大涡模拟(LES)以及统一二阶矩两相湍流模型(USM-RANS),对轴对称突扩气固流动进行了数值模拟。结果表明:LES模拟的颗粒轴向平均速度和均方根脉动速度,以及USM-RANS计算得到的模拟结果都与实验结果吻合较好。LES得到的颗粒-气体轴向脉动速度关联值比USM-RANS模拟值更接近实验值,这表明USM-RANS模拟还有改进余地,并基本上得到了LES的验证。LES的瞬态模拟能显示气固湍流流动的各向异性两相湍流结构和颗粒弥散的发展过程,而USM-RANS无法实现。     

城垛形喷口圆湍射流拟序结构的流动显示

用流动显示的方法研究了城垛形喷口对圆湍射流拟序结构的影响。实验中喷口流动Reynolds数为5000,Strouhal数为0.50。实验结果表明:城垛形喷口能够影响圆湍射流中大涡结构的产生和演化,进而控制湍射流。无动态激励时,该喷口对射流中拟序结构的影响与垛口数目和垛口高度有关。有动态激励时,城垛形喷口和声波扰动一起,在促进大涡结构的产生和演化方面,总是相互加强的。因此,城垛形喷口是控制圆湍射流拟序结构的一个有效途径。     

声激励对圆射流流场结构控制的大涡模拟

为研究声激励对圆射流流场结构的控制作用,采用大涡模拟方法计算相锁定全局声激励下的圆射流(Re=2 020)流场.计算得出的未加激励时射流的优势频率与实验符合得很好.从多角度描述声激励对射流速度场和涡量场的影响,分析流场对声激励响应的频率选择特性.通过速度场的平均值、均方根值,概率密度函数,偏度,峰度,以及动量厚度的分布,显示声激励引起速度场和混合特性的改变.涡量和Q准则揭示流场拟序结构的演化,发现激励控制流场的主导涡结构是希尔球涡.研究表明,声激励是流场控制非常准确和有效的手段,当激励频率在优势频率附近时影响尤其明显,很小的能量输入便可以引起流场结构的显著改变.     

湍流横向射流的大涡模拟及其涡结构特性

采用大涡模拟方法数值模拟了流向和展向椭圆喷嘴的湍流横向射流, 重点研究了其中旋涡结构的产生、发展等动力学演化过程. 结果表明文献中所报道的横向射流基本涡结构, 如反向旋转涡对、前缘涡、后缘涡、悬涡、肾涡、反肾涡等等是分别对应于新发现的横向射流中的基本涡结构——起始于喷嘴的三维拉伸涡环的局部结构, 因此, 在湍流横向射流中真正占主导作用的是拉伸、扭曲、沿展向摆动和沿流向扭动的三维涡环. 研究还发现: 涡环的脱落频率比流场信号分析得到的脉动频率小得多.     

可压缩混合层流动的三维大涡模拟

为了加深对复杂流动混合现象的认识,为化工设备设计提供理论基础,用大涡模拟方法对可压缩气体混合层进行了三维数值模拟.计算结果显示了混合层流动的三维瞬态发展过程以及流动过程中涡的产生、生长、破裂以及涡列的形成与各涡的合并过程,揭示了可压混合层三维演变过程中同样具有自相似性,计算结果与实验结果相符.通过对混合层初始状态条件下不同扰动和不同初始速度比的计算表明:随机扰动量对混合层的失稳演变过程影响最大,而大混合层速度比同样可使混合层失稳加剧并提前.     

离心泵叶轮内部清水湍流的动态大涡模拟

提出旋转坐标下的单相液体湍流的二阶双系数动态亚格子应力大涡模拟模型并加以验证。考虑亚格子应力的对称性和量纲一致性 ,在基于应变率张量和旋转率张量的不可约量及 Smagorinsky模型和亥姆霍兹定理的基础上 ,提出亚格子应力应表示为可解的应变率张量和旋转率张量的函数。在贴体坐标系下 ,利用交错网格系统和有限体积法离散湍流控制方程 ,采用 SIMPL EC方法求解方程。水泵叶轮中湍流的流速分布规律和压力分布等计算结果都与试验结果基本吻合 ,从而证明了此模型的可行性 ,以及计算方法和程序的可靠性     

四阶紧致格式有限体积法湍流大涡模拟

为准确预测不可压复杂湍流,提出了一种应用于大涡模拟的高精度有限体积法。该方法在非交错网格上数值求解大涡模拟方程,空间离散采用有限体四阶紧致格式,时间推进用四阶Runge-Kutta法,压力速度耦合应用四阶紧致格式的动量插值。通过直接求解顶盖驱动方腔流动和振荡平板上方的流动,证实了该方法具有近四阶精度;并在此基础上,采用动力Smagorinsky亚格子应力模式,成功地完成了充分发展槽道湍流的大涡模拟计算,所得结果与直接数值模拟结果吻合良好。结果表明,该方法是实现高精度湍流大涡数值模拟的一个有效途径。     

淹没条件下水射流涡旋特性大涡模拟及实验研究

采用大涡模拟方法对淹没条件下水射流的涡量场进行数值模拟,分析流场中涡旋的产生与扩散机制,并通过相同条件下粒子图像测速仪测量射流的涡量场,对模拟结果和方法进行验证。模拟研究泵压和围压对淹没射流涡旋特性的影响。结果表明:在射流流场中,由喷嘴出口产生一系列涡量集中的点涡旋,随着射流的前进涡旋逐渐扩散,卷吸周围介质并传递能量,卷吸范围逐渐扩大,而卷吸能力沿射流轴向呈指数衰减;随着泵压升高,整个流场中涡旋的涡量值明显增大,涡旋扩散长度直线上升;围压对涡量基本没有影响,围压的增加会使涡旋扩散区长度直线下降,减小卷吸作用范围。     

气粒湍流流动中考虑重力时颗粒的扩散

对气固两相后台阶流动,利用大涡模拟方法模拟气相场,轨道法模拟颗粒相运动。计算中气相亚格子模式(SGS)采用了标准的Smagorinsky模式。不考虑重力条件,研究了不同入口滑移条件下颗粒的扩散行为,表明颗粒扩散行为与连续相大涡运动行为存在差异。在此基础上深入研究了考虑重力条件时颗粒在流场中的扩散行为。该文工作为在细观框架下进一步研究气固两相相互作用奠定了基础。     

大涡模拟理论及其应用综述

对大涡模拟(Large Eddy Simulation-LES)理论的发展过程、基本理论及其应用现状和前景进行了系统的论述,并概要地介绍了运用LES理论,数值模拟方柱绕流的流动分离现象和动量射流与线性波相互作用的数值模拟成果,表明了LES理论对紊流细部结构具有较强的数值模拟能力。可望用以解决大量的工程实际问题。     

矩形射流中的流向涡分布特性及作用

为研究矩形射流中的流向涡形态及其作用,对矩形射流流动过程进行了三维大涡模拟。计算结果表明在矩形射流的发展区中流向涡才开始逐渐增强,在射流中远场流向涡和展向涡在涡强度、前后涡间距等方面大致相当。与规则匀称的展向涡相比流向涡比较杂乱,在射流的横截面上分布很不均匀,旋向相反的流向涡比邻存在。分析结果表明通过流向涡的旋转运动中心射流和环境流体互相掺混,流体在横截面上不同流向涡之间的游走极大地强化了射流的混合过程。     

获得空间发展槽道流动大涡模拟入口条件的最不稳定波方法

在空间发展流动的大涡模拟中,构造有效的入口扰动是关键。该文提出了利用Fourier分析获得流动最不稳定波的方法来确定入口扰动。结果表明,与同均方差的G auss分布扰动相比,最不稳定波扰动能够显著缩短流动的起涡距离。同时,要达到同样的起涡位置,G auss分布扰动的均方差是最不稳定扰动的350~500倍。可见,提出的最不稳定波入口扰动能够迅速诱发流动本身的K elv in-Helm holtz不稳定,从而提高模拟空间区域的有效性,减小计算量,节约计算时间。