• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    有機(jī)朗肯循環(huán)中的工質(zhì)熱源轉(zhuǎn)折溫度及其特征

    2017-09-06 03:00:38李新國(guó)王偉翟哲南來(lái)福王競(jìng)逸
    化工進(jìn)展 2017年9期
    關(guān)鍵詞:臨界溫度工質(zhì)熱源

    李新國(guó),王偉,翟哲,南來(lái)福,王競(jìng)逸

    (天津大學(xué)機(jī)械工程學(xué)院,天津 300350)

    有機(jī)朗肯循環(huán)中的工質(zhì)熱源轉(zhuǎn)折溫度及其特征

    李新國(guó),王偉,翟哲,南來(lái)福,王競(jìng)逸

    (天津大學(xué)機(jī)械工程學(xué)院,天津 300350)

    有機(jī)朗肯循環(huán)中存在有極值功循環(huán)與無(wú)極值功循環(huán),即存在一種凈功的轉(zhuǎn)折線與所對(duì)應(yīng)的工質(zhì)熱源轉(zhuǎn)折溫度,它對(duì)工質(zhì)的選擇與循環(huán)性能評(píng)價(jià)具有指標(biāo)性的作用?;谔菪窝h(huán)及其理論模型,本文建立工質(zhì)飽和線形狀(飽和液線斜率、飽和氣線斜率及其斜率比與臨界溫度)及物性為自變量的數(shù)學(xué)關(guān)系,并提出和定義線性飽和線工質(zhì)模型,建立循環(huán)性能與工質(zhì)物性之間的數(shù)學(xué)關(guān)系。分別推導(dǎo)出有機(jī)朗肯循環(huán)中工質(zhì)熱源轉(zhuǎn)折溫度的經(jīng)驗(yàn)公式,以及梯形循環(huán)下線性飽和線工質(zhì)的熱源轉(zhuǎn)折溫度理論公式及其修正。研究熱源轉(zhuǎn)折溫度及其對(duì)應(yīng)的優(yōu)化工況的特征與一般性規(guī)律,得到21種工質(zhì)熱源轉(zhuǎn)折溫度的有機(jī)朗肯循環(huán)數(shù)值計(jì)算與經(jīng)驗(yàn)公式之間的偏差小于1.97%。線性飽和線工質(zhì)模型下,工質(zhì)的熱源轉(zhuǎn)折溫度為臨界溫度和蒸發(fā)器窄點(diǎn)溫差之和;熱源轉(zhuǎn)折溫度所對(duì)應(yīng)的最高優(yōu)化工況(最高優(yōu)化蒸發(fā)溫度與最大凈功極值)隨斜率比單調(diào)增加,至干工質(zhì)時(shí),最高優(yōu)化工況接近工質(zhì)的臨界點(diǎn)。

    有機(jī)朗肯循環(huán);梯形循環(huán);線性飽和線工質(zhì);工質(zhì)熱源轉(zhuǎn)折溫度;最高優(yōu)化工況

    有機(jī)朗肯循環(huán)(organic Rankine cycle,ORC)采用低沸點(diǎn)有機(jī)工質(zhì)發(fā)電,可以回收利用300℃以下的中低溫?zé)崮?,如太?yáng)能、地?zé)崮堋⒐I(yè)余熱等。目前ORC的研究一般是針對(duì)具體工質(zhì)的數(shù)值計(jì)算或?qū)嶒?yàn)研究。徐榮吉等[1]以R245fa為工質(zhì),提出基于對(duì)數(shù)傳熱溫差的內(nèi)回?zé)崞餍阅苡?jì)算方法,并分析過(guò)熱溫度、過(guò)冷溫度對(duì)內(nèi)回?zé)酧RC性能的影響,表明內(nèi)回?zé)釡p少了循環(huán)的蒸發(fā)負(fù)荷和冷凝負(fù)荷,提高了循環(huán)效率。韓中合等[2]針對(duì)煙氣余熱構(gòu)建無(wú)回?zé)岷陀袃?nèi)回?zé)岬腛RC系統(tǒng),通過(guò)熱源參數(shù)變化引起工質(zhì)吸熱量的變化,分析工質(zhì)在兩種系統(tǒng)中的初溫、凈功量、熱耗率及?損的變化規(guī)律,得出較優(yōu)工質(zhì)和各工質(zhì)對(duì)內(nèi)回?zé)崞鞯钠ヅ湫?。楊緒飛等[3]在給定熱源條件下,探討進(jìn)口過(guò)熱度對(duì)膨脹機(jī)性能和系統(tǒng)性能的影響,表明隨膨脹機(jī)進(jìn)口過(guò)熱度遞減,膨脹機(jī)機(jī)械效率遞增,等熵效率遞減,膨脹機(jī)軸功和實(shí)際運(yùn)行效率呈先增后減的變化趨勢(shì)。王明濤等[4]建立能量平衡方程和?方程,研究蒸發(fā)壓力對(duì)不同烷烴類工質(zhì)的熱效率、單位質(zhì)量煙氣凈功、?效率等參數(shù)的影響規(guī)律,表明在保證膨脹機(jī)體積流量比不大于50情況下,環(huán)戊烷具有最高的熱效率和?效率。張軍輝等[5]以最大做功能力和?效率為目標(biāo)函數(shù),對(duì)10種工質(zhì)進(jìn)行分析,表明每種工質(zhì)均存在一優(yōu)化蒸發(fā)溫度使循環(huán)凈輸出功最大,而且工質(zhì)臨界溫度越高,對(duì)應(yīng)的優(yōu)化蒸發(fā)溫度也越高。賀超等[6]采用不同優(yōu)化目標(biāo),對(duì)亞臨界ORC進(jìn)行蒸發(fā)溫度和冷凝溫度的優(yōu)化分析。在排煙溫度423.15K條件下,對(duì)干工質(zhì),不同優(yōu)化目標(biāo)下蒸發(fā)溫度和冷凝溫度優(yōu)化值差異較大;濕工質(zhì)R134a與R152a臨界溫度低于熱源初始溫度20K±2K時(shí),系統(tǒng)存在優(yōu)化蒸發(fā)溫度。LI等[7]在熱源溫度120~200℃下,對(duì)5種有機(jī)工質(zhì)從安全、環(huán)保、經(jīng)濟(jì)效益方面進(jìn)行分析,表明R123非常不安全,但它的熱經(jīng)濟(jì)效率最高,尤其熱源溫度較高時(shí)。YU等[8]研究表明,當(dāng)ORC系統(tǒng)在近臨界溫度區(qū)域內(nèi)運(yùn)行,即熱源進(jìn)口溫度與工質(zhì)臨界溫度之差在適當(dāng)范圍內(nèi)時(shí),系統(tǒng)可獲得最大的輸出功。

    也有學(xué)者從理論角度對(duì)有機(jī)朗肯循環(huán)進(jìn)行分析。最早由嚴(yán)家祿[9]對(duì)兩種熱源形式的ORC,分別推導(dǎo)了工質(zhì)優(yōu)化蒸發(fā)溫度和冷卻水溫升的數(shù)學(xué)關(guān)系式和相應(yīng)的修正系數(shù),結(jié)果表明這些計(jì)算式精確性比較高。WANG等[10-11]建立了循環(huán)凈功的理論模型,推導(dǎo)理論計(jì)算公式,計(jì)算結(jié)果與ORC的數(shù)值計(jì)算結(jié)果也非常一致。HE等[12-13]建立了循環(huán)凈功與優(yōu)化蒸發(fā)溫度的解析表達(dá)式,表明優(yōu)化蒸發(fā)溫度的理論值與ORC數(shù)值計(jì)算之間的偏差僅為–0.86%~2.3%,平均偏差小于1%。CHEN等[14]綜述了ORC和超臨界朗肯循環(huán)以及工質(zhì)篩選標(biāo)準(zhǔn),討論了工質(zhì)汽化潛熱、密度、液態(tài)比熱容對(duì)循環(huán)的影響和關(guān)于過(guò)熱的經(jīng)濟(jì)性,推導(dǎo)出工質(zhì)在膨脹機(jī)等熵膨脹過(guò)程中比焓變的關(guān)系式。另外,MAIZZA等[15]提出汽化潛熱越大、密度越高、液態(tài)比熱容越小的工質(zhì)越適合于ORC,而YAMAMOTO等[16]認(rèn)為低汽化潛熱的工質(zhì)更好,因?yàn)榕蛎洐C(jī)進(jìn)口處工質(zhì)為飽和狀態(tài)時(shí)才達(dá)到優(yōu)化工況。LIU等[17]提出工質(zhì)的飽和氣線斜率公式,但未分析飽和氣線斜率對(duì)循環(huán)性能的影響。

    綜上所述,常規(guī)有機(jī)朗肯循環(huán)一般是針對(duì)實(shí)際工質(zhì)的數(shù)值計(jì)算、實(shí)驗(yàn)研究或應(yīng)用,由于實(shí)際工質(zhì)的種類有限,且工質(zhì)物性是定量、非變量,使得有機(jī)朗肯循環(huán)的研究結(jié)果或結(jié)論往往是具體、局限的,不具有普遍性?;谇捌谘芯刻岢龅哪M有機(jī)朗肯循環(huán)的梯形循環(huán)及其理論模型[18],本文建立工質(zhì)飽和線形狀(飽和液線斜率、飽和氣線斜率與臨界溫度)及物性為自變量的數(shù)學(xué)關(guān)系,即將實(shí)際工質(zhì)的物性變量化,應(yīng)用梯形循環(huán)及其理論模型,進(jìn)行梯形循環(huán)(或有機(jī)朗肯循環(huán))的理論研究,從理論角度研究有機(jī)朗肯循環(huán)的普遍性熱力學(xué)規(guī)律。重點(diǎn)針對(duì)有機(jī)朗肯循環(huán)中有極值功向無(wú)極值功轉(zhuǎn)化的凈功轉(zhuǎn)折線,所對(duì)應(yīng)的熱源轉(zhuǎn)折溫度及其最高優(yōu)化工況等特征進(jìn)行研究。

    1 工質(zhì)飽和線的數(shù)學(xué)模型

    前期研究提出并建立了梯形循環(huán)(trapezoidal cycle,TPC)及其理論模型[18]。梯形循環(huán)是將基本的內(nèi)可逆ORC簡(jiǎn)化為T-s圖中的梯形循環(huán),如圖1所示。其中將ORC中的兩個(gè)過(guò)程作如下簡(jiǎn)化:①用過(guò)程3–4代替ORC中的過(guò)程3–3′–4,這樣造成了液相區(qū)比功“Δw1”的減少(圖中3–3′–4–3圍成的面積);②采用過(guò)程2′–2–2′代替ORC中干膨脹后的過(guò)熱冷卻過(guò)程2′–2′,這樣,對(duì)于干工質(zhì)或“干膨脹”,增加了比功“Δwsup.h”(圖1中2′–2–2′–2′圍成的面積),但對(duì)濕工質(zhì)與等熵工質(zhì)或“濕膨脹”,則無(wú)此功增量。計(jì)算表明[18]:TPC與ORC之間的偏差為:功率偏差為0.3%~–3.3%,熱效率偏差為0.59%~–2.94%,相對(duì)偏差的絕對(duì)值均小于5%,表明梯形循環(huán)完全可以模擬有機(jī)朗肯循環(huán)。

    圖1 有機(jī)朗肯循環(huán)與梯形循環(huán)

    基于梯形循環(huán)及其理論模型,建立工質(zhì)飽和線的數(shù)學(xué)模型。如圖2(a)、(b)、(c),線段4–5為工質(zhì)飽和液線上蒸發(fā)溫度處的切線,斜率為kl,線段1–5為工質(zhì)飽和氣線上蒸發(fā)溫度處的切線,斜率為kg;并定義kl與kg之比為斜率比rk,rk=kl/kg。5點(diǎn)溫度定義為斜率切線的交點(diǎn)溫度Tcr',其值可反映工質(zhì)臨界溫度的高低。6點(diǎn)溫度定義為頂點(diǎn)溫度Tm。Tcr'與Tm之間的差值為ΔTm。

    (1)對(duì)濕工質(zhì),如圖2(a),由三角形相似得式(1) 。

    式中,Δsr為蒸發(fā)溫度Te下兩相區(qū)的熵變。

    (2)對(duì)干工質(zhì),如圖2(b),由三角形相似得式(3)。

    (3)對(duì)于等熵工質(zhì),Tcr'=Te+klΔsr,即公式(2)或公式(4)中kg→∞。

    因此,干、濕及等熵工質(zhì)的切線交點(diǎn)溫度Tcr'表達(dá)式一致,公式(2)或公式(4)為通式。

    圖2 頂點(diǎn)溫度與飽和線斜率及斜率切線交點(diǎn)溫度的關(guān)系

    由此可推導(dǎo)出工質(zhì)的頂點(diǎn)溫度Tm表達(dá)式如式(5)。

    2 工質(zhì)熱源轉(zhuǎn)折溫度及其經(jīng)驗(yàn)公式與偏差

    研究表明:工質(zhì)循環(huán)的優(yōu)化工況(即凈輸出功極大值)只存在于蒸發(fā)器傳熱窄點(diǎn)位于工質(zhì)泡點(diǎn)處的情況[18]?;谔菪文P停瑐鳠嵴c(diǎn)位于工質(zhì)泡點(diǎn)時(shí)的凈輸出功W表達(dá)式為式(6)。

    如圖3所示,熱源轉(zhuǎn)折溫度Th,shift所對(duì)應(yīng)的凈功為有極值功向無(wú)極值功的轉(zhuǎn)折;當(dāng)熱源溫度低于該Th,shift時(shí),工質(zhì)循環(huán)存在優(yōu)化工況,即存在凈功極大值與對(duì)應(yīng)的優(yōu)化蒸發(fā)溫度;當(dāng)熱源溫度高于該Th,shift時(shí),工質(zhì)循環(huán)凈功不再具有極大值,而是隨蒸發(fā)溫度單調(diào)上升。表明亞臨界有機(jī)朗肯循環(huán)中存在著有極值功循環(huán)與無(wú)極值功循環(huán),或功率由有極值向無(wú)極值的發(fā)展,如圖4所示。明顯地,熱源轉(zhuǎn)折溫度下的優(yōu)化工況為優(yōu)化工況的最高工況,即熱源轉(zhuǎn)折溫度Th,shift所對(duì)應(yīng)的優(yōu)化蒸發(fā)溫度Te,opt_M與凈輸出功極值Wmax_M是優(yōu)化工況的最高值。

    研究表明,除凈功的轉(zhuǎn)折線與對(duì)應(yīng)的工質(zhì)熱源轉(zhuǎn)折溫度,還存在另外兩條凈功特征線與對(duì)應(yīng)的工質(zhì)熱源特征溫度[19],如圖3所示:①熱源轉(zhuǎn)變點(diǎn)溫度Th,turn,蒸發(fā)溫度趨近于工質(zhì)臨界溫度時(shí),凈功由一直下降出現(xiàn)向上增加趨勢(shì),對(duì)應(yīng)的熱源溫度定義為熱源轉(zhuǎn)變點(diǎn)溫度Th,turn;②熱源上限溫度Th,up,蒸發(fā)溫度趨近于工質(zhì)臨界溫度時(shí),凈功由一直上升出現(xiàn)下降趨勢(shì),對(duì)應(yīng)的熱源溫度定義為熱源上限溫度Th,up。

    表1給出了21種工質(zhì)在蒸發(fā)器傳熱窄點(diǎn)溫差5℃和冷凝溫度35℃下,有機(jī)朗肯循環(huán)中工質(zhì)熱源轉(zhuǎn)折溫度,其中文獻(xiàn)[19]的計(jì)算精確度為0.1℃,本文為0.01℃。表2給出了R227ea的熱源轉(zhuǎn)折溫度隨窄點(diǎn)溫差(0~20℃)的變化。結(jié)果表明:①工質(zhì)熱源轉(zhuǎn)折溫度隨蒸發(fā)器傳熱窄點(diǎn)溫差成近似的線性關(guān)系;②熱源轉(zhuǎn)折溫度所對(duì)應(yīng)的最高優(yōu)化工況:最高優(yōu)化蒸發(fā)溫度Te,opt_M與最大凈輸功極值Wmax_M不隨蒸發(fā)器窄點(diǎn)溫差而變化。

    圖3 工質(zhì)R227ea循環(huán)凈功特征線與對(duì)應(yīng)的熱源水特征溫度

    圖4 一定熱源水溫度,不同臨界溫度工質(zhì)的循環(huán)凈功比較

    表1 工質(zhì)的熱源轉(zhuǎn)折溫度(ΔTp=5℃,Tc=35℃)

    表2 R227ea在不同窄點(diǎn)溫差下的熱源轉(zhuǎn)折溫度及其最高優(yōu)化工況值

    通過(guò)對(duì)大量工質(zhì)熱源轉(zhuǎn)折溫度的計(jì)算與擬合,可推導(dǎo)出有機(jī)朗肯循環(huán)中工質(zhì)熱源轉(zhuǎn)折溫度的經(jīng)驗(yàn)公式,如式(9)。

    表1給出了熱源轉(zhuǎn)折溫度經(jīng)驗(yàn)公式與有機(jī)朗肯循環(huán)數(shù)值計(jì)算之間的偏差,二者的偏差較小,其中最大偏差為R600的1.97%。表明該經(jīng)驗(yàn)公式可用于有機(jī)朗肯循環(huán)工質(zhì)熱源轉(zhuǎn)折溫度的模擬計(jì)算。

    3 線性飽和線工質(zhì)模型與熱源轉(zhuǎn)折溫度及其理論公式

    為簡(jiǎn)化分析,將工質(zhì)的飽和液線與飽和氣線簡(jiǎn)化為直線,構(gòu)成線性飽和線工質(zhì),如圖5。這樣,線性飽和線工質(zhì)的飽和液線斜率kl與飽和氣線斜率kg不隨溫度而變化,為常量。

    圖5 線性飽和線工質(zhì)模型

    此時(shí),線性飽和線工質(zhì)的臨界溫度Tcr即為交點(diǎn)溫度Tcr'。所以,公式(2)可寫為式(10)。

    蒸發(fā)溫度處的潛熱r=Te×Δsr,代入工質(zhì)流量得到式(11)。

    當(dāng)蒸發(fā)溫度Te達(dá)到頂點(diǎn),即臨界溫度Tcr處,則為工質(zhì)的熱源轉(zhuǎn)折溫度(Th,shift_id),如式(14)。

    式(14)表明:工質(zhì)的熱源轉(zhuǎn)折溫度與工質(zhì)的臨界溫度及蒸發(fā)器窄點(diǎn)溫差分別成線性關(guān)系,與前述的ORC分析結(jié)論相一致。

    因此,由式(15)計(jì)算出的優(yōu)化工況則為最高優(yōu)化工況,即最高優(yōu)化蒸發(fā)溫度與最大凈功極值。

    但是由于mw并不是Te的一元函數(shù),還是rk和cp的函數(shù),所以式(14)會(huì)有偏差。取線性飽和線工質(zhì)的臨界溫度Tcr范圍90~200℃,斜率比rk為–2~0.35,對(duì)公式(14)的偏差進(jìn)行計(jì)算如下。

    當(dāng)rk=0(即等熵工質(zhì))時(shí),偏差范圍為0~00.15℃,可以不用修正,仍采用式(14)。

    當(dāng)rk>0(即干工質(zhì))時(shí),偏差范圍為0.15~– 0.05℃,可以不用修正,仍采用式(14)。

    當(dāng)rk<0(即濕工質(zhì))時(shí),修正為式(16)。

    表3給出了臨界溫度Tcr150℃下,線性飽和線工質(zhì)熱源轉(zhuǎn)折溫度Th,shift_id的梯形循環(huán)數(shù)值計(jì)算與理論公式(14)或修正式(16)的計(jì)算結(jié)果。大量計(jì)算表明:在斜率比rk為–2~0.35范圍內(nèi),線性飽和線工質(zhì)熱源轉(zhuǎn)折溫度的梯形循環(huán)數(shù)值計(jì)算與上述理論公式(14)或修正式(16)之間的最大偏差分別為:Tcr=100℃時(shí)為0.25%,Tcr=150℃時(shí)為0.24%,Tcr=200℃時(shí)為0.68%。

    表3 線性飽和線工質(zhì)熱源轉(zhuǎn)折溫度下的優(yōu)化工況與斜率比的關(guān)系(Tcr=150℃;ΔTp=5℃)

    4 線性飽和線工質(zhì)熱源轉(zhuǎn)折溫度下的優(yōu)化工況及其特征

    由前述分析,工質(zhì)熱源轉(zhuǎn)折溫度Th,shifit_id下的優(yōu)化工況為最高優(yōu)化工況,即最高優(yōu)化蒸發(fā)溫度Te,opt_M與最大凈功極值Wmax_M。

    式(15)工質(zhì)熱源轉(zhuǎn)折溫度下凈功計(jì)算式中無(wú)窄點(diǎn)溫差 ΔTp,表明熱源轉(zhuǎn)折溫度下的凈功與窄點(diǎn)溫差無(wú)關(guān),由此可推論:熱源轉(zhuǎn)折溫度下的最高優(yōu)化工況參數(shù)(Te,opt_M、Wmax_M)也與窄點(diǎn)溫差無(wú)關(guān),與前述ORC的分析結(jié)論相一致。

    下面具體計(jì)算與分析熱源轉(zhuǎn)折溫度下的優(yōu)化工況隨斜率比的變化。假設(shè)計(jì)算條件為:熱源水進(jìn)口溫度為工質(zhì)熱源轉(zhuǎn)折溫度Th,shifit_id,熱源水流量mh為1kg/s,冷凝溫度Tc為35℃,蒸發(fā)器窄點(diǎn)溫差ΔTp為5℃;線性飽和線工質(zhì)的臨界溫度Tcr設(shè)為150℃,斜率比rk為–2~0.35,計(jì)算結(jié)果見表3。表明最高優(yōu)化工況隨斜率比rk單調(diào)增加,至干工質(zhì)時(shí),最高優(yōu)化工況達(dá)到工質(zhì)的臨界點(diǎn)。由于蒸發(fā)溫度為Tcr時(shí),熱源進(jìn)出口已無(wú)溫差,可取比Tcr稍低的溫度,如Tcr–0.5進(jìn)行計(jì)算。

    5 結(jié)論

    有機(jī)朗肯循環(huán)中存在有極值功循環(huán)與無(wú)極值功循環(huán),即存在一種凈功的轉(zhuǎn)折線與所對(duì)應(yīng)的工質(zhì)熱源轉(zhuǎn)折溫度,它對(duì)工質(zhì)的選擇與循環(huán)性能評(píng)價(jià)具有指標(biāo)性的作用。基于梯形循環(huán)及其理論模型,本文建立工質(zhì)飽和線形狀(飽和液線斜率、飽和氣線斜率及其斜率比與臨界溫度)及物性為自變量的數(shù)學(xué)關(guān)系,并提出和定義了線性飽和線工質(zhì)模型,建立循環(huán)性能與工質(zhì)物性之間的數(shù)學(xué)關(guān)系。分別推導(dǎo)出有機(jī)朗肯循環(huán)中工質(zhì)熱源轉(zhuǎn)折溫度的經(jīng)驗(yàn)公式,以及梯形循環(huán)中線性飽和線工質(zhì)的熱源轉(zhuǎn)折溫度理論公式及其修正。研究熱源轉(zhuǎn)折溫度及其對(duì)應(yīng)的優(yōu)化工況的特征與一般性規(guī)律。得到如下結(jié)論。

    (1)計(jì)算表明,21種工質(zhì)熱源轉(zhuǎn)折溫度的有機(jī)朗肯循環(huán)數(shù)值計(jì)算與經(jīng)驗(yàn)公式之間的偏差小于1.97%。

    (2)有機(jī)朗肯循環(huán)中,工質(zhì)熱源轉(zhuǎn)折溫度與蒸發(fā)器窄點(diǎn)溫差成線性關(guān)系,熱源轉(zhuǎn)折溫度所對(duì)應(yīng)的最高優(yōu)化工況(最高優(yōu)化蒸發(fā)溫度與最大凈功極值)不隨蒸發(fā)器窄點(diǎn)溫差而變化。

    (3)線性飽和線工質(zhì)模型下,推導(dǎo)出工質(zhì)熱源轉(zhuǎn)折溫度Th,shift_id的理論公式,為臨界溫度與蒸發(fā)器窄點(diǎn)溫差之和,其中濕工質(zhì)要進(jìn)行一定的修正。Th,shift_id下的最高優(yōu)化工況(最高優(yōu)化蒸發(fā)溫度Te,opt_M與最大凈功極值Wmax_M)隨斜率比單調(diào)增加,至干工質(zhì)時(shí),最高優(yōu)化工況接近工質(zhì)的臨界點(diǎn)。

    符號(hào)說(shuō)明

    c—— 比熱容,J/(kg·K)

    Ja—— 雅克比數(shù),J/J

    k—— 斜率,K2·kg/J

    m—— 質(zhì)量流量,kg/s

    q—— 比吸熱量,J/kg

    r—— 斜率比,kl/kg,量綱為1

    s—— 比熵,J/(kg·K)

    T—— 溫度,K或℃

    W—— 凈輸出功,J

    w—— 比功,J/kg

    η—— 熱效率,%

    上、下角標(biāo)

    c —— 冷凝

    cr —— 臨界

    e —— 蒸發(fā)

    h —— 熱源

    i —— 進(jìn)口

    m —— 頂點(diǎn)

    max —— 極大值

    opt —— 優(yōu)化

    ORC —— 有機(jī)朗肯循環(huán)

    p —— 壓力

    shift —— 轉(zhuǎn)折

    TPC—— 梯形循環(huán)

    w—— 工質(zhì)

    [1] 徐榮吉,張曉暉,閆美玉,等.過(guò)熱/過(guò)冷對(duì)內(nèi)回?zé)嵊袡C(jī)朗肯循環(huán)影響的熱力學(xué)分析[J].化工進(jìn)展,2016,35(12):3783-3792.XU Rongji,ZHANG Xiaohui,YAN Meiyu,et al.Thermal dynamics analysis on the organic Rankine cycle (ORC)with internal heat regenerator at superheated and subcooling conditions[J].Chemical Industry and Engineering Progress,2016,35(12):3783-3792.

    [2] 韓中合,潘歌,范偉,等.內(nèi)回?zé)崞鲗?duì)低溫有機(jī)朗肯循環(huán)熱力性能的影響及工質(zhì)選擇[J].化工進(jìn)展,2016,35(1):40-47.HAN Zhonghe,PAN Ge,F(xiàn)AN Wei,et al.Effect of internal heat exchanger on thermodynamic performance of low temperature organic Rankine cycle and working fluid selection[J].Chemical Industry and Engineering Progress,2016,35(1):40-47.

    [3] 楊緒飛,戚風(fēng)亮,劉秀龍,等.有機(jī)朗肯循環(huán)膨脹機(jī)入口過(guò)熱度實(shí)驗(yàn)[J].化工進(jìn)展,2016,35(7):2007-2014.YANG Xufei,QI Fengliang,LIU Xiulong,et al.Experiment on expander inlet superheat of organic Rankine cycle[J].Chemical Industry and Engineering Progress,2016,35(7):2007-2014.

    [4] 王明濤,方箏,劉啟一.渦輪增壓柴油機(jī)余熱利用的有機(jī)郎肯循環(huán)烴類高溫工質(zhì)熱力學(xué)分析[J].化工進(jìn)展,2016,35(9):2721-2727.WANG Mingtao,F(xiàn)ANG Zheng,LIU Qiyi.Thermodynamic analysis of an organic Rankine cycle for waste heat recovery of a turbo-charged diesel engine based on working fluids of alkanes[J].Chemical Industry and Engineering Progress,2016,35(9):2721-2727.

    [5] 張軍輝,劉娟芳,陳清華.有機(jī)朗肯循環(huán)系統(tǒng)最佳蒸發(fā)溫度和?分析[J].化工學(xué)報(bào),2013,64(3):820-826.ZHANG Junhui,LIU Fangjuan,CHEN Qinghua.Optimal evaporating temperature and exergy analysis for original Rankine cycle[J].CIESC Journal,2013,64(3):820-826.

    [6] 賀超,劉朝,高虹,等.工業(yè)余熱有機(jī)朗肯循環(huán)參數(shù)優(yōu)化和分析[J].工程熱物理學(xué)報(bào),2012,33(12):2042-2046.HE Chao,LIU Chao,GAO Hong,et al.Parameters optimization and performance analysis of organic Rankine cycle for industrial waste heat recovery[J].Journal of Engineering Thermophysics,2012,33(12):2042-2046.

    [7] LI H,MA X L,WEI X L,et al.Selection of working fluids for low-temperature waste heat recovery using organic Rankine cycle[J].Advanced Materials Research,2012,512-515:1217-1222.

    [8] YU H S,F(xiàn)ENG X,WANG Y F.A new pinch based method for simultaneous selection of working fluid and operating conditions in an ORC (Organic Rankine Cycle) recovering waste heat[J].Energy,2015,90:36-46.

    [9] 嚴(yán)家祿.低溫?zé)崮馨l(fā)電方案中選擇工質(zhì)和確定參數(shù)的熱力學(xué)原則和計(jì)算[J].工程熱物理學(xué)報(bào),1982,3(1):1-7.YAN Jialu.Thermodynamic principles and formulas for choosing working fluids and parameters in desgning power plant of low temperature heat[J].Journal of Engineering Thermophysics,1982,3(1):1-7.

    [10] WANG D X,LING X,PENG H,et al.Efficiency and optimal performance evaluation of organic Rankine cycle for low grade waste heat power generation[J].Energy,2013,50(1):343-352.

    [11] WANG D X,LING X,PENG H.Cost–effectiveness performance analysis of organic Rankine cycle for low grade heat utilization coupling with operation condition[J].Applied Thermal Engineering,2013,58(1-2):571-584.

    [12] HE C,LIU C,GAO H,et al.The optimal evaporation temperature and working fluids for subcritical organic Rankine cycle[J].Energy,2012,38(1):136-143.

    [13] HE C,LIU C,ZHOU M T,et al.A new selection principle of working fluids for subcritical organic Rankine cycle coupling with different heat sources[J].Energy,2014,68(15):283-291.

    [14] CHEN H,GOSWAMI D Y,STEFANAKOS E K.A review of thermodynamic cycles and working fluids for the conversion of low-grade heat[J].Renewable & Sustainable Energy Reviews,2010,14(9):3059-3067.

    [15] MAIZZA V,MAIZZA A.Working fluids in non-steady flows for waste energy recovery systems[J].Applied Thermal Engineering,1996,16:579-90.

    [16] YAMAMOTO T,F(xiàn)URUHATA T,ARAI N,et al.Design and testing of the organic Rankine cycle[J].Energy,2001,26(3):239-251.

    [17] LIU B T,CHIEN K H,WANG C C.Effect of working fluids on organic Rankine cycle for waste heat recovery[J]. Energy,2004,29(8):1207-1217.

    [18] LI X G.A trapezoidal cycle with theoretical model based on organic Rankine cycle[J].International Journal of Energy Research,2016,40(12):1624-1637.

    [19] LI X G,ZHAO W J,LIN D D,et al.Working fluid selection based on critical temperature and water temperature in organic Rankine cycle[J].Science China Technological Sciences,2015,58(1):138-146.

    Shift-temperature of heating fluid and its characteristics for working fluid in organic Rankine cycle

    LI Xinguo,WANG Wei,ZHAI Zhe,NAN Laifu,WANG Jingyi
    (School of Mechanical Engineering,Tianjin University,Tianjin 300350,China)

    The existence of the maximum power cycle and no maximum power cycle in organic Rankine cycle(ORC)means that there is a kind of shift-curve of net power output and its corresponding shift-temperature of heating fluid for the working fluids; which played an index function in the selection of the working fluids and the evaluation of the cycle performance. Based on trapezoidal cycle and its theoretical model,the relation and formulas of the saturated line shape(saturated liquid line slope,saturated gas line slope and its ratio and critical temperature)and thermal properties of the working fluids as variables were established. The model of working fluid with linear saturated line was proposed and defined. The mathematical model and relation between the cycle performance and properties of the working fluids were built. The empirical equations of the shift-temperature of heating fluid for the working fluid in ORC and the theoretical equations of the shift-temperature of heating fluid in working fluid with linear saturated line were derived,respectively. The characteristics and general principles of the shift-temperature of heating fluid and its corresponding optimum condition were investigated. Results showed that the deviation of shift-temperature of heating fluid with 21 working fluids between empirical equations and numerical simulations in ORC was less than 1.97%.The shift-temperature of heating fluid in model of working fluid with linear saturated line is the additionof the critical temperature of working fluid and the pinch point temperature difference of evaporator.The maximum optimal conditions(maximum optimal evaporation temperature and maximum net power output)corresponding to the shift-temperature of heating fluid increases with the ratio of the saturated line slope, which will be close to the critical point of the working fluid when it becomes the dry working fluid.

    organic Rankine cycle;trapezoidal cycle;working fluid with linear saturated line;shift-temperature of heating fluid for working fluid;maximum optimum condition

    TK123

    :A

    :1000-6613(2017)09-3223-08

    10.16085/j.issn.1000-6613.2017-0056

    2017-01-10;修改稿日期:2017-05-08。

    國(guó)家自然科學(xué)基金項(xiàng)目(51276122)。

    及聯(lián)系人:李新國(guó)(1965—),男,博士,教授,主要從事工程熱力學(xué)研究。E-mail:xgli@tju.edu.cn。

    猜你喜歡
    臨界溫度工質(zhì)熱源
    海洋溫差能發(fā)電熱力循環(huán)系統(tǒng)的工質(zhì)優(yōu)選
    橫流熱源塔換熱性能研究
    煤氣與熱力(2021年3期)2021-06-09 06:16:20
    Bogoliubov-Tolmachev-Shirkov模型臨界溫度和能隙解的數(shù)值方法
    采用R1234ze(E)/R245fa的非共沸混合工質(zhì)有機(jī)朗肯循環(huán)系統(tǒng)實(shí)驗(yàn)研究
    采用二元非共沸工質(zhì)的有機(jī)朗肯循環(huán)熱力學(xué)分析
    基于啟發(fā)式動(dòng)態(tài)規(guī)劃的冷熱源優(yōu)化控制
    電子制作(2017年19期)2017-02-02 07:08:31
    若干低GWP 純工質(zhì)在空調(diào)系統(tǒng)上的應(yīng)用分析
    中部槽激光-MAG復(fù)合熱源打底焊焊接工藝研究
    焊接(2015年8期)2015-07-18 10:59:13
    多類型熱源點(diǎn)共存下的區(qū)域熱力統(tǒng)籌供應(yīng)探討
    RDX基炸藥熱起爆臨界溫度的測(cè)試及數(shù)值計(jì)算
    久久热精品热| 国产aⅴ精品一区二区三区波| 一个人免费在线观看电影| 亚洲国产精品999在线| 激情在线观看视频在线高清| 直男gayav资源| 亚洲av熟女| 国产一区二区激情短视频| 久久精品国产自在天天线| 在线观看免费视频日本深夜| 淫妇啪啪啪对白视频| 日韩欧美在线二视频| 国产单亲对白刺激| 51国产日韩欧美| 深爱激情五月婷婷| 精品国产三级普通话版| 九九在线视频观看精品| 国产午夜精品论理片| 99精品在免费线老司机午夜| 一级av片app| 一边摸一边抽搐一进一小说| 日本免费a在线| 我要看日韩黄色一级片| www.999成人在线观看| 国产伦精品一区二区三区视频9| 午夜福利高清视频| 欧美三级亚洲精品| 成人高潮视频无遮挡免费网站| 人人妻人人看人人澡| 精品久久久久久成人av| 午夜激情欧美在线| 午夜福利视频1000在线观看| 日本 欧美在线| 91久久精品国产一区二区成人| 欧美日韩中文字幕国产精品一区二区三区| www.色视频.com| 人妻夜夜爽99麻豆av| 国内精品美女久久久久久| 欧美精品啪啪一区二区三区| 国产伦精品一区二区三区四那| 天堂影院成人在线观看| 又爽又黄无遮挡网站| 亚洲成人久久爱视频| 午夜免费激情av| 国产精品女同一区二区软件 | 老女人水多毛片| 亚洲av不卡在线观看| 久久久国产成人免费| 男女视频在线观看网站免费| 深爱激情五月婷婷| 黄色视频,在线免费观看| 欧美黄色片欧美黄色片| 国产午夜福利久久久久久| 国产精品亚洲美女久久久| 国产综合懂色| 成熟少妇高潮喷水视频| 国产人妻一区二区三区在| 亚洲av免费高清在线观看| 男女视频在线观看网站免费| 搡老妇女老女人老熟妇| 日韩欧美三级三区| 99久久精品热视频| 国产一区二区在线观看日韩| 国内揄拍国产精品人妻在线| 免费看a级黄色片| 嫩草影视91久久| 亚洲精品在线观看二区| 两个人的视频大全免费| 波野结衣二区三区在线| 久久久久久久久久成人| 免费av观看视频| 欧美日本视频| 最近最新中文字幕大全电影3| 日韩欧美在线二视频| 亚洲欧美精品综合久久99| 国产成人欧美在线观看| 午夜激情欧美在线| 可以在线观看的亚洲视频| 别揉我奶头~嗯~啊~动态视频| 毛片女人毛片| 欧美色视频一区免费| 亚洲中文字幕日韩| 老熟妇乱子伦视频在线观看| 老司机深夜福利视频在线观看| 国产精品99久久久久久久久| 亚洲天堂国产精品一区在线| 欧洲精品卡2卡3卡4卡5卡区| 精品99又大又爽又粗少妇毛片 | 亚洲aⅴ乱码一区二区在线播放| 欧美在线一区亚洲| 久久亚洲真实| 亚洲18禁久久av| 亚洲精品色激情综合| 欧美午夜高清在线| 国产伦精品一区二区三区视频9| 色在线成人网| 精华霜和精华液先用哪个| 丰满的人妻完整版| 亚洲综合色惰| 久久精品国产自在天天线| 美女免费视频网站| 搡老岳熟女国产| 美女大奶头视频| 精品久久久久久,| 男女做爰动态图高潮gif福利片| 深夜精品福利| 黄片小视频在线播放| 国产欧美日韩一区二区三| 国产精品av视频在线免费观看| 一区福利在线观看| 色在线成人网| 精品久久久久久,| 国内揄拍国产精品人妻在线| 国产91精品成人一区二区三区| 免费av毛片视频| 国产高清视频在线观看网站| 国产精品一及| 男女床上黄色一级片免费看| 在线免费观看的www视频| 欧美精品国产亚洲| 99在线人妻在线中文字幕| 国产av不卡久久| 久久久久久大精品| 国产精品日韩av在线免费观看| 丰满乱子伦码专区| 国产三级黄色录像| 免费看美女性在线毛片视频| 毛片一级片免费看久久久久 | 久久国产乱子免费精品| 国产真实伦视频高清在线观看 | 自拍偷自拍亚洲精品老妇| 美女大奶头视频| 欧美bdsm另类| 国产精品久久久久久精品电影| av在线老鸭窝| 精品无人区乱码1区二区| 女人被狂操c到高潮| 欧美极品一区二区三区四区| 欧美成人免费av一区二区三区| av黄色大香蕉| 又黄又爽又刺激的免费视频.| 精品久久久久久久久久免费视频| 啦啦啦观看免费观看视频高清| 国产欧美日韩精品一区二区| 看片在线看免费视频| 99国产精品一区二区三区| 国产乱人伦免费视频| 亚洲专区国产一区二区| 国产亚洲精品久久久久久毛片| 九九在线视频观看精品| 亚洲专区中文字幕在线| 日韩中字成人| av欧美777| 免费电影在线观看免费观看| 国产精品一区二区三区四区久久| av在线天堂中文字幕| 嫩草影院入口| 久9热在线精品视频| xxxwww97欧美| 久久久久精品国产欧美久久久| 国产色爽女视频免费观看| 国产aⅴ精品一区二区三区波| av专区在线播放| 色5月婷婷丁香| 黄色日韩在线| 亚洲中文字幕一区二区三区有码在线看| 成年女人看的毛片在线观看| 亚洲综合色惰| 久久性视频一级片| 国产av一区在线观看免费| 变态另类成人亚洲欧美熟女| 亚洲精品成人久久久久久| 啪啪无遮挡十八禁网站| 欧美极品一区二区三区四区| 内射极品少妇av片p| 亚洲国产精品sss在线观看| 国内毛片毛片毛片毛片毛片| 亚洲五月天丁香| 一边摸一边抽搐一进一小说| 欧美成狂野欧美在线观看| 精品熟女少妇八av免费久了| 在线免费观看不下载黄p国产 | 欧美色视频一区免费| 国产69精品久久久久777片| 亚洲久久久久久中文字幕| 中文字幕人妻熟人妻熟丝袜美| 看片在线看免费视频| 亚洲欧美日韩东京热| 国产免费一级a男人的天堂| 亚洲精品亚洲一区二区| 亚洲成人中文字幕在线播放| 亚洲av美国av| 观看美女的网站| 国产精华一区二区三区| 日韩 亚洲 欧美在线| 一卡2卡三卡四卡精品乱码亚洲| 免费看日本二区| 99久久成人亚洲精品观看| 日日夜夜操网爽| 一个人看视频在线观看www免费| 久久精品91蜜桃| 最好的美女福利视频网| www日本黄色视频网| 一二三四社区在线视频社区8| 超碰av人人做人人爽久久| 亚洲黑人精品在线| 欧美成狂野欧美在线观看| 亚洲欧美精品综合久久99| 久久伊人香网站| 免费一级毛片在线播放高清视频| 一个人看视频在线观看www免费| 一级黄色大片毛片| 日韩免费av在线播放| 成人三级黄色视频| 宅男免费午夜| 噜噜噜噜噜久久久久久91| 久久久久久久精品吃奶| 给我免费播放毛片高清在线观看| 最新在线观看一区二区三区| 亚洲avbb在线观看| 中文字幕高清在线视频| 精品福利观看| 精品一区二区三区人妻视频| 亚洲成人中文字幕在线播放| 国产精品一及| 色综合站精品国产| 亚洲av熟女| 男插女下体视频免费在线播放| 国产美女午夜福利| av黄色大香蕉| 日韩大尺度精品在线看网址| 国产熟女xx| 国产老妇女一区| 欧美日本亚洲视频在线播放| 91麻豆精品激情在线观看国产| 国产精品女同一区二区软件 | 国产精品自产拍在线观看55亚洲| 久久久久国内视频| 亚洲一区高清亚洲精品| 国产精品人妻久久久久久| 国产亚洲精品久久久久久毛片| 日本五十路高清| 亚洲精品色激情综合| 国产精品女同一区二区软件 | 一边摸一边抽搐一进一小说| 亚洲人成伊人成综合网2020| 99精品在免费线老司机午夜| 在线观看66精品国产| 2021天堂中文幕一二区在线观| 波多野结衣高清作品| 岛国在线免费视频观看| 别揉我奶头~嗯~啊~动态视频| 又爽又黄a免费视频| 国产精品99久久久久久久久| 午夜福利在线观看免费完整高清在 | 一本久久中文字幕| 人妻久久中文字幕网| 啦啦啦韩国在线观看视频| 看十八女毛片水多多多| 老熟妇仑乱视频hdxx| 五月伊人婷婷丁香| 久久久久久大精品| 一本一本综合久久| 中文在线观看免费www的网站| 在线十欧美十亚洲十日本专区| 日韩欧美精品v在线| 中文字幕熟女人妻在线| 欧美日韩福利视频一区二区| 国产av一区在线观看免费| 色综合婷婷激情| 久久热精品热| 免费电影在线观看免费观看| 亚洲aⅴ乱码一区二区在线播放| 久久99热6这里只有精品| 午夜日韩欧美国产| 日韩欧美 国产精品| 高清日韩中文字幕在线| 久久久精品欧美日韩精品| 欧美成狂野欧美在线观看| 日韩欧美在线二视频| 99久久99久久久精品蜜桃| 国产乱人视频| 99国产极品粉嫩在线观看| 最新中文字幕久久久久| 如何舔出高潮| 怎么达到女性高潮| 在线观看av片永久免费下载| 欧美日韩中文字幕国产精品一区二区三区| 国产亚洲av嫩草精品影院| 又爽又黄a免费视频| 国产黄片美女视频| 日本 av在线| 性色av乱码一区二区三区2| 99在线视频只有这里精品首页| 亚洲不卡免费看| 亚洲欧美日韩东京热| 午夜免费成人在线视频| 久久精品影院6| 久久久色成人| 亚洲在线自拍视频| 网址你懂的国产日韩在线| 一级a爱片免费观看的视频| 黄色视频,在线免费观看| 51国产日韩欧美| 久久久久久久亚洲中文字幕 | 国产黄a三级三级三级人| 欧美xxxx黑人xx丫x性爽| 日本在线视频免费播放| 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产精品sss在线观看| 一级毛片久久久久久久久女| 国产成人欧美在线观看| 亚洲一区二区三区不卡视频| av在线观看视频网站免费| 亚洲av电影在线进入| 三级毛片av免费| 午夜视频国产福利| 此物有八面人人有两片| 露出奶头的视频| 中文字幕精品亚洲无线码一区| 99热精品在线国产| 亚洲中文字幕日韩| 国产私拍福利视频在线观看| 日本 欧美在线| 亚洲天堂国产精品一区在线| 99视频精品全部免费 在线| 色吧在线观看| 嫩草影院入口| 欧美高清性xxxxhd video| 成人高潮视频无遮挡免费网站| 欧美bdsm另类| av在线天堂中文字幕| 美女 人体艺术 gogo| 国产在线男女| 又粗又爽又猛毛片免费看| 欧美日韩中文字幕国产精品一区二区三区| 能在线免费观看的黄片| 中文字幕高清在线视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲片人在线观看| 久久久精品欧美日韩精品| 免费在线观看成人毛片| 午夜免费男女啪啪视频观看 | 久久久久久久久大av| 在线观看美女被高潮喷水网站 | 日韩亚洲欧美综合| 又爽又黄a免费视频| 成人高潮视频无遮挡免费网站| 亚洲 国产 在线| 男女床上黄色一级片免费看| 亚洲精华国产精华精| 亚洲男人的天堂狠狠| 婷婷亚洲欧美| 国产一区二区在线观看日韩| 51午夜福利影视在线观看| 午夜激情福利司机影院| 国产伦精品一区二区三区四那| 一区二区三区激情视频| 亚洲激情在线av| 18禁裸乳无遮挡免费网站照片| 精品久久久久久久久久久久久| 三级男女做爰猛烈吃奶摸视频| 国产男靠女视频免费网站| 女生性感内裤真人,穿戴方法视频| 女人十人毛片免费观看3o分钟| 在线免费观看不下载黄p国产 | 一个人免费在线观看电影| 99热6这里只有精品| 麻豆一二三区av精品| 麻豆久久精品国产亚洲av| 亚洲av免费在线观看| 亚洲人成伊人成综合网2020| 人妻制服诱惑在线中文字幕| 久久久成人免费电影| 欧美乱妇无乱码| 成人一区二区视频在线观看| 精品人妻1区二区| 最近中文字幕高清免费大全6 | 亚洲av日韩精品久久久久久密| 午夜激情福利司机影院| 搞女人的毛片| 91av网一区二区| 69av精品久久久久久| 午夜久久久久精精品| 精品一区二区三区视频在线观看免费| 长腿黑丝高跟| 露出奶头的视频| x7x7x7水蜜桃| 午夜视频国产福利| 日本五十路高清| 69人妻影院| 欧美高清性xxxxhd video| 国产老妇女一区| 琪琪午夜伦伦电影理论片6080| 99久久精品国产亚洲精品| 两性午夜刺激爽爽歪歪视频在线观看| 精品久久久久久久久亚洲 | 久久人人爽人人爽人人片va | 男人狂女人下面高潮的视频| 在线观看66精品国产| 亚洲 国产 在线| 极品教师在线免费播放| 欧美国产日韩亚洲一区| 一区二区三区高清视频在线| 黄色视频,在线免费观看| 亚洲美女搞黄在线观看 | 日韩高清综合在线| 亚洲欧美日韩高清在线视频| 国产精品久久久久久亚洲av鲁大| 69人妻影院| 中文在线观看免费www的网站| 1000部很黄的大片| 国产一区二区三区在线臀色熟女| 亚洲国产精品久久男人天堂| 12—13女人毛片做爰片一| 噜噜噜噜噜久久久久久91| 成人亚洲精品av一区二区| 全区人妻精品视频| 哪里可以看免费的av片| 国产日本99.免费观看| 成人国产综合亚洲| 国产成年人精品一区二区| 女人被狂操c到高潮| 亚洲成人免费电影在线观看| 精品久久国产蜜桃| 亚洲av不卡在线观看| 亚洲中文字幕一区二区三区有码在线看| 欧美日韩中文字幕国产精品一区二区三区| 免费看日本二区| 日本黄色视频三级网站网址| av在线天堂中文字幕| 校园春色视频在线观看| 身体一侧抽搐| 在线观看美女被高潮喷水网站 | 欧美日韩瑟瑟在线播放| 搞女人的毛片| 国产精品,欧美在线| 成人鲁丝片一二三区免费| 亚洲欧美日韩东京热| 国产精品一区二区三区四区久久| 午夜激情欧美在线| 婷婷亚洲欧美| 精品乱码久久久久久99久播| 很黄的视频免费| .国产精品久久| 国产亚洲欧美在线一区二区| 婷婷丁香在线五月| 又粗又爽又猛毛片免费看| 美女cb高潮喷水在线观看| 日日摸夜夜添夜夜添av毛片 | 国产三级在线视频| 18禁裸乳无遮挡免费网站照片| 国产伦精品一区二区三区四那| 久久国产乱子免费精品| 性色avwww在线观看| 搞女人的毛片| 欧美黑人欧美精品刺激| 老司机午夜十八禁免费视频| 黄色一级大片看看| 男人舔女人下体高潮全视频| 欧美乱妇无乱码| 日韩国内少妇激情av| 床上黄色一级片| 又爽又黄a免费视频| 午夜福利18| 欧美乱妇无乱码| 欧美精品啪啪一区二区三区| 亚洲精品色激情综合| 成熟少妇高潮喷水视频| 久久久久久九九精品二区国产| 欧美成人a在线观看| 人人妻人人看人人澡| 97超视频在线观看视频| 又爽又黄无遮挡网站| 国内少妇人妻偷人精品xxx网站| 久久亚洲真实| 在线观看舔阴道视频| 97超视频在线观看视频| 欧美成人免费av一区二区三区| 91字幕亚洲| 国产亚洲精品久久久com| 免费一级毛片在线播放高清视频| 一区福利在线观看| 99久久精品国产亚洲精品| 欧美三级亚洲精品| 日日干狠狠操夜夜爽| 淫妇啪啪啪对白视频| 搡老岳熟女国产| 午夜精品久久久久久毛片777| 日本成人三级电影网站| 日本a在线网址| 亚洲中文日韩欧美视频| 在线观看一区二区三区| 国产69精品久久久久777片| 12—13女人毛片做爰片一| 天天躁日日操中文字幕| 亚洲中文日韩欧美视频| 热99在线观看视频| 性欧美人与动物交配| 日本免费一区二区三区高清不卡| 欧美xxxx黑人xx丫x性爽| 在线观看舔阴道视频| 国产成年人精品一区二区| 亚洲精品乱码久久久v下载方式| av中文乱码字幕在线| 99在线视频只有这里精品首页| 婷婷精品国产亚洲av| 男女之事视频高清在线观看| 精品久久久久久久人妻蜜臀av| 久久中文看片网| 国产久久久一区二区三区| 成人三级黄色视频| 国内精品久久久久久久电影| 三级国产精品欧美在线观看| 天天躁日日操中文字幕| 久久精品国产亚洲av涩爱 | 99国产精品一区二区蜜桃av| 国产探花极品一区二区| 琪琪午夜伦伦电影理论片6080| 天堂影院成人在线观看| 日韩中字成人| 美女高潮的动态| 成年女人看的毛片在线观看| 哪里可以看免费的av片| 欧美另类亚洲清纯唯美| 国产高清视频在线播放一区| 国产精品,欧美在线| 亚洲人成伊人成综合网2020| 三级国产精品欧美在线观看| 亚洲五月天丁香| 伊人久久精品亚洲午夜| 乱人视频在线观看| 看片在线看免费视频| 国产高清视频在线观看网站| 国内久久婷婷六月综合欲色啪| av中文乱码字幕在线| 亚洲人成电影免费在线| 国产午夜福利久久久久久| 午夜福利在线观看吧| 黄色配什么色好看| 直男gayav资源| 少妇人妻一区二区三区视频| 最后的刺客免费高清国语| 国产午夜精品论理片| 琪琪午夜伦伦电影理论片6080| 99热这里只有是精品50| 香蕉av资源在线| 精品久久久久久久久久久久久| 最新中文字幕久久久久| 男人狂女人下面高潮的视频| 亚洲国产精品999在线| 国产日本99.免费观看| 老熟妇仑乱视频hdxx| 久久久色成人| 亚洲成人免费电影在线观看| 午夜老司机福利剧场| 波野结衣二区三区在线| 两个人视频免费观看高清| 两性午夜刺激爽爽歪歪视频在线观看| 九九热线精品视视频播放| 天堂av国产一区二区熟女人妻| 国产精品日韩av在线免费观看| 看免费av毛片| 少妇的逼好多水| 麻豆一二三区av精品| 长腿黑丝高跟| 高潮久久久久久久久久久不卡| 看免费av毛片| 国产综合懂色| 日韩精品青青久久久久久| 成人毛片a级毛片在线播放| 我要搜黄色片| 国产色爽女视频免费观看| 99久久久亚洲精品蜜臀av| 亚洲狠狠婷婷综合久久图片| 啦啦啦观看免费观看视频高清| 国产精品久久久久久久电影| 日本免费a在线| 一级av片app| 1000部很黄的大片| 69人妻影院| 中文字幕人妻熟人妻熟丝袜美| 精品久久久久久,| 亚洲国产色片| 婷婷精品国产亚洲av| 女人十人毛片免费观看3o分钟| 91在线观看av| 亚洲自偷自拍三级| 色av中文字幕| 18禁黄网站禁片免费观看直播| 最后的刺客免费高清国语| av专区在线播放| 亚洲自偷自拍三级| 啦啦啦观看免费观看视频高清| 男女床上黄色一级片免费看| 久久伊人香网站| 国产伦精品一区二区三区视频9| 亚洲人成网站在线播| 2021天堂中文幕一二区在线观| 国产成人aa在线观看| 欧美国产日韩亚洲一区| 3wmmmm亚洲av在线观看| 日本黄色片子视频| 精品人妻视频免费看| 久久性视频一级片| av天堂中文字幕网| 欧美xxxx黑人xx丫x性爽| 欧美日韩瑟瑟在线播放| 极品教师在线免费播放| 亚洲av中文字字幕乱码综合| 亚洲熟妇熟女久久| 欧美潮喷喷水| 少妇的逼好多水| 午夜精品一区二区三区免费看| 88av欧美|