两相端浮力分层流k—ε模型

将研究精细的两相湍浮力分层流预报模型作为目标,广泛地涉及了两相湍浮力分层流精细模拟了的理论和方法,以通用的ｋ－ε湍流模型为基础,考虑了浮力对紊动的影响,建立了两相湍浮力分层流的ｋ－ε模型,对存在速度和密度分层的湍浮力流动进行了数值模拟,计算结果正确地反映了分层及紊动特征,与实测结果吻合较好。     

利用能量最小多尺度模型预测圆射流的速度分布

将能量最小多尺度(EMMS)模型拓展到圆射流的研究.利用圆射流的稳定性条件,预测圆射流的流速分布.分析表明,湍流圆射流的速度分布指数为1.94.通过与实验结果的对比分析,EMMS模型的平均误差和最大误差分别是-0.69%和-10.57%,低于传统的Gauss分布解、Tollmien解和G?rtler解.     

两相湍浮力分层流k-ε模型

将研究精细的两相湍浮力分层流预报模型作为目标,广泛地涉及了两相湍浮力分层流精细模拟的理论和方法．以通用的ｋ－ε湍流模型为基础,考虑了浮力对紊动的影响,建立了两相湍浮力分层流的ｋ－ε模型,对存在速度和密度分层的湍浮力流动进行了数值模拟．计算结果正确地反映了分层及紊动特征,与实测结果吻合较好     

球形汇聚-反射激波作用下湍动能方程的模化

球形激波与湍流间的相互作用在惯性约束聚变(ICF)、超新星爆发、冲击波碎石等许多工程问题中有重要意义.本文通过直接数值模拟,研究了球形激波的汇聚与反射诱导的湍流及其相互作用问题,分析了这一流动的特点和湍流的产生机制.重点讨论了该流动的湍动能系综平均方程中各项的作用和起主导作用的物理机制,发现传统的剪切湍流的产生机制和黏性项量级普遍偏小,激波才是诱发湍流的主要因素,因而在湍动能方程中起主导作用.由于激波的主导作用,压力脉动对湍动能方程有重要影响,与压力脉动有关的各个项在湍动能方程中均起重要作用.一些在一般湍流中可以忽略的项,例如质量流量项、压力体胀项等,在球形汇聚-反射激波与湍流相互作用问题中起重要作用.本文还探讨了涡黏性模式的有效性,发现满足可实现性原理的涡黏性模式可以较好地模拟雷诺应力.同时,本文揭示了现有湍动能方程封闭模式在本问题的应用中存在很大误差,并尝试对这些项提出新的封闭模式.     

大气边界层湍流积分尺度的分析方法

分析了大气边界层风洞中模拟湍流和实测大气湍流的积分尺度，通过风洞中的模拟湍积分尺度分析，检验了Taylor假设的合理性，并比较了五种常用方法的计算结果，确定出合适的分析方法；利用自相关函数直接积分和自拟合模型（AR模型）的方法分析了实际大气湍流的积分尺度，结果证明自相关函数直接积分的方法最为简便合理。     

风环境数值模拟中模拟植被的数学模型与应用

在建筑风环境数值模拟中,植被的影响通常不可忽略,而且流经植被覆盖区的湍流流动非常复杂.为此,采用在动量和湍流输运方程中附加源/汇项的方法,结合粗糙壁面率修正,模拟植被绕流效应.首先基于剪切应力输运(shear-stress transport,SST)k(湍动能)-ω(耗散频率)模型,从理论上推导了湍动能耗散频率ω方程的新的附加源项模型;其次,采用数值模拟方法,以一植被挡风墙试验为例,对该模型进行了验证;最后,将模型应用到一高层建筑群风环境的数值模拟优化设计中.结果表明,该模型可以较准确有效地模拟建筑风环境中的植被效应,因而具有一定的理论和实用参考价值.     

大差异双分散循环流化床多场耦合数值模拟

