两相平板混合层流动的双向耦合离散涡法数值模拟

提出了一种新的双向耦合离散涡方法,其能方便地将颗粒对流场的作用结合进离散涡模拟中,同时又不破坏离散涡方法固有的优点。采用这种方法对两相平板混合层流动进行了数值模拟。模拟结果表明,在颗粒作用存在的情况下,混合层流动仍能保持良好的自相似性,其动量厚度沿流向仍呈线性发展,说明颗粒不影响大涡结构间配对与合并的随机性。同时,颗粒对混合层流动高速端自由流有阻滞作用,并会使连续相湍流衰减。     

空间发展湍流气粒两相平面混合层的大涡模拟与统计

本文对空间发展的湍流气固两相平面混合层流动进行了大涡模拟研究,其中气相亚网格尺度（SGS）使用结构函数模型,气相控制方程组采用SIMPLE方法求解,固体颗粒运动用拉格朗日方法计算。计算结果正确重现了流体涡结构的卷起、合并和破碎过程,以及小尺寸颗粒在涡边缘（低涡度区）的局部富集现象。对直径分别为42μm、72μm和135μm分别进行了模拟,并将统计结果和实验测量结果（Hishida　et　al[1]）比较,表明两者的平均速度吻合很好,但颗粒数密度和脉动速度存在较明显的差异,因此有必要对亚网格应力和颗粒之间的耦合作用以及拟序结构的三维性对颗粒运动的影响开展深入研究。     

提升管内稠密气粒两相流动大涡模拟

本文采用双流体模型,引入颗粒动力学理论,对提升管内的稠密气粒两相流动进行了大涡模拟。采用改进的分步投影法对滤波后的方程进行显式求解,小尺度量采用Smagorinsky亚格子模式模拟。模拟结果给出的颗粒相速度分布、浓度分布与实验值基本吻合,气固两相存在速度滑移。模拟结果合理预报出了提升管内的环-核流动结构。     

空间发展平板混和层流动的大涡模拟-3D与2D模拟的比较

本文利用大涡模拟方法研究了空间发展平板混合层流动,分别在三维模型和二维模型下,比较了大尺度结构涡的演化过程,以及流向速度、流向脉动速度、横向脉动速度和Reynolds应力的统计时均结果,并与实验结果对照,指出3D模拟无论在物理真实性还是预报结果准确性方面都要明显优于2D模拟。     

乙醇-空气液雾两相流动和燃烧的大涡模拟

采用动态Smagorinsky亚网格应力模型,EBU(ED,eddy dissipation)燃烧模型和欧拉-拉氏两相流模型,对乙醇-空气湍流两相流动和燃烧进行了大涡模拟(LES).瞬态模拟结果与实验结果定性一致.在火焰瞬态结构的高温区域,旋涡强度较大,在高温区边缘附近存在的拟序结构有脱落的趋势.乙醇液滴通过喷雾器喷出...     

混合层流动的格子涡方法数值模拟

采用格子涡方法对二维平板混合层流动进行了数值模拟,重点考察了格子涡方法定量模拟混合层流动流场下游不同流向位置各湍流统计平均量的能力。模拟结果表明:除横向脉动速度均方根比实验值大外,混合层流动在下游不同流向位置的流向平均速度,流向脉动速度均方根和雷诺应力的模拟结果都与实验结果定量符合得很好。此外,模拟结果还很好地反应出混合层流动的自相似性质。从而表明,格子涡方法具有对空间发展混合层流动进行较精确的定量模拟的能力。     

两相混合层中颗粒运动的数值模拟

本文采用离散涡方法对平板混合层流动进行了数值模拟,得到了与实验完全定量符合的速度场。再用单向耦合方法模拟了混合层流场中颗粒的运动。分析了混合层流动中大尺度涡结构及Ｓｔｏｋｅｓ数对颗粒扩散的影响。与前人工作中所采用的每个时间步一个颗粒在固定的位置进入计算域的方法不同,本文中每个时间步有多个颗粒在入口处以随机的横向位置进入计算域。因此,在不需增加太多计算量的基础上,计算域中可以包含足够多的颗粒以获得较精确的统计结果。采用本文方法得到的颗粒速度场与实验结果定量符合得很好。     

气固两相平面混合层的直接数值模拟

本文采用直接数值模拟(DNS)方法研究了气固两相混合层的涡量场和颗粒扩散,分析了湍流混合层中涡结构的卷起和配对过程,并讨论了大尺度涡结构的配对过程对平均速度、雷诺应力的影响。同时分析了混合层中不同Stokes数的颗粒在涡结构的作用下的混合和扩散。结果显示Stokes数为1的颗粒主要分布在流场大尺度涡结构的外边界,而Stokes数为0.01在涡结构的作用下,在流场中充分混合。     

后台阶气固两相流动的大涡模拟

本文基于开源计算流体力学软件OpenFOAM,开发了气固两相流动双向耦合求解器,并对Re=13800、颗粒直径为150μm的玻璃球的后台阶两相流动进行了计算。其中气相采用大涡模拟方法,用盒式过滤函数对流体相控制方程过滤求解。颗粒相采用拉格朗日方法对颗粒进行跟踪。将计算结果和实验数据进行对比,以此来检验颗粒受力计算模型和亚格子气固两相求解器的准确性。     

流速比对颗粒在液固两相平板混合层中沉降的影响

