氯化镁中试热解炉热工过程数值计算研究.pdf

第 55 卷 第 12 期2023 年 12 月Vol.55 No.12Dec.，2023无机盐工业INORGANIC CHEMICALS INDUSTRY氯化镁中试热解炉热工过程数值计算研究王锦吉1，赵亮1，张孟辉1，徐寒露1，董辉1，2（1.东北大学冶金学院，辽宁沈阳 110819；2.辽宁省流程工业节能与绿色低碳技术工程研究中心，辽宁沈阳 110819）摘要：为了探究MgCl2喷雾热解工艺中的核心设备热解炉内的复杂热过程，提高工艺产品质量和设备使用寿命，采用欧拉-拉格朗日方法建立了气固两相流传热与热解反应耦合的数学模型，研究了热解炉稳定运行时的温度场、颗粒分解过程和HCl分布情况，并探讨了进口气体温度和流量对MgO转化率和HCl分布的影响。结果表明，在研究范围内，随着进口烟气温度和流量的增加，产品转化率逐渐提高。当烟气入口流量为1.5 m3/s、入口温度为1 423 K时和烟气入口温度为1 273 K、入口流量为2.25 m3/s时，产品转化率均达到100%。随着颗粒向下运动，HCl浓度逐渐增大，在出料口处HCl达到最大质量分数为6.5%。此外，随着进口烟气温度和流量的增加，出料口处的HCl浓度呈现先增后减的变化规律。对应HCl浓度最大值时的进口烟气温度和流量分别为1 373 K和1.75 m3/s。关键词：氯化镁；热解；EulerLagrange方法；计算流体力学；数值模拟中图分类号：TQ132.2 文献标识码：A 文章编号：1006-4990（2023）12-0140-06Numerical simulation of thermal characteristics in magnesium chloride pilot pyrolysis furnaceWANG Jinji1，ZHAO Liang1，ZHANG Menghui1，XU Hanlu1，DONG Hui1，2（1.School of Metallurgy，Northeastern University，Shenyang 110819，China；2.Liaoning Provincial Engineering Research Center for Energy Conservation and Green and Lowcarbon Technology of Process Industry，Shenyang 110819，China）Abstract：In order to explore the complex thermal process in the pyrolysis furnace，the core equipment of MgCl2 spray pyrolysis process，and improve the quality of process products and equipment service life，Euler-Lagrange method was used to establish a mathematical model of gasparticle twophase flow heat transfer coupled with pyrolysis reaction.The temperature field，particle decomposition process and HCl distribution during the stable operation of pyrolysis furnace were studied.And the influence of gas inlet temperature and flow rate on MgO conversion and HCl distribution was explored.The results showed that the product conversion rate was gradually increased with the increase of flue gas inlet temperature and inlet flow in the scope of the study.When the flue gas inlet was kept at 1.5 m3/s，the product conversion rate reaches 100%at the inlet temperature of 1 423 K.When the flue gas inlet temperature was kept at 1 273 K，the product conversion rate reaches 100%at the inlet flow rate of 2.25 m3/s.As the particles moving downward，the HCl concentration was radually increased and reached a maximum concentration of 6.5%at the discharge port.Moreover，with the increase of flue gas inlet temperature and inlet flow，the HCl concentration at the discharge port showed a change law of firstly increasing and then decreasing.The flue gas inlet temperature and inlet flow rate corresponding to the maximum value were 1 373 K and 1.75 m3/s，respectively.Key words：MgCl2；pyrolysis；EulerLagrange method；CFD；numerical simulation中国盐湖镁资源储量极其丰富1，但是缺乏有效的利用手段。盐湖中钾、锂等元素被提取利用后，会导致盐湖镁离子的富集，对盐湖自然生态造成破坏，而且也会增加其他元素的提取难度2。为高效 化工装备与设计 引用格式：王锦吉，赵亮，张孟辉，等.氯化镁中试热解炉热工过程数值计算研究 J.无机盐工业，2023，55（12）：140-145.Citation：WANG Jinji，ZHAO Liang，ZHANG Menghui，et al.Numerical simulation of thermal characteristics in magnesium chloride pilot pyrolysis furnace J.Inorganic Chemicals Industry，2023，55（12）：140-145.收稿日期：2023-03-27作者简介：王锦吉（1997），男，硕士，主要研究方向为水氯镁石喷雾热解炉内流动传热特性；E-mail：。通讯作者：董辉（1969），博士，教授，主要研究方向为气流输送床传热与流动特性、烧结余热资源回收与利用；E-mail：。Doi:10.19964/j.issn.1006-4990.2023-0166 1402023 年 12 月王锦吉等：氯化镁中试热解炉热工过程数值计算研究利用盐湖镁资源，同时提高中国高纯MgO产能，学者们开发了先沉淀、后煅烧的“两步法”制备氧化镁的工艺3，然而“两步法”需要大量使用沉淀剂，成本较高。如果采用喷雾热解的方式直接将MgCl2溶液热解成MgO，不仅可以将制备高纯氧化镁成本降低50%以上，而且制备的氧化镁产品粒度更细、品质更高。目前，MgCl2喷雾热解工艺由于存在HCl腐蚀和产品质量不稳定等问题，在国内尚未实现工程化应用，造成这些问题的主要原因在于工艺核心设备热解炉内部的热工过程尚不明确，热解炉内各物理场和组分分布规律缺乏理论依据，因此在生产过程中不能准确地把握热解炉的操作参数。MgCl2喷雾热解炉的研究主要集中在工艺参数的计算和小型的实验上4-6，然而实验手段难以获得热解炉内详细的热工过程，因此也有少部分学者采用简化的挥发模型模拟计算了气固逆流式实验炉的结构参数对炉内流场及温度场分布的影响7。但是未将真实的化学反应动力参数嵌入数值模型，也未对热解炉内颗粒的运动、分解过程进行探究。基于此，本文采用EulerLagrange方法8-9，以某公司设计的原料处理量为200 kg/h的中试热解炉为研究对象，建立了一种嵌入了氯化镁多步反应动力学参数的气体-颗粒流动传热模型。分析了热解炉内的场分布规律及颗粒分解过程，并探究了烟气入口温度及入口流量对产品转化率和HCl分布的影响，为喷雾热解工艺的产业化应用提供理论指导与技术支持。1模型建立1.1物理模型与基本假设热解炉结构如图1a所示，热烟气由顶部沿切向进入热解炉，在炉内呈旋流向下运动，原料（MgCl2溶液）由主体段顶部中心注入，随烟气螺旋向下并不断分解，最终大部分产品会在缩口段底部堆积，定期排出，少部分产品随气流从排气管排出，进入气固分离装置。采用非结构化网格对热解炉进行网格划分，对壁面边界层网格进行局部加密，设计了5种网格划分方案，网格数量分别为63万、91万、114万、182万和213万，监测不同网格划分方案下的出口温度和产品转化率，确定网格数为114万为合适的网格数量，具体网格划分如图1bd所示。根据此热解炉工艺特征，对模型进行如下简化：1）炉内运行工况稳定；2）雾化后的液滴近似为球形颗粒；3）炉壁安装有保温层，散热损失被忽略；4）原料颗粒注射简化为实心锥喷射；5）考虑热解炉内物料颗粒为稀疏相，忽略颗粒与颗粒间的碰撞等相互作用；6）由于颗粒密度远大于气体密度，因此假设颗粒只受曳力和重力作用。1.2数学模型颗粒由注射点注入热解炉内，与从烟气入口而来的高温烟气发生热量交换，颗粒内部的自由水首先受热转变为水蒸气，自由水蒸发完毕后颗粒温度继续升高，剩余的MgCl26H2O受热分解逐步转化为MgO并释放H2O和HCl气体。在整个过程中颗粒与烟气之间不断发生质量、动量、能量和组分交换。将热解炉内烟气作为连续相处理，采用Euler方法描述，颗粒视为离散相，采用Lagrange方法描述，采用牛顿第二定律描述颗粒运动过程。将气体与颗粒之间的质量、动量、能量及组分的交换以源项的形式嵌入控制方程中。选择Realizable k-模型作为湍流模型，因为其在强流线曲率、涡和旋转的模拟上表现较好10，采用离散随机游走模型描述颗粒与湍流涡的相互作用。采用汽化模型和沸腾模型11-12描述颗粒中自由水的相变过程。MgCl26H2O分解速率可以表示为：km=Amexp(-Em/RTp)式中：km为颗粒单位面积上的表面组分m的反应速率，kg/s；Am为指前因子，s-1；Em为反应活化能，kJ/mol；R为通用气体常数，取8.314 J/（molK）；Tp为颗粒温度，K。MgCl26H2O的分解是一个多步反应的过程，将其简化为5步，每一步的反应方程式及动力学参数如表1所示13。a热解炉几何结构；b整体网格划分；c截面网格分布图；d局部网格放大图。图 1热解炉几何结构及网格划分Fig.1Geometric structure and mesh division ofpyrolysis furnace 无机盐工业第 55 卷第 12 期1.3边界条件与模型验证入口边界条件设置为速率入口条件，入口速率为5 m/s，出口边界条件设置为压力出口条件，表压为0。颗粒与烟气进出口和排料口的作用设置为逃逸（Escape）条件，与炉壁的作用设置为反弹（Reflect）条件，具体设置如表2所示。在生产中，为了取得较好的雾化效果，将MgCl2配制成饱和溶液，常温下饱和MgCl2溶液中自由水含量过高，将其升温至353 K，此时MgCl2饱和溶液MgCl26H2O 质量分数为 90%，自由水质量分数为10%。进入热解炉的颗粒初始条件如表3所示。在中试现场进行不同温度下的煅烧实验，测量气体出口处温度和产品MgO的转化率，分别与数值模拟计算结果进行比较，结果如表4所示。由表4可知，误差范围在5%以内，表明模型可信度较高。2结果与分析分析热解炉内温度分布、颗粒热解过程、HCl分布等热工过程；探究烟气入口温度及入口流量对产品转化率和HCl分布的影响；揭示热解炉内发生的复杂的热工过程规律。2.1温度分布热解炉内的温度分布是最典型的热工参数，不仅直接决定了产品质量，还影响着热解炉受到的热应力大小，关系到热解炉的寿命。图2为热解炉内温度分布云图。由图2可知，在炉高911 m处，热解炉横截面上的温度呈现出中间低、炉壁高的规律，随着烟气携带着颗粒继续向下，横截面处的温度迅速趋于均匀，且横截面温度随炉高的降低而降低。在炉高911 m处温度的变化速率最快，说明颗粒在此段反应最为剧烈；在炉高09 m段，温度变化平缓，颗粒的平均出口温度为838.2 K，气体的平均出口温度为887.96 K。热解炉整体温度分布较为均匀，热应力较小。2.2颗粒热分解过程研究热解炉内颗粒的运动及分解过程，有助于从理论层面理解影响产品质量的因素，实现对热解炉操作、结构参数的精确把控。图3为颗粒分解过程。热解炉内颗粒轨迹如图3a所示，着色表示颗粒中MgO的质量分数，用以表示颗粒分解的程度。图表1MgCl26H2O热解动力学参数Table 1Kinetic parameters of MgCl26H2Othermal decomposition编号12345方程式MgCl26H2OMgCl24H2O+2H2OMgCl24H2OMgCl22H2O+2H2OMgCl22H2OMgCl2H2O+H2OMgCl2H2OMgOHCl+HClMgOHClMgO+HCl指前因子Am/s-13.61098.810174.610978.61.2103活化能E/（kJmol-1）66.8138.077.2135.692.2表 2数值模型边界条件设置Table 2Setting of boundary conditions for numerical models边界烟气入口烟气出口排料口 其余炉壁边界类型 Velocity inlet Pressure outlet Wall WallDPM设置EscapeEscapeEscapeReflect热边界类型绝热绝热表 3颗粒注射条件设置Table 3Parameters of injection model参数参数设置温度/K353颗粒注射类型实心锥流量/（kgs-1）0.056颗粒分布Rosin-Rammler分布14最小粒径/m50最大粒径/m100平均粒径/m75w（自由水）/%10喷射角度/（）60表4数值模拟结果与试验结果比较Table 4Comparison of test value and simulation value入口温度/K1 2231 2731 3231 373气体出口温度及误差实验值/K890.83897.50998.641 020.10模拟值/K906.96937.54967.601 002.12误差/%-1.78-4.273.211.79产品转化率/%实验值74.5089.1490.6999.38模拟值77.9087.9794.2597.58误差-4.361.33-3.781.84注：误差=（实验值-模拟值）/模拟值100%。图 2热解炉内温度分布云图Fig.2Cloud diagram of temperature distributionin pyrolysis furnace 1422023 年 12 月王锦吉等：氯化镁中试热解炉热工过程数值计算研究3b显示了不同高度各组分的平均转化率，各组分转化率上升的速率有所不同。从图3b可以看出，前3步反应的产物MgCl24H2O、MgCl22H2O及MgCl2H2O的转化率曲线基本重合，在炉高10 m处时就已经达到了97.34%，而MgOHCl和MgO的转化率分别仅有40.32%和48.31%。从化学反应动力学参数的计算中发现（表1），相同温度下前3步反应的反应速率比后2步反应大6个数量级以上，因此后2步反应是制约转化率的主要因素。在颗粒出口处MgO的转化率达到87.97%，MgOHCl的转化率达到97.72%。出口的产品主要杂质为MgOHCl，与现场实验得到产品检测结果相同。2.3HCl分布由于MgCl2热解会产生腐蚀性的HCl气体，研究热解炉内HCl的分布，可以更好地指导热解炉的防腐工作。热解炉内HCl质量分数分布如图4a所示，由图4a可以看出，由于旋流的作用，除911 m段外，HCl的质量分数在各横截面上的分布均匀。在气流和颗粒下降的过程中，颗粒不断分解，因此HCl的浓度随着炉高的降低而升高，如图4b所示，在排料口处最高，达到了6.5%。因此热解炉的防腐工作应当重点考虑缩口段。2.4烟气入口温度、流量的影响烟气入口温度和流量是生产过程中最直接控制的操作参数，直接影响着产品质量。本文分别探究了烟气入口温度和流量对热解炉稳定运行时的MgO产品转化率及HCl分布的影响。2.4.1烟气入口温度的影响选取烟气入口流量为1.5 m3/s，烟气入口温度为1 223、1 273、1 323、1 373、1 423、1 473 K，探究烟气入口温度对MgO转化率和HCl分布的影响。图5为烟气入口温度与MgO转化率的关系图。由图5可知，MgO的转化率随着烟气温度的升高而升高，且升高的幅度不断减小。这是由于温度升高，颗粒与气体的温差变大，导致两者之间的对流和辐射换热增强，因此颗粒分解更加完全，但是由于投料量有限，随着温度的进一步升高，颗粒分解“后劲不足”，MgO转化率呈现出较强的边际效应。当烟气入口温度由1 223 K升高至1 423 K时，产品转化率逐渐由77.90%升高至100%，增加了22.10%，再提高入口温度转化率基本不变。图6为不同烟气入口温度热解炉内HCl分布。不同入口温度下HCl质量分数随炉高的变化如图6a所示，随着烟气入口温度的升高，除了排料口附近，不同高度的HCl浓度均有所升高。图6b为排料口处HCl质量分数随烟气入口温度变化。由图6b显示，热解炉排料口处HCl浓度随烟气入口温度的升a颗粒轨迹；b不同高度颗粒各组分平均转化率。图3颗粒分解过程Fig.3Decomposition process of particle图5烟气入口温度与MgO转化率的关系图Fig.5Relationship between inlet temperatureof flue gas and conversion rate of MgOaHCl分布云图；bHCl浓度随炉高变化。图4热解炉内HCl分布情况Fig.4Distribution of HCl in pyrolysis furnace 无机盐工业第 55 卷第 12 期高呈现出先增大后减小的趋势，当温度为1 373 K时，达到最大值7.51%。这是由于当入口温度小于1 373 K时，颗粒分解程度随入口温度的升高而增大，导致排料口处的HCl浓度增大；温度进一步升高时，颗粒还未至排料口处就已经分解完全，再向下运动的过程中不会产生新的HCl气体，HCl浓度不再升高，反而使得排料口处HCl浓度呈现出随入口温度的升高而下降的现象。2.4.2烟气入口流量的影响选取烟气入口温度为1 273 K，烟气入口流量为1.25、1.50、1.75、2.00、2.25、2.50 m3/s，探究烟气入口流量对热解炉内MgO转化率和HCl分布的影响。MgO 转化率随烟气入口流量的变化如图7所示。由图7可知，MgO的转化率随着烟气入口流量的增大而增大，且增长率不断变小，这是由于入口流量的增大使得单位时间内进入热解炉的热量增加，单位质量颗粒吸收的热量也随之增加，产品转化率提高。由于投料量有限，MgO转化率随烟气入口流量的增大也呈现出较强的边际效应。烟气入口流量为1.25 m3/s时，MgO的转化率仅为76.84%，随着烟气入口流量增大至2.25 m3/s时，MgO的转化率提高至100%，转化率增加了23.16%，再增大进气量，转化率不再提高。图 8 为不同烟气入口流量下热解炉内 HCl 分布。烟气入口流量的增大，一方面会使得进入热解炉内热量增大，促进颗粒分解，引起HCl浓度增大；另一方面，入口流量的增大会对HCl有稀释作用，会使得 HCl 浓度降低，这就导致了如图 8a 所示的结果，不同进气量下，热解炉横截面HCl平均质量分数随炉高变化线图互相交错，无明显规律。由图8b可以看出，热解炉排料口处HCl浓度随烟气入口流量的增大呈现出先增大后减小的趋势，在入口流量为1.75 m3/s时最大（6.79%）。这是由于进气量较低时，颗粒分解还不完全，随着烟气入口流量的增大，颗粒分解程度增大，导致了排料口处的HCl浓度增大，然aHCl浓度随炉高的变化；b排料口处HCl浓度随烟气入口流量的变化。图8不同烟气入口流量下热解炉内HCl分布Fig.8Mass fraction distribution of HCl in pyrolysisfurnace at different inlet gas flow ratesa不同烟气入口温度下HCl质量分数随炉高的变化；b排料口处HCl质量分数随烟气入口温度的变化。图6不同烟气入口温度热解炉内HCl分布Fig.6Distribution of HCl mass fraction in pyrolysisfurnaces at different flue gas inlet temperatures图7MgO转化率随烟气入口流量的变化Fig.7Changes in MgO conversion rate withflue gas inlet flow rate 1442023 年 12 月王锦吉等：氯化镁中试热解炉热工过程数值计算研究而随着烟气入口流量的进一步增大，颗粒分解趋于完全，烟气对HCl的稀释作用占据了主导地位，使得排料口处的HCl浓度不断降低。3结论采用 Euler-Lagrange 方法对 MgCl2在中试热解炉内流动及传热传质现象进行数值研究，并在模型中嵌入了氯化镁多步反应的动力学参数。探究了烟气入口温度和入口流量对MgO的转化率和HCl分布情况的影响，主要结论如下。1）通过对颗粒分解过程和HCl分布进行分析，结果表明氯化镁发生分解的5步反应中，第13步反应发生迅速，第45步反应进行缓慢，是制约产品转化率的主要原因，出口产品杂质主要是MgOHCl。在颗粒下行过程中HCl不断堆积，在排料口处HCl浓度达到最大。2）提高烟气入口温度和入口流量均对MgO的转化率有正向影响。烟气入口温度由1 223 K升至1 423 K时，MgO的转化率提高了22.10%，烟气入口流量由1.25 m3/s增加至2.25 m3/s时，MgO的转化率提高了23.16%。提高烟气入口温度和入口流量均会使排料口处的HCl浓度呈现出先增大后减小的趋势。参考文献：1 张苏江，张琳，姜爱玲，等.中国盐湖资源开发利用现状与发展建议 J.无机盐工业，2022，54（10）：13-21.ZHANG Sujiang，ZHANG Lin，JIANG Ailing，et al.Current situation and development suggestions of development and utilization of salt lake resources in China J.Inorganic Chemicals Industry，2022，54（10）：13-21.2 蒋晨啸，陈秉伦，张东钰，等.我国盐湖锂资源分离提取进展 J.化工学报，2022，73（2）：481-503.JIANG Chenxiao，CHEN Binglun，ZHANG Dongyu，et al.Progress in isolating lithium resources from China salt lake brine J .CIESC Journal，2022，73（2）：481-503.3 刘家辉，宁志强，谢宏伟，等.氧化镁的制备方法及发展趋势综述 J.有色金属科学与工程，2021，12（3）：12-19.LIU Jiahui，NING Zhiqiang，XIE Hongwei，et al.Preparation methods and development trend of magnesium oxideJ.Nonferrous Metals Science and Engineering，2021，12（3）：12-19.4 王超超.水氯镁石制备高纯氧化镁和碱式硫酸镁水泥研究 D.太原：山西大学，2020.WANG Chaochao.Study on preparation of high purity magnesium oxide and basic magnesium sulfate cement with bischofite D.Taiyuan：Shanxi University，2020.5 刘源滔，刘富舟，杜玮，等.水氯镁石喷雾热解制备氧化镁 J.材料科学与工程学报，2018，36（6）：1010-1015，903.LIU Yuantao，LIU Fuzhou，DU Wei，et al.Preparation of magnesium oxide from bischofite by spray pyrolysis J.Journal of Materials Science and Engineering，2018，36（6）：1010-1015，903.6 都永生，孙庆国，火焱，等.氯化镁喷雾热解制备氧化镁的研究 J.无机盐工业，2013，45（4）：13-14，25.DU Yongsheng，SUN Qingguo，HUO Yan，et al.Study on preparation of magnesium oxide by spray pyrolysis of magnesium chloride J.Inorganic Chemicals Industry，2013，45（4）：13-14，25.7 DU Wei，SUN Ze，LU Guimin，et al.CFD aided design of integrated spray pyrolysis furnace for liquid ore exploitation J.Chemical Engineering and Processing：Process Intensification，2017，122：245-257.8 KARIM M R，BHUIYAN A A，SARHAN A A R，et al.CFD simulation of biomass thermal conversion under air/oxyfuel conditions in a reciprocating grate boiler J.Renewable Energy，2020，146：1416-1428.9 DIBA M F，KARIM M R，NASER J.Numerical modelling of a bubbling fluidized bed combustion：A simplified approachJ.Fuel，2020，277：118170.10 SHAVER D，PODOWSKI M.Modeling of interfacial forces for bubbly flows in subcooled boiling conditions J.Transactions of the American Nuclear Society，2015，113（10）：1368-1371.11 MILLER R S，HARSTAD K，BELLAN J.Evaluation of equilibrium and nonequilibrium evaporation models for manydroplet gas liquid flow simulationsJ.International Journal of Multiphase Flow，1998，24（6）：1025-1055.12 SAZHIN S S.Advanced models of fuel droplet heating and evaporationJ.Progress in Energy and Combustion Science，2006，32（2）：162-214.13 HUANG Qiongzhu，LU Guimin，WANG Jin，et al.Mechanism and kinetics of thermal decomposition of MgCl26H2O J.Metallurgical and Materials Transactions B，2010，41（5）：1059-1066.14 祝杰，刘振华，杨云峰，等.运用RosinRammler函数研究喷淋塔内液滴粒径分布规律 J.高校化学工程学报，2015，29（5）：1059-1064.ZHU Jie，LIU Zhenhua，YANG Yunfeng，et al.Applying the rosinrammler function to study on drop size distribution in spray tower J.Journal of Chemical Engineering of Chinese Universities，2015，29（5）：1059-1064.145