为探讨颗粒双分散效应对气固流化体系综合性能的影响,对大差异双分散循环流化床进行了多场耦合数值模拟。采用欧拉双流体模型和能量最小多尺度原理(EMMS)介尺度曳力模型,对Geldart A类与B类颗粒组成的气固体系流动特性进行模拟研究;引入丙烷脱氢(PDH)反应动力学模型,对体系的传热和反应性能进行了相应探究;基于两类颗粒在循环流化床内的快速传热与分离,进一步提出流化催化裂化(FCC)与丙烷脱氢工艺耦合思路,初步论证了其可行性。数值模拟结果显示:因大差异双分散特性影响,两颗粒相在流化稳定时体现为分离态,存在稳定的分离层,B类颗粒处于鼓泡流化状态,A类颗粒处于湍动流化状态;丙烷入口温度为30℃、入口速度为0.1 m·s^-1、催化剂初始床层温度及循环温度为700℃时,该耦合工艺可将催化剂温度降低至669℃,也可进一步实现丙烷到丙烯的良好转化,转化率可达95%以上。研究表明了该耦合工艺的可行性,可在优化提升管内反应条件的同时,进一步获得更高附加值的轻烃油气产品。该研究可为深入理解气固体系内颗粒双分散效应的影响提供一定参考,所提耦合工艺可为FCC工艺优化提供新思路。     

基于格子Boltzmann方程的大涡模拟对湍流时空关联性的研究

将格子Boltzmann方程和大涡模拟(LBE-LES)相结合,提出适应于格子Boltzmann方法(LBM)的涡黏性亚格子尺度模型,开展均匀各向同性湍流时空关联性的研究.采用D3Q19格式计算湍流的三维能谱、湍动能耗散率和其它高阶统计量,与实验和直接数值模拟结果的比较表明,该模型比传统涡黏模型有明显改进.考察了不同亚格子模型预测湍流频率波数能量谱的能力,结果表明,尺度涡产生的横扫作用是造成小尺度涡时间去关联的主要因素,不同波数的频率能量谱之间有一定的相似性,横扫速度是描述湍流频率波数能量谱的特征量.     

分区多流域提升管EMMS曳力模型的构建

考虑CFB提升管快速床中因颗粒团聚物的存在而发生的气固相间曳力下滑的现象,建立基于能量最小多尺度方法(EMMS)的曳力模型,并将其与CFD两相流模型相耦合模拟传统等径提升管和新型分区多流域提升管两种不同结构的提升管。结果表明:将EMMS曳力模型耦合CFD双流体模型,可以合理地预测提升管内轴向与径向上的流场分布,实现对提升管内气固两相流动的准确表达;当CFB提升管的结构与操作工况发生变化时须重新寻求更为合适的曳力模型。     

雷诺时均模拟的两种转捩模型比较

对两种基于RANS方程的转捩模型进行分析讨论,这两个模型都只需采用局部变量计算,其中SST-γ-Reθt模型基于实验经验公式,层流动能模型基于转捩现象。两个模型都被耦合到RANS低雷诺数湍流模型中,并通过不同工况下T3系列平板的数值计算对模型进行评估。结果表明转捩模型较原湍流模型需要消耗更多计算资源;相对层流动能模型,SST-γ-Reθt模型对y+更为敏感,在无压力平板算例中两个模型能够较为准确地预测转捩,只有在湍流度增大时,SST-γ-Reθt模型预测的精度才会下降;在有压力梯度情况下,两个模型预测转捩起始点都较实验值延后,当来流雷诺数较大时层流动能模型预测的转捩长度较实验值偏小;从平板上不同位置(层流区、转捩区和湍流区)湍动能的分布情况能够看出SST-γ-Reθt模型只是在数值上模拟转捩过程,并不考虑转捩内在规律,层流动能模型预测结果与实验值较为吻合。     

基于k-ε模型的剪切流中湍流动能衰减规律研究

