式中:P为流体密度，k为湍动能。为时间，二为速度矢量，二为空间坐标声为动力茹度，二为湍动能耗散项，从为湍动茹度，G、是平均速度梯度引起的湍动能k的产生项G。是浮力引起的湍动能k的产生项，瑞是可压湍流中脉动扩张对湍动能k的贡献，C,e ,矶二、几二为经验常数，o-、和o-二分别为湍动能k和湍动耗散率二对应的Prandtl数，凡和S、是用户定义的源项。式中为位置向量，为方向向量，为散射方向，a为吸收系数，n为折射系数，o-为散射系数，I(r,s)为位置:处沿方向‘的辐射强度，o-为Stefan-Boltzmann常数(5.67x10-8W·m一zm·K一“),T为温度，动为相位函数，仔为空间立体角，s为介质的光学深度。式中:E为玻璃的辐射系数，QCM火焰空间向玻璃硅碳棒窑炉的传热量，双，为玻璃液面温度。计算火焰空间对配合料表面的传热时，将玻璃液的辐射系数、温度相应替换为配合料的辐射系数和温度。物料性质参数及边界条件硅碳棒窑炉熔制玻璃为E玻璃纤维，建模所需的玻璃液性质参数见表3所示。玻璃纤维硅碳棒窑炉模型中边界条件设置见表4所示图3为模拟得到顶烧玻璃硅碳棒窑炉底部温度与硅碳棒窑炉底部热电偶测量温度对比图。其中池底温度误差最大值位于距离硅碳棒窑炉前端池壁位置7. 5 m处，模拟温度为1486℃，实际测得温度1 437℃，模拟温度较实际测量值偏高，主要是由于模拟得到的池底温度受玻璃液有效导热系数以及池底散热情况的影响，建模采用的玻璃液有效导热测量困难，只能采用经验值，且池底耐火材料散热复杂，模型中散热边界条件与实际硅碳棒窑炉散热有差异。图3表明模拟硅碳棒窑炉温度与实际测得温度变化趋势相符，模拟结果精度较高，能准确模拟硅碳棒窑炉内部璃液温度场。http://www.zbqunqiang.cn/