CN117688819B - 一种碳-氧反应作用下炼钢转炉熔池流场仿真方法及仿真系统 - Google Patents
一种碳-氧反应作用下炼钢转炉熔池流场仿真方法及仿真系统 Download PDFInfo
- Publication number
- CN117688819B CN117688819B CN202410139965.2A CN202410139965A CN117688819B CN 117688819 B CN117688819 B CN 117688819B CN 202410139965 A CN202410139965 A CN 202410139965A CN 117688819 B CN117688819 B CN 117688819B
- Authority
- CN
- China
- Prior art keywords
- carbon
- molten pool
- oxygen reaction
- converter
- oxygen
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 131
- CSJDCSCTVDEHRN-UHFFFAOYSA-N methane;molecular oxygen Chemical compound C.O=O CSJDCSCTVDEHRN-UHFFFAOYSA-N 0.000 title claims abstract description 130
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000009628 steelmaking Methods 0.000 title claims abstract description 36
- 238000004088 simulation Methods 0.000 title claims abstract description 28
- 230000000694 effects Effects 0.000 title description 4
- 238000007664 blowing Methods 0.000 claims abstract description 88
- 239000007791 liquid phase Substances 0.000 claims abstract description 41
- 238000003723 Smelting Methods 0.000 claims abstract description 32
- 238000004364 calculation method Methods 0.000 claims abstract description 31
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 23
- 230000009471 action Effects 0.000 claims abstract description 19
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims abstract description 17
- 239000003546 flue gas Substances 0.000 claims abstract description 17
- 230000008878 coupling Effects 0.000 claims abstract description 8
- 238000010168 coupling process Methods 0.000 claims abstract description 8
- 238000005859 coupling reaction Methods 0.000 claims abstract description 8
- 239000012071 phase Substances 0.000 claims description 31
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 29
- 239000001301 oxygen Substances 0.000 claims description 29
- 229910052760 oxygen Inorganic materials 0.000 claims description 29
- XKRFYHLGVUSROY-UHFFFAOYSA-N Argon Chemical compound [Ar] XKRFYHLGVUSROY-UHFFFAOYSA-N 0.000 claims description 28
- 229910000831 Steel Inorganic materials 0.000 claims description 21
- 239000010959 steel Substances 0.000 claims description 21
- 239000007789 gas Substances 0.000 claims description 20
- 239000000243 solution Substances 0.000 claims description 19
- 230000000704 physical effect Effects 0.000 claims description 17
- 238000004519 manufacturing process Methods 0.000 claims description 15
- 239000002184 metal Substances 0.000 claims description 15
- 229910052751 metal Inorganic materials 0.000 claims description 15
- 229910052786 argon Inorganic materials 0.000 claims description 14
- 239000002131 composite material Substances 0.000 claims description 11
- 239000002893 slag Substances 0.000 claims description 11
- 238000012805 post-processing Methods 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 10
- 239000012530 fluid Substances 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 8
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 6
- 229910052799 carbon Inorganic materials 0.000 claims description 6
- 230000008676 import Effects 0.000 claims description 6
- 238000005261 decarburization Methods 0.000 claims description 5
- 230000008034 disappearance Effects 0.000 claims description 5
- 230000008602 contraction Effects 0.000 claims description 3
- 238000004134 energy conservation Methods 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 238000002347 injection Methods 0.000 claims description 2
- 239000007924 injection Substances 0.000 claims description 2
- 239000007788 liquid Substances 0.000 claims description 2
- 238000007670 refining Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 abstract description 4
- 238000009851 ferrous metallurgy Methods 0.000 abstract 1
- 239000010931 gold Substances 0.000 abstract 1
- 229910052737 gold Inorganic materials 0.000 abstract 1
- 230000008859 change Effects 0.000 description 5
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 4
- 238000003756 stirring Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 101001121408 Homo sapiens L-amino-acid oxidase Proteins 0.000 description 2
- 102100026388 L-amino-acid oxidase Human genes 0.000 description 2
- 101100012902 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) FIG2 gene Proteins 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011423 initialization method Methods 0.000 description 2
- 229910052742 iron Inorganic materials 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000009865 steel metallurgy Methods 0.000 description 2
- 101000827703 Homo sapiens Polyphosphoinositide phosphatase Proteins 0.000 description 1
- 102100023591 Polyphosphoinositide phosphatase Human genes 0.000 description 1
- 101100233916 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) KAR5 gene Proteins 0.000 description 1
- 230000003542 behavioural effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000004945 emulsification Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 239000011261 inert gas Substances 0.000 description 1
- 238000011112 process operation Methods 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Fluid Mechanics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Carbon Steel Or Casting Steel Manufacturing (AREA)
Abstract
本发明提供一种碳‑氧反应作用下炼钢转炉熔池流场仿真方法及仿真系统，涉及钢铁冶金领域。包括：获取转炉结构参数，构建转炉三维几何模型；对几何模型进行网格划分；设定模型基本假设，选择计算模型，设定物性参数、边界条件以及求解算法，进行初始化操作和迭代计算，得到复吹条件下熔池流场；获取转炉烟气信息，求解熔池碳‑氧反应比例，计算得到熔池液相部分碳‑氧反应气泡量；构建碳‑氧反应气泡生成点，注入碳‑氧反应气泡，对碳‑氧反应气泡在渣‑金界面行为进行编译；耦合计算得到碳‑氧反应作用下炼钢转炉熔池模拟结果，后处理得到不同吹炼时刻碳‑氧反应作用下熔池流场图。本发明方法能够分析不同冶炼时刻碳‑氧反应作用下的熔池流场。
Description
技术领域
本发明属于钢铁冶金领域，尤其涉及一种碳-氧反应作用下炼钢转炉熔池流场仿真方法及仿真系统。
背景技术
转炉炼钢法是当前世界上最主要的炼钢方法，顶吹供氧和底吹惰性气体是复吹转炉冶炼的主要工艺方式。冶炼过程中驱动熔池搅拌和流动的能量来源分为三个方面：顶吹超音速氧气射流、底吹惰性气泡羽流以及熔池内部碳-氧反应生成的CO气泡群流。其中，顶吹超音速氧气射流有利于熔池气-渣-金的乳化，促进脱磷和脱碳反应，极大地影响冶炼过程各种复杂现象的发生；底吹惰性气泡羽流有利于熔池的搅拌和混匀，提升熔池的动力学条件；碳-氧反应生成的CO气泡群流动在不同冶炼时期对熔池形成不同程度的搅拌与混匀作用。明晰上述三种能量源作用下的转炉熔池流场变化规律对于转炉反应器的设计以及高效稳定地运行至关重要。
现有技术中，炼钢转炉冶炼过程中碳-氧反应呈现出典型的三阶段模式，脱碳速率呈先升高后平稳最后降低的趋势，整个过程碳含量逐渐降低，在不同吹炼阶段，熔池内部碳-氧反应生成的气泡对熔池有着复杂的搅拌影响，此外，顶吹氧枪枪位、流量以及底吹流量会随时间而变化。然而，目前现有的技术大多集中于针对顶吹、底吹以及顶底复合吹炼作用下的熔池流场进行研究，有关碳-氧反应作用下炼钢转炉熔池流场的行为特征尚缺乏研究。
发明内容
为了克服现有技术存在的问题，本发明提供一种碳-氧反应作用下炼钢转炉熔池流场仿真方法及仿真系统，用于克服目前该领域研究中存在的缺陷。
本发明的技术方案如下：
一种碳-氧反应作用下炼钢转炉熔池流场仿真方法，包括如下步骤：
步骤S1：获取实际生产中炼钢转炉的结构参数，进行合理简化，建立转炉三维几何模型；
步骤S2：将所述转炉三维几何模型进行结构化网格划分处理，得到转炉网格模型；
步骤S3：将所述转炉网格模型导入仿真软件（具体可以采用ANSYS Fluent仿真软件），设定模型求解的基本假设，选择计算模型，设定物性参数、设定边界条件以及设定求解算法，随后进行初始化操作，最后进行迭代计算，得到顶吹和底吹复合吹炼条件下的转炉熔池流场；
步骤S4：获取转炉实际生产过程中的烟气信息数据，根据碳守恒计算得到不同冶炼时刻熔池气相和液相部分碳-氧反应气泡总量，计算得到转炉不同冶炼时刻熔池液相部分碳-氧反应比例（具体可以通过FactSage软件计算得到），得到熔池液相部分碳-氧反应气泡量；
步骤S5：建立熔池内部碳-氧反应气泡生成点，将所述熔池液相部分碳-氧反应气泡量设定为碳-氧反应的输入流量，对碳-氧反应气泡在渣-金界面消失行为进行编译（具体可以调用UDF程序进行编译），得到熔池内部碳-氧反应气泡行为；
步骤S6：将所述熔池内部碳-氧反应气泡行为与步骤S3得到的顶吹和底吹复合吹炼条件下的转炉熔池流场进行耦合计算，求解得到顶吹、底吹以及碳-氧反应作用下的转炉熔池模拟结果；
步骤S7：对所述碳-氧反应作用下的转炉熔池模拟结果进行后处理操作，得到不同吹炼时刻下碳-氧反应作用下的炼钢转炉熔池流场。
进一步地，所述步骤S2中，将所述转炉三维几何模型进行结构化网格划分处理的过程中，对氧枪拉瓦尔喷头区域、底部元件区域以及气-渣-金界面区域进行局部网格加密处理，所述加密处理具体为：计算域采用六面体结构化网格划分，整体网格质量大于0.5，网格数量大于50万。
进一步地，所述步骤S1中，所述结构参数包括转炉炉体的结构参数、氧枪拉瓦尔喷头的结构参数和底部元件的结构参数；所述进行合理简化的方法为：
对于转炉炉体，简化后的考虑因素包括炉膛内部几何形状、熔池钢液区域及熔池钢液区域上部700-2000 mm范围处的炉渣和气相区域；
对于氧枪拉瓦尔喷头，简化后的考虑因素包括收缩段、喉口以及扩张段内部区域；
将底部元件简化成圆孔形状。
进一步地，所述步骤S3中，所述设定模型求解的基本假设包括：
（1）氧气为可压缩牛顿流体；氩气、钢液和炉渣为不可压缩牛顿流体；
（2）氧气满足粘性萨瑟兰定律，其它流体的物性参数保持恒定。
进一步地，所述步骤S3中，所述选择计算模型，具体为：
选择VOF多相流模型对气-渣-金界面进行求解，选择DPM离散相模型对底吹氩气泡进行求解，选择能量守恒方程对高温熔池内部传热行为进行求解，选择标准k-ε湍流模型对熔池流场进行求解。
进一步地，所述步骤S3中，所述设定物性参数包括：设定氧气和氩气的物性参数值为300 K温度条件下所对应的物性参数值，炉渣和钢液的物性参数值为1873 K温度条件下所对应的物性参数值；
所述设定边界条件时包括：顶吹入口采用Profile边界条件，氧气相体积分数为1，顶吹入口温度为300 K；顶面出口采用压力出口边界条件，顶面出口温度为1873 K；底吹入口采用壁面边界条件同时DPM选项处的边界采用escape类型，氩气相体积分数为1，氩气的注入温度为300 K；转炉炉体除了顶吹入口、顶面出口、底吹入口的其余部位采用壁面边界条件，并将所述其余部位设定为绝热壁面；
所述步骤S3中，所述设定求解算法包括：压力-速度耦合采用PISO算法，压力插值采用PRESTO!算法，自由界面插值采用Geo-Reconstruct算法，控制方程对流项的离散使用二阶迎风格式。其中，PRESTO!算法（PREssure STaggering Option，PRESTO!）为交错压力算法，用于对压力项进行离散，以使压力更好地同时满足动量方程和连续性方程；PRESTO!算法是fluent软件中一个常规选项，为常规技术。
进一步地，所述步骤S3中，所述进行初始化操作包括：采用混合初始化方式对计算域进行初始化操作，同时通过Patch功能将计算域中钢液区域、炉渣区域以及氧气区域所对应的相体积分数设置为1。
进一步地，所述步骤S4，具体包括：
步骤S41，获取转炉实际生产过程中的烟气信息数据，所述烟气信息数据包括不同冶炼时刻转炉烟气总流量、CO和CO2体积比例，采用公式（1）和（2）计算得到不同冶炼时刻熔池气相和液相部分碳-氧反应气泡总量，公式如下：
（1）
式（1）中，为脱碳速率，单位为：%·s-1；Qgas为转炉烟气总流量，单位为：Nm3·s-1；φ (CO+ CO2)为CO和CO2的体积比例，单位为：%；Wsteel为钢液质量，单位为：t；
（2）
式（2）中，GCO为熔池气相和液相部分碳-氧反应气泡总量，单位为：kg·s-1；
步骤S42，通过FactSage软件中Multiple Interconnected Equilibrium和Adiabatic Stoichiometric Reactor方法结合Macro programming功能计算得到不同冶炼时刻转炉熔池液相部分碳-氧反应比例λ；
步骤S43，采用公式（3）计算得到熔池液相部分碳-氧反应气泡量，公式如下：
QCO=λ·GCO（3）
式（3）中，QCO为熔池液相部分碳-氧反应气泡量，单位为：kg·s-1。
进一步地，所述步骤S5，具体包括：
步骤S51：通过ANSYS Fluent中的Bounded和Sample Points功能建立熔池内部碳-氧反应气泡生成点，保证熔池上部碳-氧反应气泡生成点数量高于熔池下部区域；
步骤S52：UDF程序中调用DEFINE_DPM_SCALAR_UPDATE宏，当熔池液相碳-氧反应气泡上升到渣-金界面相体积分数为0.5的区域时，气泡被自动去除。
一种碳-氧反应作用下炼钢转炉熔池流场仿真系统，包括：
转炉三维几何模型构建模块：获取实际生产中炼钢转炉的结构参数，进行合理简化，建立转炉三维几何模型；
结构化网格划分模块：将所述转炉三维几何模型进行结构化网格划分处理，得到转炉网格模型；
顶吹、底吹以及碳-氧反应作用下的转炉熔池模拟仿真模块：
用于将所述转炉网格模型导入仿真软件，设定模型求解的基本假设，选择计算模型，设定物性参数、设定边界条件以及设定求解算法，随后进行初始化操作，最后进行迭代计算，得到顶吹和底吹复合吹炼条件下的转炉熔池流场；获取转炉实际生产过程中的烟气信息数据，根据碳守恒计算得到不同冶炼时刻熔池气相和液相部分碳-氧反应气泡总量，计算得到转炉不同冶炼时刻熔池液相部分碳-氧反应比例，得到熔池液相部分碳-氧反应气泡量；建立熔池内部碳-氧反应气泡生成点，将所述熔池液相部分碳-氧反应气泡量设定为碳-氧反应的输入流量，对碳-氧反应气泡在渣-金界面消失行为进行编译，得到熔池内部碳-氧反应气泡行为；以及，将所述熔池内部碳-氧反应气泡行为与得到的所述顶吹和底吹复合吹炼条件下的转炉熔池流场进行耦合计算，求解得到顶吹、底吹以及碳-氧反应作用下的转炉熔池模拟结果；
后处理模块：对所述碳-氧反应作用下的转炉熔池模拟结果进行后处理操作，得到不同吹炼时刻下碳-氧反应作用下的炼钢转炉熔池流场。
与现有技术相比，本发明的有益效果如下：
本发明提供的一种碳-氧反应作用下炼钢转炉熔池流场仿真方法，通过考虑熔池液相部分碳-氧反应生成的CO气泡上浮行为，模拟了转炉冶炼过程中三种能量源：顶吹超音速氧气射流、底吹惰性气泡羽流以及熔池内部碳-氧反应生成的CO气泡群流作用下的熔池流场，所得模拟结果更加符合实际转炉熔池内部流动行为，可以获得碳-氧反应作用下熔池动力学特性的变化规律，上述结果对于实际转炉的设计优化、底吹工艺的优化、吹炼制度的优化以及转炉高效稳定的生产具有重要的指导意义。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图进行简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被当作是对本发明范围的限定。
图1为本发明实施例中碳-氧反应作用下炼钢转炉熔池流场仿真方法流程图；
图2为本发明实施例中转炉三维几何模型；
图3为本发明实施例中转炉网格模型；
图4为本发明实施例中顶底复合吹炼条件下转炉熔池流场；
图5为本发明实施例中转炉烟气流量、CO和CO2气体比例随时间变化图；
图6为本发明实施例中熔池液相和气相部分碳-氧反应气泡量随时间变化图；
图7为本发明实施例中熔池液相部分碳-氧反应比例随时间变化图；
图8为本发明实施例中熔池内部碳-氧反应气泡生成点图；
图9为本发明实施例中碳-氧反应作用下冶炼时刻为100 s、200 s、400 s、700 s时的炼钢转炉熔池流场图；
图10为本发明实施例中模型准确性验证图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
本实施例提供了一种碳-氧反应作用下炼钢转炉熔池流场仿真方法，如图1所示，该方法具体包括以下步骤：
步骤S1：根据某炼钢厂提供的100吨转炉炉体及其配套的4孔氧枪拉瓦尔喷头和底部元件设计图纸，使用SolidWorks软件构建了转炉炉膛内部包括钢液区域及其上部800 mm范围的区域，氧枪拉瓦尔喷头只考虑了收缩段、喉口以及扩张段内部区域，底部元件简化成了圆孔形状的转炉三维几何模型，如图2所示，保存为.STEP格式文件；
步骤S2：将格式为.STEP的转炉三维几何模型导入ICEM软件中，首先对每个面进行命名，然后进行Edge关联和O-Block划分处理，如图3所示，在氧枪拉瓦尔喷头区域、底部元件区域以及气-渣-金界面区域通过增加节点数量的方式进行了网格加密处理，整个模型的网格质量为0.55，网格数量为100万，保存为.msh格式文件;
步骤S3：将.msh格式的转炉网格模型文件导入ANSYS Fluent软件中，其中，模型求解的基本假设包括：
（1）氧气为可压缩牛顿流体；氩气、钢液和炉渣为不可压缩牛顿流体；
（2）氧气满足粘性萨瑟兰定律，其它流体物性参数保持恒定。
根据模型内部流场行为的具体情况，选择VOF多相流模型对气-渣-金界面进行求解，选择DPM离散相模型对底吹氩气泡进行求解，开启能量守恒方程对高温熔池内部传热行为进行求解，选择标准k-ε湍流模型对熔池流场进行求解。
设定物性参数，其中，氧气为主项，温度为300 K，密度为可压缩，比热容为919.31J·kg-1·K-1，导热系数为0.0246 W·m-1·K-1，粘度符合萨瑟兰定律；钢液为第二相，温度为1873 K，密度为7100 kg·m-3，比热容为670 J·kg-1·K-1，导热系数为40 W·m-1·K-1，粘度为0.0065 kg·m-1·s-1；炉渣为第三项，温度为1873 K，密度为3500 kg·m-3，比热容为1200J·kg-1·K-1，导热系数为1.7 W·m-1·K-1，粘度为0.1 kg·m-1·s-1；氩气为DPM离散相，温度为300 K，密度为1.6228 kg·m-3，比热容为520.64 J·kg-1·K-1；
设定边界条件，顶吹入口采用profile边界条件，氧气相体积分数为1，温度为300K；顶面出口采用压力出口边界条件，温度为1873 K；底吹入口采用壁面边界条件同时DPM采用escape类型，氩气相体积分数为1，温度为300 K；其余部位采用壁面边界条件，并设定为绝热壁面；
设定求解算法，压力-速度耦合采用PISO算法，压力插值采用PRESTO!算法，自由界面插值采用Geo-Reconstruct算法，控制方程对流项的离散使用二阶迎风格式；
进行初始化操作，采用混合初始化方式对计算域进行初始化操作，同时通过Patch功能将计算域中钢液区域、炉渣区域以及氧气区域所对应的相体积分数设置为1；进行迭代计算，时间步长采取自适应模型；每0.1 s时间间隔对计算结果进行保存，并导出.cas和.dat格式的文件，如图4所示，文件通过Fluent自带的后处理功能得到顶底复合吹炼条件下的转炉熔池流场；
步骤S4：针对100吨转炉生产过程，通过原位式激光分析仪获取转炉冶炼周期内烟气流量、CO和CO2气体的体积比例，如图5所示；根据公式（1）和（2）计算得到不同冶炼时刻熔池气相和液相部分碳-氧反应气泡总量，如图6所示；结合100吨转炉实际冶炼过程操作特征，通过FactSage软件中Multiple Interconnected Equilibrium和AdiabaticStoichiometric Reactor方法结合Macro programming功能计算得到不同冶炼时刻熔池液相部分碳-氧反应比例λ，如图7所示；最后通过公式（3）计算得到不同冶炼时刻熔池液相部分碳-氧反应气泡量；
（1）
式（1）中为脱碳速率，%·s-1；Qgas为转炉烟气总流量，Nm3·s-1；φ (CO)和φ
(CO2)为CO和CO2的气体体积比例，%；Wsteel为钢液质量，t；
（2）
式（2）中GCO为熔池气相和液相部分碳-氧反应气泡总量，kg·s-1。
QCO=λ·GCO（3）
式（3）中QCO为熔池液相部分碳-氧反应气泡量kg·s-1。
步骤S5：在步骤S3的结果模型中，通过Bounded和Sample Points功能建立了熔池内部碳-氧反应气泡生成点，如图8所示；而后选择特征冶炼时刻为100 s、200 s、400 s、700s时熔池液相部分碳-氧反应气泡流量通过DPM模型将其作为输入流量进行注入；开启UDF程序，调用DEFINE_DPM_SCALAR_UPDATE宏功能，使熔池液相碳-氧反应气泡上升到渣-金界面相体积分数为0.5区域时，气泡被自动去除；
步骤S6：将所述熔池内部碳-氧反应气泡行为与步骤S3得到的顶底复合吹炼条件下的转炉熔池流场进行耦合计算，计算时间步长采用自适应模型，每0.1 s时间间隔对模拟计算结果进行保存，并导出.cas和.dat格式的文件；
步骤S7：将步骤S6计算得到的.cas和.dat格式的文件导入Tecplot软件进行后处理操作，提取纵截面，得到冶炼时刻为100 s、200 s、400 s、700 s时的碳-氧反应作用下的转炉熔池流场，如图9所示，并将图片导出.GIF格式的文件。
本发明还提供一种碳-氧反应作用下炼钢转炉熔池流场仿真系统实施例，包括：
转炉三维几何模型构建模块：获取实际生产中炼钢转炉的结构参数，进行合理简化，建立转炉三维几何模型；
结构化网格划分模块：将所述转炉三维几何模型进行结构化网格划分处理，得到转炉网格模型；
顶吹、底吹以及碳-氧反应作用下的转炉熔池模拟仿真模块：用于将所述转炉网格模型导入仿真软件，设定模型求解的基本假设，选择计算模型，设定物性参数、设定边界条件以及设定求解算法，随后进行初始化操作，最后进行迭代计算，得到顶吹和底吹复合吹炼条件下的转炉熔池流场；获取转炉实际生产过程中的烟气信息数据，根据碳守恒计算得到不同冶炼时刻熔池气相和液相部分碳-氧反应气泡总量，计算得到转炉不同冶炼时刻熔池液相部分碳-氧反应比例，得到熔池液相部分碳-氧反应气泡量；建立熔池内部碳-氧反应气泡生成点，将所述熔池液相部分碳-氧反应气泡量设定为碳-氧反应的输入流量，对碳-氧反应气泡在渣-金界面消失行为进行编译，得到熔池内部碳-氧反应气泡行为；以及，将所述熔池内部碳-氧反应气泡行为与得到的所述顶吹和底吹复合吹炼条件下的转炉熔池流场进行耦合计算，求解得到顶吹、底吹以及碳-氧反应作用下的转炉熔池模拟结果；
后处理模块：对所述碳-氧反应作用下的转炉熔池模拟结果进行后处理操作，得到不同吹炼时刻下碳-氧反应作用下的炼钢转炉熔池流场。
为了确保本发明所建立模型的准确性，依据本实施例所得到的模拟结果首先与水力学实验所得到的混匀时间进行了比较，混匀时间是表征转炉熔池流动与混匀效果的特征参数，如图10所示，顶吹、底吹以及复吹工况条件下，模拟结果与实验结果的误差在8.05%-15.86%范围内；此外，纯顶吹工况下，熔池动能WTop为1552 J，在考虑冶炼中期400 s时熔池碳-氧反应所产生的CO气泡时，熔池动能WTop+CO为7351 J，WTop/WTop+CO=21.1%，根据陈家祥在《钢铁冶金学》一书中报道，针对100-200 t转炉而言，通过理论计算WTop/WTop+CO处于19.8~24.3%范围内，这与本实施例得到的模拟结果吻合良好。综上而言，本发明所建立的模型是准确且更加贴近于实际生产情况。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种碳-氧反应作用下炼钢转炉熔池流场仿真方法，其特征在于，包括如下步骤：
步骤S1：获取实际生产中炼钢转炉的结构参数，进行合理简化，建立转炉三维几何模型；
步骤S2：将所述转炉三维几何模型进行结构化网格划分处理，得到转炉网格模型；
步骤S3：将所述转炉网格模型导入仿真软件，设定模型求解的基本假设，选择计算模型，设定物性参数、设定边界条件以及设定求解算法，随后进行初始化操作，最后进行迭代计算，得到顶吹和底吹复合吹炼条件下的转炉熔池流场；
步骤S4：获取转炉实际生产过程中的烟气信息数据，根据碳守恒计算得到不同冶炼时刻熔池气相和液相部分碳-氧反应气泡总量，计算得到转炉不同冶炼时刻熔池液相部分碳-氧反应比例，得到熔池液相部分碳-氧反应气泡量；
步骤S5：建立熔池内部碳-氧反应气泡生成点，将所述熔池液相部分碳-氧反应气泡量设定为碳-氧反应的输入流量，对碳-氧反应气泡在渣-金界面消失行为进行编译，得到熔池内部碳-氧反应气泡行为；
步骤S6：将所述熔池内部碳-氧反应气泡行为与步骤S3得到的顶吹和底吹复合吹炼条件下的转炉熔池流场进行耦合计算，求解得到顶吹、底吹以及碳-氧反应作用下的转炉熔池模拟结果；
步骤S7：对所述碳-氧反应作用下的转炉熔池模拟结果进行后处理操作，得到不同吹炼时刻下碳-氧反应作用下的炼钢转炉熔池流场。
2.根据权利要求1所述一种碳-氧反应作用下炼钢转炉熔池流场仿真方法，其特征在于，所述步骤S2中，将所述转炉三维几何模型进行结构化网格划分处理的过程中，对氧枪拉瓦尔喷头区域、底部元件区域以及气-渣-金界面区域进行局部网格加密处理，所述加密处理具体为：计算域采用六面体结构化网格划分，整体网格质量大于0.5，网格数量大于50万。
3.根据权利要求1所述一种碳-氧反应作用下炼钢转炉熔池流场仿真方法，其特征在于，所述步骤S1中，所述结构参数包括转炉炉体的结构参数、氧枪拉瓦尔喷头的结构参数和底部元件的结构参数；所述进行合理简化的方法为：
对于转炉炉体，简化后的考虑因素包括炉膛内部几何形状、熔池钢液区域及熔池钢液区域上部700-2000 mm范围处的炉渣和气相区域；
对于氧枪拉瓦尔喷头，简化后的考虑因素包括收缩段、喉口以及扩张段内部区域；
将底部元件简化成圆孔形状。
4.根据权利要求1所述一种碳-氧反应作用下炼钢转炉熔池流场仿真方法，其特征在于，所述步骤S3中，所述设定模型求解的基本假设包括：
（1）氧气为可压缩牛顿流体；氩气、钢液和炉渣为不可压缩牛顿流体；
（2）氧气满足粘性萨瑟兰定律，其它流体的物性参数保持恒定。
5.根据权利要求1所述一种碳-氧反应作用下炼钢转炉熔池流场仿真方法，其特征在于，所述步骤S3中，所述选择计算模型，具体为：
选择VOF多相流模型对气-渣-金界面进行求解，选择DPM离散相模型对底吹氩气泡进行求解，选择能量守恒方程对高温熔池内部传热行为进行求解，选择标准k-ε湍流模型对熔池流场进行求解。
6.根据权利要求1所述一种碳-氧反应作用下炼钢转炉熔池流场仿真方法，其特征在于，所述步骤S3中，所述设定物性参数包括：设定氧气和氩气的物性参数值为300 K温度条件下所对应的物性参数值，炉渣和钢液的物性参数值为1873 K温度条件下所对应的物性参数值；
所述设定边界条件时包括：顶吹入口采用Profile边界条件，氧气相体积分数为1，顶吹入口温度为300 K；顶面出口采用压力出口边界条件，顶面出口温度为1873 K；底吹入口采用壁面边界条件同时DPM选项处的边界采用escape类型，氩气相体积分数为1，氩气的注入温度为300 K；转炉炉体除了顶吹入口、顶面出口、底吹入口的其余部位采用壁面边界条件，并将所述其余部位设定为绝热壁面；
所述步骤S3中，所述设定求解算法包括：压力-速度耦合采用PISO算法，压力插值采用PRESTO!算法，自由界面插值采用Geo-Reconstruct算法，控制方程对流项的离散使用二阶迎风格式；其中PRESTO!算法为交错压力算法，用于对压力项进行离散。
7.根据权利要求1所述一种碳-氧反应作用下炼钢转炉熔池流场仿真方法，其特征在于，所述步骤S3中，所述进行初始化操作包括：采用混合初始化方式对计算域进行初始化操作，同时通过Patch功能将计算域中钢液区域、炉渣区域以及氧气区域所对应的相体积分数设置为1。
8.根据权利要求1所述一种碳-氧反应作用下炼钢转炉熔池流场仿真方法，其特征在于，所述步骤S4，具体包括：
步骤S41，获取转炉实际生产过程中的烟气信息数据，所述烟气信息数据包括不同冶炼时刻转炉烟气总流量、CO和CO2体积比例，采用公式（1）和（2）计算得到不同冶炼时刻熔池气相和液相部分碳-氧反应气泡总量，公式如下：
（1）
式（1）中，为脱碳速率，单位为：%·s-1；Qgas为转炉烟气总流量，单位为：Nm3·s-1；φ(CO+ CO2)为CO和CO2的体积比例，单位为：%；Wsteel为钢液质量，单位为：t；
（2）
式（2）中，GCO为熔池气相和液相部分碳-氧反应气泡总量，单位为：kg·s-1；
步骤S42，通过FactSage软件中Multiple Interconnected Equilibrium和AdiabaticStoichiometric Reactor方法结合Macro programming功能计算得到不同冶炼时刻转炉熔池液相部分碳-氧反应比例λ；
步骤S43，采用公式（3）计算得到熔池液相部分碳-氧反应气泡量，公式如下：
QCO=λ·GCO（3）
式（3）中，QCO为熔池液相部分碳-氧反应气泡量，单位为：kg·s-1。
9.根据权利要求1所述一种碳-氧反应作用下炼钢转炉熔池流场仿真方法，其特征在于，所述步骤S5，具体包括：
步骤S51：通过ANSYS Fluent中的Bounded和Sample Points功能建立熔池内部碳-氧反应气泡生成点，保证熔池上部碳-氧反应气泡生成点数量高于熔池下部区域；
步骤S52：UDF程序中调用DEFINE_DPM_SCALAR_UPDATE宏，当熔池液相碳-氧反应气泡上升到渣-金界面相体积分数为0.5的区域时，气泡被自动去除。
10.一种碳-氧反应作用下炼钢转炉熔池流场仿真系统，其特征在于，包括：
转炉三维几何模型构建模块：获取实际生产中炼钢转炉的结构参数，进行合理简化，建立转炉三维几何模型；
结构化网格划分模块：将所述转炉三维几何模型进行结构化网格划分处理，得到转炉网格模型；
顶吹、底吹以及碳-氧反应作用下的转炉熔池模拟仿真模块：用于将所述转炉网格模型导入仿真软件，设定模型求解的基本假设，选择计算模型，设定物性参数、设定边界条件以及设定求解算法，随后进行初始化操作，最后进行迭代计算，得到顶吹和底吹复合吹炼条件下的转炉熔池流场；获取转炉实际生产过程中的烟气信息数据，根据碳守恒计算得到不同冶炼时刻熔池气相和液相部分碳-氧反应气泡总量，计算得到转炉不同冶炼时刻熔池液相部分碳-氧反应比例，得到熔池液相部分碳-氧反应气泡量；建立熔池内部碳-氧反应气泡生成点，将所述熔池液相部分碳-氧反应气泡量设定为碳-氧反应的输入流量，对碳-氧反应气泡在渣-金界面消失行为进行编译，得到熔池内部碳-氧反应气泡行为；以及，将所述熔池内部碳-氧反应气泡行为与得到的所述顶吹和底吹复合吹炼条件下的转炉熔池流场进行耦合计算，求解得到顶吹、底吹以及碳-氧反应作用下的转炉熔池模拟结果；
后处理模块：对所述碳-氧反应作用下的转炉熔池模拟结果进行后处理操作，得到不同吹炼时刻下碳-氧反应作用下的炼钢转炉熔池流场。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410139965.2A CN117688819B (zh) | 2024-02-01 | 2024-02-01 | 一种碳-氧反应作用下炼钢转炉熔池流场仿真方法及仿真系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410139965.2A CN117688819B (zh) | 2024-02-01 | 2024-02-01 | 一种碳-氧反应作用下炼钢转炉熔池流场仿真方法及仿真系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117688819A CN117688819A (zh) | 2024-03-12 |
| CN117688819B true CN117688819B (zh) | 2024-04-26 |
Family
ID=90135660
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410139965.2A Active CN117688819B (zh) | 2024-02-01 | 2024-02-01 | 一种碳-氧反应作用下炼钢转炉熔池流场仿真方法及仿真系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117688819B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118657039B (zh) * | 2024-08-16 | 2024-10-29 | 北京科技大学 | 一种连铸结晶器内颗粒运动及相互作用行为的仿真方法 |
| CN118682090B (zh) * | 2024-08-28 | 2024-11-05 | 北京科技大学 | 一种连铸机钢包至中间包钢水温度的稳定控制方法和系统 |
| CN118734659B (zh) * | 2024-09-02 | 2024-11-12 | 东北大学 | 一种电弧炉熔池温度的动态软测量方法及系统 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6284016B1 (en) * | 1997-03-21 | 2001-09-04 | Nippon Steel Corporation | Pressure converter steelmaking method |
| CN202519289U (zh) * | 2012-04-10 | 2012-11-07 | 中冶赛迪工程技术股份有限公司 | 脉动增强型金属熔池脉动气体搅拌装置 |
| CN114970391A (zh) * | 2022-05-25 | 2022-08-30 | 安徽工业大学 | 一种炼钢过程碳氧反应气泡动力学水模型实验方法及装置 |
| CN115261551A (zh) * | 2022-06-29 | 2022-11-01 | 北京科技大学 | 一种转炉底吹工艺优化方法及优化系统 |
| CN115270654A (zh) * | 2022-07-05 | 2022-11-01 | 北京科技大学 | 一种转炉炼钢出钢过程的数值模拟方法 |
| CN117010287A (zh) * | 2023-05-18 | 2023-11-07 | 中国恩菲工程技术有限公司 | 熔池熔炼流固耦合仿真计算方法及装置、存储介质、终端 |
| CN117171983A (zh) * | 2023-08-21 | 2023-12-05 | 北京科技大学 | 基于关键场耦合的富氧底吹熔炼炉建模仿真方法及装置 |
-
2024
- 2024-02-01 CN CN202410139965.2A patent/CN117688819B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6284016B1 (en) * | 1997-03-21 | 2001-09-04 | Nippon Steel Corporation | Pressure converter steelmaking method |
| CN202519289U (zh) * | 2012-04-10 | 2012-11-07 | 中冶赛迪工程技术股份有限公司 | 脉动增强型金属熔池脉动气体搅拌装置 |
| CN114970391A (zh) * | 2022-05-25 | 2022-08-30 | 安徽工业大学 | 一种炼钢过程碳氧反应气泡动力学水模型实验方法及装置 |
| CN115261551A (zh) * | 2022-06-29 | 2022-11-01 | 北京科技大学 | 一种转炉底吹工艺优化方法及优化系统 |
| CN115270654A (zh) * | 2022-07-05 | 2022-11-01 | 北京科技大学 | 一种转炉炼钢出钢过程的数值模拟方法 |
| CN117010287A (zh) * | 2023-05-18 | 2023-11-07 | 中国恩菲工程技术有限公司 | 熔池熔炼流固耦合仿真计算方法及装置、存储介质、终端 |
| CN117171983A (zh) * | 2023-08-21 | 2023-12-05 | 北京科技大学 | 基于关键场耦合的富氧底吹熔炼炉建模仿真方法及装置 |
Non-Patent Citations (1)
| Title |
|---|
| 基于"炉-机匹配"的炼钢-连铸过程多工序运行仿真;杨建平 等;《中南大学学报（自然科学版）》;20200331;第第51卷卷(第第3期期);第1-13页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117688819A (zh) | 2024-03-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN117688819B (zh) | 一种碳-氧反应作用下炼钢转炉熔池流场仿真方法及仿真系统 | |
| CN105930585A (zh) | 基于CFD的Shell气化炉流场和温度场的模拟方法 | |
| CN115270654B (zh) | 一种转炉炼钢出钢过程的数值模拟方法 | |
| Ding et al. | Energy efficiency assessment of hydrogen recirculation ejectors for proton exchange membrane fuel cell (PEMFC) system | |
| Hashimoto et al. | Transient model-based operation guidance on blast furnace | |
| Zhou et al. | Toward the bath flow interaction of a 250 t combined top and bottom blowing converter based on a mathematical modeling | |
| Samdani et al. | A process model for underground coal gasification–Part-III: Parametric studies and UCG process performance | |
| Lv et al. | Behavior of gas–slag–metal emulsion with nozzle‐twisted lance in converter steelmaking process | |
| CN114970383B (zh) | 一种co2和h2o气氛下渣/金熔滴脱碳行为的数值模拟方法 | |
| Huang et al. | Comprehensive numerical modeling of raceways in blast furnaces | |
| CN115612767B (zh) | 一种热交换区温度确定方法、系统、电子设备及存储介质 | |
| CN105045994B (zh) | 一种优化带式焙烧机热工制度的方法 | |
| CN106929615A (zh) | 一种高炉炉缸数值模拟与解析系统及其控制方法 | |
| CN117972999A (zh) | 一种基于Fluent-UDF的转炉熔池火点区传热行为模拟方法 | |
| Liu et al. | Numerical Study of the Kinetics of Dephosphorization with Single-Flow Post-combustion Oxygen Lance in Converter Steelmaking Process | |
| Dong et al. | Numerical study on the combustion phenomena in the gas region of a top-blown converter with single-flow post-combustion oxygen lance | |
| Puviani et al. | A novel Ansys CFX–RELAP5 coupling tool for the transient thermal-hydraulic analysis of liquid metal systems | |
| Jiang et al. | Numerical simulation of heat transfer behavior in hot spot zone of converter molten bath | |
| CN117171983A (zh) | 基于关键场耦合的富氧底吹熔炼炉建模仿真方法及装置 | |
| CN109753702A (zh) | 一种燃气轮机化学回热系统的仿真方法 | |
| CN116050174A (zh) | 一种烟气外循环下焦炉结焦过程的数值模拟方法 | |
| CN101576939A (zh) | 一种发动机系统模块化建模仿真方法 | |
| CN111599415B (zh) | 一种基于计算机仿真的高炉数字化系统实现方法 | |
| CN115048827A (zh) | 一种能够同时计算燃烧室污染物CO和NOx数值的方法 | |
| Hellberg et al. | A model of gas injection into a blast furnace tuyere |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |