创刊于1987年, 双月刊
主管:

江西理工大学

主办:

江西理工大学
江西省有色金属学会

ISSN:1674-9669
CN:36-1311/TF
CODEN YJKYA9

小方坯碳偏析模拟研究

董志成, 张炯明, 马海涛

董志成, 张炯明, 马海涛. 小方坯碳偏析模拟研究[J]. 有色金属科学与工程, 2019, 10(5): 40-45, 112. DOI: 10.13264/j.cnki.ysjskx.2019.05.007
引用本文: 董志成, 张炯明, 马海涛. 小方坯碳偏析模拟研究[J]. 有色金属科学与工程, 2019, 10(5): 40-45, 112. DOI: 10.13264/j.cnki.ysjskx.2019.05.007
DONG Zhicheng, ZHANG Jiongming, MA Haitao. Simulation of carbon segregation in billet[J]. Nonferrous Metals Science and Engineering, 2019, 10(5): 40-45, 112. DOI: 10.13264/j.cnki.ysjskx.2019.05.007
Citation: DONG Zhicheng, ZHANG Jiongming, MA Haitao. Simulation of carbon segregation in billet[J]. Nonferrous Metals Science and Engineering, 2019, 10(5): 40-45, 112. DOI: 10.13264/j.cnki.ysjskx.2019.05.007

小方坯碳偏析模拟研究

基金项目: 

国家自然科学基金资助项目 51834002

详细信息
    通讯作者:

    董志成(1994-),男,主要从事钢的洁净化与控制方向的研究,E-mail:1563426690@qq.com

  • 中图分类号: TF777.3

