線性與非線性傳熱過程的Cu rzon-Ah lborn熱機(jī)在任意功率時(shí)的效率?

李倩文 李瑩 張榮 盧燦燦 白龍

(中國礦業(yè)大學(xué)物理學(xué)院,徐州 221116)

線性與非線性傳熱過程的Cu rzon-Ah lborn熱機(jī)在任意功率時(shí)的效率?

李倩文 李瑩 張榮 盧燦燦 白龍?

(中國礦業(yè)大學(xué)物理學(xué)院,徐州 221116)

(2017年1月31日收到;2017年4月28日收到修改稿)

熱機(jī)性能的優(yōu)化是熱力學(xué)領(lǐng)域的一個(gè)重要問題,而工質(zhì)與熱源之間的傳熱過程是熱機(jī)工作時(shí)產(chǎn)生不可逆的主要來源.本文在引入功率增益和效率增益兩個(gè)重要參數(shù)的基礎(chǔ)上,基于一個(gè)簡(jiǎn)化的Curzon-Ahlborn熱機(jī)模型并利用合比分比原理,給出了線性與非線性傳熱過程的熱機(jī)在任意功率輸出時(shí)的效率表達(dá)式,結(jié)合數(shù)值計(jì)算詳細(xì)討論了熱機(jī)在任意功率輸出時(shí)的特性.研究表明,參數(shù)ξ作為功率增益δP的函數(shù)存在兩個(gè)分支：在第一分支上(不利情形),效率呈現(xiàn)出單調(diào)變化特征;在第二分支上(有利情形),效率隨著的δP變化是非單調(diào)的且有最大值.隨著傳熱指數(shù)的增加,熱機(jī)的工作區(qū)域減小,這源于非線性傳熱過程包含熱輻射所致.進(jìn)一步發(fā)現(xiàn)功率-效率關(guān)系曲線存在權(quán)衡工作點(diǎn),熱機(jī)在該點(diǎn)附近工作能夠?qū)崿F(xiàn)最有效的熱功轉(zhuǎn)換.研究結(jié)果有助于深入理解具有不同傳熱過程熱機(jī)的優(yōu)化執(zhí)行.

功率增益,效率增益,任意功率,效率

1 引 言

熱機(jī)是將外界輸入的熱能轉(zhuǎn)變成有用功的裝置,熱機(jī)效率的改進(jìn)成為第一次工業(yè)革命的引擎,這是自然科學(xué)指導(dǎo)人類生產(chǎn)實(shí)踐的典型范例.迄今為止,熱功轉(zhuǎn)換一直是熱力學(xué)領(lǐng)域中的一個(gè)重要問題,熱機(jī)性能的優(yōu)化成為熱力學(xué)發(fā)展的驅(qū)動(dòng)力之一.當(dāng)前,隨著經(jīng)濟(jì)的全球化與社會(huì)的快速發(fā)展,能源的需求導(dǎo)致不可再生能源(例如：石油,煤炭等)日益枯竭,尋找可再生能源和提高熱裝置的熱功轉(zhuǎn)換效率是人類得以可持續(xù)發(fā)展的根本保證.一個(gè)不爭(zhēng)的事實(shí)是：所有熱機(jī)的工作效率的上限是卡諾效率ηC=1-TL/TH,這里TH和TL分別是高低溫?zé)嵩吹臏囟?雖然基于平衡態(tài)熱力學(xué)的卡諾定理是熱力學(xué)第二定律的基石,但熱機(jī)要想達(dá)到卡諾效率,要求其工作過程為準(zhǔn)靜態(tài)過程,這導(dǎo)致熱機(jī)無功率輸出,致使卡諾熱機(jī)的實(shí)際應(yīng)用受到極大的限制.眾所周知,所有實(shí)際熱機(jī)都是在有限時(shí)間條件下工作的,Curzon和Ahlborn[1]在1975年基于內(nèi)部可逆假設(shè)提出了一個(gè)有限時(shí)間條件下的熱機(jī)模型,給出了熱機(jī)在最大功率輸出時(shí)的效率,即CA效率但在更早些時(shí)候,Reitlinger, Yvon,Chambadal和Novikov已經(jīng)做出了類似的結(jié)果,上述研究工作導(dǎo)致了“有限時(shí)間熱力學(xué)”成為熱力學(xué)的一個(gè)嶄新領(lǐng)域[2-4].

