崔洪江, 張春利, 陳秉智

(大連交通大學(xué) 機(jī)車車輛工程學(xué)院， 遼寧大連 116028)

機(jī)車變流柜中的整流器將發(fā)電機(jī)組的交流電整流成直流電壓，逆變器將蓄電池和直流電壓逆變成交流電，以供給機(jī)車使用[1]。5 kW逆變單元散熱器和9 kW整流單元散熱器將工作部件熱量吸收通過風(fēng)道利用風(fēng)機(jī)以強(qiáng)迫對流的方式把熱量發(fā)散到機(jī)箱外，保證工作部件溫度正常[2-4]。

WANG等[5]應(yīng)用傳熱學(xué)理論對肋片散熱器進(jìn)行熱阻計(jì)算、仿真和試驗(yàn)，并進(jìn)行了對比研究證明。ZHANG等[6]建模仿真模擬發(fā)熱塊代替IGBT功率原件的方式對液冷板的流速、壓力分布和溫度分布等性能參數(shù)進(jìn)行試驗(yàn)驗(yàn)證。ZHANG等[7]運(yùn)用理論計(jì)算、數(shù)值模擬和試驗(yàn)研究對HXD1C電力機(jī)車輔助變流器進(jìn)行散熱設(shè)計(jì)。從WANG、ZHANG等文獻(xiàn)中可得出熱阻與流阻是散熱器的兩個重要參數(shù)，對散熱器進(jìn)行仿真計(jì)算時建立模擬發(fā)熱塊也是從熱阻與流阻這兩個角度進(jìn)行出發(fā)建立的。但是試驗(yàn)得到5 kW逆變單元散熱器熱阻R=0.005 4(m2·K/W)和9 kW整流單元散熱器熱阻R=0.008 6(m2·K/W)，可見散熱器熱阻影響小可忽略不計(jì)。所以阻尼在不改變散熱器耗熱值和不考慮熱阻的情況下，加載瞬態(tài)變化流阻曲線，在CFD計(jì)算時有利于誤差的減小。

在Fluent軟件建立變流柜模型時，用Ⅰ型和Ⅱ型阻尼替換機(jī)車變流柜關(guān)鍵部分5 kW逆變單元散熱器和9 kW整流單元散熱器，其中把阻尼耗熱值為5 kW定義為Ⅰ型阻尼及耗熱值為9 kW定義為Ⅱ型阻尼。由于實(shí)驗(yàn)得到兩組Ⅰ型、Ⅱ型阻尼阻力曲線，所以通過給Ⅰ型和Ⅱ型阻尼進(jìn)出口處設(shè)置監(jiān)測點(diǎn)，得到Ⅰ型和Ⅱ型阻尼進(jìn)出口處湍動流耗散率及熵產(chǎn)率，建立折線圖應(yīng)用流體力學(xué)及熱力學(xué)第二定律進(jìn)行對比分析，得到Ⅰ型阻尼和Ⅱ型阻尼合適的加載阻力曲線。

1 模 型

1.1 物理模型

圖1是變流柜物理模型的左視圖，由3部分組成：上部有4個百葉窗大小一樣對稱分布此為變流柜的進(jìn)口，中間部分有4個5 kW逆變單元散熱器和2個整流單元9 kW散熱器呈對稱分布，其中5 kW逆變單元散熱器和9 kW整流單元散熱器的具體結(jié)構(gòu)如圖2所示。下部為主風(fēng)道與風(fēng)機(jī)相連，此為變流柜的出口。

圖1 變流柜的物理模型

1.2 計(jì)算模型

圖3是變流柜計(jì)算模型的左視圖，與物理模型的主要區(qū)別在于中間部分用Ⅰ型、Ⅱ型阻尼替換了5 kW逆變單元散熱器和9 kW整流單元散熱器。

圖2 5 kW逆變單元散熱器和9 kW整流單元散熱器的具體結(jié)構(gòu)

