CN112464543A - 一种计算vim冶炼过程中的夹杂物运动的方法 - Google Patents
一种计算vim冶炼过程中的夹杂物运动的方法 Download PDFInfo
- Publication number
- CN112464543A CN112464543A CN202110114283.2A CN202110114283A CN112464543A CN 112464543 A CN112464543 A CN 112464543A CN 202110114283 A CN202110114283 A CN 202110114283A CN 112464543 A CN112464543 A CN 112464543A
- Authority
- CN
- China
- Prior art keywords
- alloy
- calculation domain
- domain
- setting
- calculation
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 117
- 230000008569 process Effects 0.000 title claims abstract description 68
- 238000003723 Smelting Methods 0.000 title claims abstract description 59
- 230000033001 locomotion Effects 0.000 title claims abstract description 55
- 239000000956 alloy Substances 0.000 claims abstract description 153
- 239000002245 particle Substances 0.000 claims abstract description 87
- 230000006698 induction Effects 0.000 claims abstract description 41
- 239000000463 material Substances 0.000 claims abstract description 38
- 238000004088 simulation Methods 0.000 claims abstract description 34
- 230000000704 physical effect Effects 0.000 claims abstract description 9
- 238000004364 calculation method Methods 0.000 claims description 187
- 229910045601 alloy Inorganic materials 0.000 claims description 145
- 230000005484 gravity Effects 0.000 claims description 30
- 239000007788 liquid Substances 0.000 claims description 27
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 25
- 230000001681 protective effect Effects 0.000 claims description 24
- 230000008878 coupling Effects 0.000 claims description 18
- 238000010168 coupling process Methods 0.000 claims description 17
- 238000005859 coupling reaction Methods 0.000 claims description 17
- 239000012774 insulation material Substances 0.000 claims description 17
- 238000011160 research Methods 0.000 claims description 13
- 238000004458 analytical method Methods 0.000 claims description 12
- 238000009826 distribution Methods 0.000 claims description 11
- 239000012530 fluid Substances 0.000 claims description 11
- 239000011810 insulating material Substances 0.000 claims description 11
- 230000001133 acceleration Effects 0.000 claims description 9
- 230000001419 dependent effect Effects 0.000 claims description 6
- 239000000395 magnesium oxide Substances 0.000 claims description 6
- CPLXHLVBOLITMK-UHFFFAOYSA-N magnesium oxide Inorganic materials [Mg]=O CPLXHLVBOLITMK-UHFFFAOYSA-N 0.000 claims description 6
- AXZKOIWUVFPNLO-UHFFFAOYSA-N magnesium;oxygen(2-) Chemical compound [O-2].[Mg+2] AXZKOIWUVFPNLO-UHFFFAOYSA-N 0.000 claims description 6
- 230000035699 permeability Effects 0.000 claims description 6
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 claims description 5
- YKTSYUJCYHOUJP-UHFFFAOYSA-N [O--].[Al+3].[Al+3].[O-][Si]([O-])([O-])[O-] Chemical compound [O--].[Al+3].[Al+3].[O-][Si]([O-])([O-])[O-] YKTSYUJCYHOUJP-UHFFFAOYSA-N 0.000 claims description 5
- PNEYBMLMFCGWSK-UHFFFAOYSA-N aluminium oxide Inorganic materials [O-2].[O-2].[O-2].[Al+3].[Al+3] PNEYBMLMFCGWSK-UHFFFAOYSA-N 0.000 claims description 5
- 229910052802 copper Inorganic materials 0.000 claims description 5
- 239000010949 copper Substances 0.000 claims description 5
- 238000009413 insulation Methods 0.000 claims description 5
- 230000001052 transient effect Effects 0.000 claims description 5
- 239000010425 asbestos Substances 0.000 claims description 4
- 239000011449 brick Substances 0.000 claims description 4
- ODINCKMPIJJUCX-UHFFFAOYSA-N calcium oxide Inorganic materials [Ca]=O ODINCKMPIJJUCX-UHFFFAOYSA-N 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 4
- 230000000694 effects Effects 0.000 claims description 4
- 230000005284 excitation Effects 0.000 claims description 4
- 239000002657 fibrous material Substances 0.000 claims description 4
- 230000002706 hydrostatic effect Effects 0.000 claims description 4
- 230000000737 periodic effect Effects 0.000 claims description 4
- 229910052895 riebeckite Inorganic materials 0.000 claims description 4
- BRPQOXSCLDDYGP-UHFFFAOYSA-N calcium oxide Chemical compound [O-2].[Ca+2] BRPQOXSCLDDYGP-UHFFFAOYSA-N 0.000 claims description 3
- 239000000292 calcium oxide Substances 0.000 claims description 3
- 230000008014 freezing Effects 0.000 claims description 3
- 238000007710 freezing Methods 0.000 claims description 3
- TWNQGVIAIRXVLR-UHFFFAOYSA-N oxo(oxoalumanyloxy)alumane Chemical compound O=[Al]O[Al]=O TWNQGVIAIRXVLR-UHFFFAOYSA-N 0.000 claims description 3
- 238000004321 preservation Methods 0.000 claims 1
- 238000000746 purification Methods 0.000 abstract description 5
- 230000000007 visual effect Effects 0.000 abstract description 5
- 238000013461 design Methods 0.000 abstract description 3
- 239000007789 gas Substances 0.000 description 24
- 230000021715 photosynthesis, light harvesting Effects 0.000 description 7
- 238000002474 experimental method Methods 0.000 description 4
- 230000008018 melting Effects 0.000 description 4
- 238000002844 melting Methods 0.000 description 4
- 239000000523 sample Substances 0.000 description 3
- 238000003756 stirring Methods 0.000 description 3
- PXHVJJICTQNCMI-UHFFFAOYSA-N Nickel Chemical compound [Ni] PXHVJJICTQNCMI-UHFFFAOYSA-N 0.000 description 2
- 239000002826 coolant Substances 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000012535 impurity Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000007769 metal material Substances 0.000 description 2
- 229910001069 Ti alloy Inorganic materials 0.000 description 1
- 229910000797 Ultra-high-strength steel Inorganic materials 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000010261 cell growth Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 229910052593 corundum Inorganic materials 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000002500 effect on skin Effects 0.000 description 1
- 230000005674 electromagnetic induction Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 230000012010 growth Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052759 nickel Inorganic materials 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000005295 random walk Methods 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000010935 stainless steel Substances 0.000 description 1
- 229910001220 stainless steel Inorganic materials 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
- 229910001845 yogo sapphire Inorganic materials 0.000 description 1
Images
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/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)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Continuous Casting (AREA)
Abstract
本发明涉及一种计算VIM冶炼过程中的夹杂物运动的方法。所述方法:建立真空感应炉几何模型;设置真空感应炉的材料物性参数;在模型中添加并设置磁场;在模型中添加并设置湍流k‑ ε流场;在模型中添加并设置粒子跟踪物理场;对真空感应炉模型进行网格划分;配置求解器并求解;结果分析。本发明首次提供了一种基于COMSOL Multiphysics模拟计算VIM冶炼过程中的夹杂物运动轨迹的方法,可对密闭复杂的VIM冶炼过程中的夹杂物的运动轨迹和去除情况进行可视化仿真模拟,获得实验难以观测的结果,为VIM工艺方案优化设计提供理论指导和技术支持,对实现合金材料的高纯净化目标具有重要意义。
Description
技术领域
本发明属于金属的生产与精炼技术领域,尤其涉及一种计算VIM冶炼过程中的夹杂物运动的方法。
背景技术
随着现代工业技术的迅速发展,对金属材料的性能要求越来越高,纯净化已成为金属材料开发与应用的一个重要方向。夹杂物严重影响合金材料的力学性能,尤其在冶炼高品质合金材料时更需对其进行严格控制。真空感应熔炼(VIM)在真空条件下进行,具有强烈的搅拌功能,能够有效去除合金中的杂质和气体,是生产镍基高温合金、钛合金、不锈钢、超高强度钢等特种合金材料的重要冶炼方法。对VIM冶炼过程中夹杂物的运动进行研究,从而促进夹杂物的去除,对实现合金材料的高纯净化目标具有重要意义。
真空感应熔炼具有强烈的搅拌功能,夹杂物在熔池搅拌作用下主要受质量力、曳力和Saffman升力的影响。因此若想探究VIM冶炼过程中夹杂物的运动情况,除考虑夹杂物和合金的材料属性外,还需对熔池的流场分布进行求解,最后在流场分布结果上对夹杂物进行受力分析和运动行为分析。
在VIM冶炼过程中,真空感应炉处于真空密闭条件下,且熔池温度很高,很难通过观察或实验的方法直接获得炉内的冶炼情况,因此无法准确对冶炼过程进行合理有效的控制。而数值模拟方法能够对复杂问题进行精确定量的描述,通过建立合理的数学模型,可对密闭复杂的VIM冶炼过程进行可视化仿真模拟。
为了实现VIM冶炼过程中的夹杂物运动的模拟仿真,需要综合考虑到真空感应炉中的磁场分布、流场分布和夹杂物受力情况,并对上述物理场进行耦合和求解。目前常见的仿真模拟软件很难针对磁场进行求解,且多物理场的耦合也很难实现,一般通常通过开发两个或多个仿真软件的耦合程序的方法实现磁场和流场的耦合。但耦合程序过于繁琐,灵活性较差,且模拟仿真结果往往并不理想。COMSOL Multiphysics(简记为COMSOL)是以有限元法为基础,通过求解偏微分方程或偏微分方程组来实现任意多物理场的直接耦合,并以高效的计算性能和杰出的多场双向直接耦合分析能力实现高度精确的数值仿真,广泛应用于声学、化学反应、电磁学、流体动力学、热传导、结构力学等领域。中国专利申请CN201811525642.8公开了一种基于COMSOL软件的面内型热电器件结构设计方法,该专利申请通过COMSOL对热电器件的磁场和温度场进行了耦合,对不同结构热电器件工作时的温度场分布进行了数值模拟,为热电器件的结构优化提供了理论指导和技术支持。中国专利申请CN202010038956.6公开了一种反应堆-回路冷却剂流场、温度场和应力场耦合计算方法,实现了反应堆-回路冷却剂流场、温度场和应力场的耦合数值模拟计算。
COMSOL软件拥有简洁友好的操作界面,可对磁场和流场进行耦合求解,并可添加粒子跟踪模型对夹杂物的运动轨迹进行计算分析。然而,基于COMSOL适用于模拟计算VIM冶炼过程中的夹杂物运动轨迹的方法却未见报道。
发明内容
为了解决现有技术存在的问题,本发明提供了一种计算VIM冶炼过程中的夹杂物运动的方法。本发明方法基于COMSOL模拟计算VIM冶炼过程中的夹杂物运动,从而促进夹杂物的去除,从而实现合金材料的高纯净化目标。
为了实现上述目的,本发明提供了一种计算VIM冶炼过程中的夹杂物运动的方法,所述方法包括如下步骤:
(1)利用COMSOL建立真空感应炉的几何模型;所述几何模型的计算域包括合金计算域、坩埚计算域、保温材料计算域、线圈计算域、冷凝水计算域和保护气体计算域,并将所述几何模型的最外层定义为无限元域;
(2)设置真空感应炉的材料物性参数;
(3)在所述几何模型中添加并设置磁场;所述磁场设置方式为线圈电流激励,所述线圈为单导线线圈组;设定所述无限元域的外边界为磁绝缘;对所述几何模型的磁场初始值进行设置,所述几何模型内的所有计算域的磁势均设置为零;在所述合金计算域内添加力计算节点,设置为在频域同时计算瞬时变量和周期平均变量,并对VIM冶炼过程中的合金计算域的计算单元上的洛伦兹力进行计算;
(4)在所述几何模型中添加并设置湍流流场;所述湍流流场选择湍流k-ε模型对所述合金计算域的流场进行模拟分析,湍流类型采用平均纳维-斯托克斯方程描述,采用标准壁函数处理近壁面,利用动网格的方法对所述合金计算域的液面波动现象进行模拟;将所述坩埚计算域的内壁面设置为无滑移壁面;将所述合金计算域内的合金熔液设置为包含重力的不可压缩流体,设置所述合金熔液的初始速度和初始压力为零,并启用静水压力补偿;在所述湍流流场内添加重力节点,对所述合金计算域所受的重力进行计算;并利用洛伦兹力多物理场对所述合金计算域所受的时均电磁力进行计算,且将时均电磁力作为体积力输入到湍流k-ε模型中;
(5)在所述几何模型中添加并设置粒子跟踪物理场;所述粒子为随机生成的夹杂物粒子,所述夹杂物粒子的初始速度为所在位置的合金熔液的流动速度;对所述夹杂物粒子的密度、直径、类型及电荷数进行设置,设置所述坩埚计算域的内壁面和所述合金计算域的液面的壁条件为冻结或反弹,以实现夹杂物粒子去除效果的仿真模拟;然后在所述合金计算域中添加重力节点、曳力节点和萨夫曼升力节点,对所述合金计算域中的夹杂物粒子在运动过程中的受力情况进行仿真模拟;
(6)对所述几何模型进行网格划分;
(7)配置求解器并求解;求解的过程依次包括求解模型流场和求解夹杂物粒子的运动轨迹及去除情况;
(8)结果分析;所述结果分析包括绘制VIM冶炼过程中的合金计算域流场分布云图和速度随时间的变化曲线、绘制VIM冶炼过程中的夹杂物粒子的运动轨迹图以及对合金计算域的液面和坩埚计算域的内壁面的夹杂物粒子去除率进行统计。
优选地,在步骤(1)中,建立所述几何模型所需的参数包括:炉料高度;坩埚的高度、宽度和厚度;保温材料的高度、宽度和厚度;线圈的内径和外径,线圈间的间隙及匝数。
优选地,在步骤(2)中,所述真空感应炉的材料包括合金、坩埚、保温材料、线圈、冷凝水和保护气体;所述材料物性参数包括所述合金、所述坩埚、所述保温材料、所述线圈、所述冷凝水和所述保护气体的相对磁导率、电导率、相对介电常数以及所述合金的密度和动力粘度。
优选地,所述坩埚采用氧化镁、氧化钙和氧化铝中的一种或多种材料制成;所述保温材料为硅酸铝纤维材料、石棉材料、高铝砖材料中的一种或多种;所述线圈采用紫铜材料制成;所述冷凝水为水;和/或所述保护气体为Ar气。
优选地,在步骤(3)中,所述几何模型内的所有计算域符合安培定律。
优选地,在步骤(4)中,将所述重力节点的加速度设置为当地重力加速度;在步骤(4)中,在利用洛伦兹力多物理场对所述合金计算域所受的时均电磁力进行计算时,所述洛伦兹力多物理场的接口为磁场和湍流流场,在所述湍流流场中添加体积力节点,体积力大小为时均电磁力,体积力的轴向力为z方向时均电磁力,体积力的径向力为r方向时均电磁力。
优选地,在步骤(4)中,为模拟所述合金计算域的液面波动现象,对所述合金计算域的液面进行以下设置:
(a)在湍流流场中添加自由表面节点,区域设置为合金计算域液面;
(b)将所述合金计算域设置为变形域动网格,初始变形量设置为零,网格滑移类型设置为超弹性滑移。
优选地,在步骤(6)中,对所述几何模型进行自由三角形网格划分;所述合金计算域的网格相比所述坩埚计算域、所述保温材料计算域、所述线圈计算域、所述冷凝水计算域和所述保护气体计算域的网格较细,并将所述坩埚计算域的内壁面和所述线圈计算域的内壁面设置为边界层网格。
优选地,所述合金计算域的网格为流体动力学极细化网格,所述坩埚计算域、所述保温材料计算域、所述线圈计算域、所述冷凝水计算域和所述保护气体计算域的网格为普通物理超细化网格。
优选地,在步骤(7)中:求解模型流场的过程为频域-瞬态研究,物理场的接口为磁场、湍流流场和动网格,利用洛伦兹力多物理场进行耦合,并且在求解过程中自动对所述几何模型重新进行网格划分;求解夹杂物粒子的运动轨迹及去除情况的过程为瞬态研究,物理场的接口为粒子跟踪物理场,将因变量值中的不求解变量值设置为求解模型流场得出的解。
本发明与现有技术相比至少具有如下的有益效果:
(1)目前关于VIM冶炼过程中的夹杂物运动的仿真模拟研究很少,并且现有模拟仿真结果也并不理想;本发明首次提供了一种基于COMSOL Multiphysics模拟计算VIM冶炼过程中的夹杂物运动轨迹的方法,弥补了这方面的不足。本发明方法利用洛伦兹力多物理场耦合磁场和湍流流场,对VIM合金计算域的流场分布进行求解,最后在流场求解结果的基础上通过粒子跟踪算法仿真模拟夹杂物的运动轨迹。
(2)本发明方法利用COMSOL Multiphysics软件对磁场、湍流流场和粒子跟踪物理场进行耦合,在保证其耦合的正确性下,对VIM冶炼过程中的夹杂物的运动轨迹进行仿真计算,可对密闭复杂的VIM冶炼过程进行可视化仿真模拟,获得实验难以观测的结果,从而对VIM冶炼过程中夹杂物的去除情况进行预测,本发明方法能够有效保证数值模拟结果的准确性,为VIM工艺方案优化设计提供理论指导和技术支持,对实现合金材料的高纯净化目标具有重要意义。
附图说明
图1是为本发明方法的仿真计算流程图。
图2是本发明实施例1的真空感应炉几何模型图。
图3是本发明实施例1的真空感应炉几何模型的网格划分结果图。
图4是本发明实施例1在不同熔炼时间下的VIM合金计算域流场分布云图。图中(a)的熔炼时间为100s;(b)的熔炼时间为200s;(c)的熔炼时间为300s;(d)的熔炼时间为400s。
图5是本发明实施例1的合金计算域中三个探测点的速度变化曲线图。
图6是本发明实施例1在不同运动时间下的VIM合金计算域中的夹杂物粒子位置图。图中(a)的运动时间为25s;(b)的运动时间为50s;(c)的运动时间为75s,(d)的运动时间为100s。
图7是本发明实施例1在VIM冶炼过程中的坩埚内壁面和合金计算域液面对于夹杂物粒子的去除率曲线图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合本发明中的实施例,对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明提供了一种计算VIM冶炼过程中的夹杂物运动的方法,所述方法包括如下步骤:
(1)利用COMSOL建立真空感应炉的几何模型(也记作真空感应炉几何模型或真空感应炉模型);所述几何模型的计算域包括合金计算域、坩埚计算域、保温材料计算域、线圈计算域、冷凝水计算域和保护气体计算域,并将所述几何模型的最外层定义为无限元域;在本发明中,合金计算域也可用熔池表示,熔池即合金熔化后熔液所在区域;具体地,真空感应炉模型的计算域包括:合金、坩埚、保温材料、线圈、冷凝水和保护气体;建立真空感应炉模型所需参数包括:炉料高度(即合金在坩埚中的高度),坩埚的高度、宽度和厚度,保温材料的高度、宽度和厚度,线圈的内径和外径,线圈间的间隙及匝数;在本发明中,根据上述参数,利用COMSOL建立真空感应炉模型,并将所述真空感应炉模型的最外层定义为无限元域;特别说明的是,无限元是一种几何上可以趋于无限远处的单元,通常应用于无界域问题的研究;无限元可将合理的坐标缩放并应用于要研究的物理区域外部的一个虚拟域上,从而实现有限元的延伸,反映近场的边界条件和特征,加快有限元的计算和提高求解收敛性。
(2)设置所述真空感应炉的材料物性参数;在本发明中,具体地,所述真空感应炉的材料包括合金、坩埚、保温材料、线圈、冷凝水和保护气体;所述材料物性参数包括所述合金、所述坩埚、所述保温材料、所述线圈、所述冷凝水和所述保护气体的相对磁导率、电导率、相对介电常数以及所述合金的密度和动力粘度;在本发明中,所述坩埚可选择MgO、CaO或Al2O3等材料,保温材料可选择硅酸铝纤维、石棉或高铝砖等材料,线圈可选择紫铜材料,冷凝水选择水,保护气体选择Ar气,本发明根据材料属性对上述材料的相对磁导率、电导率、相对介电常数进行设置;根据熔炼合金型号,对其相对磁导率、电导率、相对介电常数、密度和动力粘度进行设置。
(3)在所述几何模型中添加并设置磁场;所述磁场设置方式为线圈电流激励,所述线圈为单导线线圈组;设定所述无限元域的外边界为磁绝缘;对所述几何模型的磁场初始值进行设置,所述几何模型内的所有计算域的磁势均设置为零;在所述合金计算域内添加力计算节点,设置为在频域同时计算瞬时变量和周期平均变量,并对VIM冶炼过程中的合金计算域的计算单元上的洛伦兹力进行计算;特别说明的是,仿真模拟基于的理论为有限元计算理论,计算过程中会将需要分析的模型进行离散化,即将连续的弹性体划分成有限个单元组成的离散体,并把作用力按等效原则移置到各节点上,其中每个单元可称为计算单元;在本发明的步骤(3)中,根据实际情况设置电流大小,其余参数皆与材料属性保持一致;假设无限元域外边界为磁绝缘,即无穷远边界处磁力线平行,磁势为零;在步骤(3)中,所述几何模型内的全部计算域符合安培定律,并对合金计算域的计算单元上的洛伦兹力进行计算;
其中,安培定律为:
洛伦兹力计算公式如下:
(4)在所述几何模型中添加并设置湍流流场(即湍流k-ε流场(spf));所述湍流流场选择湍流k-ε模型对所述合金计算域的流场进行模拟分析,湍流类型采用平均纳维-斯托克斯方程描述,采用标准壁函数处理近壁面,利用动网格(变形域动网格)的方法对所述合金计算域的液面波动现象进行模拟;将所述坩埚计算域的内壁面设置为无滑移壁面;将所述合金计算域内的合金熔液设置为包含重力的不可压缩流体,设置所述合金熔液的初始速度和初始压力为零也即指的是设置所述合金计算域的流场的初始速度为零,并启用静水压力补偿,设置湍流动能为spf.kinit m2/s2,湍流动能耗散率(简记为湍流耗散率)为spf.epinitm2/s2;在所述湍流流场内添加重力节点,对所述合金计算域所受的重力进行计算;并利用洛伦兹力多物理场对所述合金计算域所受的时均电磁力进行计算,且将时均电磁力作为体积力(即洛伦兹力)输入到湍流k-ε模型中;特别说明的是,所述近壁面即所述合金计算域与所述坩埚计算域的内壁面的交界处;在本发明中,也将合金计算域的流场记作熔池流场;在本发明的步骤(4)中,设置所述合金计算域流场初始速度为零,合金计算域所受力为自身重力和洛伦兹力;在计算湍流流场时,需要在湍流流场中加入洛伦兹力;为计算合金计算域的计算单元上的洛伦兹力(即体积力),利用洛伦兹力多物理场对合金计算域所受的时均电磁力进行计算,并将时均电磁力作为体积力输入到湍流k-ε模型中,即输入到合金计算域的计算单元上;在本发明中,步骤(4)中的洛伦兹力设置的目的在于设置一个接收项,从而能够使步骤(3)中计算得出(输出)的洛伦兹力输入到步骤(4)中的湍流模型中。
在本发明中,选择湍流k-ε(spf)模型对VIM合金计算域流场进行模拟分析,湍流类型采用平均纳维-斯托克斯方程(RANS)描述;涉及到的公式如下所示:
动量守恒方程:
其中,ρ--合金熔液密度,kg/m3;
连续性方程:
湍流动能k方程:
其中,k--湍流动能,m2/s2;
T--数学运算符号,表示矩阵的转置;
湍流动能耗散率ε方程:
在本发明中,利用变形域动网格的方法对合金计算域液面波动现象进行模拟,涉及到的公式如下所示:
(5)在所述几何模型中添加并设置粒子跟踪物理场;所述粒子为随机生成的夹杂物粒子,所述夹杂物粒子的初始速度为所在位置的合金熔液的流动速度;对所述夹杂物粒子的密度、直径、类型及电荷数进行设置,设置所述坩埚计算域的内壁面和所述合金计算域的液面的壁条件为冻结或反弹,以实现夹杂物粒子去除效果的仿真模拟;然后在所述合金计算域中添加重力节点、曳力节点和萨夫曼升力(Saffman升力)节点,对所述合金计算域中的夹杂物粒子在运动过程中的受力情况进行仿真模拟;特别说明的是,在所述合金计算域中添加重力节点,COMSOL软件会自动将重力修改为质量力,即实际在所述合金计算域中添加的是质量力;所述质量力指的是夹杂物粒子所受重力与浮力之差;
在本发明中,将夹杂物粒子视为宏观粒子;夹杂物粒子所受质量力由夹杂物和合金熔液的密度差所引起,计算公式如下:
其中,F g --质量力,N;
m p --夹杂物粒子质量,kg;
g--重力加速度,m/s2;
ρ p --夹杂物粒子密度,kg/m3;
ρ--合金熔液密度,kg/m3。
夹杂物粒子在运动过程中受合金熔液曳力影响,根据斯托克斯定律对合金熔液曳力进行计算,公式如下:
其中,F D --曳力,N;
d p --夹杂物粒子直径,m;
v--夹子物粒子运动速度,m/s;
m p --夹杂物粒子质量,kg;
ρ p --夹杂物粒子密度,kg/m3;
μ--合金熔液粘度,Pa·s;
在本发明中,除质量力和曳力外,夹杂物粒子在运动过程中还受Saffman升力影响,其计算公式如下:
其中F s --Saffman升力,N;
d p --夹杂物粒子直径,m;
μ--合金熔液粘度,Pa·s;
v--夹子物粒子运动速度,m/s;
ρ--合金熔液密度,kg/m3;
特别说明的是,在公式(1)至(16)中,相同的字母代表的物理量一致。
(6)对所述几何模型进行网格划分;在本发明中,具体地,对所述几何模型进行自由三角形网格划分,其中合金计算域的网格较细,其余计算域网格较粗,并将坩埚内壁面和线圈内壁面设置为边界层网格;在本发明中,合金计算域的网格较细,可以保证数值模拟结果的准确性;其余计算域网格较粗,可以缩短运算时间;电磁感应加热过程中,在坩埚内壁面和线圈内壁面存在较为明显的趋肤效应,且该效应会显著影响模型计算结果,因此将坩埚内壁面和线圈内壁面设置为边界层网格,且优选为边界层网格大小要小于电流的趋肤深度,可以有效保证数值模拟结果的准确性。
(7)配置求解器并求解;求解的过程依次包括求解模型流场(记作研究1)和求解夹杂物粒子的运动轨迹及去除情况(记作研究2);在本发明中,具体地,对高温合金VIM冶炼过程中夹杂物运动轨迹求解模型进行分步计算,其中研究1求解模型流场,研究2求解夹杂物粒子的运动轨迹及去除情况;研究1为频域-瞬态研究,根据VIM冶炼的实际参数设置时间步及频率,并利用洛伦兹力多物理场耦合磁场和湍流流场;研究2为瞬态研究,根据VIM冶炼的实际情况设置时间步,因变量初始值由研究1计算得出(即将因变量值中的不求解变量值设置为研究1计算得出的解);特别说明的是,在本发明中,求解模型流场即为求解所述几何模型流场;特别地,在步骤(7)配置求解器并求解之后,步骤(8)结果分析之前,还包括判断是否收敛的步骤,即判断是否能够成功计算出所述几何模型的解;在求解几何模型时,网格会影响所述几何模型是否能顺利计算出几何模型的结果;若网格质量较差或与所述几何模型不匹配,则所述几何模型不能进行顺利运算,需要针对所述几何模型对网格重新进行划分。
(8)结果分析;所述结果分析包括绘制VIM冶炼过程中的合金计算域流场分布云图和速度随时间的变化曲线、绘制VIM冶炼过程中的夹杂物粒子的运动轨迹图以及对合金计算域的液面和坩埚计算域的内壁面的夹杂物粒子去除率进行统计;在本发明中,夹杂物粒子即VIM冶炼过程中的夹杂物。
本发明首次提供了一种基于COMSOL Multiphysics模拟计算VIM冶炼过程中的夹杂物(夹杂物粒子)运动轨迹的方法;本发明方法利用COMSOL Multiphysics软件对磁场、流场和粒子跟踪物理场进行耦合,在保证其耦合的正确性下,对VIM冶炼过程中的夹杂物的运动轨迹进行仿真计算,可对密闭复杂的VIM冶炼过程进行可视化仿真模拟,获得实验难以观测的结果,从而对VIM冶炼过程中夹杂物的去除情况进行预测,本发明方法能够有效保证数值模拟结果的准确性,为VIM工艺方案优化设计提供理论指导和技术支持,对实现合金材料的高纯净化目标具有重要意义。
根据一些优选的实施方式,在步骤(1)中,建立所述几何模型所需的参数包括:炉料高度;坩埚的高度、宽度和厚度;保温材料的高度、宽度和厚度;线圈的内径和外径,线圈间的间隙及匝数。
根据一些优选的实施方式,在步骤(2)中,所述真空感应炉的材料包括合金、坩埚、保温材料、线圈、冷凝水和保护气体;所述材料物性参数包括所述合金、所述坩埚、所述保温材料、所述线圈、所述冷凝水和所述保护气体的相对磁导率、电导率、相对介电常数以及所述合金的密度和动力粘度。
根据一些优选的实施方式,所述坩埚采用氧化镁(MgO)、氧化钙(CaO)和氧化铝(Al2O3)中的一种或多种材料制成;所述保温材料为硅酸铝纤维材料、石棉材料、高铝砖材料中的一种或多种;所述线圈采用紫铜材料制成;所述冷凝水为水;和/或所述保护气体为Ar气。
根据一些优选的实施方式,在步骤(3)中,所述几何模型内的所有计算域符合安培定律。
根据一些优选的实施方式,在步骤(4)中,将所述重力节点的加速度设置为当地重力加速度;在步骤(4)中,在利用洛伦兹力多物理场对所述合金计算域所受的时均电磁力进行计算时,所述洛伦兹力多物理场的接口为磁场和湍流流场,在所述湍流流场中添加体积力节点,体积力大小为时均电磁力,体积力的轴向力为z方向时均电磁力,体积力的径向力为r方向时均电磁力。在本发明中,具体地,例如在COMSOL软件中,将所述重力节点的加速度设置为-g_const m/s2,将轴向力设置为mf.FLtzavZ,径向力设置为mf.FLtzavR。
根据一些优选的实施方式,在步骤(4)中,为模拟所述合金计算域的液面波动现象,对所述合金计算域的液面进行以下设置:(a)在湍流流场中添加自由表面节点,区域设置为合金计算域液面;(b)将所述合金计算域设置为变形域动网格,初始变形量设置为零,网格滑移类型设置为超弹性滑移。
根据一些优选的实施方式,在步骤(6)中,对所述几何模型进行自由三角形网格划分;所述合金计算域的网格相比所述坩埚计算域、所述保温材料计算域、所述线圈计算域、所述冷凝水计算域和所述保护气体计算域的网格较细,并将所述坩埚计算域的内壁面和所述线圈计算域的内壁面设置为边界层网格。
根据一些优选的实施方式,所述合金计算域的网格为流体动力学极细化网格,所述坩埚计算域、所述保温材料计算域、所述线圈计算域、所述冷凝水计算域和所述保护气体计算域的网格为普通物理超细化网格;特别说明的是,COMSOL软件可对几何模型自动进行网格划分,其中可将网格校准为普通物理网格、流体动力学网格、等离子体网格和半导体网格,并预定义了9种网格大小,分别为极粗化网格、超粗化网格、较粗化网格、粗化网格、常规网格、细化网格、较细化网格、超细化网格、极细化网格,网格大小逐渐减小;对上述两个参数进行设置后,即可对目标区域进行网格划分;优选的是,所述合金计算域的网格为流体动力学极细化网格,自动划分后,其最大单元网格大小为0.00469m,最小单元网格大小为1.4×10-5m,最大单元增长率为1.05,曲率因子为0.2,狭窄区域分辨率为1;优选的是,所述坩埚计算域、所述保温材料计算域、所述线圈计算域、所述冷凝水计算域和所述保护气体计算域的网格为普通物理超细化网格,自动划分后,其最大单元网格大小为0.028m,最小单元网格大小为1.05×10-4m,最大单元增长率为1.2,曲率因子为0.25,狭窄区域分辨率为1。
根据一些优选的实施方式,在步骤(7)中:求解模型流场的过程为频域-瞬态研究,物理场的接口为磁场、湍流流场和动网格,利用洛伦兹力多物理场进行耦合,并且在求解过程中自动对所述几何模型重新进行网格划分;求解夹杂物粒子的运动轨迹及去除情况的过程为瞬态研究,物理场的接口为粒子跟踪物理场,将因变量值中的不求解变量值设置为求解模型流场得出的解。
下文将通过举例的方式对本发明进行进一步的说明,目的是有助于本领域技术人员对本发明的构思、技术方案有更完整、准确和深入的理解,并有助于其实施,但是本发明的保护范围不限于这些实施例。
实施例1
如图1所示,本实施例提供了一种基于COMSOL模拟计算VIM冶炼过程中的夹杂物运动的方法,所述方法包括以下步骤:
(1)建立真空感应炉几何模型。
如图2所示,建立真空感应炉几何模型,模型几何尺寸如下:炉料高度为195mm;坩埚尺寸为φ340mm×H280mm,壁厚20mm;保温材料尺寸φ400mm×H310mm,壁厚30mm;线圈内径为20mm,外径为28mm,线圈间间隙为8mm,线圈匝数为8匝。根据上述参数,利用COMSOL绘制真空感应炉几何模型,并将所述真空感应炉几何模型的最外层定义为无限元域,无限元域的坐标系选定为球面坐标系。
(2)设置真空感应炉的材料物性参数。
真空感应炉涉及到的材料包括:合金、坩埚、保温材料、线圈、冷凝水和保护气体。坩埚选择MgO材料,保温材料选择硅酸铝纤维材料,线圈选择紫铜材料,冷凝水选择水,保护气体选择Ar气,根据材料属性对上述材料的物性参数进行设置,具体数值如下表1所示:
(3)在真空感应炉几何模型中添加并设置磁场(mf)。
真空感应炉几何模型内的全部计算域符合安培定律,安培定律涉及到的物性参数与材料保持一致。无限元域外边界为磁绝缘。对真空感应炉几何模型的磁场初始值进行设置,真空感应炉几何模型内的所有计算域磁势为零。添加线圈节点,线圈设置为单导线线圈组,激励方式为电流激励,电流大小为250A,其余物性参数皆与材料属性保持一致。在合金计算域内添加力计算节点,设置为在频域同时计算瞬时变量和周期平均变量,并对VIM冶炼过程中的合金计算域的计算单元上的洛伦兹力进行计算。
(4)在真空感应炉几何模型中添加并设置湍流k-ε流场(spf)。
流场选择湍流k-ε(spf)模型,湍流类型采用平均纳维-斯托克斯方程(RANS)描述,采用标准壁函数处理近壁面,并将所述合金计算域内的合金熔液设置为包含重力的不可压缩流体。将合金计算域选定为流体,设定其初始速度和初始压力为零(即设置所述合金熔液的初始速度和初始压力为零),并启用静水压力补偿,湍流动能为spf.kinit m2/s2,湍流动能耗散率为spf.epinit m2/s2。将坩埚壁设置为无滑移条件(即将所述坩埚计算域的内壁面设置为无滑移壁面)。在所述湍流流场内添加重力节点,重力加速度设置为-g_const m/s2,对所述合金计算域所受的重力进行计算。在VIM冶炼过程中,除重力外,合金熔液还应受洛伦兹力影响,为计算合金计算域内的洛伦兹力,对真空感应炉几何模型进行以下设置:①添加洛伦兹力多物理场,接口为磁场和湍流流场;②在湍流物理场接口(即湍流流场)中添加体积力节点,体积力大小为洛伦兹力多物理场计算得出的时均电磁力,体积力的轴向力为mf.FLtzavZ,体积力的径向力为mf.FLtzavR。利用动网格的方法对所述合金计算域的液面波动现象进行模拟,为模拟VIM冶炼过程中合金计算域液面波动现象,对合金计算域液面进行以下设置:(a)在湍流物理场接口下,添加自由表面节点,区域为合金计算域液面;(b)将合金计算域设置为变形域动网格,初始变形量为零,网格滑移类型为超弹性滑移。
(5)在真空感应炉几何模型中添加并设置粒子跟踪物理场(fpt)。
为计算VIM冶炼过程中夹杂物粒子在合金计算域中的运动轨迹,在真空感应炉几何模型中添加流体流动颗粒跟踪(fpt)物理场接口(即粒子跟踪物理场),并将粒子视为宏观粒子。对粒子属性进行设置,模拟粒子为高温合金中的Al2O3夹杂物粒子,粒子密度为3500kg/m3,粒子直径为20μm,粒子类型为固体颗粒,电荷数为0。粒子释放时间为0s,2s,4s,6s,8s,10s,初始位置设置为随机,每次释放粒子数为60,初始速度为速度场(spf)。在合金计算域中添加重力节点以模拟粒子运动过程中重力对其造成的影响,重力加速度为-g_const m/s2。在合金计算域中添加曳力节点以模拟熔液流动对粒子运动的影响,曳力定律选择斯托克斯定律,速度为速度场(spf),动力粘度为动力粘度(spf/fpt),湍流弥散模型选择离散随机游走,湍流动能为湍流动能(spf),湍流动能耗散率为湍流动能耗散率(spf),在合金计算域中添加Saffman升力节点,速度为速度场(spf)。合金计算域液面和坩埚壁面壁条件为冻结。在本实施例中,将将曳力和Saffman升力节点处的速度全部设为速度场(spf)的原因在于,曳力和升力都与流体流动速度相关,因此需要将湍流k-ε(spf)模型计算得出的速度场带入到粒子跟踪模型中。特别说明的是,在本发明中,mf为磁场的简称,spf为流场的简称,fpt为粒子跟踪物理场的简称,特别地,对于速度场(spf)、湍流动能耗散率(spf)、湍流动能(spf)等,表示的是流场(spf)中的速度场、湍流动能耗散率、湍流动能;动力粘度为动力粘度(spf/fpt),表示的是动力粘度值受流场和粒子跟踪物理共同影响。
(6)对真空感应炉几何模型进行网格划分。
对真空感应炉几何模型进行自由三角形网格划分,其中合金计算域的网格为流体动力学极细化网格,其余计算域网格为普通物理超细化网格。将坩埚内壁面设置为边界层网格,边界层数为2层,边界层拉伸因子为1.2,厚度为0.0015m。将线圈内壁面设置为边界层网格,边界层数为2层,边界层拉伸因子为1.2,厚度为0.004m。网格划分划分结果如图3所示,完整网格包括13187个域单元,最小单元质量为0.2869,平均单元质量为0.8637;在本发明中,单元质量用来表征网格质量好坏,无单位,其值越接近于1代表网格质量越好。
(7)求解器配置并求解(计算)。
对高温合金VIM冶炼过程中夹杂物粒子(夹杂物)运动轨迹求解模型进行分步计算,其中研究1对模型流场进行求解,研究2对夹杂物粒子(夹杂物)运动轨迹和去除情况进行求解。研究1为频域-瞬态研究,时间步为rang(0, 1, 400),频率为4kHz,物理场接口选择磁场、湍流场和动网格,利用洛伦兹力多物理场进行耦合,并且在求解过程中自动对所述几何模型重新进行网格划分。研究2为瞬态研究,时间步为rang(0, 0.1, 100),物理场接口选择流体流动颗粒跟踪(即粒子跟踪物理场(fpt)),将因变量值中的不求解变量值设置为研究1计算得出的解。
(8)结果分析。
绘制熔炼时间为100s,200s,300s,400s时VIM合金计算域流场分布云图,如图4所示;采集3个探测点(在图5中记作点1、点2、点3)的速度数据并绘制探测点速度变化曲线,如图5所示。绘制VIM冶炼400s时夹杂物粒子运动25s,50s,75s和100s时的夹杂物粒子位置图,如图6所示;并对合金计算域液面和坩埚内壁面的夹杂物粒子(夹杂物)去除率进行统计,统计结果如图7所示。
本发明首次提供了一种基于COMSOL Multiphysics模拟计算VIM冶炼过程中的夹杂物运动轨迹的方法,本发明方法利用COMSOL Multiphysics软件对磁场、湍流流场和粒子跟踪物理场进行耦合,在保证其耦合的正确性下,对VIM冶炼过程中的夹杂物的运动轨迹进行仿真计算,可对密闭复杂的VIM冶炼过程中的夹杂物的运动轨迹和去除情况进行可视化仿真模拟,获得实验难以观测的结果;本发明方法能够有效保证数值模拟结果的准确性,为VIM工艺方案优化设计提供理论指导和技术支持,对实现合金材料的高纯净化目标具有重要意义。
本发明未详细说明部分为本领域技术人员公知技术。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (10)
1.一种计算VIM冶炼过程中的夹杂物运动的方法,其特征在于,所述方法包括如下步骤:
(1)利用COMSOL建立真空感应炉的几何模型;所述几何模型的计算域包括合金计算域、坩埚计算域、保温材料计算域、线圈计算域、冷凝水计算域和保护气体计算域,并将所述几何模型的最外层定义为无限元域;
(2)设置真空感应炉的材料物性参数;
(3)在所述几何模型中添加并设置磁场;所述磁场设置方式为线圈电流激励,所述线圈为单导线线圈组;设定所述无限元域的外边界为磁绝缘;对所述几何模型的磁场初始值进行设置,所述几何模型内的所有计算域的磁势均设置为零;在所述合金计算域内添加力计算节点,设置为在频域同时计算瞬时变量和周期平均变量,并对VIM冶炼过程中的合金计算域的计算单元上的洛伦兹力进行计算;
(4)在所述几何模型中添加并设置湍流流场;所述湍流流场选择湍流k-ε模型对所述合金计算域的流场进行模拟分析,湍流类型采用平均纳维-斯托克斯方程描述,采用标准壁函数处理近壁面,利用动网格的方法对所述合金计算域的液面波动现象进行模拟;将所述坩埚计算域的内壁面设置为无滑移壁面;将所述合金计算域内的合金熔液设置为包含重力的不可压缩流体,设置所述合金熔液的初始速度和初始压力为零,并启用静水压力补偿;在所述湍流流场内添加重力节点,对所述合金计算域所受的重力进行计算;并利用洛伦兹力多物理场对所述合金计算域所受的时均电磁力进行计算,且将时均电磁力作为体积力输入到湍流k-ε模型中;
(5)在所述几何模型中添加并设置粒子跟踪物理场;所述粒子为随机生成的夹杂物粒子,所述夹杂物粒子的初始速度为所在位置的合金熔液的流动速度;对所述夹杂物粒子的密度、直径、类型及电荷数进行设置,设置所述坩埚计算域的内壁面和所述合金计算域的液面的壁条件为冻结或反弹,以实现夹杂物粒子去除效果的仿真模拟;然后在所述合金计算域中添加重力节点、曳力节点和萨夫曼升力节点,对所述合金计算域中的夹杂物粒子在运动过程中的受力情况进行仿真模拟;
(6)对所述几何模型进行网格划分;
(7)配置求解器并求解;求解的过程依次包括求解模型流场和求解夹杂物粒子的运动轨迹及去除情况;
(8)结果分析;所述结果分析包括绘制VIM冶炼过程中的合金计算域流场分布云图和速度随时间的变化曲线、绘制VIM冶炼过程中的夹杂物粒子的运动轨迹图以及对合金计算域的液面和坩埚计算域的内壁面的夹杂物粒子去除率进行统计。
2.根据权利要求1所述的方法,其特征在于:
在步骤(1)中,建立所述几何模型所需的参数包括:炉料高度;坩埚的高度、宽度和厚度;保温材料的高度、宽度和厚度;线圈的内径和外径,线圈间的间隙及匝数。
3.根据权利要求1所述的方法,其特征在于:
在步骤(2)中,所述真空感应炉的材料包括合金、坩埚、保温材料、线圈、冷凝水和保护气体;
所述材料物性参数包括所述合金、所述坩埚、所述保温材料、所述线圈、所述冷凝水和所述保护气体的相对磁导率、电导率、相对介电常数以及所述合金的密度和动力粘度。
4.根据权利要求3所述的方法,其特征在于:
所述坩埚采用氧化镁、氧化钙和氧化铝中的一种或多种材料制成;
所述保温材料为硅酸铝纤维材料、石棉材料、高铝砖材料中的一种或多种;
所述线圈采用紫铜材料制成;
所述冷凝水为水;和/或
所述保护气体为Ar气。
5.根据权利要求1所述的方法,其特征在于:
在步骤(3)中,所述几何模型内的所有计算域符合安培定律。
6.根据权利要求1所述的方法,其特征在于:
在步骤(4)中,将所述重力节点的加速度设置为当地重力加速度;
在步骤(4)中,在利用洛伦兹力多物理场对所述合金计算域所受的时均电磁力进行计算时,所述洛伦兹力多物理场的接口为磁场和湍流流场,在所述湍流流场中添加体积力节点,体积力大小为时均电磁力,体积力的轴向力为z方向时均电磁力,体积力的径向力为r方向时均电磁力。
7.根据权利要求1所述的方法,其特征在于:
在步骤(4)中,为模拟所述合金计算域的液面波动现象,对所述合金计算域的液面进行以下设置:
(a)在湍流流场中添加自由表面节点,区域设置为合金计算域液面;
(b)将所述合金计算域设置为变形域动网格,初始变形量设置为零,网格滑移类型设置为超弹性滑移。
8.根据权利要求1所述的方法,其特征在于:
在步骤(6)中,对所述几何模型进行自由三角形网格划分;所述合金计算域的网格相比所述坩埚计算域、所述保温材料计算域、所述线圈计算域、所述冷凝水计算域和所述保护气体计算域的网格较细,并将所述坩埚计算域的内壁面和所述线圈计算域的内壁面设置为边界层网格。
9.根据权利要求8所述的方法,其特征在于:
所述合金计算域的网格为流体动力学极细化网格,所述坩埚计算域、所述保温材料计算域、所述线圈计算域、所述冷凝水计算域和所述保护气体计算域的网格为普通物理超细化网格。
10.根据权利要求1所述的方法,其特征在于,在步骤(7)中:
求解模型流场的过程为频域-瞬态研究,物理场的接口为磁场、湍流流场和动网格,利用洛伦兹力多物理场进行耦合,并且在求解过程中自动对所述几何模型重新进行网格划分;
求解夹杂物粒子的运动轨迹及去除情况的过程为瞬态研究,物理场的接口为粒子跟踪物理场,将因变量值中的不求解变量值设置为求解模型流场得出的解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110114283.2A CN112464543B (zh) | 2021-01-28 | 2021-01-28 | 一种计算vim冶炼过程中的夹杂物运动的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110114283.2A CN112464543B (zh) | 2021-01-28 | 2021-01-28 | 一种计算vim冶炼过程中的夹杂物运动的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112464543A true CN112464543A (zh) | 2021-03-09 |
CN112464543B CN112464543B (zh) | 2021-04-06 |
Family
ID=74802437
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110114283.2A Active CN112464543B (zh) | 2021-01-28 | 2021-01-28 | 一种计算vim冶炼过程中的夹杂物运动的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112464543B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113779906A (zh) * | 2021-11-12 | 2021-12-10 | 北京科技大学 | 镍基高温合金真空感应熔炼过程中脱氮动力学的计算方法 |
CN114850465A (zh) * | 2022-06-15 | 2022-08-05 | 北京科技大学 | 一种钢水可浇性预测系统和方法 |
CN115125396A (zh) * | 2022-07-28 | 2022-09-30 | 中国科学院金属研究所 | 采用var工艺对特殊钢枝晶细化及夹杂物去除的控制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1395196A (zh) * | 2002-08-08 | 2003-02-05 | 上海交通大学 | 过程冶金中计算流体力学的数据可视化方法 |
US20150213163A1 (en) * | 2012-12-20 | 2015-07-30 | Institute Of Modern Physics, Chinese Academy Of Sciences | Particle flow simulation system and method |
CN107301291A (zh) * | 2017-06-21 | 2017-10-27 | 河北工业大学 | 一种基于电渣重熔过程的电磁场、温度场和流动场耦合有限元计算方法 |
CN109583066A (zh) * | 2018-11-22 | 2019-04-05 | 南京工程学院 | 一种直流架空线路绝缘子表面污秽沉积模拟方法 |
CN111950137A (zh) * | 2020-07-28 | 2020-11-17 | 北京科技大学 | 用于固体夹杂物在钢液中运动轨迹的模拟方法及系统 |
-
2021
- 2021-01-28 CN CN202110114283.2A patent/CN112464543B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1395196A (zh) * | 2002-08-08 | 2003-02-05 | 上海交通大学 | 过程冶金中计算流体力学的数据可视化方法 |
US20150213163A1 (en) * | 2012-12-20 | 2015-07-30 | Institute Of Modern Physics, Chinese Academy Of Sciences | Particle flow simulation system and method |
CN107301291A (zh) * | 2017-06-21 | 2017-10-27 | 河北工业大学 | 一种基于电渣重熔过程的电磁场、温度场和流动场耦合有限元计算方法 |
CN109583066A (zh) * | 2018-11-22 | 2019-04-05 | 南京工程学院 | 一种直流架空线路绝缘子表面污秽沉积模拟方法 |
CN111950137A (zh) * | 2020-07-28 | 2020-11-17 | 北京科技大学 | 用于固体夹杂物在钢液中运动轨迹的模拟方法及系统 |
Non-Patent Citations (1)
Title |
---|
汪瑞婷: "电渣重熔过程夹杂物运动行为的数值模拟", 《钢铁研究学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113779906A (zh) * | 2021-11-12 | 2021-12-10 | 北京科技大学 | 镍基高温合金真空感应熔炼过程中脱氮动力学的计算方法 |
CN114850465A (zh) * | 2022-06-15 | 2022-08-05 | 北京科技大学 | 一种钢水可浇性预测系统和方法 |
CN115125396A (zh) * | 2022-07-28 | 2022-09-30 | 中国科学院金属研究所 | 采用var工艺对特殊钢枝晶细化及夹杂物去除的控制方法 |
CN115125396B (zh) * | 2022-07-28 | 2023-12-01 | 中国科学院金属研究所 | 采用var工艺对特殊钢枝晶细化及夹杂物去除的控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112464543B (zh) | 2021-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112464543B (zh) | 一种计算vim冶炼过程中的夹杂物运动的方法 | |
Kharicha et al. | Review on modeling and simulation of electroslag remelting | |
Li et al. | 3D numerical simulation of successive deposition of uniform molten Al droplets on a moving substrate and experimental validation | |
Tang et al. | Numerical and experimental investigation of laser additive manufactured Ti2AlNb-based alloy | |
Karimi-Sibaki et al. | A parametric study of the vacuum arc remelting (VAR) process: Effects of arc radius, side-arcing, and gas cooling | |
CN109165469B (zh) | 方坯连铸结晶器电磁搅拌参数模拟设定及优化方法 | |
Li et al. | Modeling of electromagnetic field and liquid metal pool shape in an electroslag remelting process with two series-connected electrodes | |
Pericleous et al. | A multiscale 3D model of the vacuum arc remelting process | |
Yuan et al. | Development of an improved CBR model for predicting steel temperature in ladle furnace refining | |
Melissari et al. | Development of a heat transfer dimensionless correlation for spheres immersed in a wide range of Prandtl number fluids | |
Wang et al. | Numerical analysis of inclusion motion behavior in electroslag remelting process | |
Ghazal et al. | On the dissolution of nitrided titanium defects during vacuum arc remelting of ti alloys | |
Reddy et al. | Evaluation of glass forming ability of Zr–Nb alloy systems through liquid fragility and Voronoi cluster analysis | |
Xu et al. | Tree-ring formation during vacuum arc remelting of INCONEL 718: Part II. Mathematical modeling | |
Tsai et al. | Numerical simulation of solidification morphologies of Cu-0.6 Cr casting alloy using modified cellular automaton model | |
Nikulin | MATHEMATICAL MODELLING OF AMF GEOMETRY AND FREQUENCY IMPACTS ON VOLUME AND SURFACE MELT FLOWS AT INDUCTION MELTING. | |
CN115828697A (zh) | 一种电子束冷炉床熔炼过程中电子束对铸锭凝固影响数值模拟方法 | |
Wang et al. | Numerical investigation on species transport in electroslag remelting dual alloy ingot | |
Fallah-Mehrjardi et al. | Investigation of freeze lining in copper-containing slag systems: part III. high-temperature experimental investigation of the effect of bath agitation | |
Li et al. | Migration behavior of solidification nuclei in pure Al melt under effect of electric current pulse | |
Liu et al. | Comprehensive model for a slag bath in electroslag remelting process with a current-conductive mould | |
Wang et al. | Numerical investigation on the effect of slag thickness on metal pool profile in electroslag remelting process | |
Zhang et al. | Break the superheat temperature limitation of induction skull melting technology | |
Shi et al. | Insight into droplet formation in electroslag remelting process by numerical simulation | |
Song et al. | An electromagnetic and thermal analysis of a cold crucible melting |
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 |