電子膨脹閥射流噪聲數值模擬

黃 皓，王志毅，任夫磊，程德威

（1.浙江理工大學，a.機械與自動控制學院；b.建筑工程學院，杭州310018；2.奧克斯電氣有限公司，浙江寧波315191）

電子膨脹閥射流噪聲數值模擬

黃 皓1a，王志毅1b，任夫磊1a，程德威2

（1.浙江理工大學，a.機械與自動控制學院；b.建筑工程學院，杭州310018；2.奧克斯電氣有限公司，浙江寧波315191）

應用Realizable k-ε湍流模型和Mixture兩相流模型，對電子膨脹閥小開度下的二維流場進行了數值模擬，根據計算結果，使用寬頻噪聲計算方法得到了膨脹閥內流體噪聲分布。對電子膨脹閥的改進模型進行了數值模擬，并對流體的速度場及其噪聲輻射結果進行了分析，結果表明，閥口長度加長和閥口改為喇叭口結構可以減小射流噪聲。

電子膨脹閥；射流噪聲；FLUENT；數值模擬

0 引 言

電子膨脹閥是變頻空調中必不可少的節流元件，通過與壓縮機變容量技術的有效結合，可以有效提高變頻空調系統的舒適性與節能性［1-2］。然而，當系統在特定工況運行時，電子膨脹閥會產生較大的噪聲，有時甚至在室內都能聽到。目前，對于電子膨脹閥噪聲的研究很少，研究者主要是膨脹閥生產廠家以及空調制造商。珠海格力電器股份有限公司對冷媒流經電子膨脹閥時產生的節流噪聲進行了研究，指出電子膨脹閥內氣泡破裂噪聲一般發生在1 000～2 000 Hz的中頻區和5 000～6 000 Hz的高頻區，并針對各種噪音提出了相應的控制措施［3］。廣東美的制冷設備有限公司根據毛細管穩定性和準確性好的優點，在電子膨脹閥裝機時，在其前后增加一定長度的毛細管，降低了電子膨脹閥閥口上下游冷媒的壓力差，減小了冷媒的湍流強度，從而減小湍流噪聲［4］。但是這些改進措施均基于特定的空調系統進行，沒有提升電子膨脹閥自身的工作性能。電子膨脹閥產生的噪聲主要為流體噪聲，因此，從改進流道入手來解決電子膨脹閥的噪聲問題是企業掌握核心競爭力的必然選擇。

本文針對電子膨脹閥節流噪聲問題，通過FLUENT軟件開展流場和噪聲分布的數值模擬，分析了膨脹閥的流道幾何結構對流場和噪聲的影響，對電子膨脹閥降噪具有一定的工程指導意義。

1 數學模型

由于電子膨脹閥內流體存在相變過程，本文采用Mixture多相流模型，同時開啟空化模型。Mixture模型及空化模型的控制方程為：

a）連續性方程

b）動量方程

式中，ρ為流體介質密度，vi為介質在i方向上的速度，Xi為相應的空間i方向，Q為質量的產生率，Tij為廣義應力張量，Fi為其他源項。而是引起射流噪聲的主要原因，Tij包括三部分組成：

其中右端第一項為轉移張量，第二項為粘滯應力張量，第三項為熱傳導引起的應力張量。在射流中，力的變化成為偶極子輻射源，應力變化則成為四極子輻射源。在亞聲速和跨聲速的自由噴射流中，粘滯應力與熱傳導的影響可以忽略不計，故穩態噴射流中唯一的噪聲源來自于轉移動量的梯度張量，這種聲源一般為四極子聲源［5-6］。因此，可以使用FLUENT噪聲模型中的寬頻噪聲模型對電子膨脹閥射流噪聲進行數值分析［7］。

c）紊流動能k及其耗散率ε的輸運方程，

d）第二相的體積分數方程

關于流體噪聲的理論研究，Lighthill最早給出了自由射流中湍流噪聲總功率與噴射速度及出口直徑的關系式：

其中，W 為湍流總聲功率（W），ρ為流體密度（kg/m3），D為噴口孔徑（m），v為流體速度（m/s），ρ0和c0分別為聲傳播介質密度（kg/m3）和聲速（m/ s），K為Lighthill系數，實驗值約為0.3×10-4～1.8×10-4。對于自由射流以及固體邊界不起主要作用的噴射流，射流噪聲受噴射速度的影響最大，閥口噴射速度越大，產生的湍流噪聲也就越大。

Proudman［8］在Lighthill的噪聲理論基礎上，推導出一個各向同性湍流引起的聲功率公式。Lilley［9］通過計算Proudman的最初理論中所忽略的變化較小的時間差分項推導出另一公式，得到單位體積的聲功率為：

聲功率級計算公式為：

式中，u為湍流速度（m/s），l為特征長度（m），a0為當地聲速（m/s），α為一常數，Pref為參考聲功率（W）。

