線性與非線性傳熱過程的Cu rzon-Ah lborn熱機在任意功率時的效率?

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

(中國礦業大學物理學院,徐州 221116)

線性與非線性傳熱過程的Cu rzon-Ah lborn熱機在任意功率時的效率?

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

(中國礦業大學物理學院,徐州 221116)

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

熱機性能的優化是熱力學領域的一個重要問題,而工質與熱源之間的傳熱過程是熱機工作時產生不可逆的主要來源.本文在引入功率增益和效率增益兩個重要參數的基礎上,基于一個簡化的Curzon-Ahlborn熱機模型并利用合比分比原理,給出了線性與非線性傳熱過程的熱機在任意功率輸出時的效率表達式,結合數值計算詳細討論了熱機在任意功率輸出時的特性.研究表明,參數ξ作為功率增益δP的函數存在兩個分支：在第一分支上(不利情形),效率呈現出單調變化特征;在第二分支上(有利情形),效率隨著的δP變化是非單調的且有最大值.隨著傳熱指數的增加,熱機的工作區域減小,這源于非線性傳熱過程包含熱輻射所致.進一步發現功率-效率關系曲線存在權衡工作點,熱機在該點附近工作能夠實現最有效的熱功轉換.研究結果有助于深入理解具有不同傳熱過程熱機的優化執行.

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

1 引 言

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

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

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

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

本文所研究的CA熱機及其周圍環境如圖1(a)所示,這里TH和TL分別是高低溫熱源的溫度,THC和TLC分別是工質(工作流體)和高低溫熱源等溫接觸時的溫度(TH＞THC＞TLC＞TL).當工質在與高低溫熱源等溫接觸過程中,其傳熱過程滿足如下規律：

這里α為傳熱系數,n為傳熱指數.當n=1時,對應的是線性傳熱規律;當n=5/4時,對應的是非線性傳熱過程的Dulong-Petit定律.積分方程(1)和(2),可以獲得工質與高低溫熱源接觸過程中吸收和放出的熱量：

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

這里面th和tc分別是工質與高低溫熱源的接觸時間.相比于兩個等溫過程,熱機的兩個絕熱過程持續的時間較短,可以忽略不計,這種近似滿足CA熱機的內部可逆假設要求.根據熱力學第二定律,兩個等溫過程有如下關系,

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

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

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

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

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

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

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

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

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

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

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

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

當n=1時(即傳熱過程滿足線性規律),可以得到方程(18)的2個解析解(兩個分支)：

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

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

3 數值計算結果與討論

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

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

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

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

圖5 (a)線性傳熱情況下效率η作為TC/TH和δP函數的等高線圖,ξ取第1分支; (b)線性傳熱情況下效率η作為TC/TH和δP函數的等高線圖,ξ取第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函數的等高線圖,ξ取第1分支; (b)線性傳熱情況下效率η作為TC/TH和δP函數的等高線圖,ξ取第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.

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

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

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

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

圖7 (網刊彩色)任意功率時的效率η與功率增益δP的關系,對應卡諾效率η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 (網刊彩色)效率增益δη與功率增益δP的關系,對應卡諾效率η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 結 論

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

附錄A 合比分比定理

在此我們對正文中用到的合比分比定理給予簡要介紹.假設

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

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

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

這個恒等式還可以寫為

方程(A 4)就是著名的合比分比定理,在正文中推導方程(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物理學報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

?中央高校基本科研業務費專項資金(批準號: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