圖3 變流柜的建模圖

2 邊界條件設(shè)定

該變流柜模型計(jì)算域的長×寬×高具體尺寸為2 700 mm×2 500 mm×3 000 mm，在出口處加載速度分別是35.16 m/s、39.07 m/s、42.98 m/s，其中對Ⅰ型阻尼和Ⅱ型阻尼所加載的阻力曲線如表1～表4所示。網(wǎng)格類型采用的是非結(jié)構(gòu)化網(wǎng)格，生成網(wǎng)格數(shù)目是200萬，計(jì)算采用κ-ε湍流兩方程進(jìn)行收斂計(jì)算。

表1 Ⅰ型阻尼加載第1組阻力曲線

表3 Ⅱ型阻尼加載第1組阻力曲線

表4 Ⅱ型阻尼加載第2組阻力曲線

3 流動性分析

3.1 理論分析

流體在阻尼內(nèi)流動過程中，分別滿足連續(xù)性方程、動量方程以及能量守恒方程。

連續(xù)性方程[8-10]如式(1)所示：

(1)

式中：ux，uy為流體速度在二維平面內(nèi)x,y方向上的分量；

x，y為與流體速度相對應(yīng)的坐標(biāo)軸方向；

ρ為流體密度。

動量方程[8-10]如式(2)所示：

(2)

式中：u為流動速度；p為流動過程中的壓力；S為流動過程中的廣義源項(xiàng)；i，j為不同的流動方向。

流體湍動能[8-10]計(jì)算公式如式(3)所示：

(3)

式中：k為湍動能；ε為湍流耗散率；αk為流體Prandtl數(shù)；CK為速度梯度引起的湍動能項(xiàng)；Sk為流體湍動能源項(xiàng)。

流體的湍流耗散率[8-10]計(jì)算公式如式(4)所示：

(4)

3.2 性能分析

對阻尼流道內(nèi)的流動情況進(jìn)行數(shù)值模擬，流動過程遵循連續(xù)性方程、動量守恒方程以及能量守恒方程。為了更好地分析阻尼內(nèi)部流體的流動情況提出如下假設(shè)：(1)流體連續(xù)介質(zhì)的假定；(2)流體的黏性可以忽略，把流體看成是無黏性的理想流體；(3)流體的膨脹性可忽略不計(jì)。

該變流柜計(jì)算模型總共有4個Ⅰ型阻尼和2個Ⅱ型阻尼，并且在其進(jìn)出口處設(shè)置12個監(jiān)測點(diǎn)，監(jiān)測點(diǎn)設(shè)置位置如圖4所示,設(shè)置的目的是為了監(jiān)測Ⅰ型和Ⅱ型阻尼在進(jìn)出口處的速度、壓力、湍動能和湍流耗散率。其中湍動能和湍流耗散率的數(shù)值越大越有助于阻尼散熱，這是此次仿真在流動性能方面上選用湍動能和湍流耗散率來做評價指標(biāo)的目的。根據(jù)邊界條件設(shè)定的相關(guān)參數(shù)，應(yīng)用Fluent軟件得到Ⅰ型和Ⅱ型阻尼計(jì)算結(jié)果，如圖5～圖10所示，橫軸表示監(jiān)測點(diǎn)位置，縱軸表示阻尼的湍動能和湍流耗散率數(shù)值大小。

