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

    基于納米有機(jī)液體的高性能溫差電池

    2016-11-22 09:49:08孫浩宇濮金歡唐桂華
    物理化學(xué)學(xué)報 2016年10期
    關(guān)鍵詞:熱電電勢氯化鈉

    孫浩宇 濮金歡 唐桂華

    (西安交通大學(xué)能源與動力工程學(xué)院,熱流科學(xué)與工程教育部重點實驗室,西安710049)

    基于納米有機(jī)液體的高性能溫差電池

    孫浩宇濮金歡唐桂華*

    (西安交通大學(xué)能源與動力工程學(xué)院,熱流科學(xué)與工程教育部重點實驗室,西安710049)

    傳統(tǒng)液體溫差電池較低的熱電轉(zhuǎn)換性能一直無法得到有效改善,亟需尋找新的熱電轉(zhuǎn)換機(jī)制來提升熱電轉(zhuǎn)換效率。本文采用分子動力學(xué)(MD)方法,數(shù)值模擬研究了不同溫度下以不同配比的甘油-水為溶劑的氯化鈉溶液在碳納米管(CNT)內(nèi)離子、分子分布情況。結(jié)果表明:離子、分子的分布受溫度影響較大,近壁面凈電荷、電勢分布隨溫度升高出現(xiàn)明顯的分層。根據(jù)模擬結(jié)果提出以CNT為電極材料,甘油氯化鈉溶液或甘油水氯化鈉溶液為電解質(zhì)溶液組成溫差電池。其熱電轉(zhuǎn)換性能遠(yuǎn)優(yōu)于大多數(shù)溫差電池,同時溫度適用范圍也顯著增加。以多孔碳為電極材料,甘油氯化鈉溶液為電解質(zhì)溶液組成的熱電轉(zhuǎn)換裝置實驗驗證了可行性。

    分子動力學(xué);溫差電池;碳納米管;低溫能源利用;熱電轉(zhuǎn)換

    本文提出采用甘油氯化鈉溶液和甘油水二元混合氯化鈉溶液作為電解質(zhì)溶液,可明顯降低溶液的表面張力,使溶液更易于浸潤C(jī)NT,穩(wěn)定性也相應(yīng)地提高,同時隨著甘油在溶劑中的比例增大,相同溫差下產(chǎn)生的電勢也隨之增大,熱電轉(zhuǎn)換性能顯著增加。除此之外甘油與水相互混合增加了熱電轉(zhuǎn)換的溫度范圍,不再局限于0-100°C內(nèi),顯著增加了液體熱電轉(zhuǎn)換的適用范圍。

    2 數(shù)值方法及物理模型

    采用分子動力學(xué)方法,基于開源軟件Lammps26進(jìn)行計算。物理模型如圖1所示,計算中選用CNT(20,20)扶手椅型碳納米管并固定,其半徑R=1.356 nm,計算長度l=9.838 nm,不同模擬工況其電解質(zhì)溶液組分見表1,初始溫度T=300 K,壓力為p=0.1 MPa。本文采用Class 2力場中的PCFF力場(force field for polycarbonate)以及宏觀正則系綜(NVT系綜,其中N代表體系中原子個數(shù),V代表體積,T代表溫度),CNT與溶液中分子、離子之間的勢能函數(shù)為E=εij[2(σij/rij)9-3(σij/ rij)6],交互作用參數(shù)采用Lorentz-Berthelot混合計算法則,截斷半徑設(shè)為1 nm,盒子軸向采用周期性邊界條件,時間步長為1 fs。模擬1 ns使溶液達(dá)到平衡階段,平衡之后繼續(xù)模擬1 ns并進(jìn)行數(shù)據(jù)輸出,總模擬時間為2 ns。

    圖1  模擬的物理模型Fig.1 Physical model of simulation

    表1 甘油/甘油-水電解質(zhì)溶液各組分含量Table 1 Content of each component in glycerine/glycerolaqueous electrolyte solution

    3 結(jié)果與分析

    3.1CNT內(nèi)分子和離子分布特性

    圖2為溫度300 K時CNT內(nèi)電解質(zhì)溶液中離子、分子的分布情況,橫坐標(biāo)表示距離管徑中心的距離。圖2(A)為不同工況下Na+分布情況,圖2 (B)為H(甘油分子中與碳原子直接成健的氫原子)分布情況。從圖2(A)中可知,不同工況下近壁面處Na+具有相似分布特性,從圖2(B)可知,不同工況下近壁面處H的分布特性同樣相似。事實上不同工況下溶液中其他離子、分子均表現(xiàn)出相似的分布特性,本文不再贅述。由于不同工況下CNT內(nèi)離子、分子分布特性相似,下文選取溶劑中甘油含量為60.10%的氯化鈉電解質(zhì)溶液(簡稱60.10%甘油電解質(zhì)溶液)進(jìn)行分析。圖2(C)為Na+和Cl-在CNT中徑向分布情況;圖2(D)為甘油分子和水分子在CNT內(nèi)分布情況,左邊縱坐標(biāo)是甘油分子中各原子(C1表示分子鏈中間碳原子,C2表示碳鏈兩端的碳原子,H表示直接與碳原子連接的氫原子,HO表示與氧原子相連的氫原子,OH表示與氫原子相連的氧原子)密度在CNT中徑向分布數(shù)值,右邊縱坐標(biāo)是水分子中各原子(H1表示水中的氫原子,O1表示水中的氧原子)密度在CNT中徑向分布數(shù)值。由圖3可知,陰陽離子以及分子密度在近壁面處出現(xiàn)大的波動。溶液中的H原子距離壁面比O和C原子近,Na+距離壁面比Cl-近。前者是因為H受到CNT壁面的斥力比溶液中的O、C弱;后者主要因為Na+的半徑比Cl-的半徑小,所以Na+受到來自H原子的引力較強(qiáng),受到O、C原子的斥力較弱22。

    圖2 300 K時CNT中電解質(zhì)溶液離子、分子分布Fig.2 Number density distribution of the ions and molecules in electrolyte solution of the CNT at 300 K

    根據(jù)Poisson方程可以得到相應(yīng)的電勢φ分布,Poisson方程如下

    式中,ε0表示真空介電常數(shù);ε表示溶液的介電常數(shù);ρne(r)表示凈電荷密度,可由下式獲得:

    式中,n表示溶液中粒子總數(shù);ci表示第i個離子或原子的電荷量。

    根據(jù)二元無規(guī)混合物有效介電常數(shù)的計算公式,可得到甘油水氯化鈉溶液的有效介電常數(shù)ε,其計算公式如下27

    式中,ω表示溶質(zhì)質(zhì)量分?jǐn)?shù);ε1b=2ω/(3-ω)ε21,εb=2(1-ω)/(2+ω)ε1,ε2表示氯化鈉的介電常數(shù),ε2表示甘油或甘油水混合溶液的介電常數(shù),如表2所示28。

    圖3(A)為T=300 K時60.10%甘油電解質(zhì)溶液在CNT內(nèi)凈電荷密度分布情況??梢姡贑NT管壁附近的凈電荷密度波動較大,這是由于CNT管壁附近離子、原子密度分布不均勻,并且相對集中;同時在CNT內(nèi)近壁面凈電荷密度分布隨著溫度的升高而下降,這是因為溫度越高,近壁面附近的離子越易獲得能量,脫離近壁面的束縛,導(dǎo)致凈電荷密度降低,即溫度越高,近壁面內(nèi)凈電荷密度越低。圖3(B)表示不同溫度時60.10%甘油電解質(zhì)溶液在CNT內(nèi)電勢分布??梢?,近壁面電勢隨溫度升高而增大,這是由于CNT內(nèi)近壁面凈電荷密度分布隨溫度升高而下降所導(dǎo)致。圖3(C)為T=300 K時不同工況CNT內(nèi)凈電荷密度分布。觀察到以純水為溶劑的氯化鈉溶液與以甘油為溶劑的氯化鈉溶液相比波動的峰值比較平均,純甘油氯化鈉溶液的凈電荷主要集中在近壁面附近;不同比例的甘油水混合飽和氯化鈉溶液凈電荷密度分布具有相同的波動規(guī)律,同時隨甘油含量的增加和氯化鈉含量的減少而降低。圖3(D)給出了各工況處于不同溫度CNT相互連接(如圖3(E))輸出電勢φ隨溫差ΔT變化,同時與Chen等22模擬結(jié)論進(jìn)行對比。可知,輸出電勢φ隨溫差ΔT增大而增加;純甘油氯化鈉溶液的輸出電勢φ明顯大于純水氯化鈉溶液的輸出電勢;甘油水二元混合飽和氯化鈉溶液的輸出電勢φ隨甘油含量的增加及氯化鈉含量的減少而升高,而根據(jù)Chen等22研究,純水氯化鈉溶液的輸出電勢φ應(yīng)該隨純水氯化鈉濃度減少而降低,所以輸出電勢φ逐漸增大主要是由甘油含量增加所導(dǎo)致。由圖3(D)可知,本文模擬的氯化鈉水溶液的結(jié)果與Chen等模擬結(jié)論相符合。

    圖3 CNT內(nèi)電解質(zhì)溶液凈電荷密度及電勢分布情況Fig.3 Distribution of net charge and electric potential in the CNT

    表2 不同比例的甘油和水混合物在不同溫度下介電常數(shù)(ε2)Table 2 Dielectric constants(ε2)of glycerol-water solutions under different temperatures

    3.2熱電轉(zhuǎn)換效率與Seebeck系數(shù)

    熱電轉(zhuǎn)換效率公式為η=Winput/Woutput,其中Winput表示輸入功,Woutput表示輸出功。液體溫差電池?zé)犭娹D(zhuǎn)換過程類似于電容放電,輸出功的過程也相當(dāng)于電極放電做功,所以輸出功Woutput可由下式給出22

    其中U=Q/CT,則:

    式中CT表示電容,ΔQ表示最大電荷變量,Q為近壁面總的初始電荷量。由于電容是由陰、陽極電容組成,故:

    式中,Ca表示陽極電容;Cc表示陰極電容29。本文所選取的CNT屬于介孔(2 nm

    系統(tǒng)的輸入功Winput為:

    式中,Ci表示各原子/離子的比熱容。對于甘油分子、水分子、C、Na+和Cl-的比熱分別為202.58、74.54、8.52、27.60和33.95 J·mol-1。為了更好地描述系統(tǒng)的性能,通常采用其相對于卡諾循環(huán)的效率ηr: Seebeck系數(shù)計算公式為α=?V/?T,V表示電勢,T表示溫度;α表示電勢隨溫度變化快慢。

    圖5(A)為不同工況下熱電轉(zhuǎn)換Seebeck系數(shù)α與溫差ΔT的關(guān)系,發(fā)現(xiàn)純甘油氯化鈉溶液和甘油水氯化鈉溶液熱電轉(zhuǎn)換的Seebeck系數(shù)隨溫差ΔT增加而增加,而純水氯化鈉溶液的Seebeck系數(shù)α隨溫差ΔT升高變化不明顯;當(dāng)溫差ΔT=40,60 K時,甘油水氯化鈉溶液的α隨甘油含量的增加而增大。這說明在較大溫差時,隨著甘油含量的增加CNT內(nèi)電解質(zhì)溶液對溫度的敏感性越高。圖5(B)表示不同工況下熱電轉(zhuǎn)換效率η隨溫差變化,發(fā)現(xiàn)相同氯化鈉濃度時,純甘油氯化鈉電解質(zhì)溶液的效率η要遠(yuǎn)高于純水氯化鈉溶液;甘油水二元混合飽和氯化鈉溶液時,熱電轉(zhuǎn)換效率隨甘油含量的增加(氯化鈉濃度減小)而增加,根據(jù)Chen等22研究可知,氯化鈉濃度減小,熱電轉(zhuǎn)換效率減少,所以甘油含量的增加是引起熱電轉(zhuǎn)換效率升高的主要原因。圖5(C)為不同工況下該機(jī)制相對于卡諾循環(huán)的熱電轉(zhuǎn)換效率ηr隨溫差的變化。純甘油氯化鈉溶液和甘油水氯化鈉溶液熱電轉(zhuǎn)換ηr均隨溫差增加而減小,說明該機(jī)制的熱電轉(zhuǎn)換效率溫差越小越接近卡諾循環(huán)效率,溫差越大越偏離卡諾循環(huán)。

    綜上所述,在CNT內(nèi)氯化鈉濃度相同時氯化鈉甘油溶液比氯化鈉水溶液性能更佳;甘油水二元混合飽和氯化鈉電解質(zhì)溶液的熱電轉(zhuǎn)換性能在溫差較小、甘油含量較多情況下更加優(yōu)異。本文熱電轉(zhuǎn)換機(jī)制在ΔT=60 K時,相對于卡諾循環(huán)的熱電轉(zhuǎn)換效率ηr高達(dá)~11.30%,而大多數(shù)液體溫差電池相對于卡諾循環(huán)的熱電轉(zhuǎn)換效率ηr僅為~1.40%13,16,20,30。

    圖4 電解質(zhì)溶液中離子在介孔內(nèi)分布示意圖Fig.4 Schematic of distribution of ions in the mesopore in electrolyte solution

    圖5 不同濃度甘油溶劑電解質(zhì)溶液的α、η、ηr隨溫差ΔT的變化情況Fig.5 Changes of α,η,and ηrwith ΔT in different electrolyte solutions

    3.3實驗驗證

    當(dāng)電解質(zhì)溶液限制在CNT中,由于陰陽離子徑向分布的差異,以及溶劑分子中各原子分布不同而產(chǎn)生偶極子的有序取向31,32,會導(dǎo)致其內(nèi)部電荷分布不均勻,產(chǎn)生凈電荷及電勢。通過模擬分析可知,CNT內(nèi)凈電荷及電勢分布受溫度影響較大,由此將處于不同溫度下的CNT連接形成回路,能夠輸出電勢φ,實現(xiàn)熱能向電能的轉(zhuǎn)換。本實驗采用碳納米多孔材料炭黑代替CNT驗證該機(jī)制可行性。

    實驗以工業(yè)炭黑(純度>99.5%,天津金秋實化工有限公司)為電極材料,甘油水-氯化鈉溶液(甘油純度>99%,氯化鈉純度>99.5%)為電解質(zhì)溶液組成熱電轉(zhuǎn)換裝置。首先,將炭黑用丙酮(純度>98.5%)清洗4 h晾干;然后將其壓制成直徑為10 mm的圓盤電極,每個圓盤電極的質(zhì)量電極大約為90 mg;最后用導(dǎo)電銀膠將電極與邊長為5 mm的銅箔(作為集流器)連接并放入DZ-3AII真空干燥箱(奧環(huán)電子產(chǎn)品有限公司)內(nèi),120°C烘干4 h,并將連接集流器側(cè)的電極表面進(jìn)行防水處理。實驗系統(tǒng)示意圖如圖6所示。實驗采用F26-ME恒溫水浴系統(tǒng)(德國Julabo公司)保證兩端溫度恒定,冷、熱兩端采用直徑為5 mm的鹽橋連接。冷端溫度保持25°C不變,通過調(diào)節(jié)高溫端恒溫水浴的溫度來改變溫差ΔT。用Keithley 2700數(shù)據(jù)采集器(美國Keithley公司)采集溫度T以及不同溫差ΔT下輸出的電壓φ。

    圖6 納米有機(jī)液體溫差電池實驗示意圖Fig.6 Schematic of the thermogalvanic cell based on organic nanofluids

    圖7 電極兩端電勢隨溫差ΔT變化情況Fig.7 Output electric potential between two electrodes

    通過實驗得到甘油氯化鈉電解質(zhì)溶液溫差電池電勢隨溫差變化情況,如圖7所示。可知,電勢隨著溫差ΔT的增大而增大,相同氯化鈉濃度下氯化鈉甘油電解質(zhì)溶液產(chǎn)生的電勢高于氯化鈉水電解質(zhì)溶液,甘油/水混合氯化鈉電解質(zhì)溶液電勢隨甘油含量的增加而增大,與數(shù)值模擬結(jié)果一致。同時發(fā)現(xiàn),當(dāng)氯化鈉濃度相同時,溶劑中甘油含量下降到40%,電勢隨溫差ΔT升高而增加的量明顯下降,與純水氯化鈉電解質(zhì)溶液接近,說明此時電解質(zhì)溶液中甘油對離子分布的影響下降,本文認(rèn)為是因為氯化鈉在水中的溶解度大于在甘油中的溶解度導(dǎo)致。實驗結(jié)果與數(shù)值模擬結(jié)果定量上的差異是由于實驗所用的多孔電極的孔徑及導(dǎo)電、導(dǎo)熱性能和表面粗糙度均與碳納米管存在差異。

    4 結(jié)論

    通過MD方法模擬CNT-甘油氯化鈉溶液及CNT-甘油水氯化鈉溶液的納米有機(jī)液體溫差電池的熱電轉(zhuǎn)換過程以及實驗驗證研究,獲得以下主要結(jié)論:

    (1)CNT內(nèi)溶液中的離子、分子分布不均,近壁面凈電荷密度ρne隨溫度升高而下降,近壁面電勢φ隨溫差增加而增加。

    (2)CNT內(nèi)近壁面凈電荷密度ρne隨甘油含量增加以及氯化鈉含量的降低而降低。

    (3)將處于不同溫度下的CNT相互連接形成回路能夠?qū)崿F(xiàn)熱能向電能的轉(zhuǎn)換,熱電轉(zhuǎn)換Seebeck系數(shù)及熱電轉(zhuǎn)換效率隨溫差的增加而增大,甘油/水二元混合飽和氯化鈉溶液時,熱電轉(zhuǎn)換效率隨甘油含量增加而增加;低溫差時具有較高的ηr,其性能優(yōu)于傳統(tǒng)溫差電池。

    (4)實驗驗證發(fā)現(xiàn),氯化鈉濃度相同時,氯化鈉甘油溶液熱電轉(zhuǎn)換性能要優(yōu)于氯化鈉水溶液,甘油水氯化鈉電解質(zhì)溶液電勢隨甘油含量增加而增加。

    另外,CNT孔徑、氯化鈉的濃度、初始溫度等均會影響熱電轉(zhuǎn)換性能。該熱電轉(zhuǎn)換機(jī)制具有以下優(yōu)點:電極與溶液不存在化學(xué)反應(yīng),所以對電極材料的損耗極低,碳的排放量為零;電極材料可由工業(yè)炭黑代替,成本低;不同比例甘油水混合溶劑的凝固點和沸點不再局限于0和100°C,熱電轉(zhuǎn)換溫度范圍顯著增加。該熱電轉(zhuǎn)換機(jī)制具有較高的轉(zhuǎn)換效率,對將來的低溫?zé)犭娹D(zhuǎn)換研究具有指導(dǎo)意義。

    References

    (1) Heremans,J.P.;Dresselhaus,M.S.;Bell,L.E.;Morelli,D.T. Nat.Nanotech.2013,8,471.doi:10.1038/nnano.2013.129

    (3) Snyder,J.G.;Toberer,E.S.Nat.Mater.2008,7,105. doi:10.1038/nmat2090

    (5) Poudel,B.;Hao,Q.;Ma,Y.;Lan,Y.C.;Minnich,A.;Yu,B.; Yan,X.;Wang,D.Z.;Muto,A.;Vashaee,D.;Chen,X.Y.;Liu, J.M.;Dresselhaus,M.S.;Chen,G.;Ren,Z.F.Science 2008, 320,634.doi:10.1126/science.1156446

    (6) Kongtragool,B.;Wongwises,S.Renew.Sust.Energy Rev.2003, 7,131.doi:10.1016/S1364-0321(02)00053-9

    (7)Kong,X.O.;Wang,R.Z.;Huang,X.H.Energ.Convers. Manage.2004,45,433.doi:10.1016/j.enconman.2003.09.009

    (8) Kongtragool,B.;Wongwises,S.Renew.Energy 2007,32,547. doi:10.1016/j.renene.2006.03.003

    (9) Thombare,D.G.;Verma,S.K.Renew.Sust.Energy Rev.2008, 12,1.doi:10.1016/j.rser.2006.07.001

    (11) Mancini,T.;Heller,P.;Butler,B.;Osborn,B.;Schiel,W.; Goldberg,V.;Buck,R.;Diver,R.;Andraka,C.;Moreno,J. ASME J.Sol.Energy 2003,125,135.doi:10.1115/1.1562634

    (12) Kuzminskii,Y.V.;Zasukha,V.A.;Kuzminskaya,G.Y. J.Power Sources 1994,52,231.doi:10.1016/0378-7753(94) 02015-9

    (13) Quickenden,T.I.;Mua,Y.J.Electrochem.Soc.1995,142, 3985.doi:10.1149/1.2048446

    (15) Debethune,A,J.;Licht.T.S.;Swendeman,N.J.Electrochem. Soc.1959,106,616.doi:10.1149/1.2427448

    (16) Gunawa,A.;Lin,C.H.;Buttry,D.A.;Mujica,V.;Taylor,R.A.; Prasher,R.S.;Phelan,P.E.Nano.Micro.Therm.2013,17(4), 304.doi:10.1080/15567265.2013.776149

    (17) Quickenden,T.I.;Vernon,C.F.Sol.Energy 1986,36,63. doi:10.1016/0038-092X(86)90061-7

    (18) Kang,T.J.;Fang,S.;Kozlov,M.E.;Haines,C.S.;Li,N.;Kim, Y.H.;Chen,Y.S.;Baughman,R.H.Adv.Funct.Mater.2012, 22,477.doi:10.1002/adfm.201101639

    (19)Yu,Q.;Punyamurtual,V.K.;Han,A.;Lim,H.J.Power Sources 2008,183,403.doi:10.1016/j.jpowsour.2008.05.008

    (20) Hu,R.C.;Cola,B.A.;Haram,N.;Barisci,J.N.;Lee,S.; Stoughton,S.;Wallace,G.;Too,C.;Thomas,M.;Gestos,A.; Cruz,M.E.D.;Ferraris,J.P.;Zakhidow,A.A.;Baughman,R. H.Nano Lett.2010,10,838.doi:10.1021/nl903267n

    (21) Chen,X.;Xu,B.X.;Liu,L.Appl.Mech.Rev.2014,66,050803. doi:10.1115/1.4026913

    (22) Xu,B.X.;Liu,L.;Lim,H.;Qiao,Y.;Chen,X.Nano Energy2012,1,805.doi:10.1016/j.nanoen.2012.07.013

    (23) Xu,B.X.;Yu,Q.;Park,T.;Tak,M.;Zhou,Q.L.;Chen,X. Energy Environ.Sci.2011,4,3632.doi:10.1039/c1ee01405f

    (24)Li,Q.;Yang,D.F.;Wang,J.H.;Wu,Q.;Liu,Q.Z.Acta Phys.-Chim.Sin.2016,32(3),691.[李清,楊登峰,王建花,武琪,劉清芝.物理化學(xué)學(xué)報,2016,32(3),691.] doi:10.3866/PKU.W HXB201512182

    (25) Shen,Z.L.;He,G.H.;Zhang,N.;Hao,C.Acta Phys.-Chim. Sin.2015,31(6),1025.[諶莊琳,賀高紅,張寧,郝策.物理化學(xué)學(xué)報,2015,31(6),1025.]doi:10.3866/PKU. WHXB201504141

    (27) Xue,Q.Z.;Li,W.Y.J.University Petroleum 1999,23(4),102. [薛慶忠,李文瀛.石油大學(xué)學(xué)報,1999,23(4),102.]doi:1000-5870(1999)04-0102-03

    (28) Association,G.P.Physical Properties of Glycerine and Its Solutions;Glycerine Producers′Association:New York 1963; p 17.

    (29) Zhang,L.L.;Zhao,X.S.Chem.Soc.Rev.2009,38,2520. doi:10.1039/b813846j

    (30) Mua,Y.;Quickenden,T.I.J.Electrochem.Soc.1996,143(8), 2558.doi:10.1149/1.1837047

    (31) Mann,D.J.;Halls,M.D.Phys.Rev.Lett.2003,90,195503. doi:10.1103/PhysRevLett.90.195503

    High-Performance Thermogalvanic Cell Based on Organic Nanofluids

    SUN Hao-YuPU Jin-HuanTANG Gui-Hua*
    (Key Laboratory of Thermal-fluid Science and Engineering of Ministry of Education,School of Energy and Power Engineering,Xi′an Jiaotong University,Xi′an 710049,P.R.China)

    The thermoelectric performance of traditional thermogalvanic cells is relatively low and a more efficient conversion mechanism is required.In this paper,the distribution of glycerol/glycerin in an aqueous sodium chloride solution in a carbon nanotube(CNT)is investigated by molecular dynamics(MD)simulation. The distributions of ions,molecule net charge,and electrical potential of the system are markedly affected by temperature.We propose a novel nanofluid thermoelectric conversion method based on the CNT and glycerol/ glycerin aqueous sodium chloride solution.The thermoelectric performance of the proposed system is much higher than that of most of current liquid thermogalvanic cells,and the application temperature range is also widened considerably.Apreliminary thermal-to-electrical energy conversion experiment based on nanoporous carbon with mixtures of sodiumchloride and glycerol is also conducted to qualitatively verify the numerical results.

    Molecular dynamics;Thermogalvanic cell;Carbon nanotube;Low-temperature energy utilization;Thermoelectric conversion

    1 引言

    隨著不可再生能源的日益枯竭以及能源需求增大,低溫能源(工業(yè)廢水、汽車尾氣、地?zé)嵋约疤柲艿?的開發(fā)與利用受到越來越多的關(guān)注,其對改善能源環(huán)境問題及可持續(xù)發(fā)展具有深遠(yuǎn)的意義。針對應(yīng)用固體熱電材料1-5和斯特林發(fā)動機(jī)的低溫能源利用進(jìn)行了大量的研究6-9,雖然取得了較大的進(jìn)步,但是仍存在成本高、經(jīng)濟(jì)性低以及可靠性差10,11等問題。因此有學(xué)者將注意力轉(zhuǎn)移到液體溫差電池上12,13,利用電解質(zhì)溶液氧化還原對以實現(xiàn)熱電轉(zhuǎn)換,例如Cu-CuSO4和Zn-ZnSO4、FeCN3--FeCN4-以及Fe3+-Fe2+等14-18。然而依靠氧化還原對實現(xiàn)熱電轉(zhuǎn)換的經(jīng)濟(jì)性和熱電轉(zhuǎn)換性能仍不能滿足工業(yè)需求。為了突破液體溫差電池的技術(shù)瓶頸,部分學(xué)者開始研究碳納米管(CNT)中電解質(zhì)溶液的性質(zhì),并將其成功應(yīng)用于低溫?zé)犭娹D(zhuǎn)換過程。Yu等19以多孔碳為電極,氯化鈉水溶液為電解質(zhì)溶液組成的裝置有效地將熱能轉(zhuǎn)換成電能。Hu等20以CNT為電極,F(xiàn)eCN3--FeCN4-水溶液為電解質(zhì)溶液制作成液體溫差電池,得到較高的電勢和熱電轉(zhuǎn)換效率。Chen等21-23研究發(fā)現(xiàn)當(dāng)氯化鈉水溶液被限制在CNT內(nèi)時,溶液中的分子、離子分布會發(fā)生各向異性變化,導(dǎo)致CNT管壁產(chǎn)生電勢,并且該電勢大小受溫度影響較大,利用這個機(jī)理能夠?qū)崿F(xiàn)熱電轉(zhuǎn)換且具有較高的轉(zhuǎn)換效率;同時CNT在熱能、電能、機(jī)械能相互轉(zhuǎn)換及海水淡化24,25等方面也具有優(yōu)良的性能。但是水溶液存在表面張力較大、不易浸入CNT以及穩(wěn)定性差、易受環(huán)境影響等缺點。

    May 13,2016;Revised:June 27,2016;Published online:June 28,2016.

    .Email:ghtang@mail.xjtu.edu.cn;Tel:+86-29-82665319.

    O646

    10.3866/PKU.WHXB201606281

    The project was supported by the National Natural Science Foundation of China(51576156).國家自然科學(xué)基金(51576156)資助項目?Editorial office ofActa Physico-Chimica Sinica

    (2) Tritt,T.M.Annu.Rev.Mater.Res.2011,41,433.10.1146/ annurev-matsci-062910-100453

    (4) Bell,L.E.Science 2008,321,1457.10.1126/ science.1158899

    (10) Vining,C.B.Nat.Mater.2009,8,83.10.1038/nmat2361

    (14) Bouty,E.J.Phys.1880,9,229.10.1051/jphystap: 018800090022901

    (26) Plimpton,S.J.Comput.Phys.1995,117,1.10.1006/ jcph.1995.1039

    (32) Zhu,F.;Schulten,K.Biophys.J.2003,85,236.10.1016/ S0006-3495(03)74469-5

    猜你喜歡
    熱電電勢氯化鈉
    福州熱電兩臺660MW熱電聯(lián)產(chǎn)工程核準(zhǔn)獲批
    氯化鈉溶液的配制
    同步機(jī)和異步機(jī)之磁勢和電勢的比較
    氯化鈉水溶液結(jié)構(gòu)的研究
    熱電轉(zhuǎn)換材料的開發(fā)與應(yīng)用
    電力與能源(2017年6期)2017-05-14 06:19:41
    熱油中的食鹽為何不溶化?
    奧秘(2016年12期)2016-12-17 15:48:57
    新型熱電制冷裝置的實驗開發(fā)
    氯化鈉
    熱泵在熱電聯(lián)產(chǎn)中的應(yīng)用
    河南科技(2015年15期)2015-03-11 16:25:52
    靜電場中計算電勢的幾種方法
    精品人妻1区二区| 99国产精品一区二区三区| 热99在线观看视频| 中文字幕av在线有码专区| 国产久久久一区二区三区| 免费大片18禁| 小说图片视频综合网站| 亚洲av成人精品一区久久| 精品99又大又爽又粗少妇毛片 | 美女高潮喷水抽搐中文字幕| 成人三级黄色视频| av在线蜜桃| 国内精品久久久久久久电影| 免费搜索国产男女视频| 亚洲最大成人av| 精品国产亚洲在线| 亚洲最大成人中文| 美女被艹到高潮喷水动态| 91字幕亚洲| 亚洲av免费高清在线观看| www.999成人在线观看| 此物有八面人人有两片| 午夜两性在线视频| 又爽又黄a免费视频| 90打野战视频偷拍视频| 丁香欧美五月| 成人特级黄色片久久久久久久| 国产成+人综合+亚洲专区| 1024手机看黄色片| 日本精品一区二区三区蜜桃| 亚洲精品亚洲一区二区| 91九色精品人成在线观看| 亚洲一区高清亚洲精品| www.www免费av| 琪琪午夜伦伦电影理论片6080| 午夜福利18| 国产视频一区二区在线看| 成人欧美大片| 亚洲人成网站在线播| 成人av在线播放网站| 久久亚洲真实| 国产一区二区亚洲精品在线观看| 亚洲成人精品中文字幕电影| 亚洲 欧美 日韩 在线 免费| 国产精品乱码一区二三区的特点| 人妻久久中文字幕网| 狂野欧美白嫩少妇大欣赏| 久久久久久久午夜电影| 日韩欧美精品v在线| 麻豆成人av在线观看| 麻豆一二三区av精品| 精品99又大又爽又粗少妇毛片 | 丰满的人妻完整版| 在线观看舔阴道视频| av黄色大香蕉| 欧美高清性xxxxhd video| 日韩 亚洲 欧美在线| 简卡轻食公司| 欧美成人性av电影在线观看| 婷婷色综合大香蕉| 精品国产亚洲在线| 伦理电影大哥的女人| 色尼玛亚洲综合影院| 两个人视频免费观看高清| 十八禁网站免费在线| 国产视频一区二区在线看| 小说图片视频综合网站| 亚洲第一区二区三区不卡| 日韩欧美国产在线观看| 嫁个100分男人电影在线观看| 欧美+亚洲+日韩+国产| 国产免费男女视频| 97超级碰碰碰精品色视频在线观看| 国产精品嫩草影院av在线观看 | 中文字幕高清在线视频| 两人在一起打扑克的视频| 欧美色视频一区免费| www.999成人在线观看| 婷婷六月久久综合丁香| 精品久久久久久,| 久久这里只有精品中国| 亚洲成av人片在线播放无| 中文字幕熟女人妻在线| 国内揄拍国产精品人妻在线| 色噜噜av男人的天堂激情| 精品国产三级普通话版| 国产精品一区二区性色av| 亚洲在线观看片| 首页视频小说图片口味搜索| 俄罗斯特黄特色一大片| 深夜精品福利| www日本黄色视频网| ponron亚洲| 国产黄片美女视频| 色综合婷婷激情| 级片在线观看| 国产成人影院久久av| 亚洲av电影在线进入| 国内精品美女久久久久久| 国产免费男女视频| 人人妻人人看人人澡| 国产 一区 欧美 日韩| 欧美黄色淫秽网站| 久久午夜亚洲精品久久| 久久亚洲真实| 最近中文字幕高清免费大全6 | 欧美激情久久久久久爽电影| 精品久久久久久久久久免费视频| 精品99又大又爽又粗少妇毛片 | 欧美色视频一区免费| 亚洲成av人片免费观看| 日韩av在线大香蕉| 亚洲国产欧美人成| 美女xxoo啪啪120秒动态图 | 毛片女人毛片| 欧美成人性av电影在线观看| 少妇的逼好多水| 免费电影在线观看免费观看| 一a级毛片在线观看| 亚洲乱码一区二区免费版| 亚洲中文字幕日韩| 看十八女毛片水多多多| 久久久久免费精品人妻一区二区| 国产在视频线在精品| 久久草成人影院| 免费在线观看成人毛片| 18禁裸乳无遮挡免费网站照片| 欧美成人性av电影在线观看| 国产精品自产拍在线观看55亚洲| 国产成年人精品一区二区| 99久久精品热视频| 亚洲乱码一区二区免费版| 亚洲av成人不卡在线观看播放网| 精品一区二区三区视频在线观看免费| 免费无遮挡裸体视频| 亚洲va日本ⅴa欧美va伊人久久| 中文字幕高清在线视频| 国产精品电影一区二区三区| 国产大屁股一区二区在线视频| 欧美成狂野欧美在线观看| 黄色视频,在线免费观看| 亚洲乱码一区二区免费版| 国产精品国产高清国产av| 天天一区二区日本电影三级| h日本视频在线播放| 免费av不卡在线播放| 国产一区二区三区视频了| 白带黄色成豆腐渣| 国产伦在线观看视频一区| 人妻丰满熟妇av一区二区三区| 午夜精品一区二区三区免费看| 99精品在免费线老司机午夜| 99久久精品国产亚洲精品| 亚洲精品久久国产高清桃花| 国产精品爽爽va在线观看网站| 波多野结衣高清作品| 久久久精品大字幕| 久久久久久久精品吃奶| 亚洲最大成人av| 综合色av麻豆| 亚洲av日韩精品久久久久久密| 一个人免费在线观看的高清视频| 国产精品av视频在线免费观看| 欧美在线黄色| 免费高清视频大片| 深爱激情五月婷婷| 亚洲av.av天堂| 久久亚洲精品不卡| 亚洲美女视频黄频| 亚洲国产精品999在线| 国产午夜精品久久久久久一区二区三区 | 国产一区二区三区视频了| 国产精品一区二区三区四区免费观看 | 国产精品野战在线观看| 在线播放国产精品三级| 欧美国产日韩亚洲一区| 精品久久久久久久久久久久久| 男人狂女人下面高潮的视频| 女人被狂操c到高潮| 久久99热6这里只有精品| 特大巨黑吊av在线直播| 成年女人毛片免费观看观看9| 99热这里只有精品一区| 成人精品一区二区免费| 欧美激情国产日韩精品一区| 日韩欧美在线乱码| 国产成+人综合+亚洲专区| 日本熟妇午夜| 亚洲经典国产精华液单 | 日本a在线网址| 69av精品久久久久久| 国产成人影院久久av| 麻豆成人午夜福利视频| x7x7x7水蜜桃| 亚洲第一欧美日韩一区二区三区| 久久伊人香网站| 精品国内亚洲2022精品成人| 亚洲中文日韩欧美视频| 欧美激情国产日韩精品一区| 精品福利观看| 一区二区三区激情视频| ponron亚洲| 少妇人妻一区二区三区视频| 黄色女人牲交| 99在线人妻在线中文字幕| 嫩草影院精品99| 免费人成视频x8x8入口观看| 国产亚洲av嫩草精品影院| 热99re8久久精品国产| 国产单亲对白刺激| 欧美国产日韩亚洲一区| 欧美最黄视频在线播放免费| 亚洲,欧美,日韩| 一个人看的www免费观看视频| 亚洲中文字幕日韩| 亚洲av成人不卡在线观看播放网| 精华霜和精华液先用哪个| 国产精品久久久久久精品电影| 99在线人妻在线中文字幕| 免费大片18禁| 中亚洲国语对白在线视频| 亚洲欧美日韩高清专用| 在线看三级毛片| 亚洲精品乱码久久久v下载方式| 欧美一级a爱片免费观看看| netflix在线观看网站| 国产精品国产高清国产av| 亚洲第一区二区三区不卡| 人人妻,人人澡人人爽秒播| 久久亚洲精品不卡| 有码 亚洲区| 超碰av人人做人人爽久久| 亚洲五月天丁香| 久9热在线精品视频| 首页视频小说图片口味搜索| 级片在线观看| 色综合婷婷激情| 亚洲av二区三区四区| 麻豆久久精品国产亚洲av| 美女大奶头视频| 亚洲精品色激情综合| 亚洲性夜色夜夜综合| 97热精品久久久久久| 非洲黑人性xxxx精品又粗又长| 内地一区二区视频在线| 欧美最黄视频在线播放免费| 757午夜福利合集在线观看| 亚洲电影在线观看av| 婷婷精品国产亚洲av在线| 色综合亚洲欧美另类图片| 亚洲无线在线观看| 久久午夜福利片| 很黄的视频免费| 午夜久久久久精精品| 午夜日韩欧美国产| 午夜福利高清视频| 精品无人区乱码1区二区| 亚洲久久久久久中文字幕| 99久久精品一区二区三区| 欧美不卡视频在线免费观看| 国产高清视频在线播放一区| 最新中文字幕久久久久| 女生性感内裤真人,穿戴方法视频| 51午夜福利影视在线观看| 国产亚洲av嫩草精品影院| 老司机午夜福利在线观看视频| 人妻丰满熟妇av一区二区三区| 在线十欧美十亚洲十日本专区| 日本黄大片高清| 成人三级黄色视频| 九色成人免费人妻av| 99国产精品一区二区蜜桃av| 成人性生交大片免费视频hd| 国产v大片淫在线免费观看| 成人一区二区视频在线观看| 老熟妇乱子伦视频在线观看| 成人国产一区最新在线观看| 欧美午夜高清在线| 性色av乱码一区二区三区2| 久久久国产成人免费| 老司机福利观看| 美女被艹到高潮喷水动态| 欧美最新免费一区二区三区 | 哪里可以看免费的av片| 最近最新免费中文字幕在线| 97超级碰碰碰精品色视频在线观看| 国产在线精品亚洲第一网站| 又黄又爽又免费观看的视频| 亚洲精品亚洲一区二区| 国产精品一及| 天美传媒精品一区二区| bbb黄色大片| 亚洲人成网站高清观看| 中文字幕av成人在线电影| 欧美精品国产亚洲| 一级黄片播放器| 在线播放国产精品三级| 国产亚洲精品av在线| 麻豆久久精品国产亚洲av| 高清日韩中文字幕在线| 欧美乱妇无乱码| 啦啦啦观看免费观看视频高清| 在线a可以看的网站| 亚洲精品一区av在线观看| 欧美日韩中文字幕国产精品一区二区三区| 色精品久久人妻99蜜桃| 免费在线观看成人毛片| 精品一区二区三区av网在线观看| 久久久久国产精品人妻aⅴ院| 91久久精品电影网| 久久久久久久久大av| 国内毛片毛片毛片毛片毛片| 一级作爱视频免费观看| www.熟女人妻精品国产| 两性午夜刺激爽爽歪歪视频在线观看| 999久久久精品免费观看国产| 亚洲成人中文字幕在线播放| 久久精品综合一区二区三区| 日韩欧美精品免费久久 | 成人毛片a级毛片在线播放| 日韩欧美精品v在线| 色综合站精品国产| 桃红色精品国产亚洲av| 中文字幕免费在线视频6| 国产伦在线观看视频一区| 我的老师免费观看完整版| 97超视频在线观看视频| 我要搜黄色片| 午夜精品久久久久久毛片777| 男女那种视频在线观看| 女同久久另类99精品国产91| 久久久成人免费电影| 精品久久久久久久末码| 尤物成人国产欧美一区二区三区| 欧美日韩黄片免| 国产亚洲欧美在线一区二区| 国产淫片久久久久久久久 | 久久久久久久亚洲中文字幕 | 99久久无色码亚洲精品果冻| 老司机福利观看| 国产精品国产高清国产av| 午夜福利在线观看免费完整高清在 | 九色国产91popny在线| 亚洲色图av天堂| 日本 av在线| 久久亚洲真实| 99久久成人亚洲精品观看| 白带黄色成豆腐渣| www.色视频.com| 欧美激情在线99| 国产aⅴ精品一区二区三区波| 色综合婷婷激情| 成人毛片a级毛片在线播放| 国产亚洲欧美在线一区二区| 97精品久久久久久久久久精品| 亚洲精华国产精华液的使用体验| av网站免费在线观看视频| 久久久久久国产a免费观看| 大香蕉97超碰在线| 亚洲人成网站高清观看| 亚洲成色77777| 亚洲欧洲国产日韩| 亚洲精品影视一区二区三区av| 久久精品国产鲁丝片午夜精品| 全区人妻精品视频| 久久久久久伊人网av| 99久久精品国产国产毛片| 高清毛片免费看| 纵有疾风起免费观看全集完整版| 有码 亚洲区| 日本黄色片子视频| 日韩国内少妇激情av| 我的女老师完整版在线观看| 亚洲在线观看片| 交换朋友夫妻互换小说| 少妇的逼好多水| 熟女人妻精品中文字幕| 亚洲人成网站在线观看播放| 亚洲精品视频女| 亚洲av电影在线观看一区二区三区 | 又爽又黄无遮挡网站| 国产一区二区三区综合在线观看 | 国产一级毛片在线| 国产成人精品福利久久| 少妇人妻 视频| 成人无遮挡网站| 国产在线一区二区三区精| 欧美最新免费一区二区三区| 国产成人免费观看mmmm| 久久久色成人| 国产在线男女| 26uuu在线亚洲综合色| 王馨瑶露胸无遮挡在线观看| 又黄又爽又刺激的免费视频.| 欧美xxⅹ黑人| 高清av免费在线| 国产美女午夜福利| 看黄色毛片网站| 亚洲精品久久久久久婷婷小说| 舔av片在线| 干丝袜人妻中文字幕| 亚洲最大成人中文| 国内揄拍国产精品人妻在线| 国产免费一区二区三区四区乱码| av黄色大香蕉| 国产精品一区二区性色av| 国产亚洲午夜精品一区二区久久 | 国产男女超爽视频在线观看| 丰满人妻一区二区三区视频av| 男人和女人高潮做爰伦理| 精品视频人人做人人爽| 欧美精品人与动牲交sv欧美| 国产av国产精品国产| 免费人成在线观看视频色| 2022亚洲国产成人精品| 啦啦啦在线观看免费高清www| 2021少妇久久久久久久久久久| 国产69精品久久久久777片| 成人午夜精彩视频在线观看| av播播在线观看一区| 91在线精品国自产拍蜜月| 久久精品国产亚洲av天美| 精品熟女少妇av免费看| 国产国拍精品亚洲av在线观看| 啦啦啦啦在线视频资源| 不卡视频在线观看欧美| 精品国产三级普通话版| 色婷婷久久久亚洲欧美| 日日摸夜夜添夜夜添av毛片| av免费在线看不卡| 日本黄色片子视频| av在线蜜桃| 亚洲国产精品成人综合色| 国产探花在线观看一区二区| 97人妻精品一区二区三区麻豆| 美女cb高潮喷水在线观看| 免费大片黄手机在线观看| 免费电影在线观看免费观看| 99久国产av精品国产电影| 国产精品国产av在线观看| 亚洲天堂av无毛| 男人和女人高潮做爰伦理| 不卡视频在线观看欧美| 久久久国产一区二区| 18禁动态无遮挡网站| 亚洲欧美精品自产自拍| 一级毛片久久久久久久久女| 极品教师在线视频| 九九在线视频观看精品| 久久久久久久午夜电影| 蜜桃亚洲精品一区二区三区| av在线蜜桃| 久久6这里有精品| 免费少妇av软件| 亚洲精品日韩av片在线观看| 欧美国产精品一级二级三级 | 午夜福利视频1000在线观看| 国产免费又黄又爽又色| 午夜福利在线观看免费完整高清在| 噜噜噜噜噜久久久久久91| 免费黄网站久久成人精品| 精品人妻偷拍中文字幕| 欧美成人一区二区免费高清观看| 精品久久久久久久人妻蜜臀av| 免费播放大片免费观看视频在线观看| videossex国产| 久久久欧美国产精品| 视频区图区小说| 丝袜脚勾引网站| 三级经典国产精品| 久久精品综合一区二区三区| 成年av动漫网址| 国产乱人视频| 国产免费一区二区三区四区乱码| 波多野结衣巨乳人妻| 草草在线视频免费看| 麻豆成人午夜福利视频| 欧美激情久久久久久爽电影| 国产亚洲91精品色在线| 晚上一个人看的免费电影| 人妻夜夜爽99麻豆av| 久久6这里有精品| 国产欧美日韩精品一区二区| 免费播放大片免费观看视频在线观看| 国产人妻一区二区三区在| 国产精品一区二区三区四区免费观看| 亚洲内射少妇av| 欧美高清性xxxxhd video| 男女边摸边吃奶| 一区二区三区乱码不卡18| 国产成人一区二区在线| 亚洲自偷自拍三级| 在线a可以看的网站| 国产精品秋霞免费鲁丝片| 一级黄片播放器| 熟女电影av网| 国产综合精华液| 乱码一卡2卡4卡精品| 日韩欧美精品v在线| 久久久久网色| 亚洲av二区三区四区| 在线观看国产h片| 国产成人精品一,二区| 国产一区二区在线观看日韩| 亚洲人与动物交配视频| 插逼视频在线观看| 黄色怎么调成土黄色| 99久久精品一区二区三区| 少妇人妻久久综合中文| 超碰97精品在线观看| 91久久精品电影网| av国产久精品久网站免费入址| 亚洲国产精品999| 又粗又硬又长又爽又黄的视频| 99热这里只有精品一区| 只有这里有精品99| 狂野欧美白嫩少妇大欣赏| 亚洲av欧美aⅴ国产| 日本欧美国产在线视频| 在线天堂最新版资源| 国产免费视频播放在线视频| 中国美白少妇内射xxxbb| 国产男女内射视频| 国产久久久一区二区三区| 免费高清在线观看视频在线观看| 春色校园在线视频观看| www.av在线官网国产| 嫩草影院精品99| 成人综合一区亚洲| 国产精品秋霞免费鲁丝片| 色视频www国产| 亚洲精品aⅴ在线观看| 国模一区二区三区四区视频| 国产伦在线观看视频一区| 有码 亚洲区| 真实男女啪啪啪动态图| 久久久久精品久久久久真实原创| 免费看日本二区| 天天躁夜夜躁狠狠久久av| 好男人在线观看高清免费视频| 国产老妇女一区| 嫩草影院新地址| 熟女av电影| 亚洲综合色惰| 亚洲图色成人| 亚洲天堂av无毛| 久久精品久久久久久久性| 国产极品天堂在线| 亚洲国产精品成人久久小说| 身体一侧抽搐| 久久久色成人| 国内少妇人妻偷人精品xxx网站| 97精品久久久久久久久久精品| 国产成人精品久久久久久| 日韩欧美精品免费久久| 人妻制服诱惑在线中文字幕| 国产精品精品国产色婷婷| 午夜免费观看性视频| 久久人人爽人人爽人人片va| 涩涩av久久男人的天堂| 在线观看一区二区三区| 狠狠精品人妻久久久久久综合| 亚洲自拍偷在线| 久久久久性生活片| 亚洲精品久久久久久婷婷小说| 国产 精品1| 国产成人精品久久久久久| 国产在视频线精品| 国产黄片美女视频| 国产高清有码在线观看视频| 搞女人的毛片| 精品久久久久久久久av| 乱码一卡2卡4卡精品| 国产成人免费观看mmmm| 极品少妇高潮喷水抽搐| 国产精品久久久久久精品电影| 亚洲不卡免费看| 亚洲怡红院男人天堂| 久久久久久久国产电影| 91aial.com中文字幕在线观看| 麻豆精品久久久久久蜜桃| 一区二区三区精品91| 久久精品久久久久久久性| 国产真实伦视频高清在线观看| 欧美丝袜亚洲另类| 婷婷色麻豆天堂久久| 国产精品一二三区在线看| 国产黄色免费在线视频| 五月天丁香电影| 国产精品一区二区性色av| 中文乱码字字幕精品一区二区三区| 国产免费福利视频在线观看| 一本一本综合久久| 1000部很黄的大片| 国产精品久久久久久精品电影小说 | 亚洲四区av| 国产有黄有色有爽视频| 美女cb高潮喷水在线观看| 大陆偷拍与自拍| 日本一二三区视频观看| 肉色欧美久久久久久久蜜桃 | 免费大片黄手机在线观看| 一区二区三区四区激情视频| 色综合色国产| 色播亚洲综合网| 99精国产麻豆久久婷婷| 婷婷色综合www| 亚洲欧美一区二区三区黑人 | 国产精品99久久99久久久不卡 | 极品教师在线视频| 精品国产乱码久久久久久小说| 九九爱精品视频在线观看| 国产成人91sexporn| 日本一本二区三区精品| 一边亲一边摸免费视频|