趙吉松

摘要：針對壁面摩擦力測量問題，建立了一種基于剪切敏感液晶（SSLC）涂層技術(shù)的平板表面摩擦力矢量場全局測量方法。該方法利用SSLC涂層在摩擦力作用下的顏色變化特性（不同方向顯示不同顏色）并結(jié)合其顏色變化與摩擦力大小之間的校準(zhǔn)關(guān)系解算摩擦力矢量的方向和大小，能夠測量整個待測區(qū)域的摩擦力矢量分布。應(yīng)用所述方法測量了平板表面凸起物繞流的摩擦力矢量場。試驗結(jié)果表明，該方法不僅能夠高分辨率測量出平板表面凸起物繞流的摩擦力矢量場，而且能夠研究凸起物的尾跡區(qū)隨著流動速度增加的發(fā)展過程。

關(guān)鍵詞：壁面摩擦力;測量;剪切敏感液晶;凸起物;干擾區(qū)

中圖分類號：O355 文獻(xiàn)標(biāo)識碼：A

流體與固體壁面之間的摩擦阻力是流體動力學(xué)領(lǐng)域的一個重要參數(shù)。在空氣動力學(xué)領(lǐng)域，許多重要信息可以通過顯示或者測量壁面摩擦阻力信息而獲取。流過飛行器表面的氣流產(chǎn)生的摩擦阻力會顯著影響飛行器的性能。內(nèi)部流動產(chǎn)生的摩擦阻力（如噴氣發(fā)動機(jī)壓縮空氣產(chǎn)生的摩擦阻力）同樣會對發(fā)動機(jī)的推阻性能產(chǎn)生重要影響。準(zhǔn)確測量摩擦阻力無論是在理論研究中還是工程實際中都具有重要意義。然而，壁面摩擦阻力的測量一直是一個難題，一直沒有很好的測量方法和技術(shù)。參考文獻(xiàn)[1]和參考文獻(xiàn)[2]對壁面摩擦力的主要測量方法進(jìn)行了綜述和評論。傳統(tǒng)的摩擦力測量方法主要采用機(jī)械或者電子的方法，如機(jī)械式天平、侵入式探針和傳感器等。這些方法屬于局部方法，只能測量單點摩擦力信息，并且會給流動帶來干擾，甚至?xí)o壁面帶來破壞。如果能夠高分辨率測量出壁面摩擦力的矢量分布，顯然具有重要的科學(xué)意義和工程價值。

剪切敏感液晶（Shear-Sensitive Liquid Crystal，SSLC）涂層技術(shù)是一種非接觸式壁面摩擦力矢量分布的全局測量方法。SSLC是一種固醇類液晶，其分子在物面上會形成一種螺旋軸垂直于壁面的螺旋結(jié)構(gòu)，螺旋軸的長度與可見光的波長量級相同。這種螺旋結(jié)構(gòu)在視覺上非?；钴S，對白光（光譜連續(xù)分布的光，如太陽光、鎢燈光等）選擇性反射，而且反射波長與螺旋軸長度成比例關(guān)系。在壁面摩擦力的作用下，螺旋結(jié)構(gòu)會變形，螺旋軸會傾斜，綜合效果是對人射光的反射具有很強(qiáng)的方向性，并且這種變化具有快速可逆性，隨摩擦力變化而變化。如果能將SSLC涂層的顏色隨這些參數(shù)變化的規(guī)律進(jìn)行校準(zhǔn)，那么便可以應(yīng)用SSLC涂層測量壁面摩擦力的矢量場?；谶@一思路，SSLC涂層技術(shù)已經(jīng)成功用于測量一些典型流動的壁面摩擦力矢量場[3～10]。除SSLC涂層技術(shù)，目前國際上提出的全局摩擦力場測量方法主要還有面應(yīng)力敏感薄膜技術(shù)[11，12]、微柱切應(yīng)力傳感器[13]、油膜干涉法[14]、熒光油膜法[15，16]、摩擦力診斷技術(shù)[17～19]以及軟基質(zhì)薄膜法[20]等。與這些方法相比，SSLC涂層技術(shù)的主要特色之處在于其對摩擦力的反應(yīng)是彩色的、視覺可見的，并且具有很高的時間分辨率（1kHz）和空間分辨率（像素級別）。