2 電子膨脹閥流場數值模擬

2.1 幾何建模及網格無關性分析

電子膨脹閥三維剖面圖如圖1所示，當流體從彎管流入、直管流出時為正向流動；當從直管流入、彎管流出時為反向流動，本文所有模擬均基于反向流動。為了縮短計算時間又能夠體現流場特征，將此模型進行了簡化，建立二維模型用于流場及噪聲分析［10］。

在開度為200脈沖工況下，針對網格無關性驗證，分別做了1萬、2萬和6萬3種網格數量的模型，如圖2所示。從圖2中可見，閥口附近網格經過了加密處理，最小尺寸均小于10-5m。通過在流道出口處設置一個質量流量監測點，得到換算為圓管道后的流量分別為0.033 0、0.033 9、0.035 4 kg/s，其相對誤差值均在5%以內，表明網格密度對計算結果的精度影響不大，本研究選擇2萬的計算網格。數值模擬所需的邊界條件設置為給定進、出口壓力條件，圖形左側顯示為黃色的邊界線設置為軸對稱邊界條件，其他均為無滑移壁面條件。

圖1 電子膨脹閥三維剖面

圖2 3種不同網格密度對比

2.2 流場模擬

數值模擬首先采用穩態求解器進行計算，然后以收斂的穩態解作為瞬態計算的初始值再次計算。計算過程采用的兩相流模型為Mixture模型、湍流模型為Realizable k-ε模型、空化模型為Schnerr-Sauer模型。設定邊界條件為進口壓力2.5 MPa，出口壓力0.8 MPa。流體物性參數為R410A兩相參數，2.5 MPa條件下對應的R410A飽和液相密度為967.5 kg/m3、動力粘度為9.72×10-5kg/m·s，飽和汽相密度為105.67 kg/m3、粘度為1.4×10-5kg/m·s，定義液相為主相，汽相為第二相。求解方法設定為SIMPLE算法，松弛因子均設置在0.3～0.5之間。

根據計算結果，得到流體速度分布如圖3所示。從圖3可以看出噴口處速度最大，在噴口下游閥針壁面帶拐角處速度梯度較其他地方大，此處湍流強度最高。從與其對應的湍流云圖圖4可以看出，此處即為射流的過渡區，其輻射出的噪聲強度亦為最大，達到46.5 dB，如圖5所示。將噪聲計算結果與文獻［11-12］對比后發現，模擬結果基本滿足國標要求，因此，可認為該計算模型適用于電子膨脹閥的噪聲模擬。

圖3 電子膨脹閥流體等速線

圖4 電子膨脹閥原始結構湍流動能分布

圖5 電子膨脹閥流體噪聲云圖

2.3 流道模型改進及對比分析

電子膨脹閥產生射流噪聲的主要原因是閥口前后壓差大，過渡區內流體速度梯度大、湍流強度高，而湍流強度又與噴管形狀和流體參數有關［13］。張勃等［14］通過實驗，研究了不同寬高比的圓轉矩形收斂噴管的射流寬、窄兩個對稱面上的湍流強度，分析結果表明，當寬高比大于8時，在射流中心線上距離噴口同一距離處，寬高比越大，湍流強度越大。這是因為噴口高度的減小，使得其對出口截面附近渦系發展的促進作用減弱。因此，本文嘗試通過加長上游閥口長度，來促進渦系發展以減弱湍流強度。改進前、后模型的參數對比如表1所示。

表1 不同結構模型參數對比

將原始結構電子膨脹閥的數學計算模型，應用到改進后的電子膨脹閥上，得到了改進后的電子膨脹閥的模擬結果。閥口長度由2 mm加長為3 mm的湍流動能分布如圖6所示，上游管道連接處由直角變為喇叭口的湍流動能分布（如圖7所示）。與圖4進行對比后發現，3種不同結構對應的混合區與過渡區的湍流分布均相似，改進后模型的湍流強度最大值均有不同程度的減弱，完全湍流區的湍流強度變大且占據更大的空間。當噴口加長后，出口處的渦系更容易控制在射流軸線范圍內并向下游傳播，射流結構中3個分區的連接更加緊密，湍流變化率趨于平緩，因此產生的湍流噪聲也相應減小，從圖8和圖9可以明顯發現這一變化。閥口加長后射流聲功率減小了0.3 dB，喇叭口結構對應的聲功率減小了1.2 dB，由此可以推斷，加長閥口長度和閥口改為喇叭口結構可以減小射流噪聲。

圖6 加長閥口結構湍流動能分布

圖7 喇叭口結構湍流動能分布

圖8 加長閥口結構聲功率分布

圖9 喇叭口結構聲功率分布

