質子交換爐溫度場均勻性分析與優化

2019-05-05 02:17:34伏娜張晞

北京航空航天大學學報 2019年4期

伏娜，張晞

(北京航空航天大學儀器科學與光電工程學院, 北京 100083)

鈮酸鋰(LiNbO3)晶體具有優異的電光和聲光等特性，被廣泛用于制備各種光波導器件[1-2]。LiNbO3光波導是高精度光纖陀螺系統的核心器件，LiNbO3光波導的大批量生產是實現光纖陀螺產業化的關鍵之一。質子交換爐是用質子交換法制備LiNbO3光波導的主要設備，為了提高生產效率，往往將批量LiNbO3基底同時放到質子交換爐中進行質子交換以制備光波導。若爐內軸向溫度分布不均勻，將導致不同溫度梯度下進行質子交換時LiNbO3光波導的薄膜厚度會有所不同，會對同一批次制備的光波導的性能一致性造成影響，因此提高質子交換爐的溫度場均勻性很有必要[3]。目前對質子交換爐的要求為：質子交換爐軸向均勻溫區長度≥200 mm，均勻溫區內的最大溫差在5 min內≤1℃。

目前國內外對質子交換爐內部溫度的研究重點大多在于溫度的精確控制[4-5]，而對爐內溫度分布均勻性的相關研究較少，導致現有的質子交換爐很難用于大批量LiNbO3光波導的制備。質子交換爐內存在熱傳導、熱對流等多種傳熱方式，十分復雜，難以建立精確的數學模型。面對復雜的質子交換爐系統設計，一般采用模仿或經驗改進等方法對質子交換爐進行設計，不僅成本高而且研究周期長，生產出來的質子交換爐往往不能達到理想的溫度分布。隨著計算機技術的發展，計算流體力學(CFD)的出現給質子交換爐的設計提供了極大的方便[6-7]。FLUENT是國際上流行的CFD軟件包，它提供了用戶自定義函數(UDF)作為FLUENT的二次開發，當標準的FLUENT界面不能滿足用戶需求時，用戶可以編寫UDF程序并動態地連接到FLUENT求解器上以實現期望的控制算法、模型方程等[8]。1994年，Chicatelli等[9]首次將CFD和控制方法相結合，用于高速推進系統的模擬和控制。Yang等[10]針對垃圾焚燒爐難以測量內部溫度的問題，通過CFD建模研究前饋控制效果，并將數值模擬結果作為控制系統信息數據庫。Gao等[11]針對PID控制器參數整定時難以評估整體控制效果的問題，提出了一種基于FLUENT的PID控制器參數整定評估方法，驗證了該方法的有效性。這些成果為本文模擬質子交換爐的控溫奠定了一定的基礎。本文采用基于FLUENT的PID控制方法，利用FLUENT模擬控制對象，編寫UDF控溫程序，使用FLUENT UDF進行閉環控制，模擬質子交換爐爐溫的PID控制。與傳統的PID控制器不同之處在于基于FLUENT的PID控制器不需要建立被控對象的傳遞函數模型，而是建立被控對象的網格模型，利用FLUENT對網格計算求解得到最終的控制效果。

本文基于質子交換爐提出了多種加熱控溫方案，利用FLUENT對質子交換爐加熱過程的溫度場進行仿真，分析此過程中溫度場分布特點，研究質子交換爐控溫方式、傳感器位置布置、加熱絲布置對爐內溫度場均勻性的影響，找到最佳加熱方案，提高爐內溫度場均勻性。

1 模型建立

1.1 物理模型

本文研究的質子交換爐采用圓筒結構，立式放置，爐體結構由外至內分別有保溫層、加熱絲、爐壁，如圖1所示，爐體高H=0.8 m，爐膛內徑為0.12 m。加熱絲均勻纏繞在爐壁四周，將產生的熱量傳遞給爐壁，爐壁與爐膛內的空氣主要以傳導和對流的方式進行熱量傳遞，由于冷熱空氣密度不同，爐膛內的熱氣流會向上流動。爐體壁面和底部有保溫結構，為方便工藝過程中打開爐口送料，爐頂一般沒有保溫結構，故爐頂處主要以對流

