大氣濕度對柴油機缸內(nèi)工質(zhì)熱物性參數(shù)的影響

倪培永，王向麗

(1.南通大學(xué) 機械工程學(xué)院，江蘇 南通 226019；2.南通大學(xué) 電氣工程學(xué)院，江蘇 南通 226019)

大氣濕度對柴油機缸內(nèi)工質(zhì)熱物性參數(shù)的影響

倪培永1，王向麗2

(1.南通大學(xué) 機械工程學(xué)院，江蘇 南通 226019；2.南通大學(xué) 電氣工程學(xué)院，江蘇 南通 226019)

將有關(guān)熱物性參數(shù)計算表達式嵌入KIVA-3軟件中，這些參數(shù)包括壓力、溫度、比熱容、比熱比、比焓、導(dǎo)熱系數(shù)、動力粘度、音速、壓縮系數(shù)和普朗特數(shù)。通過改變進氣成分，計算進氣含濕量分別為0，5，10和20 g/kg時柴油機燃燒過程中缸內(nèi)工質(zhì)熱物性參數(shù)的變化規(guī)律。計算結(jié)果表明，大氣濕度對燃燒溫度、比熱容和普朗特數(shù)影響較大，對壓力、比熱比、比焓和音速影響較小，對壓縮系數(shù)影響很小。含濕量增加，在高溫時動力粘度和導(dǎo)熱系數(shù)兩者稍微增加。

柴油機；含濕量；熱物性參數(shù)；模擬

0 引 言

由于天氣、季節(jié)等客觀因素會造成大氣濕度的變化。已有研究結(jié)果表明，大氣濕度尤其是高濕度對柴油機排放有重要的影響[1-4]。 這說明濕空氣會對燃燒過程產(chǎn)生一定的影響。 若能研究濕空氣在柴油機燃燒過程中熱物性參數(shù)變化規(guī)律，一方面可獲得濕度和熱物性參數(shù)之間的定量關(guān)系，另一方面可加深了解濕空氣在燃燒過程中所起的物理作用，對燃燒現(xiàn)象給以更詳盡的解釋。

目前，濕空氣熱物性參數(shù)的研究多集中于濕空氣透平(HAT)循環(huán)方面，具體的研究包括：計算高濕空氣在高溫高壓區(qū)域內(nèi)的熱物性參數(shù)[5-7]，介紹熱物性參數(shù)的計算方法[8]，研究濕空氣擴散燃燒特性[9-10]。然而,對柴油機濕空氣燃燒過程中熱物性參數(shù)研究報道并不多?；诖?，本文采用KIVA-3軟件模擬不同含濕量時柴油機燃燒過程熱物性參數(shù)規(guī)律。研究結(jié)果可為應(yīng)用進氣加濕方法改善柴油機氮氧排放提供研究依據(jù)。

1 計算模型的建立

1.1 幾何模型

計算用柴油機主要技術(shù)參數(shù)如表1所示。燃燒室結(jié)構(gòu)已經(jīng)簡化，簡化后的計算網(wǎng)格結(jié)構(gòu)如圖1所示。柴油機噴油器噴孔數(shù)為5，為節(jié)省計算時間，故選取燃燒室的1/5作為計算區(qū)域。

大氣濕度有2種表示方式：一種是絕對濕度，即含濕量，另一種是相對濕度，兩者可相互轉(zhuǎn)換。本文采用含濕量來定義大氣濕度，其定義為1 kg干空氣同時并存的水蒸氣量，單位是g/kg。本文以含濕量為0，5，10和20 g/kg為例進行計算[11]。

表1 柴油機主要技術(shù)參數(shù)Tab.1 Main parameters of diesel engine

圖1 上止點時計算網(wǎng)格Fig.1 Computational mesh at TDC

1.2 熱物性參數(shù)表達式

流體熱物性包括流體的熱力學(xué)性質(zhì)參數(shù)和遷移(輸運)性質(zhì)參數(shù)，其中熱力學(xué)性質(zhì)參數(shù)除了壓力、溫度外，還包括比熱容、比熱比、比焓壓縮系數(shù)、普朗特數(shù)、聲速等，而遷移性質(zhì)參數(shù)則包括粘度、導(dǎo)熱系數(shù)等。

