彭　壽，陳淑勇，陳　凱，王　蕓，彭　程

(蚌埠玻璃工業(yè)設(shè)計(jì)研究院，浮法玻璃新技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室，蚌埠　233000)

?

玻璃粉料熔融球化過(guò)程的模擬研究

彭壽，陳淑勇，陳凱，王蕓，彭程

(蚌埠玻璃工業(yè)設(shè)計(jì)研究院，浮法玻璃新技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室，蚌埠233000)

采用數(shù)值模擬方法，對(duì)單個(gè)熔融玻璃粉料顆粒的恒溫球化過(guò)程進(jìn)行了模擬研究。應(yīng)用FLUENT軟件，結(jié)合流體體積法追蹤顆粒球化過(guò)程的液滴界面，模擬得到了熔融玻璃粉料顆粒在表面張力作用下自然球化的演化過(guò)程。提出了以熔融液滴的表面積變化來(lái)表征球化進(jìn)程，分析了恒溫條件下，熔融液滴的初始球形度、體積、粘度、表面張力因素對(duì)液滴球化時(shí)間的影響，并根據(jù)模擬結(jié)果擬合得到了特征球化時(shí)間的經(jīng)驗(yàn)公式。

球化； 數(shù)值模擬； 球化時(shí)間； VOF方法

1　引　言

玻璃微珠是直徑為數(shù)微米至數(shù)百微米的玻璃質(zhì)球形粉體材料，分為實(shí)心、空心、多孔玻璃微珠[1-3]。玻璃微珠材料具有強(qiáng)度高、流動(dòng)性好、光學(xué)性能優(yōu)異等特點(diǎn)，被廣泛用作光反射材料、凈化拋光材料、填充材料、載體材料等。玻璃微珠的生產(chǎn)方法分為物理法和化學(xué)法，采用物理法制備的玻璃微珠具有更優(yōu)的性能[4-6]。生產(chǎn)中較為常用的粉末法是一種二次成形工藝，其主要工藝流程為：首先將大塊狀玻璃通過(guò)物理方法粉碎至要求的粒徑并進(jìn)行分級(jí)，然后將分級(jí)合格的玻璃粉料加入球化爐中，使玻璃顆粒加熱熔融，在表面張力作用下自然球化，最后將球化后的玻璃液滴冷卻收集，得到玻璃微珠產(chǎn)品。

玻璃微珠的粉末法制備工藝，其優(yōu)點(diǎn)為球徑易于控制、得球率高、可生產(chǎn)硬質(zhì)玻璃微珠，但也具有產(chǎn)量低、成本高、能耗大的缺點(diǎn)[6]。粉末法工藝中的關(guān)鍵設(shè)備為高溫球化爐，因此，需對(duì)粉料在球化爐中的熔融球化過(guò)程進(jìn)行深入研究，為工藝的優(yōu)化改進(jìn)指明方向。由于球化爐內(nèi)為封閉的高溫環(huán)境，實(shí)驗(yàn)觀測(cè)較為困難，因此，研究者主要采用數(shù)值模擬手段來(lái)對(duì)球化爐進(jìn)行研究，可獲得球化爐內(nèi)的流場(chǎng)、溫度場(chǎng)，進(jìn)而可計(jì)算不同粒徑的顆粒在爐內(nèi)運(yùn)動(dòng)軌跡及顆粒溫度，并定性分析顆粒的球化效率[7-10]。

玻璃顆粒在球化爐中加熱熔融后，其自然球化在很短時(shí)間內(nèi)即可完成，一般小于0.1 s，難以設(shè)計(jì)實(shí)驗(yàn)進(jìn)行觀測(cè)。文獻(xiàn)中關(guān)于液滴球化的研究集中于噴霧過(guò)程的研究。Nichiporenko與Naida[11,12]分析了噴霧法制備金屬粉末過(guò)程，得到了金屬液滴球化時(shí)間的經(jīng)驗(yàn)公式，發(fā)現(xiàn)可通過(guò)調(diào)整液滴表層的固化速度來(lái)控制粉末的形貌。

