王銳，周致富，白飛龍，陳斌，王國祥,2

（1 西安交通大學(xué)動力工程多相流國家重點實驗室，陜西 西安 710049；2 美國Akron 大學(xué)機械工程系，Ohio Akron 44325）

引 言

葡萄酒色斑（port wine stain，PWS）是一種無法自愈的先天性血管畸形，表現(xiàn)為真皮淺層毛細血管數(shù)量過多，常出現(xiàn)在患者面部或頸部等暴露位置，給患者帶來容貌損害及心理障礙[1]。PWS 的傳統(tǒng)治療方法有同位素放射、冷凍、手術(shù)切除等[2]，但治療效果不佳。目前激光已成為治療PWS 首選方法，其治療原理是病灶血管中的血紅蛋白優(yōu)先吸收特定波長激光能量并發(fā)生不可逆熱破壞[3-4]。臨床廣泛使用的是波長585 nm 和595 nm 的脈沖染料激光，然而表皮中的黑色素對可見光波段的激光都有較強的吸收，易導(dǎo)致皮膚的爆皮、瘢痕、硬化、色素沉著等并發(fā)癥，且使激光攜帶的能量難以到達病灶血管，使得治療效率低、療效差[5-6]。

1995年，Nelson 等[7]提出了以R134a 為工質(zhì)的制冷劑噴霧冷卻（cryogen spray cooling，CSC）表皮的方式。在激光作用于皮膚之前，將飽和制冷劑液體噴射到皮膚表面，低沸點高揮發(fā)性的低溫制冷劑液滴沖擊皮膚表面，可將皮膚表面溫度迅速降低50℃以上。由于CSC 作用時間只有數(shù)十毫秒，可快速冷卻病變處的表皮組織，而對病變血管組織影響很小，從而達到空間選擇性的冷卻效果[8]。

當(dāng)前臨床使用的CSC 制冷劑為R134a，其常壓下沸點溫度為-26.07℃，提供的冷量有限，不足以對黑色素含量較高的深膚人種PWS 患者提供有效的表皮冷保護[9]。相較于R134a，R404a 具有更低的沸點（常壓下-46.5℃）和更強的揮發(fā)性，且具有環(huán)境友好性及對人體無害性，有望成為改善深膚色人種皮膚病患者治療效果的新型制冷劑[10]。本課題組在前期開展了R404a 制冷劑瞬態(tài)噴霧冷卻表面動態(tài)傳熱規(guī)律的研究[11]，并分析了噴霧距離與噴霧時間對表面瞬態(tài)溫度及熱通量的影響。制冷劑噴霧冷卻表面溫度及熱流表現(xiàn)出強烈的不均勻分布，但目前針對噴霧冷卻過程中表面?zhèn)鳠岵痪鶆蚍植继匦缘难芯枯^少[12]。本文在前期研究的基礎(chǔ)上，對R404a 瞬態(tài)噴霧冷卻表面?zhèn)鳠岬臅r間和空間不均勻性進行實驗研究，旨在對實際臨床激光手術(shù)中的噴霧冷卻提供指導(dǎo)。

1 實驗系統(tǒng)和方法

圖1 實驗系統(tǒng)示意Fig.1 Schematic of experiment system

圖2 噴嘴結(jié)構(gòu)Fig.2 Nozzle structure

表1 表皮與環(huán)氧樹脂的熱物性比較Table 1 Comparison of thermal properties of epidermis and epoxy resin

如圖1所示，制冷劑閃蒸瞬態(tài)噴霧冷卻的實驗系統(tǒng)由制冷劑儲罐、快速反應(yīng)電磁閥、電動平移臺、噴嘴、高壓連接軟管、冷卻基質(zhì)、熱電偶、高速攝 像儀及數(shù)據(jù)采集系統(tǒng)組成。制冷劑儲罐置于室溫環(huán)境，罐內(nèi)絕對壓力1.45 MPa。該壓力高于R404a 室溫條件下的飽和壓力，可有效防止R404a 在連接軟管中發(fā)生相變。噴射持續(xù)時間由控制快速反應(yīng)電磁閥的開閉時間決定，電磁閥反應(yīng)時間小于5 ms。噴嘴結(jié)構(gòu)如圖2所示，膨脹腔的長寬均為10 mm，直管出口段內(nèi)徑與長度分別為1 mm 及60 mm。采用這種帶膨脹腔的噴嘴，可使噴嘴出口制冷劑具有更好的霧化效果，減小制冷劑噴霧的擴散范圍，增強傳熱表面的冷卻效果[13]。冷卻基體采用幾何尺寸為50 mm×50 mm×5 mm 的環(huán)氧樹脂，其熱物性與人體表皮相近，如表1所示。

