1、自上世紀(jì)末期開(kāi)始，計(jì)算機(jī)技術(shù)的發(fā)展取得了巨大進(jìn)步，目前柴油機(jī)領(lǐng)域越來(lái)越多的依賴于CFD，數(shù)值模擬可以提供實(shí)驗(yàn)難以給出的流動(dòng)、傳熱細(xì)節(jié)，在發(fā)動(dòng)機(jī)設(shè)計(jì)、性能測(cè)試和優(yōu)化等方面體現(xiàn)了重要價(jià)值。柴油機(jī)缸內(nèi)輻射傳熱數(shù)值模擬的關(guān)鍵包括輻射傳輸方程的求解、流動(dòng)與輻射耦合求解、高溫燃?xì)庵蠧O2和H2O的物性求解。目前，對(duì)缸內(nèi)輻射換熱的研究認(rèn)識(shí)相對(duì)較少。本文在自主開(kāi)發(fā)的計(jì)算多物理軟件GTEA基礎(chǔ)上，添加基于格心型有限體積法的輻射模型，實(shí)現(xiàn)流動(dòng)與輻射耦合計(jì)算

2、，并模擬了缸內(nèi)輻射換熱，主要內(nèi)容如下:
　　（1）建立輻射傳輸方程求解計(jì)算程序。針對(duì)吸收、發(fā)射、散射性介質(zhì)的輻射換熱問(wèn)題，基于格心型有限體積法(FVM)求解輻射傳遞方程，采用FORTRAN語(yǔ)言編寫輻射換熱求解程序。輻射計(jì)算可以與流場(chǎng)計(jì)算采用同一套網(wǎng)格，簡(jiǎn)化了計(jì)算程序，提高了計(jì)算的通用性。當(dāng)網(wǎng)格和立體角的離散個(gè)數(shù)滿足一定的對(duì)應(yīng)關(guān)系時(shí)，才能得到更精確的計(jì)算結(jié)果。當(dāng)壁面為冷邊界時(shí)，壁面輻射熱流隨著吸收系數(shù)的增加而增大，隨著散射反照率的增

3、加逐漸減小。
　?。?）實(shí)現(xiàn)軟件GTEA已有的流場(chǎng)模塊與添加的輻射模塊的耦合。針對(duì)封閉方腔內(nèi)含吸收、發(fā)射及各向同性散射介質(zhì)的輻射與自然對(duì)流耦合傳熱問(wèn)題，采用有限體積法進(jìn)行了詳細(xì)研究。并分析普朗克數(shù)（0.0001≤Pl≤10）、衰減系數(shù)（0.2≤β0≤5）、溫度比（0.25≤γ≤0.9）、散射反照率（0≤ω≤1）、傾斜角（-60°≤α≤90°）和長(zhǎng)寬比（1≤AR≤5）等參數(shù)對(duì)等溫線、流線、溫度、速度和Nusselt數(shù)的影響。

4、　　（3）采用合并寬窄譜帶關(guān)聯(lián)k分布模型(CWNBCK)計(jì)算非灰氣體物性。對(duì)于缸內(nèi)輻射換熱問(wèn)題，不需要知道譜帶詳細(xì)的光譜數(shù)據(jù)，主要關(guān)心譜帶的輻射換熱量。根據(jù)文獻(xiàn)中合并寬窄譜帶的方法，本文進(jìn)一步將窄譜帶合并成兩個(gè)寬帶。采用合并后的寬帶k分布模型計(jì)算CO2和H2O混合氣體總的吸收系數(shù)。本方法兼顧了計(jì)算精度與計(jì)算效率，適用于缸內(nèi)燃?xì)廨椛湮镄缘挠?jì)算。
　?。?）考慮缸內(nèi)氣相對(duì)液滴的輻射換熱。軟件GTEA已有的液滴蒸發(fā)模塊忽略了熱輻射的影響

5、，本文考慮輻射能以類似于對(duì)流換熱的形式作用于液滴表面。采用B.Abramzon與S. S. Sazhin簡(jiǎn)化的輻射吸收模型及Dombrovsky和Sazhin給出的平均吸收因子的最佳近似表達(dá)式，編寫液滴輻射換熱計(jì)算程序。通過(guò)單液滴算例研究表明：當(dāng)氣相溫度較低（673K）時(shí)，輻射傳熱對(duì)液滴直徑和溫度的影響較小。
　　（5）實(shí)現(xiàn)燃燒室內(nèi)兩相燃?xì)饬髦械妮椛渑c對(duì)流耦合傳熱數(shù)值模擬。運(yùn)用軟件GTEA中已有的傳熱、流動(dòng)、噴霧、化學(xué)反應(yīng)等數(shù)值求