馬英杰,郝繼光

(北京理工大學(xué) 宇航學(xué)院，北京 100081)

【基礎(chǔ)理論與應(yīng)用研究】

水液滴撞擊導(dǎo)流面成膜過程數(shù)值模擬

馬英杰,郝繼光

(北京理工大學(xué) 宇航學(xué)院，北京 100081)

對液滴撞擊親水導(dǎo)流壁面后的液膜形成過程進行了理論分析，用基于CLSVOF方法的數(shù)值模型，分析了液滴撞擊親水表面后的鋪展特性，研究了不同表面張力系數(shù)和液滴直徑對液滴撞擊親水表面潤濕直徑的影響。結(jié)果表明， CLSVOF方法可以用于液滴撞擊親水導(dǎo)流壁面的模擬；液滴撞擊壁面達到最大潤濕直徑的時間隨著表面張力系數(shù)的增大而減?。灰旱巫矒粼谟H水壁面上的最大潤濕直徑隨著液滴初始直徑的增大而增大，液滴完成一次鋪展收縮過程需要的時間也隨液滴初始直徑的增大而增加；為進一步研究液滴撞擊導(dǎo)流壁面成膜對導(dǎo)流壁面的保護作用提供了理論基礎(chǔ)。

液滴撞擊；CLSVOF；表面張力系數(shù)；液滴直徑

航天發(fā)射的很多場合都有注水到高溫超音速氣體射流減少其伴生危害的應(yīng)用，如航天發(fā)射場注水降低燃氣射流溫度和噪聲，導(dǎo)彈熱發(fā)射時注水降低燃氣射流及反濺流輻射的紅外信號，燃氣蒸汽式彈射器內(nèi)注水降低燃氣射流溫度等。研究液滴撞擊親水表面過程對于了解液膜降溫效果具有重要意義，液滴撞擊導(dǎo)流或發(fā)射筒壁面，這些壁面均為親水表面，研究液滴撞擊導(dǎo)流面即研究液滴撞擊親水壁面的鋪展過程。液滴撞擊壁面鋪展是多相流的重要基礎(chǔ)過程，是自Worthington[1]首次研究液滴撞擊壁面一個多世紀以來，研究者感興趣的研究對象。影響液滴撞擊壁面后的運動和形態(tài)的因素有很多，主要有液滴的物性(密度、黏度、表面張力系數(shù)、尺寸)、壁面條件、液滴與壁面間的相互作用、液滴碰撞速度、撞擊方向、環(huán)境壓力和溫度等因素。本文主要對液滴的物性參數(shù)中的表面張力系數(shù)和液滴直徑對液滴撞擊壁面的影響進行了數(shù)值模擬研究，闡明了液滴撞壁鋪展過程的機理，為其在航天發(fā)射領(lǐng)域的應(yīng)用提供了理論依據(jù)和研究方法。

近年來，國內(nèi)外研究者對液滴撞壁鋪展過程進行了大量研究。實驗研究方面，1876年Worthinton[1]首先實驗研究了液滴撞擊金屬盤子的過程，后來又研究了牛奶和水銀液滴從不同高度落下撞擊金屬盤子的過程；Mao[2]研究了撞擊參數(shù)對最大鋪展直徑和反彈的影響，得出最大鋪展直徑與Re數(shù)、We數(shù)和靜態(tài)接觸角的關(guān)系；Sikalo等[3-4]實驗研究了不同液滴速度、直徑、黏性和表面張力系數(shù)的液滴撞擊石蠟、玻璃壁面鋪展直徑、速度和高度的變化規(guī)律；Wang等[5]實驗研究了液滴撞壁鋪展高度和鋪展直徑隨時間的變化，得出了不同液滴屬性、速度和壁面濕潤性對液滴鋪展情況不同的影響程度；Radulovic等[6]、Hamraoui等[7]和Fischer等[8]分別研究了pH值、高濃度和低濃度表面活性劑在固體壁面或預(yù)置液膜表面上的鋪展現(xiàn)象；畢菲菲等[9]用蒸餾水、無水乙醇和甘油3種材料的液滴做了液滴碰撞水平固體壁面的實驗，研究了液滴屬性對液滴碰壁鋪展不同時期的影響；李維仲等[10]實驗研究了不同液滴半徑、壁面親疏水性和粗糙度對液滴撞壁鋪展過程的影響。數(shù)值模擬研究方面，M.Pasandideth- Fard等[11]用實驗并基于VOF方法進行了數(shù)值模擬，研究了水液滴撞擊不銹鋼表面的鋪展過程，結(jié)果與實驗結(jié)果有很好的一致性；Lunkad等[12]數(shù)值模擬研究了碰撞表面的傾斜度、潤濕特性、液滴屬性和碰撞速度對液滴撞擊平壁面和傾斜壁面的鋪展動態(tài)過程的影響；Gunjal[13]應(yīng)用VOF方法數(shù)值模擬了不同液滴速度、黏度、表面張力和尺寸等對液滴鋪展系數(shù)和鋪展高度的影響，與實驗結(jié)果基本一致；Kim等[14]用VOF方法數(shù)值模擬研究了液滴撞擊速度和表面張力等參數(shù)對非牛頓流體撞擊固體壁面的動態(tài)特性的影響。劉邱祖等[15]用格子Boltzmann方法模擬研究了液滴最大鋪展直徑、壁面親疏水性與液滴表面張力和鋪展速度之間的關(guān)系；周龍玉等[16]用實驗方法結(jié)合數(shù)值模擬研究了液滴碰壁鋪展和反彈的現(xiàn)象，得出了液滴表面張力系數(shù)和液滴半徑對其鋪展過程的影響；陳石等[17]通過理論分析和數(shù)值模擬相結(jié)合對液滴受到的慣性力、黏性力和表面張力進行了理論分析，分析了液滴物性參數(shù)對液滴撞壁鋪展的影響。

