炼镁还原罐内传热及其强化换热分析.pdf

第 卷第 期 年 月西安科技大学学报 张小艳，柯亚萍，刘浪，等 炼镁还原罐内传热及其强化换热分析 西安科技大学学报，（）：，（）：收稿日期：基金项目：国家自然科学基金项目（）通信作者：张小艳，女，陕西西安人，博士，教授，：炼镁还原罐内传热及其强化换热分析张小艳，柯亚萍，刘浪，侯东壮，韩子怡，刘清江 西安科技大学 能源学院，陕西 西安 ；西安科技大学 西部矿井开采及灾害防治教育部重点实验室，陕西 西安 ；西安科技大学 地质与环境学院，陕西 西安 ）摘要：在碳达峰碳中和背景下，皮江法炼镁工艺的高能耗问题不容忽视，而造成该工艺高能耗的原因之一是还原罐内物料传热效率低下。为提高还原罐的传热效率，对还原罐内的传热过程进行数值模拟，分析球团直径、球团初始温度、还原罐内径和肋片参数对还原罐径向传热性能的影响规律。结果表明：球团直径越大、初始温度越高、还原罐内径越小，越有利于热量向还原罐内部传递，其中球团直径的变化对还原罐中心温度变化影响最大，直径为 ，和 的球团依次加热 ，和 后，达到还原反应所需温度（大于 ）；肋片长度越长、宽度越宽、数量越多，对还原罐中心温度的提升效果越明显，其中肋片长度的变化对其影响最大，肋片长度由 以 增量增至 ，加热 时，还原罐中心温度依次增大 、。研究结果为皮江法炼镁还原罐的换热及其强化提供理论参考，也为缩短炼镁周期提供很好的解决思路。关键词：还原罐；温度场；还原程度；强化换热；数值模拟中图分类号：文献标志码：文章编号：（）：开放科学（资源服务）标识码（）：，（，；，；，）：，（：，（），：；引言中国丰富的镁矿储量为镁产业的发展奠定了基础 ，自 年以来，中国的产镁总量占世界生产总量的 ，是全球第一产镁大国 。皮江法炼镁工艺（硅热还原法）是中国最主要的金属镁生产工艺，但该工艺存在还原周期长导致的能耗大问题 ，是单位能耗最高的有色冶金行业之一 。在碳达峰碳中和背景下，为进一步降低皮江法炼镁的能耗，重点应放在优化产业结构及提高能源利用效上 ，而提高能源利用率，应该从皮江法炼镁工艺的传热机理等方面深入研究寻找突破口。目前对皮江法炼镁工艺的研究主要集中在还原反应条件的优化方面。唐祁峰等重点介绍碳热还原炼镁技术的动力学研究进展，指出目前碳热还原法炼镁要实现工业化生产还需解决的问题 。穆晓辉等介绍国内外不同还原剂的热法炼镁研究进展，并分析不同还原剂对炼镁过程的影响 。等研究用铝热法从低品位菱镁矿中直接提取镁的还原过程，结果表明该方法的镁提取率高于常规炼镁工艺 。谢兴以单质镁生产过程产生的熵?作为主要研究手段应用，分析热过程和热设备中有效能损失的动力学原因 。邓军平等通过对水泥中掺加适量镁渣不仅提高了复合水泥的强度，还一定程度上解决了镁渣粉尘污染问题 。等研究还原剂类型和还原剂量对皮江法炼镁工艺的影响 。等开发了一种新型的微波驱动皮江法炼镁工艺 。等采用 压力下耦合电热场制备金属镁，该方法有效地降低了通过常规方法提取镁所需的反应装置温度 。等提出一种在流动氩气中进行的非真空生产镁的新工艺，结果表明新工艺可缩短炼镁 。乔锦华通过对球团中硅铁的不同配比进行试验，分析硅含量对粗镁产量的影响 。大量学者对皮江法炼镁工艺的研究中发现，该工艺存在还原罐内球团温升速度过慢或者存在冷芯使罐内球团反应不充分等问题，从而影响单位时间产镁率和产镁量。为探究其中的原因，李波分析高硅白云石产镁率低的原因，并通试验验证皮江法炼镁工艺的原材料采用高硅白云石是可行的 。等对硅热还原过程的化学动力学机理进行试验研究，建立包含化学反应的三维传热非定常数值方法，并对此进行了验证 。等提出将白云石、硅铁、萤石和粘结剂等粉末材料混合制成球团的一种生产镁的新工艺 。除此之外也有不少学者对炼镁还原罐结构等进行研究。傅大学等为解决皮江法还原罐内传热慢的问题，采用数值方法研究预制球团在还原罐内的传热规律 。等开发利用高温碳化硅颗粒增强还原罐内传热的循环利用新工艺，并对此进行数值模拟，结果表明该技术可显著提高传热效率，缩短生产周期，降低生产成本 。邓军平等设计一种有利于加快物料还原反应速率的新型内热式炼 西 安科技大学学报 年第 卷第 期张小艳，等：炼镁还原罐内传热及其强化换热分析镁反应器，该反应器具有传热速度快，内部温度梯度小的特点 。在对皮江法炼镁工艺进行降耗研究过程中，主要目标是大力提高资源和能源的利用率。目前造成皮江法炼镁工艺高能耗的原因之一是真空热还原阶段还原罐内物料传热效率低下，从而导致还原周期（一般为 ）较长、产镁率低等问题。为缩短炼镁周期，提高还原罐的传热速率，以还原罐为研究对象，基于 仿真模拟软件与传热学基本理论，对还原罐内传热过程进行研究，利用温度和还原反应程度变化，获知各影响因素下还原罐传热性能的变化特性。为进一步降低皮江法炼镁工艺的能耗，缩短还原周期，实现炼镁绿色化、低碳化目标提供理论参考。数值模拟方法 物理模型皮江法炼镁过程是将还原罐内抽成真空，然后在 高温下进行固相反应还原出镁单质，得到粗镁。还原罐内部传热过程如图所示，主要包括：（）罐壁与球团间的辐射；（）球团与球团之间的导热；（）罐壁与球团间的导热；（）球团之间的辐射和球团自身的化学反应吸热 种形式。还原罐为壁厚 ，长度 的圆筒，核桃状球团填满罐内。在大量文献研究中，学者们多将还原罐内球团简化为一个整体，由于还原罐内球团与球团之间存在辐射传热，将球团视作一个整体去建模会引起较大误差，所以将球团简化为球状球团，通过等效体积法可计算得到球团直径。由于还原罐长度较长，直径较小且还原炉内温度均匀，认为热量仅沿还原罐径向传递，故可在轴向取单层球团厚度建立三维物理模型，如图 所示，并且为了后续表述简洁化，模型几何参数尺寸符号见表 。由于传热过程十分复杂，为了便于模拟计算结合实际做出如下假设。）罐壁向球团传热、球团由外向内的传热均可视为一维径向传热。）在整个炼镁周期内，未发生化学反应之前无内热源，达到反应温度后球团持续发生还原反应并忽略镁蒸气对传热过程的影响。）还原罐外侧为定壁温条件。?图 还原罐内部传热过程 表 物理模型参数符号 名称符号还原罐内直径 还原罐内半径 还原罐壁厚 还原罐外径 名称符号球团直径 还原罐肋片长度 还原罐肋片宽度 还原罐肋片数量 个 数学模型 控制方程还原罐外壁向内壁的导热过程为恒壁温条件下无内热源的三维非稳态导热过程，控制方程为（）()()()（）式中为换热时间，；为还原罐的温度，；为还原罐的密度，；为还原罐的比热容，（）；为还原罐导热系数，（）。对于球团的传热过程为三维非稳态有内热源的导热过程，该三维数值模型在极坐标系中的控制方程为（）()()()（）式中为球团的温度，；为球团的密度，图 还原罐物理模型 ；为球团的比热容，（）；为球团导热系数，（）。为由化学反应引起的吸热热源，随后通过用户定义的函数计算并作为源项添加到控制方程中。热源 表达式 为 （）（）式中为生成单位质量镁所需热量，其值取为 ；为单位体积球料的产镁速率，；为单位体积球料的理论产镁量，；为反应转化率。转化率 的微分方程为 ()（）（）式中为比例系数，。边界条件设置一般情况下，初始条件的设置需符合模型的实际条件，边界条件为 提供的壁面条件、内部边界条件和内部单元区域。炉温波动范围不大，为简化计算假设炉温为定值，还原罐模型在传热学中的初始条件及边界条件设置见下式初始条件还原罐初始温度（，）球团初始温度（，）边界条件还原罐外壁温（，）还原罐内壁面以及球团表面都为内壁面，在 中设置为“”耦合面。还原罐外壁面通过导热作用对罐内球团进行加热，球团与球团之间存在导热和辐射换热，加热时，需开启能量方程和辐射模型。（）（）（）（）式中为罐体内壁面和罐内球团表面间的辐射换热量，；为还原罐外壁与内壁热传导换热量，；为斯特藩常量，（）；为壁面辐射换热综合因子；为加热面平均辐射温度，；为球团外表面平均辐射温度，；为加热面或降温面的平均辐射温度，；为材料的导热系数，（）；为某点所在处的温度梯度。计算方法及网格划分采用 模拟软件对还原罐传热过程进行数值模拟，采用 软件对还原罐及球团（还原罐外径为 ，内径为 ，球团直径 ）进行建模，完成建模后利用 的 模块进行网格划分（图 ），对罐内球团进行布尔减操作，画出流体域部分，整个模型的网格划分采用结构性四面体和六面体网格，为计算精度要求对罐内球团部分进行加密处理。同时，在还原罐传热过程的数值模拟中，采用压力求解器和非稳态计算法，开启能量方程和 辐射模型，为提高收敛速度与计算精度，采用有限元体积法，压力耦合器采用 ，松弛因子保持不变。网格无关性验证在数值模拟中，网格数量对计算精度有一定影响，为了提高计算精度，针对还原罐内 的温升过程，分别进行了网格无关性验证及步长独立性验证，以确保模拟结果的准确性。图为罐中心温度随时间变化。从图 可以看出，网格数为 时比 的温度略微升高，和 对比几乎没有区别，所以选择网格数为 进行数值模拟计算；时间步长为 ，的计算结果基本一致，所以选取 的时间步长。模型验证为验证模拟计算结果和源项 编译的准确性，采用文献 的结果进行验证。所验证还原罐模型的尺寸为内径 ，外径 ，壁厚 ，罐内球团直径 ，罐体和球团的物 西 安科技大学学报 年第 卷?图 还原罐传热物理模型及网格划分 第 期张小艳，等：炼镁还原罐内传热及其强化换热分析性参数及边界条件的设置与文献一致，其中罐体加热温度为 ，球团初始温度为 。根据模拟结果计算还原反应程度与文献模拟结果进行对比，图是数值模拟与文献模拟结果的对比情况。其中还原反应程度定义为已反应球团体积和总球团体积之比。通过对比可知，还原罐内温度模拟结果与文献模拟结果接近，误差在 内，最大误差在 时为 ，验证了还原罐模型和模拟方法的准确性和可靠性。?图 还原罐中心温度随时间变化?图 还原罐中还原反应程度随时间变化 结果与分析 球团直径对还原罐内传热的影响在还原罐外径为 ，内径为 ，球团初始温度为 ，球团直径分别为 ，的情况下，对还原罐内的温度场进行数值模拟。从图 可以看出，不同球团直径时，还原罐中心温度在不同时间段升温速率差别明显，在 内直径为 的球团温升速度最快，在 以后中心温度逐渐趋于稳定。当球团直径由 以 增量依次增加至 时，加热的中心温度由 分别增至 ，增幅分别为 ，和 ，可见随着球团直径的增加，还原罐中心温度明显升高，但增加幅度逐渐减小，其中 增长至 时的增幅最大。从加热时间来看，的球团在加热 后才达到还原反应所要求的温度（大于 ），而 、和 的球团分别需要加热，和，可以看出球团直径越大加热达到化学反应温度的时间越短。在达到反应条件后，球团就会发生还原反应进行吸热，从而影响热量向罐内的传递，还原罐内中心温度的升幅逐渐减小最终趋于平稳。?图 不同球团直径还原罐中心温度随时间变化 从图 可以看出，在相同加热时间下，随着球团直径的增加，温度分布更加均匀。在加热 时，球团直径由 每增加 ，罐内最大温差（罐内 与 处球团的温差）由 分别减小至 ，和 ，可以看出球团直径增加会使罐内最大温差减小，但减幅是逐渐降低的，其中 增长至 时的减幅最大。此外图中可以明显看出，直径为 ，的球团在加热 后的温度分布曲线分别与直径为 ，的球团在加热 后的温度分布曲线十分接近，可见增大球团直径可以有效提高传热速率。?图 不同球团直径还原罐内温度径向分布 球团初始温度对还原罐内传热的影响在还原罐外径为 ，内径为 ，球团直径为 时，球团初始温度分别为 ，和 的情况下对还原罐内的温度场进行数值模拟。从图可以看出，在相同加热条件下，提高球团初始温度，可使还原罐中心温度更快达到反应温度。当球团初始温度为 时，加热 还原罐中心温度才能达到 ，从而使最内部球团开始反应；当球团初始温度提高到 时加热，与初温 加热 的温度进行对比，只是略有升高；当球团初始温度提高到 时，加热只需 就可使罐中心温度可提升至 ；当球团初始温度提高到 时，使罐中心温度升至 仅需加热 。从以上数据可以看出，提高球团初始温度，可以缩短还原罐内部球团达到化学反应所需温度的时间，但对于球团初始温度的提升越大则效果越好，并且在初始温度 范围内对缩短达到化学反应条件时间的作用可以忽略不计，但升至 以上会有明显效果。造成该现象的原因是球团升温过程中在达到反应条件后，球团就会发生还原反应而吸热，则向内部传递的热量会低于未发生化学反应时，只有使球团初始温度越接近化学反应温度，才能更快缩短炼镁周期。图 显示不同球团初始温度下还原罐内的温度分布云图，从图 可以看出，不同初始温度下还原罐内的温度分布规律，在加热 时，初始温度为 的还原罐内高温球团数量明显多于初始温度为 ；在加热 时，初始温度为 的还原罐内依然存在代表 的绿色球团，而初始温度为 的还原罐内最低温度已经达到 ，在加热 时，初始温度为 的还原罐内球团均已达到反应温度，与初始温度为 的对比很明显可以缩短还原罐的加热周期。?图 不同球团初始温度还原罐中心温度随时间变化?还原罐内径对还原罐内传热的影响球团初始温度为 ，球团直径为 ，在还原罐内径分别为 ，和 ，壁厚 的情况下对还原罐内的温度场进行数值模拟。从图 可以看出，还原罐内径的增大虽然增加了还原罐受热面积，但同时因球团装料量的增加使得热阻变大，不利于热量向罐内的传递。在相同加热时长下，还原罐直径越小罐中心温度越高。当还原罐内径由 每增加 ，加热 时，还原罐中心温度由 依次减小至 ，加热 时，还原罐中心温度由 依次减小至 ，。以还原罐中心温度加热到 为节点，从加热时间来看，内径为 ，和 所需要的时间分别为 ，而当还原罐内径增加至 时，加热 仍未达到 。从以上数据可明显看出增加还原罐内径严重影响还原罐内传热。西 安科技大学学报 年第 卷图 不同球团初始温度下还原罐内的温度分布 第 期张小艳，等：炼镁还原罐内传热及其强化换热分析?图 不同直径还原罐中心温度随时间变化 从图 可以看出，在相同加热时间下，随着还原罐直径的增加，还原反应程度逐渐降低，且还原罐内径越大罐中心温度达到反应温度所需的时间越长。加热时，内径为 ，和 的还原罐内还原反应程度分别为 、和 ，在这一时间段内还原罐内球团处于温升阶段，所以随着还原罐内径的增加还原程度变化差异不是很明显；还原罐内径由 每增加 ，加热 时，还原罐内还原反应程度由 依次减小至 、，由此可明显看出增大还原罐内径会严重影响炼镁周期。?图 不同直径还原罐还原反应程度随时间变化 肋片参数对加肋片的还原罐内传热的影响扩大还原罐直径可提高还原罐球团填充量，但会严重影响炼镁周期，若能解决传热慢的问题则能提高单罐产镁量。针对大型还原罐传热慢的问题借鉴文献 提出的适用于竖罐皮江法炼镁的肋片罐结构设计方案，将其中空罐及肋片分布和数量进行了修改，肋片罐结构示意图如图 所示。通过数值模拟的方法研究肋片参数对还原罐内传热的影响规律。与内径 的还原罐对比，内径 罐体的球团装料量提升 倍。选择直径 的还原罐，球团直径为 ，初始温度为 ，分别对含有不同肋长、肋宽以及数量的还原罐内温度场进行数值模拟，运用控制变量法分析肋长、肋宽以及数量对还原罐传热的影响规律。?：为肋片数；为肋片长度；为肋片宽度图 肋片罐结构示意?图 不同肋长还原罐中心温度随时间变化 从图 可以看出，在同一加热时长下，加肋片的罐体中心温度均大于无肋片罐，且肋长越长罐中心温度越高，肋长等量增加对应中心温度也近似等量升高。与无肋片罐相比，肋宽为 、肋长分别为 ，的肋片罐，加热 时，还原罐中心温度由 依次增加至 ，增量分别为 ，。随着加热时间的延长，各肋长的还原罐中心温度差值逐渐减小，加热 时，还原罐中心温度依次为 ，增量均小于 ，分析认为，随着罐内温度的升高，与热源的温差逐渐减小，传热速率下降，罐内温度也逐渐趋于平缓。从图 可以看出，同一加热时长下，有肋片的还原罐罐中心温度均高于无肋片还原罐，且肋片越宽，罐体中心温度越高。与无肋片罐相比，肋长为 ，宽度分别为 ，的肋片罐，加热 时，还原罐中心温度由 依次增至 ，增幅分别为 和 ，增幅略有增大。随着加热时间的延长，各肋宽的还原罐中心温度差值逐渐减小，在 时，差值均小于 。?图 肋长不同肋宽还原罐中心温度随时间变化 从图 可以看出，在同一加热时长下，肋片数越多，罐体中心温度越高，且随着加热时间的延长，不同数量的肋片罐中心温度差值逐渐减小。肋片数量由 每增加时，加热 时的中心温度由 分别增至 ，增幅分别为 和 ，可见随着肋片数量的增加，还原罐中心温度明显升高，但增加幅度略微减小，其中 增长至 时的增幅最大。?图 不同肋片数下还原罐中心温度随时间变化 虽然肋片数的增加可有效提高罐中心温度，但与增加肋长和肋宽相比，肋片数的增加对罐内球团的填充量有很大影响。图 为不同肋片数时还原罐内的温度分布云图，从球团填充量来看，肋片数越多，罐内球团填充量则越少，在肋片数 时，罐内球团可以最大程度上充满罐体，在 时，球团数量减小，在 时，球团数量最少。从温度分布来看，在同一加热时长下，肋片数越多则罐内温度分布越均匀。在加热 ，肋片数 时，处在 的低温球团数量（蓝色区域）明显多于肋片数 和 ，随着加热时间的延长，这一规律保持不变，均为肋片数少的罐中心温度较低的球团数量多于肋片数较多的。在加热 时，肋片数 的还原罐内温度大于 的高温球团（红色区域）数量明显少于肋片数 和 ，尤其是 时，除了肋片罐中心的球团外其余球团温度均处于 以上。?图 不同肋片数的还原罐内温度分布（，）（，）各影响因素下还原罐传热性能综合评估表 为各影响因素下加热 （炼镁周期内温升阶段）时还原罐中心温度的变化情况，罐中心温度变化值及变化率均为一个炼镁周期（）计算所得。从表 可以看出，球团直径越大、初始温度越高、还原罐内径越小，越有利于热量向罐体内部的传递。在球团直径以 增量由 增加至 的过程中，加热 时，还原罐中心温度依次升高 、。在球团初始温度由 升高至 的过程中，由于 至 范围内数据差异不大，所以直接对比了 升高至 的结果。当球团初始温度由 依次增至 ，加热时，还原罐中心温度依次增大了 、西 安科技大学学报 年第 卷第 期张小艳，等：炼镁还原罐内传热及其强化换热分析 。在还原罐直径由 增加至 的过程中，加热 时，直径每增加 ，还原罐中心温度依次减小 、。表 各影响因素下还原罐中心温度变化（）（）影响因素变化范围罐中心温度 变化值 变化率 变化值 变化率 变化值 变化率 球团直径 球团初始温度 罐体内径 图 更直观地显示出球团直径、球团初始温度与还原罐中心温度呈正相关，还原罐内径与还原罐中心温度呈负相关，并且随着加热时间的延长，还原罐中心温度随球团直径、球团初始温度、还原罐内径变化而变化的幅度逐渐减小。其中对还原罐中心温度变化影响最大的是球团直径的变化，其次是还原罐内径，最后是球团初始温度。此处需对此进行说明，图中的球团初始温度影响因素的数据并非等量递增的，由于球团初始温度 的差值较大，所以使球团初始温度影响因素下 的数据高于球团直径影响因素下 处的数据。表 为肋片罐各影响因素下加热 时的还原罐中心温度变化情况，从表 可以看出，肋片长度越长、宽度越宽、数量越多，越有利于热量向罐体内部的传递。在肋片长度由 增加至 的过程中，每增加 ，加热 时，还原罐中心温度依次增大 、。在肋片宽度由 增加至 的过程中，每增加 ，加热 时，还原罐中心温度依次增大 、。在肋片数由 个增加至 个的过程中，每增加 个，加热 时，还原罐中心温度依次增大 、。图 更直观地显示出肋片长度、肋片宽度以及肋片数量与还原罐中心温度呈正相关，并且同表 各影响因素下肋片罐中心温度变化（）（）影响因素变化范围罐中心温度 变化值 变化率 变化值 变化率 变化值 变化率 肋片长度 肋片宽度 肋片数量 个?图 各影响因素下还原罐中心温度变化 样呈现出随着加热时间的延长，还原罐中心温度随肋片长度、宽度、数量变化而变化的幅度逐渐减小。其中肋片长度的变化对还原罐中心温度的影响最大，其次是肋片数量，肋片宽度对还原罐中心温度的影响最小。?图 肋片罐的各影响因素下还原罐中心温度变化 优化分析研究的目的是提高还原罐的传热效率，缩短炼镁周期从而提高炼镁厂的经济效益。通过分析结果可以看出，增大球团直径、提高球团初始温度、减小还原罐内径以及还原罐内增生肋片的方式均能有效提高还原罐传热效率，在实际工业生产的设备改造过程中，还要考虑经济性。在现有条件下，要提高单罐产镁量，首先应该选择增大罐体直径，在大直径的还原罐内壁增设肋片并选择大直径的球团进行还原反应，原因是在提高了单罐产镁量的基础上，增大还原罐传热效率，从而使炼镁周期不变甚至缩短，这大大提高了能源利用率。结论）当球团直径增加至 时，加热 可达到还原反应所需温度，此外将球团初始温度提升至 时，加热 可达到 ，表明提高球团直径和初始温度，可以使还原罐中心更快达到反应温度。）还原罐内径越大还原反应程度越低，达到完全反应时间也越长。还原罐内径由 增加至 ，在加热 时，还原反应程度由 减小为 。通过在还原罐内壁加肋片可以有效解决大型还原罐传热慢的问题，加肋片的罐体中心温度在同一加热时长下均大于无肋片的，肋片长度越长、宽度越宽、数量越多，还原罐内传热效率越高。）综合分析表明球团直径的变化对还原罐中心温度变化影响最大，其次是还原罐直径，球团初始温度对还原罐中心温度影响最小。还原罐内壁加肋片有利于热量向罐中心传递，肋片长度的变化对还原罐中心温度影响最大，其次是肋片数量，肋片宽度对还原罐中心温度影响最小。）在实际工业生产的设备改造过程中，考虑经济性的前提下，要提高单罐产镁量，首先应该选择增大罐体直径，对其内壁增设肋片并选择大直径的球团进行还原反应。参考文献（）：于旭光，邱竹贤 镁工业生产及应用的现状和展望 材料与冶金学报，（）：，（）：，（）：马登民 镁及镁合金产业发展研究 世界有色金属，：，：林如海，孙前 走向新时代的中国镁工业 中国有色金属，（）：，（）：朱祖武，曹黎华 国内皮江法炼镁技术的缺陷与改进途径 南昌高专学报，（）：，（）：，（）：任虎奎，任建勋 我国金属镁冶炼技术现状与发展趋势 陕西煤炭，（）：，（）：，西 安科技大学学报 年第 卷第 期张小艳，等：炼镁还原罐内传热及其强化换热分析 ，（）：尚梅，李强 区域经济发展与能源消费间关系研究 西安科技大学学报，（）：，（）：唐祁峰，高家诚，陈小华，等 真空碳热还原制镁工艺的研究进展 材料导报，（）：，（）：，穆晓辉，罗洪杰，王耀武，等 热法炼镁还原剂研究进展 轻金属，（）：，（）：，（）：谢兴 新法铝热真空还原镁工艺及热力学分析 沈阳：东北大学，：，邓军平，陈新年，郭一萍 镁渣和矿渣对复合水泥性能的影响 西安科技大学学报，（）：，（）：，（）：，（）：，：，（）：乔锦华，梁国军 球团中硅含量对粗镁产量的影响 冶金管理，（）：，（）：李波，褚丙武 高硅白云石皮江法炼镁实验 有色矿冶，（）：，（）：，：，（）：傅大学，张廷安，豆志河，等 球团预制 硅热还原炼镁还原罐内传热 东北大学学报（自然科学版），（）：，（），（）：，（）：邓军平，王晓刚 内热式炼镁工艺实验研究 材料导报，（）：，（）：，杨康定，陈群，任建勋 炼镁还原罐内传热与化学反应的耦合特性 清华大学学报（自然科学版），（）：，（），（）：冯俊小，张志远 镁还原罐内强化换热研究 北京科技大学学报，（）：，（）：（责任编辑：刘洁）