鑒于對(duì)玻璃粉料熔融球化過(guò)程實(shí)驗(yàn)觀測(cè)的困難性，本文采用數(shù)值模擬手段，對(duì)液滴球化過(guò)程進(jìn)行了研究。通過(guò)FLUENT軟件，結(jié)合流體體積法(Volume of Fluid, VOF)，對(duì)粒徑為10～30 μm的單個(gè)熔融玻璃粉料顆粒的恒溫球化過(guò)程進(jìn)行了模擬研究，不但可補(bǔ)充實(shí)驗(yàn)方法的不足，模擬研究所得到的數(shù)據(jù)也遠(yuǎn)遠(yuǎn)超過(guò)實(shí)驗(yàn)觀測(cè)。通過(guò)對(duì)熔融玻璃顆粒球化過(guò)程的模擬，本文得到了熔融玻璃顆粒在表面張力作用下的自然球化演化過(guò)程，并分析了顆粒球化時(shí)間的影響因素。

2　數(shù)值模型

2.1物理模型和簡(jiǎn)化假設(shè)

圖1　模擬計(jì)算域示意圖Fig.1　Sketch of the computational domain

本文模擬所采用的三維物理模型如圖1所示。選擇50×50×50 μm3的立方體計(jì)算域，中心區(qū)域放置一長(zhǎng)方體形熔融玻璃液滴，為平衡模擬精度與計(jì)算耗時(shí)，本文取網(wǎng)格空間步長(zhǎng)Δx=1.0 μm，即12.5萬(wàn)控制體單元。應(yīng)用FLUENT軟件，對(duì)熔融玻璃液滴在表面張力作用下的自然球化過(guò)程進(jìn)行研究，并根據(jù)模擬結(jié)果分析各因素對(duì)液滴球化過(guò)程的影響。

為簡(jiǎn)化模型，降低模擬計(jì)算耗時(shí)，本文引入如下簡(jiǎn)化假設(shè)：(a)過(guò)程恒溫；(b)玻璃液與空氣為不可壓牛頓流體；(c)忽略重力影響；(d)壓力為標(biāo)準(zhǔn)大氣壓。

2.2控制方程

根據(jù)引入的模型簡(jiǎn)化假設(shè)，本文所研究的玻璃液滴球化過(guò)程主要涉及流場(chǎng)求解、液滴界面捕捉。其中，流場(chǎng)控制方程包括質(zhì)量守恒方程(連續(xù)性方程)、動(dòng)量守恒方程(Navier-Stokes equations)，其不可壓縮形式如下：

▽·u=0

(1)

(2)

式中：u為流體速度，m·s-1；t為時(shí)間，s；ρ為流體密度，kg·m-3；p為壓力表壓，Pa；μ為流體粘度，Pa·s；F為動(dòng)量源項(xiàng)，N·m-3。

液滴界面捕捉采用VOF方法，定義α為液相體積分?jǐn)?shù)，表征液相在某控制體單元內(nèi)所占的體積分?jǐn)?shù)。根據(jù)質(zhì)量守恒定律，液相體積分?jǐn)?shù)α的控制方程如下[13]：

(3)

在VOF方法中，上述方程中的相關(guān)流體性質(zhì)可由下式計(jì)算：

ρ=(1-α)ρG+αρL

(4)

μ=(1-α)μG+αμL

(5)

式中：下標(biāo)L、G分別代表液相、氣相。

本文所模擬液滴球化過(guò)程中，動(dòng)量方程(2)中的源項(xiàng)F為表面張力在液滴界面兩側(cè)產(chǎn)生的附加壓差。根據(jù)離散理論，將該壓差表示為作用于控制體的體積力源項(xiàng)，其形式為：

(6)

式中：σ為液體表面張力，N·m-1；κ為氣液界面曲率，m-1。對(duì)于表面張力造成的源項(xiàng)，本文采用基于VOF方法的連續(xù)表面張力模型(continue surface force, CSF)，則界面曲率可表示為[14]：

n=▽?duì)?/p>

(7)

式中：n為氣液界面法向矢量。

2.3模型設(shè)置

模擬中所采用的物性參數(shù)如表1所示，其數(shù)值為球化爐實(shí)際操作溫度下(取為1600 ℃)的估算值，其中熔融玻璃液粘度、表面張力為一虛擬范圍，用以研究玻璃液物性對(duì)球化過(guò)程的影響。