本文采用CLSVOF方法對液滴撞擊導(dǎo)流壁面成膜過程進行了數(shù)值模擬研究，分析了表面張力系數(shù)和液滴初始直徑對液滴撞壁變形情況和潤濕直徑的影響，揭示了潤濕直徑隨時間的變化規(guī)律。

1 理論基礎(chǔ)

1.1 數(shù)學(xué)模型

本文采用CLSVOF方法即VOF方法和Level Set方法的耦合，該方法即能保證流體的質(zhì)量守恒又能精確計算界面的曲率和法向，克服了兩種方法的不足之處，能對相界面進行更加精確追蹤。

連續(xù)性方程為

(1)

動量方程為

(2)

表面張力模型表達式為

(3)

能量方程為

(4)

1.2 計算模型的建立

采用二維軸對稱模型進行計算，所建立的模型如圖1所示。

圖1 初始時刻模型

圖1為初始時刻的計算模型，采用純水液滴，溫度為常溫，以一定碰撞速度U0=0.77m/s垂直撞擊固體壁面，靜態(tài)接觸角為87.4°。采用有限體積法對控制方程進行離散，壓力速度耦合采用SIMPLE方法，壓力求解采用PRESTO!方法，計算單元液相體積分數(shù)采用Geo-Reconstruct方法離散求解，時間步長Δt=10-6s。

2 計算模型驗證

前文確立了基于CLSVOF方法的數(shù)值計算模型，在進行數(shù)值模擬之前，要進行網(wǎng)格劃分，網(wǎng)格數(shù)量對計算精度和計算時間有非常關(guān)鍵的影響，因此進行了網(wǎng)格無關(guān)性和數(shù)值模擬結(jié)果的驗證。

2.1 無量綱化

量綱定理是E.Buckingham[18]于1951年提出的。根據(jù)量綱π定理進行無量綱化。

第一步：寫出本問題涉及的所有物理量的量。

T為時間量綱，M為質(zhì)量量綱，L為長度量綱，表1所示為各物理量量綱中時間、質(zhì)量和長度量綱的指數(shù)，則t和Db可由ρ,D,σ線性表示。

表1 各物理量量綱中時間、質(zhì)量和長度量綱的指數(shù)

第二步，尋找剩余物理量對應(yīng)的無量綱量。

ln(t)=x1lnρ+x2lnD+x3lnσ

(5)

解出x1,x2,x3，即可得到時間t的無量綱量為

(6)

同理可得潤濕直徑的無量綱量

(7)

2.2 網(wǎng)格無關(guān)性驗證

圖2 無量綱化后不同網(wǎng)格數(shù)量計算的液滴潤濕壁面直徑隨時間的變化

圖2為無量綱化后不同網(wǎng)格數(shù)量計算的液滴潤濕壁面直徑隨時間的變化。比較了4種不同的網(wǎng)格，其計算結(jié)果變化趨勢比較相似，只是數(shù)值上略有不同，隨著網(wǎng)格數(shù)量的減少，其計算結(jié)果與網(wǎng)格數(shù)為240的網(wǎng)格的計算結(jié)果差距越大，而且考慮到隨著網(wǎng)格數(shù)量的增加，需要消耗更多的計算時間，綜合以上考慮，本文最終選用網(wǎng)格數(shù)為180的網(wǎng)格為計算網(wǎng)格。

