Effects of Adiabatic Exponent on Richtmyer-Meshkov Instability

We present a systematical numerical study of the effects of adiabatic exponent γ on Richtmyer-Meshkov instability (RM1) driven by cylindrical shock waves, based on the γ model for the multi-component problems and numerical simulation with high-order and high-resolution method for compressible Euler equations. The results show that the RM1 of different γ across the interface exhibits different evolution features with the case of single γ. Moreover, the large γ can hold back the development of nonlinear structures, such as spikes and bubbles.     

R-T界面不稳定性及湍流混合的大涡模拟

采用大涡模拟(LES)方法研究R-T(Rayleigh Taylor)不稳定性,显示了不稳定发展的小扰动阶段、变形阶段、规则非线性阶段、不规则非线性阶段和湍流混合阶段界面演化与特征,分析了湍流混合层厚度增长的规律,同时还研究了湍流二阶相关量如动量、质量和能量通量在大尺度和亚格子尺度所占分量,验证了LES方法是处理R-T界面不稳定性和湍流混合的一种有效方法.     

RICHTMYER-MESHKOV不稳定性的数值模拟

采用了自适应的非结构网格和基于有限体积法的二阶Godunov格式,数值模拟了在激波作用下两种不同密度流体的交界面的演化过程.着重讨论了Richtmyer Meshkov不稳定性以及斜压效应在交界面演化过程中的作用,并给出了交界面的扰动增长率.     

不同界面组分分布对Richtmyer-Meshkov不稳定性的影响

基于二维非定常Euler方程,对平面激波与不同界面组分分布下氦气气柱作用过程所引起的Richtmyer-Meshkov不稳定性现象进行了数值模拟,探讨了激波冲击轻质气柱后气柱界面形态的演变及流场波系结构,定量分析了气柱特征尺度(气柱长度、高度和中轴宽度)和气柱体积压缩率随时间变化.此外,结合流场压强、速度、环量和气体混合率,多角度分析了激波驱动界面气体混合的流动机制,获得了不同界面组分分布对界面不稳定性的影响.结果表明,随着气柱界面从完全扩散界面向间断界面的过渡,界面两侧的声反射系数随之增大,使入射激波与气柱界面的作用由常规透射转变为非常规透射,反射激波逐渐加强,透射激波逐渐减弱,使得Richtmyer-Meshkov不稳定性随之增强;同时,界面两侧阿特伍德数的增大,加强了Rayleigh-Taylor不稳定性和Kelvin-Helmholtz不稳定性的发展.此外,界面不稳定性的加强使得流场环量增大,导致气体混合率的增长速率随之升高.     

激波冲击V形界面重气体导致的壁面与旋涡作用及其对湍流混合的影响

基于多组分混合物质量分数模型,采用色散最小耗散可控的高分辨率有限体积方法,数值模拟了弱激波冲击V形空气/SF_6界面后,界面不稳定性生成的旋涡与固体壁面作用问题.激波冲击V形界面之后,因斜压效应诱导涡量沉积在界面附近,形成沿界面规则排列的多个涡对结构.旋涡的诱导作用使界面不断变形和卷起,同时旋涡之间不断发生相互并对,诱导更多更小尺度的旋涡产生.旋涡诱导作用的叠加效应,使界面尖端处的初始涡对向上下壁面发展.随后,涡结构开始与壁面发生复杂的相互作用.旋涡与壁面作用后沿壁面加速,使得物质界面沿壁面伸展,随后,旋涡从壁面回弹,并诱导二次旋涡产生.旋涡与壁面相互作用的过程,能够明显加剧物质混合.本文从物质混合的角度研究了该过程的机理,分析了旋涡与壁面作用对物质混合的影响.     

铝电解槽内湍流流动与界面波动的数值模拟

建立了模拟铝电解槽内电解液和铝熔液湍流流及两层熔液界面波动的数学模型，利用交错网格上的一种基于ＳＩＭＰＬＥ算法的改形格式，计算了２８０ｋＡ电解槽内两层熔液的流动状况与界面形状。     

湍流模型在内流流场数值模拟中的应用

