鉑金通道溫度控制工藝的數(shù)值模擬?

李 青 汪偉軍 王麗紅 姜文成 郝 藝

（1.東旭集團(tuán)有限公司 石家莊 050021）（2.平板顯示玻璃技術(shù)和裝備國家工程實驗室 石家莊 050035）

鉑金通道溫度控制工藝的數(shù)值模擬?

李 青1，2汪偉軍1，2王麗紅1，2姜文成1，2郝 藝1，2

（1.東旭集團(tuán)有限公司 石家莊 050021）（2.平板顯示玻璃技術(shù)和裝備國家工程實驗室 石家莊 050035）

應(yīng)用流體力學(xué)專業(yè)的數(shù)值分析軟件Ansys-Fluent，對鉑金通道進(jìn)行溫度場與流場的耦合計算；應(yīng)用軟件的“用戶定義函數(shù)（udf）”，編寫自動調(diào)節(jié)鉑金通道的加熱功率的程序，使各段鉑金通道的溫度達(dá)到預(yù)定目標(biāo)。從計算結(jié)果分析得到結(jié)論：在溫度上升段鉑金通道內(nèi)，對壁面附近的氣泡上升有利，但當(dāng)氣泡上升到中間或接近頂部位置時，上升速度會降低甚至停止；在溫度下降段鉑金通道內(nèi)，對壁面附近的氣泡上升不利，但對已經(jīng)上升到中間或接近頂部位置的氣泡溢出影響不大；在工藝設(shè)計時應(yīng)將溫度上升段管道放在前面，下降段放在后面。通過計算加熱功率與實際使用加熱功率比較，兩者基本一致，驗證了加熱功率數(shù)值模型的可靠性。

液晶基板玻璃；數(shù)值模擬；鉑金通道；溫度控制；加熱功率

1 引言

在液晶基板玻璃生產(chǎn)中，從窯爐出來的熔融玻璃里面含有大量的氣泡，這些氣泡如果被帶到最終的產(chǎn)品中，可能形成缺陷。光學(xué)玻璃、汽車玻璃以及IT行業(yè)使用的玻璃對氣泡的要求非常嚴(yán)格，直徑大于0.2mm的氣泡為缺陷氣泡，要求每公斤玻璃中直徑小于 0.2mm的氣泡不得大于 0.5個［1～2］。熔融玻璃的澄清溫度高、過程復(fù)雜，溢流下拉法制備基板玻璃的澄清過程主要發(fā)生在鉑金通道中。鉑金通道通過分段調(diào)節(jié)控制溫度使熔融玻璃中的氣泡減少或排除。實際生產(chǎn)中一般通過調(diào)節(jié)鉑金管道的加熱功率來控制鉑金通道的溫度。

近年來，隨著液晶基板玻璃在國內(nèi)的發(fā)展，國內(nèi)對鉑金通道的研究也逐步增多。科技期刊公開發(fā)表的鉑金通道相關(guān)文獻(xiàn)數(shù)量較少，2012年王庚、李桂玲、張濤等［3］對鉑金通道的熱電偶的波動進(jìn)行了研究并提供了對策。2015年吳愛軍、張濤、程秦川、王群偉等［4］研究了鉑金通道耐火材料的問題，找出了符合鉑金通道工藝要求的耐火材料組分比例。2016年，李青、張春喜等［5］探討了延長鉑金通道壽命的問題。

鉑金通道研究多集中在專利方面，2009年，宋金虎、鄭權(quán)等［6］發(fā)明了一種鉑金通道電加熱法蘭結(jié)構(gòu)；2010年，王保平、劉文泰、李兆廷等［7］發(fā)明的一種鉑金通道溫度調(diào)節(jié)裝置；2010年，李兆廷、王保平、斯沿陽、楊亞彬等［8］發(fā)明的一種鉑金通道電加熱結(jié)構(gòu)。以上發(fā)明內(nèi)容大都集中在加熱裝置或結(jié)構(gòu)方面。在鉑金通道加熱工藝方面的發(fā)明有如：韓俊、宋金虎等［9］在2007年申請了專利，該專利發(fā)明了一種鉑金通道電加熱工藝設(shè)計的方法，該方法通過計算鉑金管壁以及法蘭的電阻率等來計算加熱功率，但是沒有考慮通道的熱量散失；2012年，程英、樊宏剛［10］發(fā)明的一種基于相序設(shè)置的鉑金通道加熱方法，此方法重點電路方面的設(shè)計，實現(xiàn)對法蘭及銅排溫度的控制。歸納國內(nèi)的文獻(xiàn)和專利等，多在研究鉑金通道的結(jié)構(gòu)或者是加熱裝置，僅有的韓俊、宋金虎等申請的專利是研究鉑金通道加熱功率的計算方法，可以為工藝設(shè)計作指導(dǎo)，但美中不足的是此方法只是考慮鉑金管壁以及法蘭的電阻率，沒有考慮通道的熱量損失。

