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

    耗氧型惰化系統(tǒng)反應器性能理論

    2019-12-02 10:49:22謝輝輝馮詩愚彭孝天潘俊王洋洋
    北京航空航天大學學報 2019年11期
    關鍵詞:催化劑

    謝輝輝, 馮詩愚,*, 彭孝天, 潘俊, 王洋洋

    (1. 南京航空航天大學航空學院飛行器環(huán)境控制與生命保障工業(yè)和信息化部重點實驗室, 南京 210016;2. 中國航空工業(yè)集團有限公司南京機電液壓工程研究中心航空機電系統(tǒng)綜合航空科技重點實驗室, 南京 211106)

    燃油箱的防火抑爆能力不僅關系到飛機生存能力和易損性,同時也關系到飛機利用率、成本和乘客生命安全[1]。因此,無論對于軍機還是民機,都必須采用有效措施來防止油箱燃爆事故發(fā)生。從20世紀70年代后期,通過機載設備制取惰化氣體的機載惰化技術發(fā)展十分迅速,相對于攜帶惰化氣體方式而言,后勤保障要求低,且能實現(xiàn)飛行過程中始終保持低于支持燃油燃燒的極限氧濃度,但該系統(tǒng)需要大量引氣,膜分離組件成本昂貴、壽命低。

    作為世界航空技術的最前沿陣地,美國的相關科研機構和廠商一直致力于研究新型的燃油箱防爆技術,以克服中空纖維膜應用的不足,例如冷卻惰化技術[2]、吸附惰化[3]、耗氧型惰化[4],而耗氧型惰化由于其流程簡單、效率高,被認為是下一代最可能得到實際應用的新型惰化系統(tǒng),其代表性的技術有派克宇航公司的綠色催化惰化技術和霍尼韋爾公司的蒸發(fā)催化惰化技術[5-8]。2013年,美軍成功地在UH-60(黑鷹)上進行了綠色惰化技術的試飛驗證,也在美國聯(lián)邦航空管理局(FAA)大西城實驗室的環(huán)境模擬倉進行了全飛行包線模擬[9-10],并計劃于2020年后向市場推廣。

    為了突破國外技術壁壘,南京機電液壓工程研究中心和南京航空航天大學從2014年開始,率先在國內(nèi)對新型惰化技術展開了研究,在國外綠色惰化技術的基礎上,提出了一種3CIS惰化技術(low temperature Controllable oxygen Consumed Catalytic Inerting System),其原理是將油箱上部燃油蒸汽混合物預熱后通入反應器,在催化劑的作用下,燃油蒸汽被氧化成水和二氧化碳,氧氣被消耗,反應后的低氧含量氣體再次被送入油箱上部,從而降低燃油箱氣體中的氧含量,達到惰化的目的。

    反應器作為3CIS系統(tǒng)核心部件,對惰化效率有很大影響。催化燃燒的反應器常見有固定床、蜂窩陶瓷載體、金屬基載體幾種[11],后兩種的催化劑為一塊整體,在載體中有許多微小通道,催化劑活性成分負載在通道上,接觸面積大、催化效率高、壓降低[12],但其設計難度大。固定床式反應器中填充有催化劑顆粒,結構簡單、可靠性強,因此也被廣泛使用,但是由于導熱性差的非金屬氧化物構成的催化劑載體之間為點接觸,所以床層間的導熱性差,可能會產(chǎn)生飛溫現(xiàn)象。

    與化學工業(yè)中反應器不同的是,后者通常工作在穩(wěn)定工況下,高效率是其追求的主要指標,而產(chǎn)物的溫度、設備的重量體積等則為次要指標。但是3CIS系統(tǒng)中,由于氧氣在惰化過程中不斷被消耗,燃油蒸汽濃度又與油箱溫度和載油量有關,所以進口反應物組成在整個惰化過程中變化極大,為了消除對燃油箱的潛在安全隱患,因此保證反應器出口氣體溫度在燃油自燃溫度以下是設計的主要指標,而重量和尺寸的要求也十分嚴格,但是對效率的追求則是次要指標。

    由于反應器中存在催化反應、氣體流動和傳熱等復雜因素的耦合,且3CIS系統(tǒng)中反應器尺寸較小,如果通過實驗測量內(nèi)部參數(shù)變化和溫度分布難度和代價較大。隨著計算機的發(fā)展,CFD作為一種有效科學的手段早已經(jīng)廣泛用于反應器的仿真。Benedettoa等[13]建立了二維的甲烷催化燃燒CFD模型,并耦合了流體流動與化學反應過程來討論不同壓力對反應速率、甲烷轉(zhuǎn)化率以及反應器內(nèi)溫度的影響。Smith等[14]利用Fluent建立了低溫水氣變換反應固定床模型,對4種不同的化學反應動力學進行模擬,但是這些模擬中認為催化劑固相和氣體始終處于熱平衡狀態(tài),這與實際過程存在一定差異。Gao等[15]對固定床反應器中流動行為及一氧化碳生成草酸二乙酯(DEO)的催化反應進行模擬研究,并采用兩溫度多孔介質(zhì)模型,結果表明預測結果與實驗重合度較熱平衡模型高。

    本文首先通過實驗測量了催化劑反應速率隨溫度和物料濃度的變化關系,基于阿倫尼烏斯關系式回歸出相關系數(shù),然后通過用戶自定義變量(UDS)添加固相能量方程,用戶自定義函數(shù)(UDF)添加催化反應,采用多孔介質(zhì)兩溫度模型對3CIS系統(tǒng)中的反應器在不同燃油蒸汽濃度、不同流速下的內(nèi)部溫度分布、出口產(chǎn)物組成、反應器效率等進行了模擬分析。

    1 催化劑性能測試實驗

    催化劑的性能需要通過實驗的方式進行測定,圖1給出了實驗流程圖,采用微量注射泵將燃油注入混合器中,通過調(diào)整3種氣體的流量,來調(diào)節(jié)反應氣體的濃度,并在混合器中與燃油混合,在加熱器中將其溫度調(diào)整至所需的溫度后進入石英反應器,反應后的產(chǎn)物通過離子火焰檢測器和氣相色譜來檢定其成分。

    在石英反應器內(nèi),催化劑的裝填量為5 mL,反應氣體流量為5 000 mL/min,空速為60 000 h-1,壓力為常壓,RP-3濃度保持0.48%,改變氧氣濃度為2%、4.5%、7.5%、10%、12.5%和15%,分別測試反應溫度為170、180、190和200℃時的反應速率。由于催化劑用量非常少,反應器直徑很小,可認為不存在縱向和橫向的溫度梯度和反應物的濃度梯度,測試結果如圖2所示,當氧氣濃度低于10%時,反應速度與氧氣濃度呈近似線性關系,當氧濃度高于10%時,反應速率與氧氣濃度幾乎無關。

    控制氧氣濃度為10%,調(diào)整RP-3濃度為0.13%、0.26%、0.48%、0.7%、0.9%和1.5%,測試結構如圖3所示,從圖中可見,在實驗溫度范圍內(nèi)反應溫度為170、 180、190和200℃時的反應速率隨燃油濃度增加而線性上升,屬一級依賴關系。

    采用阿倫尼烏尼斯公式,聯(lián)立圖2和圖3數(shù)據(jù)求出活化能、指數(shù)前因子,得到反應動力學方程為

    (1)

    圖1 催化劑性能測試實驗Fig.1 Catalyst performance test experiment

    圖2 不同氧氣濃度下反應速率隨溫度變化Fig. 2 Variation of reaction rate with temperature at different oxygen concentrations

    圖3 不同溫度下反應速率隨RP-3濃度變化Fig.3 Variation of reaction rate with RP-3 concentration at different temperatures

    式中:CRP-3為RP-3的物質(zhì)的量濃度;CO2為氧氣的物質(zhì)的量濃度;XO2為氧氣體積分數(shù);Rr為RP-3的反應速率;R為理想氣體常數(shù);Tf為氣相溫度。

    2 物理模型及控制方程

    2.1 物理模型

    反應器結構和尺寸如圖4所示,內(nèi)部空間分為3個部分,其中中間為顆粒狀催化劑隨機堆積而成的催化床層。

    從燃油箱中抽出的燃油蒸汽混合物被加熱后進入反應器在催化床層中發(fā)生化學反應生成水蒸氣和二氧化碳,并釋放熱量。RP-3燃油其氫質(zhì)量分數(shù)為14.46%,分子式可簡化為C7.12H14.16,分子量為141 g/mol,其化學反應方程為

    圖4 催化反應器結構Fig.4 Catalytic reactor structure

    C7.12H14.16+10.74O2=7.12CO2+7.23H2O

    (2)

    燃燒熱值為43.53 MJ/kg。反應器內(nèi)部主要參數(shù)由表1列出。

    表1 反應器催化床主要參數(shù)Table 1 Main parameters of reactor catalytic bed

    2.2 控制方程和邊界條件

    本研究中,對反應器模型提出以下幾點假設。

    1) 由于催化劑顆粒的直徑很小,將反應器催化床孔隙率視為常值,并且催化劑內(nèi)化學性質(zhì)均一,將模型簡化為擬均相模型。

    2) 催化床為多孔介質(zhì),添加動力阻力項。

    3) 反應氣體為不可壓縮氣體,適用理想氣體方程。

    4) 由于化學反應發(fā)生在催化劑顆粒表面以及內(nèi)部孔隙中,認為化學反應熱全部由固相吸收,通過氣固相傳熱達到平衡。

    5) 采用氣-固相局部熱平衡,在多孔介質(zhì)中添加固體能量方程和熱傳遞方程。

    6) 反應物的體積力對流場分布影響小,因此可忽略不計。

    在仿真中采用Fluent 17.0多孔介質(zhì)模型來簡化反應器床層復雜無規(guī)律的空隙通道,實際上是給氣體添加額外的動量源項來表征催化劑顆粒給氣體帶來的阻礙作用。采用UDS在多孔介質(zhì)區(qū)域添加固相能量方程,化學反應通過UDF給對應的組分添加質(zhì)量源項表征化學反應中各組分濃度的變化,化學反應熱以能量源項形式添加給固相能量方程,并通過傳熱作用將熱量傳遞給氣相。以下為多孔介質(zhì)區(qū)域主要控制方程。

    氣相質(zhì)量守恒:

    (3)

    式中:γ為催化床孔隙率;ρg為氣相密度;v為氣體速度。

    氣相動量方程為

    (4)

    式中:τ為剪切應力;P為氣相壓力;S為動量源項,表達式為

    (5)

    式中:μg為氣體動力黏度;dp為催化劑顆粒直徑。

    氣相狀態(tài)方程為

    (6)

    式中:Yi為氣相第i組分質(zhì)量分數(shù);Mi為氣相第i組分摩爾質(zhì)量。

    組分運輸方程為

    (7)

    式中:Ri為第i組分的質(zhì)量源項;Nr為組分的序號;Ji為氣相第i組分的質(zhì)量擴散:

    (8)

    其中:Di,m為氣相第i組分質(zhì)量擴散系數(shù);μt為湍流黏度;Sct為湍流施密特數(shù)。

    氣相能量方程為

    (9)

    式中:Ef為氣體能量;kf為氣相導熱系數(shù);Ts為固相溫度;Cp.g為氣相比熱容;hi為氣相第i組分焓變;hv為氣固相對流傳熱系數(shù)。

    固相能量方程為

    hv(Ts-Tf)+ΦR

    (10)

    式中:ρs為固相密度;ks為固相導熱系數(shù);ΦR為化學反應放熱;Cp·s為固定相比熱容。

    氣固相對流傳熱系數(shù)為

    (11)

    式中:Db為反應器管徑;Re為雷諾數(shù);Nur為努塞爾數(shù);Pr為普朗特數(shù);kg為氣相導熱系數(shù)。

    Guardo等[16]在Fluent 6.0平臺上使用5種不同的湍流模型:Spalart-Allmaras 模型,k-ω模型, 標準k-ε模型, Realizablek-ε模型和RNGk-ε模型來研究湍流模型對固定床壁面?zhèn)鳠岬挠绊?,研究表明使用Spalart-Allmaras 模型更真實反映氣體在復雜的床層間流動情況。湍流方程為

    (12)

    在入口處使用質(zhì)量流量入口邊界條件,規(guī)定了入口氣體的質(zhì)量流量、溫度、各組分摩爾分數(shù)、壓力,出口采用壓力出口邊界條件。

    壁面條件采用對流熱傳導邊界條件,認為惰化系統(tǒng)中環(huán)境溫度為40℃,反應器外壁面對流傳熱系數(shù)為50 W·m-2·s-1。

    3 仿真結果與分析

    3.1 網(wǎng)格無關性驗證

    反應器采用軸對稱模型,將三維模型簡化為兩維。為了證明網(wǎng)格的數(shù)量滿足仿真要求,分別采用5套不同網(wǎng)格來計算氣體通過反應器產(chǎn)生的壓降,結果如圖5所示。75×150網(wǎng)格數(shù)量已滿足要求,因此采用此數(shù)量網(wǎng)格進行仿真。

    圖5 網(wǎng)格數(shù)量與反應器壓降關系Fig.5 Relationship between number of grids and reactor pressure drop

    3.2 反應器靜態(tài)特性

    改變反應器進口氧氣濃度與RP-3濃度,固定氣體流速為1 m/s,反應溫度為210℃,得到不同操作條件下的RP-3轉(zhuǎn)化率,如圖6所示。由圖6可知,在相同RP-3濃度下,增大氧氣濃度,轉(zhuǎn)化率也隨之增加,當氧氣濃度大于某個值時,轉(zhuǎn)化率不再隨著氧氣濃度增加而增加(將這個值命名為氧濃度飽和值)。在RP-3濃度為0.13%~0.48%條件下,氧濃度飽和值大約在12%左右,而隨著RP-3濃度增加,氧濃度飽和值會增加。在小于不同RP-3濃度對應飽和值的氧氣濃度區(qū)間內(nèi),如圖中氧氣濃度為6%~10%時,對于相同氧濃度值,增加RP-3濃度會降低轉(zhuǎn)化率;而在大于RP-3濃度對應飽和值的氧氣濃度區(qū)間內(nèi),有相反的特性。

    固定RP-3濃度為0.26%,氧氣濃度為15%,改變溫度與氣體流速得到圖7數(shù)據(jù)??芍诓煌瑲怏w流速下,隨著溫度增加轉(zhuǎn)化率增加率不同。流速較小時溫度對轉(zhuǎn)化率的促進作用要強于流速較大時;RP-3轉(zhuǎn)化率隨反應溫度增加而增加,當反應溫度在240℃以上時,RP-3轉(zhuǎn)化率增加緩慢甚至幾乎不變,由于催化劑設計時希望最佳工作溫度不超過200℃左右且反應動力學方程主要由170~200℃工況下的實驗數(shù)據(jù)擬合,對于阿倫尼烏斯公式擬合的反應動力學方程其仿真結果在溫度較高時誤差有所增大,因此其表現(xiàn)為反應當溫度遠離200℃時,溫度對化學反應的促進作用衰減。

    圖6 RP-3轉(zhuǎn)化率與反應物濃度關系Fig.6 Relationship between RP-3 conversion rate and reactant concentration

    3.3 反應器內(nèi)溫度及各組分濃度動態(tài)特性

    假設進入反應器的燃油蒸汽混合物被預熱到150℃,賦予反應器初始溫度150℃,RP-3濃度為0.26%,氧氣濃度為15%,氣體流速為1 m/s,研究反應器被充分預熱后其內(nèi)部氣-固相溫度及RP-3濃度變化過程。圖8和圖9分別給出了反應器在80 s時刻的氣體溫度云圖以及RP-3濃度云圖。圖10出了t=40,80,120 s時刻對應的催化床軸向氣固溫度。圖11給出了3個時刻的RP-3濃度沿反應器軸向變化曲線。圖12、圖13分別給出了0~120 s時間內(nèi)催化床3個截面的氣相平均溫度及RP-3濃度變化曲線。

    由圖10和圖11可以看出,由于反應熱被固相吸收,固相溫度大于氣相溫度。40 s時刻由于催化床溫度較低,反應較為緩慢,RP-3轉(zhuǎn)化率低;到了80 s時刻反應加劇,溫度沿著床層迅速增加,軸線上的RP-3被完全消耗,熱點位于催化床后端;120 s時刻,由于反應器催化床的導熱性差,催化床后端溫度持續(xù)上升,催化床前后溫差巨大。由圖12得出反應器出口溫度在80 s后增加迅速,出現(xiàn)飛溫現(xiàn)象。由圖13可以看出,反應器在X=0.075m處轉(zhuǎn)化率一直處于較低水平,即反應器只有后半段工作效率高。

    圖8 t=80 s時刻反應器內(nèi)氣體溫度云圖Fig.8 Gas temperature contour in reactor at t=80 s

    圖9 t=80 s時刻反應器內(nèi)RP-3濃度云圖Fig 9 RP-3 concentration contour in reactor at t=80 s

    圖10 不同時刻催化床軸線處溫度分布Fig.10 Temperature distribution at axis of catalytic bed at different moments

    圖11 不同時刻反應器軸線處RP-3濃度分布Fig.11 RP-3 concentration distribution at reactor axis at different moments

    圖12 不同截面處氣相溫度隨時間變化Fig.12 Variation of gas phase temperature at different cross-sections with time

    圖13 不同截面處RP-3平均濃度隨時間變化Fig.13 Variation of average concentration of RP-3 at different cross-sections with time

    圖14 t=80 s時刻催化床不同軸向處RP-濃度Fig.14 RP-3 concentration in different axial positions of catalytic bed at t=80 s

    圖15 t=80 s時刻催化床不同軸向處徑向溫度Fig.15 Radial temperature in different axial positions of catalytic bed at t=80 s

    圖14和圖15給出了t=80 s時刻不同軸向處的RP-3濃度和徑向溫度分布。可以發(fā)現(xiàn)靠近反應器外壁處,由于邊界的黏性作用以及外壁自然對流對熱量的損失,加之催化床的導熱性差,溫度難以向壁面?zhèn)鬟f,造成中心處的溫度遠高于靠近外壁處;不均勻的溫度又造成RP-3濃度在徑向上有巨大差異。由圖15可知,化學反應主要發(fā)生在靠近反應器軸線的位置。

    4 結 論

    通過UDS添加固相能量方程,以UDF源項形式添加化學反應的方式,以多孔介質(zhì)模型來代替復雜的顆粒催化床,建立了氣-固耦合傳熱的反應器,對耗氧惰化系統(tǒng)中的反應器進行仿真,得出以下結論。

    1) 當達到氧濃度飽和值,增加氧氣濃度將不對RP-3轉(zhuǎn)化起促進作用;在小于氧濃度飽和值的氧氣濃度區(qū)間內(nèi),RP-3濃度越大,轉(zhuǎn)化率越低,在大于氧濃度飽和值的氧氣濃度區(qū)間內(nèi),RP-3濃度越小,轉(zhuǎn)化率越低。

    2) 反應過程中,熱點位于靠近催化床出口位置,隨著時間推移,RP-3轉(zhuǎn)化率迅速增加,因此不采取額外的冷卻措施,反應器將飛溫。

    3) 化學反應主要發(fā)生在反應器的后半段,且靠近反應器軸線處,使催化床的溫度均勻分布對提高反應器效率有重要意義。

    本文將反應器作為獨立的對象,未充分與惰化系統(tǒng)聯(lián)系且反應器有飛溫現(xiàn)象,因此后續(xù)工作將重點解決飛溫現(xiàn)象,并與系統(tǒng)聯(lián)系起來展開進一步研究。

    猜你喜歡
    催化劑
    走近諾貝爾獎 第三種催化劑
    大自然探索(2023年7期)2023-11-14 13:08:06
    直接轉(zhuǎn)化CO2和H2為甲醇的新催化劑
    鋁鎳加氫催化劑在BDO裝置運行周期的探討
    如何在開停產(chǎn)期間保護克勞斯催化劑
    智富時代(2018年3期)2018-06-11 16:10:44
    新型釩基催化劑催化降解氣相二噁英
    掌握情欲催化劑
    Coco薇(2016年2期)2016-03-22 02:45:06
    碳包覆鐵納米催化劑的制備及其表征
    V2O5-WO3/TiO2脫硝催化劑回收研究進展
    負載型催化劑(CuO/TUD-1,CuO/MCM-41)的制備及其在一步法氧化苯合成苯酚中的應用
    合成化學(2015年4期)2016-01-17 09:01:27
    復合固體超強酸/Fe2O3/AI2O3/ZnO/ZrO2催化劑的制備及其催化性能
    應用化工(2014年3期)2014-08-16 13:23:50
    久久久国产精品麻豆| 久久青草综合色| 中文字幕人妻熟女乱码| 国产午夜精品久久久久久| 国产高清视频在线播放一区| 窝窝影院91人妻| 欧美成狂野欧美在线观看| 欧美黑人欧美精品刺激| 精品亚洲乱码少妇综合久久| 午夜福利视频精品| 婷婷丁香在线五月| 一级片免费观看大全| 精品亚洲乱码少妇综合久久| 精品卡一卡二卡四卡免费| 国产精品亚洲av一区麻豆| 国产精品影院久久| 久久午夜综合久久蜜桃| 69av精品久久久久久 | 色94色欧美一区二区| 丰满饥渴人妻一区二区三| 精品一区二区三区av网在线观看 | 热99re8久久精品国产| 欧美老熟妇乱子伦牲交| 另类精品久久| av国产精品久久久久影院| 三级毛片av免费| 肉色欧美久久久久久久蜜桃| 黑人巨大精品欧美一区二区蜜桃| 国产老妇伦熟女老妇高清| 国产男女内射视频| a在线观看视频网站| 欧美黑人欧美精品刺激| 精品亚洲成国产av| 日韩中文字幕欧美一区二区| 国产成人精品久久二区二区免费| 婷婷丁香在线五月| 美女主播在线视频| 免费人妻精品一区二区三区视频| 久久久水蜜桃国产精品网| 国产xxxxx性猛交| 999精品在线视频| 色综合婷婷激情| 久久国产精品人妻蜜桃| 婷婷丁香在线五月| 精品一区二区三区四区五区乱码| 中亚洲国语对白在线视频| 肉色欧美久久久久久久蜜桃| 天堂动漫精品| 巨乳人妻的诱惑在线观看| 午夜视频精品福利| 欧美日韩av久久| 国产在视频线精品| 久久久久精品国产欧美久久久| 久久香蕉激情| 欧美精品一区二区大全| 久久久久精品人妻al黑| 亚洲欧美色中文字幕在线| 亚洲色图综合在线观看| 99国产精品一区二区蜜桃av | 精品国产一区二区久久| 久久精品国产亚洲av香蕉五月 | 最近最新免费中文字幕在线| 婷婷成人精品国产| 十分钟在线观看高清视频www| 国产主播在线观看一区二区| 黄色怎么调成土黄色| 中文字幕另类日韩欧美亚洲嫩草| 变态另类成人亚洲欧美熟女 | 啦啦啦中文免费视频观看日本| 国产精品成人在线| 人成视频在线观看免费观看| 亚洲中文av在线| 99国产精品免费福利视频| 国产精品久久久av美女十八| 五月天丁香电影| 国产在视频线精品| 青青草视频在线视频观看| 国产高清videossex| 夜夜夜夜夜久久久久| 在线观看免费午夜福利视频| 国产成人精品无人区| 性高湖久久久久久久久免费观看| 国产淫语在线视频| 亚洲精品中文字幕在线视频| 丝袜美足系列| 一级毛片女人18水好多| www日本在线高清视频| 悠悠久久av| 人妻 亚洲 视频| 国产不卡一卡二| 在线观看免费视频网站a站| 日韩欧美一区二区三区在线观看 | 色综合婷婷激情| 女性被躁到高潮视频| 大片免费播放器 马上看| 一本大道久久a久久精品| 757午夜福利合集在线观看| 亚洲色图av天堂| 国产精品1区2区在线观看. | 黑人猛操日本美女一级片| 久久久久久久久免费视频了| 欧美精品av麻豆av| 女警被强在线播放| 交换朋友夫妻互换小说| cao死你这个sao货| 欧美日韩亚洲综合一区二区三区_| 国产午夜精品久久久久久| 午夜福利一区二区在线看| 精品亚洲成a人片在线观看| 少妇粗大呻吟视频| 成人特级黄色片久久久久久久 | 国产真人三级小视频在线观看| 一级黄色大片毛片| av天堂久久9| 国产精品偷伦视频观看了| 老熟女久久久| 我的亚洲天堂| 久久九九热精品免费| 性高湖久久久久久久久免费观看| 免费黄频网站在线观看国产| 午夜福利视频精品| 免费久久久久久久精品成人欧美视频| cao死你这个sao货| 亚洲成人国产一区在线观看| 国产成人系列免费观看| 精品一区二区三区视频在线观看免费 | 免费看a级黄色片| 午夜视频精品福利| 国产成人欧美| 国产野战对白在线观看| 欧美变态另类bdsm刘玥| 国产精品久久久人人做人人爽| 高清视频免费观看一区二区| 久久国产精品大桥未久av| 日本a在线网址| 亚洲av日韩在线播放| 亚洲五月婷婷丁香| 亚洲人成电影免费在线| 无人区码免费观看不卡 | 欧美 亚洲 国产 日韩一| 手机成人av网站| 老汉色∧v一级毛片| 精品亚洲乱码少妇综合久久| av在线播放免费不卡| 18禁国产床啪视频网站| 999久久久国产精品视频| 在线观看免费视频日本深夜| 亚洲中文字幕日韩| 999精品在线视频| 精品亚洲乱码少妇综合久久| 久久狼人影院| av在线播放免费不卡| 99国产精品免费福利视频| 国产精品一区二区在线观看99| 蜜桃在线观看..| 少妇的丰满在线观看| 国产国语露脸激情在线看| 两人在一起打扑克的视频| 两个人看的免费小视频| 丁香六月欧美| 国产老妇伦熟女老妇高清| 丝瓜视频免费看黄片| 国产主播在线观看一区二区| 亚洲,欧美精品.| 久久国产精品大桥未久av| 电影成人av| 美女主播在线视频| av网站在线播放免费| 久久久久精品国产欧美久久久| 男男h啪啪无遮挡| 777久久人妻少妇嫩草av网站| 一区二区三区激情视频| 亚洲欧美激情在线| 无遮挡黄片免费观看| 国产不卡av网站在线观看| 午夜福利,免费看| 日韩欧美三级三区| 日韩欧美一区二区三区在线观看 | 9191精品国产免费久久| 亚洲精品国产色婷婷电影| 91老司机精品| 美国免费a级毛片| 一二三四在线观看免费中文在| 另类精品久久| 色在线成人网| 亚洲视频免费观看视频| 亚洲人成77777在线视频| 超色免费av| 成人国语在线视频| 搡老乐熟女国产| 日韩免费高清中文字幕av| 亚洲欧美日韩另类电影网站| 亚洲精品一卡2卡三卡4卡5卡| 亚洲五月色婷婷综合| 丝袜在线中文字幕| 国产不卡一卡二| 岛国毛片在线播放| 国产有黄有色有爽视频| 老鸭窝网址在线观看| 别揉我奶头~嗯~啊~动态视频| 中文字幕色久视频| 搡老乐熟女国产| 色老头精品视频在线观看| 不卡一级毛片| 国产野战对白在线观看| 18禁国产床啪视频网站| 久久久国产精品麻豆| 啪啪无遮挡十八禁网站| 亚洲国产av影院在线观看| 久久影院123| 美女福利国产在线| 亚洲精品成人av观看孕妇| 国产欧美日韩一区二区三区在线| 亚洲少妇的诱惑av| 久久毛片免费看一区二区三区| 成在线人永久免费视频| 一级,二级,三级黄色视频| 国产国语露脸激情在线看| 亚洲国产中文字幕在线视频| 怎么达到女性高潮| 一级a爱视频在线免费观看| 久热爱精品视频在线9| 精品国产国语对白av| 色综合欧美亚洲国产小说| 99热国产这里只有精品6| 精品久久久久久电影网| 亚洲 欧美一区二区三区| 777久久人妻少妇嫩草av网站| 日本五十路高清| 精品乱码久久久久久99久播| 色94色欧美一区二区| 久久精品成人免费网站| 91麻豆av在线| 亚洲精品av麻豆狂野| 不卡av一区二区三区| www.自偷自拍.com| 黄片小视频在线播放| 欧美激情极品国产一区二区三区| 欧美日韩亚洲综合一区二区三区_| 亚洲av日韩在线播放| 精品国产乱码久久久久久小说| 国产人伦9x9x在线观看| 咕卡用的链子| 每晚都被弄得嗷嗷叫到高潮| 欧美av亚洲av综合av国产av| 国产在线视频一区二区| 日韩人妻精品一区2区三区| 日本一区二区免费在线视频| av福利片在线| 久久人人爽av亚洲精品天堂| 又黄又粗又硬又大视频| 18禁观看日本| 美女高潮喷水抽搐中文字幕| 国产成人精品久久二区二区91| 久久九九热精品免费| 久久久久精品国产欧美久久久| 国产精品亚洲一级av第二区| 一本色道久久久久久精品综合| 丝袜美足系列| 麻豆乱淫一区二区| 久久国产精品男人的天堂亚洲| 日本a在线网址| 久久久久国产一级毛片高清牌| 成年人午夜在线观看视频| 国产xxxxx性猛交| av网站免费在线观看视频| 成人18禁在线播放| 亚洲欧美色中文字幕在线| 亚洲九九香蕉| 久久精品国产99精品国产亚洲性色 | 99re6热这里在线精品视频| 超色免费av| 我的亚洲天堂| 少妇粗大呻吟视频| 大香蕉久久网| bbb黄色大片| 国产成人免费无遮挡视频| 日韩制服丝袜自拍偷拍| 人人妻人人澡人人爽人人夜夜| 亚洲国产av新网站| 热re99久久国产66热| 国产亚洲精品第一综合不卡| 日本黄色日本黄色录像| 日韩人妻精品一区2区三区| www.精华液| 亚洲五月色婷婷综合| 亚洲成人国产一区在线观看| 又大又爽又粗| 最新美女视频免费是黄的| 在线观看www视频免费| 国产主播在线观看一区二区| a在线观看视频网站| 人人妻人人添人人爽欧美一区卜| 亚洲精品国产区一区二| 日本黄色日本黄色录像| 国产欧美日韩综合在线一区二区| 午夜福利视频精品| 亚洲专区中文字幕在线| 国产av精品麻豆| 亚洲综合色网址| 在线观看免费视频日本深夜| 亚洲久久久国产精品| 亚洲一卡2卡3卡4卡5卡精品中文| 国产1区2区3区精品| 女人高潮潮喷娇喘18禁视频| 精品国产一区二区三区久久久樱花| 成人国产av品久久久| 久久精品aⅴ一区二区三区四区| 色精品久久人妻99蜜桃| 国产人伦9x9x在线观看| 老汉色∧v一级毛片| 久久久国产精品麻豆| 免费不卡黄色视频| 国产老妇伦熟女老妇高清| 国产不卡av网站在线观看| 成人手机av| 91精品国产国语对白视频| tube8黄色片| 最近最新免费中文字幕在线| 国产精品欧美亚洲77777| 一区二区三区激情视频| 国产在线免费精品| 极品人妻少妇av视频| 两性午夜刺激爽爽歪歪视频在线观看 | 成年人午夜在线观看视频| 另类亚洲欧美激情| 99香蕉大伊视频| 久久99热这里只频精品6学生| 午夜视频精品福利| 亚洲九九香蕉| 夫妻午夜视频| 黄网站色视频无遮挡免费观看| 国产欧美亚洲国产| 国产成人免费观看mmmm| 搡老熟女国产l中国老女人| 欧美精品高潮呻吟av久久| 欧美精品人与动牲交sv欧美| 天堂8中文在线网| 嫩草影视91久久| 国产有黄有色有爽视频| 国产精品欧美亚洲77777| 超碰97精品在线观看| 精品少妇黑人巨大在线播放| 精品国产一区二区久久| 不卡av一区二区三区| 91麻豆精品激情在线观看国产 | 亚洲欧美色中文字幕在线| 中国美女看黄片| 日日爽夜夜爽网站| 动漫黄色视频在线观看| 视频区欧美日本亚洲| 精品午夜福利视频在线观看一区 | 日本五十路高清| 黑丝袜美女国产一区| 美女扒开内裤让男人捅视频| 久久国产亚洲av麻豆专区| 欧美av亚洲av综合av国产av| 久久天躁狠狠躁夜夜2o2o| 午夜免费鲁丝| 十八禁高潮呻吟视频| 亚洲精品中文字幕一二三四区 | avwww免费| a在线观看视频网站| 精品少妇内射三级| 国产又色又爽无遮挡免费看| 一边摸一边抽搐一进一小说 | 一区二区三区国产精品乱码| 精品国产一区二区三区四区第35| 中国美女看黄片| 亚洲色图综合在线观看| 无遮挡黄片免费观看| 丝袜人妻中文字幕| 夫妻午夜视频| 欧美日韩亚洲高清精品| 国产伦人伦偷精品视频| 国产在线精品亚洲第一网站| 蜜桃国产av成人99| 天天操日日干夜夜撸| 啦啦啦视频在线资源免费观看| 黄色 视频免费看| 精品人妻熟女毛片av久久网站| 桃花免费在线播放| 黑人巨大精品欧美一区二区mp4| 久久精品人人爽人人爽视色| 欧美中文综合在线视频| 超碰成人久久| 精品人妻熟女毛片av久久网站| 亚洲国产欧美在线一区| 精品国产亚洲在线| av有码第一页| 国产亚洲午夜精品一区二区久久| 成人av一区二区三区在线看| 中文字幕人妻丝袜制服| 久久久国产成人免费| 国产精品久久久久成人av| www.999成人在线观看| 成人三级做爰电影| 超碰97精品在线观看| 露出奶头的视频| 成人国产av品久久久| 国产欧美日韩一区二区三| 欧美成狂野欧美在线观看| 久久 成人 亚洲| 伊人久久大香线蕉亚洲五| 大陆偷拍与自拍| 亚洲国产中文字幕在线视频| 亚洲一区中文字幕在线| 在线永久观看黄色视频| 久久香蕉激情| 亚洲专区字幕在线| 国产91精品成人一区二区三区 | 亚洲精品中文字幕在线视频| 成人亚洲精品一区在线观看| 国产精品免费大片| 每晚都被弄得嗷嗷叫到高潮| 国产一区二区激情短视频| 热99re8久久精品国产| 欧美精品亚洲一区二区| 十八禁网站网址无遮挡| 女人精品久久久久毛片| 女性被躁到高潮视频| 国产欧美日韩精品亚洲av| 国产av又大| 黑人操中国人逼视频| 亚洲精品一二三| 99热网站在线观看| 水蜜桃什么品种好| 色94色欧美一区二区| svipshipincom国产片| 黄色成人免费大全| 国产成人精品在线电影| 色综合欧美亚洲国产小说| 不卡一级毛片| 国产成人精品在线电影| 成人影院久久| 丝瓜视频免费看黄片| 丝袜人妻中文字幕| 亚洲欧洲日产国产| 国产精品影院久久| 久久免费观看电影| 亚洲国产欧美网| 国产一卡二卡三卡精品| 国产成人免费观看mmmm| 婷婷成人精品国产| 99国产精品一区二区蜜桃av | 人人妻人人添人人爽欧美一区卜| 欧美在线黄色| 一区二区av电影网| 中文亚洲av片在线观看爽 | 99精国产麻豆久久婷婷| 青青草视频在线视频观看| 丁香欧美五月| 亚洲一卡2卡3卡4卡5卡精品中文| 成年版毛片免费区| 国产成人精品在线电影| 欧美国产精品一级二级三级| 国产亚洲精品第一综合不卡| 青草久久国产| a级毛片黄视频| 亚洲av电影在线进入| 精品第一国产精品| videos熟女内射| 嫁个100分男人电影在线观看| 无限看片的www在线观看| 国产一区二区激情短视频| 成人影院久久| 搡老岳熟女国产| 亚洲熟女精品中文字幕| 国产精品久久久久久精品电影小说| 久久中文看片网| 老司机福利观看| 精品国产亚洲在线| 国产精品电影一区二区三区 | 十八禁人妻一区二区| 国产精品电影一区二区三区 | 少妇精品久久久久久久| 日韩人妻精品一区2区三区| 亚洲 欧美一区二区三区| 91麻豆精品激情在线观看国产 | 国产视频一区二区在线看| 国产av精品麻豆| 午夜激情久久久久久久| 日韩欧美国产一区二区入口| 一区二区三区精品91| 精品国产国语对白av| 啦啦啦 在线观看视频| 久久久精品94久久精品| 丝瓜视频免费看黄片| 免费观看a级毛片全部| 国产日韩欧美视频二区| 亚洲国产精品一区二区三区在线| 色播在线永久视频| 日韩三级视频一区二区三区| 精品第一国产精品| 亚洲国产欧美日韩在线播放| 啦啦啦 在线观看视频| 少妇的丰满在线观看| 99九九在线精品视频| 极品少妇高潮喷水抽搐| 视频区欧美日本亚洲| 怎么达到女性高潮| 欧美午夜高清在线| 亚洲第一欧美日韩一区二区三区 | 波多野结衣av一区二区av| 欧美日韩福利视频一区二区| 在线观看66精品国产| 国产又色又爽无遮挡免费看| 后天国语完整版免费观看| avwww免费| 又大又爽又粗| 男男h啪啪无遮挡| 国产高清激情床上av| 亚洲人成伊人成综合网2020| 丰满饥渴人妻一区二区三| 日韩视频一区二区在线观看| 免费久久久久久久精品成人欧美视频| 天天躁夜夜躁狠狠躁躁| 十八禁人妻一区二区| 老司机深夜福利视频在线观看| 久久av网站| 成人国语在线视频| 999久久久国产精品视频| 欧美精品人与动牲交sv欧美| 啦啦啦视频在线资源免费观看| 久久天堂一区二区三区四区| 亚洲成a人片在线一区二区| 黑人巨大精品欧美一区二区mp4| 叶爱在线成人免费视频播放| 国产精品免费大片| 不卡一级毛片| 亚洲精品中文字幕一二三四区 | 新久久久久国产一级毛片| 最新的欧美精品一区二区| 老司机亚洲免费影院| 国产精品一区二区在线观看99| 丰满人妻熟妇乱又伦精品不卡| 丝袜喷水一区| 久久久久国内视频| 成人永久免费在线观看视频 | 免费观看av网站的网址| 久久热在线av| 人妻 亚洲 视频| 最新在线观看一区二区三区| 麻豆乱淫一区二区| 国产精品九九99| 青青草视频在线视频观看| 1024香蕉在线观看| 91大片在线观看| 男女无遮挡免费网站观看| 久久久久久亚洲精品国产蜜桃av| 午夜免费成人在线视频| 一进一出抽搐动态| 女人被躁到高潮嗷嗷叫费观| 最近最新中文字幕大全免费视频| 在线观看免费高清a一片| 一本一本久久a久久精品综合妖精| 中文字幕最新亚洲高清| 两个人免费观看高清视频| 精品午夜福利视频在线观看一区 | 高清视频免费观看一区二区| 亚洲欧洲精品一区二区精品久久久| 少妇被粗大的猛进出69影院| 国产亚洲欧美精品永久| 久久国产精品男人的天堂亚洲| 狂野欧美激情性xxxx| 999久久久精品免费观看国产| 免费人妻精品一区二区三区视频| 水蜜桃什么品种好| 欧美精品一区二区大全| 免费在线观看视频国产中文字幕亚洲| 亚洲av日韩精品久久久久久密| 久久 成人 亚洲| 性高湖久久久久久久久免费观看| 人人妻人人添人人爽欧美一区卜| 亚洲视频免费观看视频| 伦理电影免费视频| 久久天堂一区二区三区四区| 正在播放国产对白刺激| 久热这里只有精品99| 老司机福利观看| 久久狼人影院| 自拍欧美九色日韩亚洲蝌蚪91| 丁香六月天网| 亚洲av国产av综合av卡| 国产精品国产高清国产av | 黄色视频不卡| 中文字幕av电影在线播放| 十八禁网站网址无遮挡| 亚洲 欧美一区二区三区| www.999成人在线观看| 亚洲精品久久午夜乱码| 国产精品二区激情视频| 日本av手机在线免费观看| 捣出白浆h1v1| 久久99热这里只频精品6学生| 美女高潮到喷水免费观看| av电影中文网址| 日韩人妻精品一区2区三区| 妹子高潮喷水视频| 在线观看免费日韩欧美大片| 在线观看66精品国产| 一区在线观看完整版| 女人高潮潮喷娇喘18禁视频| 高清av免费在线| 国产欧美日韩一区二区三| 精品国内亚洲2022精品成人 | 日本欧美视频一区| 久久久久精品国产欧美久久久| 亚洲精品美女久久久久99蜜臀| 午夜日韩欧美国产| 99国产综合亚洲精品| 精品福利观看| 国产av又大| 新久久久久国产一级毛片| 熟女少妇亚洲综合色aaa.|