圖3 不同時刻R404a 噴霧對環(huán)氧樹脂表面的沖擊以及液膜殘留 Fig.3 Snapshots of droplet impacting on Plexiglas surface of R404a spray

為研究表面?zhèn)鳠岬臅r間和空間的不均勻性，通過電動平移臺（調(diào)節(jié)精度8 μm）調(diào)節(jié)噴嘴與冷卻表面的相對位置，使用熱電偶測量冷卻表面不同測點（測點與噴霧中心徑向距離分別為0、2、4、6、8、10 mm）處的瞬時溫度隨時間的變化。本文未采用常規(guī)的外部嵌入式熱電偶測溫，而是通過磁控濺射技術(shù)直接在環(huán)氧樹脂表面噴涂厚度約2 μm 的T 型薄膜熱電偶。薄膜熱電偶測溫范圍-200～350℃，溫度測量相對誤差小于1%，熱電偶的靜態(tài)溫度校核與動態(tài)響應(yīng)特性詳見本課題組之前的工作[14]。此方法具有測溫響應(yīng)時間快，且對被測物體表面的溫度場無干擾的優(yōu)點。為研究制冷劑噴霧形態(tài)對表面?zhèn)鳠岵痪鶆蛐缘挠绊?，實驗利用高速攝像儀（MotionXtra HG-100）對噴嘴外部噴霧形態(tài)進行可視化拍攝，并使用大功率LED 連續(xù)光源小角度仰視補光，拍攝頻率1000 fps，拍攝分辨率1504×1128。

由本課題組對R404a 瞬態(tài)噴霧冷卻表面動態(tài)傳熱規(guī)律的研究結(jié)果[11]可知，噴霧時間對表面動態(tài)傳熱影響很小，噴霧距離為30 mm 時噴霧表面最低溫度最小。所以本文研究R404a 瞬態(tài)噴霧冷卻傳熱空間不均勻性時，固定噴霧距離為30 mm、噴霧時間為50 ms，分別測量噴霧中心軸線及距離中軸線2、4、6、8 及10 mm 的徑向位置處的瞬態(tài)溫度。測得表面溫度后，應(yīng)用杜哈梅爾定理求解表面瞬態(tài)熱通量，其表達式為[17]

式中，q為基質(zhì)表面熱通量，ρ、c、λ分別為環(huán)氧樹脂基質(zhì)的密度、比熱容及熱導(dǎo)率，T為表面瞬時測量溫度。由于基質(zhì)表面溫度Ti是在離散的時間點ti（i=1，2，…，M）上進行測量的，假設(shè)溫度在每個相鄰時間段內(nèi)隨時間線性變化，積分式可展開為

研究傳熱表面溫度及熱通量的時間變化及徑向分布規(guī)律，即可獲得噴霧冷卻表面?zhèn)鳠岵痪鶆蛐蕴卣鳌?/p>

2 實驗結(jié)果與討論

2.1 R404a 噴霧的時間和空間的不均勻性

圖3示出了不同時刻R404a 瞬態(tài)噴霧沖擊環(huán)氧樹脂表面以及噴霧結(jié)束后液膜殘留圖片。從圖中可以看出，R404a 噴霧對冷卻表面的覆蓋面積隨時間強烈變化。因此，R404a 瞬態(tài)噴霧冷卻表面單點位置傳熱特性不僅與時間有關(guān)，而且與距噴霧中心的距離相關(guān)，噴霧冷卻表面?zhèn)鳠峋哂袝r間和空間的不均勻性。此外，在噴霧結(jié)束后，可以看到明顯的不均勻液膜殘留現(xiàn)象。

圖4 不同位置、不同時刻表面熱通量與溫度的變化曲線Fig.4 Surface heat flux as a function of surface temperature for different time and locations