由圖5、圖6可見，(1)Ⅰ型和Ⅱ型阻尼加載第1組阻力曲線所得湍動能在體積流量為4.62 m3/s監(jiān)測點(diǎn)為2的位置處最大是9.6 m2/s2，在體積流量為3.78 m3/s監(jiān)測點(diǎn)為1的位置處最小是2.65 m2/s2。(2)Ⅰ型和Ⅱ型阻尼加載第2組阻力曲線所得湍動能在體積流量為4.62 m3/s監(jiān)測點(diǎn)為2的位置處最大是9.97 m2/s2，在體積流量為3.78 m3/s監(jiān)測點(diǎn)為1的位置處最小是3.12 m2/s2。(3)在外部條件相同的情況下，隨著出口處加載的速度增大，在阻尼進(jìn)出口處的湍動能數(shù)值也隨著增大，原因是出口處加載的速度增大，在流經(jīng)阻尼的氣流中湍流的比例增多，進(jìn)而導(dǎo)致了阻尼進(jìn)出口處湍動能的數(shù)值增大。(4)監(jiān)測點(diǎn)1、2、3和4湍動能數(shù)值差異最大，原因是包含監(jiān)測點(diǎn)的阻尼距離下部主風(fēng)道的距離是最近的，氣體中湍流的比例最多導(dǎo)致的。

圖4 散熱器的監(jiān)測點(diǎn)

圖5 Ⅰ型和Ⅱ型阻尼加載第1組阻力曲線所得湍動能的折線圖

如圖7所示，在體積流量是4.20 m3/s條件下，Ⅰ型和Ⅱ型阻尼加載第二組阻力曲線所得湍動能數(shù)值大于Ⅰ型和Ⅱ型阻尼加載第一組阻力曲線所得湍動能數(shù)值，說明對于阻尼在耗熱值不變的情況下，隨著相同流速下靜壓力的的增大，會導(dǎo)致湍動能數(shù)值也隨著變小，這是符合實(shí)際情況的。因?yàn)樯崞魇峭ㄟ^系統(tǒng)中風(fēng)機(jī)的風(fēng)量來進(jìn)行散熱的，因此在流經(jīng)散熱器的風(fēng)量會有一定的阻力，阻力不能太大，太大會減小風(fēng)量通過散熱器，從而達(dá)不到散熱的要求。當(dāng)然阻力也不能太小，太小會使風(fēng)量流速加快不足以帶走足夠的熱量，也達(dá)不到預(yù)想的散熱的要求。由此可得Ⅰ型和Ⅱ型阻尼加載第二組阻力曲線是符合加載條件的。

圖6 Ⅰ型和Ⅱ型阻尼加載第2組阻力曲線所得湍動能的折線圖

圖7 在體積流量4.20 m3/s兩組阻力曲線的湍動能的比較

由圖8～圖9可見，(1)Ⅰ型和Ⅱ型阻尼加載第1組阻力曲線所得湍流耗散率在體積流量為4.62 m3/s，監(jiān)測點(diǎn)7處的湍流耗散率的數(shù)值最大是887 m2/s2，在體積流量為3.78 m3/s，監(jiān)測點(diǎn)1處的湍流耗散率的數(shù)值最小是116 m2/s2。(2)Ⅰ型和Ⅱ型阻尼加載第2組阻力曲線所得湍流耗散率在體積流量為4.62 m3/s，監(jiān)測點(diǎn)7處的湍流耗散率的數(shù)值最大是1 355 m2/s2，在體積流量為3.78 m3/s，監(jiān)測點(diǎn)1處的湍流耗散率的數(shù)值最小是114 m2/s2。(3)在外部條件相同的情況下，隨著體積流量的增大，阻尼進(jìn)出口處的湍流耗散率的數(shù)值也不斷的增大。(4)湍流耗散率的數(shù)值大小與湍動能的數(shù)值大小相比較，在湍動能數(shù)值大的阻尼，其湍流耗散率的數(shù)值小。

圖8 Ⅰ型和Ⅱ型阻尼加載第1組阻力曲線所得湍流耗散率的折線圖

圖9 Ⅰ型和Ⅱ型阻尼加載第2組阻力曲線所得湍流耗散率的折線圖