事實(shí)上,CA效率是線性非平衡熱力學(xué)框架下的熱機(jī)效率的一個(gè)基本結(jié)果[5].值得注意的是,CA效率的獲得是基于一個(gè)線性熱轉(zhuǎn)移過程,即熱源和工質(zhì)之間的熱傳導(dǎo)滿足牛頓熱傳導(dǎo)規(guī)律.然而,實(shí)際熱機(jī)在工作過程中工質(zhì)和熱源之間的傳熱過程通常是復(fù)雜的,這使得研究具有非線性傳熱過程熱機(jī)的性能優(yōu)化更加具有理論和工程實(shí)踐意義.1993年,Angu lo-Brown和Paez-Hernandez[6]考慮了熱源和工質(zhì)之間的傳熱過程滿足一個(gè)非線性傳熱定律,即Dulong-Petit定律,并詳細(xì)討論了滿足非線性傳熱過程的熱機(jī)特性.Huleihil和Andresen[7]研究了對(duì)流型傳熱過程對(duì)于一個(gè)內(nèi)部可逆熱機(jī)性能的影響,給出了該熱機(jī)優(yōu)化的功率和效率的近似解析公式.文獻(xiàn)[8-11]較為系統(tǒng)地研究了工質(zhì)與熱源之間的熱傳遞過程對(duì)于不可逆Carnot-like熱機(jī)功率與效率的影響,給出了優(yōu)化功率(效率)與系統(tǒng)參數(shù)的關(guān)系.Liang等[12]基于最小熵原則分析了多種傳熱過程條件下熱機(jī)的熱功轉(zhuǎn)換特征,給出了優(yōu)化熱機(jī)性能的一般性方法.這些工作極大地促進(jìn)了熱源與工質(zhì)傳熱過程(包括線性或非線性)如何影響熱機(jī)性能的相關(guān)問題研究[13-16].

另一方面,實(shí)際的熱電廠和真實(shí)熱機(jī)并不是在最大功率時(shí)工作的,這主要是因?yàn)樽畲蠊β瘦敵鰰r(shí)熱機(jī)的效率相對(duì)較小,而在比最大功率稍微小的功率工作時(shí),熱機(jī)能夠獲得較大的熱功轉(zhuǎn)換效率. W hitney[17,18]著重分析了量子器件在任意功率輸出時(shí)的性能,給出了量子熱機(jī)在有限功率輸出時(shí)效率和功率的上限.受W hitney工作的啟發(fā),Ryabov和Holubec[19]討論了宏觀低耗散熱機(jī)模型在任意功率時(shí)的優(yōu)化問題,給出了最大效率和功率增益因子的關(guān)系.因此,研究任意功率(非最大功率)輸出時(shí)熱機(jī)性能的優(yōu)化問題更具有工程實(shí)踐意義[20,21].然而,熱機(jī)的不可逆工作過程主要源自于熱源和工作流體之間的傳熱過程,這種熱傳遞過程如何影響熱機(jī)在任意功率輸出時(shí)的性能卻很少有人關(guān)注.因此,本文系統(tǒng)地研究了具有線性和非線性傳熱過程的熱機(jī)在任意功率輸出條件下的性能優(yōu)化問題.基于一個(gè)簡(jiǎn)化的CA熱機(jī)模型并且結(jié)合代數(shù)中的等比和分比原則,解析性地推導(dǎo)出了任意功率輸出時(shí)熱機(jī)的效率與功率增益之間關(guān)系,并進(jìn)一步討論了熱機(jī)性能優(yōu)化與系統(tǒng)參數(shù)的關(guān)系,研究結(jié)果有助于人們深入理解真實(shí)熱機(jī)的執(zhí)行過程.

2 線性與非線性傳熱Curzon-Ahlborn熱機(jī)和理論公式