目前，根據(jù)SSLC涂層的顏色變化特性解算摩擦力矢量場的方法有兩種。一種本文稱之為多視角法，其原理是基于Reda等[21]的研究發(fā)現(xiàn)：在垂直光照射下，從不同方向觀測的SSLC涂層顏色的波長（或者色調(diào)）可以由Gauss曲線擬合，曲線的對稱軸為摩擦力方向，曲線峰值與摩擦力大小存在對應(yīng)關(guān)系。另一種本文稱之為兩視角法151。該方法采用兩個同步相機(jī)，借助于特定的試驗裝置提前校準(zhǔn)每個相機(jī)觀測的SSLC涂層的顏色隨摩擦力大小和方向的變化規(guī)律，在測量未知流場時能夠根據(jù)兩個相機(jī)采集的待測表面每個點的顏色信息在校準(zhǔn)曲線之間插值解算摩擦力大小和方向。兩視角法的優(yōu)勢是只需要從兩個方向觀測SSLC涂層顏色，但是其校準(zhǔn)過程比較復(fù)雜。

在測量精度方面，根據(jù)Reda等1211研究，兩視角方法的精度低于多視角方法。因為兩視角法只利用了兩個方向觀測SSLC涂層顏色，而多視角法利用5個或者更多方向觀測SSLC涂層顏色。盡管多視角法在原理上具有更高的測量精度，但是Reda等[3，4]在其研究中采用單個相機(jī)依次拍攝SSLC涂層不同方向的顏色信息。由于流動本身的非定常性，單個相機(jī)在從不同視角拍攝SSLC涂層的顏色時對應(yīng)的流動狀態(tài)并不完全一致，因而容易引入額外噪聲。Zhao等[7，8]將多視角法推廣至風(fēng)洞試驗測量，但是由于只采用單個相機(jī)拍攝不同方向的SSLC顏色，為了降低噪聲，需要多次測量取平均值。顯然，采用多個相機(jī)同步采集SSLC涂層不同方向的顏色可以解決這一問題[10]。此外，采用多臺相機(jī)同步采集SSLC涂層不同方向的顏色具有測量瞬時摩擦力場的潛力。

本文基于多臺同步相機(jī)和多視角方法建立了一種應(yīng)用SSLC涂層測量平板表面摩擦力矢量場的方法。應(yīng)用該方法測量了平板表面凸起物繞流的摩擦力矢量場，并且研究了凸起物尾跡區(qū)隨流動速度增加的發(fā)展過程。試驗結(jié)果展示了SSLC涂層技術(shù)測量壁面摩擦力場的能力。

1 試驗裝置

試驗研究在小型開口射流風(fēng)洞中開展，試驗裝置如圖1所示。該風(fēng)洞的出口為亞聲速噴管，噴管出口尺寸為寬4cm，高2.2cm。風(fēng)洞出口的氣流速度通過調(diào)整噴管壓比（Nozzle Pressure Ratio，NPR）進(jìn)行調(diào)節(jié)。試驗平臺固定于噴管出口處，平板表面與噴管下唇口平齊，平板尺寸為25cm×20cm。為了提高SSLC涂層顏色變化的對比度，在試驗平臺中間嵌入一個10cm×10cm的黑色電鍍鋁塊。電鍍鋁塊的前緣距離試驗平板前緣5cm。六面體凸起物（厚度0.5cm，寬度1.0cm）固定于試驗平臺，位于電鍍鋁塊的上游，距離其前緣1.0cm。

試驗中采用鹵鎢小燈泡（20W）提供法向照射光。燈泡置于距離測量區(qū)域正上方120cm處，以減小測量區(qū)域的光線照射方向的差異。該光源具有較好的光線平行度（見參考文獻(xiàn)[8]）。采用6臺相機(jī)（Canon EOS 80D）從不同的周向角φ同時拍攝SSLC涂層的顏色信息。相機(jī)俯視角（相機(jī)視線方向與測量表面之間的夾角）設(shè)置為28.5°。使用快門控制器控制相機(jī)同步拍照，測試結(jié)果表明相機(jī)的最高同步拍照速度最高可達(dá)250Hz。但是，由于照射光源的亮度不夠，本試驗中無法采用如此高的拍照速度（否則會曝光不足）。在本試驗中，6臺相機(jī)均設(shè)置為手動模式，其中光圈半徑F=11，曝光時間T=1/4s，感光度ISO=3200。