表1　數(shù)值模擬中采用的流體物性參數(shù)

邊界條件設(shè)置：如圖1所示的計(jì)算域外邊界采用壓力入口邊界，p=0 Pa。

初始條件設(shè)置：如圖1所示，初始化一長(zhǎng)方體狀靜止熔融玻璃液滴于計(jì)算域中心，初始流速及壓力為：u0=0 m·s-1，p0=0 Pa。

數(shù)值方法設(shè)置：采用PISO(Pressure-Implicit with Splitting of Operators)方法處理速度—壓力耦合求解，采用CICSAM(Compressive Interface Capturing Scheme for Arbitrary Meshes)方法處理VOF方程，以獲取銳利的相界面。

本文所模擬的液滴球化過(guò)程為非穩(wěn)態(tài)過(guò)程，在初始時(shí)刻，液滴球形度小，表面張力源項(xiàng)最大，將導(dǎo)致較大的流體速度，需要較小的時(shí)間步長(zhǎng)以保證數(shù)值收斂性及計(jì)算精度；球化過(guò)程后期，液滴球形度增大，表面張力源項(xiàng)變小，流體速度變小，采用較大的時(shí)間步長(zhǎng)仍可保證數(shù)值收斂性，并可顯著降低計(jì)算耗時(shí)。因此，本文采用可變時(shí)間步長(zhǎng)設(shè)置：全局庫(kù)朗特?cái)?shù)(Courant number)設(shè)為0.5，最小時(shí)間步長(zhǎng)為Δt=1×10-8s，最大時(shí)間步長(zhǎng)為Δt=1×10-5s。

3　模擬結(jié)果與討論

3.1液滴球化過(guò)程模擬結(jié)果

采用圖1所示物理模型，按照前文所述模型假設(shè)及設(shè)置，對(duì)液滴球化過(guò)程進(jìn)行模擬研究，模擬條件為：玻璃液滴初始形狀為20×20×20 μm3的立方體；玻璃液物性為μL=1 Pa·s，σL=0.25 N·m-1。模擬結(jié)果如圖2所示。

由圖2所示模擬結(jié)果可知，在初始階段(圖2a)，立方體形液滴的邊角處，具有較大的界面曲率κ，因此表面張力在邊角處產(chǎn)生較大的附加壓差，即方程(2)中的動(dòng)量源項(xiàng)F較大，使得邊角處的液體具有較大的指向液滴內(nèi)部的速度；隨著液滴在表面張力作用下逐漸球化(圖2b)，液滴表面各處的曲率趨于一致，液滴表面處的壓力趨于一致，液體流速在粘性耗散作用下變小，液滴繼續(xù)球化，但球化速度變緩；隨著時(shí)間繼續(xù)，液滴最終將變?yōu)榍蛐?，液滴?nèi)部壓力一致，液體流速趨于0，液滴球化完成。

熔融玻璃液滴球化過(guò)程中，液滴內(nèi)外平均壓力差以及液滴內(nèi)平均速度隨時(shí)間變化如圖3所示。如前文所分析，初始階段液滴內(nèi)部具有較高的壓力及流速；隨著球化進(jìn)行，液滴內(nèi)部平均壓力及速度均呈指數(shù)衰減趨勢(shì)。當(dāng)達(dá)到穩(wěn)定時(shí)，液滴內(nèi)部流速趨于0，壓力趨于定值。由圖3中結(jié)果可發(fā)現(xiàn)，液滴球化完成后，液滴內(nèi)外平均壓力差模擬值與理論計(jì)算值存在一定誤差，這是由于數(shù)值解存在一定的誤差導(dǎo)致。

圖2　液滴外形、x軸中心面壓力及速度矢量分布在不同時(shí)刻的模擬結(jié)果Fig.2　Simulation results of droplet shape, pressure and velocity vector at the middle

圖3　玻璃液滴球化過(guò)程內(nèi)外平均壓力差及平均速度變化趨勢(shì)圖Fig.3　Charts of mean pressure difference andmean velocity versus time of the spheroidization process

圖4　熔融玻璃液滴球化過(guò)程表面積隨時(shí)間變化趨勢(shì)圖Fig.4　Chart of surface area versus time of the moltendroplet spheroidization process