定壓比熱容

(1)

式中：γ為氣體比熱比；p為氣體壓力；T為氣體溫度；ρ為氣體密度。

比熱比

(2)

式中f(eth)為與比焓有關(guān)的函數(shù)。

根據(jù)Sutherland方程，氣體動力粘度

(3)

式中A1和A2為常數(shù)。

導(dǎo)熱系數(shù)

(4)

式中為A3和A4為常數(shù)。

壓縮系數(shù)

(5)

式中v為氣體體積。

比焓

(6)

式中u為比內(nèi)能。

普朗特數(shù)

(7)

當(dāng)?shù)匾羲?/p>

(8)

式中：k為多變指數(shù)；R為氣體常數(shù)。

1.3 數(shù)學(xué)模型和初始、邊界條件

本研究中油滴破碎模型選擇TAB模型，湍流模型為RNGk-ε模型，燃燒模型為Characteristictime模型[11]。計算從進氣門關(guān)閉開始(-150°CAATDC)到排氣門打開(124°CAATDC)結(jié)束。進氣門關(guān)閉時刻氣缸壓力為0.095MPa，缸內(nèi)溫度為320K。 計算工況：標(biāo)定轉(zhuǎn)速；75%負荷。在KIVA-3軟件輸入?yún)?shù)文件itape5中，通過改變H2O的質(zhì)量分數(shù)調(diào)節(jié)含濕量。

2 計算結(jié)果與分析

2.1 含濕量對熱力學(xué)參數(shù)的影響

圖2給出了不同含濕量時柴油機缸內(nèi)工質(zhì)熱力學(xué)參數(shù)的變化歷程?？梢钥闯?，含濕量對缸內(nèi)壓力影響較小，即對柴油機動力性影響不大，但對燃燒溫度有一定的影響。與含濕量為0時相比，含濕量增大到20g/kg時，缸內(nèi)最高燃燒溫度可降低30K。因此，濕度的增加有利于降低柴油機NOx排放。由于水的比熱容較大，因此含濕量增加時，混合氣比熱容增加，如圖2(c)所示。這就定性、定量解釋了進氣中水蒸汽含量增加后燃燒溫度降低的現(xiàn)象。

比熱比隨曲軸轉(zhuǎn)角變化情況如圖2(d)所示。比熱比不僅與氣體分子微觀自由度有關(guān)，隨自由度增加而減少，而且與溫度有關(guān)，隨著溫度的增加而減小。溫度增加，自由度隨之增加，因此比熱比減小。當(dāng)空氣中含濕量增加時，一方面三原子分子H2O含量增加，雙原子分子氮氣和氧氣含量減少，引起分子自由度增加，導(dǎo)致比熱比減小，另一方面缸內(nèi)溫度下降，也會導(dǎo)致比熱比減小。

圖2(e)給出了比焓隨曲軸轉(zhuǎn)角變化的情況。由圖可看出，在壓縮和做功過程中，比焓先增后減，隨著含濕量增加，比焓變化很小。由于比焓與濕度、溫度有關(guān)，濕度增加時，比焓增加，但是溫度也減小，因此比焓變化不大。

圖2(f)給出了壓縮系數(shù)隨曲軸轉(zhuǎn)角變化的情況。由圖可見，含濕量對壓縮系數(shù)幾乎沒有影響，在壓縮過程中，壓縮系數(shù)減小，壓縮難度加大，在上止點時壓縮難度最大，在最大爆發(fā)壓力后，壓縮系數(shù)逐漸增加。

普朗特數(shù)是流體力學(xué)中表征流體流動中能量交換中一個重要參數(shù)。普朗特數(shù)與動力粘度和比熱成正比，與導(dǎo)熱系數(shù)成反比。圖2(g)給出了普朗特數(shù)隨曲軸轉(zhuǎn)角變化的情況。在活塞開始上行時，溫度傳遞速度快，溫度邊界層厚度發(fā)展得快，普朗特數(shù)呈下降趨勢，在壓縮階段后期(-50°CA之后)，流動邊界層發(fā)展較快，因此普朗特數(shù)又呈增加趨勢。隨著濕度增加，定壓比熱增加較多，導(dǎo)熱系數(shù)下降和動力粘度下降幅度相當(dāng)，因此普朗特數(shù)仍增加。

