稠密气固两相流双流体大涡数值模拟研究.pdf

Classified Index: O359, TQ051.13 U.D.C.: 621 Dissertation for the Doctoral Degree in Engineering EULER-EULER LARGE EDDY APPROACH AND NUMERICAL SIMULATION OF DENSE GAS-SOLID TWO-PHASE FLOW Candidate ： Yin Lijie Supervisor： Prof. Lu Huilin Academic Degree Applied for ： Doctor of Engineering Speciality： Thermal Energy Engineering Affiliation ： School of Energy Science and Engineering th Date of Defence ： December, 23 , 2009 Degree-Conferring-Institution ： Harbin Institute of Technology 摘 要 摘 要 随着计算机技术和计算方法的不断发展，数值模拟方法已成为研究气固两 相流动的重要方法之一。不断完善数值模拟方法，对深入研究气固两相流动的 复杂机理以及影响因素有重要的意义。 大涡模拟（LES ）同时包含了直接数值模拟（DNS ）和雷诺平均法 （RANS ）的思想，与 DNS 相比，LES 可以用于模拟较高的 Reynold 数和较 复杂的几何结构的流动，与 RANS 相比可获得更多的脉动信息，在湍流的数 值模拟研究中有着重要的发展前景。目前，大涡数值模拟方法在单相流中已有 较成熟的应用，在气固两相流中的应用尚需进一步的研究。 本文采用体积分数加权平均方法，对气固两相流的输运方程进行过滤。考 虑气相不可解尺度湍流和可解尺度湍流之间的动量输运，建立了考虑颗粒相浓 度影响的气相亚格子粘度模型；基于颗粒动理学理论，根据 Smagorinsky 涡粘 模型，考虑过颗粒相不可解尺度湍流和可解尺度湍流之间的动量输运，推导了 颗粒相亚格子粘度模型；采用涡扩散模式，提出了颗粒相亚格子热传导系数方 程，封闭了可解尺度拟颗粒温度方程；考虑了高颗粒浓度下介观尺度的影响， 修正了颗粒相亚格子压力方程。建立了稠密气固两相流双流体大涡数值模拟 （LESg-θ-LESp ）模型。 应用稠密气固两相流双流体大涡数值模拟模型，研究了鼓泡流化床内气固 两相流动特性。模拟观察到了鼓泡流化床中气泡的长大、合并和破碎的运动过 程。模拟结果得到了床内颗粒相时均可解尺度速度和脉动速度等分布，与 Yuu 等实测结果相吻合；得到了床内颗粒相时均浓度分布，与 Taghipour 等实测结 果相吻合。研究发现不考虑颗粒相的亚格子尺度时，在床层中形成了较大的气 泡，随着 Cs 的增加，气泡数量逐渐增加，平均气泡直径逐渐减小，时均可解 尺度拟颗粒温度减小，颗粒浓度在床层中心处逐渐增加而在边壁处降低，时均 滑移速度在床层中心处减小而在边壁处增加；随着颗粒浓度的增加，颗粒相亚 格子尺度粘度和热传导系数减小，而颗粒相亚格子尺度压力增大。 应用稠密气固两相流的双流体大涡模拟模型，研究了提升管内的气固两相 流动特性。通过模拟结果可观察到提升管内的环-核流动结构，得到提升管内 的时均颗粒相浓度、速度等分布，模拟结果与 Knowlton 等实测结果基本吻 合。随着颗粒浓度的增加，可解尺度和亚格子尺度下的时均颗粒相的粘度和热 传导