隨著計算機的快速發(fā)展，特別是近年來CAE（Computer Aided Engineer，CAE）技術(shù)的廣泛應(yīng)用，為鉑金通道的加熱工藝設(shè)計提供了一種切實有效的手段。本文利用流體力學(xué)專業(yè)分析軟件Ansys-Fluent的流動-熱耦合模塊，實現(xiàn)對熔融玻璃流場和熱場的耦合模擬計算。應(yīng)用“用戶定義函數(shù)（udf）”編寫一個簡單的程序，使軟件根據(jù)出口計算所得溫度與目標(biāo)控制溫度比較，自動調(diào)節(jié)加熱功率大小。

本文通過對鉑金通道的模擬計算，分析鉑金通道不同的溫度段對氣泡澄清的作用，以解決鉑金通道的溫度工藝布置問題。同時，將模擬計算所得電加熱功率與實際生產(chǎn)電加熱功率進(jìn)行比較，驗證計算模型的可靠性，為電加熱工藝參數(shù)提供參考。

2 幾何尺寸與材料參數(shù)及邊界條件

根據(jù)熔融玻璃澄清工藝建立的溫度制度，為使熔融玻璃在短時間、短距離內(nèi)實現(xiàn)高效澄清，設(shè)計鉑金通道的結(jié)構(gòu)如圖1所示，鉑金通道內(nèi)壁為鉑金管壁，熔融玻璃在鉑金管壁內(nèi)流動，鉑金管壁外面包覆著五層耐火材料，這些材料的特性均不相同，耐火材料外為生產(chǎn)環(huán)境。實際生產(chǎn)中，需要通過對鉑金管壁進(jìn)行電加熱來控制熔融玻璃的溫度。在工藝設(shè)計階段，由于鉑金通道結(jié)構(gòu)的復(fù)雜性為加熱功率工藝設(shè)計帶來了極大的難度。

圖1 通道結(jié)構(gòu)三維圖

2.1 鉑金通道各段尺寸和材料參數(shù)

鉑金通道按照不同的功能可分為四段（鉑金通道1、2、3、4），每一段鉑金通道包括鉑金壁管和外部的五層耐火材料，它們的幾何尺寸如表1。

表1 鉑金通道幾何尺寸

鉑金通道包括鉑金壁管和外部的五層耐火材 料，它們的材料參數(shù)如表2。

表2 耐火材料參數(shù)

2.2 熔融玻璃的參數(shù)

各段鉑金通道中熔融玻璃的參數(shù)如表3。

表3 熔融玻璃參數(shù)

2.3 鉑金通道各段的邊界條件

鉑金通道各段的邊界條件如表4。

表4 通道各段邊界條件

3 建模、網(wǎng)格劃分及計算模型設(shè)定

3.1 建模和網(wǎng)格劃分

如圖1，利用三維CAD軟件，根據(jù)鉑金通道幾何尺寸數(shù)據(jù)生成三維實體模型。計算時利用對稱條件以減小計算量，所以只生成二分之一模型。

利用Fluent前置處理軟件gambit對模型進(jìn)行網(wǎng)格劃分，網(wǎng)格全部采用六面體單元，如圖2，網(wǎng)格數(shù)1098762。

圖2 鉑金通道網(wǎng)格圖

3.2 計算模型設(shè)定

生產(chǎn)過程中鉑金通道的傳熱現(xiàn)象錯綜復(fù)雜，主要包括以下六個方面的熱交換：熔融玻璃內(nèi)部之間的熱交換、熔融玻璃與鉑金管壁之間的熱交換、鉑金管壁的電加熱熱源、鉑金管壁與耐火材料之間的熱交換、各層耐火材料之間的熱交換、外圍耐火材料與環(huán)境之間的熱交換，在正常生產(chǎn)中。這是一個穩(wěn)定的熱交換過程，計算時也采用穩(wěn)態(tài)計算。在建立數(shù)學(xué)模型時，做如下簡化：

1）一般中堿玻璃在1200℃的溫度時電阻率約為1000Ω·m；而無堿玻璃約為10000～50000Ω·m［11］，而鉑金在1200℃的溫度時電阻率約為0.00465Ω·m［12］，它們的電阻率相差好幾個量級，因此認(rèn)為電加熱功率都作用在鉑金管壁上，忽略熔融玻璃的電加熱效應(yīng)。

