稠密气固两相流双流体大涡数值模拟研究
随着计算机技术和计算方法的不断发展，数值模拟方法已成为研究气固两相流动的重要方法之一。不断完善数值模拟方法，对深入研究气固两相流动的复杂机理以及影响因素有重要的意义。 大涡模拟（ LES）同时包含了直接数值模拟（ DNS）和雷诺平均法（RANS）的思想，与DNS相比，LES可以用于模拟较高的Reynold数和较复杂的几何结构的流动，与RANS相比可获得更多的脉动信息，在湍流的数值模拟研究中有着重要的发展前景。目前，大涡数值模拟方法在单相流中已有较成熟的应用，在气固两相流中的应用尚需进一步的研究。 本文采用体积分数加权平均方法，对气固两相流的输运方程进行过滤。考虑气相不可解尺度湍流和可解尺度湍流之间的动量输运，建立了考虑颗粒相浓度影响的气相亚格子粘度模型；基于颗粒动理学理论，根据 Smagorinsky涡粘模型，考虑过颗粒相不可解尺度湍流和可解尺度湍流之间的动量输运，推导了颗粒相亚格子粘度模型；采用涡扩散模式，提出了颗粒相亚格子热传导系数方程，封闭了可解尺度拟颗粒温度方程；考虑了高颗粒浓度下介观尺度的影响，修正了颗粒相亚格子压力方程。建立了稠密气固两相流双流体大涡数值模拟（LESg-θ-LESp）模型。 应用稠密气固两相流双流体大涡数值模拟模型，研究了鼓泡流化床内气固两相流动特性。模拟观察到了鼓泡流化床中气泡的长大、合并和破碎的运动过程。模拟结果得到了床内颗粒相时均可解尺度速度和脉动速度等分布，与Yuu等实测结果相吻合；得到了床内颗粒相时均浓度分布，与Taghipour等实测结果相吻合。研究发现不考虑颗粒相的亚格子尺度时，在床层中形成了较大的气泡，随着 Cs的增加，气泡数量逐渐增加，平均气泡直径逐渐减小，时均可解尺度拟颗粒温度减小，颗粒浓度在床层中心处逐渐增加而在边壁处降低，时均滑移速度在床层中心处减小而在边壁处增加；随着颗粒浓度的增加，颗粒相亚格子尺度粘度和热传导系数减小，而颗粒相亚格子尺度压力增大。 应用稠密气固两相流的双流体大涡模拟模型，研究了提升管内的气固两相流动特性。通过模拟结果可观察到提升管内的环-核流动结构，得到提升管内的时均颗粒相浓度、速度等分布，模拟结果与Knowlton等实测结果基本吻合。随着颗粒浓度的增加，可解尺度和亚格子尺度下的时均颗粒相的粘度和热传导系数逐渐增加，可解尺度下的时均颗粒相压力先增加，后减小，亚格子尺度下的时均颗粒相压力增加；随着 Cs的增加，在床层中心区域，时均颗粒相浓度增加，气相和颗粒相可解尺度速度、颗粒脉动速度以及颗粒相可解尺度拟颗粒温度降低，在边壁区域变化趋势相反；有效时均颗粒相粘度和颗粒相热传导系数随 Cs的增加而增加，有效时均颗粒相压力和颗粒湍动能随 Cs的增加而减小；对曳力系数进行修正后时均颗粒相浓度、可解尺度拟颗粒温度、床层压降升高，气相和颗粒相可解尺度速度有所降低，但幅度很小，有效时均颗粒相粘度、颗粒相热传导系数和有效时均颗粒相压力增加。 采用稠密气固两相流双流体大涡模拟模型，研究了化学链燃烧（CLC）反应器内的气固两相流动特性，研究发现，空气反应器和燃料反应器中的颗粒浓度呈现中间稀、边壁浓，下浓上稀的分布趋势，形成局部环-核流动的非均匀结构，其流动结构分别与单个提升管和鼓泡流化床类似。采用异相化学反应模型，研究了无载氧体循环的燃料反应器中CH4和NiO之间的催化燃烧反应，研究发现，随着时间的推移，燃料反应器中的NiO不断被还原为 Ni，Ni、CO2和H2O的含量逐渐增多，但是由于反应器温度的降低和NiO含量的减少，使得化学反应速率逐渐减慢，导致Ni以及CO2和H2O增加的速度减慢，料层中CH4的含量逐渐增多。此外，由于此反应为体积增大的反应，在流化过程中，湍动更加剧烈，由于气泡的携带作用，燃料反应器顶部自由空间区域的CH4含量略高于料层上部，而CO2和H2O的含量低于料层上部。
大涡模拟;数值模拟;气固两相流;双流体模型;异相反应
哈尔滨工业大学
博士
热能工程
陆慧林
2009
中文
O359.1
141
2013-03-20（万方平台首次上网日期，不代表论文的发表时间）