如圖10所示，在體積流量是4.20 m3/s情況下，Ⅰ型和Ⅱ型阻尼加載第2組阻力曲線所得湍流耗散率數(shù)值大于Ⅰ型和Ⅱ型阻尼加載第1組阻力曲線所得湍流耗散率數(shù)值，其中在監(jiān)測點(diǎn)5、7位置處，Ⅰ型和Ⅱ型阻尼加載第2組阻力曲線所得湍流耗散率數(shù)值明顯大于Ⅰ型和Ⅱ型阻尼加載第1組阻力曲線所得湍流耗散率數(shù)值，由此可得Ⅰ型和Ⅱ型阻尼加載第2組阻力曲線更有利于散熱器的散熱。

圖10 在體積流量4.20 m3/s兩組阻力曲線的湍流耗散率的比較

4 熵產(chǎn)分析

4.1 理論分析

熱力學(xué)第二定律(熵增加原理)其表述為：不可能把從低溫物體傳到高溫物體而不產(chǎn)生其他影響，或不可能從單一熱源取熱使之完全轉(zhuǎn)換為有用的功而不產(chǎn)生其他影響，或不可逆熱力過程中熵的微增量總是大于零。又稱“熵增定律”，表明了自然過程中，一個孤立系統(tǒng)的總混亂度(即“熵”)不會減小。

在熱力學(xué)中，熵是系統(tǒng)狀態(tài)參數(shù)，它的物理表達(dá)式為：

或dS=dQ/T

式中：S表示熵，Q表示熱量，T表示溫度。

(1)進(jìn)出口壓降

ΔP=P1-P2

式中：P1，P2分別為監(jiān)測進(jìn)出口截面監(jiān)測點(diǎn)的靜壓力，Pa；

ΔP為散熱器的風(fēng)阻，Pa。

(2)傳熱過程中的不可逆損失即熵產(chǎn)。考慮氣流及阻尼構(gòu)成的系統(tǒng)，在該系統(tǒng)中，氣流為低溫?zé)嵩?，熵變由下式?jì)算：

式中：T1，T2分別為監(jiān)測進(jìn)出口截面監(jiān)測點(diǎn)的溫度，K;qm為空氣質(zhì)量流量，kg/s;cp為空氣比定壓熱容，J/(kg·K);R為空氣的氣體常數(shù)，J/(kg·K);P1，P2分別為監(jiān)測進(jìn)出口截面監(jiān)測點(diǎn)的靜壓力，Pa。

4.2 性能分析

熵產(chǎn)最小法：從熱力學(xué)第一定律和第二定律的角度分析，自然界的一切實(shí)際過程都是不可逆的。散熱器作為一種散熱過程是典型的不可逆過程的裝置，熱量以導(dǎo)熱、對流換熱和輻射換熱3種方式進(jìn)行傳遞。散熱器中的不可逆損失分為有限溫差傳熱引起的不可逆損失和流體流動過程中克服摩擦阻力引起的不可逆損失，兩者之和即是散熱器中總的不可逆性損失。基于最小熵產(chǎn)原理，熵產(chǎn)越小，散熱器不可逆損失越小，其綜合性能越好[11-14]。

如圖11、圖12可見，(1)Ⅰ型和Ⅱ型阻尼加載第1組阻力曲線所得熵產(chǎn)率在體積流量是4.62 m3/s，阻尼6的熵產(chǎn)率最大值是233 J/(kg·K)。在體積流量是3.78 m3/s，阻尼1的熵產(chǎn)率的最小值是11 J/(kg·K)。(2)Ⅰ型和Ⅱ型阻尼加載第2組阻力曲線所得熵產(chǎn)率在體積流量是4.62 m3/s，阻尼6的熵產(chǎn)率最大值是91 J/(kg·K)，在體積流量是3.78 m3/s，阻尼2的熵產(chǎn)率的最小值是4 J/(kg·K)。(3)在Ⅰ型和Ⅱ型阻尼加載同一組阻力曲線，Ⅱ型阻尼熵產(chǎn)率明顯大于Ⅰ型阻尼熵產(chǎn)率，原因是Ⅰ型阻尼耗熱值低于Ⅱ型阻尼耗熱值，其中有限溫差傳熱引起的不可逆損失大于流體流動過程中克服摩擦阻力引起的不可逆損失，所以在熵產(chǎn)率中有限溫差傳熱引起的不可逆損失占主導(dǎo)地位，導(dǎo)致了Ⅱ型阻尼熵產(chǎn)率大于Ⅰ型阻尼熵產(chǎn)率，當(dāng)然也不可忽視流體流動過程中克服摩擦阻力引起的不可逆損失。