本文所研究的CA熱機(jī)及其周圍環(huán)境如圖1(a)所示,這里TH和TL分別是高低溫?zé)嵩吹臏囟?THC和TLC分別是工質(zhì)(工作流體)和高低溫?zé)嵩吹葴亟佑|時(shí)的溫度(TH＞THC＞TLC＞TL).當(dāng)工質(zhì)在與高低溫?zé)嵩吹葴亟佑|過程中,其傳熱過程滿足如下規(guī)律：

這里α為傳熱系數(shù),n為傳熱指數(shù).當(dāng)n=1時(shí),對(duì)應(yīng)的是線性傳熱規(guī)律;當(dāng)n=5/4時(shí),對(duì)應(yīng)的是非線性傳熱過程的Dulong-Petit定律.積分方程(1)和(2),可以獲得工質(zhì)與高低溫?zé)嵩唇佑|過程中吸收和放出的熱量：

圖1 (a)CA熱機(jī)模型;(b)簡(jiǎn)化的CA熱機(jī)模型Fig.1.(a)The CA heat enginemodel;(b)a sim p lified version of the CA heat engine.

這里面th和tc分別是工質(zhì)與高低溫?zé)嵩吹慕佑|時(shí)間.相比于兩個(gè)等溫過程,熱機(jī)的兩個(gè)絕熱過程持續(xù)的時(shí)間較短,可以忽略不計(jì),這種近似滿足CA熱機(jī)的內(nèi)部可逆假設(shè)要求.根據(jù)熱力學(xué)第二定律,兩個(gè)等溫過程有如下關(guān)系,

熱機(jī)的功率(P)和效率(η)可以分別表示為：

為了研究熱機(jī)在任意功率時(shí)的性能,我們引入了功率增益因子(δP)和效率增益因子(δη),分別表示為：

其中Pmax為熱機(jī)的最大功率,ηmaxP為熱機(jī)最大功率時(shí)的效率,且-1≤δP≤0.

為了優(yōu)化CA熱機(jī)的性能,我們定義工質(zhì)與高溫?zé)嵩唇佑|時(shí),工質(zhì)和熱源之間的溫差為x=TH-THC;工質(zhì)與低溫?zé)嵩唇佑|而放熱時(shí),工質(zhì)與低溫?zé)嵩吹臏夭顬閥=TLC-TL.參數(shù)(x,y)為深入研究熱機(jī)的執(zhí)行提供了有效途徑.如果令x=y,經(jīng)典的CA熱機(jī)模型可以約化為一個(gè)簡(jiǎn)化的CA熱機(jī)模型[22],這個(gè)簡(jiǎn)化的CA熱機(jī)模型能夠在優(yōu)化熱機(jī)性能時(shí)避免繁雜的數(shù)學(xué)計(jì)算,更為重要的是它能夠較好地反映不可逆熱機(jī)的工作過程,并且已被證明是一個(gè)有效的熱機(jī)模型[23],簡(jiǎn)化的CA熱機(jī)模型如圖1(b)所示.這樣 QHC=αxnth和QLC=αxntc,基于方程(5)可以得到如下公式：

利用代數(shù)中的合比分比原理,可以得到

關(guān)于合比分比原理,可以詳見附錄A.利用方程(6)和(10),非線性傳熱熱機(jī)的功率可以表示為

將功率P對(duì)變量x求導(dǎo),并且滿足極值條件,即d P/d x=0.這樣,滿足最大功率輸出時(shí)

將方程(12)代入方程(11),最大功率為

將方程(12)代入方程(7)并結(jié)合方程(9),最大功率時(shí)的效率為

結(jié)合方程(14)和(16),任意功率時(shí)的效率可以表示為

圖2 (網(wǎng)刊彩色)P/Pmax與ξ的關(guān)系曲線Fig.2.(color online)The relation curve of P/Pmax vsξ.

理解任意功率時(shí)熱機(jī)執(zhí)行的關(guān)鍵是建立效率與功率增益的關(guān)系.基于方程(15),功率增益因子δP與ξ的關(guān)系為

當(dāng)n=1時(shí)(即傳熱過程滿足線性規(guī)律),可以得到方程(18)的2個(gè)解析解(兩個(gè)分支)：

當(dāng)取n=5/4(即傳熱過程滿足非線性規(guī)律)時(shí),解析求解(18)是非常困難的.利用數(shù)值計(jì)算,可以給出方程(18)的兩個(gè)數(shù)值解(兩個(gè)分支),這里假設(shè)數(shù)值解可以寫為

