EMMS-based Meso-Scale Model and Its Application in Simulating Gas-Solid Two-Phase Flows

气固两相流是个典型的非线性非平衡系统，呈现出复杂的时空多尺度结构。有关它的连续介质模拟方法，当前主流观点认为：在网格足够细的情况下，不考虑结构影响的双流体模型（Two-Fluid Model, TFM）也称为一般双流体模型（General Two-Fluid Model, GTFM），能准确预测非均匀流动的行为。而我们认为，必须将结构因素考虑到双流体模型中，建立介尺度结构模型，实现对曳力系数的修正，才能准确预测非均匀两相流动行为。上述争论实际上涉及两相流模拟的基本问题，有必要对它进行明确回答。 针对上述基本问题，我们开展了系统的数值实验（论文第二章）。首先针对与粗网格模拟中单个网格尺寸相当的双周期边界Geldart A类颗粒系统，采用GTFM模拟，发现计算结果随着网格的减小逐渐变化，约在10个颗粒直径尺寸时趋于稳定。将此结果应用于实际的提升管体系，不断细化网格直至网格足够小，发现模拟结果仍与实验值相差甚远；对于Gelart B类颗粒系统，也发现上述类似现象。相比较而言，采用基于EMMS的介尺度模型（EMMS/Matrix模型），针对周期边界系统的模拟表明计算结果与网格基本无关。而模拟实际的提升管系统时，发现即使采用较粗网格（只需保证满足一定的差分精度要求），预测的颗粒通量、轴径向浓度分布、径向速度分布以及介观结构变化都与实验观察吻合较好。基于上述结果可知，单纯细化网格并不可行，对于复杂的非均匀气固两相流动，必须建立合理的介尺度模型，而其中，EMMS/Matrix模型则是目前可达到网格无关的较理想模型。 论文第三章进一步简化了原EMMS/Matrix介尺度模型。基于严格的推导，非均匀结构因子HD的三个自变量（ug, us, εg）被简化成两个变量（uslip, εg）, 由此基础上，提出二步拟合法，对HD曲面进行了拟合。采用拟合公式（3%相对误差）取代原来的矩阵插值法，简化了与双流体模型的耦合。 EMMS/Matrix模型可以成功预测A类和B类颗粒的提升管系统的复杂流动，但却不适用于模拟典型的循环流化床锅炉（Circulating Fluidized Bed Bolier, CFBB），其主要问题在于团聚物方程。论文第四章分析了原团聚物方程的适用范围以及失效原因。通过调整输入能量的定义，使得改进后的团聚物方程也能适用于CFBB系统的计算，而原方程则是它的一个特例。方程预测的团聚物直径在固定通量下随着气速的增加而减小，在某个临界值后，又随着气速的增加而增大，体现了颗粒流体协调（particle-fluid compromise, PFC）和气体控制（fluid dominating, FD）的两种极值趋势的交替。采用新团聚物方程的模拟表明，出口通量、流型状态、各区浓度分布等主要监测数据都与实验值吻合较好。 对于循环流化床中典型的状态多值性行为―噎塞，论文第五章利用双流体模型结合基于EMMS的介尺度模型成功预测了实验室和工业尺度反应器中的A类颗粒流化系统中的这种典型现象。在此基础上，论文利用该模拟方法考察了提升管高度对噎塞的影响。模拟表明，噎塞区域随着提升管高度的增加而逐渐增大，而饱和夹带量和噎塞平台起点则维持不变。 论文第六章总结了本论文获得的主要成果，展望了本文的气固介尺度模型的前景以及进一步开展研究的方向

TiO_2/SiO_2复合薄膜的形成历程与超亲水性关系

考察了溶胶-凝胶法制备TiO2／SiO2复合薄膜过程中，正硅酸乙酯(TEOS)的水解时间与膜的超亲水性的关系，通过XRD分析了SiO2添加量对TiO2锐钛矿晶型的影响，并研究复合薄膜的表面处理对其亲水性的影响。实验发现，当无因次反应时间为0．4时，复合薄膜的亲水性最好；SiO2的引入能够抑制金红石晶型的生长，最佳含量在10～20％

TiO2／SiO2复合薄膜的形成历程与超亲水性关系

考察了溶胶-凝胶法制备TiO2／SiO2复合薄膜过程中，正硅酸乙酯(TEOS)的水解时间与膜的超亲水性的关系，通过XRD分析了SiO2添加量对TiO2锐钛矿晶型的影响，并研究复合薄膜的表面处理对其亲水性的影响。实验发现，当无因次反应时间为0．4时，复合薄膜的亲水性最好；SiO2的引入能够抑制金红石晶型的生长，最佳含量在10～20％

TiO2／SiO2复合薄膜的形成历程与超亲水性关系

考察了溶胶-凝胶法制备TiO2／SiO2复合薄膜过程中，正硅酸乙酯(TEOS)的水解时间与膜的超亲水性的关系，通过XRD分析了SiO2添加量对TiO2锐钛矿晶型的影响，并研究复合薄膜的表面处理对其亲水性的影响。实验发现，当无因次反应时间为0．4时，复合薄膜的亲水性最好；SiO2的引入能够抑制金红石晶型的生长，最佳含量在10～20％

