浮法玻璃熔窑液流的三维数学模型.docxVIP

浮法玻璃熔窑液流的三维数学模型

1玻璃熔窑的研究意义

浮法玻璃生产工艺是目前平板玻璃生产工艺中最先进的生产工艺。熔窑则是浮法生产工艺中最重要的热工设备之一。因此,加强对玻璃熔窑的科学研究,深入了解掌握窑内液流的温度场、速度场分布,提高玻璃熔窑设计水平,改进熔窑结构和作业条件,对降低玻璃生产成本,提高玻璃质量,节约能源都具有重要意义。

本文运用SIMPLE算法,采用三维数学模型计算了卡脖开度变化时的玻璃液温度场和速度场,并对其产生的影响作了分析和研究。

2玻璃熔制过程数学模型

本文采用一个能模拟浮法熔窑内的主要物理输运过程,反映玻璃生产中的一些主要影响因素的三维数学模型,该模型主要由配合料熔化模型、与上部火焰空间耦合的半经验模型和玻璃液运动模型组成,从而为工程界理解和优化玻璃生产工艺提供一个指导工具,为将来可以发展成为玻璃熔制过程综合模型奠定基础。

2.1熔融玻璃液的特性

在玻璃液流运动模型中,作如下假设:

(1)作为一个稳定运行的窑炉,其平均运动状态是相对稳定的,处于稳态;

(2)熔融玻璃液充满整个熔化池,配合料的作用只表现在液面温度分布中;

(3)熔融玻璃液是均质的牛顿粘性流体;

(4)熔融玻璃液处于蠕动状态,可以忽略摩擦升温;

(5)忽略气泡和化学反应对流动和温度分布的影响;

(6)玻璃液定压比热Cp和膨胀系数β不随温度变化;

(7)玻璃液密度、粘度、导热系数均为温度T的函数。

2.2显色剂为理想观念的通用方程

作了上述简化之后,玻璃液在窑池内状态可用流体连续性方程,运动方程,动量守恒方程及能量守恒方程来描述,对这些方程联立求解,即可获得熔窑中玻璃液流的速度场和温度场。为了便于编制通用计算机程序,将上述方程写成通用微分方程(1),式中各符号的含义见表1。

??τ(?)+div(→U?)=div(Γ?grad?)+S?(1)??τ(?)+div(U??)=div(Γ?grad?)+S?(1)

该方程中的四项依次分别是不稳态项(由于假设玻璃液是稳态流动,故略去此项)、对流项、扩散项以及源项。其中ρ表示密度,τ表示时间,Γ?表示扩散系数,S?表示源项,?表示因变量。对于特定意义的?,具有特定Γ?的量和S?。

2.3边境条件

2.3.1热边界条件

(1)辐射传热系数的计算

在玻璃熔窑表面配合料单位面积内得到的净热量可由下式表示:

QΓ=Cfm?[(Τf100)4-(Τm100)4]+αfm?(Τf-Τm)(2)QΓ=Cfm?[(Tf100)4?(Tm100)4]+αfm?(Tf?Tm)(2)

式中,Cfm-火焰与物料之间的辐射传热系数;Tf-火焰温度/K;Tm-物料温度/K;αfm-火焰与物料之间的对流换热系数/W·m-2·K-1。

(2)温度b土壤散热系数t

在卡脖和冷却部,玻璃液的热量主要是通过窑墙、大碹和胸墙内壁传至外壁,再通过自然对流和辐射散热散到大气环境中去,其热流方程由下式表示:

Q=ε0?C0?[(Τg100)4-(Τin100)4]?F1=Κw?(Τin-Τout)?F2/δ=α?(Τout-Τair)?F3=α′?(Τg-Τair)?F1(3)ε0=11ε1+F1F2?(1ε1-1)(4)

式中,Q-卡脖至出料口液面总散热量/W;ε0-玻璃液面与窑墙、大碹与胸墙内壁之间的传热黑度;ε1-玻璃液面黑度;C0-黑体辐射系数,等于5.69W·m-2·K-1;Tg-卡脖至出料口液面平均温度/K;Tin-窑墙内壁平均温度/K;Tout-窑墙外壁平均温度/K;Tair-环境温度/K;Kw-窑墙平均导热系数/W·m-2·K-1;F1-卡脖至出料口液面面积/m2;F2-窑墙内壁表面积/m2;F3-窑墙外壁表面积/m2;δ-窑墙厚度/m;α-窑墙外壁平均换热系数/W·m-2·K-1;α′-等效换热系数/W·m-2·K-1。

(3)导热系数的测定

池底和池壁热流强度可表述为:

Qw=Κw?(Τin-Τout)/δ=α?(Τout-Τair)=α′?(Τin-Τair)(5)

式中,Kw-池底或池壁窑墙的导热系数/W·m-2·K-1;α-池壁外侧换热系数/W·m-2·K-1;α′-等效换热系数/W·m-2·K-1,其余符号解释同上。

2.3.2速度边界条件

(1)配合料前进速度

在升温段:U=U0,V=0,W=0(6)

其中U0为配合料前进速度;

在熔化段:U=U0,V=0,W=Wbij(7)

Wbij为第i,j单元熔化的配合料进入玻璃液速度。

(2)玻璃溶液的自由表面

?U?z=0,?V?z=0,W=0(8)

(3)进料口

U=Uout(9)

Uout—出料口平均速度,其大小随出料量变化。

(4)池壁界面

U=0?V=0?W=0(10)

3卡充填热性能计算

针对杭州新型建