Tracking simulation of ladle thermal state based on inverse heat conduction problem
-
摘要: 现代化高效连续炼钢的关键是钢水温度控制，钢包热状态是钢包周转过程中钢水温降的重要影响因素，因而烘烤结束时的钢包热状态至关重要.由于钢包内部温度较高，生产场地布局复杂，电子原件的使用寿命较短等原因，传统的实测法和数值模拟等方法都无法实现包衬内部实时接触式测温.针对上述情况，文章运用有限差分法正向求解包衬温度场，再建立非稳态的钢包传热一维数学模型，采用Fluent软件模拟火焰温度场，用传热反问题研究方法以钢包易测量的包壳温度为已知条件，对包衬温度分布进行数学反演，并通过计算机C#语言对210t钢包烘烤过程进行智能化模拟追踪，编写了烘烤过程的温度场实时监控软件，为钢包调度和编制合理的烘烤制度提供了一个切实可行的新方法.Abstract: Temperature control is the key to realize modernized continuous-steelmaking, just like ladle thermal state greatly influence the temperature loss of molten steel in the ladle turnover process, so the ladle thermal state is significant especially in the endpoint of ladle preheating. But the traditional measurement method and numerical simulation can't realize real-time non-contact temperature measuring of ladle lining inside because of the high ladle lining temperature, complex production site layout and short service life of electronic component. To address the problem, using finite difference method to solve the ladle lining temperature field and using Fluent simulate flame temperature field, with the help of established mathematical ladle one-dimensional heat transfer model and adopted inverse heat conduction problem, mathematic inversion method was used to calculate the temperature field of ladle lining inside and implement the visual monitor of ladle thermal state during the preheating procedure by means of computer language, which contained in the real-time monitoring software compiled by c#, provides a feasible new method to managing ladle and making reasonable preheating rules.
-
Keywords:
- ladle thermal state /
- model /
- inverse heat conduction problem /
- visual monitor
-
随着核反应堆服役时间不断增长，功率不断增大，堆中包壳和结构材料用奥氏体不锈钢的辐照肿胀问题日益凸显[1-5].辐照肿胀是由材料内空位聚集形成空洞所引起，会使材料体积发生膨胀，造成紧固件断裂，危害反应堆安全.现有的抑制肿胀方法有冷加工和添加合金元素[6-8]，但研究表明，经这2种方法处理的奥氏体不锈钢所能经受的辐照剂量无法超过100 dpa[9-10]，远未达到先进核反应堆对奥氏体不锈钢性能的要求.
核反应堆是相对封闭的高危环境，无法原位观察反应堆中奥氏体不锈钢的辐照肿胀行为[11].奥氏体不锈钢在核反应堆中应用广泛，其服役温度也相对较宽，如压水堆中工作温度约为300 ℃[12]，而快堆中工作温度却在500 ℃左右[13].利用超高压电子显微镜，Sekio等[14]对奥氏体不锈钢在450 ℃下的辐照肿胀行为进行原位观察.考虑到工作温度不同，奥氏体不锈钢在核反应堆中的辐照肿胀行为必然也存在差异，因此，有必要对奥氏体不锈钢在其他温度下的辐照肿胀行为进行系统原位观察和分析.
以SUS316L奥氏体不锈钢为研究对象，原位观察其在300~500 ℃区间内电子辐照过程中的辐照肿胀行为，对比分析不同辐照温度下空洞形核及长大的差异，为预测不同工作温度下奥氏体不锈钢辐照肿胀情况提供实验依据.
1 材料及方法
本研究所用材料为SUS316L奥氏体不锈钢，其化学成分如表 1所示.
表 1 SUS316L奥氏体不锈钢化学成分/（质量分数，%）Table 1. Chemical composition of SUS316L austenitic stainless steel /(mass fraction, %)元素 C Si Mn P S Ni Cr Mo V Al N Fe 含量 0.013 0.20 1.28 0.024 0.001 0 13.32 17.24 2.04 0.04 0.014 0.039 6 余量 利用超高压电子显微镜 (型号为JEOL，JEM-ARM1300) 对直径为3 mm的透射电镜薄膜试样进行电子束辐照.辐照电压1 250 kV，辐照剂量3.6 dpa (30 min)，辐照温度分别为300 ℃、350 ℃、400 ℃、450 ℃和500 ℃.辐照过程中，对辐照区形貌进行原位观察.辐照完成后，利用透射电镜 (Tecnai G2 F20) 观察辐照后空洞分布，并利用Image-pro plus软件对空洞尺寸进行测量统计.
2 结果及讨论
2.1 电子束辐照条件下辐照区形貌变化
电子束辐照过程中，试样内不断引入过饱和点缺陷，即空位和间隙原子[15-16].为趋于能量平衡，这些过饱和点缺陷会不断扩散迁移，聚集形成位错环、空洞等二次缺陷[17-18].以500 ℃电子束辐照过程中原位观察辐照区形貌变化为例，如图 1所示.辐照前由于取向差异，晶界两侧视场衬度不同，但并无点缺陷簇、位错环等缺陷组织 (图 1(a)).辐照5 min后，晶界附近产生大量点缺陷簇，而且越靠近晶界，点缺陷簇的密度越高 (见图 1(b)).受晶界的尾闾作用[19]，辐照引入的过饱和点缺陷不断向晶界扩散，导致晶界附近点缺陷相互聚集几率增加，点缺陷簇密度也随之增大.随着辐照时间增长，形成点缺陷簇和位错环的区域也不断增大.辐照30 min后，整个辐照区都可以明显观察到点缺陷簇 (如图 1(d)).
![]()
图 1 500 ℃电子辐照下辐照区形貌Figure 1. Micrographs of irradiated area under electron irradiation at 500 ℃电子束辐照条件下，不仅有点缺陷簇形成，还可以观察到位错环.利用透射电镜对辐照30 min后的位错环形貌进行观察，如图 2所示.同样以500 ℃为例，可以看出电子束辐照30 min后，试样内形成形态不一的位错环，尺寸也不尽相同，所观察到的最大尺寸将近300 nm.
2.2 不同温度下空洞的演变
除了点缺陷簇及位错环外，辐照引入的过饱和空位还可聚集形成空洞.空洞的产生可分为形核和长大2个过程[20-21].即辐照引入的过饱和空位先相互聚集形成空洞核心，然后空洞核不断吸收空位而长大.
对500 ℃电子辐照条件下空洞的产生及演变进行原位观察，如图 3所示.辐照前，晶界两侧并无空洞 (图 3(a)).当辐照5 min，辐照剂量达到0.6 dpa时，晶界附近已有少量小尺寸空洞形成 (见图 3(b))，空洞在此辐照剂量下已完成形核；随着辐照时间不断增加，空洞数量和尺寸不断增大；辐照30 min时 (辐照剂量为3.6 dpa)，辐照区内清晰可见大量尺寸较大的空洞 (见图 3(d)).
450 ℃、400 ℃、350 ℃及300 ℃电子辐照条件下空洞的演变分别如图 4、图 5、图 6和图 7所示.这4种辐照温度下样品内也都有空洞生成，但空洞的产生时间及尺寸却存在较大差异.450 ℃时，电子辐照5 min后，晶界附近有少量小尺寸空洞形成 (如图 4(b))，表明此时的空洞也已完成形核.但随着辐照温度降低，空位扩散能力下降，导致空洞形核所需时间增加.400 ℃时，电子辐照10 min后晶界附近仍未观察到空洞，15 min后晶界附近出现少量小尺寸空洞；而350 ℃和300 ℃下分别辐照15 min和25 min后才观察到空洞.另外，随着温度降低，空洞尺寸也在逐渐下降.300 ℃时，即使辐照30 min，晶界附近也只能观察到一些尺寸非常小的空洞.
![]()
图 4 450 ℃电子辐照过程中空洞的产生及演变过程Figure 4. Void evolution under electron irradiation at 450 ℃![]()
图 5 400 ℃电子辐照过程中空洞的产生及演变过程Figure 5. Void evolution under electron irradiation at 400 ℃![]()
图 6 350 ℃电子辐照过程中空洞的产生及演变过程Figure 6. Void evolution under electron irradiation at 350 ℃![]()
图 7 300 ℃电子辐照过程中空洞的产生及演变过程Figure 7. Void evolution under electron irradiation at 300 ℃利用透射电镜对辐照30 min后的空洞形貌进行观察，如图 8所示.对比图 8(a)~图 8(d)可知，在300~500 ℃范围内，辐照温度越高，空洞尺寸越大，这与原位观察的结果相符.经统计300 ℃、350 ℃、400 ℃、450 ℃和500 ℃时，空洞的平均尺寸分别为3.6 nm、5.5 nm、7.9 nm、13.3 nm和14.5 nm.
![]()
图 8 不同温度下电子辐照30 min后空洞形貌Figure 8. Void micrographs after 30 min electron irradiation at different temperatures对比300~500 ℃区间内电子辐照过程中空洞的演变情况 (图 3~图 7) 可知，450 ℃和500 ℃时，空洞形核所需时间最短；同时，这2种温度下空洞尺寸也远大于其他温度 (见图 8).空洞是空位在三维空间通过扩散而聚集形成的体缺陷.空位的扩散能力随温度增加而增加，因而在本研究中，高温下 (如≥450 ℃)，空位具有较高扩散能力，易于聚集形成空洞核心并促进空洞长大.低温下 (如≤400 ℃)，空位扩散能力相对较低，因而空洞形核所需时间增长，空洞长大也受到了限制.因此，通过降低奥氏体不锈钢中空位扩散能力，可以抑制空洞长大，减小辐照肿胀.
3 结论
1) 300~500 ℃区间内，3.6 dpa辐照剂量以内，随着电子辐照时间增长，辐照区内点缺陷簇的数量不断增加，尺寸也不断增大.位错环呈现不同形貌，最大尺寸可达300 nm左右.
2) 300~500 ℃区间内，经30 min电子辐照，试样内都有空洞形成.随着辐照温度升高，空洞尺寸不断增大，450 ℃和500 ℃时尺寸最大，分别为13.3 nm和14.5 nm.该研究结果可为不同工作环境下奥氏体不锈钢中辐照肿胀情况的评估提供实验依据.
-
![]()
图 6 测温点3处包衬随时间的温度场
Fig 6. Ladle lining temperature field of temperature measuring points 3
表 1 包衬的物性参数
Table 1 Physical parameters of ladle lining
钢包
结构材料
名称厚度/mm 温度/℃ 导热系数λ
/[W·(m·k)-1]比热C/[J·
(kg·℃)-1]包壁
工作
层尖晶石
砖180 100 0.264 780 200 0.485 862 400 0.915 931 600 1.182 977 800 1.428 1013 包壁
永久
层尖晶石
浇注料100 100 0.995 848 200 1.074 906 400 1.155 978 600 1.251 1024 800 1.281 1064 包壁
包壳钢板 32 / 46.667 502.8 注:“/”表示无数据 
下载: 导出CSV
-
[1] 刘占增, 郭鸿志.钢包传热研究的发展与现状[J].钢铁研究, 2007, 35(1):59-62. http://www.cnki.com.cn/Article/CJFDTOTAL-GTYJ200701017.htm [2] 季乐乐, 贺东风, 徐安军, 等.蓄热式钢包烘烤的数值模拟[J].钢铁, 2013, 48(4):76-81. [3] 吴鹏飞, 徐安军, 贺东风, 等.绝热层对钢包热行为的影响[J].北京科技大学学报, 2012, 5:0-14. http://www.cnki.com.cn/Article/CJFDTOTAL-BJKD201205014.htm [4] 潘学峰, 周海斌, 朱苗勇.连铸中间包内钢液温度变化规律的研究[J].中国冶金, 2007, 17(10):45-49. http://www.cnki.com.cn/Article/CJFDTOTAL-ZGYE200710015.htm [5] FREDMAN T P, TORRKULLA, SAXEN H. Two-dimensional dynamic simulation of the thermal state of ladles[J]. Metallurgical and Materials Transactions B, 1999, 30(2):323-330. doi: 10.1007/s11663-999-0061-2
[6] GLASER B, GÖRNERUP M, SICHEN D. Thermal modelling of the ladle preheating process[J]. steel research international, 2011, 82(12):1425-1434. doi: 10.1002/srin.v82.12
[7] 褚福强, 吴晓敏, 朱毅, 等.氧化铝纤维传热方式及反问题研究[J].工程热物理学报, 2015, 8:029. http://www.cnki.com.cn/Article/CJFDTOTAL-GCRB201508029.htm [8] 王旭东, 臧欣阳, 杜凤鸣, 等.基于实测温度的结晶器瞬态传热反问题计算方法[J].铸造技术, 2014, 7:34-39. http://www.cnki.com.cn/Article/CJFDTOTAL-ZZJS201407034.htm [9] 梅宁, 焦思.基于工程反问题的圆管湍流流体热物性参数反演[J].中国海洋大学学报:自然科学版, 2014, 44(10). http://www.cnki.com.cn/Article/CJFDTOTAL-QDHY201410018.htm [10] 高强权, 王堃, 陈林海.利用共轭梯度法反演钢包内钢液温度分布[J].世界科技研究与发展, 2013, 35(005):578-581. http://www.cnki.com.cn/Article/CJFDTOTAL-SJKF201305004.htm [11] 唐光星.基于传热反问题方法的连铸凝固传热数学建模研究[D].大连理工大学, 2010:45-60. http://cdmd.cnki.com.cn/article/cdmd-10141-2010111120.htm [12] BOZZOLI F, CATTANI L, RAINIERI S, et al.Estimation of local heat transfer coefficient in coiled tubes under inverse heat conduction problem approach[J]. Experimental Thermal and Fluid Science, 2014, 59:246-251. doi: 10.1016/j.expthermflusci.2013.11.024
[13] JIN T, HONG J, ZHENG H, et al. Measurement of boiling heat transfer coefficient in liquid nitrogen bath by inverse heat conduction method[J]. Journal of Zhejiang University SCIENCE A, 2009, 10(5):691-696. doi: 10.1631/jzus.A0820540
[14] LÖHLE S, FUCHS U, DIGEL P, et al. Analysing inverse heat conduction problems by the analysis of the system impulse response[J]. Inverse Problems in Science and Engineering, 2014, 22(2):297-308. doi: 10.1080/17415977.2013.780170
[15] MOVAHEDIAN B, BOROOMAND B. The solution of direct and inverse transient heat conduction problems with layered materials using exponential basis functions[J]. International Journal of Thermal Sciences, 2014, 77:186-198. doi: 10.1016/j.ijthermalsci.2013.10.021
[16] BECK J V, Blackwell B, Clair Jrcrs. Inverse heat conduction Ill-posed problems[M]. James Beck, 1985:88-94. [17] 吴鹏飞.大型预制块钢包热行为过程数值模拟与优化研究[D].北京科技大学, 2012:46-54. -
期刊类型引用(3)
1. 李连奇，杨占兵. 奥氏体不锈钢辐照肿胀和偏析的研究进展. 材料导报. 2021(05): 5122-5129 . 
百度学术
2. 于开元，李彦，孙祎帆，张干，张展瑄. Cu-Nb合金薄膜的微观结构与辐照损伤. 稀有金属. 2021(03): 372-377 . 百度学术
3. 晋克勤. 高校冶金工程实验室安全管理初探. 科技视界. 2019(12): 230-231 . 百度学术
其他类型引用(3)
计量
- 文章访问数: 51
- HTML全文浏览量: 50
- PDF下载量: 3
- 被引次数: 6