2）鉑金材料和耐火材料之間以及各層耐火材料之間，存在接觸熱阻，但是在安裝工藝中已經(jīng)在這些縫隙之間灌注了耐火泥，因此認(rèn)為各個材料之間沒有熱阻。

3）認(rèn)為熔融玻璃、鉑金材料、各個耐火材料的熱學(xué)性能為各向同性。

4）認(rèn)為在同一段鉑金通道中熔融玻璃、鉑金材料、各個耐火材料的熱學(xué)性能不隨溫度發(fā)生變化。

5）在鉑金通道最外層，它通過對流轉(zhuǎn)熱和輻射傳熱，向周圍環(huán)境散失熱量。根據(jù)斯特藩-玻爾茲曼定律［13］，耐火材料的輻射力為

其中：ε為法向發(fā)射率，耐火磚大約為0.5左右；σ為斯特藩-玻爾茲曼常數(shù)，值為5.67W/m^2；T1為外層耐火材料溫度，t2為廠房墻壁內(nèi)面溫度。

正常生產(chǎn)中，由于鉑金通道完全處于廠房內(nèi)，基本上是一個封閉的系統(tǒng)，當(dāng)生產(chǎn)達(dá)到穩(wěn)定狀態(tài)時，t1與t2可以認(rèn)為是相等的，即是說耐火材料的輻射散熱可以忽略不計。在本文中只計算它的對流散熱。

在設(shè)定計算模型時，采用壓力速度耦合的simple算法來計算控制方程的傳輸項；動量方程、能量方程采用二階離散迎風(fēng)格式；湍流模型采用標(biāo)準(zhǔn)的k-ε模型［14］。

用“用戶定義函數(shù)”（簡稱：UDF）編制程序（程序原代碼見附錄）來使計算程序自動調(diào)節(jié)鉑金管壁電加熱功率。此程序?qū)⒃谟嬎闶諗亢螅x取并計算出口的平均溫度，將此溫度與出口欲控制的目標(biāo)溫度進(jìn)行比較，再根據(jù)比較結(jié)果調(diào)整鉑金管壁的加熱功率源項，直至出口溫度達(dá)到所要求的溫度為止。

4 計算結(jié)果分析

4.1 出口溫度分布

如圖3是出口橫截面上的溫度分布圖，其中a、b、c、d分別對應(yīng)鉑金通道1、鉑金通道2、鉑金通道3和鉑金通道4。

圖3 鉑金通道出口溫度分布圖

由計算結(jié)果得到，圖3中a、b、c接近鉑金管壁的熔融玻璃溫度較中間的溫度高，而d則與之相反。由設(shè)定的條件知道，鉑金通道1、2、3的出口溫度高于進(jìn)口溫度，在熔融玻璃流經(jīng)這些管道時，鉑金管壁在對熔融玻璃逐漸加溫，因此接近管壁的熔融玻璃溫度首先升高，使之高于中間熔融玻璃的溫度；而鉑金通道4的情況恰好與之相反。在鉑金通道1、2、3中，由于接近管壁的熔融玻璃溫度升高，在管壁附近的氣泡受浮力將增大，同時熔融玻璃的粘度將減小，氣泡所受的黏滯力將減小，這些都有助于氣泡的上升；但是當(dāng)氣泡到達(dá)中間及頂部位置后，由于這里的溫度降低，氣泡上升速度會降低甚至停在。而在鉑金通道4中，由于管壁附近的溫度低、中間及頂部溫度高，而整個管道內(nèi)從進(jìn)口到出口縱向的溫度在逐漸降低，特別是管壁附近的溫度降低最為明顯，所以，這段映鉑金通道中不利于管壁附近的氣泡上升，對于已經(jīng)到達(dá)中間位置的氣泡溢出影響不大。所以，有利于氣泡溢出的角度出發(fā)，應(yīng)該將溫度上升段的映鉑金通道布置在前面，而溫度降低段的映鉑金通道應(yīng)該布置在其后。

4.2 計算結(jié)果與討論

根據(jù)“用戶定義函數(shù)”自編程序計算結(jié)果數(shù)據(jù)如表5，其中“單位體積熱功率”是指鉑金壁面單位體積的加熱功率，“Pt體積”是指鉑金管壁的體積，“熱功率”是指鉑金管壁加熱總功率。