圖2(h)給出了當(dāng)?shù)匾羲匐S曲軸轉(zhuǎn)角變化的情況。由于在壓縮和燃燒過程缸內(nèi)工質(zhì)狀態(tài)是變化的，因此當(dāng)?shù)匾羲僖彩亲兓摹穸仍黾?，缸?nèi)溫度和比熱比均減少，因此當(dāng)?shù)匾羲贉p小。

圖2 缸內(nèi)工質(zhì)熱力學(xué)參數(shù)變化規(guī)律Fig.2 Thermodynamic parameters of in-cylinder working fluids with crank angle

2.2 濕度對遷移參數(shù)的影響

柴油機缸內(nèi)工質(zhì)遷移參數(shù)變化規(guī)律如圖3所示。圖3(a)給出了導(dǎo)熱系數(shù)隨曲軸轉(zhuǎn)角變化的情況。由于水的導(dǎo)熱系數(shù)小于氮氣和氧氣的導(dǎo)熱系數(shù)，因此當(dāng)含濕量增加時，導(dǎo)熱系數(shù)減小。當(dāng)含濕量增加時，缸內(nèi)溫度減小，因此氣體導(dǎo)熱系數(shù)隨溫度減小而減小，在低溫時影響不明顯，在高溫時影響明顯。圖3(b)給出了動力粘度隨曲軸轉(zhuǎn)角變化的情況。動力粘度變化趨勢與導(dǎo)熱系數(shù)類似，均隨曲軸轉(zhuǎn)角變化是先增后減。相對于導(dǎo)熱系數(shù)，含濕量對動力粘度影響更小。

圖3 缸內(nèi)工質(zhì)遷移參數(shù)變化規(guī)律Fig.3 Migratory parameters of in-cylinder working fluids with crank angle

3 結(jié) 語

1)含濕量對柴油機缸內(nèi)熱效應(yīng)影響明顯，含濕量增加，缸內(nèi)溫度、比熱容和普朗特數(shù)下降較多。

2)含濕量對缸內(nèi)工質(zhì)壓力、比熱比和比焓影響較小，對壓縮系數(shù)影響很小。

3)缸內(nèi)工質(zhì)動力粘度與導(dǎo)熱系數(shù)變化規(guī)律類似。與低溫時相比，在高溫時隨著含濕量的增加，兩者增加較多。

[1]TOBACKAT,HEARNEJS,KURITZB,etal.Theeffectofambienttemperatureandhumidityonmeasuredidlingemissionsfromdieselschoolbuses[C]//SAEPaper:Detroit,USA, 2004,2004-01-1087.

[2]LINCherng-yuan,CHENWei-cheng.Effectsofpotassiumsulfidecontentinmarinedieselfueloilonemissioncharacteristicsofmarinefurnacesundervaryinghumidityofinletair[J].OceanEngineering,2006,33(8/9):1260-1270.

[3]CASATIR,SCHEERV,VOGTR,etal.Measurementofnucleationandsootmodeparticleemissionfromadieselpassengercarinrealworldandlaboratoryinsitudilution[J].AtmosphericEnvironment,2007,41(10):2125-2135.

[4]STRAWAAW,KIRCHSTETTERTW,HALLARAG,etal.Opticalandphysicalpropertiesofprimaryon-roadvehicleparticleemissionsandtheirimplicationsforclimatechange[J].JournalofAerosolScience,2010,41(1):36-50.

[5] 丁皓,吉曉燕,秦建華，等.高溫高壓下飽和濕空氣焓與濕度的預(yù)測[J].化工學(xué)報,2002,53(8):879-882.

DINGHao,JIXiao-yan,QINJian-hua,etal.Predictionofsaturatedhumidityandenthalpyofhumidairathightemperatureandpressure[J].JournalofChemicalIndustryandEngineering,2002,53(8):879-882.

