潘保國，汪衛(wèi)華,，儲德林，王 鵬

計算機(jī)科學(xué)與自動化技術(shù)

基于橢圓檢測的液態(tài)PbLi與RAFM鋼的潤濕角的快速測算

潘保國1，汪衛(wèi)華1,2，儲德林2，王 鵬1

（1. 陸軍炮兵防空兵學(xué)院 基礎(chǔ)部，安徽 合肥 230031；2. 安徽大學(xué) 物質(zhì)科學(xué)與信息技術(shù)研究院，安徽 合肥 230039）

在惰性氣體保護(hù)下，利用高溫靜態(tài)座滴法測量并采用基于弧支撐線段的橢圓檢測方法計算液態(tài)PbLi在RAFM鋼表面的潤滑角，得到在255～425 ℃間不同溫度時液態(tài)PbLi在RAFM鋼表面的潤濕角均超過98°。由于不需要人工選擇邊界點，該橢圓檢測方法與橢圓擬合相比，速度更快，實現(xiàn)了實驗結(jié)果的快速測算。

霍夫變換；橢圓檢測；液態(tài)PbLi；RAFM鋼；潤濕角

在核聚變反應(yīng)堆中，包層是能量轉(zhuǎn)換的關(guān)鍵部件[1]，目前受到國際上普遍關(guān)注的、最有發(fā)展?jié)摿Φ木圩兌寻鼘邮且簯B(tài)金屬包層，而液態(tài)金屬與結(jié)構(gòu)材料之間的潤濕性能關(guān)系到熱交換的效率以及結(jié)構(gòu)材料的抗等離子體輻照，進(jìn)而決定反應(yīng)堆結(jié)構(gòu)材料的服役壽命。在材料科學(xué)中，金屬的液態(tài)與固態(tài)兩相之間的潤濕性能是一項基礎(chǔ)性的研究，屬于基本物理化學(xué)性質(zhì)，是反映固液界面潤濕性能大小的重要參數(shù)，對焊接、冶煉、鑄造、粉末冶金等過程都具有重要意義。

圖1 潤濕角示意圖

如圖1所示，當(dāng)液滴在固體基座表面上達(dá)到平衡狀態(tài)時，夾角即為潤濕角。通過對潤濕角的研究，可以獲得固、液、氣分子互相作用的諸多信息，如表面潤濕性、表面粗糙度、固液界面張力等[2-4]。

低活化馬氏體鋼（Reduced Activation Ferritic/ Martensitic Steel，RAFM鋼）是聚變反應(yīng)堆環(huán)境中優(yōu)先考慮的重要候選材料，本文采用高溫靜態(tài)座滴法測算液態(tài)PbLi[5]在RAFM鋼表面潤濕角。

目前潤濕角測量技術(shù)以實驗研究為主，主要分為角度測量法、力測量法、長度測量法和透射測量法，但以上方法都存在測量裝置復(fù)雜、人為操作因素影響較大等缺點。因此，近年的文獻(xiàn)較多地采用拍攝液滴在基座表面的靜態(tài)數(shù)字圖片，然后采用數(shù)字圖像處理技術(shù)測算出液滴在基座表面的潤濕角。王曉輝等人[6]采用了多項式擬合的方法計算潤濕角，多項式擬合中液滴輪廓曲線點數(shù)選取的多少，以及多項式階數(shù)的選取都會對最終計算結(jié)果產(chǎn)生較大的影響。徐志鈕等人[7]使用橢圓擬合方法可以更加有效地利用液滴邊界點信息，其擬合的結(jié)果要比多項式擬合的結(jié)果更加穩(wěn)定。但在使用橢圓擬合時必須依靠人工選取液滴的輪廓曲線并刪除其它無關(guān)的點線。這種采取人工干預(yù)的方法將會增加主觀因素造成的誤差。

本文首先將獲得的數(shù)字圖像進(jìn)行增強處理，使用霍夫變換定位底座表面所在的直線，計算出該直線方程。然后利用橢圓檢測的方法定位液滴輪廓所在的橢圓曲線，并計算出該橢圓曲線方程。最后依據(jù)底座表面直線方程和輪廓曲線的橢圓方程，計算出底座直線與輪廓曲線的交點，在交點處計算該橢圓切線的斜率，得到潤濕角。