的方式與外界環境進行換熱。圖1中爐膛內的箭頭代表熱量流動。

質子交換爐的爐溫是通過改變加熱絲的控制電流來調節的，綜合考慮控制精度、爐溫均勻性等多方面因素，可將加熱絲分為多段，每段加熱絲對應一個溫區，構成多個加熱回路，并在每個區安裝溫度傳感器(通常采用熱電偶)用于反饋各區溫度，各區的控制回路根據目標設定溫度值與實時反饋的溫度值之差調整相應的控制器輸出以控制加熱絲加熱, 如圖1所示，ri為設定溫度值，yi為反饋溫度值，ei為溫度偏差，PIDi代表控制器，ui為控制器的輸出值(其中i為控制回路的個數，i=1,2,…,n)。

由于本文的研究重點在于質子交換爐內部的溫度場分布，故對實際模型作如下簡化：

1) 爐膛內的溫度是本文的研究重點，故將模型簡化為內徑0.12 m，高0.8 m的圓柱體。

2) 加熱絲均勻纏繞在爐膛周圍，可假設壁面熱流量是均勻的。

質子交換工藝的目標溫度為240℃，安裝在爐壁外表面的溫度傳感器將采集的溫度信號發送給控制器，控制器將采集溫度值與目標溫度值進行比較后經PID計算輸出控制量，并由繼電器控制加熱絲的電流通斷來調節發熱量，使爐膛的溫度穩定在目標溫度值。

1.2 數學模型

1.2.1 質子交換爐內空氣的傳熱模型

爐膛內的空氣和爐壁之間的熱量傳遞主要靠熱傳導和熱對流的方式進行，由能量守恒原理和傅里葉定律可得柱坐標系下空氣的傳熱微分方程[12]：

(1)

為簡化式(1)，現做如下合理假設：①由于空氣是完全透明體，因此可以忽略輻射換熱；②空氣是不可壓縮的流體，空氣的流速是固定的；③徑向空氣流量可忽略不計；④內熱源生成熱為0。

實際上空氣的流速在徑向上不為零，但是空氣在軸向上的速率比徑向上快的多，因此可忽略不計。上述假設下的空氣瞬態熱傳導模型如下：

(2)

對于軸對稱的爐體，變量T與φ無關，因此上述模型可簡化為

(3)

初始條件和邊界條件為

(4)

式中：h為表面對流傳熱系數；TH為爐體頂部的溫度；Tf為周圍流體的溫度，為已知量；T0為初始溫度；Tm為空氣溫度。

1.2.2 CFD模型

CFD模型是基于微元內的質量守恒、動量守恒和能量守恒確定的一組偏微分方程。基本控制方程包括質量守恒方程、動量守恒方程和能量守恒方程，數學表達式為[13]

質量守恒方程：

(5)

x方向動量守恒方程：

(6)

y方向動量守恒方程：

z方向動量守恒方程：

(7)

能量守恒方程：

(8)

式中：U=(u,v,w)，u、v、w分別為U在x、y、z方向的分速度；Su、Sv、Sw、ST均為廣義源項；k為傳熱系數；μ為動力黏度系數;p為壓強。

2 FLUENT仿真

2.1 網格劃分及邊界條件設置

進行CFD分析的第1步是劃分網格[14]。其主要思想是將空間連線的計算區域分割成足夠小的計算區域，然后在每一計算區域里應用流體控制方程，求解所有區域的流體計算方程，最終獲得整個計算區域的物理量分布[15]。本文針對不同的方案分別采用Gambit生成圓柱體的結構化網格。

邊界條件的設置如表1所示。