試驗采用的剪切敏感液晶為Hallerest BCN/192。本研究將液晶溶解于丙酮中，使用空氣刷均勻噴涂到待測表面。丙酮快速蒸發(fā)，留下一層紅色的 SSLC涂層。SSLC涂層厚度約為10μm（根據(jù)質(zhì)量守恒并考慮噴涂損失估算）。該液晶的清色溫度為49℃±1℃，即當(dāng)SSLC涂層溫度低于該溫度時，SSLC涂層對溫度不敏感，當(dāng)其溫度高于該溫度時，SSLC涂層變成無色。

2 摩擦力測量方法

本研究應(yīng)用Reda等[3]提出的多視角法解算摩擦力矢量場，不同之處是采用6臺相機(jī)同時拍攝SSLC涂層在不同方向的顏色變化。圖2給出摩擦力矢量場測量方法，共分為4步：

（1）在相機(jī)周向φ和摩擦力方φ一致的情況下，校準(zhǔn)SSLC涂層的顏色變化與摩擦力大小之間的關(guān)系。本試驗中，在測量區(qū)域的射流中心線投影線上剛好滿足φ=φ=0。用于顏色校準(zhǔn)的摩擦力大小，通過測量邊界層速度和采用修正Coles-Fernholz公式[22，23]計算得到。該計算公式已經(jīng)

（2）對于任意的待測流場，采用6臺相機(jī)同時記錄不同周向角觀測的SSLC涂層顏色。

（3）對于測量平面的每個點，采用Gauss曲線擬合hue-φ數(shù)據(jù)，擬合曲線的峰值對應(yīng)的周向角為摩擦力矢量的方經(jīng)過大量試驗數(shù)據(jù)的驗證。為了便于描述SSLC徐層的顏色，采用參考文獻(xiàn)[24]中的三色模型將相機(jī)拍攝的RGB信息轉(zhuǎn)換為色調(diào)（hue）信息。本試驗的照射光源和成像設(shè)備與參考文獻(xiàn)[10]相同，因此可以直接采用參考文獻(xiàn)[10]中的顏色校準(zhǔn)曲線，這里不再重復(fù)給出。向φ。

（4）將摩擦力方向?qū)?yīng)的hue值與步驟（1）得到的顏色校準(zhǔn)曲線相結(jié)合，解算出摩擦力矢量的大小。對于待測表面的每個點，重復(fù)步驟（3）和步驟（4）便可以得到全表面的摩擦力矢量場。

3 試驗結(jié)果與分析

圖3給出平板表面凸起物繞流試驗的一張原始照片（φ=17.7°）?？梢?，SSLC涂層通過綠色清晰地顯示出凸起物的尾跡區(qū)。

圖4給出相機(jī)從不同方向采集的SSLC涂層顏色。圖4中只給出了2號、4號和6號相機(jī)（相機(jī)編號參見圖2）采集的圖片，其他位于對稱位置的相機(jī)采集的圖片為鏡像圖片。圖4中的圖片已經(jīng)被變換成正視圖（采用測量區(qū)域周圍的正方形4個頂點作為校準(zhǔn)點），并且不需要部分已經(jīng)被裁剪。進(jìn)行這種視角變換是因為多視角方法要求用于Gauss曲線擬合的不同方向觀測的SSLC涂層顏色取自相同的物理點。每張圖片對應(yīng)的測量區(qū)域為9cm×9cm。圖片對稱軸與射流中心線在測量平面的投影重合?？梢姡琒SLC涂層在不同方向顯示不同的顏色。Reda等[3]研究發(fā)現(xiàn)，當(dāng)視線方向與摩擦力方向一致時，觀測到的SSLC涂層的顏色變化最大，并且摩擦力越大對應(yīng)的SSLC涂層的顏色變化越大。根據(jù)這一結(jié)論，可以定性分析不同區(qū)域的摩擦力信息（每張圖片的觀測方向是已知的）。如圖4所示的三張圖片中，圖4（a）中的SSLC涂層顏色變化較大（深綠色），因而摩擦力的主要方向與圖4（a）對應(yīng)的觀測方向較為接近（相對于其他兩張圖片而言）;凸起物尾跡區(qū)的SSLC涂層顏色變化較大，因而該區(qū)域的摩擦力較大。