將(19)-(22)式代入(17)式,任意功率時(shí)的效率(η)與功率增益因子(δP)的關(guān)系可以被確立,這為進(jìn)一步認(rèn)識(shí)熱機(jī)以有限功率輸出時(shí)的性能提供了有力工具.

3 數(shù)值計(jì)算結(jié)果與討論

圖3 (網(wǎng)刊彩色).與ξ的關(guān)系曲線.Fig.3.(color on line)The relation cu rve of vsξ.

為了進(jìn)一步研究任意功率輸出時(shí)熱機(jī)的效率,認(rèn)識(shí)ξ與功率增益δP的關(guān)系是非常重要的.我們?cè)诘诙糠纸o出了兩種傳熱過程條件下ξ與功率增益δP的關(guān)系[見方程(19)-(22)],基于這些函數(shù)關(guān)系式,圖4給出了 ξ隨著δP的變化特征曲線.兩種傳熱過程條件下ξ-δP的特征曲線均包括兩個(gè)分支：熱機(jī)工作在第1分支(也稱不利分支)上時(shí),當(dāng)x相對(duì)于x?減小時(shí),功率的減小伴隨著效率的降低;熱機(jī)操作在第2分支(也稱有利分支)上時(shí),當(dāng)x相對(duì)于x?增加時(shí),功率的增加伴隨著效率的增強(qiáng),甚至大于最大功率時(shí)的效率ηmaxP,參見圖2和圖3.非線性傳熱過程的ξ-δP特征曲線包含在線性傳熱過程的ξ-δP特征曲線之中.基于方程(17)并結(jié)合關(guān)系式(19)-(22),圖5和圖6分別給出了線性和非線性傳熱情況下熱機(jī)效率η作為溫度比(TL/TH)和效率增益(δP)函數(shù)的等高線圖.從圖5(a)和圖6(a)可以看出,對(duì)于ξ取第1分支時(shí),η隨著TL/TH或者|δP|的增加單調(diào)減少;當(dāng)ξ取第2分支時(shí),η隨著TL/TH增加而減小,而隨著|δP|的增加先達(dá)到最大值,然后逐漸減小甚至消失.值得注意的是,雖然兩種傳熱過程熱機(jī)效率的等高曲線是相似的,但是線性傳熱熱機(jī)工作區(qū)域大于非線性傳熱熱機(jī)的工作區(qū)域.線性傳熱過程僅由熱傳導(dǎo)支配工質(zhì)和熱源之間的熱交換,而非線性傳熱過程還包含熱輻射過程,這導(dǎo)致熱機(jī)工作時(shí)的不可逆程度和熱損失增加.因此,熱機(jī)工作在非線性傳熱過程時(shí),工作區(qū)域以及效率會(huì)相應(yīng)減小.

圖4 (網(wǎng)刊彩色)ξ與δP的關(guān)系曲線Fig.4.(color on line)The relation cu rve ofξvsδP.

為了進(jìn)一步闡釋不同傳熱過程對(duì)于熱機(jī)效率的影響,在固定TL/TH=0.5的情況下,圖7給出了兩種傳熱條件下任意功率時(shí)的效率η與功率增益δP的關(guān)系曲線.從圖7可以看出,無論是線性或者非線性傳熱過程,η均是δP的雙值函數(shù),這是ξ-δP特征曲線的兩分支結(jié)構(gòu)在效率η與功率增益δP關(guān)系曲線中的反映.發(fā)現(xiàn)η-δP的曲線是閉合的“環(huán)”狀結(jié)構(gòu),除了重合之處外,線性傳熱過程的效率大于非線性傳熱過程的效率,這源于較復(fù)雜的非線性傳熱過程引起熱損耗增加所致.進(jìn)一步發(fā)現(xiàn),ξ第2分支(虛線)上的效率值明顯大于第1分支(實(shí)線)上的效率,因?yàn)闊釞C(jī)工作在有利分支上時(shí),能實(shí)現(xiàn)更大效率的輸出.另外,η隨著δP變化的最大值出現(xiàn)在|δP|較小的值附近而不是δP=0處,即熱機(jī)以比最大功率稍小的功率輸出時(shí)能夠獲得最大效率.這意味著熱機(jī)工作性能曲線存在著權(quán)衡工作點(diǎn),在該點(diǎn)附近工作的熱機(jī)能實(shí)現(xiàn)最有效的熱功轉(zhuǎn)換.從應(yīng)用角度來看,這些結(jié)果為如何優(yōu)化熱機(jī)的執(zhí)行提供了有效途徑.