分别从理论推导和数值模拟的角度,分析了剪切流中湍流动能的衰减规律.通过简化湍流变量输运方程,采用Matlab求解了使用标准k-ε模型计算剪切流时湍动能在不同高度处的数值解,通过fluent数值模拟,验证了理论数值解,并与均匀流中的数值模拟结果进行对比,分析了在剪切流中入口湍流强度和湍流黏性比、来流速度梯度和输运方程扩散项等对湍流变量衰减的影响规律.结果表明:剪切流中存在湍流动能衰减;入口湍流强度对湍流变量输运方程生成项的影响较大,湍流黏性比对湍流变量输运方程生成项的影响较小;当入口湍流强度相同时,入口湍流黏性比越大,湍动能衰减的越快,来流速度梯度对湍动能衰减的影响越小;当入口湍流黏性比相同时,入口湍流强度越大,湍动能衰减的越慢,来流速度梯度对湍动能衰减的影响越大.     

Reynolds应力模型在热分层流场中的应用及优化

对Reynolds应力模型(RSM)在数值模拟热分层流场方面进行了优化,即在Reynolds应力方程和湍动能耗散方程的基础上加入标量通量方程和温度扰动方程,使原来的7方程成为了11方程。优化后的模型可以通过精确计算湍流通量输运方程来求解加热湍流中重要的热浮力项。采用7方程RSM、11方程RSM,并用k-ε湍流模型作为参考,对平板间非稳定热分层流动和模拟池式快堆堆心上方的腔室内受迫对流和热分层之间的混合流动进行数值模拟和比较分析,并将计算结果与现有实验数据相比较。结论是:优化后的11方程RSM比其他的湍流模型更适合于计算各向异性的热分层流动。     

基于VOF模型的导弹低速入水数值模拟方法

基于商用软件Fluent框架,利用动网格技术和UDF方法实现了回转体低速0°攻角直入水过程的数值模拟。采用水、气2相VOF方法捕捉界面,分别采用湍流与层流建立了导弹入水冲击的数值模拟方法,并与实验获得的照片进行了对比分析。分析结果表明:层流模型下导弹入水冲击的数值模拟和试验结果吻合较好,湍流模型下的数值模拟结果误差较大,因此,在低速条件下,湍流模型不适用,应该用层流模型。     

群体互扰下的风场绕流特性和风压分布分析

研究群体互扰下的风场绕流和风压分布问题.改进了流场模拟大涡模型的亚格子应力的湍动黏度,以改善大涡模型对壁面湍流的描述,在壁面处得到与直接数值模拟(DNS)相近的结果.基于此种改进的大涡湍流模型,在单物体和3个物体绕流情形下,分别选择结构刚性模型和弹性模型开展了流场和风压的数值模拟,得出相应的绕流特性和风压分布并与运用其他湍流模型的计算结果进行对比分析,以验证本文提出的改进技术.同时,计算表明,单物体和3个物体情形下的流场特性和风压分布显著不同.     

基于网格数值耗散修正的圆湍射流模型

基于民用机舱内个性通风口,以圆湍射流作为研究对象,以标准κ-ε湍流模型较精细数值模拟与实验对比结果为出发点,通过对粗网格下数值结果的分析,针对网格尺度增大导致计算中数值耗散相应增大、射流速度沿轴向和径向衰减过快问题,从湍流耗散率方程出发,对影响湍流耗散的生成项和耗散项进行修正,平衡因网格尺度产生流动衰减过快的因素,进而提出了耗散率修正的κ-ε湍流模型.采用该修正模型验证了不同雷诺数和非等温条件下射流的适用性.应用修正的湍流模型可以避免机舱内大尺度空间流动对网格的苛刻要求,明显减少网格总数,同时达到对舱内复杂结构和边界条件下流场特征的正确模拟.     

螺旋桨初生空化湍流的多相流数值模拟