圖4示出了表面熱通量（q）與溫度（T）的變化關(guān)系，圖4(a)自上而下各條實線分別表征距噴霧中心距離為0、2、4、6、8、10 mm 表面各位置處 的q-T關(guān)系，圖4(b)繪出30、60、90、120 ms 各時刻下的q-T關(guān)系，其中同位置不同時刻點通過虛線連接。圖4(a)中，相同噴霧半徑處，q-T曲線沿逆時針方向時間遞增；噴霧開始后，隨時間的變化，表面制冷劑液滴逐漸增多，T降低，q升高；至基質(zhì)表面剛好完全潤濕時，q達到峰值；T繼續(xù)下降后，表面制冷劑蒸發(fā)符合核態(tài)沸騰機制[18]，q-T曲線波動較大，過熱度的降低致使氣泡核減少，q隨之降低；表面制冷劑完全蒸發(fā)后，q逐漸減小至零，T上升。此外注意到，僅靠近噴霧中心（r=0、2 mm）時，表面溫度可達制冷劑沸點溫度Tboiling以下，噴霧邊緣區(qū)域溫度無法靠近沸點溫度。圖4(b)中，相同時間下，q-T曲線沿箭頭所示方向徑向距離遞增，T遞增，q基本呈下降趨勢。噴霧開始30 ms 及60 ms后，表面噴霧中心最低溫度可低至-40℃左右；90ms及120ms 時，最低溫度可低至Tboiling以下。由上述分析易知，不同時間、不同噴霧半徑處的熱通量與溫度變化曲線均不同，R404a 噴霧冷卻表面?zhèn)鳠峋哂袝r間和空間的不均勻性。

2.2 噴霧表面?zhèn)鳠犭S時間的變化

圖5展示了6 個噴霧半徑下R404a 瞬態(tài)噴霧冷卻環(huán)氧樹脂表面溫度T隨時間t的變化，溫度變化曲線自下而上分別對應(yīng)距噴霧中心距離為0、2、4、6、8、10 mm 表面各位置。環(huán)氧樹脂表面各位置處初始溫度均等于室溫，電磁閥開啟后開始計時，經(jīng)10 ms 左右的遲滯（包括電磁閥開啟時間、制冷劑在噴嘴內(nèi)流動時間及制冷劑噴霧液滴飛行30 mm距離的時間），噴霧液滴開始撞擊環(huán)氧樹脂，表面溫度均快速下降，50 ms 時噴霧結(jié)束。之后幾十至幾百毫秒內(nèi)表面溫度仍會繼續(xù)下降直至最小值Tmin，之后緩慢升溫直至室溫。從圖中可以看出，測量點越接近噴霧中心，溫度曲線下降速度越快，噴霧中 心附近區(qū)域表面溫度會降至R404a 沸點溫度以下，這是由環(huán)氧樹脂表面噴霧中心區(qū)域大量制冷劑液滴相變換熱引起的；而噴霧外緣區(qū)域，如本實驗中的半徑10 mm 處，表面制冷劑液滴量較少，溫度曲線變化和緩，最低溫度與室溫相近。

圖5 不同噴霧半徑處表面溫度隨時間變化Fig.5 Time variation for surface temperature

圖6 不同噴霧半徑處表面熱通量隨時間變化Fig.6 Time variation for surface heat flux

對應(yīng)地，圖6展示了6 個噴霧半徑下表面熱通量q隨時間t的變化，熱通量變化曲線自上而下分別對應(yīng)6 組半徑0～10 mm。與表面溫度一樣，表面熱通量也隨時間強烈變化。當(dāng)制冷劑液滴撞擊表面時，q快速上升至峰值qmax；之后制冷劑液滴在環(huán)氧樹脂表面形成液膜，q變化較為緩慢；當(dāng)液膜完全消失時，q下降至0。從圖中可以看出，測量點越接近噴霧中心，熱通量上升速度越快，且達到的熱通量峰值越大；而噴霧外緣區(qū)域，如本實驗中的半徑10 mm 處，表面制冷劑液滴量較少，表面熱通量值較小。