圖5 (a)線性傳熱情況下效率η作為TC/TH和δP函數(shù)的等高線圖,ξ取第1分支; (b)線性傳熱情況下效率η作為TC/TH和δP函數(shù)的等高線圖,ξ取第2分支Fig.5.(a)Forξ=ξ1,the contou r-line p lot ofηas functions of TC/TH andδP with the linear heat transfer process; (b)forξ=ξ2,the contou r-line p lot ofηas functions of TC/TH andδP with the linear heat transfer process.

圖6 (a)非線性傳熱情況下效率η作為TC/TH和δP函數(shù)的等高線圖,ξ取第1分支; (b)線性傳熱情況下效率η作為TC/TH和δP函數(shù)的等高線圖,ξ取第2分支Fig.6.(a)Forξ=ξ1,the contou r-line p lot ofηas functions of TC/TH andδP with the nonlinear heat transfer process;(b)forξ=ξ2,the contour-line p lot ofηas functions of TC/TH andδP with the nonlinear heat transfer process.

最后,根據(jù)方程(16),功率增益δη表示為

方程(19)-(22)提供了線性與非線性傳熱情況下ξ與功率增益δP的公式,結(jié)合(23)式與方程(19)-(22),可以獲得δη-δP的關(guān)系曲線.對(duì)于線性傳熱過程,分別對(duì)應(yīng)于ξ第1分支和第2分支的功率增益與效率增益的關(guān)系為

而對(duì)于非線性傳熱過程,分別對(duì)應(yīng)于ξ第1分支[方程(21)]和第2分支[方程(22)]的δη-δP的關(guān)系為

基于方程(24)-(27),圖8給出了功率增益δη與效率增益δP的關(guān)系曲線.從圖8可以看出,兩種傳熱情況下的 δη-δP曲線均為“環(huán)”形狀的閉合曲線,在第1分支上曲線幾乎重合(實(shí)黑線和實(shí)紅線所示),而在第2分支上線性和非線性傳熱過程的δη-δP曲線略有差別(黑虛線和紅虛線所示),這種微小的差異源于不同傳熱過程對(duì)于熱機(jī)工作的影響.

圖7 (網(wǎng)刊彩色)任意功率時(shí)的效率η與功率增益δP的關(guān)系,對(duì)應(yīng)卡諾效率ηC=1-TC/TH=0.5Fig.7.(color on line)The relation curve of the efficiency at arbitrary power(η)asa function of the power gain(δP).Here,the correspond ing Carnot efficiency ηC=1-TC/TH=0.5.

圖8 (網(wǎng)刊彩色)效率增益δη與功率增益δP的關(guān)系,對(duì)應(yīng)卡諾效率ηC=1-TC/TH=0.5Fig.8.(color on line)The relation cu rve of the efficiency gain(δη)as a function of the power gain (δP).Here,the corresponding Carnot efficiencyηC= 1-TC/TH=0.5.

4 結(jié) 論