1 實驗部分

1.1 實驗裝置和實驗方法

如圖2所示，將平板加熱盤①置于手套箱②內(nèi)放置平穩(wěn)，利用水平儀將其保持水平，連接好電源。將RAFM鋼板③置于平板加熱盤上中心位置，連接好溫度測試儀的溫度探頭。將PbLi連同包裝以及處于敞開狀態(tài)的石英罩置于手套箱內(nèi)。連接好真空泵抽氣管道④、保護(hù)氣體進(jìn)氣管道⑤。

圖2 實驗裝置圖

啟動真空泵抽真空，在手套箱內(nèi)氣體壓強減小到2 Pa左右時，充入高純氬氣（99.999%），并重復(fù)1次。借助手套箱的手套打開PbLi貯存裝置，取出PbLi顆粒⑥置于RAFM鋼上，然后關(guān)閉PbLi貯存裝置，罩上石英罩⑦。

打開平板加熱盤電源，對PbLi合金進(jìn)行加熱使其熔化，分別在255 ℃、280 ℃、306 ℃、332 ℃、357 ℃、380 ℃、395 ℃、411 ℃、425 ℃時，使用可變焦CCD相機(jī)近距離拍攝液滴圖片，并以RAW圖片格式保存。

1.2 潤濕角的計算

圖像在成像時，光學(xué)系統(tǒng)的失真、相對運動以及在傳輸過程中的噪聲污染等原因，都會造成圖像的某些降質(zhì)，這給分析測算潤濕角帶來困難，因此需要對圖像進(jìn)行處理，將圖像中感興趣的特征即液滴邊界進(jìn)行突出，而衰減其它信息。采取的步驟如圖3所示。

圖3 圖像處理及潤濕角測算步驟

(a)原始圖像(b)圖像銳化(c)灰度圖像(d)閾值轉(zhuǎn)換

中值濾波是一種非線性平滑技術(shù)，一般采用一個含有奇數(shù)個點的滑動窗口，將窗口中各點的灰度值的中值來代替指定點（一般是該窗口的中心點）的灰度值。以255 ℃時的PbLi液滴在RAFM鋼表面的潤濕角圖片（大小為138×78像素）為例，中值濾波處理圖像的關(guān)鍵步驟結(jié)果如圖4所示。

提取輪廓的算法就是掏空內(nèi)部點：如果原圖中有一點為黑（或白），且它的8個相鄰點都是同色時（即該點是內(nèi)部點），則將該點刪除。

輪廓提取的結(jié)果如圖5所示。

圖5 液滴輪廓曲線

霍夫變換是在一個參數(shù)空間中通過計算累計結(jié)果的局部最大值，得到一個符合某種幾何形狀的集合作為霍夫變換的結(jié)果。

圖6 基于霍夫變換的直線檢測結(jié)果

圖7 基于弧支撐線段的橢圓檢測結(jié)果

2 結(jié)果與討論

采用前述的實驗方法，分別在熔化的PbLi液滴被加熱至255 ℃、280 ℃、306 ℃、332 ℃、357 ℃、380 ℃、395 ℃、411 ℃、425 ℃時采集液滴圖片，并分別對這些圖像進(jìn)行處理、分析和測算，最后計算可得在上述各溫度時，基于霍夫變換和橢圓檢測計算后所得的潤濕角結(jié)果數(shù)據(jù)如表1所示。

表1 基于橢圓檢測計算所得潤濕角

由表1中的實驗結(jié)果可知在255 ℃～425 ℃下，液態(tài)PbLi合金與RAFM鋼的潤濕角均大于98°，潤濕性能較差。

在圖8中，小方框表示原始的實驗數(shù)據(jù)，虛線表示對實驗數(shù)據(jù)進(jìn)行直線擬合所得的結(jié)果。使用Matlab軟件中的corrcoef函數(shù)計算出潤濕角與溫度之間的相關(guān)系數(shù)為-0.446 3，負(fù)相關(guān)表示潤濕角隨著溫度的升高而減小，與圖8中擬合結(jié)果一致。這是因為，隨著環(huán)境溫度的升高，液相中的分子間作用力增大，表面張力減小，潤濕性增大，從而導(dǎo)致潤濕角減小[10-11]。