根據模擬的結果，得到3種模型噴口處的速度變化曲線，如圖10所示。從圖10中可以看出，兩種改進后的電子膨脹閥，出口速度均比原始結構大，這是因為上游流道的改變為流體流動提供了較長時間的動力，所以，流體動量增加，速度變大。根據Lighthill理論，膨脹閥出口速度越大，射流噪聲強度越大，噪聲聲壓級越大，這與模擬結果正好相反。產生這種現象的原因，是Lighthill理論只適用于自由射流和固體邊界不起主要作用的半受限射流，它無法對電子膨脹閥的射流噪聲進行解釋，這是由于電子膨脹閥的管道口徑較小，內部流動很容易受到壁面的影響。因此，在進行理論計算時，必須將邊界層的作用考慮在內。

圖10 不同結構閥口速度變化對比

3 結 論

本文采用Realizable k-ε模型和Mixture兩相流模型對小開度下3種不同結構的電子膨脹閥的流場進行了數值模擬，得到結論如下：

a）電子膨脹閥射流流場的湍動能分布與聲功率級分布極為相似，驗證了射流噪聲強度主要取決于速度梯度這一理論。

b）Lighthill理論無法解釋流場噪聲分布和閥口處的速度分布對比結果，這說明對自由射流和固體邊界不起主要作用的噴射流均適用的Lighthill理論不適用于口徑較小的電子膨脹閥噪聲計算。

c）通過寬頻噪聲模型計算得到的聲功率分布圖發現，加長閥口長度和閥口改為喇叭口結構的電子膨脹閥產生的射流噪聲均比原始結構要小，即閥口長度加長和閥口改為喇叭口結構可以減小射流噪聲。

［1］商萍君，易佳婷.電子膨脹閥的優勢和發展趨勢［J］.制冷與空調，2007，21（2）：72-76.

［2］閆長林.電子膨脹閥在家用變頻空調器中的應用［J］.家電科技，2013（3）：52-54.

［3］何 珍.淺談電子膨脹閥節流噪聲的產生與控制［J］.科技創業家，2013（4）：79-81.

［4］林 竹，吳 空.變頻電子膨脹閥空調冷媒流動異音分析研究［J］.家電科技，2012（11）：78-79.

［5］韓陽泉.噴注噪聲的數值模擬［D］.哈爾濱：哈爾濱工程大學，2007：45-48.

［6］馬大猷.現代聲學理論基礎［M］.北京：科學出版社，2004：94-95.

［7］安世亞太流體產品業務部.ANSYS Fluent航空氣動噪聲解決方案［J］.航空制造技術，2012（9）：103-104.

［8］Proudman I.The generation of noise by isotropic turbulence［J］.Proc Roy Soc，1952，A214：119.

［9］Lilley G M.The radiated noise from isotropic turbulence revisited［R］//NASA Contract Report.Hampton：NASA Langley Research Center，1993：75.

［10］王志剛，王冬青，柳成文，等.電子膨脹閥內制冷劑兩相流動過程的數值計算與實驗驗證［J］.流體機械，1999，27（8）：51-55.

［11］中華人民共和國國家經濟貿易委員會.JB/T 10386-2002，家用和類似用途空調電子膨脹閥［S］.2003.

［12］中華人民共和國工業和信息化部.JB/T 10212-2011，制冷空調用直動式電子膨脹閥［S］.2011.

［13］Shih T H，Chen J Y，Lumley J L.Second order modeling of boundary-free turbulent shear flows［J］.AIAA Journal，1992，30（6）：1553-1560.

［14］張 勃，吉洪湖，曹廣州，等.寬高比對矩形噴管射流湍流強度影響試驗［J］.航空動力學報，2010，25（10）：2244-2248.

NumericaI SimuIation of Jet Noise of EIectronic Expansion VaIve

HUANG Hao1a，WANG Zhi-yi1b，REN Fu-lei1a，CHENG De-wei2
（1.Zhejiang Sci-Tech University，a.School of Mechanical Engineering&Automation；b.School of Civil Engineering and Architecture，Hangzhou 310018，China；2.AUX Electrical Appliance Co.，Ltd.，Ningbo 315191，China）

Numerical simulation of two-dimensional flow field was carried out for electronic expansion valve with Realizable k-εturbulent model and mixture two-phase flow model.Based on the calculation result，flow noise distribution in the expansion valve was gained with broadband noise calculation method. Numerical simulation of the improved model of the electronic expansion valve was conducted.Besides，velocity field of the flow and noise radiation results were analyzed.The results show longer valve port and trumpet structure of the valve port can reduce jet noise.

electronic expansion valve；jet noise；FLUENT；numerical simulation

O358

A

（責任編輯：康 鋒）

1673-3851（2014）05-0512-05

2014-03-01

國家科技支撐計劃項目（2012BAF01B05）

黃 皓（1988-），男，河南信陽人，碩士研究生，研究方向為流體機械理論及應用。

王志毅，E-mail：zywang@zstu.edu.cn