為了降低圖片噪聲，同時保留顏色變化的梯度信息，對每張圖片的色調(diào)信息進(jìn)行中值濾波處理，濾波窗口取1.4mm×1.4mm。采用本文第2節(jié)的方法將不同周向角觀測的SSLC涂層顏色轉(zhuǎn)換為摩擦力矢量場。圖5給出在不同流動條件下（NPR=1.03～1.08）測量的平板凸起物繞流的摩擦力矢量場。為了便于顯示，圖中只在若干個v為常值的截面上，每隔1.2mm給出一個摩擦力矢量。圖中彩色云圖表示摩擦力大小，箭頭表示摩擦力方向，箭頭長短與摩擦力大小成正比。與參考文獻(xiàn)[5]、參考義獻(xiàn)[8]相比，圖5所示的摩擦力場噪聲較小，一方面因為多視角法具有較高的精度，另一方面因為本文采用多臺相機(jī)同步測量不同方向的SSLC涂層顏色，避免了由于流動不穩(wěn)定導(dǎo)致的噪聲。

需要說明的是，圖5給出的不同速度下的摩擦力場是由同一份SSLC涂層測量得到的?？梢?，同一份SSLC涂層可以用于測量不同來流速度下的摩擦力場。

觀察圖5可知，測量的摩擦力場基本對稱分布，這對于穩(wěn)態(tài)流動是比較合理的，并且也側(cè)面反映了SSLC涂層技術(shù)的可復(fù)現(xiàn)性（左右兩側(cè)的測量結(jié)果相當(dāng)于鏡像復(fù)現(xiàn)）。對比圖5（a）～圖5（d）所示的摩擦力矢量場可知，SSLC涂層技術(shù)成功捕獲了摩擦力場隨速度增加的發(fā)展過程。對于NPR=1.03到NPR =1.06，如圖5（a）～圖5（c）所示，凸起物的擾動效應(yīng)使得其尾跡區(qū)內(nèi)的摩擦力大小顯著增加（相對于其周圍區(qū)域），并且尾跡區(qū)的范圍以及尾跡區(qū)內(nèi)的摩擦力大小隨著來流速度增加而增加。但是，當(dāng)NPR進(jìn)一步增加至1.08時，凸起物尾跡區(qū)內(nèi)的摩擦力相對于周圍區(qū)域并沒有顯著增大，甚至有所減小，如圖5（d）所示。此時，凸起物的主要影響不再是擾動流場而是阻擋氣流并將氣流推向兩側(cè)。由此可見，SSLC涂層技術(shù)除了能夠高分辨率測量壁面摩擦力矢量場之外，還能夠提供流場發(fā)展過程的重要信息。

4 結(jié)論

本文基于多視角SSLC涂層技術(shù)建立了一種測量平板表面摩擦力矢量場的方法。該方法使用多臺同步相機(jī)從不同方向同時觀測SSLC涂層的顏色，與使用單臺相機(jī)相比，測量效率更高，能夠避免由于流動不穩(wěn)定導(dǎo)致的噪聲，并且具有測量瞬時摩擦力場的潛力。該方法的原理是采用Gauss曲線擬合不同方向的SSLC涂層顏色，解算摩擦力的方向和大小。應(yīng)用該方法測量了平板表面凸起物繞流的摩擦力矢量場，研究了凸起物尾跡區(qū)隨來流速度增加的變化過程。試驗結(jié)果表明，該方法高分辨率測量出平板凸起物繞流的摩擦力矢量場，并且捕獲了凸起物尾跡區(qū)隨流動速度增加的發(fā)展過程：當(dāng)NPR從1.03增加至1.06時，凸起物的擾動作用使其尾跡區(qū)的摩擦力增加;當(dāng)NPR進(jìn)一步增加至1.08時，凸起物的主要影響是阻擋氣流，其尾跡區(qū)的摩擦力相對周圍區(qū)域有所減小。在未來研究中，可提高照射光源的亮度，將所述方法推廣至測量瞬時摩擦力場。

參考文獻(xiàn)

[1]Naughton J W，Sheplak M.Modem developments in shear-stress measurement[J].Progress in Aerospace Sciences，2002，38（6）：515-570.

[2]Vinuesa R，orlii R.Measurement of wall shear stress[M].U.S.：CRC Press Taylor&Francis Group，2017.

[3]Reda D C，Wilder M C，F(xiàn)arina D J，et al.New methodology forthe measurement of surface shear stress vector distributions[J].AIAAJoumal，1997，35（4）：608-614.