圖8 潤濕角與溫度的關(guān)系圖

3 結(jié)論

采用基于弧支撐線段的橢圓檢測方法定位液態(tài)PbLi合金在RAFM鋼表面的液滴輪廓，得到它在255～425 ℃下的潤濕角均大于98°。該方法不需要人為選取輪廓點，減少了由此產(chǎn)生的誤差，且具有重現(xiàn)性好、檢測速度快、自動化程度高的特點。

[1] 汪衛(wèi)華,陳俠.聚變堆液態(tài)包層MHD效應(yīng)的研究現(xiàn)狀與發(fā)展[J].艦船電子工程,2012,32(4):126-128.

[2] Wang X D, Tian Y, Peng X F. Self-aggregation of vapor- liquid phase transition[J]. Progress in Natural Science, 2003, 13(6): 451-456.

[3] Kwok D Y, Neumann A W. Contact angle measurement and contact angle interpretation[J]. Advances in Colloid and Interface Science, 1999, 81(3): 167-249.

[4] Emil C, Rafael P C. Problems of contact angle and solid surface energy determination[J]. Advances in Colloid and Interface Science, 2002, 98(2): 245-246.

[5] Yoshitaka U, Keiichi N, Tomoaki K, et al. Contact angle measurement of molten lead–lithium on silicon carbide surfaces[J]. Fusion Engineering and Design, 2011, 86(9- 11): 2297-2300.

[6] 王曉輝,李軍建,楊威,等.接觸角的圖像處理與檢測[J].光電子技術(shù),2011,31(1):14-19.

[7] 徐志鈕,律方成,趙麗娟,等.憎水性材料接觸角算法的選擇及在硅橡膠憎水性檢測中的應(yīng)用[J].高電壓技術(shù),2012, 38(4):814-823.

[8] Ballard D H. Generalizing the Hough transform to detect arbitrary shapes[J]. Pattern Recognition: The Journal of the Pattern Recognition Society, 1981, 13(2): 111-122.

[9] Lu C, Xia S, Shao M, et al. Arc-Support Line Segments Revisited: An Efficient High-Quality Ellipse Detection[J]. IEEE Transactions on Image Processing, 2020, 29: 768- 781.

[10] 張斌,李墨林,呂新華,等.溫度和合金元素對SiC與Al、AlSi7、AlSi7Mg1潤濕性的影響[J].大連鐵道學(xué)院學(xué)報, 1997,18(4):80-84.

[11] 錢健行.潤濕角測定裝置設(shè)計及其應(yīng)用研究[D].沈陽:東北大學(xué),2014:43-45.

Rapid Measurement of Wetting Angle of Molten PbLi on RAFM Steel Surfaces Based on Ellipse Detection

PAN Bao-guo1, WANG Wei-hua1,2, CHU De-lin2, WANG Peng1

(1. Department of Basic Course, Army Academy of Artillery and Air Defense, Hefei 230031, China; 2. Institutes of Physical Science and Information Technology, Anhui University, Hefei 230039, China)

Under the protection of inert gas, the wetting angle of liquid PBLI on the surface of RAFM steel was measured by the high-temperature static seat drop method and calculated by the ellipse detection method based on the arc support line segment. It was obtained that the wetting angle of liquid PBLI on the surface of rafm steel exceeded 98° at different temperatures from 255 ℃ to 425 ℃. Compared with ellipse fitting, this ellipse detection method can reduce the cost and error because there is no need to manually select the boundary points and the speed is faster. Based on this, the fast calculation of experimental results can be realized by programming.

Hough transform; ellipse detection; liquid PbLi; RAFM steel; wetting angle

TP391.41

A

1009-9115(2022)03-0059-04

10.3969/j.issn.1009-9115.2022.03.016

國家自然科學(xué)基金項目（51576208）

2021-08-11

2021-12-13

潘保國（1977-），男，安徽懷寧人，碩士，副教授，研究方向為科學(xué)計算與可視化。

（責(zé)任編輯、校對：田敬軍）