由表5所得計算加熱功率可以看出鉑金通道4所需功率明顯小于鉑金通道1、2、3，這是由于鉑金通道1、2、3所控制的出口溫度大于進(jìn)口溫度，所加載的熱功率除了一部分用來提高熔融玻璃溫度外，還有一部分用來補償鉑金管道的熱量散失。而鉑金通道4中，所控制的出口溫度小于進(jìn)口溫度，所加載的熱功率只需補償一部分熱散失即可。

表5 計算結(jié)果

如表6所列，為各鉑金通道計算所需加熱功率與實際所加載的功率比較。

表6 功率比較

由表6各個鉑金通道比較可以看出，實際使用熱功率與計算所得熱功率非常接近。鉑金通道1、2、3實際使用熱功率比計算所得熱功率略大，是因為計算時，熱輻射以及溢出氣泡帶走的熱量沒有計算；而鉑金通道4實際實際使用熱功率比計算所得熱功率略小，原因同上所述。

5 結(jié)語

應(yīng)用CAE技術(shù)，對鉑金通道的熱場和流場耦合分析，根據(jù)分析可以得到如下結(jié)論：

1）鉑金通道出口溫度大于進(jìn)口溫度的管道內(nèi)，管壁附近的熔融玻璃溫度高于中間及頂部的熔融玻璃溫度，有利于管道壁面附近的氣泡上升，但是氣泡上升到中間位置后，上升速度會逐漸降低甚至停止；在鉑金通道出口溫度小于進(jìn)口溫度的管道內(nèi)，管壁附近的熔融玻璃溫度低于中間及頂部的熔融玻璃溫度，由于管道內(nèi)壁面附近縱向的溫度在逐漸降低，所以這段管道不利于壁面附近的氣泡上升，而對于已經(jīng)到達(dá)中間或接近頂部的氣泡的溢出影響不大。在設(shè)計管道溫度工藝時，應(yīng)將溫度上升的管道設(shè)計在前，而將溫度下降的管道設(shè)計在后。

2）計算所得加熱功率與實際使用基本一致，可以應(yīng)用CAE技術(shù)為加熱功率工藝參數(shù)設(shè)計作指導(dǎo)。在實際生產(chǎn)中，溫度上升段管道所用功率應(yīng)比計算所得功率略大；在溫度下降段管道，反之。

［1］董振波.TFT-LCD玻璃基板生產(chǎn)中氣泡缺陷分析與對策［C］//中國硅酸鹽學(xué)會電子玻璃分會：2014年光電子玻璃技術(shù)研討會，2014.131-131.DONG Zhenbo.Analysis and Countermeasure of Bubble Defect in Production of TFT-LCD Glass Substrate［C］//China Silicate Society Electronic Glass Branch：2014 Optoelectronic Glass Technology Seminar，2014.131-131.

［2］王俊明，董振波，付冬偉，等.電子玻璃基板生產(chǎn)中氣泡缺陷分析與對策［J］.電子技術(shù)與軟件工程，2013（16）：160-160.WANG Junming，DONG Zhenbo，F(xiàn)U Dongwei，et al.Analysis and Countermeasure of Bubble Defect in Production of Electronic Glass Substrate［J］.Electronic Technology and Software Engineering，2013（16）：160-160.

［3］.王庚，李桂玲，張濤.TFT基板玻璃通道鉑銠熱電偶溫度波動現(xiàn)象分析和對策［J］.電子世界，2012（13）：79-80.WANG Geng，LI Guiling，ZHANG Tao.Analysis and CountermeasureofTemperatureFluctuation ofPlatinum-rhodium Thermocouple on Glass Substrate Glass Channel［J］.Electronic World，2012（13）：79-80.

［4］吳愛軍，張濤，程秦川，等.TFT-LCD玻璃基板熔爐用鋯質(zhì)填充料的開發(fā)及應(yīng)用［C］//2015耐火材料綜合學(xué)術(shù)年會暨全國不定形耐火材料學(xué)術(shù)會議：2015耐火原料學(xué)術(shù)交流會，2015.56-57.WU Aijun，ZHANG Tao，CHENG Qinchuan，et al.Development and Application of Zirconium Filler for TFT-LCD Glass Substrate［C］//2015 Conference on Refractory Materials and National Conference on Unshaped Refractories：2015 Insulation Seminar for Refractory Raw Materials，2015.56-57.

［5］李青，張春喜，斯沿陽，等.延長鉑金通道使用壽命的措施探討［J］.電子玻璃技術(shù)，2016（48）：93-95.LI Qing，ZHANG Chunxi，SI Yanyang，et al.Discussion on the Measures to Extend the Service Life of Platinum Channel［J］.Electronic Glass Technology，2016（48）：93-95.