2.3 數(shù)值模擬結(jié)果驗證

為了驗證所建立的液滴撞擊固體壁面模型的準確性，本文將數(shù)值模擬結(jié)果與實驗值進行了對比。實驗驗證基準采用的是H.Y.Kim等[19]的親水表面的液滴沉積實驗結(jié)果。

圖3為網(wǎng)格數(shù)為180時液滴撞擊壁面數(shù)值模擬結(jié)果與實驗結(jié)果對比圖?？梢砸姷?，在整個過程中數(shù)值模擬結(jié)果與實驗結(jié)果變化趨勢有很好的一致性，但有一定的誤差。因為數(shù)值模擬計算時，所有初始條件均可瞬時實現(xiàn)，而具體的實驗操作有實驗條件限制，并且實驗裝置的操作誤差和外界環(huán)境條件也是造成誤差的原因。數(shù)值模擬誤差在允許范圍之內(nèi)，通過對液滴撞擊壁面過程所建模型的驗證，說明所建模型設(shè)置合理，所得結(jié)果較為準確，因此可以模擬液滴撞壁鋪展的實際過程。

圖3 網(wǎng)格數(shù)為180時液滴撞擊壁面數(shù)值模擬結(jié)果與實驗結(jié)果對比

3 結(jié)果分析

液滴撞擊親水導(dǎo)流壁面的影響因素很多，本項目研究的是不同表面張力系數(shù)和不同液滴直徑對液滴撞擊親水壁面的影響。

3.1 不同表面張力系數(shù)的影響

要研究表面張力系數(shù)的影響，取液滴速度v=0.77 m/s，初始液滴直徑為D=3.6 mm，表面張力系數(shù)分別取σ=0.022 N/m、0.048 N/m、0.072 8 N/m。

圖4為不同表面張力系數(shù)的液滴撞壁后形狀隨時間的變化?？梢钥闯?，在撞擊開始階段，3種表面張力下的液滴形狀差別不大，到23 ms時，3種液滴都處于收縮狀態(tài)，液滴潤濕直徑大小順序為Db(σ=0.022 N/m)>Db(σ=0.048 N/m)>Db(σ=0.072 8 N/m)的液滴處于收縮狀態(tài)，另兩種液滴處于鋪展狀態(tài)，液滴潤濕直徑大小順序為Db(σ=0.022 N/m)

圖5為不同表面張力系數(shù)液滴撞擊壁面潤濕直徑隨時間的變化。可以看出，液滴撞擊導(dǎo)流面第一次鋪展過程中達到最大潤濕直徑的時間隨著表面張力系數(shù)的增大而減小，相應(yīng)的液滴對導(dǎo)流壁面的降溫效果也會更好。將液滴的一個運動周期(如0～26 ms)分成鋪展階段和收縮階段。因為有表面張力的存在，液滴在鋪展過程中會克服表面張力做功，在相同的速度下，液滴初始能量是相同的，表面張力越大，液滴運動會為克服表面張力消耗更多的能量。因此液滴在鋪展狀態(tài)時，液滴最大鋪展直徑越小，用的時間越少；在收縮階段，表面張力越大，液滴收縮越快，而且在鋪展階段的最大鋪展直徑越大，在收縮時會消耗更多的時間。液滴在導(dǎo)流壁面上鋪展的速度越快，到達高溫區(qū)域的速度越快，可以更高效地達到降溫效果，更好地保護導(dǎo)流壁面。

圖4 不同表面張力系數(shù)的液滴撞壁后形狀隨時間的變化

圖5 不同表面張力系數(shù)液滴撞擊壁面潤濕直徑Db隨時間t的變化

圖6為無量綱化后液滴撞擊壁面潤濕直徑隨時間的變化?？梢钥闯觯?條曲線有很好的同步性，說明無量綱化后3條曲線頻率一致。

3.2 不同液滴直徑的影響

要研究液滴直徑的影響，取液滴速度v=0.77 m/s，表面張力系數(shù)σ=0.072 8 N/mm，液滴直徑分別取3.6 mm、2 mm和1 mm。