Simulation of carbon segregation in billet

  • 摘要: 基于Fluent建立了150 mm×150 mm结晶器的三维模型,模拟计算了结晶器内流场、温度场及溶质分布的变化,并对二冷区宏观偏析进行了模拟.结果发现,结晶器角部传热方式为二维传热,与表面一维传热相比凝固速度较快.结晶器角部钢液存在回流,同时弯月面处钢液也存在小的回流.受回流及凝固的影响,碳元素在结晶器内会重新分配,上部表现为正偏析、回流通道表现为负偏析.并且发现,由于固液相扩散系数的不同,直到凝固终点,铸坯冷却过程中都会存在环形负偏析.
    Abstract: Based on Fluent software, a three-dimensional model of 150 mm×150 mm mold was built. The variation in flow field, temperature field and solute distribution in the mold were calculated and the macro-segregation in the secondary cooling zone was simulated. It was found that the solute at the corner of the mold solidified faster than that on the surface as heat transfer of the former is two-dimensional not one-dimensional. The molten steel at the corner of the mold recirculated, and that at the meniscus of the mold refluxed slightly. Under the influence of recirculation and solidification, carbon was redistributed in the mold so that there appeared positive carbon segregation at the upper part of the mold and negative one in the recirculation channel. It was also found that there was annular negative carbon segregation during the entire cooling of the casting billet as a result of the difference in solid-liquid phase diffusion coefficient.
  • 超声波是一种振动频率较高的声波,其频率范围在20 kHz~10 MHz[1].当超声波的能量足够高(超过该液体的“空化阈”)时,就会产生“超声空化”现象[2].由于声波是一种纵波,所以超声波在液体中传播的过程中会在液体中交替产生压缩相和稀疏相,稀疏相所形成的高负压的条件可以使液体中的气体过饱和从而析出,产生大量的微气泡[3].超声空化就是指液体中的微气核随着声压的变化产生周期性的生长、收缩,直至崩溃的动力学过程[4].

    超声波空化作用的强弱与其声学参数及液体的物理化学性质有关.对于液体来说,空化作用的强度与液体的表面张力、黏滞系数、温度、密度以及其中所溶解的气体的量[5]等性质有关;对于超声波本身来说,其空化作用强度则主要与超声波的频率和功率等参数有关[6-7].此外,超声波空化作用的强度还与液体所处外环境的压力大小存在一定关系[8].

    对于充气体系的液体,由于超声波的空化作用,会在气泡壁上产生交替的压力和拉力,并且空化泡崩溃过程中会在周围的液体产生微射流和湍流,这都有利于将液体中的气泡粉碎成更小的气泡.气泡尺寸变小,上升速度会降低,从而使液体的气含率增加[9].浮选过程中,气含率的增加不但可以提高精矿产率和回收率,还可以增加浮选速率[10-11].

    王凡[12]的研究发现,使用超声波强化煤泥浮选时,浮选柱中的气泡尺寸明显变小,可以增加浮选柱的气含率.邹华生[13]利用自主研发的超声鼓泡塔测定了不同频率的超声场对气泡直径的影响,发现在一定条件的超声场中,气泡直径会显著减小.Ozkan[14]在富集菱镁矿时,发现超声处理可以使浮选泡沫变小,并增加了泡沫的稳定性.

    本文通过使用自主设计的超声波气泡发生装置,用排液收集气体的方法测定了在不同条件下超声空化作用所产生气体的量,使用液位上升法探索了超声波对充气条件下溶液气含率的影响.

    实验装置如图 1所示,超声波发生器的型号为JY98-IIID,超声波变幅杆分为20 kHz和28 kHz 2种,气泡发生器上部为标有刻度的玻璃管,可以记录液位的变化.气泡收集器可视为一个倒置的底部封口的漏斗,颈部标有刻度,可以读出所收集的气体的体积.微孔滤片为一圆形微孔材料片,安装在设备底部,可以将氮气瓶冲入的气体进行分散.实验所用药剂MIBC为化学纯,溶液均使用蒸馏水配制.

    图  1  实验设备装置示意图
    1.超声波发生器 2.超声波变幅杆 3.气泡反应器 4.氮气瓶 5.气泡收集器 6.微孔滤片
    Figure  1.  Schematic diagram of experimental equipment

    实验前先调节好超声波发生器的功率为实验所需功率,将气泡发生器和气泡收集器都充满溶液,并使液面位于气泡收集器底部以上的位置,实验过程中氮气瓶的阀门处于关闭状态,不充气.当打开超声波发生器时,超声波变幅杆的前端会产生大量的微气泡,一部分会在前端随着声波不断震荡,一部分会上升至顶部的气泡收集器中,上升至顶部的气泡会将气泡收集器细颈段的液体排出,从而可以通过标有的刻度读出液面高度,得到进入细颈段的气体的体积.超声波发生器设置为间歇工作制度,单段工作时间为10 s,间歇时间为3 s.

    使用GBX公司生产的界面张力仪测定不同浓度MIBC溶液的表面张力,用来表示超声波空化作用产生的气泡量与溶液表面张力的关系.测试方法如下:

    1) 配置溶液:取一定量的蒸馏水加入起泡剂,使用欧洲之星搅拌器搅拌3 min,搅拌均匀配成所需浓度的起泡剂溶液后,在避光处静置3 h用来测试表面张力.

    2) 仪器的校正:把测试所用的铂金片和测量杯用无水乙醇润洗后,使用酒精灯将铂金片和测量杯点燃.测量杯冷却以后,加入30 mL蒸馏水后用表面张力仪测试其表面张力,如果表面张力仪在71~73 mN/m,则校正值正确,校正完毕.

    3) 表面张力的测量:在室温25 ℃左右的条件下进行测试,将MIBC溶液浓度由低到高排列,每次取30 mL加入测量杯中进行测量,每个样品测量3次.如果测量值相差不大,则取算数平均值作为结果,如果相差较大,则再混匀静置后继续测量.

    实验前需先移除气泡收集器.首先配置一定浓度的MIBC溶液,体积为v,加入气泡反应器中,使液面处于标有刻度的中部玻璃管段,记录此时的液面高度.气泡反应器中部玻璃管段的半径为r,所以截面积为πr2.打开氮气罐,调节充气速度为v,待液面稳定后,记录充气速度为v时的气液界面高度2,打开超声波发生器,待液面稳定后记录此时的气液界面高度h2.

    未启动超声波发生器时,溶液气含率的计算公式为:

    $${{\varepsilon }_{1}}\text{=}\frac{\pi {{r}^{2}}\times \left( {{h}_{2}}-{{h}_{1}} \right)}{V+\pi {{r}^{2}}\times \left( {{h}_{2}}-{{h}_{1}} \right)}$$

    超声波发生器开启条件下,溶液气含率的计算公式为:

    $${{\varepsilon }_{2}}\text{=}\frac{\pi {{r}^{2}}\times \left( {{h}_{3}}-{{h}_{1}} \right)}{V+\pi {{r}^{2}}\times \left( {{h}_{3}}-{{h}_{1}} \right)}$$

    其中,为溶液的气含率;R为玻璃管段的半径;h为气液界面距玻璃管底部的高度;V为加入MIBC溶液的体积.

    实验使用蒸馏水,超声波变幅杆的频率为20 kHz,每个功率条件都收集3个工作段的气体,结果如图 2所示.从结果可以看出,超声波空化作用产生的气泡的速率并不是很快,相比一般的充气式发泡,体积基本可以忽略.也可以得出超声波空化作用产生气体的体积是随着超声波功率的增加而增加的,且基本呈线性关系.

    图  2  超声功率对气泡体积的影响
    Figure  2.  Influence of power on gas volume

    实验使用溶液仍为蒸馏水,超声波变幅杆的频率为20 kHz.调节超声波发生器的功率为1 000 W,不更换气泡反应器内的液体.单段超声时间为10 s,3段为1组.打开超声波发生器,记录随着时间增长,每组超声空化产生气体的体积变化,共记录6组,实验结果如图 3所示.由实验结果可知,随着时间的增长,超声波空化产生气体的体积不断减少,但会逐渐趋于平缓.这是由于随着超声空化的不断进行,水中溶解气体的量不断减小,空化作用产生的气体体积也相应减小.而减小到一定程度后,由于外环境空气溶于水的速度与空化作用产生气体速度不断接近,溶液中溶解气体的量趋于平衡,空化作用产生气体的体积也就趋于一定值.

    图  3  气泡体积随时间的变化
    Figure  3.  Change of gas volume over time

    实验使用不同浓度的MIBC溶液,超声波变幅杆频率为20 kHz,超声波发生器的功率为1 000 W.配制浓度分别为5 mg/L、10 mg/L、15 mg/L、20 mg/L、25 mg/L的MIBC溶液进行实验,先测定不同浓度的MIBC溶液的表面张力,再进行发泡实验,实验结果如图 4所示. 从图 4可以看出,MIBC溶液的表面张力是随着MIBC溶液浓度的升高而逐渐降低的,但在不同MIBC浓度条件下,超声波空化产生的气体体积相近,而且与2.1.1实验中同等条件下使用蒸馏水的结果也基本一致.所以可以得出结论: MIBC浓度的改变,会引起溶液表面张力的改变,MIBC浓度越高,溶液的表面张力就越小,但溶液的表面张力的改变对超声波空化作用产生的气泡的体积没有太大影响.

    图  4  不同MIBC浓度下溶液表面张力和产生气泡体积变化
    Figure  4.  Change of gas volume and surface tension under different concentration of MIBC

    把变幅杆更换为28 kHz的变幅杆,使用蒸馏水进行实验,对比20 kHz和28 kHz 2种不同频率的变幅杆超声空化作用产生气泡体积的差异.由于超声波发生器的限制,28 kHz的变幅杆功率调节区间较低,故只能在200~600 W的条件下进行实验,2种频率的实验结果如图 5所示.在使用28 kHz变幅杆的条件下,气泡体积仍然是随着功率的增大而有显著增长,也基本呈线性关系.对比2种频率的实验结果,发现使用28 kHz的变幅杆,空化作用产生的气泡体积要高于使用20 kHz的变幅杆.这是由于超声波的频率越高,溶液中交替产生正负压的速率越快,空化作用发生的概率就越高,所以产生气泡的体积就越大.

    图  5  28 kHz和20 kHz 2种变幅杆空化产生气体体积
    Figure  5.  Gas volume under different power by using the booster of 28kHz and 20kHz

    实验在MIBC浓度分别为10 mg/L 、20 mg/L 、30 mg/L的条件下进行,使用超声波变幅杆为20 kHz,频率为1 000 W,探索在不同充气速率条件下超声波对溶液气含率的影响,实验结果如图 6图 7图 8所示.由实验结果可以得知,在超声波作用的条件下,溶液的气含率有所增加,但并不是随着充气速率的增大而增大.充气速率在125 m3/h时,超声波对溶液气含率的影响最大,充气速率达到150 m3/h时,超声波对溶液气含率的影响反而有所降低.通过对实验过程的观察,分析得出超声波增大溶液的气含率的原因主要有以下2个方面:

    图  6  MIBC浓度为10 mg/L时,超声波对溶液气含率的影响
    Figure  6.  Influence of ultrasonic on gas hold-up when the concentration is 10 mg/L
    图  7  MIBC浓度为20 mg/L时,超声波对溶液气含率的影响
    Figure  7.  Influence of ultrasonic on gas hold-up when the concentration is 20mg/L
    图  8  MIBC浓度为30 mg/L时,超声波对溶液气含率的影响
    Figure  8.  Influence of ultrasonic on gas hold-up when the concentration is 30mg/L

    1) 超声波可以减小溶液中气泡的粒径[15].在实验过程中可以明显观察到超声波对于气泡有粉碎作用,在打开超声波发生器以后,进入超声场中的气泡会明显变小.并且在充气速率较低的时候,由于气泡的上升速度较慢,超声波作用时间充分,对于气泡的粉碎作用更加明显.充气速率增大时,由于气泡的上升速度很快,所以使得超声波对于上升速率大的气泡作用时间不足,从而减弱了对溶液气含率的影响.

    2) 在超声波作用下,气泡发生器内部溶液的流场会发生变化.由于超声波声场的压力,在超声场中的溶液内部会产生震荡.气泡通过超声场时,会受到溶液流场的干扰,改变原本直线上浮的运动轨迹,这就使得气泡在溶液中的滞留时间增长,上升速度减慢,从而增加了溶液的气含率.对于上升速率较慢的气泡,流场的干扰效果会更明显,所以表现为充气速率过快时,超声波对于溶液的气含率影响会有所减弱.实验过程中,可以明显看出反应器内流场的变化.

    图 9为超声波影响气含率的机制示意图.左为未加超声波的充气条件,右为打开超声波发生器,附加超声场后充气条件的变化.

    图  9  超声波影响气含率机制示意图
    Figure  9.  Schematic diagram of the influence of ultrasonic on gas hold-up

    1) 超声波通过空化作用可以产生气泡,产生气泡的速率跟超声波的功率、频率及溶液中溶解气体的量有关,并且都呈正相关关系.

    2) 超声波空化作用产生气泡的速率与溶液的表面张力关系不大,溶液表面张力的减小,空化作用产生的气泡体积并没有太大变化.

    3) 超声波作用于充气溶液中的气泡,通过对气泡的粉碎作用和对流场的改变,可以使溶液的气含率增加,并且在一定的充气速率下,超声波的强化效果最好.

  • 图  1   几何模型及网格划分

    Fig  1.   Geometric model and meshing

    图  2   铸坯低倍及计算结果对比

    Fig  2.   Slab low organization and comparison of calculation results

    图  3   结晶器内钢液流线及速度云图

    Fig  3.   Flow and velocity cloud diagram of molten steel in the mould

    图  4   结晶器内温度云图及温度曲线

    Fig  4.   Temperature cloud and temperature curve in the mould

    图  5   结晶器中液相分数

    Fig  5.   Liquid fraction fraction in the mould

    图  6   结晶器内C浓度分布曲线

    Fig  6.   C concentration curve in the mould

    图  7   结晶器内C元素分布云图

    Fig  7.   C element distribution map in the mould

    图  8   铸坯横截面C元素分布

    Fig  8.   Distribution of element C in cross section of slab

    表  1   二冷区各段长度及冷却水流量

    Table  1   The length of each section of the cold zone and the cooling water flow

    下载: 导出CSV

    表  2   C元素含量、液相线斜率及分配系数

    Table  2   C element content, liquidus slope and partition coefficient

    下载: 导出CSV

    表  3   固液相中C的扩散系数

    Table  3   Diffusion coefficient of C in solid-liquid phase

    下载: 导出CSV

    表  4   钢种物性参数

    Table  4   Physical properties of steel

    下载: 导出CSV
  • [1] 王泽鹏, 肖鹏程, 朱立光, 等.方坯连铸结晶器铜管温度场分析[J].钢铁, 2018, 53(3):38-43. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gt201803007
    [2] 张炯明, 张立, 杨会涛, 等.板坯结晶器钢水凝固的数值模拟[J].北京科技大学学报, 2004(2):130-134. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=bjkjdxxb200402005
    [3] 毛学, 肖鹏程, 刘增勋, 等.方坯连铸结晶器三维热力耦合分析[J].炼钢, 2019, 35(1):47-53. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=lg201901008
    [4] 张兴中, 任素波, 周超, 等.新型结晶器内流场及温度场数值模拟[J].钢铁, 2015, 50(9):53-57. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gt201509009
    [5] 文艳梅, 韩延申, 刘青, 等.大方坯结晶器内的多物理场模拟优化[J].钢铁, 2018, 53(6):53-60. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gt201806009
    [6] 车晓梅, 陈伟庆, 赵江, 等.水口类型对大方坯连铸结晶器内流场和温度场的影响[J].有色金属科学与工程, 2012, 3(1):32-36. http://ysjskx.paperopen.com/oa/DArticle.aspx?type=view&id=201201008
    [7] 李凤善, 陈伟庆.矩形坯连铸结晶器流场和温度场的数值模拟[J].上海金属, 2016, 38(3):73-78. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=shjs201603015
    [8] 董其鹏.方坯连铸疏松及宏观偏析的模拟研究[D].北京: 北京科技大学, 2018.
    [9] 安航航.高碳钢大方坯凝固机理与宏观偏析关键控制技术研究[D].北京: 北京科技大学, 2018.
    [10] 陈桂彬, 贺东风, 袁飞, 等.基于传热反问题的钢包热状态跟踪模拟[J].有色金属科学与工程, 2016, 7(2):25-31. http://ysjskx.paperopen.com/oa/DArticle.aspx?type=view&id=201602005
    [11] 郭学益, 王松松, 王亲猛, 等.氧气底吹炼铜模拟软件SKSSIM开发与应用[J].有色金属科学与工程, 2017, 8(4):1-6. http://ysjskx.paperopen.com/oa/DArticle.aspx?type=view&id=201704001
    [12] 李茂旺, 张国锋, 安航航.150 mm×150 mm断面小方坯结晶器电磁搅拌磁场与流场耦合的数值模拟研究[J].工业加热, 2018, 47(2):49-53. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gyjr201802013
    [13] 苌真真, 张炯明, 王顺玺, 等.150t钢包自由表面旋涡的数值模拟与湍流模型的选择[J].有色金属科学与工程, 2014, 5(4):37-43. http://ysjskx.paperopen.com/oa/DArticle.aspx?type=view&id=201405021
    [14] 朱苗勇, 娄文涛, 王卫领.炼钢与连铸过程数值模拟研究进展[J].金属学报, 2018, 54(2):131-150. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=jsxb201802002
    [15] 佟志芳, 肖成, 魏战龙.分子动力学模拟及其在冶金炉渣中的应用研究[J].有色金属科学与工程, 2016, 7(3):15-20. http://ysjskx.paperopen.com/oa/DArticle.aspx?type=view&id=2016030003
    [16] 任兵芝, 朱苗勇, 王宏丹, 等.大方坯连铸结晶器电磁搅拌三维电磁场与流场的数值模拟[J].金属学报, 2008(4):507-512. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=jsxb200804023
    [17] 干勇, 仇圣桃.连续铸钢过程数学物理模[M].北京:冶金工业出版社, 2001.
    [18] 汪成义. 180 mm×180 mm方坯连铸结晶器内流场及温度场模拟研究[A]//中国金属学会.第十一届中国钢铁年会论文集——SO2.炼钢与连铸[C], 北京: 中国金属学会, 2017: 10.
图(8)  /  表(4)
计量
  • 文章访问数:  118
  • HTML全文浏览量:  30
  • PDF下载量:  5
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-05-06
  • 发布日期:  2019-10-30
  • 刊出日期:  2019-09-30

目录

/

返回文章
返回