熱機(jī)性能的優(yōu)化對(duì)于工程實(shí)踐和國民經(jīng)濟(jì)的可持續(xù)發(fā)展有著重要意義.實(shí)際熱機(jī)在工作過程中輸出最大功率時(shí)其效率并非最大,而以略小于最大功率輸出時(shí)能夠獲得最大效率.本文基于一個(gè)簡(jiǎn)化的CA熱機(jī)模型,分析了熱機(jī)工作過程中不同的傳熱過程對(duì)于熱機(jī)性能的影響.基于兩個(gè)引入的參數(shù)：功率增益和效率增益,并且結(jié)合分比合比定理,給出了線性和非線性傳熱過程熱機(jī)的最大功率、最大功率時(shí)的效率以及任意功率時(shí)效率的普遍公式.進(jìn)一步發(fā)現(xiàn)參數(shù)ξ作為功率增益δP的函數(shù)存在2個(gè)分支.在此基礎(chǔ)上,我們?cè)敿?xì)地給出了不同傳熱過程的熱機(jī)在任意功率輸出時(shí)的效率與系統(tǒng)參數(shù)的關(guān)系.發(fā)現(xiàn)含非線性傳熱過程的熱機(jī)的工作區(qū)域小于線性傳熱過程的熱機(jī),這主要因?yàn)榉蔷€性傳熱過程存在著熱輻射,而且隨著傳熱指數(shù)的增加, η/ηmax的最大值略有增加且左移.任意功率時(shí)效率η作為TL/TH和δP函數(shù)的等高線圖表明：對(duì)應(yīng)于ξ-δP的第1分支,η隨著TL/TH和|δP|增加呈現(xiàn)出單調(diào)減少的趨勢(shì);對(duì)應(yīng)于ξ-δP的第2分支,η隨著TL/TH增加而單調(diào)變化,但η隨著|δP|的增加呈現(xiàn)出非單調(diào)特征,在TL/TH與|δP|均較小的區(qū)域存在效率最大值.非線性傳熱過程的熱機(jī)的執(zhí)行區(qū)域略有減小,這主要源于非線性傳熱過程的復(fù)雜性所致.兩種傳熱過程的熱機(jī)的η-δP的曲線具有閉合“環(huán)”狀特征,當(dāng)|δP|→0和|δP|→1時(shí),熱機(jī)的η-δP的曲線存在一定的重疊,而在其他區(qū)域,傳熱指數(shù)較小(線性傳熱過程)的熱機(jī)效率大于傳熱指數(shù)較大(非線性傳熱過程)的熱機(jī)效率;在熱機(jī)以比最大功率稍小的功率輸出工作時(shí),可以實(shí)現(xiàn)最大效率的執(zhí)行,這表明熱機(jī)實(shí)現(xiàn)最佳性能工作時(shí)存在著權(quán)衡工作點(diǎn).另外,兩種傳熱過程的熱機(jī)效率增益和功率增益的關(guān)系曲線也呈現(xiàn)出閉合“環(huán)”形狀,在ξ-δP的第2分支上略有差別.本文的研究結(jié)果不僅能夠幫助我們深入理解不同傳熱過程如何影響實(shí)際熱機(jī)的性能,而且對(duì)優(yōu)化受控于不同傳熱過程實(shí)際熱機(jī)也是有價(jià)值的.

附錄A 合比分比定理

在此我們對(duì)正文中用到的合比分比定理給予簡(jiǎn)要介紹.假設(shè)

這里面α/=β,β/=0,ν/=0.令方程(A 1)兩邊同時(shí)加上1,得到如下等式

再令方程(A 1)兩邊同時(shí)減去1,得到如下等式

讓方程(A 3)除以方程(A 2),則

這個(gè)恒等式還可以寫為

方程(A 4)就是著名的合比分比定理,在正文中推導(dǎo)方程(10)的過程中,我們用到了此定理.

[1]Curzon F,Ah lborn B 1975 Am.J.Phys.43 22

[2]Vaud rey A,Lanzetta F,Feid t M 2014 J.Non-Equil. Therm.39 199

[3]And resen B 2011 Angew.Chem.In t.Ed.50 2690

[4]Chen L G,Sun F R 1998 Acta Phys.Sin.18 395(in Chinese)[陳林根,孫豐瑞1998物理學(xué)報(bào)18 395]

[5]van den B roeck C 2005 Phys.Rev.Lett.95 190602

[6]Angu lo-B row n F,Phez-Hernhndez F 1993 J.App l.Phys. 74 2216

[7]Hu leihil M,And resen B 2006 J.Appl.Phys.100 014911

[8]Chen L,Sun F,W u C 1999 J.Phys.D:App l.Phys.32 99

[9]Zhou S,Chen L,Sun F,W u C 2005 App l.Energy 81 376

[10]Chen L,Li J,Sun F 2008 Appl.Energy 85 52

[11]Chen L,Sun F,W u C 2006 Appl.Energy 83 71