［6］宋金虎，鄭權(quán)，劉文泰，等.一種對鉑金通道電加熱的法蘭結(jié)構(gòu)：中國，CN101708946A［P］，2009.11.02.SONG Jinhu，ZHENG Quan，LIU Wentai，et al.A flange structure for platinum channel heating： China，CN101708946A［P］，2009.11.02.

［7］王保平，劉文泰，李兆廷.一種鉑金通道溫度調(diào)節(jié)裝置：中國，CN201660542U［P］，2010.4.9.WANG Baoping，LIU Wentai，LI Zhaoting.A platinum channel temperature regulator：China，CN201660542U［P］，2010.4.9.

［8］李兆廷，王保平，斯沿陽，等.一種鉑金通道電加熱結(jié)構(gòu)：中國，CN201660544U［P］，2010.4.22.LI Zhaoting，WANG Baoping，SI Yanyang，et al.A platinum channel electric heating structure： China，CN201660544U［P］，2010.4.22.

［9］韓俊，宋金虎.鉑金通道電加熱設(shè)計方法：中國，CN101144925［P］，2007.9.30.HAN Jun，SONG Jinhu.Design method of platinum channel electric heating：China，CN101144925［P］，2007.9.30.

［10］程英，樊宏剛.一種基于相序設(shè)置的鉑金通道加熱方法：中國，CN103011555A［P］，2012.12.20.CHENG Ying，F(xiàn)AN Honggang.A platinum channel heating method based on phase sequence setting：CN103011555A［P］，2012.12.20.

［11］佚名.熔融玻璃電阻率的測定［J］.玻璃纖維，1977（6）：11-13.YI Ming.Determination of resistivity of molten glass［J］.Glass fiber，1977（6）：11-13.

［12］Callister W D，Rethwisch D G.Fundamentals of Materials Science and Engineering［M］.北京：化學(xué)工業(yè)出版社，2004：66-73.Callister W D，Rethwisch D G.Fundamentals of Materials Science and Engineering［M］.Beijing：Chemical Industry Press，2004：66-73.

［13］楊世銘，陶文銓.傳熱學(xué)［M］.北京：高等教育出版社，1998：272-284.YANG Shiming，TAO Wenquan.Heat transfer［M］.Beijing：Higher Education Press，1998：272-284.

［14］任志安，郝點，謝紅杰.幾種湍流模型及其在Flute中的應(yīng)用［J］.化工裝備技術(shù)，2009，30（2）：38.REN Zhian，HAO Dian，XIE Hongjie.There are several turbulence models and their application in Flute［J］.Chemical Equipment Technology，2009，30（2）：38.

Numerical Simulation of Platinum Channel Temperature Control Process

LI Qing1，2WANG Weijun1，2WANG Lihong1，2JIANG Wencheng1，2HAO Yi1，2
（1.Tunghsu Group Co.，LTD，Shijiazhuang 050021）
（2.National Engineering Laboratory of Flat Panel Display Glass Technology and Equipment，Shijiazhuang 050035）

Ansys-Fluent，a numerical analysis software for hydrodynamics，is used to calculate the temperature field and flow field of the platinum channel.The“user-defined function（udf）”of the software is used to write the program to automatically adjust the heating power of the platinum channel.The temperature of the platinum channel reaches the intended target.From the analysis of the results，it is concluded that the rise of the bubble near the wall is favorable in the platinum channel of the temperature rise，but the rising speed will decrease or even stop when the bubble rises to the middle or near the top position.In the platinum channel of the temperature decrease，The bubble rise near the wall is unfavorable，but has little effect on the bubble overflow that has risen to the middle or near the top position.In the process design，the temperature rise section pipe should be placed in front and the descending section should be placed.By comparing the heating power with the actual use of heating power，the two are basically the same，verify the reliability of the heating power numerical model.

liquid crystal substrate glass，numerical simulation，platinum channel，temperature control，heating power

TP392

10.3969/j.issn.1672-9722.2017.10.034

Class Number TP392

2017年4月8日，

2017年5月26日

國家重點研發(fā)計劃項目（編號：2016YFB0303705）資助。

李青，女，碩士，教授級高級工程師，研究方向：機械設(shè)計。汪偉軍，男，工程師，研究方向：機械設(shè)計。姜文成，男，工程師，研究方向：機械設(shè)計。王麗紅，女，高級工程師，研究方向：精密儀器設(shè)計。郝藝，女，碩士，助理工程師，研究方向：光電子材料與器件。