對(duì)于液滴的球化過(guò)程，本文提出采用液滴表面積S隨時(shí)間t的變化來(lái)表征球化的進(jìn)度，模擬所得S-t結(jié)果如圖4所示，可發(fā)現(xiàn)球化過(guò)程中S隨t呈指數(shù)衰減趨勢(shì)。為了進(jìn)一步對(duì)液滴球化過(guò)程進(jìn)行量化分析，本文采用如下指數(shù)衰減函數(shù)對(duì)模擬所得S-t數(shù)據(jù)進(jìn)行擬合：

S=ΔSe-t/τ+S∞

(8)

式中：S∞、ΔS、τ為待定擬合系數(shù)。針對(duì)本文所研究液滴球化過(guò)程，S∞為玻璃液滴的球化后表面積，m2；ΔS為液滴初始狀態(tài)與球化后的表面積差，m2；τ為特征球化時(shí)間，s，表征液滴球化過(guò)程的快慢。物理破碎得到的粉料顆粒球形度一般大于0.5，故若以液滴外表面積與等體積球體表面積相對(duì)偏差小于0.1%作為液滴球化完全的標(biāo)準(zhǔn)，則根據(jù)公式，液滴完全球化所需時(shí)間約為10τ。

根據(jù)S-t數(shù)據(jù)擬合公式，得到圖4中模擬數(shù)據(jù)的擬合結(jié)果如下：

S=3.3496×10-10×e-t/9.511×10-6+1.9223×10-9

(9)

即，粘度為1 Pa·s、形狀為20×20×20 μm立方體的熔融玻璃液滴，其在表面張力作用下完全球化所需時(shí)間約為1×10-4s。

3.2球化過(guò)程影響因素分析

根據(jù)S-t數(shù)據(jù)擬合公式可知，參數(shù)S∞、ΔS表征了液滴的形狀特征，而參數(shù)τ表征了液滴球化過(guò)程的速度，因此可以通過(guò)研究τ來(lái)分析液滴球化過(guò)程中各影響因素的作用。

熔融液滴球化過(guò)程的主要影響因素為液體性質(zhì)、液滴初始形狀，其中，液體性質(zhì)因素包括粘度μL、表面張力σL，液滴初始形狀因素包括液滴初始球形度φ0、液滴體積VL。上述各影響因素相互獨(dú)立，為簡(jiǎn)化分析的復(fù)雜性，本文進(jìn)行了上述各參數(shù)的單因素影響分析。

(1)粘度影響分析

對(duì)粘度進(jìn)行單因素影響分析，固定其他因素參數(shù)：σL=0.25 N·m-1，VL=8.0×10-15m3，φ0=0.806。根據(jù)模擬所得S-t數(shù)據(jù)，擬合得到不同粘度下的S∞、ΔS、τ數(shù)據(jù)如表2所示。

表2　不同粘度條件下模擬數(shù)據(jù)擬合結(jié)果

由表2數(shù)據(jù)可知，在液滴形狀參數(shù)VL、φ0固定的情況下，參數(shù)S∞、ΔS的擬合值基本一致，也證明了公式的可靠性。根據(jù)表中結(jié)果，可發(fā)現(xiàn)擬合所得特征球化時(shí)間τ隨液體粘度變化顯著，且呈線性趨勢(shì)，如圖5所示。

采用圖5中所示的線性函數(shù)，對(duì)τ-μL數(shù)據(jù)進(jìn)行擬合，所得到的τ-μL關(guān)系擬合函數(shù)如下：

τ=1.0322×10-5μL

(10)

(2)表面張力影響分析

對(duì)表面張力進(jìn)行單因素影響分析，固定其他因素參數(shù)：μL=10 Pa·s，VL=8.0×10-15m3，φ0=0.806。為簡(jiǎn)化篇幅，未給出不同表面張力條件下模擬數(shù)據(jù)擬合結(jié)果表。擬合所得特征球化時(shí)間τ與液體表面張力σL關(guān)系如圖6所示。

圖5　特征球化時(shí)間τ與液體粘度μL關(guān)系圖Fig.5　Relationship of characteristic spheroidizingtime τ versus liquid viscosity μL