[12]Cheng X T,W angW H,Liang X G 2012 Chin.Sci.Bu ll. 55 2847

[13]A rias-Hernandez L A,A res de Prarga G,Angu lo-B rown F 2003 Open Sys.Inform at.Dyn.10 351

[14]Li J,Chen L,Sun F 2007 Appl.Energy 84 944

[15]Shu LW,Chen L,Sun F 2009 Sci.China Ser.B:Chem. 52 1154

[16]Chim al-Eguia J C,Barranco-Jim enez M A 2006 Open Sys.Inform at.Dyn.13 43

[17]W hitney R S 2014 Phys.Rev.Lett.112 130601

[18]W hitney R S 2015 Phys.Rev.B 91 115425

[19]Holubec V,Ryabov A 2016 J.Stat.M ech.2016 073204

[20]Long R,Liu W 2016 Phys.Rev.E 90 052114

[21]Ryabov A,Holubec V 2016 Phys.Rev.E 93 050101

[22]Agrawal D C 2009 Eur.J.Phys.30 1173

[23]Paez-Hernandz R T,Portillo-D iaz P,Ladino-Luna D 2016 J.Non-Equil.Therm.41 19

(Received 31 January 2017;revised manuscript received 28 April 2017)

Efficiency at arbitrary power for the Curzon-Ahlborn heat engine in linear and nonlinear heat transferprocesses?

LiQian-Wen Li Ying Zhang Rong Lu Can-Can Bai Long?

(School of Physics,China University ofM ining and Technology,Xuzhou 221116,China)

The optimal performance of heat engine is an im portantissue in thermodynam ics,but the heat transfer between the working medium and two heat reservoirs induces the irreversibility during the operation of heat engine.Based on two im portant param eters introduced in this paper(nam ely,the power gain and the efficiency gain),for heat engine operating in the linear and nonlinear heat transfer processes,the formula for the efficiency at arbitrary power is achieved in terms of a simplified Curzon-Ah lborn heat enginemodel and the“componendo and dividendo”rule.The features of heat engine at arbitrary power output are also discussed in detail based on the num erical calculations.It is indicated that the parameterξas a function of the power gainδP contains two branches:the efficiency shows themonotonous variation on the first branch(the favorable case);the efficiency exhibits the non-m onotonous characteristics and has the maximum value on the second branch(the un favorable case).The working region of the heat engine is reduced as the heat transfer exponent increases,which results from the radiative contribution in the nonlinear heat transfer process. For the fi rst branch,the contour-line p lot ofηversus TL/THandδP clearly dem onstrates thatηhas the decreasing trend with increasing TL/THand|δP|;for the second branch,ηmonotonically deceases as TL/THincreases,butηshows the non-monotonic behaviors as|δP|increases.The efficiency has themaximum value in the region where TL/THand|δP| have the sm all values,and the working regim e of heat engines in the nonlinear heat transfer process is relatively sm all due to the com plexity of the nonlinear heat transfer process.The curves of the efficiency in two heat transfer processes are loop-shaped,when|δP|→ 0 and|δP|→ 1,the curves ofη～δP in two heat transfer processes are sam e.But in other regim es,the efficiency of the heat engine with the linear heat transfer process is bigger than in the nonlinear heat transfer process.Furthermore,it is found that a considerably larger efficiency can be obtained when heat engine working close to them aximum power.This im p lies that there exists the trade-off working point where the heat engine can perform the m ost effective heat-work conversion.In addition,the curves of the power gain vs.the efficiency gain also display the loop-shaped characteristics,but there is the weak difference on the second branch.Our results are very conducive to understanding the op tim al perform ance of heat engines in different heat transfer processes.

power gain,efficiency gain,arbitrary power,efficiency

PACS:05.70.Ln,05.20.-y,05.70.-a DO I:10.7498/aps.66.130502

?中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金(批準(zhǔn)號(hào):2015XKM S082)資助的課題.

?通信作者.E-m ail:bailong2200@163.com

PACS:05.70.Ln,05.20.-y,05.70.-a DO I:10.7498/aps.66.130502

*Pro ject supported by the Fundam ental Research Funds for the Central Universities,China(G rant No.2015XKM S082).?

.E-m ail:bailong2200@163.com