圖7為不同直徑液滴撞擊壁面的潤濕直徑隨時間的變化?？梢钥闯?，隨著液滴直徑的增大，液滴在壁面上的潤濕直徑變化幅度增大，液滴一次鋪展收縮的時間也會相應(yīng)增加，相應(yīng)的液滴降溫效果會更好。3.6 mm的液滴一次鋪展收縮的時間約為26 ms，2 mm的液滴一次鋪展收縮的時間約為9.7 ms，1 mm的液滴一次鋪展收縮的時間約為4 ms。隨著時間增加，曲線趨于平緩，說明液滴在逐漸停止鋪展收縮，趨于靜止。因為液滴在壁面上鋪展達到的最大潤濕直徑越大，可以覆蓋更大的區(qū)域，相同數(shù)量的液滴可以鋪展更大的范圍，就可以保護更大范圍的導(dǎo)流壁面；而且液滴鋪展收縮的速度越慢，與導(dǎo)流壁面的接觸時間越長，足夠的相同時間內(nèi)也會使接觸面積更大，對導(dǎo)流壁面的保護作用更好。

圖6 無量綱化后不同表面張力系數(shù)液滴撞擊壁面潤濕直徑隨時間t*的變化

圖7 不同直徑液滴撞擊壁面的潤濕直徑Db隨時間t的變化

圖8為無量綱化后不同直徑液滴撞擊壁面的潤濕直徑隨時間的變化?？梢钥闯觯?條曲線有很好的同步性，說明無量綱化后3條曲線頻率一致。

圖8 無量綱化后不同直徑液滴撞擊壁面的潤濕直徑隨時間t*的變化

隨著液滴直徑的增加，液滴達到穩(wěn)定狀態(tài)就需要更多的時間，因為液滴直徑越大初始動能越大，液滴表面張力需要做更多的功才能抵消掉這部分動能；液滴直徑越大，鋪展收縮一次需要更多的時間，也是這個原理。圖7中曲線的每一次振蕩都比前一次振蕩幅值減小，而且振蕩的時間增大，是液滴克服表面張力做功，能量逐漸降低的結(jié)果。圖7中可以觀察到曲線的變化和彈簧系統(tǒng)的變化曲線非常相像，可以把液滴撞擊壁面過程類比彈簧系統(tǒng)，液滴有不同的初始速度相當于給彈簧加了不同大小的初始擾動，表面張力就相當于彈簧張力，可以類比彈簧的工作原理解釋一些現(xiàn)象。

4 結(jié)論

1) 通過與其他文獻中的實驗結(jié)果比較，證明CLSVOF方法可以用于液滴撞擊導(dǎo)流面模擬；

2) 液滴撞擊導(dǎo)流面形成液膜達到最大潤濕直徑的時間隨著表面張力系數(shù)的增大而減小，液滴對導(dǎo)流面的降溫效果更好；液滴撞擊在親水導(dǎo)流壁面上的最大潤濕直徑隨著液滴初始直徑的增大而增大，形成的液膜對導(dǎo)流壁面的保護作用更好；液滴的完成一次鋪展收縮過程需要的時間也隨液滴初始直徑的增大而增加，液滴對導(dǎo)流壁面的降溫效果也更好；

3) 無量綱化后，不同表面張力系數(shù)和不同初始直徑的液滴潤濕直徑隨時間的變化曲線有很好的同步性，頻率一致。

[1] WORTHINGTON A M.On the Forms Assumed by Drops of Liquids Falling Vertically on a Horizontal Plate[J].Royal Society of London Proceedings,1876,25(171/178):261-271.

[2] MAO T,KULMD C S,TRAN H.Spread and Rebound of Liquid Drops upon Impact on Flat Surface[J].ALCHE J,1997,43(9):2169-2179.

[3] SIKALO S,MARENGO M,TROPEA C,et al.Analysis of impact of droplets on horizontal surfaces[J].Experimental Thermal and Fluid Science,25(2002):503-510.

[4] SIKALO S,GANIC E N.Phenomena of Droplet-Surface Interactions[J].Experimental Thermal and Fluid Science,2006,31:97-110.

[5] WANG M J,HUNG Y L,LIN F H.Dynamic behaviors of droplet impact and spreading:a universal relationship study of dimensionless wetting diameter and droplet height[J].Experimental Thermal and Fluid Science,2009,33(7):1112-1118.

[6] RADULOVIC J,SEFIANE K,MARTIN E R，et al.On the effect of PH on spreading of surfactant solutions on hydrophobic surfaces[J].Journal of Colloid and Interface Science,2009,332(2):497-504.

