D0I:10.13374/1.issnl00103.2009.07.036 第31卷第7期 北京科技大学学报 Vol.31 No.7 2009年7月 Journal of University of Science and Technology Beijing Jd.2009 多流中间包流动特征的数学模型 潘宏伟程树森 北京科技大学治金与生态工程学院，北京100083 摘要提出了针对多流中间包流动特征的数学模型，给出了绘制停留时间分布(RTD)曲线的标准及其相关计算的规范：明 确了“死区“的特性，并计算了其体积分数；模型中考虑到了短路流，指出计算短路流的方法·模型可以更好地刻画活塞流的流 动特性 关键词中间包：停留时间分布：短路流：死区：水模型：数值模拟 分类号TF777.2 Mathematical model of flow characterization in multi-strand continuous casting tundishes PAN Hong"wei,CHENG Shu-sen School of Metallurgical and Ecological Engineering University of Science and Technology Beijing.Beijing 100083.China ABSTRACT A mathematical model of flow characterization in multi-strand continuous casting tundishes was attempted to put for- ward,the standard of plotting residual time distribution (RTD)curves was definitely presented and a relative calculation method was proposed.The feature of dead region was carefully studied and a method to calculate its volume fraction was given in the model. What's more,the model considers bypass flow and proposes a method to calculate its volume fraction.Experimental results show that this new model can better capture the difference between the reality flow pattern and the ideal plug flow pattern and reflects the fea- ture of plug flow· KEY WORDS tundish:residual time distribution (RTD):bypass flow dead region:water model:mathematical simulation 中间包是连续操作的反应器，也是高质量钢材 1O0%.文献[9]在Sahai提出的模型基础上，将单流 生产流程中关键一环可].按流动方式的不同，中 中间包的数学模型引用到多流中间包，除此之外， 间包内可分为短路流区、活塞流区、返混流区和死 多流中间包流动特征的数学模型研究很少见到报 区，四个区域的体积分数反映了中间包内的流动特 道。经过仔细分析，发现现有数学模型仍存在亟待 征，本文研究中间包流动特征的数学模型的目的就 完善之处：①模型中没有提出与之紧密相关的停留 是要建立合理的数学表达式，根据表达式计算出中 时间分布(RTD)曲线的绘制标准，使冶金工作者在 间包内各区域的体积分数.Sahai在文献[G-7]提出 这一问题上认识不统一；②模型中没有涉及短路流， 了基于单流中间包的流动特征数学模型：此后国内 而短路流显著影响着中间包治金工艺：③在计算死 外治金工作者常常借用该模型计算多流中间包问 区体积分数时将死区视为绝对的静止，这不符合实 题，但计算结果不尽如人意，各流股对应的相同区域际情况，在一些情况下利用Saai的模型计算得到 体积分数之和往往偏大，甚至大于100%，这是不符 死区的体积分数为负，文献[10]中对此进行了陈述； 合实际情况的，且流股越多，和的数值越大，如文献 ④模型不能全面体现活塞流的流动特性，因此，进 [8]中各流股对应的返混区体积分数之和大于 一步研究多流中间包流动特征的数学模型，完善多 收稿日期：2008-07-31 基金项目：国家自然科学基金资助项目(N。.60872147) 作者简介：潘宏伟(1984-)，男，博士研究生，E mail:pan5025@163.com:程树森(1964一)，男，教授，博士生导师
多流中间包流动特征的数学模型 潘宏伟 程树森 北京科技大学冶金与生态工程学院北京100083 摘 要 提出了针对多流中间包流动特征的数学模型给出了绘制停留时间分布（RTD）曲线的标准及其相关计算的规范；明 确了“死区”的特性并计算了其体积分数；模型中考虑到了短路流指出计算短路流的方法．模型可以更好地刻画活塞流的流 动特性． 关键词 中间包；停留时间分布；短路流；死区；水模型；数值模拟 分类号 TF777∙2 Mathematical model of flow characterization in mult-i strand continuous casting tundishes PA N Hong-weiCHENG Shu-sen School of Metallurgical and Ecological EngineeringUniversity of Science and Technology BeijingBeijing100083China ABSTRACT A mathematical model of flow characterization in mult-i strand continuous casting tundishes was attempted to put forwardthe standard of plotting residual time distribution （RTD） curves was definitely presented and a relative calculation method was proposed．T he feature of dead region was carefully studied and a method to calculate its volume fraction was given in the model． What’s morethe model considers bypass flow and proposes a method to calculate its volume fraction．Experimental results show that this new model can better capture the difference between the reality flow pattern and the ideal plug flow pattern and reflects the feature of plug flow． KEY WORDS tundish；residual time distribution （RTD）；bypass flow；dead region；water model；mathematical simulation 收稿日期：2008-07-31 基金项目：国家自然科学基金资助项目（No．60872147） 作者简介：潘宏伟（1984—）男博士研究生E-mail：pan5025＠163．com；程树森（1964—）男教授博士生导师 中间包是连续操作的反应器也是高质量钢材 生产流程中关键一环［1—5］．按流动方式的不同中 间包内可分为短路流区、活塞流区、返混流区和死 区．四个区域的体积分数反映了中间包内的流动特 征．本文研究中间包流动特征的数学模型的目的就 是要建立合理的数学表达式根据表达式计算出中 间包内各区域的体积分数．Sahai 在文献［6—7］提出 了基于单流中间包的流动特征数学模型；此后国内 外冶金工作者常常借用该模型计算多流中间包问 题但计算结果不尽如人意各流股对应的相同区域 体积分数之和往往偏大甚至大于100％这是不符 合实际情况的且流股越多和的数值越大如文献 ［8］ 中各流股对应的返混区体积分数之和大于 100％．文献［9］在 Sahai 提出的模型基础上将单流 中间包的数学模型引用到多流中间包．除此之外 多流中间包流动特征的数学模型研究很少见到报 道．经过仔细分析发现现有数学模型仍存在亟待 完善之处：①模型中没有提出与之紧密相关的停留 时间分布（RTD）曲线的绘制标准使冶金工作者在 这一问题上认识不统一；②模型中没有涉及短路流 而短路流显著影响着中间包冶金工艺；③在计算死 区体积分数时将死区视为绝对的静止这不符合实 际情况在一些情况下利用 Sahai 的模型计算得到 死区的体积分数为负文献［10］中对此进行了陈述； ④模型不能全面体现活塞流的流动特性．因此进 一步研究多流中间包流动特征的数学模型完善多 第31卷 第7期 2009年 7月 北 京 科 技 大 学 学 报 Journal of University of Science and Technology Beijing Vol．31No．7 Jul．2009 DOI:10．13374／j．issn1001－053x．2009．07．036
,816 北京科技大学学报 第31卷 流中间包内各区域体积分数的计算十分必要， 1200 820 1水模型实验 直接对高温下的钢液流场进行测量十分困难， 而用水模型来模拟钢液流动被证明是可行的，并能 中部水口边部水口 够正确反映实际钢液流动的规律.根据相似理论， 1990200 980 当中间包内钢液流动与模型中流体流动处于同一自 2240 5860 模化区时，只要保证原型和模型的F相等，就可保 证模型和原型的相似，所以模型采用F,准数做定 图1中间包几何尺寸示意图（单位：mm) 性准则，根据Fr准数相等可以得到Fr。=Frm;即 Fig.I Schematic diagram of a tundish (in mm) 中间包控流参数的优劣，本次水模型实验示踪剂采 上-正取比例因子=乞= 。3则有： 用饱和的KC溶液，用电导仪测量中间包出水口的 Vm=JλV。=0.577Vp, 电压随时间变化曲线，由于中间包左、右两侧对称， 实验中仅需测得中间包右半部分的出水口电压随时 Qm=λ2Qp=0.06415Qp 间的变化，如图2所示.引用单流中间包流动特征 式中，Vp和Vm分别为原型和模型的拉速，Frp和 的数学模型计算出死区体积，结果如表1所示，从计 Frm分别为原型和模型的弗劳德准数 算结果中可以看出，该中间包右侧两流股（中部和边 基于以上分析以及板坯连铸中间包的实际尺 部)的死区体积分数之和大致为30%~54%，如果 寸，按模型与实物1：3的比例制作了有机玻璃材质 将中间包左、右两侧各流股视为对称，整个中间包内 的模型.中间包原形状和尺寸如图1所示， 各流股对应的死区体积分数之和将为60%~ 在水模型实验中，利用“刺激一响应”法来测量 108%,表明一大半的中间包处于“死寂”状态，这在 注入的示踪剂在中间包模型内的停留时间分布曲 实际中是不可能的，死区体积超过100%更加是不 线，定量地描述中间包内总体的流动特性，最终评价 可能的 16 20 (b) 12 10 200 400 60080010001200 200 400 60080010001200 时间s 时间s 图2水模型实验电压随时间变化结果()中部水口电压变化：(b)边部水口电压变化 Fig2 Water model experimental resu for time varied voltage:(a)time varied volage for the inner strand:(b)time varied voltage for the outer strand 表1不同入口速度条件下水模型计算结果 模型引用到多流中间包后，计算出的死区体积分数 Table 1 Calculation results for different speeds at entrance in water 明显偏大，各流对应的体积分数之和甚至超过 model experiment 100%,明显不符合实际情况 0/ 水口 类别 Lmin/s 7ls t/s (mmi凸位置 /% 2数值模拟结果分析 中部 39 499.6 558.6 10.6 优化后 0.85 边部 40 449.9 558.6 19.5 针对某厂另外两种不同设计理念的中间包进行 中部 36 407 478.5 14.9 了数值模拟，分别命名为1#中间包和2#中间包 优化后 1.0 边部 40 407.1478.5 14.9 为减小数值模拟计算量并便于与水模型相互验证， 中部 14 374.9558.6 32.9 数值模拟中中间包尺寸为原模型的1/3，流体为水. 原型 0.85 边部 75 437.3558.621.7 常用的微分方程和边界条件广泛见于各类文献，不 再赘述，本文着重讨论计算结果， 由此可见，直接将单流中间包流动特征的数学 图3(a)描述了1中间包流场中的流线（顶部
流中间包内各区域体积分数的计算十分必要． 1 水模型实验 直接对高温下的钢液流场进行测量十分困难 而用水模型来模拟钢液流动被证明是可行的并能 够正确反映实际钢液流动的规律．根据相似理论 当中间包内钢液流动与模型中流体流动处于同一自 模化区时只要保证原型和模型的 Fr 相等就可保 证模型和原型的相似．所以模型采用 Fr 准数做定 性准则．根据 Fr 准数相等可以得到 Frp＝ Frm；即 V 2 p gLp ＝ V 2 m gL m ；取比例因子 λ＝ L m Lp ＝ 1 3 ；则有： V m＝ λV p＝0∙577V p Qm＝λ 5 2 QP＝0∙06415QP． 式中V p 和 V m 分别为原型和模型的拉速Frp 和 Frm 分别为原型和模型的弗劳德准数． 基于以上分析以及板坯连铸中间包的实际尺 寸按模型与实物1∶3的比例制作了有机玻璃材质 的模型．中间包原形状和尺寸如图1所示． 在水模型实验中利用“刺激—响应”法来测量 注入的示踪剂在中间包模型内的停留时间分布曲 线定量地描述中间包内总体的流动特性最终评价 图1 中间包几何尺寸示意图（单位：mm） Fig．1 Schematic diagram of a tundish （in mm） 中间包控流参数的优劣．本次水模型实验示踪剂采 用饱和的 KCl 溶液．用电导仪测量中间包出水口的 电压随时间变化曲线．由于中间包左、右两侧对称 实验中仅需测得中间包右半部分的出水口电压随时 间的变化如图2所示．引用单流中间包流动特征 的数学模型计算出死区体积结果如表1所示从计 算结果中可以看出该中间包右侧两流股（中部和边 部）的死区体积分数之和大致为30％～54％如果 将中间包左、右两侧各流股视为对称整个中间包内 各流股对应的死区体积 分 数 之 和 将 为 60％ ～ 108％表明一大半的中间包处于“死寂”状态这在 实际中是不可能的死区体积超过100％更加是不 可能的． 图2 水模型实验电压随时间变化结果．（a）中部水口电压变化；（b）边部水口电压变化 Fig．2 Water-model experimental results for time-varied voltage：（a） time-varied voltage for the inner strand；（b） time-varied voltage for the outer strand 表1 不同入口速度条件下水模型计算结果 Table1 Calculation results for different speeds at entrance in watermodel experiment 类别 v／ （m·min —1） 水口 位置 tmin／s t／s τ／s V d V ／％ 优化后 0．85 中部 39 499．6 558．6 10．6 边部 40 449．9 558．6 19．5 优化后 1．0 中部 36 407 478．5 14．9 边部 40 407．1 478．5 14．9 原型 0．85 中部 14 374．9 558．6 32．9 边部 75 437．3 558．6 21．7 由此可见直接将单流中间包流动特征的数学 模型引用到多流中间包后计算出的死区体积分数 明显偏大各流对应的体积分数之和甚至超过 100％明显不符合实际情况． 2 数值模拟结果分析 针对某厂另外两种不同设计理念的中间包进行 了数值模拟分别命名为1＃ 中间包和2＃ 中间包． 为减小数值模拟计算量并便于与水模型相互验证 数值模拟中中间包尺寸为原模型的1／3流体为水． 常用的微分方程和边界条件广泛见于各类文献不 再赘述本文着重讨论计算结果． 图3（a）描述了1＃ 中间包流场中的流线（顶部 ·816· 北 京 科 技 大 学 学 报 第31卷
第7期 潘宏伟等：多流中间包流动特征的数学模型 ,817 视角)从中可以看到，流体流出挡墙上的溢流孔 图4(a)描述了2#中间包内流线图（前方视 后，左、右两侧各产生了一大一小两个漩涡.图3(b) 角)·从中可以看到，流体流出溢流孔后，斜向上流 是流体穿过溢流孔水平面的速度矢量图，箭头方向 动，再沿着上部的自由液面流动，然后向下流向出 表示速度方向，箭头长度表示其大小，就该图来看， 口，是一种强烈的活塞流流动趋势，图4(b)显示其 中间包内并无大面积死区，存在较大范围的返混区， 穿过溢流孔后的水平速度场，该图显示整体流动比 (a) 浸人式 (b) 挡墙上的 0.0093m+s 入口 挡墙上的 溢流孔 溢流孔 第1流出口 第2流出口 第3流出口第4流出口 图31°中间包内流线示意图()和穿过溢流孔水平速度矢量图(b) Fig.3 Schematic diagram of streamlines in Tundish No.1(a)and velocity field at the horizontal plane through the spout hole (b) a 挡墙 一浸人式水口 (b) 挡墙上的 0.0093m.s ,导流杯 滋流孔 第1流出口第2流出口第3流出口第4流出口 图42中间包内流线示意图(a)和穿过溢流孔水平速度失量图(b) Fig.4 Schematic diagram of streamlines for Tundish No.2 (a)and velocity field at the horizontal plane through the spout hole No.2 (b) 较均匀，没有明显死区， 算结果如图5和表2所示.从中可以看到两个缺 鉴于上述两种中间包的对称性，仅讨论中间包 点：其一，2#中间包内死区体积分数计算结果出现 的右半部分，即第3流（靠近中部）和第4流（靠近边 负值（一3.6%），在实际中这是不可能的；其二，各流 部) 对应的返混区体积分数之和超过100%，这也是不 200 合理的 ☑☑☑1'中间包 160 因2”中间包 表21产和2*中间包内各区域体积分数计算结果 Table 2 Calculation results of different regions'volume fraction for 20 第3流对应的死区 Tundishes No.I and No.2 体积分数的计算结果 为负值，导致第3流 中间包 体积分数 第3流 第4流 总计 80 和第4流死区体积 分数之和几乎为零 Valv 2.35 1.49 3.8 1*中间包 Vi/v 79.3 80.4 159.7 VV 18.3 18.1 36.4 V/V Valv -3.6 3.62 0.02 2*中间包 Vi/v 77.6 73.6 151.1 图51*和2#中间包内各区域体积分数柱状图 V/V 26.1 22.8 48.9 Fig.5 Histogram of different regions'volume fraction for both Tundishes No.I and No.2 另外，根据Sahail提出的数学模型中假设死区 如果直接引用Sahai提出的单流中间包数学模 是绝对静止的，则中间包实际有效容积为V一V, 型来计算此多流中间包内各区域的体积分数，其计 进而可得：
视角）．从中可以看到流体流出挡墙上的溢流孔 后左、右两侧各产生了一大一小两个漩涡．图3（b） 是流体穿过溢流孔水平面的速度矢量图箭头方向 表示速度方向箭头长度表示其大小．就该图来看 中间包内并无大面积死区存在较大范围的返混区． 图4（a）描述了2＃ 中间包内流线图（前方视 角）．从中可以看到流体流出溢流孔后斜向上流 动再沿着上部的自由液面流动然后向下流向出 口是一种强烈的活塞流流动趋势．图4（b）显示其 穿过溢流孔后的水平速度场该图显示整体流动比 图3 1＃中间包内流线示意图（a）和穿过溢流孔水平速度矢量图（b） Fig．3 Schematic diagram of streamlines in Tundish No．1（a） and velocity field at the horizontal plane through the spout hole （b） 图4 2＃中间包内流线示意图（a）和穿过溢流孔水平速度矢量图（b） Fig．4 Schematic diagram of streamlines for Tundish No．2（a） and velocity field at the horizontal plane through the spout hole No．2（b） 较均匀没有明显死区． 鉴于上述两种中间包的对称性仅讨论中间包 的右半部分即第3流（靠近中部）和第4流（靠近边 部）． 图5 1＃和2＃中间包内各区域体积分数柱状图 Fig．5 Histogram of different regions’volume fraction for both Tundishes No．1and No．2 如果直接引用 Sahai 提出的单流中间包数学模 型来计算此多流中间包内各区域的体积分数其计 算结果如图5和表2所示．从中可以看到两个缺 点：其一2＃中间包内死区体积分数计算结果出现 负值（—3∙6％）在实际中这是不可能的；其二各流 对应的返混区体积分数之和超过100％这也是不 合理的． 表2 1＃和2＃中间包内各区域体积分数计算结果 Table2 Calculation results of different regions’volume fraction for Tundishes No．1and No．2 中间包 体积分数 第3流 第4流 总计 V d／V 2∙35 1∙49 3∙8 1＃中间包 V ba／V 79∙3 80∙4 159∙7 V p／V 18∙3 18∙1 36∙4 V d／V —3∙6 3∙62 0∙02 2＃中间包 V ba／V 77∙6 73∙6 151∙1 V p／V 26∙1 22∙8 48∙9 另外根据 Sahai ［7］提出的数学模型中假设死区 是绝对静止的则中间包实际有效容积为 V — V d 进而可得： 第7期 潘宏伟等： 多流中间包流动特征的数学模型 ·817·
,818 北京科技大学学报 第31卷 0.=vo,0'=0-0m 数输入，其输入总量为： v/Q M=Qco△t'. 模型中也不考虑短路流，则： 式中，M为示踪剂在入口的输入总量；c0为示踪剂 Q=0,Q'=Q 输入浓度；△t为示踪剂注入时间，为满足6函数， 所以在文献[7]中得到死区的体积分数为： 相对整个计算时间该值应该非常小；Q为流体在入 1-0=1-w-ya/0=ya v/Q 口的体积流量， V 反应器 但是，在实际冶金过程中，死区的定义是一个相对的 死区 概念：超过平均停留时间2倍以上的区域，并不是绝 人口 出口 对静止的.因而，中间包实际有效容积并不能等效 流体 活塞区 返混区 为V一Va,文献[7]中计算死区的体积分数也就不 短路流 合理， 综上所述，水模型实验和数值模拟结果均显示， 采用Sahai的针对单流的数学模型来计算多流中间 图6本文模型示意图 包内各区域体积分数存在不足，因而利用其结果指 Fig.6 Schematic diagram of the proposed model 导多流中间包优化设计和多流中间包冶金是不可靠 示踪剂在中间包出口的输出总量为： 的 w= Qici(ti)dt 3提出新模型 根据质量守恒，总的输出量等于总的输入量： 一个合理的计算模型应该计算得到合理的计算 Qa(a=oaai 结果，并且能够捕捉到不同流动方式之间的显著差 异 对于多流中间包，Q:为每一流对应的流量，为 方差。常用于计算一组数据的离散度，也就可 便于讨论，本文假设各流流量均相等（如果各流流量 以衡量停留时间分布曲线上数据的离散程度，结合 不相等则代入实际流量)，示踪剂注入浓度c0等于 中间包RTD曲线，方差表示如下： 1,示踪剂注入时间△t等于1s·归一化后可以得到 般的量纲1形式，离散后得： =J0(4一)'E()d= 0(9-1(9a0=6所E(9a0-1, 22(9a=1. g=6/x, =2 △0=△t/t, =1 该方差可以描述实际流动模式与理想的活塞流 E(91=.00s-a(9p) Qco△t 动之间的偏差.当流动模式仅仅呈现出理想活塞流 t=V/Q: 时，方差取最小值0：当流动呈现完全混合流时（完 全混合流动是最大程度的返混流动，返混区域在化 [6a 学反应工程学中定义为：包含不同停留时间流体的 此处，系数α应该引起关注，它是由测量误差 区域)，方差将达到最大值1；当流动模式中存在短 或计算误差引起的，通常在实验过程中，可以检测 路流区、死区或部分返混区时，方差的取值将在0和 某一点的浓度，但是基于质量守恒方程计算流出的 1之间，返混区越大，RTD曲线越平展，方差取值越 示踪剂总量时，需要出口平面上的平均浓度.因此， 大 根据检测到的点浓度，需要通过α系数等效成出口 当前并没有明确的针对多流中间包流动特征的 平面上的平均浓度， 数学模型，提出新模型需要考察多流中间包内活塞 如图7所示，如果E()在0达到某一个趋向 区、返混区、死区和短路流的特性，并规范停留时间 于0的数时，其值不为0，且达到一个较高的值后 分布(RTD)的相关计算，本文尝试提出多流中间包 突然下降，这表明存在短路流： 流动特征模型以供讨论，模型示意图如图6所示， -5以1af-0, 示踪剂在中间包入口的输入方式为近似δ函
θav＝ （ V — V d）／Q′ V／Q Q′＝ Q— Qby． 模型中也不考虑短路流则： Qby＝0Q′＝ Q． 所以在文献［7］中得到死区的体积分数为： 1—θav＝1— （ V — V d）／Q V／Q ＝ V d V ． 但是在实际冶金过程中死区的定义是一个相对的 概念：超过平均停留时间2倍以上的区域并不是绝 对静止的．因而中间包实际有效容积并不能等效 为 V — V d文献［7］中计算死区的体积分数也就不 合理． 综上所述水模型实验和数值模拟结果均显示 采用 Sahai 的针对单流的数学模型来计算多流中间 包内各区域体积分数存在不足因而利用其结果指 导多流中间包优化设计和多流中间包冶金是不可靠 的． 3 提出新模型 一个合理的计算模型应该计算得到合理的计算 结果并且能够捕捉到不同流动方式之间的显著差 异． 方差 σ2 常用于计算一组数据的离散度也就可 以衡量停留时间分布曲线上数据的离散程度结合 中间包 RTD 曲线方差表示如下： σ2 i＝∫ ∞ 0 （ tj—tj） 2E（ tj）d t＝ ∫ ∞ 0 （θj—1） 2E（θj）dθ＝∫ ∞ 0 θ2 jE（θj）dθ—1 σ2＝ ∑ n i＝1 σ2 i． 该方差可以描述实际流动模式与理想的活塞流 动之间的偏差．当流动模式仅仅呈现出理想活塞流 时方差取最小值0；当流动呈现完全混合流时（完 全混合流动是最大程度的返混流动返混区域在化 学反应工程学中定义为：包含不同停留时间流体的 区域）方差将达到最大值1；当流动模式中存在短 路流区、死区或部分返混区时方差的取值将在0和 1之间返混区越大RTD 曲线越平展方差取值越 大． 当前并没有明确的针对多流中间包流动特征的 数学模型提出新模型需要考察多流中间包内活塞 区、返混区、死区和短路流的特性并规范停留时间 分布（RTD）的相关计算．本文尝试提出多流中间包 流动特征模型以供讨论模型示意图如图6所示． 示踪剂在中间包入口的输入方式为近似 δ函 数输入其输入总量为： M＝ Qc0Δt′． 式中M 为示踪剂在入口的输入总量；c0 为示踪剂 输入浓度；Δt′为示踪剂注入时间为满足 δ函数 相对整个计算时间该值应该非常小；Q 为流体在入 口的体积流量． 图6 本文模型示意图 Fig．6 Schematic diagram of the proposed model 示踪剂在中间包出口的输出总量为： M＝∫ n∫1 ∞ 0 Qici（ tj）d t． 根据质量守恒总的输出量等于总的输入量： ∫ n∫1 ∞ 0 Qici（ tj）d t＝ Qc0Δt′ 对于多流中间包Qi 为每一流对应的流量．为 便于讨论本文假设各流流量均相等（如果各流流量 不相等则代入实际流量）示踪剂注入浓度 c0 等于 1示踪剂注入时间Δt′等于1s．归一化后可以得到 一般的量纲1形式离散后得： ∑ n i＝1 ∑ ∞ j＝0 ［ Ei（θj）Δθ］＝1 θj＝tj／τ Δθ＝Δt／τ Ei（θj）＝α （1／n） Qci（ tj）τ Qc0Δt′ ＝α 1 n τci（θτj ） τ＝ V／Q α＝c0Δt′ ∑ n i＝1 ∑ ∞ j＝0 1 n ci（ tj）Δt ． 此处系数 α应该引起关注它是由测量误差 或计算误差引起的．通常在实验过程中可以检测 某一点的浓度但是基于质量守恒方程计算流出的 示踪剂总量时需要出口平面上的平均浓度．因此 根据检测到的点浓度需要通过 α系数等效成出口 平面上的平均浓度． 如图7所示如果 E（θ）在 θ达到某一个趋向 于0的数θ′时其值不为0且达到一个较高的值后 突然下降这表明存在短路流： V byi V ＝∫ θ′ 0 Ei（θ）dθθ′→0 ·818· 北 京 科 技 大 学 学 报 第31卷
第7期 潘宏伟等：多流中间包流动特征的数学模型 ,819 3.5 3.0 =-, 2.5 RTD曲线初始阶段，存在 a=之yad 2.0 一个迅速的急剧上升和一个 =1 15 并发的急剧下降过程，表明 根据质量守恒，可以得到活塞区的体积分数： 存在短路流 1.0 Vpi1 Vhai Vdi Vbyi VR VV 0.5 =2片 基于新的多流中间包流动特征的数学模型，根 图7E曲线上短路流示意图 据原有的数值模拟得到的数据，重新计算如图4所 Fig.7 Schematic diagram of the character of bypass flow in the E- 示中间包内各区域体积分数，如表3所示.计算结 curve 果显示：①1和2*中间包内不存在短路流；②新的 模型计算的死区体积分数为正值且数据比较合理， V 各区域体积分数之和为100%：③捕捉到了1#中间 0可以在I曲线上进行定义、I曲线描述的是， 包内强烈的返混特点和2中间包内强烈的活塞流 从0到0时间段内留在容器内的示踪剂分布函数， 流动趋势 可以表示如下： 为了验证数学模型的可靠性，对1#中间包和 I(0=1- E(0d0. 2#中间包分别进行了水模型实验·基于中间包流场 的对称性，水模型实验测得了中间包右侧第3流和 如图8所示，在0达到0以前有个急剧的初始 第4流出口电压随时间的变化曲线.由于出口电压 下降阶段，在【曲线上，急剧下降的线段长度就是 值是浓度的线性函数，为了便于利用本文模型，将出 短路流的体积分数 口电压值近似视为出口浓度进行数据处理，其误差 1.0 可利用本文模型中的修正系数α进行调整，处理结 /曲线上，急剧下降的 果如图9所示，这样既能保证水模型数据的真实 线段长度即是短路流 0.6 体积分数 性，又便于利用本文模型进行计算，中间包各区域的 体积分数计算结果如表3所示. 04 0.14 0.2 0.12 2”中间包 0.10 0.08 图8I曲线上短路流示意图 0.06 Fig.8 Schematie diagram of the character of bypass flow in the I- 0.04 1“中间包 curve 0.02 另外，根据死区的定义（流体的停留时间超过平 0.5 1.0 15 2 均停留时间2倍以上的区域)可得： -2a0, 图91云和2“中间包水模型实验结果对比 Fig.9 Water model experimental results for Tundishes No.I and va=vai No-2 V 台v 表3中的数据表明，对于2中间包，水模型实 此处，积分下限可以根据实际需求进行调整 验的计算结果与数值模拟的计算结果吻合较好，证 分析证实，死区和返混区使得实际流动形式在 实了本文模型的可靠性.但是，对于1中间包，水 偏离理想活塞流的过程中具有共性，通过上述计算 模型实验结果与数值模拟结果间有一定差异.其原 得到的死区体积分数，从而可以得到返混区的体积 因主要是：本文提出的数学模型基于质量守恒，要求 分数： 中间包各流RTD曲线下的积分面积之和为1；总体
图7 E 曲线上短路流示意图 Fig．7 Schematic diagram of the character of bypass flow in the Ecurve V by V ＝ ∑ n i＝1 V byi V ． θ′可以在 I 曲线上进行定义．I 曲线描述的是 从0到θ时间段内留在容器内的示踪剂分布函数 可以表示如下： I（θ）＝1—∫ θ 0 E（θ）dθ． 如图8所示在θ达到θ′以前有个急剧的初始 下降阶段在 I 曲线上急剧下降的线段长度就是 短路流的体积分数． 图8 I 曲线上短路流示意图 Fig．8 Schematic diagram of the character of bypass flow in the Icurve 另外根据死区的定义（流体的停留时间超过平 均停留时间2倍以上的区域）可得： V di V ＝∫ ∞ 2 Ei（θ）dθ V d V ＝ ∑ n i＝1 V di V ． 此处积分下限可以根据实际需求进行调整． 分析证实死区和返混区使得实际流动形式在 偏离理想活塞流的过程中具有共性．通过上述计算 得到的死区体积分数从而可以得到返混区的体积 分数： V bai V ＝σ2 i— V di V — V byi V V ba V ＝ ∑ n i＝1 V bai V ． 根据质量守恒可以得到活塞区的体积分数： V pi V ＝ 1 n — V bai V — V di V — V byi V V p V ＝ ∑ n i＝1 V pi V ． 基于新的多流中间包流动特征的数学模型根 据原有的数值模拟得到的数据重新计算如图4所 示中间包内各区域体积分数如表3所示．计算结 果显示：①1＃和2＃中间包内不存在短路流；②新的 模型计算的死区体积分数为正值且数据比较合理 各区域体积分数之和为100％；③捕捉到了1＃中间 包内强烈的返混特点和2＃中间包内强烈的活塞流 流动趋势． 为了验证数学模型的可靠性对1＃ 中间包和 2＃中间包分别进行了水模型实验．基于中间包流场 的对称性水模型实验测得了中间包右侧第3流和 第4流出口电压随时间的变化曲线．由于出口电压 值是浓度的线性函数为了便于利用本文模型将出 口电压值近似视为出口浓度进行数据处理其误差 可利用本文模型中的修正系数 α进行调整处理结 果如图9所示．这样既能保证水模型数据的真实 性又便于利用本文模型进行计算中间包各区域的 体积分数计算结果如表3所示． 图9 1＃和2＃中间包水模型实验结果对比 Fig．9 Water-model experimental results for Tundishes No．1 and No．2 表3中的数据表明对于2＃中间包水模型实 验的计算结果与数值模拟的计算结果吻合较好证 实了本文模型的可靠性．但是对于1＃ 中间包水 模型实验结果与数值模拟结果间有一定差异．其原 因主要是：本文提出的数学模型基于质量守恒要求 中间包各流 RTD 曲线下的积分面积之和为1；总体 第7期 潘宏伟等： 多流中间包流动特征的数学模型 ·819·