计算流体动力学和动态耦合热力学毕业论文中英文资料外文翻译文献

计算流体动力学和动态耦合热力学软件

在顶吹转炉中的应用

Mikael ERSSON, Lars HÖGLUND, Anders TILLIANDER,

Lage JONSSON and Pär JÖNSSON

应用程序冶金部，皇家技术学院（KTH）

SE-10044，瑞典首都斯德哥尔摩。

（2007年11月8日收到，2007年12月10日接受）

一种新的建模方法被提出，这种建模方法是使用计算流体力学软件相连结热力学数据库以获取动态模拟冶金过程的现象。这种建模方法已被应用在一个基本的氧气顶吹转炉模型。通过各种气体之间的反应研究。结果表明,大量的表面气体的流通是完全受对流控制的。此外，在这个过程中大量产生的CO脱碳可能会放慢从浴缸喷出的液滴的脱碳率。在目前的模拟反映实验室的实验条件下，这点也被证实在这个过程中所产生的炉渣（FeO和/或SiO2）接近于零，即只产生的气体（二氧化碳CO2）就好比是氧气射流击中钢液。它也说明如何从几秒钟的采样推算脱碳速度，只要是含碳量足够的高可以在后期的时候做含碳量的模拟，从而得到的碳含量的粗略估计。总的结论是，通过Thermo-Calc的数据库和CFD软件的动态耦合来达到冶金动态模拟是有可能的。

关键词：转炉;计算流体力学，热力学建模;炉渣和动态模拟。

介绍

在许多涉及氧气喷射撞击到钢液面的冶金过程中，为了优化涉及动力学的部分，如脱碳，底层流体动力学是需要的。现在对于这个问题已经几个实验报告和一些数值或计算流体动力学（CFD）的报告。Szekely and Asai已经介绍了一种将液体的冲击射流表面的计算模型。Ngyen and Evans通过使用这种方法计算溶池喷嘴直径比液体表面所造成变形的冲击射流的影响，张等人模仿了一种同时使用顶吹和底吹复合吹炼的情况。Odenthal等人展示了一种顶吹转炉

多相CFD模型，在这种顶吹转炉中由于冲击射流以及底部和顶部转换混合时存在飞溅现象。Nakazono等人描述了铁液表面的超音速氧气喷射冲击时的含碳量的两阶段的数值分析。通过计算表明在真空和表面处理条件下天然气和钢铁之间的变化。该模型采用稳态方法无飞溅等。在其他文献中还其他非顶吹CFD模型介绍，琼森等。提出了硫精炼耦合计算流体力学和热力学模型

热力学是成立的在CFD程序中作为一个自定义的子程序，尤其是专门为调查系统所写的子程序。图中1可以看到这样的做法一个示意图。

图1。合并计算流体力学和热力学建模与数值模型方法的示意图。由文献（14）

最近包含顶吹系统CFD模型已经出现并和实验数据进行比较，这里提出，这个模型是包括气/液/渣的流体体积（VOF）等反应扩展的一个多相流模型，这个模型适用于顶吹系统。这个模型以及各个阶段的扩展已经获得批准使用这一方法所建的装备使得终于有一个可以说明基本的顶吹转炉模型结果。

2。数值模型

目前的建模方法如图2所示。

图2。合并数值模型方法的示意图

为了解决新的研究可以轻易的合并到这个模型的问题建立了以模块化的方式，这也意味着，当从一个系统变到另一个系统是只需要一个很小的重编程程序，只要是热力学数据在目前的热力学数据库中是存在的以及不超过CFD软件的运算能力。这个模拟过程已经在包含6个节点的集群的Linux PC上执行。现实生活中的模拟时间是和运算系统的数量相关的，所以没有固定和典型的模拟时间;它可能会在一小时和10天的实时变化亦或是10秒的时间。在图3中可以看出其物理域和数值域的原理。

图3。顶吹转炉示意图。a）物理域b）数值域。

所涉及的边界条件为入口速度，出口压力，无滑墙壁和对称轴。所有的墙壁都用标准壁面函数表示，合理的k–e模型已使用在所有的例子中。域的宽度为0.075米，高度是0.13米。最开始，1500°C的15.6公斤距顶枪0.01米以上的钢熔液用来试验，大量的气体从入口吹入相当于25升/分钟的纯体积流量氧气。用于模拟的表1中的各种参数可以看出和表2中的初始浓度的不同。

表1.不同阶段密度和扩散系数

表2.初始浓度

2.1。差价

ANSYS软件同时也被使用，这是一个为计算求解流体体积、质量的有限元分析软件，

这样动量守恒方程中所涉及的质量、体积都得到了解决。根据湍流模型同时使用了一些额外的守恒方程作为补充,例如湍流动能守恒,k,和动荡能量耗散,e,都是用标准的k-e 模型，下面的公式用来计算任何形式的

这里r是密度,u是指速度矢量,当基于雷诺兹平均使用湍流模型,G是扩散系数,正如从表3中看到的一样方程(1)和表3描述了质量、动量、湍流动能,动荡能量耗散、能量、物质和体积分数。

表3守恒方程参数

2.2热计算

为了获得准确描述了热力学特征，软件Thermo-Calc被使用使用。这是一个通用的软件方案的多组分的相平衡计算。它使用一种技术,它允许非常灵活的设置条件的平衡状态,从而适用过程模拟。

问题的解决方法是以下。第一,质量和热量含量分别计算每个阶段。然后,总质量和热量的内容是总结。该系统是氢原子核装备中的。最后,程序将计算各阶段温度,新的成分和含量。

能与热力学的软件应用程序编程接口使用TQ操作。这个接口是一个接口

Thermo-Calc中可用的软件包,使其能要推广实施不同系统(例如组件和元素,该系统由),而无需更改代码。

2.3. 耦合和假设

其目在于CFD-package和Thermo-Calc数据库软件两个软件的耦合,是创建一个通用的计算模型,包括化学反应冶金系统。在下面的文本描述的耦合。主要的假设是局部均衡、每个时间步可以到达每个计算单元的过程中,。软件之间的接口CFD-package和热力学数据库分别被编码在C和FORTRAN。

2.4．多相考虑

所有的散装阶段的建模为不可压缩;有统一的密度,见表1。与一个精化的模型应该有可能使用一个理想气体定律假设对气态和一些温度对钢密度模型的依赖。简要地说明一些细节的功能假设一摩尔的氧和碳反应形成的融化2摩尔的一氧化碳。在计算单元在考虑将会有一个扩张的气相到约两倍于原体积。这反过来将很可能意味着气体向外扩张的计算单元。如果再加上CFD,扩张将发生在以下步骤中,作为额外的气体质量被添加到细胞源项。经过计算细胞已经达到平衡它也将会有一个特定的温度平衡温度。

3.结论

3.1.基本顶级吹转炉

在图5显示气体射流向量的情况,正如水落在钢铁溶液中所显示的那样。这是看到射流失去其轴向冲击就好像其滴落在刚溶液表面。从射流中来的低流量给出了一个相对较小的渗透率在钢铁溶液中。

图6说明了流场的钢铁浴引起碰撞射流。目前的最高速度的射流冲击的钢液面和面积。流体措施阻止了渗透区向墙壁,然后进入到钢液中。一个大型的循环,循环中形成的冲击流(即轴和墙之间)和几个小的循环回路。应该指出的是,在更大的循环级速度是非常小的。这是由于低流量来自高层的射击流。更高的流率呈现较大的渗透和较强的循环。

图5.在气体速度矢量。矢量有固定长度;他们是由色彩反映速度级(m / s)。