為更好顯示表面溫度及傳熱的不均勻性，定義表面不同區(qū)域與噴霧中心位置的溫度差 ( ,)T r tΔ=T(r,t) -T(0,t)及熱通量差 Δq(r,t)=q(0,t) -q(r,t)，ΔT與Δq隨時間變化關(guān)系示于圖7。如圖7(a)所示，各位置處ΔT在噴霧10 ms 后急速上升，均于15 ms左右時達到峰值，此時空間溫度分布不均勻性達到最大。半徑2 mm 位置處的ΔT取值始終接近零，說明噴霧中心區(qū)域溫度分布較均勻；相同噴霧時間下，半徑越大，ΔT值越大；噴霧邊緣的半徑10 mm處，ΔT最大可達60℃左右。圖7(b)所示的表面熱通量差值最大可達400 kW·m-2，此外Δq有負值出現(xiàn)，說明中心區(qū)域的熱通量不是始終高于周圍 區(qū)域。

2.3 噴霧表面?zhèn)鳠嵫貜较虻淖兓?/h3>

圖7 不同噴霧半徑與中心區(qū)域處溫度差及 熱通量差隨時間變化Fig.7 Time variation for differences in temperature and heat flux between spray center and other locations

圖8示出了表面溫度及熱通量在不同時刻的徑向分布。初始時刻，各位置溫度為室溫，熱通量為 零。噴霧開始25 ms 時，表面不同位置處以不同的熱通量進行換熱，溫度分布呈現(xiàn)不均勻性。50 ms時噴霧結(jié)束，表面不同位置仍維持較大熱通量，溫度繼續(xù)降低，溫度徑向分布曲線近似成高斯分布，與文獻[12,19]結(jié)果一致。之后一段時間內(nèi)，溫度曲線依然保持相似的形狀。此外特別注意到，中心區(qū)域與半徑2 mm 處的溫度及熱通量取值始終保持相近的值，說明噴霧表面中心存在徑向分布均勻的小區(qū)域。在本實驗條件下，中心區(qū)域半徑2 mm 范圍內(nèi)表面?zhèn)鳠岬目臻g分布均勻，50～75 ms 時間段內(nèi)，半徑4 mm 范圍內(nèi)表面?zhèn)鳠岬目臻g分布均勻。

圖8 表面溫度及熱通量在特定時間的徑向分布Fig.8 Radial distribution for surface temperature and heat flux

為表征噴霧冷卻表面?zhèn)鳠岬奶匦?，將表面溫度極值Tmin和表面熱通量峰值qmax及其發(fā)生時間的徑向分布示于圖9。由圖可知，Tmin隨半徑增加而遞增，qmax隨半徑增加而遞減，說明噴霧表面冷卻效果隨半徑增加逐漸降低。Tmin和qmax的取值在半徑8～10 mm 處發(fā)生突變，說明此處為噴霧邊緣區(qū)域；8 mm 半徑處最低溫度-32.7℃，仍符合臨床上激光治療對冷卻溫度的要求，在本實驗條件下符合冷卻 溫度要求的表面有效區(qū)域半徑即為8 mm 左右。制冷劑噴霧液滴溫度在噴嘴出口處已接近大氣壓下對應(yīng)的沸點溫度（-46.5℃），液滴繼續(xù)飛行過程中，由于其高揮發(fā)性，周圍對流換熱不足以提供液滴快速蒸發(fā)需要的潛熱，還需要從其自身剩余部分吸熱，因而導(dǎo)致其溫度不斷降低（可遠低于沸點溫度）[11]。低溫液滴沖擊基質(zhì)表面后繼續(xù)蒸發(fā)，部分殘留并形成液膜，通過導(dǎo)熱、對流與沸騰多種方式與基質(zhì)換熱可快速降低基質(zhì)表面溫度，在液滴沖擊較強的冷卻中心區(qū)域可使得基質(zhì)溫度降低至工質(zhì)沸點溫度以下[20]。噴霧中心區(qū)域附近換熱強烈，20 ms 內(nèi)熱通量迅速升高并達到峰值；隨半徑增加，表面換熱速度減慢，需較長時間達到熱通量峰值?；|(zhì)表面殘留液膜可維持表面溫度持續(xù)降低，經(jīng)較長時間（約250 ms）達到溫度最低值；隨半徑增加，表面溫度無法達到沸點溫度以下，表面液滴迅速蒸發(fā)而無法堆積形成液膜，遠離中心區(qū)域處的表面液滴在極短時間內(nèi)完全蒸發(fā)，幾十毫秒內(nèi)即達到溫度最低值，之后逐漸升溫。臨床上治療持續(xù)時間通常為150 ms左右，在本實驗條件下符合冷卻持續(xù)時間要求的表 面有效區(qū)域半徑為2 mm 左右。