本文将Spalart-Allmaras湍流模型方程推导为守恒形式,并应用于内流湍流流场计算.采用此-方程湍流模型分别对扩压器流场以及NASA 67转子流场进行了数值模拟.结果表明常用在外流计算中的Spalart-Allmaras湍流模型在内流计算中仍具有较高的精度、数值稳定性和效率.     

混合对流的数值模拟

分别采用(-v'2)-f模型、k-ε模型和室内零方程模型对通风空调中典型的气流组织形式混合对流进行数值模拟.研究表明:在通风入口附近存在一个分离涡,对这个分离涡预测的准确程度是整个计算成败的关键因素.将预测结果和实验进行比较,发现这个涡的尺度介于(-v'2)-f模型和零方程模型预测的结果之间,比k-ε模型预测的结果大.虽然从流场上看,零方程模型预测的结果和实验结果吻合很好,但从模型构造的角度看,它不具备普遍的意义.从对温度场的预测看,(-v'2)-f模型预测结果最好,室内零方程模型次之,k-ε模型结果最差.     

颗粒湍流扩散的数值模拟

1引言颗粒湍流扩散是气固两相流中的一个重要的研究课题,近年来,国外有大量文献报道这方面的工作。因为颗粒的湍流扩散的研究的一个基本问题就是湍流的描述,它一方面可以来检验湍流描述的正确性;另一方面,该问题是气固两相流理论中,特别是颗粒群轨道模型中一个难点。描述湍流大体上可分为两种方法。一是直接数值模拟方法（DNS）,J．B．Mclaughlin对管流内颗粒沉降的研究就基于直接数值模拟山,K．D．squires和J．K．Eaton对颗粒扩散的研究也基于DNS［2］,DNS一般只适用于低Reynolds数。另一方面,用随机方法来模拟湍流,该…     

水力旋流器湍流流动的数值模拟

本文采用雷诺应力模型计算了水力旋流器的水相湍流流场,计算结果与实验数据吻合很好。与相关文献中采用修正的κ-ε模型的计算结果比较,本文采用雷诺应力模型的计算结果更接近于实验结果。计算得到了水力旋流器内的流线图、等压线以及零速包络面。     

高温空气燃烧炉内湍流混合特性的数值研究

应用自行研发的三维流动、燃烧、传热和污染物NO_x湍流生成的数值模拟程序,对高温空气燃烧实验模型炉进行了湍流扩散燃烧混合特性的数值模拟.数值预报了燃烧室内气体燃料和空气的混合物分数及其湍流脉动的三维分布.数值研究结果表明:在一定的几何条件和气体动力学条件下,高温空气燃烧的湍流混合在更广泛的区域内以较小梯度的进行;混合物分数的脉动主要分布在燃烧区,这表明高温空气燃烧的火焰厚度更大,具有燃烧释热更趋均匀的特性.数值模拟结果与相关的实验结果有相同的规律.     

Rayleigh-Taylor和Richtmyer-Meshkov不稳定性的FCT方法数值模拟

介绍了二维流体力学不稳定性程序的FCT数值方法,数值模拟Rayleigh-Taylor和Richt-myer-Meshkov流体不稳定性,在线性阶段,与线性理论符合很好;在非线性阶段,与俄罗斯激波管实验的计算结果符合很好。计算结果表明:FCT方法有较高计算精度,给出了不稳定性发展的好的图象,适合于ICF烧蚀和内爆流体不稳定性问题的计算。     

激光湍流大气传输数值模拟中计算参量的选取

根据抽样原则、湍流介质特性、湍流效应和光传播效应等全面分析了网格间距、相屏间距和网格数目等关键计算参量选取的基本要求.通过数值计算,分析各种要求的可行性,表明在具体的数值模拟过程中,并非所有要求都能同时满足,若严格按照要求选择计算参量,所需的计算量常超出计算机负荷范围.结合具体数值计算,综合考虑计算机承受能力和计算准确度,给出了一种可靠实用的计算参量选取思路和方案.     

多模态RM不稳定性对初始扰动条件的依赖性分析北大核心CSCD

利用可压缩多介质黏性流动和湍流大涡模拟程序MVFT,对多次冲击作用下的三维多模态Richtmyer-Meshkov(RM)不稳定性发展及其对初始扰动条件的依赖性进行了数值模拟分析。湍流混合区宽度在初始冲击后以幂次律增长,在反射冲击后和第一次反射稀疏波作用后,以具有不同增长因子的指数规律增长,在第一次反射压缩波作用后近似以线性规律增长;而湍流混合区统计量则以类似的规律衰减。多模态RM不稳定性发展对初始扰动条件有很强的依赖性,主要体现在初始冲击后至反射冲击前和反射冲击后至第一次反射稀疏波作用前这两个阶段,即在第一次反射稀疏波作用后,湍流混合区的发展逐渐失去对初始扰动条件的记忆。     

液体在线混合室混合效果的数值模拟

本文根据液体在线混合室的两种不同结构归纳出两种入口界条件，相应获得二维和三维物理模型，应用有限差分法对其混合效果分别作了数值模拟。研究结果表明：设计液体在线混合室时，为了达到最佳混合效果，应使得浓度梯度与速度梯度方向一致，即让两种液体的入口管在同一根轴线上。     

二阶STC格式及其对跨音速湍流流动的数值模拟

ＳＴＣ格式是一种求解积分形式物理守恒律的计算方法,具有守恒性好,精度高的特点。目前求解Ｅｕｌｅｒ方程的ＳＴＣ格式已经建立起来,但求解Ｎ－Ｓ方程的ＳＴＣ格式尚在发展之中。为了使ＳＴＣ格式能够有效地求解粘性流动的问题,本文构造出了求解曲线坐标系下的二维Ｎ－Ｓ方程的二阶ＳＴＣ格式。用该格式对跨音速湍流流动问题进行了计算。     

粗糙壁面上湍流边界层的直接数值模拟

本文采用多重直接力内嵌边界法计算流体与粗糙壁面之间的相互作用力,直接数值模拟了粗糙壁面上湍流边界层的流动状况。粗糙壁面则通过流向和展向分别叠加正弦函数的方法产生。研究发现,峰值较大的位置优先诱导产生出涡结构,并且加大了涡结构的倾斜角度。此外,较大峰值处的粗糙度所产生出的"喷出"和"扫入"事件更加强烈,其对外层的作用范围也就更大。同时粗糙壁面的峰值大小对于雷诺应力的输运也起着至关重要的作用。     

不同湍流模型对圆射流数值模拟的讨论

采用FLUENT软件,选择Spalart-Allmaras模型,标准κ-ε模型以及RNG κ-ε模型对气体无限空间淹没湍流圆射流进行了数值模拟,并将计算结果与理论分析及实验测量方法所得结果进行比较.结果表明,三种模型都能较好反映射流规律,但计算数值与理论值和实测值有所差异,标准κ-ε模型的模拟效果最好.     

旋转槽道湍流标量场输运特性数值模拟研究

展向旋转槽道湍流中的标量场输运过程是与众多工程流动直接相关的模型问题。基于直接数值模拟工作对该问题开展系统性的研究。由壁面摩擦速度定义的流动Reynolds数固定在180,重点考察Schmidt数和旋转数的影响。结果表明,较弱旋转强度即可对主导流动结构形态产生明显的影响：此时槽道不稳定侧产生流向大尺度结构,由此导致标量场出现条带状结构。强旋转时不稳定侧出现被湍流充分混合的区域,而在稳定侧流动层流化并出现近似传导区。平均标量剖面在湍流区和层流区呈现斜率不同的线性分布。Schmidt数小于1时,湍流区标量场脉动和湍流输运随旋转数出现非单调变化,而Schmidt数大于等于1时两者都随旋转数单调下降。由此导致总标量传输率在Schmidt数小于1时随旋转数先上升后下降,而当Schmidt数大于1时单调下降且在弱旋转区域下降趋势最快。     

非光滑表面湍流减阻的数值模拟研究

本文采用CFD数值模拟的方法对二维光滑平板阻力进行计算,与经典理论公式比较,以验证CFD方法的可靠性.从已有的肋条减限的实验结论出发,对带有肋条结构的平板与光滑平板的阻力进行CFD计算,经对比带肋条平板的阻力减小18.39%.分析数值模拟结果的速度场.涡量场等,初步得到肋条结构减阻的机理.肋条顶部产生的两个方向相反的旋...