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

    空間電荷層效應(yīng)對(duì)固體氧化物燃料電池三相界面附近氧空位傳輸?shù)挠绊?

    2021-07-01 09:42:54徐晗張璐
    物理學(xué)報(bào) 2021年12期
    關(guān)鍵詞:界面效應(yīng)

    徐晗 張璐

    (西安交通大學(xué)建筑環(huán)境與能源應(yīng)用工程系, 西安 710049)

    1 引 言

    高溫運(yùn)行條件(800—1000 ℃)賦予固體氧化物燃料電池(solid oxide fuel cell, SOFC)發(fā)電效率高、污染物排放低以及燃料適用性廣等優(yōu)勢(shì), 同時(shí)也導(dǎo)致了電池成本高與耐久性差等一系列限制SOFC商業(yè)化的技術(shù)瓶頸, 因此運(yùn)行條件中低溫化(< 700 ℃)已成為SOFC商業(yè)化的必經(jīng)之路.然而運(yùn)行溫度的降低顯著影響SOFC電極與電解質(zhì)中的傳輸與反應(yīng)過(guò)程, 進(jìn)而增大極化阻抗, 導(dǎo)致電池電化學(xué)性能以指數(shù)規(guī)律下降.薄膜電解質(zhì)技術(shù)的快速發(fā)展顯著地降低了電解質(zhì)歐姆阻抗, 因此, 如何有效地降低電極極化阻抗已成為中低溫SOFC發(fā)展的關(guān)鍵瓶頸[1,2].

    在電極材料中加入電解質(zhì)材料形成復(fù)合電極,可顯著地增加電子導(dǎo)體相、離子導(dǎo)體相與氣相接觸形成的三相界面(three phase boundary, TPB)面積, 有效地降低整個(gè)電極電化學(xué)反應(yīng)引發(fā)的活化極化阻抗; 同時(shí), TPB也為電子導(dǎo)體與離子導(dǎo)體兩相材料組成的異質(zhì)界面, 可誘發(fā)空間電荷層(space charge layer, SCL)效應(yīng)[3,4]: 帶電點(diǎn)缺陷聚集使得異質(zhì)界面帶正電, 同時(shí)電荷守恒原理導(dǎo)致異質(zhì)界面附近形成帶負(fù)電的空間電荷層.帶電點(diǎn)缺陷的重新分布可使得TPB附近局部載流子濃度顯著區(qū)別于導(dǎo)體體相濃度, 改變界面附近電荷傳輸能力, 進(jìn)而促進(jìn)TPB處電化學(xué)反應(yīng)動(dòng)力學(xué), 降低反應(yīng)活化極化阻抗.尤其對(duì)于界面體積占比高的納米復(fù)合電極, SCL效應(yīng)的影響更為顯著, 例如文獻(xiàn)[5]實(shí)驗(yàn)研究發(fā)現(xiàn)電子導(dǎo)體納米薄膜修飾可顯著提升陰極性能, 而這與薄膜覆蓋阻礙氧離子傳輸理應(yīng)降低陰極性能的傳統(tǒng)認(rèn)知恰好相反, 他們猜測(cè)SCL效應(yīng)可能是該分歧產(chǎn)生的原因.可見(jiàn), 目前關(guān)于SCL效應(yīng)對(duì)納米復(fù)合電極性能的影響機(jī)理尚未達(dá)成定論,而深刻理解SCL效應(yīng)對(duì)TPB附近電荷傳輸過(guò)程的調(diào)控機(jī)理, 對(duì)降低反應(yīng)活化極化阻抗、提高中低溫條件下SOFC電極性能尤為重要.

    TPB特征尺度小且存在于電極內(nèi)部[6], 因此難以實(shí)驗(yàn)測(cè)量其附近的電荷傳輸過(guò)程, 相比之下數(shù)值模擬是研究TPB附近SCL效應(yīng)的有效手段.目前關(guān)于SOFC內(nèi)SCL效應(yīng)的模擬研究主要聚焦于摻雜釔的鋯酸鋇(BZY)等質(zhì)子導(dǎo)體[7-10]以及摻雜鈧的氧化鈰(SDC)、摻雜釓的氧化鈰(GDC)等CeO2基氧離子導(dǎo)體[11-16], 大多采用Mott-Schottky模型或Gouy-Chapman模型[3,17], 這兩個(gè)模型均源自于Poisson-Boltzmann (PB)方程.然而,PB方程基于載流子電化學(xué)平衡假設(shè), 簡(jiǎn)化SCL內(nèi)載流子濃度為Boltzmann分布, 該假設(shè)不適用于以下三種情況[18]: 1)載流子電化學(xué)勢(shì)梯度不為0(即具有宏觀運(yùn)動(dòng)導(dǎo)致的凈電流); 2)界面非均勻、離散帶電; 3)遠(yuǎn)離帶電界面處載流子濃度與其本體濃度不相等.受限于PB方程假設(shè), 現(xiàn)有針對(duì)SOFC內(nèi)SCL效應(yīng)的模擬研究普遍關(guān)注同相材料內(nèi)晶粒與晶粒之間形成的晶界, 尚未考慮TPB等不同相材料接觸形成的異質(zhì)界面; 重點(diǎn)關(guān)注燒結(jié)溫度等制備工藝的影響(凈電流密度近似為0, 電化學(xué)平衡假設(shè)近似成立), 尚未涉及運(yùn)行工況; 大多采用一維模型, 均未考慮非均勻離散帶電界面.例如文獻(xiàn)[7, 8]采用一維Mott-Schottky模型研究了不同釔的摻雜量與材料燒結(jié)溫度下BZY內(nèi)SCL效應(yīng)對(duì)晶界電阻的影響.Kim等[12,13,16]基于Mott-Schottky假設(shè)建立了預(yù)測(cè)晶界電勢(shì)的一維線性擴(kuò)散模型, 用于分析GDC內(nèi)晶界附近的SCL效應(yīng).綜上所述, 現(xiàn)有SOFC內(nèi)SCL效應(yīng)相關(guān)模擬研究尚未涉及運(yùn)行工況下SOFC電極內(nèi)TPB附近的載流子傳輸過(guò)程.

    作者前期耦合Poisson方程和載流子質(zhì)量守恒方程, 建立了適用于載流子具有宏觀運(yùn)動(dòng)時(shí)(電化學(xué)勢(shì)梯度不為零)氧離子導(dǎo)體內(nèi)晶界或者異質(zhì)界面附近載流子傳輸過(guò)程的數(shù)理模型, 基于一維模型假設(shè)研究了氧離子導(dǎo)體SCL內(nèi)氧空位傳輸機(jī)理[19].本文將基于上述數(shù)理模型, 結(jié)合SOFC模式電極, 進(jìn)一步考慮TPB分布的離散屬性, 發(fā)展適用于SOFC運(yùn)行工況的二維數(shù)理模型及格子玻爾茲曼數(shù)值模擬方法, 揭示SCL效應(yīng)對(duì)SOFC內(nèi)TPB附近電荷傳輸過(guò)程的影響規(guī)律及機(jī)理.

    2 求解問(wèn)題描述與數(shù)值方法

    借助SOFC模式電極構(gòu)筑幾何形狀與位置可控的非均勻離散TPB.圖1描述了本文研究中模式電極的幾何結(jié)構(gòu).如圖1(c)與圖1(d)中紅色實(shí)心點(diǎn)與紅色實(shí)線所示TPB位置, 此處發(fā)生電化學(xué)反應(yīng), 同時(shí)也為異質(zhì)界面, 導(dǎo)致離子導(dǎo)體內(nèi)該界面附近形成空間電荷層, 重點(diǎn)關(guān)注該區(qū)域的氧空位傳輸特性.考慮電極寬度(沿x軸方向)和高度(沿y軸方向)遠(yuǎn)小于電極長(zhǎng)度以及電解質(zhì)半徑, 可取二維x-y平面內(nèi)一個(gè)電極電解質(zhì)單元作為研究對(duì)象(如圖1(b)與圖1(c)); 同時(shí)由于每個(gè)電極電解質(zhì)單元為中心對(duì)稱(chēng)結(jié)構(gòu), 最終計(jì)算區(qū)域?yàn)閳D1(d)所示的SCL與體相(Bulk)區(qū)域.需要注意的是,本文重點(diǎn)關(guān)注TPB附近SCL內(nèi)的氧空位傳輸特性, 而TPB厚度與SCL厚度相比不可忽略, 因此不同于電池層面計(jì)算中可忽略TPB厚度, 本文計(jì)算中考慮TPB厚度lTPB(如圖1(d)所示).

    圖1 (a) SOFC模式電極幾何結(jié)構(gòu)示意圖; (b), (c), (d)本文計(jì)算區(qū)域, 包括TPB附近的SCL與體相(Bulk)區(qū)域Fig.1.(a) Schematic of a patterned SOFC electrode; (b),(c), (d) computational domain considered in the present study, including the SCL and bulk area adjacent to the TPB.

    采用文獻(xiàn)[19]中發(fā)展的泊松-載流子質(zhì)量守恒耦合方程描述TPB附近氧空位傳輸過(guò)程, 該耦合方程可準(zhǔn)確描述載流子具有宏觀運(yùn)動(dòng)時(shí)氧離子導(dǎo)體內(nèi)載流子傳輸特性.由于SOFC常用ZrO2基、CeO2基等螢石結(jié)構(gòu)氧化物陶瓷作為材料, 重點(diǎn)研究M2O3摻雜的AO2氧化物, 其內(nèi)部載流子包括陽(yáng)離子受體與氧空位(分別用下標(biāo)a與V表示),忽略僅在燒結(jié)溫度(遠(yuǎn)高于SOFC運(yùn)行溫度)下才會(huì)移動(dòng)的陽(yáng)離子受體濃度變化[3], 只考慮氧空位作為載流子.由導(dǎo)體體相電中性假設(shè)可得:ca(x,y)=因此, 控制電勢(shì)和氧空位濃度分布的耦合方程組可表示為

    式中,φ為電勢(shì),F為法拉第常數(shù),zV為氧空位攜帶電荷數(shù),c為載流子濃度,ε0與εr分別表示真空與相對(duì)介電常數(shù),D為擴(kuò)散系數(shù),t為時(shí)間.本文關(guān)注穩(wěn)態(tài)下氧空位傳輸特性, 故(2)式可進(jìn)一步簡(jiǎn)化為

    氧空位在濃度梯度與電勢(shì)梯度的共同驅(qū)動(dòng)下傳輸,且產(chǎn)生的電流密度為

    式中, 等號(hào)右邊第一項(xiàng)為濃度梯度導(dǎo)致的擴(kuò)散電流密度idif, 第二項(xiàng)為電勢(shì)梯度導(dǎo)致的遷移電流密度imig.

    SCL效應(yīng)導(dǎo)致的額外電阻RSCL(以下簡(jiǎn)稱(chēng)SCL電阻)可表示為

    式中,Rtot為考慮SCL效應(yīng)時(shí)的總電阻,Rbulk為相同計(jì)算與邊界條件下不考慮SCL效應(yīng)時(shí)的電阻.已有文獻(xiàn)基于(5)式給出了SCL電阻一維計(jì)算公式[2,19,20], 但無(wú)法適用于具有非均勻離散TPB的二維情況.參考磚塊模型中對(duì)晶界電阻的近似處理方式[20], 本文構(gòu)建圖2所示的二維電阻網(wǎng)絡(luò)圖, 據(jù)此可得以下表達(dá)式:

    圖2 二維電阻網(wǎng)絡(luò)圖Fig.2.Two-dimensional resistance network.

    式中,σV為電導(dǎo)率且Δl為y方向每一個(gè)離散電阻的長(zhǎng)度, Δx為x方向每一個(gè)離散電阻的寬度, Δz為垂直x-y平面的單位厚度.

    選取電解質(zhì)厚度l0、上下邊界電勢(shì)差φ0、氧空位體相濃度cV0分別為特征長(zhǎng)度、特征電勢(shì)與特征濃度, 對(duì)(1)式、(3)式與(4)式進(jìn)行無(wú)量化處理,獲得以下控制TPB附近氧空位傳輸過(guò)程的無(wú)量綱參數(shù):

    式中, 上標(biāo)*表示無(wú)量綱量;i0=4F2cV0Dφ0/(RTl0)為特征電流密度;λD為Debye長(zhǎng)度:

    根據(jù)圖1(d)所示計(jì)算區(qū)域與氧空位傳輸物理過(guò)程, 采用Neumann邊界條件描述左邊界與右邊界的電勢(shì)和氧空位濃度分布; 下邊界為體相區(qū)域,取該處電勢(shì)為參考電勢(shì), 氧空位濃度為體相濃度;上邊界非TPB區(qū)域采用不可滲透邊界; 上邊界TPB區(qū)域有氧空位移動(dòng), 故給定電流密度與電勢(shì).歸納本文數(shù)值模型的無(wú)量綱邊界條件列于表1.

    作者前期研究表明格子玻爾茲曼方法(lattice Boltzmann, LB)方法能夠很好地預(yù)測(cè)模式電極電化學(xué)特性以及氧離子導(dǎo)體空間電荷層內(nèi)載流子傳輸過(guò)程[19,21-24], 在文獻(xiàn)[19]中詳細(xì)推導(dǎo)了求解泊松-載流子質(zhì)量守恒耦合方程的LB模型, 并基于一維模型假設(shè)研究了氧離子導(dǎo)體SCL內(nèi)載流子傳輸特性, 該LB模型可直接用于本文二維模式電極求解, 只需修改邊界條件即可.本文采用非平衡外推格式[25]處理表1所列的邊界條件.

    表1 本文的邊界條件Table 1.Boundary conditions of the present study.

    3 結(jié)果分析

    3.1 模型驗(yàn)證

    將60 × 60, 120 × 120和240 × 240網(wǎng)格下的計(jì)算結(jié)果分別與480 × 480網(wǎng)格下計(jì)算結(jié)果進(jìn)行比較, 計(jì)算誤差分別為13.6%, 6.2%與2.6%, 因此可認(rèn)為當(dāng)網(wǎng)格數(shù)為240 × 240時(shí), 程序計(jì)算網(wǎng)格無(wú)關(guān)性已滿足.綜合考慮計(jì)算準(zhǔn)確性與效率, 后續(xù)計(jì)算網(wǎng)格數(shù)選為240 × 240.

    本文與文獻(xiàn)[19]采用相同的控制方程與LB模型, 涉及不同的邊界條件, 而邊界條件處理亦會(huì)影響模型準(zhǔn)確性.作者已在文獻(xiàn)[19]中對(duì)求解泊松-載流子質(zhì)量守恒耦合方程的LB模型進(jìn)行了驗(yàn)證, 本文進(jìn)一步驗(yàn)證模式電極相關(guān)邊界條件實(shí)施的可靠性.若不考慮SCL效應(yīng), 則?cV=0 , Poisson方程(1)可簡(jiǎn)化為L(zhǎng)aplace方程:

    作者前期發(fā)展并實(shí)驗(yàn)驗(yàn)證了求解方程(11)的LB模型, 同時(shí)應(yīng)用于模式陽(yáng)極內(nèi)電荷傳輸特性研究[21,22].采用上述經(jīng)過(guò)實(shí)驗(yàn)驗(yàn)證的求解Laplace方程的LB模型預(yù)測(cè)圖1(d)所示結(jié)構(gòu)中電勢(shì)分布,并采用計(jì)算所得的TPB區(qū)域電勢(shì)和電流密度分布作為本文泊松-載流子質(zhì)量守恒耦合方程LB模型的邊界條件, 若本文數(shù)值模型模式電極邊界條件實(shí)施正確, 則會(huì)模擬出與Laplace方程一樣的計(jì)算結(jié)果.如圖3所示, 耦合LB模型計(jì)算的電勢(shì)分布與不考慮SCL效應(yīng)的Laplace方程的計(jì)算結(jié)果完全相同, 且模擬獲得的氧空位濃度在整個(gè)計(jì)算區(qū)域?yàn)槎ㄖ?(x,y)=1 ), 有效地證明了本文模式電極邊界條件實(shí)施的正確性與可靠性.

    圖3 本文LB模型驗(yàn)證: 右邊界(x/l0 = 1)與上邊界(y/l0 =1)電勢(shì)分布Fig.3.Validation of the present LB model: Potential distributions at the right (x/l0 = 1) and top (y/l0 = 1) boundary.

    3.2 基準(zhǔn)工況分析

    從本節(jié)開(kāi)始研究不同無(wú)量綱平均電流密度、無(wú)量綱Debye長(zhǎng)度與無(wú)量綱電勢(shì)下SCL效應(yīng)對(duì)TPB附近氧空位傳輸過(guò)程的影響規(guī)律.基準(zhǔn)工況為:iav/i0= —0.4,λD/l0= 4.30 × 10—2,zVFφ0/(RT) =4.33 × 10—3.圖4描述了基準(zhǔn)工況下TPB附近氧空位傳輸特性.電勢(shì)與氧空位濃度均在TPB處具有最大值, 沿著遠(yuǎn)離TPB方向減小(如圖4(a)與圖4(b)); 電勢(shì)在靠近TPB區(qū)域放射狀減小, 在遠(yuǎn)離TPB區(qū)域其等高線逐漸平行于下邊界; 不同于不考慮SCL效應(yīng)時(shí)氧空位濃度在整個(gè)計(jì)算區(qū)域保持不變, 受SCL效應(yīng)影響, 氧空位濃度在靠近TPB區(qū)域劇烈變化, 遠(yuǎn)離TPB區(qū)域其數(shù)值為本體濃度保持不變.同時(shí)電勢(shì)、氧空位濃度等高線與上邊界非TPB區(qū)域、左邊界和右邊界垂直, 且遷移與擴(kuò)散電流均從TPB流入至底部(y/l0= 0)流出, 其流線方向與上邊界非TPB區(qū)域、左邊界和右邊界平行(如圖4(c)與圖4(d)).取idif,tot==0.001為SCL與本體區(qū)域分界線以確定SCL厚度(如圖4(b)中粉紅色實(shí)線所示),發(fā)現(xiàn)SCL近似為以TPB為中心的1/4圓形, 且厚度大于Debye厚度.

    圖4 基準(zhǔn)工況下TPB附近氧空位傳輸特性, 其中, 電勢(shì)(a)與氧空位濃度(b)分布; 遷移(c)與擴(kuò)散(d)電流密度大小及流線分布; 右邊界(x/l0 = 1)電勢(shì)、電勢(shì)梯度與氧空位濃度(e), 以及電流密度與電荷密度分布(f)Fig.4.Oxygen vacancy transport adjacent to the TPB under standard case: Potential (a) and oxygen vacancy concentration(b) distribution; migration (c) and diffusion (d) current density streamline; distributions of potential, potential gradient, oxygen vacancy concentration (e), current density and charge density (f) at x/l0 = 1.

    圖4 (e)與圖4(f)為計(jì)算區(qū)域右邊界(x/l0= 1)各物理量分布.沿著y/l0增大的方向, 電勢(shì)φ/φ0不斷增大且變化速度先增大后減小, 導(dǎo)致?(φ/φ0)呈現(xiàn)上凸的變化趨勢(shì), 且其數(shù)值在TPB附近的減小程度遠(yuǎn)大于氧空位濃度的增大程度, 由(4)式可知遷移電流密度imig由電勢(shì)梯度與氧空位濃度共同決定, 所以其絕對(duì)值先增大后減小; 隨著向TPB靠近, 氧空位濃度變化愈加劇烈, 導(dǎo)致?(cV/cV0)不斷增大, 因此擴(kuò)散電流密度idif絕對(duì)值呈現(xiàn)逐漸增大的變化趨勢(shì); 總電流度密度itot=imig+idif的絕對(duì)值也呈現(xiàn)增大的變化趨勢(shì).無(wú)量綱凈電荷密 度呈現(xiàn)和氧空位濃度類(lèi)似的變化趨勢(shì): 在體相其數(shù)值為零, 在SCL區(qū)域不斷增大, 同時(shí)可以發(fā)現(xiàn)由于F/ε0εr~ 1015, 極小的氧空位濃度變化會(huì)引起劇烈的凈電荷密度變化, 從而導(dǎo)致電勢(shì)分布顯著變化.同時(shí)基準(zhǔn)工況下, TPB中心點(diǎn)處imig與idif分別為—2.56與—5.44, SCL效應(yīng)導(dǎo)致的擴(kuò)散電流密度idif占據(jù)總電流密度的68%, 說(shuō)明此時(shí)SCL效應(yīng)對(duì)TPB附近氧空位傳輸特性具有重要影響.

    3.3 不同無(wú)量綱平均電流密度下空間電荷層效應(yīng)的影響規(guī)律

    圖5描述了不同無(wú)量綱平均電流密度(iav/i0,其中i0=4F2cV0DVφ0/(RTl0) )下SCL效應(yīng)對(duì)氧空位傳輸過(guò)程的影響規(guī)律.無(wú)量綱平均電流密度表征相同計(jì)算與邊界條件下, 考慮SCL效應(yīng)時(shí)的平均電流密度與不考慮SCL效應(yīng)時(shí)的平均電流密度的偏離程度, 可綜合考慮不同金屬氧化物和摻雜濃度引起的氧空位濃度cV0與擴(kuò)散系數(shù)DV、運(yùn)行溫度T、電解質(zhì)兩端過(guò)電勢(shì)φ0與厚度l0等參數(shù)的影響規(guī)律.圖5中iav/i0為負(fù)表示電流方向與y軸正方向相反; 反之則與y軸正方向同向.不考慮SCL效應(yīng)時(shí),iav/i0= —0.34.當(dāng)iav/i0< 0且iav/i0<—0.34時(shí), 沿著靠近TPB的方向, 電勢(shì)呈現(xiàn)先增大后減小的變化趨勢(shì)(圖5(c)), 導(dǎo)致遷移電流密度imig/i0在遠(yuǎn)離TPB區(qū)域小于0、靠近TPB區(qū)域大于0 (圖5(e)); 氧空位濃度cV/cV0則不斷增大(圖5(d)), 導(dǎo)致idif/i0< 0 (圖5(f)).隨著iav/i0的增大, TPB附近電勢(shì)梯度?(φ/φ0),cV/cV0, 以及?(cV/cV0)均呈現(xiàn)減小的變化規(guī)律, 導(dǎo)致該區(qū)域imig/i0與idif/i0絕對(duì)值均減小.當(dāng)iav/i0> 0時(shí), 呈現(xiàn)相反的物理量分布及變化規(guī)律.總體來(lái)說(shuō), 隨著iav/i0從負(fù)值增大到正值,φ/φ0的極值不斷減小,TPB中心點(diǎn)處cV/cV0與imig/i0不斷減小、idif/i0不斷增大; 由于氧空位濃度從大于本體濃度減小到小于本體濃度, TPB附近氧空位傳輸能力下降, 導(dǎo)致計(jì)算區(qū)域SCL電阻不斷增加(圖5(a)), 而SCL厚度隨iav/i0與不考慮SCL效應(yīng)時(shí)iav/i0的差值的絕對(duì)值的增大而增大(圖5(b)).其中,RSCL< 0說(shuō)明SCL效應(yīng)促進(jìn)氧空位傳輸, 反之阻礙氧空位傳輸.因此, 當(dāng)平均電流密度與不考慮SCL效應(yīng)時(shí)電流密度同向時(shí), 隨著iav/i0絕對(duì)值的增大, SCL效應(yīng)對(duì)氧空位傳輸?shù)拇龠M(jìn)作用越顯著; 當(dāng)平均電流密度與不考慮SCL效應(yīng)時(shí)電流密度反向時(shí), SCL效應(yīng)會(huì)增大氧空位傳輸阻力.

    圖5 不同無(wú)量綱平均電流密度(iav/i0)下SCL效應(yīng)的影響, 其中, 不同iav/i0下SCL電阻與TPB中心點(diǎn)氧空位濃度(a)、以及SCL厚度(b)分布; 當(dāng)iav/i0分別為—1.6, —1, 1, 1.6, 以及x/l0 = 1時(shí)的界面電勢(shì)(c)、氧空穴濃度(d)、遷移(e)與擴(kuò)散(f)電流密度分布Fig.5.Influences of SCL effect under different dimensionless average current densities (iav/i0): The SCL resistance, oxygen vacancy concentration at central TPB (a), and SCL thickness (b) under different iav/i0; distributions of potential (c), oxygen vacancy concentration (d), migration (e) and diffusion (f) current density when iav/i0 = —1.6, —1, 1 and 1.6 at x/l0 = 1.

    3.4 不同無(wú)量綱Debye長(zhǎng)度下空間電荷層效應(yīng)的影響規(guī)律

    Debye長(zhǎng)度可以表征SOFC電解質(zhì)材料AO2-M2O3中金屬氧化物類(lèi)型和摻雜濃度的影響[19].無(wú)量綱Debye長(zhǎng)度(λD/l0)表征SCL厚度與電解質(zhì)厚度相對(duì)大小, 可綜合考慮運(yùn)行溫度T、電解質(zhì)厚度l0、不同金屬氧化物和摻雜濃度引起的不同介電常數(shù)ε0εr與氧空位濃度cV0等參數(shù)的影響規(guī)律.這里0.5.如圖6(a)所示,λD/l0對(duì)TPB中心點(diǎn)氧空位濃度cV/cV0影響并不顯著, 當(dāng)λD/l0較大時(shí),cV/cV0沿著靠近TPB方向在更大幾何空間內(nèi)更為緩和地減小到1, 意味著靠近TPB區(qū)域的?(cV/cV0)更小(圖6(d)), 導(dǎo)致此處idif/i0絕對(duì)值更小, 而遠(yuǎn)離TPB區(qū)域idif/i0絕對(duì)值更大(圖6(f)).電荷密度分布呈現(xiàn)和cV/cV0相同的變化趨勢(shì), 隨著λD/l0的增大, 會(huì)沿著y/l0反方向, 在更大厚度內(nèi)更為緩慢地減小至0, 由Poisson方程(1)可知, 導(dǎo)致電勢(shì)分布曲率更小(圖6(c)).隨著λD/l0的增大, SCL厚度大于Debye厚度且不斷增大, SCL電阻不斷減小, 說(shuō)明隨著λD/l0的增大, SCL效應(yīng)對(duì)氧空位傳輸?shù)拇龠M(jìn)作用增強(qiáng).另外, 對(duì)比圖5可知, 改變無(wú)量綱Debye長(zhǎng)度時(shí), SCL效應(yīng)對(duì)氧空位傳輸過(guò)程的影響程度變化并不顯著.

    圖6 不同無(wú)量綱Debye長(zhǎng)度(λD/l0)下SCL效應(yīng)的影響, 其中, 不同λD/l0下SCL電阻與TPB中心點(diǎn)氧空位濃度(a)、以及SCL厚度(b)分布; 當(dāng)λD/l0分別為0.005, 0.01, 0.05, 以及x/l0 = 1時(shí)的界面電勢(shì)(c)、氧空穴濃度(d)、遷移(e)與擴(kuò)散(f)電流密度分布Fig.6.Influences of SCL effect under different dimensionless Debye length (λD/l0): The SCL resistance, oxygen vacancy concentration at central TPB (a), and SCL thickness (b) under different λD/l0; distributions of potential (c), oxygen vacancy concentration(d), migration (e) and diffusion (f) current density when λD/l0 = 0.005, 0.01 and 0.05 at x/l0 = 1.

    3.5 不同無(wú)量綱電勢(shì)下空間電荷層效應(yīng)的影響規(guī)律

    圖7 描述了不同無(wú)量綱電勢(shì)(zVFφ0/(RT))下SCL效應(yīng)對(duì)TPB附近氧空位傳輸特性的影響規(guī)律.zVFφ0/(RT)表征驅(qū)動(dòng)導(dǎo)體內(nèi)氧空位傳輸過(guò)程的過(guò)電勢(shì)與熱勢(shì)的相對(duì)重要性.如圖7(c), 隨著zVFφ0/(RT)的增大, 在右邊界靠近TPB附近,電勢(shì)分布曲率增大, 由Poisson方程可得, 電荷密度也會(huì)增大, 從而導(dǎo)致氧空位濃度cV/cV0增大且?(cV/cV0)顯著增加(圖7(a)與圖7(d)), 伴隨著擴(kuò)散電流流密度idif/i0絕對(duì)值不斷增大(圖7(f)),同時(shí)遷移電流密度imig/i0也呈現(xiàn)增大的變化趨勢(shì)(圖7(e)).同時(shí)由圖7(a)與圖7(b)可得, 增大zVFφ0/(RT)會(huì)導(dǎo)致SCL電阻與厚度均減小, 說(shuō)明隨著zVFφ0/(RT)的增大, SCL效應(yīng)對(duì)氧空位傳輸?shù)拇龠M(jìn)作用越顯著.同時(shí)發(fā)現(xiàn)當(dāng)zVFφ0/(RT) ≤10—3時(shí),zVFφ0/(RT)幾乎不影響氧空位傳輸過(guò)程;zVFφ0/(RT) > 10—3時(shí), SCL電阻對(duì)zVFφ0/(RT)的敏感性顯著增加.

    圖7 不同無(wú)量綱電勢(shì)(zVFφ0/(RT ))下SCL效應(yīng)的影響, 其中, 不同zVFφ0/(RT )下SCL電阻與TPB中心點(diǎn)氧空位濃度(a)、以及SCL厚度(b)分布; 當(dāng)zVFφ0/(RT )分別為0.001, 0.01, 0.1, 以及x/l0 = 1時(shí)的界面電勢(shì)(c)、氧空穴濃度(d)、遷移(e)與擴(kuò)散(f)電流密度分布Fig.7.Influences of SCL effect under different dimensionless potential (zVFφ0/(RT )): The SCL resistance, oxygen vacancy concentration at central TPB (a), and SCL thickness (b) under different zVFφ0/(RT ); distributions of potential (c), oxygen vacancy concentration (d), migration (e) and diffusion (f) current density when zVFφ0/(RT ) = 0.001, 0.01 and 0.1 at x/l0 = 1.

    4 討 論

    納米復(fù)合電極通常采用具有高電催化活性的納米結(jié)構(gòu)[4,26,27](納米顆粒、薄膜或者顆粒-薄膜復(fù)合結(jié)構(gòu))修飾傳統(tǒng)多孔電極骨架.以納米復(fù)合陰極為例, 孔相中氧氣遷移至骨架表面電化學(xué)活性位處, 與導(dǎo)體相中傳輸過(guò)來(lái)的電子發(fā)生氧氣還原反應(yīng), 生成氧離子并通過(guò)導(dǎo)體相傳輸至電解質(zhì).顯然,電極性能由氣體和載流子(電子與氧離子)傳輸與電化學(xué)反應(yīng)的相互耦合及競(jìng)爭(zhēng)決定.納米結(jié)構(gòu)修飾會(huì)顯著改變傳輸與反應(yīng)的耦合過(guò)程, 進(jìn)而影響電極性能: 一方面, 納米顆粒顯著增加骨架表面電化學(xué)活性位的面積與電催化活性, 從而提高表面電化學(xué)反應(yīng)能力, 納米顆粒尺寸越小、數(shù)量越多, 表面反應(yīng)能力的提升效果越明顯; 另一方面, 納米結(jié)構(gòu)的覆蓋可能增大載流子和氣體在反應(yīng)界面附近的傳輸阻力, 納米顆粒尺寸越小、數(shù)量越多時(shí), 導(dǎo)致電化學(xué)活性位附近孔隙尺寸越小, 其對(duì)電荷和氣體傳輸?shù)淖璧K作用越顯著.另外, 納米結(jié)構(gòu)可與多孔骨架接觸形成異質(zhì)界面, 并誘發(fā)SCL效應(yīng), 顯著改變導(dǎo)體內(nèi)異質(zhì)界面附近的載流子傳輸能力[28].對(duì)SOFC而言, SCL不僅存在于TPB附近, 也存在于電極與電解質(zhì)接觸構(gòu)筑的異質(zhì)界面處.在TPB處, SCL效應(yīng)影響穿透界面的載流子傳輸能力; 在電極與電解質(zhì)交界面, SCL效應(yīng)被認(rèn)為可能提高沿著界面方向的載流子傳輸能力[29,30].

    以SOFC常用的離子導(dǎo)體材料YSZ為例, 在典型運(yùn)行溫度下, 本文耦合模型計(jì)算的SCL厚度可達(dá)1 nm.可見(jiàn), 當(dāng)采用納米尺寸的顆?;蛘弑∧ば揎梻鹘y(tǒng)多孔電極形成納米復(fù)合電極時(shí), SCL厚度與顆粒或者薄膜尺寸在同一量級(jí), 此時(shí)異質(zhì)界面誘發(fā)的SCL效應(yīng)對(duì)整個(gè)電極性能的影響極為顯著,且需要詳細(xì)考慮SCL內(nèi)的載流子傳輸過(guò)程; 而當(dāng)采用微米尺寸顆粒修飾形成傳統(tǒng)復(fù)合電極時(shí), 約1 nm的SCL厚度在整個(gè)微米顆粒中占比小, 整個(gè)導(dǎo)體的導(dǎo)電性能仍然由本體區(qū)域載流子傳輸過(guò)程決定[31], 或者采用界面電阻描述SCL效應(yīng)影響即可.現(xiàn)有文獻(xiàn)使用的一維SCL厚度lMS預(yù)測(cè)模型為[3]:表示異質(zhì)界面處無(wú)量綱電勢(shì)).采用該公式計(jì)算典型運(yùn)行溫度下SOFC內(nèi)YSZ離子導(dǎo)體中SCL厚度約為0.3 nm, 同樣證明了SCL效應(yīng)對(duì)納米復(fù)合電極的作用相比于傳統(tǒng)復(fù)合電極更加重要.然而, 已有預(yù)測(cè)模型比本文耦合模型計(jì)算的SCL厚度小, 其根本原因在于其采用的載流子電化學(xué)平衡假設(shè)只適用于導(dǎo)體凈電流密度為零的一維問(wèn)題, 應(yīng)用于本文SOFC運(yùn)行工況下(導(dǎo)體凈電流密度不為零)的離散異質(zhì)界面二維分析會(huì)低估SCL效應(yīng)的影響.本文借助幾何結(jié)構(gòu)可控、反應(yīng)活性位和物質(zhì)傳輸路徑明確的模式電極[21,22], 構(gòu)造離散TPB異質(zhì)界面, 發(fā)展適用于導(dǎo)體凈電流密度不為零的Poisson方程與載流子質(zhì)量守恒方程耦合模型, 聚焦SOFC運(yùn)行工況下離子導(dǎo)體內(nèi)TPB附近的氧空位傳輸, 并著重分析SCL效應(yīng)的影響規(guī)律.結(jié)合本文發(fā)展的SCL數(shù)理模型與作者前期發(fā)展的SOFC多孔電極內(nèi)氣體、電荷傳輸與電化學(xué)反應(yīng)耦合模型[19,21-24], 可直接研究真實(shí)納米復(fù)合電極內(nèi)TPB、電極-電解質(zhì)交界面這兩類(lèi)異質(zhì)界面附近氧空位傳輸特性, 且所得研究結(jié)論也有助于納米復(fù)合電極的設(shè)計(jì)與運(yùn)行優(yōu)化策略研究.

    文獻(xiàn)[32]在LSC113表面修飾LSC214納米顆?;蛘呒{米薄膜, 發(fā)現(xiàn)電極氧氣還原反應(yīng)能力可提升3—4個(gè)數(shù)量級(jí), 并指出SCL效應(yīng)是電極反應(yīng)能力提升的關(guān)鍵原因之一, 而本文的研究結(jié)論恰好可指導(dǎo)相關(guān)電極材料、制備和運(yùn)行工藝的選擇和優(yōu)化.例如可通過(guò)改變電極材料、幾何結(jié)構(gòu)(例如電極厚度l0等)與運(yùn)行參數(shù)(例如運(yùn)行溫度T、電壓φ0等), 構(gòu)造不同的無(wú)量綱平均電流密度(iav/i0, 其中i0= 4F2cV0DVφ0/(RTl0))、無(wú)量綱Debye長(zhǎng)度(λD/l0, 其中與無(wú)量綱電勢(shì)(zVFφ0/(RT)), 從而達(dá)到理性設(shè)計(jì)和調(diào)控電極性能的目的, 其中電極材料的改變可通過(guò)謹(jǐn)慎選擇金屬氧化物、摻雜物和調(diào)整摻雜濃度(改變電解質(zhì)介電常數(shù)ε0εr、導(dǎo)體氧空位擴(kuò)散系數(shù)DV與氧空位體相濃度cV0等)實(shí)現(xiàn).需要說(shuō)明的是, 由于氧離子導(dǎo)體內(nèi)氧空位擴(kuò)散系數(shù)與氧離子躍遷頻率、一次躍遷距離和躍遷活化能相關(guān)[33], 選擇不同的金屬氧化物與摻雜物, 可導(dǎo)致不同的躍遷距離與活化能, 且不同的運(yùn)行溫度也會(huì)影響活化能的大小, 因此通過(guò)選擇不同的金屬氧化物和摻雜物、以及調(diào)整運(yùn)行溫度均可改變氧空位擴(kuò)散系數(shù).另外, 通過(guò)合理選擇上述電極設(shè)計(jì)與運(yùn)行參數(shù), 使得無(wú)量綱平均電流密度減小、無(wú)量綱Debye長(zhǎng)度以及無(wú)量綱電勢(shì)增大, 均可有效地利用異質(zhì)界面效應(yīng)提高SCL內(nèi)的氧離子傳輸能力.同時(shí), 由第3節(jié)的研究結(jié)論可知, 無(wú)量綱平均電流密度對(duì)TPB附近SCL內(nèi)氧空位傳輸過(guò)程影響最為顯著、無(wú)量綱電勢(shì)的影響次之、無(wú)量綱Debye的影響最小, 因此應(yīng)注意前兩個(gè)無(wú)量綱參數(shù)的合理構(gòu)造能夠更為顯著地調(diào)控納米復(fù)合電極性能.

    5 結(jié) 論

    納米復(fù)合電極是提高中低溫SOFC電化學(xué)性能的有效技術(shù)手段.兩相材料構(gòu)成復(fù)合電極可顯著地增大TPB面積, 同時(shí)可能在TPB處構(gòu)筑異質(zhì)界面引發(fā)SCL效應(yīng), 顯著影響SOFC電極內(nèi)載流子傳輸能力.現(xiàn)有SOFC運(yùn)行工況下的性能模擬分析則尚未考慮異質(zhì)界面導(dǎo)致的SCL效應(yīng); 同時(shí),關(guān)于SOFC內(nèi)SCL效應(yīng)的模擬研究均是采用基于載流子電化學(xué)平衡假設(shè)(即凈電流為零)的PB方程, 聚焦同相材料晶粒與晶粒間形成的晶界附近的SCL效應(yīng), 重點(diǎn)關(guān)注燒結(jié)溫度等制備工藝的影響規(guī)律, 尚未涉及運(yùn)行條件下(載流子存在電化學(xué)勢(shì)梯度導(dǎo)致的宏觀運(yùn)行, 即凈電流不為零) SOFC電極內(nèi)TPB (不同相材料組成的異質(zhì)界面)附近的載流子傳輸過(guò)程, 也尚未考慮異質(zhì)界面的離散性.

    基于SOFC模式電極構(gòu)筑非均勻離散TPB幾何結(jié)構(gòu), 本文建立Poisson方程與載流子質(zhì)量守恒方程耦合的數(shù)理模型, 模擬研究了運(yùn)行條件下電極TPB附近氧空位傳輸過(guò)程, 重點(diǎn)揭示了SCL效應(yīng)的影響規(guī)律及機(jī)理.模擬研究表明, 基準(zhǔn)工況下,由SCL效應(yīng)導(dǎo)致的擴(kuò)散電流密度占據(jù)總電流密度的68%, 證明了SCL效應(yīng)對(duì)TPB附近氧空位傳輸?shù)闹匾?通過(guò)研究不同無(wú)量綱平均電流密度、無(wú)量綱Debye長(zhǎng)度與無(wú)量綱電勢(shì)下SCL效應(yīng)的影響規(guī)律, 獲得以下研究結(jié)論: 1)增大無(wú)量綱平均電流密度(iav/i0), SCL效應(yīng)導(dǎo)致的氧空位傳輸阻力會(huì)增大, 且SCL厚度隨iav/i0絕對(duì)值的增大而增大;2)較大的無(wú)量綱Debye長(zhǎng)度可減小SCL電阻、增大SCL厚度, 即更加促進(jìn)氧空位傳輸; 3) SCL效應(yīng)對(duì)無(wú)量綱電勢(shì)(zVFφ0/(RT))的敏感性隨zVFφ0/(RT)的增加而增大; 當(dāng)zVFφ0/(RT) ≤ 10—3, 其數(shù)值幾乎不影響SCL電阻, 當(dāng)zVFφ0/(RT) ≥ 10—2,增大zVFφ0/(RT)可減小氧空位傳輸阻力與厚度.上述研究結(jié)論可為通過(guò)調(diào)控電極內(nèi)TPB附近的SCL效應(yīng)、以提高中低溫SOFC電極性能的納米復(fù)合電極技術(shù)的發(fā)展提供指導(dǎo)與理論依據(jù).

    猜你喜歡
    界面效應(yīng)
    鈾對(duì)大型溞的急性毒性效應(yīng)
    懶馬效應(yīng)
    場(chǎng)景效應(yīng)
    國(guó)企黨委前置研究的“四個(gè)界面”
    基于FANUC PICTURE的虛擬軸坐標(biāo)顯示界面開(kāi)發(fā)方法研究
    空間界面
    金秋(2017年4期)2017-06-07 08:22:16
    應(yīng)變效應(yīng)及其應(yīng)用
    電子顯微打開(kāi)材料界面世界之門(mén)
    人機(jī)交互界面發(fā)展趨勢(shì)研究
    手機(jī)界面中圖形符號(hào)的發(fā)展趨向
    新聞傳播(2015年11期)2015-07-18 11:15:04
    欧美日韩精品网址| 人人妻人人澡人人看| 国产黄色视频一区二区在线观看| 亚洲欧洲日产国产| 久久精品久久久久久噜噜老黄| 亚洲欧美精品自产自拍| 亚洲男人天堂网一区| 丝袜喷水一区| 黄色视频在线播放观看不卡| 亚洲精品国产色婷婷电影| av福利片在线| 满18在线观看网站| 欧美黑人欧美精品刺激| 丁香六月天网| a级毛片黄视频| 人人澡人人妻人| 欧美国产精品一级二级三级| 街头女战士在线观看网站| 十八禁网站网址无遮挡| 老司机亚洲免费影院| 国产一区亚洲一区在线观看| 欧美xxⅹ黑人| 国产又色又爽无遮挡免| 日韩视频在线欧美| 国产精品久久久久久久久免| 水蜜桃什么品种好| 欧美日韩视频高清一区二区三区二| 欧美中文综合在线视频| 午夜av观看不卡| 日韩av不卡免费在线播放| 亚洲精品aⅴ在线观看| 伊人亚洲综合成人网| 叶爱在线成人免费视频播放| 性色av一级| 日韩成人av中文字幕在线观看| 国产视频首页在线观看| 99精国产麻豆久久婷婷| 国产精品一二三区在线看| 丰满少妇做爰视频| 男女下面插进去视频免费观看| 久久精品熟女亚洲av麻豆精品| 成人影院久久| 制服诱惑二区| 色94色欧美一区二区| 日韩一区二区视频免费看| 国产av国产精品国产| 成人18禁高潮啪啪吃奶动态图| 日韩制服骚丝袜av| 国产黄色免费在线视频| 色94色欧美一区二区| 宅男免费午夜| 伊人久久国产一区二区| 亚洲欧美日韩另类电影网站| 老司机亚洲免费影院| 天美传媒精品一区二区| 爱豆传媒免费全集在线观看| 国产精品偷伦视频观看了| 亚洲第一区二区三区不卡| 久久久国产成人精品二区| 午夜亚洲福利在线播放| 黑人操中国人逼视频| 久久久久久国产a免费观看| 嫁个100分男人电影在线观看| 亚洲午夜理论影院| 又大又爽又粗| 美女免费视频网站| 国产1区2区3区精品| 1024视频免费在线观看| 精品一品国产午夜福利视频| 亚洲av成人不卡在线观看播放网| 美国免费a级毛片| 久久九九热精品免费| 18禁美女被吸乳视频| 精品久久久久久久毛片微露脸| 日本精品一区二区三区蜜桃| 亚洲欧美日韩无卡精品| 亚洲国产高清在线一区二区三 | 黄片大片在线免费观看| 黑丝袜美女国产一区| 亚洲av美国av| 国产精品久久久久久精品电影 | 不卡一级毛片| 激情在线观看视频在线高清| 国产伦一二天堂av在线观看| 欧美一级a爱片免费观看看 | 国产野战对白在线观看| 亚洲五月色婷婷综合| 国产亚洲精品av在线| 国产精品野战在线观看| 淫妇啪啪啪对白视频| 久久久久国产一级毛片高清牌| 99久久国产精品久久久| 亚洲人成电影观看| 久久精品91无色码中文字幕| 午夜福利视频1000在线观看 | 久久影院123| svipshipincom国产片| 在线观看免费午夜福利视频| 久久久久国内视频| 久久久久久久精品吃奶| 久久人妻福利社区极品人妻图片| 老汉色av国产亚洲站长工具| 两个人免费观看高清视频| 午夜福利18| 午夜影院日韩av| 韩国精品一区二区三区| 激情视频va一区二区三区| 韩国av一区二区三区四区| 亚洲精品av麻豆狂野| 精品久久久久久久久久免费视频| 精品久久久久久久人妻蜜臀av | 精品久久久久久久久久免费视频| 亚洲avbb在线观看| 亚洲第一欧美日韩一区二区三区| 黄色片一级片一级黄色片| 后天国语完整版免费观看| 久久精品成人免费网站| 人成视频在线观看免费观看| 亚洲狠狠婷婷综合久久图片| av有码第一页| 国产精品一区二区三区四区久久 | 亚洲中文日韩欧美视频| 国产高清videossex| 免费人成视频x8x8入口观看| 欧美日韩中文字幕国产精品一区二区三区 | 午夜久久久在线观看| 国产精品亚洲av一区麻豆| 可以免费在线观看a视频的电影网站| 免费不卡黄色视频| 可以在线观看的亚洲视频| 国产欧美日韩一区二区精品| 午夜免费激情av| 一个人观看的视频www高清免费观看 | 91麻豆av在线| 亚洲成a人片在线一区二区| 一进一出抽搐动态| 午夜福利视频1000在线观看 | 又黄又爽又免费观看的视频| www日本在线高清视频| 老司机深夜福利视频在线观看| 免费在线观看视频国产中文字幕亚洲| 精品久久久精品久久久| 香蕉久久夜色| 国产精品亚洲美女久久久| 国产成+人综合+亚洲专区| av电影中文网址| 美女扒开内裤让男人捅视频| 无人区码免费观看不卡| 国产精品久久视频播放| 亚洲无线在线观看| 精品国产美女av久久久久小说| 日韩成人在线观看一区二区三区| 老司机在亚洲福利影院| 精品国产国语对白av| 日韩欧美国产在线观看| 欧美色欧美亚洲另类二区 | x7x7x7水蜜桃| 后天国语完整版免费观看| 免费看十八禁软件| 国产伦人伦偷精品视频| 黄频高清免费视频| 国产精品久久久久久人妻精品电影| 久久九九热精品免费| 大陆偷拍与自拍| 满18在线观看网站| 国产成人精品久久二区二区91| 最近最新中文字幕大全免费视频| 俄罗斯特黄特色一大片| 91av网站免费观看| 两个人视频免费观看高清| 两个人视频免费观看高清| 日韩av在线大香蕉| 一级片免费观看大全| 久久亚洲真实| 两性夫妻黄色片| 十分钟在线观看高清视频www| 亚洲国产中文字幕在线视频| 制服人妻中文乱码| 久久人人爽av亚洲精品天堂| 精品第一国产精品| 亚洲色图综合在线观看| 国产亚洲精品综合一区在线观看 | 欧美午夜高清在线| 一进一出好大好爽视频| 在线天堂中文资源库| 变态另类丝袜制服| 精品国产乱码久久久久久男人| 日韩视频一区二区在线观看| 国产亚洲精品久久久久5区| 免费看十八禁软件| 国产极品粉嫩免费观看在线| 午夜久久久在线观看| 久久国产亚洲av麻豆专区| 一区二区三区激情视频| 波多野结衣高清无吗| 久久久国产欧美日韩av| 男人舔女人下体高潮全视频| 欧美激情极品国产一区二区三区| 亚洲av熟女| 男人的好看免费观看在线视频 | 欧美日韩精品网址| 老汉色∧v一级毛片| 国产野战对白在线观看| 一区二区日韩欧美中文字幕| 在线免费观看的www视频| 19禁男女啪啪无遮挡网站| 90打野战视频偷拍视频| 亚洲狠狠婷婷综合久久图片| 一级片免费观看大全| 69精品国产乱码久久久| av福利片在线| 亚洲第一欧美日韩一区二区三区| 级片在线观看| 黄色毛片三级朝国网站| 搡老熟女国产l中国老女人| 国语自产精品视频在线第100页| 欧美精品亚洲一区二区| 色哟哟哟哟哟哟| 免费看a级黄色片| 在线永久观看黄色视频| 成人亚洲精品一区在线观看| 免费人成视频x8x8入口观看| 久9热在线精品视频| 中文字幕色久视频| 熟女少妇亚洲综合色aaa.| 国产片内射在线| 国产99久久九九免费精品| 老司机深夜福利视频在线观看| 亚洲 欧美 日韩 在线 免费| 国产av在哪里看| 亚洲视频免费观看视频| 国产三级黄色录像| av网站免费在线观看视频| 真人一进一出gif抽搐免费| 亚洲成国产人片在线观看| 久久人妻av系列| 亚洲专区字幕在线| 精品不卡国产一区二区三区| 好男人电影高清在线观看| 国产私拍福利视频在线观看| av在线天堂中文字幕| 一级,二级,三级黄色视频| 欧美成狂野欧美在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 一级毛片女人18水好多| www.精华液| 夜夜躁狠狠躁天天躁| 亚洲国产精品合色在线| 精品电影一区二区在线| 黄色女人牲交| 极品教师在线免费播放| 国产亚洲av嫩草精品影院| 老司机午夜十八禁免费视频| 高清在线国产一区| 色播在线永久视频| 久热爱精品视频在线9| 亚洲伊人色综图| 亚洲视频免费观看视频| 黑人巨大精品欧美一区二区mp4| 不卡av一区二区三区| 欧美日本视频| 欧美日韩瑟瑟在线播放| 在线观看66精品国产| 亚洲成人免费电影在线观看| 侵犯人妻中文字幕一二三四区| 久热这里只有精品99| 少妇的丰满在线观看| 免费高清视频大片| 欧美日韩黄片免| 伦理电影免费视频| 国产亚洲欧美98| 琪琪午夜伦伦电影理论片6080| 伊人久久大香线蕉亚洲五| 侵犯人妻中文字幕一二三四区| 亚洲va日本ⅴa欧美va伊人久久| 99国产精品99久久久久| 国产日韩一区二区三区精品不卡| 欧美日韩亚洲国产一区二区在线观看| 久久久精品国产亚洲av高清涩受| 欧美日韩福利视频一区二区| 1024视频免费在线观看| 露出奶头的视频| 国产麻豆成人av免费视频| 欧美日韩瑟瑟在线播放| 欧美激情极品国产一区二区三区| 国产精品 欧美亚洲| 成人av一区二区三区在线看| 久久人妻熟女aⅴ| 欧美一级a爱片免费观看看 | 1024视频免费在线观看| 女性被躁到高潮视频| 国产一区二区三区综合在线观看| 亚洲国产中文字幕在线视频| 日本五十路高清| 又大又爽又粗| 亚洲va日本ⅴa欧美va伊人久久| 国产一区二区在线av高清观看| 亚洲国产中文字幕在线视频| 美女免费视频网站| 极品人妻少妇av视频| 国产亚洲精品久久久久久毛片| 在线免费观看的www视频| 欧美中文日本在线观看视频| 纯流量卡能插随身wifi吗| 久久久久久久精品吃奶| 韩国av一区二区三区四区| 免费少妇av软件| 午夜亚洲福利在线播放| 两个人看的免费小视频| 搡老岳熟女国产| 欧美乱妇无乱码| 亚洲专区字幕在线| 久热这里只有精品99| 十八禁人妻一区二区| 成人永久免费在线观看视频| 国产精品香港三级国产av潘金莲| 禁无遮挡网站| 国产精品98久久久久久宅男小说| 日韩欧美免费精品| 男男h啪啪无遮挡| 丰满的人妻完整版| 色综合欧美亚洲国产小说| 18禁观看日本| 国产av又大| 国产乱人伦免费视频| 黄色视频,在线免费观看| 精品一区二区三区av网在线观看| 国产av一区二区精品久久| 国产亚洲欧美98| 日本免费一区二区三区高清不卡 | 午夜视频精品福利| 国产精华一区二区三区| 亚洲欧美精品综合一区二区三区| 亚洲av第一区精品v没综合| 中文字幕av电影在线播放| 免费av毛片视频| 狂野欧美激情性xxxx| 精品电影一区二区在线| 真人一进一出gif抽搐免费| 97人妻天天添夜夜摸| 俄罗斯特黄特色一大片| 国产成人精品久久二区二区免费| 长腿黑丝高跟| 9色porny在线观看| 丰满的人妻完整版| 欧美国产日韩亚洲一区| 久久狼人影院| av片东京热男人的天堂| 成人国产一区最新在线观看| 国产精华一区二区三区| 亚洲熟女毛片儿| 欧美成人性av电影在线观看| 黑人操中国人逼视频| 长腿黑丝高跟| 国产99久久九九免费精品| 国产欧美日韩一区二区三| 亚洲欧美精品综合久久99| 国产精品久久久久久亚洲av鲁大| 国产国语露脸激情在线看| 亚洲 欧美 日韩 在线 免费| 又大又爽又粗| 久久香蕉精品热| 亚洲精华国产精华精| 成人国产综合亚洲| 99国产精品一区二区蜜桃av| 天天一区二区日本电影三级 | 色婷婷久久久亚洲欧美| 欧美日韩乱码在线| 精品久久久精品久久久| or卡值多少钱| 亚洲中文日韩欧美视频| 国产一区二区激情短视频| 国产精品一区二区免费欧美| 欧美激情 高清一区二区三区| 午夜久久久在线观看| 亚洲国产精品999在线| 老司机午夜十八禁免费视频| 韩国精品一区二区三区| 99久久国产精品久久久| 亚洲欧美精品综合久久99| 此物有八面人人有两片| 日韩有码中文字幕| 精品福利观看| av网站免费在线观看视频| 欧美 亚洲 国产 日韩一| 久久久精品欧美日韩精品| 女警被强在线播放| 国内久久婷婷六月综合欲色啪| av天堂久久9| 一级毛片女人18水好多| 身体一侧抽搐| 麻豆国产av国片精品| 亚洲三区欧美一区| 大香蕉久久成人网| 成人特级黄色片久久久久久久| 午夜福利在线观看吧| 757午夜福利合集在线观看| 日韩欧美免费精品| 国产熟女xx| av有码第一页| 99国产精品99久久久久| 亚洲狠狠婷婷综合久久图片| 91在线观看av| 丝袜美腿诱惑在线| 欧美国产精品va在线观看不卡| 国内毛片毛片毛片毛片毛片| 亚洲专区字幕在线| 大型av网站在线播放| 久久精品亚洲熟妇少妇任你| 狠狠狠狠99中文字幕| 亚洲自偷自拍图片 自拍| 久久中文字幕一级| 少妇熟女aⅴ在线视频| 搡老妇女老女人老熟妇| 亚洲最大成人中文| 51午夜福利影视在线观看| 国产亚洲精品一区二区www| 免费高清在线观看日韩| 亚洲国产日韩欧美精品在线观看 | 波多野结衣一区麻豆| 精品国产亚洲在线| 两性夫妻黄色片| 亚洲精品在线美女| 欧美成狂野欧美在线观看| 黄色 视频免费看| 男人操女人黄网站| 亚洲成人久久性| av福利片在线| 亚洲全国av大片| 国产极品粉嫩免费观看在线| 成人国产一区最新在线观看| 免费在线观看亚洲国产| 午夜影院日韩av| 成人国语在线视频| 一区福利在线观看| 精品午夜福利视频在线观看一区| 成人国产一区最新在线观看| 亚洲免费av在线视频| 国产亚洲精品第一综合不卡| 国产欧美日韩综合在线一区二区| ponron亚洲| 一进一出抽搐动态| 亚洲少妇的诱惑av| 成年女人毛片免费观看观看9| 美女大奶头视频| 欧美一级毛片孕妇| 一区二区日韩欧美中文字幕| 欧美国产日韩亚洲一区| 欧美黑人欧美精品刺激| 国产精品野战在线观看| 国产精品自产拍在线观看55亚洲| 国产免费男女视频| 日本a在线网址| 女性被躁到高潮视频| a级毛片在线看网站| 变态另类成人亚洲欧美熟女 | 国产精品乱码一区二三区的特点 | 亚洲欧美日韩另类电影网站| 一区二区三区精品91| 不卡一级毛片| 免费在线观看黄色视频的| 自线自在国产av| 欧美色欧美亚洲另类二区 | 国产在线精品亚洲第一网站| 久久久久久国产a免费观看| 丝袜美足系列| 日本免费a在线| 久久久久国产精品人妻aⅴ院| 啪啪无遮挡十八禁网站| 亚洲熟女毛片儿| 人人妻人人澡人人看| 美女大奶头视频| 亚洲精品国产色婷婷电影| 天堂√8在线中文| av电影中文网址| 一级,二级,三级黄色视频| 国产一级毛片七仙女欲春2 | 欧美乱码精品一区二区三区| 亚洲精华国产精华精| 国产成人欧美| 女人精品久久久久毛片| 亚洲视频免费观看视频| 黄色毛片三级朝国网站| 久久人妻熟女aⅴ| 精品少妇一区二区三区视频日本电影| 久久香蕉激情| 天天躁狠狠躁夜夜躁狠狠躁| 国产片内射在线| 在线天堂中文资源库| 搡老熟女国产l中国老女人| 69av精品久久久久久| 亚洲电影在线观看av| 亚洲人成伊人成综合网2020| 国产精品香港三级国产av潘金莲| 一级毛片女人18水好多| 制服丝袜大香蕉在线| 国产精品免费视频内射| 人人妻人人澡人人看| 久久狼人影院| 一级毛片女人18水好多| 51午夜福利影视在线观看| 亚洲欧美激情在线| 国产国语露脸激情在线看| 亚洲国产看品久久| 国产精品电影一区二区三区| 性色av乱码一区二区三区2| 一边摸一边做爽爽视频免费| 很黄的视频免费| 男女下面进入的视频免费午夜 | 国产熟女xx| 欧美成人一区二区免费高清观看 | 亚洲免费av在线视频| 一进一出好大好爽视频| 精品久久久久久久人妻蜜臀av | 国产成人免费无遮挡视频| 黄频高清免费视频| 夜夜看夜夜爽夜夜摸| 精品一区二区三区四区五区乱码| 在线观看免费视频网站a站| 欧美一级a爱片免费观看看 | 亚洲人成77777在线视频| 色老头精品视频在线观看| 真人一进一出gif抽搐免费| 波多野结衣av一区二区av| 性欧美人与动物交配| 国产真人三级小视频在线观看| 多毛熟女@视频| 岛国视频午夜一区免费看| 一级毛片高清免费大全| 久久精品国产亚洲av高清一级| 不卡一级毛片| 精品午夜福利视频在线观看一区| 亚洲精品美女久久av网站| 欧美精品啪啪一区二区三区| 如日韩欧美国产精品一区二区三区| 久久国产乱子伦精品免费另类| 亚洲欧洲精品一区二区精品久久久| 久久久久久国产a免费观看| 老司机在亚洲福利影院| 99精品久久久久人妻精品| 午夜影院日韩av| xxx96com| 欧美激情高清一区二区三区| 黄网站色视频无遮挡免费观看| 国产区一区二久久| 欧美老熟妇乱子伦牲交| 91在线观看av| 国产精华一区二区三区| 丁香欧美五月| 人妻丰满熟妇av一区二区三区| 欧美日韩乱码在线| 黄色 视频免费看| 一卡2卡三卡四卡精品乱码亚洲| 亚洲av成人av| 少妇熟女aⅴ在线视频| 久久香蕉激情| 黑人巨大精品欧美一区二区蜜桃| 久久伊人香网站| 免费高清在线观看日韩| 成人特级黄色片久久久久久久| 欧美中文日本在线观看视频| 99国产精品免费福利视频| 在线十欧美十亚洲十日本专区| 久久性视频一级片| 久久国产精品男人的天堂亚洲| 黄色a级毛片大全视频| 99国产综合亚洲精品| 国产99白浆流出| 欧美一级毛片孕妇| 亚洲中文av在线| 欧美 亚洲 国产 日韩一| 欧美人与性动交α欧美精品济南到| 美女免费视频网站| 欧美最黄视频在线播放免费| 亚洲黑人精品在线| 伊人久久大香线蕉亚洲五| 欧美中文综合在线视频| 久久精品人人爽人人爽视色| 手机成人av网站| bbb黄色大片| 日韩欧美一区二区三区在线观看| 国产亚洲精品综合一区在线观看 | 午夜亚洲福利在线播放| 成人av一区二区三区在线看| 成人18禁在线播放| 日韩视频一区二区在线观看| 夜夜看夜夜爽夜夜摸| 99久久精品国产亚洲精品| 国产精品,欧美在线| 国产成+人综合+亚洲专区| 成年人黄色毛片网站| 美女午夜性视频免费| 免费一级毛片在线播放高清视频 | 午夜精品国产一区二区电影| √禁漫天堂资源中文www| 极品教师在线免费播放| 久久久久久免费高清国产稀缺| 中出人妻视频一区二区| 黄色视频不卡| 一边摸一边做爽爽视频免费| 99久久综合精品五月天人人| 19禁男女啪啪无遮挡网站| 国产精品98久久久久久宅男小说| 国产区一区二久久| 在线天堂中文资源库| 精品不卡国产一区二区三区| 在线av久久热| 欧美日本中文国产一区发布| 久热爱精品视频在线9| 夜夜夜夜夜久久久久| 中文字幕人妻丝袜一区二区| 欧美激情久久久久久爽电影 | 国产免费av片在线观看野外av| 韩国av一区二区三区四区| 国产熟女午夜一区二区三区|