圖9 Tmin 和qmax 的峰值及發(fā)生時間的徑向分布Fig.9 Radial distribution for magnitude and occurrence time of Tmin and qmax

表面溫度及熱通量等值曲線示于圖10。由等溫線圖可知，表面溫度于較大噴霧時間及較小的噴霧半徑處取得低值；噴霧時間20 ms 之后，溫度沿徑向變化顯著，僅噴霧中心區(qū)域附近達到低溫；在噴霧中心附近區(qū)域，溫度隨時間變化不顯著，且低溫值均符合激光手術(shù)治療PWS 的要求。由等熱通量曲線圖可知，高熱通量值僅發(fā)生在幾十毫秒時的靠近噴霧中心附近處，這說明CSC 是時間尺度10-3s和空間尺度10-3m 的瞬時小尺度過程。結(jié)合圖3可以看出，表面?zhèn)鳠岬臅r空不均勻性與噴霧形態(tài)密切相關(guān)。噴霧開始后較短時間內(nèi)，如圖3中的5 ms和15 ms，噴霧錐角較小，表面制冷劑覆蓋范圍小，僅中心區(qū)域附近溫度降至-30℃以下低溫，大量制冷劑集中于中心區(qū)域使表面熱通量值較高；噴霧穩(wěn)定狀態(tài)時，如圖3中的25 ms，大噴霧錐角使表面制冷劑覆蓋范圍增大，中心半徑6 mm 區(qū)域內(nèi)表面溫度均可降至-30℃以下低溫，制冷劑分散較開致使中心區(qū)域表面熱通量明顯降低；噴霧結(jié)束后一段時間，表面殘留的制冷劑液膜使表面持續(xù)維持低溫，表面熱通量逐漸降低至趨于零。

圖10 表面溫度及傳熱的時間空間分布Fig.10 Time and space dependent surface temperature and heat flux profiles

3 結(jié) 論

R404a 制冷劑閃蒸瞬態(tài)噴霧的流動與傳熱過程非常復(fù)雜，從實驗結(jié)果發(fā)現(xiàn)冷卻基質(zhì)的表面?zhèn)鳠峋哂袝r間和空間的不均勻性。噴霧中心區(qū)域附近溫度響應(yīng)快，噴射開始十幾毫秒內(nèi)即可達到熱通量峰值，最低溫度值可達到制冷劑沸點溫度以下，且能在幾百毫秒內(nèi)維持低溫。而靠近邊緣區(qū)域的表面溫度高于中心、熱通量低于中心區(qū)域，制冷劑噴霧冷卻對表皮的冷保護作用不均勻，噴霧中心與噴霧邊緣的表皮溫度差最大可達60℃左右，表面熱通量差值最大可達400 kW·m-2。高熱通量值僅發(fā)生在幾十毫秒時的噴霧中心附近區(qū)域，說明CSC 是時間尺度10-3s 和空間尺度10-3m 的瞬時小尺度過程。在噴霧中心附近區(qū)域存在表面?zhèn)鳠峋鶆蚍植嫉男^(qū)域，

在本文的實驗條件下，噴霧中心處與半徑2 mm 位置處的溫度及熱通量幾乎始終保持一致，可認為半徑2 mm 的中心區(qū)域內(nèi)R404a 噴霧冷卻表面?zhèn)鳠峋鶆蚍植?。在實際臨床激光手術(shù)治療PWS 過程中，可參考本實驗結(jié)果保證CSC 對治療范圍內(nèi)的表皮進行均勻冷保護。

[1]Lin Xiaoxi (林曉曦),Qi Zuoliang (祁佐良),Dong Jiasheng (董佳生).Research on the mechanism of enlargement of portwine stains [J].Journal of Practical Aesthetic and Plastic Surgery(實用美容整形外科),2000,11 (3):126-129