[7] HAMRAOUI A,CACHILE M,POULARD C,et al.Fingering phenomena during spreading of surfactant solution[J].Colloids and Surfaces A Physicochemical and Engineering Aspects,2004,250(1-3):215-221.

[8] FISCHER B J,TROIAN S M.Growth and decay of localized disturbances on a surfactant-coated spreading film[J].Physical Review E,2003,6701(1):6309-6309.

[9] 畢菲菲,郭亞麗,沈勝強,等.液滴撞擊固體表面鋪展特性的實驗研究[J].物理學(xué)報,2012,61(18):295-300.

[10]李維仲,朱衛(wèi)英,權(quán)生林,等.液滴撞擊水平固體壁面的可視化實驗研究[J].熱科學(xué)與技術(shù),2008,7(2):155-160.

[11]PASANDIDEH-FARD M,AZIZ S D,CHANDRA S,et al.Cooling effectiveness of a water drop impinging on a hot surface[J].International Journal of Heat and Fluid Flow,2001,22:201-210.

[12]LUNKAD S F,BUWA V V,NIGAM K D P.Numerical simulation of drop impact and spreading on horizontal and inclined surfaces[J].Chemical Engineering Science,2007,62(24):7214-7224.

[13]GUNJAL P R,RANADE V V,CHAUDHARI R V.Dynamics of drop impact on solid surface:Experiments and VOF simulations[J].AIChE Journal,2005,51(1):59-78.

[14]EUNJEONG KIM,JEHYUN.Numerical study of the parameters governing the impact dynamics of yield-stress fluid droplets on a solid surface[J].Journal of Non-Newtonian Fluid Mechanics,2012(173):62-71.

[15]劉邱祖,寇子明,韓振南,等.基于格子Boltzmann方法的液滴沿固壁鋪展動態(tài)過程模擬[J].物理學(xué)報,2013,62(23):234701-234701.

[16]周龍玉,田瑞峰.液滴碰撞實驗與數(shù)值模擬研究[J].原子能科學(xué)技術(shù),2014,48(s1):314-319.

[17]陳石,王輝,沈勝強,等.液滴振蕩模型及與數(shù)值模擬的對比[J].物理學(xué)報,2013,62(20):312-317.

[18]BUCKINGHAM E.On Physically Similar Systems; Illustrations of the Use of Dimensional Equations[J].Physical Review,1951,4(4):345-376.

[19]KIM H Y,CHUN J H.The recoiling of liquid droplets upon collision with solid surfaces[J].PHYSICS OF FLUIDS,2001,13(3):643-659.

(責任編輯 楊繼森)

Numerical Simulation of Liquid Film Formed by a Water Droplet Impinging onto a Diversion Surface

MA Ying-jie，HAO Ji-guang

(School of Aerospace Engineering, Beijing Institute of Technology, Beijing 100081, China)

The theory analysis of the liquid film is formed by the droplet impinging on the hydrophilic diversion wall was carried out. With a theoretical model based on the CLSVOF method, the spreading characteristics of droplet impinging on a hydrophilic surface were analyzed. The effects on wetting diameter of different surface tension coefficients and different droplet diameters were researched. The results have shown that the CLSVOF method can be used for the simulation of a droplet impinging on a hydrophilic surface. The time of the droplet impinging on the diversion wall to reach to the maximum wetting diameter decreases as the increase of the surface tension coefficient of the droplet. The maximum wetting diameter of the droplet impinging on a hydrophilic surface increases as the increase of the initial droplet diameter. The time of the droplet spreading and recoiling increases as the increase of droplet diameters. It provides a theoretical basis for the further study of the protective effects of the liquid film formed by the droplet impinging on the diversion wall.

droplet impingement; CLSVOF; surface tension coefficient; droplet diameter

2016-10-20；

2016-11-20 基金項目：自然科學(xué)基金(51406012)

馬英杰(1991—)，男，碩士研究生,主要從事氣液兩相流傳熱傳質(zhì)研究。

郝繼光(1979—)，男，講師，碩士生導(dǎo)師。

10.11809/scbgxb2017.03.039

馬英杰,郝繼光.水液滴撞擊導(dǎo)流面成膜過程數(shù)值模擬[J].兵器裝備工程學(xué)報，2017(3):175-179.

format:MA Ying-jie，HAO Ji-guang.Numerical Simulation of Liquid Film Formed by a Water Droplet Impinging onto a Diversion Surface[J].Journal of Ordnance Equipment Engineering,2017(3):175-179.

TJ760.9

A

2096-2304(2017)03-0175-05