圖6　特征球化時(shí)間τ與液體表面張力σL關(guān)系圖Fig.6　Relationship of characteristic spheroidizingtime τ versus liquid surface tension σL

根據(jù)圖6所示，特征球化時(shí)間τ與液體表面張力σL呈反比趨勢(shì)，按照?qǐng)D中的擬合公式，擬合所得τ-σL結(jié)果如下：

(11)

(3)液滴初始球形度影響分析

對(duì)液滴初始球形度進(jìn)行單因素影響分析，固定其他因素參數(shù)：μL=1 Pa·s，σL=0.25 N·m-1，VL=8.0×10-15m3。所得特征球化時(shí)間τ與液滴初始球形度φ0關(guān)系如圖7所示。

根據(jù)圖7所示，特征球化時(shí)間τ與(1-φ0)呈指數(shù)衰減關(guān)系，按照?qǐng)D中的擬合公式，擬合所得τ-(1-φ0)結(jié)果如下：

τ=-7.4049×10-5e-(1-φ0)/0.1239+2.5036×10-5

(12)

對(duì)于液滴球化過(guò)程，只有當(dāng)初始球形度φ0=1時(shí)，球化特征時(shí)間參數(shù)τ才趨于0，而所得擬合公式，在φ0=0.87時(shí)，τ即趨于0。出現(xiàn)該偏差的原因可能有：(a)液滴球化過(guò)程復(fù)雜，其初始形狀的影響難以單一地采用球形度φ進(jìn)行描述；(b)模擬采用控制方程的數(shù)值解來(lái)分析液滴球化過(guò)程，具有一定的數(shù)值誤差。該結(jié)果與理論分析雖然存在一定偏差，但在所本文所研究范圍內(nèi)，可認(rèn)為該公式可靠。

(4)液滴體積影響分析

對(duì)液滴體積進(jìn)行單因素影響分析，固定其他因素參數(shù)：μL=1 Pa·s，σL=0.25 N·m-1，φ0=0.806。所得特征球化時(shí)間τ與液滴體積VL關(guān)系如圖8所示。

圖7　特征球化時(shí)間τ與液滴初始球形度φ0關(guān)系圖Fig.7　Relationship of characteristic spheroidizingtime τ versus initial sphericity φ0

圖8　特征球化時(shí)間τ與液滴體積VL關(guān)系圖Fig.8　Relationship of characteristic spheroidizingtime τ versus droplet volume VL

根據(jù)圖8所示，特征球化時(shí)間τ與液滴體積VL呈冪函數(shù)關(guān)系，采用圖中的公式，擬合所得τ-VL結(jié)果如下：

(13)

式中：VL1/3表征液滴的特征尺寸。

圖9　液滴特征球化時(shí)間τ公式擬合結(jié)果Fig.9　Fitting results of characteristic spheroidizing time τ

(14)

式中：a1、a2、a3為方程待定系數(shù)。采用非線性最小二乘擬合，得到公式中的待定參數(shù)為：a1=0.9262，a2=0.1312，a3=0.1277。

由τ公式擬合結(jié)果可知，所提出的液滴球化過(guò)程特征時(shí)間τ與影響因素μL、σL、φ0、VL的經(jīng)驗(yàn)擬合公式，在本文所研究范圍內(nèi)具有良好的可靠性。該結(jié)果對(duì)液滴的球化過(guò)程的定量研究具有指導(dǎo)意義。

4　結(jié)　論

本文采用數(shù)值模擬方法，對(duì)粉末法制備玻璃微珠工藝中，熔融玻璃粉料在表面張力作用下的自然球化過(guò)程進(jìn)行了模擬研究，得到了恒溫條件下熔融液滴自然球化過(guò)程中的速度、壓力隨時(shí)間的演化，并以液滴的表面積變化表征了液滴球化進(jìn)程；根據(jù)模擬結(jié)果，分析了液滴球化過(guò)程中相關(guān)影響因素的作用，并得到了液滴球化過(guò)程特征時(shí)間與影響因素間的經(jīng)驗(yàn)擬合公式。對(duì)深入了解高溫熔融球化過(guò)程，優(yōu)化生產(chǎn)工藝提供了參考依據(jù)。