[2]Kim B S,Lee J B,Jang H S,Kwon Y W,Oh C K.Multiple basal cell carcinomas arising in a port-wine stain with a remote history of therapeutic irradiation [J].J.Dermatol.,2004,31 (10):820-823

[3]Kercher,D M,Sheer R E,So R M C.Short duration heat transfer studies at high free-stream temperatures [J].Journal of Engineering for Gas Turbines and Power,1983,105 (1):156-166

[4]Kelly K M,Nelson J S.Update on the clinical management of port wine stains [J].Lasers in Medical Science,2000,15 (4):220-226

[5]Karapetian E,Aguilar G,Kimel S,Lavernia E J,Nelson J S.Effects of mass flow rate and droplet velocity on surface heat flux during cryogen spray cooling [J].Physics in Medicine and Biolog,2003,48:N1-N6

[6]Majaron B,Nelson J S.Optical-Thermal Response of Laser-Irradiated Tissue[M].Netherlands:Springer,2011:859-913

[7]Nelson J S,Milner T E,Anvari B,Tanenbaum B S,Kimel S,Svaasand L O,Jacques,S L.Dynamic epidermal cooling during pulsed laser treatment of port-wine stain:a new methodology with preliminary clinical evaluation [J].Arch.Dermatol.,1995,131 (6):695-700

[8]Chang C J,Anvari B,Nelson S J.Cryogen spray cooling for spatially selective photocoagulation of hemangiomas:a new methodology with preliminary clinical reports [J].Plastic and Reconstructive Surgery,1998,102 (2):459-463

[9]Tunnell J W,Chang D W,Johnstone C,Torres J H,Patrick C W,Miller M J,Anvari B.Effects of cryogen spray cooling and high radiant exposures on selective vascular injury during laser irradiation of human skin [J].Arch.Dermatol.,2003,139 (6):743-750

[10]Dai T,Yaseen M A,Diagaradjane P,Yaseen M A,Diagaradjane P.Comparative study of cryogen spray cooling with R-134a and R-404a:implications for laser treatment of dark human skin [J].J.Biomed.Opt.,2006,11 (4):041116-041127

[11]Zhou Z F,Chen B,Wang Y S,Guo L J,Wang G X.An experimental study on pulsed spray cooling with refrigerant R-404a in laser surgery [J].Applied Thermal Engineering,2012,39:29-36

[12]Franco W,Liu J,Wang G X.Radial and temporal variations in surface heat transfer during cryogen spray cooling [J].Phys.Med.Biol.,2005,50 (2):387-397

[13]Sher E,Bar-Kohany T,Rashkovan A.Flash-boiling atomization [J].Progress in Energy and Combustion Science,2008,34 (4):417-439

[14]Zhou Zhifu (周致富),Wu Weitao (吳威濤),Wang Guoxiang (王國祥),Chen Bin (陳斌),Wang Yueshe (王躍社) ,Gong Zheng (龔政).Thin-film thermocouples for rapid measurement of transient surface temperature in cryogen spray cooling [J].CIESC Journal(化工學(xué)報),2011,62 (11):3060-3064

[15]Duck F A.Physical Properties of Tissues:A Comprehensive Book [M].London:Acadamic Press,1990

[16]Pikkula B M,Torres J H,Chang D W,Lavernia E J,Nelson J S.Effects of mass flow rate and droplet velocity on surface heat flux during cryogen spray cooling [J].Phys.Med.Biol.,2001,28 (2):113-120

[17]Taler J.Theory of transient experimental techniques for surface heat thansfer [J].Heat Mass Transfer,1996,17:3733-3748

[18]Aguilar G,Wang G X,Nelson J S.Effect of spurt duration on the heat transfer dynamics during cryogen spray cooling [J].Phys.Med.Biol.,2003,48:2168-2181

[19]Choi B,Welch A J.Infrared imaging of 2-d temperature distribution during cryogen spray cooling Trans [J].ASME J.Biomech.Eng.,2002,124:669-675

[20]Zhou Z F,Wu W T,Chen B,Wang G X,Guo L J.An experimental study on the spray and thermal characteristics of R134a two-phase flashing spray [J].International Journal of Heat and Mass Transfer,2012,55 (15/16):4460-4468