表1中的UDF flux表示使用UDF編寫的熱流量輸出函數。爐壁的邊界條件分為2種情況，有加熱絲的部分設為壁面wall，熱流量通過調用UDF編寫的熱流量輸出函數來獲得；沒有加熱絲的部分設為壁面wall，由于有保溫層，熱流量設為0。

表1 邊界條件的設置Table 1 Boundary condition setting

2.2 模擬監測點設置

質子交換爐內的軸向溫度分布是影響LiNbO3光波導生產質量的關鍵。由于研究重點是質子交換爐內軸向方向的溫度場分布，爐口沒有保溫結構，受外界環境影響較大，故均勻溫區在爐底。工藝要求均勻溫區長度Δh≥200 mm，故在均勻溫區從爐底依次沿軸向方向均勻布置了5個監測點，間距l=0.05 m，其位置坐標分別為(0，0，0)m、(0，0，0.05) m、(0，0，0.1)m、(0，0，0.15)m、(0，0，0.2)m。質子交換爐內均勻溫區溫度分布用這5個測點的溫度分布來描述。

2.3 求解器的設置

激活能量模型，湍流模型選擇標準k-ε模型，壁面函數選為標準壁面函數，求解方法選擇SIMPLEC，湍流動能與湍流耗散選擇二階差分，湍流黏度和能量的松弛因子選為0.8，其他保持默認值。質子交換爐內空氣初始溫度設為300 K，預期達到的溫度設為513 K，根據溫度傳感器采集到的溫度，控制壁面熱流量的變化。對于瞬態求解，時間步長是一個重要的概念，時間步長設置太大會導致計算難以收斂，設置太小會增加迭代次數，計算時間過長。為了盡量減少計算時間而不降低計算精度，在仿真中設時間步長為0.5 s，時間步數為800，即仿真時間為400 s。始終采用上述設置完成以下仿真內容。

3 不同方案下的仿真結果與對比分析

本文提出3種控溫方式，分別為整段控溫、兩段控溫、三段控溫。整段控溫即加熱絲為一整段，傳感器位置布置在整段加熱絲上，由其反饋的溫度經過PID處理后控制整段加熱絲的熱流量；兩段控溫即將均勻纏繞在爐膛周圍的加熱絲均勻分為2段，2個傳感器位置分別布置在每段加熱絲上，其反饋的溫度經過PID處理后分別控制每段加熱絲的熱流量；三段控溫即將均勻纏繞在爐膛周圍的加熱絲均勻分為3段，3個傳感器位置分別布置在每段加熱絲上，其反饋的溫度經過PID處理后分別控制每段加熱絲的熱流量。

本文在提出的3種控溫方式下分別設計了不同的傳感器位置布置，如圖2所示。在每種傳感器位置布置下通過改變加熱絲布置高度h，即h=2Δh、h=3Δh，可得出表2中18種設計方案，利用FLUENT分別模擬這18種方案下的溫度場。

為準確描述溫度場的均勻性，采用在同一時刻所有監測點溫度的最高溫度與最低溫度之差，即最大偏差來衡量溫度場均勻性的好壞，其表達式為

ΔTmax=yimax-yimin

(9)

式中：yimin為第i時刻所有監測點的最低溫度值；yimax為第i時刻所有監測點的最高溫度值。

對以上18種方案分別仿真400 s，經過100 s后爐內溫度基本上達到動態平衡，本文選取100 ～400 s之間的最大溫差ΔTmax來分析各種方案的均勻性好壞，所得結果如圖3～5所示。