TiO_2/SiO_2复合薄膜的形成历程与超亲水性关系

考察了溶胶-凝胶法制备TiO2／SiO2复合薄膜过程中，正硅酸乙酯(TEOS)的水解时间与膜的超亲水性的关系，通过XRD分析了SiO2添加量对TiO2锐钛矿晶型的影响，并研究复合薄膜的表面处理对其亲水性的影响。实验发现，当无因次反应时间为0．4时，复合薄膜的亲水性最好；SiO2的引入能够抑制金红石晶型的生长，最佳含量在10～20％

基于EMMS介尺度模型的双分散鼓泡流化床的模拟

双分散气固鼓泡流化床中颗粒通常具有不同粒径或密度,导致产生颗粒偏析等现象,影响传递和反应行为。颗粒分离和混合与气泡运动密不可分,其中相间曳力起关键作用。最近Ahmad等提出了一种基于气泡结构的双分散介尺度曳力模型,能成功预测双分散鼓泡流化床的床层膨胀系数。本研究耦合该曳力模型与连续介质方法,模拟了两种不同的双分散鼓泡流化床,通过分析不同流化状态下的气泡运动、颗粒浓度比的轴向分布等参数,进一步检验模型的适用性。研究表明,当双分散颗粒处于完全流化状态时,耦合双分散介尺度曳力模型可合理预测不同颗粒的分离现象;而其处于过渡流化状态时,新曳力模型和传统模型均无法获得合理结果,此时调节固固曳力可改进模拟结果

基于emms介尺度模型的双分散鼓泡流化床的模拟

双分散气固鼓泡流化床中颗粒通常具有不同粒径或密度,导致产生颗粒偏析等现象,影响传递和反应行为。颗粒分离和混合与气泡运动密不可分,其中相间曳力起关键作用。最近Ahmad等提出了一种基于气泡结构的双分散介尺度曳力模型,能成功预测双分散鼓泡流化床的床层膨胀系数。本研究耦合该曳力模型与连续介质方法,模拟了两种不同的双分散鼓泡流化床,通过分析不同流化状态下的气泡运动、颗粒浓度比的轴向分布等参数,进一步检验模型的适用性。研究表明,当双分散颗粒处于完全流化状态时,耦合双分散介尺度曳力模型可合理预测不同颗粒的分离现象;而其处于过渡流化状态时,新曳力模型和传统模型均无法获得合理结果,此时调节固固曳力可改进模拟结果

基于emms介尺度模型的双分散鼓泡流化床的模拟

双分散气固鼓泡流化床中颗粒通常具有不同粒径或密度,导致产生颗粒偏析等现象,影响传递和反应行为。颗粒分离和混合与气泡运动密不可分,其中相间曳力起关键作用。最近Ahmad等提出了一种基于气泡结构的双分散介尺度曳力模型,能成功预测双分散鼓泡流化床的床层膨胀系数。本研究耦合该曳力模型与连续介质方法,模拟了两种不同的双分散鼓泡流化床,通过分析不同流化状态下的气泡运动、颗粒浓度比的轴向分布等参数,进一步检验模型的适用性。研究表明,当双分散颗粒处于完全流化状态时,耦合双分散介尺度曳力模型可合理预测不同颗粒的分离现象;而其处于过渡流化状态时,新曳力模型和传统模型均无法获得合理结果,此时调节固固曳力可改进模拟结果

基于多尺度模型的MIP提升管反应历程数值模拟

MIP（maximizing iso—paraffins process）工艺采用两个反应器串联技术，可有效改善汽油质量。MIP反应器的冷态模拟虽能揭示反应器内的流动行为及几何结构的影响，但无法考虑反应引发的变化。为更准确地揭示该反应器中的油气及颗粒运动行为，尝试了三维瞬态反应模拟。模拟采用双流体模型结合十二集总反应动力学模型，并在相间动量传递模型和传热模型中考虑了多尺度结构的影响，然后与基于均匀分布的传统模型作对比。结果表明，相比于传统模型，多尺度模型能较准确预测二反段内的流动结构、颗粒浓度以及温度分布。在预测产率方面，两种模型所得结果类似，都对油浆和柴油的预测较好，对液化气和干气的预测偏差较大。这说明，仅在动量传递及传热模型中考虑多尺度结构的影响是不够的

基于多尺度模型的MIP提升管反应历程数值模拟

MIP（maximizing iso—paraffins process）工艺采用两个反应器串联技术，可有效改善汽油质量。MIP反应器的冷态模拟虽能揭示反应器内的流动行为及几何结构的影响，但无法考虑反应引发的变化。为更准确地揭示该反应器中的油气及颗粒运动行为，尝试了三维瞬态反应模拟。模拟采用双流体模型结合十二集总反应动力学模型，并在相间动量传递模型和传热模型中考虑了多尺度结构的影响，然后与基于均匀分布的传统模型作对比。结果表明，相比于传统模型，多尺度模型能较准确预测二反段内的流动结构、颗粒浓度以及温度分布。在预测产率方面，两种模型所得结果类似，都对油浆和柴油的预测较好，对液化气和干气的预测偏差较大。这说明，仅在动量传递及传热模型中考虑多尺度结构的影响是不够的