[6] 楊文濱，蘇明.高濕空氣在高溫高壓區(qū)域內(nèi)的熱物性計算及分析[J].動力工程，2004,24(3):383-387,410.

YANGWen-bin,SUMing.Calculationandanalysisofthermodynamicpropertiesforhighhumidairinhightemperatureandhighpressure[J].PowerEngineering,2004,24(3):383-387,410.

[7] 劉朝，楊智勇，劉娟芳.高溫高壓下濕空氣的焓和熵計算[J].工程熱物理學(xué)報，2007,28(4)：557-560.

LIUChao,YANGZhi-yong,LIUJuan-fang.Thecalculationsofmoistair′senthalpyandentropyinhightemperatureandpressure[J].JournalofEngineeringThermophysics,2007,28(4)：557-560.

[8] 常勇強,曹子棟,趙振興，等.多組分氣體熱物性參數(shù)的計算方法[J].動力工程學(xué)報，2010，30(10)：772-776.

CHANGYong-qiang,CAOZi-dong,ZHAOZhen-xing,etal.Calculationmethodforthermalpropertiesofmulti-componentgas[J].JournalofChineseSocietyofPowerEngineering,2010，30(10)：772-776.

[9] 周見廣,臧述升,葛冰, 等.大濕度旋流擴散燃燒的變工況特性[J].上海交通大學(xué)學(xué)報，2004，38(4):1617-1622.

ZHOUJian-guang,ZANGShu-sheng,GEBing,etal.Thecharacterofswirldiffusionhumidaircombustionunderdifferentworkconditions[J].JournalofShanghaiJiaotongUniversity,2004，38(4):1617-1622.

[10] 顧欣，臧述升，葛冰.濕空氣擴散燃燒火焰結(jié)構(gòu)特性研究[J].工程熱物理學(xué)報，2006,27(2)：343-346.

GUXin,ZANGShu-sheng,GEBing.Characteristicsofnon-premixedflamestructureinhumidaircombustion[J].JournalofEngineeringThermophysics,2006,27(2)：343-346.

[11]VEYNANTED,VERVISCHL.Turbulentcombustionmodeling[J].ProgressinEnergyandCombustionScience,2002,28(3):193-266.

[12] 王向麗，王忠，倪培永，等.空氣濕度對柴油機NOx和碳煙排放影響的模擬分析[J].農(nóng)業(yè)工程學(xué)報，2012，28(12)：68-73.

WANGXiang-li,WANGZhong,NIPei-yong,etal.SimulationanalysisofeffectsofairhumidityonNOxandsootemissionsindieselengine[J].TransactionsoftheChineseSocietyofAgriculturalEngineering,2012，28(12)：68-73.

Effect of atmospheric humidity on thermal physical parameters of in-cylinder working fluids for diesel engine

Ni Pei-yong1,WANG Xiang-li2

(1.School of Mechanical Engineering, Nantong University,Nantong 226019,China;2.School of Electrical Engineering, Nantong University,Nantong 226019,China)

The expressions of related thermal physical parameters were embedded into KIVA-3 software.Those parameters include pressure, temperature, specific heat, specific heat ratio, specific enthalpy, thermal conductivity, kinetic viscosity, sound velocity, compression coefficient and Prandtl number.By changing inlet gas composition, the thermal physical parameters of in-cylinder working fluids were calculated during diesel combustion at different atmospheric humidity (0, 5, 10 and 10 g/kg).The simulation results show that humidity has important effect on combustion temperature, specific heat and Prandtl number, a little effect on pressure, specific heat ratio, specific enthalpy and sound velocity, and little effect on compression coefficient.With increasing humidity ratio, the thermal conductivity and kinetic viscosity increase slightly at high temperature.

diesel engine; humidity ratio; thermal physical parameters; simulation

2015-04-21；

2015-08-06

中國博士后科學(xué)基金資助項目(2013M541613)

倪培永(1977-)，男，博士，副教授，研究方向為內(nèi)燃機燃料與燃燒。

TK421.1

A

1672-7649(2015)12-0054-05

10.3404/j.issn.1672-7649.2015.12.011