[4]Reda D C，Wilder M C，Mehta R D，et al.Measurement ofcontinuous pressure and shear distributions using coating andimaging techniques[J].AIAA Journal，1998，36（6）：895-899.

[5]Fujisawa N，F(xiàn)unatani S，Kosaka S.Measurement of shear-stress distribution by liquid-crystal coating[J].MeasurementScience and Technology，2003（14）：1655-1661.

[6]Nakano T，F(xiàn)ujisawa N.Wind tunnel testing of shear-stressmeasurement by liquid-crystal coating[J].Journal ofVisualization，2006，9（2）：135-136.

[7]Zhao J，Scholtz P，Go L.Color change characteristics of twoshear-sensitive liquid crystal mixtures（BCN/192，BN/R50C）and their application in surface shear stress measurements[J]Science Bulletin，2011，56（27）：2897-2905.

[8]Zhao J，Scholtz P，Go L.Measurement of surface shear stressvector distribution using shear-sensitive liquid crystal coatings陰.Acta Mechanica Sinica，2012，28（5）：1261-1270.

[9]陳星，畢志獻(xiàn)，宮建，等.基于剪敏液晶涂層的光學(xué)摩阻測量技術(shù)研究[J].實驗流體力學(xué)，2012，26（6）：70-74.

[10]Zhao J.Measurement of wall shear stress in high speed airflow using shear-sensitive liquid crystal coating[J].Sensors，2018，18（5）：1605.

[11]Crafton J，F(xiàn)onov S，Jones E，et al.Optical measurements ofpressure and shear in a plasma[C]//35th AIAA FluidDynamics Conference and Exhibit，Toronto，Ontario Canada，2005.

[12]Fonov S，Jones G，Crafton J，et al.The development ofoptical techniques for the measurement of pressure and skinfriction[刀.Measurement Science&Technology，2006，16（6）：1-8.

[13]GroBe S，Schr6der W Mean wall-shear stress measurementsusing the micro-pillar shear-stress sensor MPS3[J].Measurement Science&Technology，2007，19（1）：015403.

[14]Vinuesa R，Bartrons E，Chiu D，et al.New insight into flowdevelopment and two dimensionality of turbulent channelflows[J].Experiments in Fluids，2014，55（6）：1759.

[15]Liu T，Montefort J，Woodiga S，et al.Global luminescent oil-film skin-friction meter[J].AUA Journal，2008，46（2）：476-485.

[16]Liu T，Woodiga S，Ma T.Skin friction topology in a regionenclosed by penetrable boundary[J].Experiments in Fluids，2011，51（6）：1549-1562.

[17]Woodiga S，Liu T.Feasibility of global skin frictiondiagnostics using temperature sensitive paint[J].Measurement Science&Technology，2011，22（11）：115402.

[18]Liu T，Woodiga S，Gregory J，et al.Global skin-frictiondiagnostics based on surface mass-transfer visualizations[J].AIAA Journal，2014，52（11）：2369-2383.

[19]Liu T，Misaka T，Asai K，et al.Feasibility of skin-frictiondiagnostics based on surface pressure gradient field[J].Measurement Science&Technology，2017，27（12）：125304.

[20]Bashirzadeh Y，Qian S，Manxthamuthu V Non-intrusivemeasurement of wall shear stress in flow channels[J].Sensors&Actuators a Physical，2018，271：118-123.

[21]Reda D C，Wilder M C，F(xiàn)arina D J，et al.Arealmeasurements of surface shear stress vector distributionsusing liquid crystal coatings[C]//34th Aerospace SciencesMeeting and Exhibit，Reno，NV，USA，1996.

[22]Nagib H M，Chauhan K A，Monkewitz P A.Approach to anasymptotic state for zero pressure gradient turbulentboundary layers[J].Philosophical Transactions aMathematical Physical Engineeging Sciences，2007，365（1852）：755-770.

[23]Samniguel Vila C，Vinuesa R，Discetti S，et al.Theidentification of well-behaved turbulent boundary layers[J].Journal of Fluid Mechanics，2017，822：109-138.

[24]Hay J L，Hollingsworth D K.A comparison of trichromiesystems for use in the calibration of polymer-dispersedthermochromic liquid crystals[J].Experimental Thermal&Fluid Science，1996，12（12）：1-12.