本文采用流动显示的方法对平板混合层中上下层流体速度比对固体颗粒在混合层中沉降的影响进行了研究。实验中分别采用粒径小于 ４０ ｕｍ,粒径 ９８～１０４ｕｍ,粒径 １５４～１６０ｕｍ的玻璃微珠以及环氧树脂作为固相颗粒,对这些颗粒在速度比分别为１：１．２、１：２和１：２．８的液相混合层中的运动进行了显示。结果表明混合层中大涡结构对固体颗粒的沉降具有迟滞作用,其作用程度取决于混合层中上下层流体速度比。速度比越大,颗粒的沉降越慢。     

三维槽道两相流颗粒运动的大涡模拟

利用基于动态粘性亚网格模型的大涡模拟方法,分别模拟了Reτ=180的三维槽道湍流内的三种颗粒,以及Reτ分别为180和395时的塑料颗粒的运动状况。其中气相控制方程的求解采用分步投影方法,泊松方程采用基于FFT的快速求解,颗粒的拉氏控制方程采用二阶龙格-库塔积分．模拟统计得到了流体和颗粒的平均速度、脉动速度等,不同的颗粒在流场中(特别是近壁区)表现出来了不同的行为,并再现了颗粒在壁面附近局部富集的现象．随Reτ提高,颗粒各方向上的速度脉动都有所增强．     

有旋和无旋气固两相湍流的不同结构

了解有旋和无旋突扩气粒两相详细湍流结构对控制燃料-空气混合、火焰稳定以及燃烧污染物生成很重要.周力行等曾经对其中的两相时平均流场和湍流特性进行过雷诺平均的RANS模拟和测量的研究,但是BANS模拟不能给出详细的两相湍流的瞬态结构.本文建立了二阶矩两相亚网格尺度应力模型,对有旋和无旋同轴突扩气粒两相流动进行了双流体大涡模...     

旋流两相流动的DSM—PDF两相湍流模型

本文提出了气相雷诺应力方程模型和颗粒概率密度输运方程模型相结合而构成的ＤＳＭ—ＰＤＦ两相湍流模型,并用该模型和ｋ－ε－Ｋｐ模型模拟了文献中测量的旋流数为０．４７的旋流突扩气粒两相流动．和测量结果的对照表明,两模型都能较好地模拟族流数不高的两相流动平均速度场,但ＤＳＭ—ＰＤＦ模型可以揭示出两相湍流的各向异性,因此有预报强旋两相流动的潜力．     

DSM-LPDF两相湍流模型及旋流两相流动的模拟

本文由流体－颗粒速度的拉氏联合概率密度函数（ＰＤＦ）输运方程出发,用Ｓｉｍｏｎｉｎ建议的Ｌａｎｇｅｖｉｎ模型封闭颗粒所遇到流体瞬时速度的条件期望项,并用Ｍｏｎｔｅ Ｃａｒｌｏ方法直接求解 ＰＤＦ输运方程,将其和求解流体雷诺应力方程模型的有限差分方法结合,建立了雷诺应力－拉氏ＰＤＦ（ＤＳＭ－ＬＰＤＦ,简称ＤＬ）两相湍流模型．用此模型模拟了旋流数为０．４７的突扩旋流气粒两相流动,并与文献中ＰＤＰＡ实验和用类似于单相流动湍流模型封闭方法的时平均统一二阶矩（ＵＳＭ）模型的预报进行了对比．     

流体机械中高雷诺数流动的大涡模拟

本文应用大涡模拟方法建立了管道流在高雷诺数条件下的数值求解方法，推导了用求和下标所表示的三维控制方程在贴体坐标下的张量形式．应用有限差分法对控制方程进行离散求解，用编制的三维湍流场LES计算程序对方弯管内的流动进行了计算验证．     

不同亚格子模式在后台阶湍流流动大涡模拟中的应用

本文用大涡模拟方法研究了湍流后台阶流场中的大涡演变过程,并在此基础上研究了目前大涡模拟中比较常用的六种亚格子模式。在相同的流动几何参数以及计算条件下,给出了不同亚格子模式下湍流流动瞬时压力场以及流场中瞬时粘性大小的分布,尤其给出了流场瞬时演变的大尺度涡结构。将六种亚格子模式从回流区长度、计算时间、计算结果的准确性和流场稳定性等不同角度进行了比较。为选取合适的亚格子模式深入研究湍流耗散机理奠定了基础。     

气固两相圆湍射流颗粒对气相湍流调制的机理分析

对空间发展的气固两相圆湍射流进行了双向耦合的大涡模拟,采用对拖拽力源项扰动进行频谱分析的方法揭示了大小颗粒对气相湍流不同调制规律的物理机制.结果表明:大颗粒的拖拽力源项扰动具有较强的低频特征,其与上游区域以大尺度涡为特征的气相低频特征一致,从而强化了气相湍流,而小颗粒却类似无特征频率的噪声,难以与气相脉动特征发生正向耦合而弱化了气相湍流.     

各向同性湍流中颗粒引起湍流变动的直接模拟

本文通过直接数值模拟对均匀各向同性湍流中颗粒对湍流的变动作用进行了研究.颗粒相的体积分数很小而质量载荷足够大,以至于颗粒之间的相互作用可以忽略不计,而重点考虑颗粒与湍流间能量的交换。颗粒对湍流的反向作用使得湍动能的耗散率增强,以至于湍动能的衰减速率增大.湍动能的衰减速率随颗粒惯性的增大而增大。三维湍动能谱显示,颗粒对湍动能的影响在不同的尺度上是不均匀的。在低波数段,流体带动颗粒,而高波数段则相反.