本文對(duì)球化液滴球化過(guò)程的影響因素的分析研究，簡(jiǎn)化為單因素分析，忽略了各影響因素間的交互作用，且實(shí)際粉料球化過(guò)程涉及粉料顆粒與環(huán)境的換熱，因此需要進(jìn)一步的深入研究，以對(duì)玻璃粉料熔融球化過(guò)程有更深入、準(zhǔn)確的認(rèn)識(shí)。

[1] 袁啟華,羅大兵,雷麗文.載體多孔玻璃微珠的制備及其應(yīng)用[J].硅酸鹽通報(bào),2000,19(3):52-55.

[2] 楊玉香,邵謙,葛圣松.玻璃微珠的應(yīng)用研究進(jìn)展[J].中國(guó)粉體技術(shù),2006,12(2):45-48.

[3] 彭壽,王蕓,彭程，等.空心玻璃微珠制備方法及應(yīng)用研究進(jìn)展[J].硅酸鹽通報(bào),2012,31(6):1508-1513.

[4] 施陽(yáng)和,鄭華.球形粉體的制備方法及應(yīng)用[J].中國(guó)粉體技術(shù),2015,21(4):71-75.

[5] 宋曉睿,楊輝.空心玻璃微球制備技術(shù)研究進(jìn)展[J].硅酸鹽學(xué)報(bào),2012,40(3):450-457.

[6] 石成利,梁忠友,李春紅.玻璃微珠的制備方法[J].山東建材,2005,26(1):33-35.

[7] Tsymbalist M M,Rudenskaya N A,Kuz’min B P,et al.Low-temperature plasma spheroidizing of polydisperse powders of refractory materials[J].ProtectionofMetals,2003,39(4):338-343.

[8] Lu X,Zhu L P,Zhang B,et al.Simulation of flow field and particle trajectory of radio frequency inductively coupled plasma spheroidization[J].ComputationalMaterialsScience,2012,65:13-18.

[9] 邵宗恒.高頻等離子體球化石英粉的數(shù)值模擬研究[D].上海：華東交通大學(xué)碩士學(xué)位論文,2012.

[10] Ramasamy R,Selvarajan V.Injected particle behavior in a thermal plasma[J].AppliedSurfaceScience,2001,169-170:617-621.

[11] Nichiporenko O S,Naida Y I.Heat exchange between metal particles and gas in the atomization process[J].SovietPowderMetallurgyandMetalCeramics,1968,7(7):509-512.

[12] Nichiporenko O S,Naida Y I.Fashioning the shape of sprayed powder particles[J].PowderMetallurgyandMetalCeramics,1968,7(10):753-755.

[13] Hirt C W,Nichols B D.Volume of fluid (VOF) method for the dynamics of free boundaries[J].JournalofComputationalPhysics,1981,39(1):201-225.

[14] Brackbill J U,Kothe D B,Zemach C.A continuum method for modeling surface tension[J].JournalofComputationalPhysics,1992,100(2):335-354.

Simulation Research of Glass Powder Particle Melt Spheroidizing Process

PENGShou,CHENShu-yong,CHENKai,WANGYun,PENGCheng

(State Key Laboratory for Advanced Technology of Float Glass,Bengbu Design & Research Institute for Glass Industry,Bengbu 233000,China)

The isothermal spheroidizing process of a single molten glass particle was investigated via numerical simulation. By utilizing the FLUENT software, the evolution of molten glass particle natural spheroidizing process induced by surface tension effect was simulated, which based on tracking the liquid interface of spheroidizing process by volume of fluid method. The spheroidizing process is characterized by the variation of the molten droplet’s surface area. Additionally, the influence factors of spheroidizing time were also analyzed, including the initial sphericity, volume, viscosity, surface tension of the molten droplet. Finally, an empirical formula of characteristic spheroidizing time was obtained with the simulation results.

spheroidization;numerical simulation;spheroidizing time;VOF method

彭壽(1960-)，男，教授級(jí)高工.主要從事浮法玻璃新技術(shù)、特種玻璃、光電顯示玻璃等方面的研究.

陳淑勇,博士研究生.

TQ170.1

A

1001-1625(2016)06-1800-07