圖11 Ⅰ型和Ⅱ型阻尼加載第1組阻力曲線所得熵產(chǎn)率的折線圖

如圖13所示，在體積流量是4.20 m3/s條件下，Ⅰ型和Ⅱ型阻尼加載第2組阻力曲線所得熵產(chǎn)率數(shù)值小于Ⅰ型和Ⅱ型阻尼加載第1組阻力曲線所得熵產(chǎn)率數(shù)值，并且可見Ⅰ型和Ⅱ型阻尼加載第2組阻力曲線所得熵產(chǎn)率數(shù)值是Ⅰ型和Ⅱ型阻尼加載第1組阻力曲線所得熵產(chǎn)率數(shù)值1/2倍，符合熵產(chǎn)最小法。所以從熵產(chǎn)最小法的角度出發(fā)，Ⅰ型和Ⅱ型阻尼加載第2組阻力曲線更適合。

圖12 Ⅰ型和Ⅱ型阻尼加載第2組阻力曲線所得熵產(chǎn)率的折線圖

圖13 在體積流量4.20 m3/s兩組阻力曲線的熵產(chǎn)率的比較

5 結(jié) 論

采用阻尼的方式對機(jī)車變流柜中的散熱器進(jìn)行仿真計(jì)算，建立變流柜模型時中間部分用Ⅰ型和Ⅱ型阻尼替換了5 kW逆變單元散熱器和9 kW整流單元散熱器，應(yīng)用流體力學(xué)及熱力學(xué)第二定律，把湍動能、湍流耗散率和熵產(chǎn)率作為評價指標(biāo)，得到Ⅰ型阻尼和Ⅱ型阻尼合適的加載阻力曲線。

(1) 在對散熱器進(jìn)行熱分析時，由于阻尼在不改變散熱器耗熱值和不考慮熱阻的情況下，加載的流阻曲線是瞬態(tài)變化，這樣在CFD計(jì)算時有利于誤差的減小，并且應(yīng)用方便，可作為對散熱器進(jìn)行CFD計(jì)算的一種新的分析方法。

(2) 散熱器在相同尺寸大小的情況下，第1組阻力曲線所對應(yīng)散熱器的翅片密度大于第2組阻力曲線所對應(yīng)散熱器的翅片密度，但是第1組阻力曲線湍動能、湍流耗散率數(shù)值均小于第2組阻力曲線湍動能、湍流耗散率數(shù)值，并且第1組阻力曲線熵產(chǎn)率大于第2組阻力曲線熵產(chǎn)率也不符合最小熵產(chǎn)率，說明散熱器的翅片密度不能太大，否則適得其反。

(3) 把湍動能和湍流耗散率作為仿真結(jié)果處理的評價指標(biāo)，而不是以往速度和壓力，在3組不同體積流量加載的情況下，Ⅰ型和Ⅱ型阻尼加載第1、2組阻力曲線得到圖3～圖7計(jì)算結(jié)果，更加說明湍流的擾亂度大使得散熱器翅片邊界層厚度變小有利于其散熱。

(4) 把熵產(chǎn)率作為仿真結(jié)果處理的評價指標(biāo)，而不是溫度，在3組不同體積流量加載的情況下，Ⅰ型和Ⅱ型阻尼加載第1、2組阻力曲線得到圖11～圖13計(jì)算結(jié)果，更加說明不可逆熱力過程中熵的微增量總是大于零的。