對比圖3(a)和(b)可以發現，采用整段控溫方式時改變傳感器的位置對加熱絲布置高度h=2Δh的均勻性影響較大，對于h=3Δh和h=4Δh的均勻性幾乎沒有影響；由圖3(a)可以看出，傳感器位置在z=0時，加熱絲布置高度h=2Δh時爐內溫度場均勻性最差，而h=4Δh時爐內溫度場均勻性最好。這是由于加熱絲布置的高度過低時，沒有布置加熱絲的爐壁熱流量為0，不能完全補償爐頂的對流換熱作用引起的熱量散失，導致爐底均勻溫區長度變短，達不到期望的均勻溫區長度。因此整段控溫方式下的最佳方案為方案3，即傳感器位置z=0，加熱絲布置高度h=4Δh。

圖2 3種控溫方式的傳感器位置布置Fig.2 Sensor position arrangement of three kinds of temperature control methods

對比圖4(a)和(b)可以看出，采用兩段控溫方式時改變傳感器位置對加熱絲布置高度h=2Δh的均勻性影響很大，與組合方式1相比，采用組合方式2時h=2Δh的均勻性提高了2～3倍，對于h=3Δh和h=4Δh的均勻性也有一定提高；由圖4(b)可以看出，傳感器位置布置均采用組合方式2時，加熱絲布置高度h=2Δh時爐內均勻溫區的均勻性最差，而h=4Δh時爐內均勻溫區的均勻性最好。因此兩段控溫方式下的最佳方案為方案12，即傳感器位置布置組合方式2，加熱絲布置高度h=4Δh。

表2 18種設計方案Table 2 18 kinds of design schemes

圖3 整段控溫時傳感器位置對質子交換爐均勻性的影響Fig.3 Effect of sensor position on uniformity of proton exchange furnace under only one temperature controller

對比圖5(a)和(b)可以看出，采用三段控溫方式時，傳感器位置布置采用組合方式2可以更好地改善均勻溫區的均勻性；圖5(b)中傳感器位置布置均采用組合方式2時，加熱絲布置高度h=4Δh時爐內溫度場均勻性最好。因此三段控溫方式下的最佳方案為方案18，即傳感器位置布置組合方式2，加熱絲布置高度h=4Δh。

圖4 兩段控溫時傳感器位置對質子交換爐均勻性的影響Fig.4 Effect of sensor position on uniformity of proton exchange furnace under two temperature controllers

圖5 三段控溫時傳感器位置對質子交換爐均勻性的影響Fig.5 Effect of sensor position on uniformity of proton exchange furnace under three temperature controllers

對以上3種控溫方式下的最佳方案進行進一步對比，取400 s時刻y=0截面的溫度分布云圖，如圖6所示，從左到右依次為方案18、方案12、方案3對應溫度范圍512.5～513.5 K的溫度分布，可見三段控溫下的最佳方案對應的均勻溫區長度最長，整段控溫下的最佳方案對應的均勻溫區最短。

在100～400 s時間段內對3種控溫方式下的最佳方案進行最大偏差的對比，如表3所示。

表3中，采用方案18爐內均勻溫區的均勻性最好，即三段控溫、傳感器位置布置采用組合方式2、加熱絲布置高度h=4Δh時最大溫度偏差最小，ΔTmax可達0.03℃，滿足生產工藝的要求。由表3可以看出，當加熱絲布置高度均為h=4Δh時，隨著控溫分段數目的增加，均勻溫區的最大溫度偏差越小。這是由于爐內溫度分布存在梯度時，整個壁面熱流量的不同位置需要不同程度的調整，而整段控溫會同時改變整個壁面熱流量，且改變程度相同，這就導致有些位置會有熱量積累而有些位置散失的流量得不到補償，爐內溫度分布均勻性較差；隨著控溫分段數目的增加，爐內溫度分布存在梯度時，壁面熱流量可以分時進行不同程度的調整，快速補償散失的熱量，因此爐內溫度分布越均勻。但是控溫段數越多，生產成本也會越高，而且各個控制回路之間的耦合效應也越來越嚴重，提高爐溫均勻性的同時有可能會降低控溫精度。因此滿足生產工藝均勻性要求的前提下，三段控溫的方式較好。