摘要：
同时采用修正剪切应力输运(SST)湍流模型和Baseline雷诺应力模型（RSM）求取了E779A螺旋桨在无空化状态和初生空化状态下的梢涡运动轨迹,分析了涡核最小压力系数、湍动能、轴向速度分量和涡核半径沿运动轨迹的变化,并从模拟得到的梢涡卷曲起始和梢涡涡束的角度阐述了梢涡形成机理.空化模型采用改进Sauer模型,考虑了非凝结性气核质量分数、体积分数和气泡初始半径以及湍流脉动的影响,并针对轻度、中度和重度空化面积进行了可信性校验.当空化数σ>初生空化数σi时,叶梢截面压力系数分布相对不再改变的判定准则来确定.涡核中心位于螺旋线垂向截面上最小压力点,涡核边界由湍流涡频率峰值决定.数值模拟结果表明,RSM模拟梢涡路径较修正SST湍流模型稍长、局部梢涡空化范围略大、叶梢最小压力系数和轴向速度分量要小,涡核湍动能分布更为合理.但两者模拟得到的涡核运动轨迹几乎重合,并且初生空化状态下的涡核运动轨迹、最小压力系数和轴向速度分布均与各自无空化状态下非常接近,表明了初生空化状态判定的正确性和改进数值模型对梢涡运动轨迹模拟的适用性.     

FBM模型在栅中翼型云状空化流动数值计算中的应用

评价滤波器湍流模型(FBM)在栅中翼型计算中的应用.分别采用修正RNG κ-ε湍流模型和FBM湍流模型,以及不同滤波器尺寸的FBM模型对栅中翼型的云状空化流动进行模拟,并与水洞实验结果进行对比.结果表明:在模拟计算中采用合适滤波器尺寸的FBM湍流模型,并且在空化模型中考虑到湍动能对空化流动的影响,可以更准确捕捉栅中翼型空化区尾部涡团交替脱落的非定常细节,升阻力与实验更为接近.     

大曲率结构化表面湍流流场特性对比研究

湍流模型选择对于小尺寸大曲率圆管内软性磨粒流湍流流场模拟计算具有重要影响.针对此问题,提出使用标准k-ε模型与考虑两相磨粒流平均流动中的旋转及旋流流动情况的重整化群(RNG)k-ε模型进行仿真对比研究的方法.基于颗粒动力学理论的欧拉模型,对软性磨粒流的控制方程进行求解.通过使用两种湍流模型,对小尺寸大曲率圆管内湍流流动进行仿真计算,得到了两湍流模型下软性磨粒流流场的速度、压力、湍动能等数值模拟结果.仿真结果对比分析表明:标准k-ε模型模拟得到的速度值与实验值存在一定误差,而RNG k-ε模型仿真的速度值误差相对较小;标准k-ε模型对湍动能的估值计算比RNG k-ε湍流模型高,难以计算有滞止点的流动过程;RNG k-ε湍流模型能更有效地模拟有大曲率带分离的湍流流动,从而证明该模型更适合小尺寸大曲率圆管内软性磨粒流湍流流动的研究.     

低磁雷诺数不可压缩磁流体湍流的非线性涡黏性k-ω封闭模型

在低磁雷诺数近似下,通过逐项模化湍动能输运方程的磁流体源项,提出一个适于任意的磁场布置形式的非线性涡黏性的磁流体湍流k-ω模型.在该模型中,使用Kenjeres-Hanjalic模型和Song的非线性涡黏性模型分别模化脉动电场-速度相关量及各向异性的雷诺正应力,并采用Wilcox低雷诺数k-ω模型对层流化进行估计.通过与直接数值模拟结果和实验结果比较,表明本文模型给出的速度型误差小于5%,摩擦阻力系数误差小于10%,能够满足工程需要.     

关于湍流多尺度统计理论方法的论证

用黑河化音测站的实测风、温脉动资料对一种湍流多尺度统计理论方法进行了分析验证,并应用此方法通过对相关系数、湍流时间尺度、湍谱和扩散参数等湍流特征量的研究,分析了黑河戈壁地区的湍流微结构和扩散特征.结果表明:多尺度分析方法正适合研究具有多尺度特征的大气湍流现象,可以比传统方法提供更多的有关湍流属性的信息,能够更直观地表现出不同尺度湍涡对扩散作用的大小,是发展湍流统计理论的一种有效工具