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

    好氧發(fā)酵過程微界面?zhèn)髻|(zhì)耦合模擬分析

    2020-11-18 01:57:24馮堯成任厲泰張鋒張志炳
    化工學(xué)報 2020年11期
    關(guān)鍵詞:傳質(zhì)氣液氣泡

    馮堯成,任厲泰,張鋒,張志炳

    (南京大學(xué)化學(xué)化工學(xué)院,江蘇南京210023)

    引 言

    2,3-丁二醇作為一種重要的化工原料和液體燃料,廣泛應(yīng)用于化工、能源及食品等領(lǐng)域[1-2],因此尋求經(jīng)濟(jì)、高效和低污染的2,3-丁二醇制備方式并最終將其產(chǎn)業(yè)化一直是科學(xué)研究的熱點[3]。相較于傳統(tǒng)化學(xué)法生產(chǎn)2,3-丁二醇,微生物發(fā)酵法因其工藝條件溫和、操作簡便、成本低廉的特點,優(yōu)勢明顯[4]。雖然目前國內(nèi)外的一系列研究都能達(dá)到相當(dāng)高的糖轉(zhuǎn)化率,但這些成果大多數(shù)都建立在消耗高葡萄糖含量的精細(xì)底物淀粉的基礎(chǔ)上,不利于2,3-丁二醇的大規(guī)模發(fā)酵生產(chǎn)。如何利用廉價糖源達(dá)到較為良好的發(fā)酵效果是當(dāng)今面臨的難題[5]。微生物好氧發(fā)酵過程是一種典型的多相反應(yīng)體系,該體系在化工、能源、環(huán)保等多個領(lǐng)域都有廣泛的運用,但無論是傳質(zhì)分離過程還是化學(xué)反應(yīng)過程,多相之間物質(zhì)傳輸?shù)乃俾识际侵萍s反應(yīng)速率和物質(zhì)利用率的關(guān)鍵因素[6-9],其根本原因在于氣液或氣液固之間的相界面積太小[10]。而在微生物發(fā)酵生產(chǎn)2,3-丁二醇的菌種生長階段,恰好需要較高的溶氧量來支持菌體的代謝[1]。微界面強化傳質(zhì)技術(shù)通過界面?zhèn)髻|(zhì)強化來實現(xiàn)過程強化?;谖⒔缑鎻娀砟睿芯空唛_發(fā)了一系列新型高效的反應(yīng)裝置,并成功應(yīng)用于工業(yè)過程[11-14]。

    微生物發(fā)酵過程中,氣泡是氣液間物質(zhì)傳輸?shù)闹匾d體,其在液層中的運動與傳質(zhì)過程直接決定了微生物發(fā)酵工藝過程的效率[15-16]。而對于單氣泡而言,其半徑?jīng)Q定了它的運動速率[17],半徑的變化受物質(zhì)傳輸與氣泡內(nèi)壓的控制,而傳質(zhì)速率又會受到運動速率與氣泡尺寸的顯著影響[18-19],顯然,氣泡的上升與傳質(zhì)是一個互相影響的耦合過程。因此,建立合適的數(shù)學(xué)模型,描述氣泡在體系內(nèi)上升與傳質(zhì)的過程,并與微界面強化技術(shù)相結(jié)合,對發(fā)酵罐的設(shè)計以及發(fā)酵工藝中的各種工況與參數(shù)的設(shè)定都有非常有意義的指導(dǎo)作用。本文在這一背景下,構(gòu)建多組分氣泡運動與傳質(zhì)耦合模型,并結(jié)合能量耗散理論,利用數(shù)值模擬的方式分析微界面體系對好氧發(fā)酵的強化過程,并選擇預(yù)設(shè)工況下的優(yōu)化參數(shù)。

    1 模型介紹

    1.1 氣泡上升與傳質(zhì)耦合模型

    田恒斗等[20-22]結(jié)合氣泡瞬態(tài)加速度模型與氣泡瞬態(tài)非平衡傳質(zhì)模型,以描述單組分氣泡在液體中的上升和傳質(zhì)作用,其核心方程組如下:

    其中第一個方程是描述氣泡瞬態(tài)加速度的微分方程。第二個是描述氣泡瞬態(tài)傳質(zhì)的微分方程,CA為液相中的氣體濃度,CI為氣液相界面處的氣體濃度。第四個是利用非平衡傳質(zhì)理論求解CI的公式,C*為氣液界面處氣體的平衡濃度,H 為氣體的溶解度系數(shù),σ為液體的表面張力。

    1.2 多組分氣相在混合溶液中的上升與傳質(zhì)耦合模型

    在工業(yè)過程中,多組分氣體構(gòu)成的反應(yīng)體系較為常見,特別是以空氣為反應(yīng)物的,如好氧發(fā)酵、PX氧化、濕式氧化等過程。多組分氣體傳質(zhì)與反應(yīng)過程中存在組分間的競爭與協(xié)調(diào)關(guān)系,并影響傳質(zhì)與反應(yīng)過程。因此,構(gòu)建氣泡運動傳質(zhì)耦合模型,以適用多組分氣泡在液體中的運動和傳質(zhì)行為,具有更為重要的工業(yè)意義。

    在含N 個組分的氣泡和K 個組分的液相體系中,令yi為氣相i 組分的摩爾分?jǐn)?shù),ρi為氣相i 組分的密度,Mg,i為氣相中i 組分的摩爾質(zhì)量,yl,i為液相i 組分的摩爾分?jǐn)?shù),μi為液相i 組分的黏度,Ml,i為液相i組分的摩爾質(zhì)量,可得到一系列數(shù)組:

    因此,混合氣體的密度可由式(8)計算:

    混合液體的黏度μl[23]可表示為:

    通過這一系列擴(kuò)展過程,可以用式(10)取代式(1)中第一行的微分方程,使描述氣泡瞬態(tài)加速度的方程適用于多組分。

    為了簡化問題,假設(shè)氣相各組分在邊界層的傳質(zhì)是互不干涉的,則單個氣泡的物質(zhì)傳輸速率可表示為:

    式中,DAB為氣體在液相中的擴(kuò)散系數(shù),I =[1 1 … 1]。

    任意時刻氣泡內(nèi)所含氣體質(zhì)量可表示為:

    式(12)兩邊對時間t求導(dǎo)可得:

    對于包含N個組分的氣泡,代入式(8),則有:

    設(shè)mg,i為氣泡中i 組分的質(zhì)量,Mg為該參數(shù)的合集數(shù)組;再設(shè)Ng=[n1n2… nN],ni是氣泡中i組分的物質(zhì)的量;因為Ng= Mg/Molg,所以:

    將式(10)、式(11)代入式(7)可得:

    用式(19)取代式(1)中第二行的微分方程,同樣使描述氣泡瞬態(tài)非平衡傳質(zhì)的方程適用于多組分氣泡。

    1.3 能量耗散模型

    能量耗散率(energy dissipation rate,ε),定義為單位質(zhì)量流體耗散能量的速率,單位為W·kg-1。基于本課題組前期對于微界面強化技術(shù)的研究[24-25],在超細(xì)氣泡控制塔中,氣泡產(chǎn)生過程中的能量耗散速率與氣泡體系的尺寸有很大的關(guān)系。可以用ε來表征微氣泡體系的能量消耗情況,其與氣泡尺寸的關(guān)系如下:

    式中,dmin和dmax分別代表體系中最小和最大的氣 泡 直 徑;d32為 氣 泡 的Satuer 直 徑;Wecrit為 臨 界Weber數(shù),通常取1.24[26]。

    2 模型計算方法

    本文利用Runge-Kutta 變階方法[27-28]進(jìn)行上述模型(ODE 系統(tǒng))數(shù)值模擬。Runge-Kutta 變階方法是基于Runge-Kutta 定階方法改進(jìn)的數(shù)值積分方法。對改進(jìn)后的式(1)使用該方法,引入初始條件,并輸入所需參數(shù),即可對任意t 時刻對應(yīng)的氣泡半徑、上浮速度、所在位置以及瞬時傳質(zhì)速率等參數(shù)實現(xiàn)快速數(shù)值求解。

    3 計算實例與結(jié)果分析

    模擬好氧發(fā)酵體系中微氣泡的上升運動及氧傳質(zhì)情況,并調(diào)節(jié)關(guān)鍵參數(shù),分析體系的傳質(zhì)與能耗情況,從而獲得較為理想的工況??諝庵饕裳鯕夂偷獨鈨煞N氣體組成,模擬計算中所用到的該兩種物質(zhì)的相關(guān)參數(shù)如表1所示。

    計算過程中用到的其他條件包括:微生物發(fā)酵過程的操作壓強為一個標(biāo)準(zhǔn)大氣壓(101325 Pa),溫度為20℃。液相密度ρL為1000 kg·m-3,液體的表面張力σ 為0.0728 N·m-1,動力黏度系數(shù)μ 為1.021×10-3Pa·s,液相氣體濃度CA=0,發(fā)酵罐液面高度為5 m。微界面強化的核心思想即利用微界面取代傳統(tǒng)的宏界面,從而成倍地提高傳質(zhì)速率,根據(jù)微界面的定義[8]選擇一系列初始?xì)馀莅霃絉0(0.1,0.3,0.5,0.8,1.0,1.5,2.0 mm)來考察微界面體系對傳質(zhì)效果的提升。初值條件設(shè)為氣泡高度h0=0,氣泡初始速度vb0=0。預(yù)設(shè)的運動時間為50 s。

    表1 體系中主要物質(zhì)的各項參數(shù)Table 1 Parameters of main substances in the system

    3.1 氣泡運動狀態(tài)

    氣泡半徑是決定其運動與傳質(zhì)狀態(tài)的關(guān)鍵參數(shù)[29],氣泡上升過程中,其大小會受到氣體傳輸和氣泡內(nèi)壓的制約。

    如圖1 所示,氣泡上升過程中,其半徑逐漸減小,這種氣泡尺寸的變化主要是氣液兩相間傳質(zhì)造成的,且初始半徑越大的氣泡尺寸變化越顯著。另外大氣泡由于運動速度較大,能夠在短時間內(nèi)迅速離開液層。氣泡上升高度以及上升速率的變化情況分別如圖2、圖3所示。

    圖1 氣泡半徑隨時間的變化情況Fig.1 Change of bubble radius with time

    圖2 氣泡高度隨時間的變化情況Fig.2 Change of bubble height with time

    圖3 氣泡速率隨時間的變化情況Fig.3 Change of bubble velocity with time

    在圖2 中,實心點所代表的氣泡在設(shè)定的上升時間內(nèi)會離開液面,而空心點表示的氣泡在預(yù)設(shè)時間內(nèi)將一直停留在液層內(nèi),由此可知,較小的氣泡所構(gòu)成的微界面體系可充分利用罐體高度,長時間停留在液層中,實現(xiàn)了充分傳質(zhì),而較大的氣泡極可能未達(dá)到充分傳質(zhì)前已從液體中逸出,導(dǎo)致物料的浪費。氣泡半徑是決定氣泡運動速率的關(guān)鍵因素,因此,氣泡上升速率變化的趨勢應(yīng)當(dāng)與氣泡半徑變化的趨勢相吻合,結(jié)合圖1 與圖3 來看,這二者之間變化的趨勢是符合預(yù)期的。

    3.2 氣泡傳質(zhì)情況

    氣泡內(nèi)的主要組分為N2與O2。由于氣泡大小對于氣液傳質(zhì)界面有較大的影響,為了評價各不同尺寸氣泡的傳質(zhì)效果,引入傳質(zhì)通量(mass transfer flux,J),定義為單個氣泡中的某組分在單位時間內(nèi)的質(zhì)量變化量與該時刻氣泡表面積的比值。此外,氧氣利用率(O2utilization,η),即傳輸至液相中的氧氣質(zhì)量占?xì)馀葜性醒鯕赓|(zhì)量的比值,也是衡量傳質(zhì)效果的重要指標(biāo)。兩種指標(biāo)的具體計算公式如下:

    式中,S為t時刻氣泡的表面積。

    以初始半徑為0.5 mm 的氣泡為例,表2 展示了氣泡內(nèi)兩種氣體的組成變化,圖4 為其傳質(zhì)通量的變化。

    由表2 可知,氣泡內(nèi)的氣體組成隨氣泡的運動與傳質(zhì)而改變,N2所占的摩爾分?jǐn)?shù)略有上升,O2的摩爾分?jǐn)?shù)下降,這是由于兩者擴(kuò)散差異所致。但O2在液相中的擴(kuò)散系數(shù)也僅是略高于N2,因此在傳質(zhì)速率上不會存在過大的差異。這在圖4 也有所反映,從整體趨勢上來看,兩種氣體的傳質(zhì)通量變化趨勢幾乎一致,但N2的傳質(zhì)通量要明顯高于O2。這是因為在空氣中,N2的物質(zhì)的量幾乎是O2的四倍,其在界面處存在更高的相界面濃度,從而導(dǎo)致更高的傳質(zhì)通量。但在實際的好氧發(fā)酵過程中,液相溶解的N2不被消耗,故存在傳質(zhì)平衡,而氧氣會被不斷消耗。

    表2 氣泡內(nèi)氣體摩爾分?jǐn)?shù)的變化Table 2 Changes of gas mole fraction in bubbles

    圖4 兩種氣體傳質(zhì)通量隨時間變化的情況Fig.4 Change of mass transfer flux of two gases with time

    對比不同大小氣泡中O2的傳輸效率是選擇合適初始?xì)馀莅霃降年P(guān)鍵步驟。

    圖5 初始半徑對O2傳質(zhì)通量的影響Fig.5 Effect of initial radius on mass transfer flux of O2

    圖6 初始半徑對O2利用率的影響Fig.6 Effect of initial radius on O2 utilization

    由圖5 可知,氣泡中氧的傳質(zhì)速率會迅速上升達(dá)到一個極值,而后緩慢下降。對于小氣泡,這種降低越慢,而且氣泡越小,能達(dá)到的傳質(zhì)通量極值就越高,這充分體現(xiàn)了微界面體系對氣液傳質(zhì)的強化特性。這得益于小氣泡較小的上升速率和較薄的傳質(zhì)邊界層。而對于整個體系而言,較小的氣泡群會擁有更大的相界面積[30],從而提高體系的氣含率和氧傳質(zhì)速率,對發(fā)酵過程同樣也能起到促進(jìn)作用。

    就物料成本考慮,大氣泡的攜帶氧氣能力強,處理量大,但傳質(zhì)通量較小。如圖6所示,其傳輸?shù)揭后w內(nèi)的O2占其原攜帶O2量的百分比明顯小于小氣泡;另外大氣泡上升速度快,停留時間短,也制約了O2在氣液間的傳輸。

    溶解氧(dissolved oxygen,DO)在發(fā)酵工藝中一般以相對含量表示,稱為DO 值,是發(fā)酵過程中的一項重要指標(biāo)。一般認(rèn)為N2吹掃后體系中的DO 值為0,而將氧氣充分溶解之后的DO 值視為100%。由于好氧發(fā)酵過程的各個階段存在不同的DO 值,研究DO 值對O2傳輸?shù)挠绊懸彩潜匾?,根?jù)非平衡理論計算所得的液相濃度CA與傳質(zhì)推動力CA-CI的關(guān)系如圖7所示。顯然,隨著CA的增大,傳質(zhì)推動力先增后減,存在最大值。

    為了便于表達(dá),設(shè)XA為液相中的氣體濃度與該條件下液相飽和濃度的比值,稱為液相飽和度,XI為界面濃度與飽和濃度的比值,稱為界面飽和度,兩者之差即為該條件下的傳質(zhì)推動力,從圖8 中可以看出,當(dāng)液相飽和度在0.2~0.3 時,傳質(zhì)推動力有最大值。以初始半徑為0.3 mm 的氣泡為例計算,考察不同DO 值下的傳質(zhì)通量與氧氣利用率,得到的圖9、圖10可以表明,計算結(jié)果與預(yù)測結(jié)果一致,當(dāng)DO值為20%時,有最好的傳質(zhì)效果; 當(dāng)DO 值達(dá)到100%之后,傳質(zhì)即停止,甚至?xí)霈F(xiàn)“倒吸”現(xiàn)象。

    圖7 不同液相濃度下的傳質(zhì)推動力Fig.7 Mass transfer driving force under different CA

    圖8 液相飽和度對傳質(zhì)推動力的影響Fig.8 Effect of XA on the mass transfer driving force

    圖9 DO值對O2利用率的影響Fig.9 Effect of DO value on O2 utilization

    3.3 微氣泡體系的能耗

    如前文所述,微氣泡在促進(jìn)氣液傳質(zhì)方面的作用顯著,但要形成小氣泡往往需要較高的能量消耗[31]。實際生產(chǎn)過程中不僅需要考慮物料成本,能量的損耗同樣也需要計算。

    圖11展示了O2傳質(zhì)效率、能量耗散系數(shù)與初始?xì)馀莅霃街g的關(guān)系。由圖可知,較高的傳質(zhì)速率與較低的能耗難以兼得,但一定存在某個區(qū)間,在能耗適中的情況下,獲得較為理想的傳質(zhì)效果。在圖11中,這個區(qū)間應(yīng)該在250 ~500 μm之間。

    圖10 DO值對O2傳質(zhì)通量的影響Fig.10 Effect of DO value on mass transfer flux of O2

    圖11 不同的初始半徑對傳質(zhì)以及能量耗散系數(shù)的影響Fig.11 Effect of different R0 on mass transfer and energy dissipation coefficient

    但是,本文的研究還是在一個高度理想化的情況下進(jìn)行的,氣泡在水平方向上是否會有運動,多氣泡之間的互相作用等一系列的問題都需要進(jìn)一步的研究,現(xiàn)有的理論成果是否可靠也有待實驗的論證。

    4 結(jié) 論

    針對好氧發(fā)酵過程,本文將原有單組分氣泡上升與傳質(zhì)模型擴(kuò)展到多組分氣泡傳質(zhì)過程,考察氣泡的不同初始半徑對不同組分傳質(zhì)的影響,結(jié)合能量耗散系數(shù)這一指標(biāo),評估體系的能耗情況,選擇出理論條件下較為合適的工況條件。

    計算結(jié)果表明,當(dāng)在20℃,1 個標(biāo)準(zhǔn)大氣壓(101325 Pa)的標(biāo)準(zhǔn)狀況下操作時,對于一個液面高度為5 m的發(fā)酵罐,當(dāng)氣泡初始半徑超過0.5 mm時,其將在較短時間內(nèi)逸出體系,造成物料浪費;較小的氣泡則擁有更好的傳質(zhì)效果,當(dāng)氣泡初始半徑小于0.1 mm 時,其對攜帶的O2的利用率可以接近50%,顯示了微界面體系對于氣液傳質(zhì)的強化特性。

    DO 值是好氧發(fā)酵工藝中的一項重要指標(biāo),通過模型的預(yù)測,在預(yù)設(shè)工況下,將體系內(nèi)的DO 值維持在20%~30%之間時,會取得較為良好的傳質(zhì)效果。

    由于產(chǎn)生微氣泡的能耗需要較高的能耗,因此一味地追求小氣泡并不能保證經(jīng)濟(jì)效益。在不同工況下,綜合考慮物料成本和能耗成本,選擇合適的分散尺寸才能帶來最佳的優(yōu)化效果。

    符 號 說 明

    C——氣體溶解于液相中的濃度,kg·m-3

    C*——氣液界面處的平衡濃度,kg·m-3

    CA——液相主體中的氣體濃度,kg·m-3

    CD——黏性阻力系數(shù)

    CI——氣液界面處的氣體濃度,kg·m-3

    DAB——氣體分子在液體中擴(kuò)散速率,m2·s-1

    dmax——氣泡群中氣泡的最大直徑,m

    dmin——氣泡群中氣泡的最小直徑,m

    d32——氣泡群的Sauter直徑,m

    E——亨利系數(shù),kPa

    g——重力加速度,N·kg-1

    H——溶解度系數(shù),kg·m-3·Pa-1

    h——氣泡上升高度,m

    J——傳質(zhì)通量,kg·s-1·m-2

    M——摩爾質(zhì)量,g·mol-1

    Mg——組分質(zhì)量合集數(shù)組,kg

    Mol——摩爾質(zhì)量的合集數(shù)組,g·mol-1

    mg——氣泡質(zhì)量,kg

    mO2——氣泡內(nèi)氧氣質(zhì)量,kg

    N——物質(zhì)的量合集數(shù)組,mol

    patm——體系外界壓強,Pa

    pg——氣泡內(nèi)部壓強,Pa

    R——氣泡半徑,m

    S——氣泡表面積,m2

    t——氣泡上升運動時間,s

    Vg——氣泡體積,m3

    vb——氣泡上升速率,m·s-1

    Wecrit——臨界Weber數(shù)

    XA——液相飽和度

    XI——界面飽和度

    Y——組分摩爾分率的合集數(shù)組

    y——組分摩爾分?jǐn)?shù)

    ε——能量耗散率,W·kg-1

    η——氧氣利用率,%

    μ——動力黏度系數(shù),Pa·s

    ρg——氣泡內(nèi)氣體密度,kg·m-3

    ρL——液體密度,kg·m-3

    ρ0——外界壓強下對應(yīng)氣體的密度,kg·m-3

    σ——表面張力,N·m-1

    下角標(biāo)

    g——氣相

    l——液相

    0——初始值

    猜你喜歡
    傳質(zhì)氣液氣泡
    檸檬氣泡水
    欣漾(2024年2期)2024-04-27 15:19:49
    SIAU詩杭便攜式氣泡水杯
    新潮電子(2021年7期)2021-08-14 15:53:12
    浮法玻璃氣泡的預(yù)防和控制對策
    微重力下兩相控溫型儲液器內(nèi)氣液界面仿真分析
    氣液分離罐液位計接管泄漏分析
    冰凍氣泡
    CO2 驅(qū)低液量高氣液比井下氣錨模擬與優(yōu)化
    氨基酸鹽吸收二氧化碳過程的傳質(zhì)特性
    PTFE膜吸收CO2的工藝參數(shù)對傳質(zhì)性能的影響
    清潔轉(zhuǎn)向酸H+表面?zhèn)髻|(zhì)行為實驗研究
    寂寞人妻少妇视频99o| 男女之事视频高清在线观看| 啦啦啦韩国在线观看视频| 久久久久久伊人网av| 亚洲久久久久久中文字幕| 久久久精品欧美日韩精品| 国产综合懂色| 婷婷六月久久综合丁香| 蜜桃久久精品国产亚洲av| 国产一区二区三区av在线 | 国产精品精品国产色婷婷| 国产精品不卡视频一区二区| 久久精品影院6| 可以在线观看毛片的网站| 99热只有精品国产| 国产高清三级在线| 国产蜜桃级精品一区二区三区| 亚洲第一电影网av| 日日撸夜夜添| 国内精品久久久久精免费| 91在线精品国自产拍蜜月| 嫩草影院入口| 淫秽高清视频在线观看| 午夜影院日韩av| 在线观看66精品国产| 国产精品电影一区二区三区| 亚洲欧美中文字幕日韩二区| 国产午夜精品久久久久久一区二区三区 | 久久精品国产清高在天天线| 精品一区二区三区av网在线观看| 亚洲国产欧美人成| 久久亚洲国产成人精品v| 日本a在线网址| 亚洲av二区三区四区| 1024手机看黄色片| 中文字幕精品亚洲无线码一区| 午夜亚洲福利在线播放| 久久久久久久久大av| 国产爱豆传媒在线观看| 日本熟妇午夜| 久久久久免费精品人妻一区二区| 毛片一级片免费看久久久久| 两个人的视频大全免费| 女人被狂操c到高潮| 国内揄拍国产精品人妻在线| 成人美女网站在线观看视频| 国产一区二区三区在线臀色熟女| 看免费成人av毛片| 亚洲va在线va天堂va国产| 精品人妻偷拍中文字幕| 久久精品国产99精品国产亚洲性色| 国产精品av视频在线免费观看| 精品久久久久久久久久免费视频| 哪里可以看免费的av片| 中文资源天堂在线| 亚洲婷婷狠狠爱综合网| 天天躁夜夜躁狠狠久久av| 美女高潮的动态| 国产成人91sexporn| 免费观看精品视频网站| 欧美一区二区亚洲| 免费无遮挡裸体视频| 亚洲av免费高清在线观看| 日韩av不卡免费在线播放| 99热这里只有精品一区| 亚洲自偷自拍三级| 中文资源天堂在线| 最后的刺客免费高清国语| 91久久精品电影网| 成人鲁丝片一二三区免费| 99久久精品热视频| 午夜精品一区二区三区免费看| 成年女人永久免费观看视频| 人人妻人人看人人澡| 校园人妻丝袜中文字幕| 午夜日韩欧美国产| 91在线精品国自产拍蜜月| 高清午夜精品一区二区三区 | 听说在线观看完整版免费高清| 亚洲色图av天堂| 能在线免费观看的黄片| 日本色播在线视频| 亚洲在线自拍视频| 亚洲性夜色夜夜综合| h日本视频在线播放| 欧美xxxx性猛交bbbb| 久久精品国产亚洲av香蕉五月| 欧美最黄视频在线播放免费| 日本a在线网址| 欧美一区二区精品小视频在线| 精品日产1卡2卡| 国产黄a三级三级三级人| 黄色一级大片看看| 热99re8久久精品国产| 久久99热这里只有精品18| 97超碰精品成人国产| 性色avwww在线观看| 国产三级在线视频| 亚洲在线自拍视频| 色播亚洲综合网| 国产三级中文精品| 午夜久久久久精精品| 天美传媒精品一区二区| 国内少妇人妻偷人精品xxx网站| 亚洲不卡免费看| 中文字幕熟女人妻在线| 精品熟女少妇av免费看| 亚洲欧美精品综合久久99| 精品久久国产蜜桃| 在线观看免费视频日本深夜| 性欧美人与动物交配| 床上黄色一级片| av天堂中文字幕网| 麻豆av噜噜一区二区三区| 精品午夜福利在线看| 亚洲色图av天堂| 啦啦啦啦在线视频资源| 国产精品国产三级国产av玫瑰| 欧美区成人在线视频| 午夜福利18| 91久久精品电影网| 精品熟女少妇av免费看| 在线观看66精品国产| 少妇高潮的动态图| 免费人成在线观看视频色| 久久午夜亚洲精品久久| 亚洲av.av天堂| 国语自产精品视频在线第100页| 国产v大片淫在线免费观看| 我要搜黄色片| 亚洲精品一区av在线观看| 国产激情偷乱视频一区二区| 日本欧美国产在线视频| 香蕉av资源在线| 亚洲经典国产精华液单| 欧美区成人在线视频| 免费av不卡在线播放| 老熟妇乱子伦视频在线观看| 日本熟妇午夜| 日本熟妇午夜| 伦精品一区二区三区| 欧美区成人在线视频| 久久人妻av系列| 亚洲精品国产av成人精品 | 国产麻豆成人av免费视频| 国产精品三级大全| 一进一出抽搐gif免费好疼| 亚洲在线自拍视频| 午夜精品在线福利| 美女黄网站色视频| 全区人妻精品视频| 级片在线观看| 美女cb高潮喷水在线观看| 国产午夜精品久久久久久一区二区三区 | 最近的中文字幕免费完整| 97在线视频观看| 日韩av不卡免费在线播放| 亚洲人成网站在线观看播放| 亚洲国产精品成人久久小说 | 2021天堂中文幕一二区在线观| 97超碰精品成人国产| 麻豆精品久久久久久蜜桃| 俄罗斯特黄特色一大片| 国产三级在线视频| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲成a人片在线一区二区| 亚洲性夜色夜夜综合| 国产男靠女视频免费网站| 欧美日韩国产亚洲二区| 夜夜夜夜夜久久久久| 国产精品女同一区二区软件| 亚洲美女黄片视频| 久久久久国产网址| 亚洲国产精品久久男人天堂| 国模一区二区三区四区视频| 久久精品国产清高在天天线| 久久精品91蜜桃| 亚洲欧美日韩高清在线视频| 1024手机看黄色片| 九九热线精品视视频播放| 国产男靠女视频免费网站| 欧美绝顶高潮抽搐喷水| 中出人妻视频一区二区| 久久韩国三级中文字幕| 丰满的人妻完整版| 精品熟女少妇av免费看| 国产欧美日韩精品一区二区| 热99re8久久精品国产| 国产精品精品国产色婷婷| 一级a爱片免费观看的视频| 成年版毛片免费区| 亚洲天堂国产精品一区在线| 国产精品嫩草影院av在线观看| 国产精品三级大全| 欧美日本亚洲视频在线播放| 国产精品乱码一区二三区的特点| 精品99又大又爽又粗少妇毛片| 亚洲精品久久国产高清桃花| 国产男人的电影天堂91| 在线观看av片永久免费下载| 尾随美女入室| 女的被弄到高潮叫床怎么办| 国模一区二区三区四区视频| 国产av不卡久久| 亚洲真实伦在线观看| 久久久久国产精品人妻aⅴ院| 久久精品国产鲁丝片午夜精品| 男人的好看免费观看在线视频| 日韩三级伦理在线观看| 大型黄色视频在线免费观看| 亚洲,欧美,日韩| 又爽又黄a免费视频| 在线免费十八禁| 少妇裸体淫交视频免费看高清| 婷婷六月久久综合丁香| 国产高潮美女av| 色5月婷婷丁香| 久久久久精品国产欧美久久久| 国产片特级美女逼逼视频| 好男人在线观看高清免费视频| 免费观看精品视频网站| 女同久久另类99精品国产91| 久久久久免费精品人妻一区二区| 亚洲久久久久久中文字幕| 亚洲欧美日韩高清专用| 免费观看人在逋| 99久久中文字幕三级久久日本| 大型黄色视频在线免费观看| 又黄又爽又免费观看的视频| 亚洲在线自拍视频| 欧美中文日本在线观看视频| 久久久久九九精品影院| 久久草成人影院| 国产色婷婷99| 亚洲无线观看免费| 国语自产精品视频在线第100页| 身体一侧抽搐| 校园人妻丝袜中文字幕| 久久综合国产亚洲精品| 久久精品影院6| 热99在线观看视频| 亚洲精品在线观看二区| 中文亚洲av片在线观看爽| 亚洲一区高清亚洲精品| 一级毛片电影观看 | 一区福利在线观看| 91午夜精品亚洲一区二区三区| 黑人高潮一二区| 舔av片在线| 欧美xxxx黑人xx丫x性爽| 成人av在线播放网站| 国产精品伦人一区二区| 男女那种视频在线观看| 国产伦在线观看视频一区| 91在线观看av| av专区在线播放| 午夜久久久久精精品| 中出人妻视频一区二区| 狂野欧美白嫩少妇大欣赏| 99久久无色码亚洲精品果冻| 乱人视频在线观看| 97超视频在线观看视频| 日韩欧美精品v在线| 日本免费一区二区三区高清不卡| 大香蕉久久网| 日韩成人伦理影院| 久久久a久久爽久久v久久| 蜜桃亚洲精品一区二区三区| 国产精品一区www在线观看| 欧美激情在线99| 乱人视频在线观看| 国产欧美日韩精品亚洲av| 22中文网久久字幕| 亚洲欧美成人精品一区二区| 最近手机中文字幕大全| 女人十人毛片免费观看3o分钟| 亚洲精品国产成人久久av| 日韩三级伦理在线观看| 亚洲欧美成人综合另类久久久 | 嫩草影院新地址| 欧美一区二区国产精品久久精品| 我的女老师完整版在线观看| 欧美日韩精品成人综合77777| 国产精品久久久久久亚洲av鲁大| 日本-黄色视频高清免费观看| 欧美不卡视频在线免费观看| 国产午夜精品论理片| 波多野结衣巨乳人妻| 成人综合一区亚洲| av黄色大香蕉| 亚洲av电影不卡..在线观看| 最新在线观看一区二区三区| 熟女人妻精品中文字幕| 男人和女人高潮做爰伦理| 97热精品久久久久久| av在线播放精品| 日本免费a在线| 91久久精品国产一区二区三区| 精品午夜福利在线看| 免费av毛片视频| 在线观看一区二区三区| 亚洲自拍偷在线| 日日啪夜夜撸| 天天躁夜夜躁狠狠久久av| 悠悠久久av| 三级国产精品欧美在线观看| 高清午夜精品一区二区三区 | 国产激情偷乱视频一区二区| 成人综合一区亚洲| 国产在线精品亚洲第一网站| 一本精品99久久精品77| 日本在线视频免费播放| 精品日产1卡2卡| 成熟少妇高潮喷水视频| 老司机午夜福利在线观看视频| 久久99热6这里只有精品| 韩国av在线不卡| 一级黄片播放器| 免费在线观看成人毛片| 国产日本99.免费观看| 精品无人区乱码1区二区| 非洲黑人性xxxx精品又粗又长| 国产黄色视频一区二区在线观看 | 少妇的逼水好多| 小蜜桃在线观看免费完整版高清| 久久99热6这里只有精品| 亚洲人与动物交配视频| 国产黄色视频一区二区在线观看 | 欧美日本视频| 国产精品不卡视频一区二区| 亚洲无线观看免费| 久久久久久久久久成人| 国产精品99久久久久久久久| 国产探花极品一区二区| 搞女人的毛片| 国产精品国产高清国产av| 精品福利观看| 三级经典国产精品| 女同久久另类99精品国产91| 在线天堂最新版资源| 国产成人精品久久久久久| 国产女主播在线喷水免费视频网站 | 国产精品久久久久久av不卡| 国产精品久久久久久av不卡| 一级毛片我不卡| 亚洲精品国产av成人精品 | 特级一级黄色大片| 成人高潮视频无遮挡免费网站| 国产毛片a区久久久久| 久久久久国产网址| 日本五十路高清| 性色avwww在线观看| 级片在线观看| 婷婷精品国产亚洲av| 看非洲黑人一级黄片| 亚洲av第一区精品v没综合| 麻豆av噜噜一区二区三区| 免费av观看视频| 亚洲av五月六月丁香网| 亚洲一区二区三区色噜噜| 亚洲激情五月婷婷啪啪| 深夜a级毛片| 国产私拍福利视频在线观看| АⅤ资源中文在线天堂| 偷拍熟女少妇极品色| 免费搜索国产男女视频| 99久国产av精品国产电影| 麻豆乱淫一区二区| 亚洲久久久久久中文字幕| 国产色爽女视频免费观看| 亚洲人成网站在线观看播放| 国产成人freesex在线 | 韩国av在线不卡| av.在线天堂| 免费搜索国产男女视频| 小蜜桃在线观看免费完整版高清| 亚洲av免费在线观看| 国产精品久久电影中文字幕| 男人舔奶头视频| 99久国产av精品国产电影| 少妇高潮的动态图| 欧美一区二区亚洲| 99热精品在线国产| 国产精品av视频在线免费观看| 亚洲国产精品成人久久小说 | 国产极品精品免费视频能看的| 听说在线观看完整版免费高清| 日韩欧美国产在线观看| av在线观看视频网站免费| 久久午夜亚洲精品久久| 一个人看的www免费观看视频| 嫩草影院精品99| 男女那种视频在线观看| 精华霜和精华液先用哪个| 简卡轻食公司| 久久久久久久久久成人| 亚洲自偷自拍三级| 欧美成人精品欧美一级黄| 久久久久久九九精品二区国产| 免费搜索国产男女视频| 午夜免费男女啪啪视频观看 | 亚洲av电影不卡..在线观看| 老熟妇仑乱视频hdxx| 一区二区三区免费毛片| 在线国产一区二区在线| 久久久久久久久久黄片| 欧美3d第一页| 亚洲一区二区三区色噜噜| 一夜夜www| 欧美潮喷喷水| 色综合亚洲欧美另类图片| 秋霞在线观看毛片| 国产老妇女一区| 美女黄网站色视频| 亚洲人与动物交配视频| 精品少妇黑人巨大在线播放 | 亚洲在线自拍视频| 欧美日韩精品成人综合77777| 久久九九热精品免费| 中文资源天堂在线| 国产高清激情床上av| 午夜激情欧美在线| 一级黄色大片毛片| 我要看日韩黄色一级片| 欧美高清成人免费视频www| 黄色视频,在线免费观看| 精品少妇黑人巨大在线播放 | 中文字幕熟女人妻在线| 日韩国内少妇激情av| 亚洲aⅴ乱码一区二区在线播放| 亚洲精品一卡2卡三卡4卡5卡| 俺也久久电影网| 色哟哟哟哟哟哟| 美女 人体艺术 gogo| 国产av在哪里看| 精品午夜福利视频在线观看一区| 亚洲综合色惰| 美女被艹到高潮喷水动态| 露出奶头的视频| 不卡视频在线观看欧美| 亚洲精品久久国产高清桃花| ponron亚洲| 欧美丝袜亚洲另类| 欧美bdsm另类| 国产综合懂色| 一级毛片电影观看 | 欧美成人免费av一区二区三区| 亚洲精品一区av在线观看| 国产片特级美女逼逼视频| 国产色婷婷99| АⅤ资源中文在线天堂| 搡老熟女国产l中国老女人| 99国产精品一区二区蜜桃av| 成人特级黄色片久久久久久久| 成人性生交大片免费视频hd| 国产真实伦视频高清在线观看| 免费黄网站久久成人精品| 国产私拍福利视频在线观看| a级一级毛片免费在线观看| 亚洲性夜色夜夜综合| 久久久久国产精品人妻aⅴ院| 日韩欧美精品免费久久| 22中文网久久字幕| 晚上一个人看的免费电影| 少妇被粗大猛烈的视频| 男女那种视频在线观看| 激情 狠狠 欧美| 欧美激情久久久久久爽电影| 插阴视频在线观看视频| 亚洲激情五月婷婷啪啪| 国产91av在线免费观看| 精品午夜福利在线看| 日韩在线高清观看一区二区三区| 免费搜索国产男女视频| 成人av一区二区三区在线看| 亚洲一区高清亚洲精品| 一进一出好大好爽视频| 亚洲av成人精品一区久久| 久久久久久大精品| 少妇高潮的动态图| 精品欧美国产一区二区三| 日产精品乱码卡一卡2卡三| 欧美另类亚洲清纯唯美| 欧美精品国产亚洲| 精品少妇黑人巨大在线播放 | 亚洲av熟女| 日本黄色视频三级网站网址| 韩国av在线不卡| 日韩精品有码人妻一区| 青春草视频在线免费观看| 久久久久久久久久黄片| 久久久久久久亚洲中文字幕| 深夜精品福利| 亚洲国产色片| 国产大屁股一区二区在线视频| 男女做爰动态图高潮gif福利片| 亚洲精品456在线播放app| 成人三级黄色视频| 麻豆精品久久久久久蜜桃| 村上凉子中文字幕在线| 性欧美人与动物交配| 深夜a级毛片| 男人狂女人下面高潮的视频| 精品人妻视频免费看| 一级av片app| 国内揄拍国产精品人妻在线| 别揉我奶头 嗯啊视频| 99热只有精品国产| 国产精品野战在线观看| 超碰av人人做人人爽久久| 欧美日本视频| 99国产精品一区二区蜜桃av| av天堂中文字幕网| 国产亚洲av嫩草精品影院| 1000部很黄的大片| 99热全是精品| 3wmmmm亚洲av在线观看| 国产乱人偷精品视频| 天堂√8在线中文| 国产亚洲精品综合一区在线观看| 18禁在线播放成人免费| 不卡视频在线观看欧美| 欧美成人一区二区免费高清观看| 在线天堂最新版资源| 国产美女午夜福利| 亚洲精品日韩在线中文字幕 | 日本撒尿小便嘘嘘汇集6| 久久国产乱子免费精品| 欧美成人免费av一区二区三区| 色噜噜av男人的天堂激情| 亚洲欧美精品综合久久99| 3wmmmm亚洲av在线观看| 精品99又大又爽又粗少妇毛片| 香蕉av资源在线| 99热精品在线国产| 最近视频中文字幕2019在线8| 成人永久免费在线观看视频| 欧美zozozo另类| 精品福利观看| videossex国产| 欧美区成人在线视频| 禁无遮挡网站| 身体一侧抽搐| 麻豆成人午夜福利视频| 三级国产精品欧美在线观看| 成人av在线播放网站| 中文字幕av成人在线电影| 久久人人爽人人片av| 欧美日韩在线观看h| 极品教师在线视频| 热99在线观看视频| 天堂影院成人在线观看| 久久久久久久久久黄片| 亚洲三级黄色毛片| 国产一区二区三区在线臀色熟女| 中国美白少妇内射xxxbb| 搡老妇女老女人老熟妇| 国产精品一区二区三区四区免费观看 | 麻豆国产av国片精品| 日韩人妻高清精品专区| 麻豆av噜噜一区二区三区| 一级a爱片免费观看的视频| 内地一区二区视频在线| 日日撸夜夜添| 日本一本二区三区精品| 综合色丁香网| 国产精品一区二区性色av| 校园春色视频在线观看| 成人特级av手机在线观看| 日本成人三级电影网站| 女人被狂操c到高潮| 国产亚洲91精品色在线| 三级经典国产精品| 久久精品国产亚洲网站| 黑人高潮一二区| 久久久久国产网址| 亚洲人与动物交配视频| 国产精品一及| 日韩国内少妇激情av| 乱系列少妇在线播放| 国产精品av视频在线免费观看| 18+在线观看网站| 国产亚洲av嫩草精品影院| 欧美日韩精品成人综合77777| 真人做人爱边吃奶动态| 一本久久中文字幕| 人人妻人人澡欧美一区二区| 精品一区二区三区av网在线观看| 欧美激情在线99| 久久欧美精品欧美久久欧美| 99热这里只有是精品50| 丝袜美腿在线中文| 蜜臀久久99精品久久宅男| 天堂网av新在线| 欧美在线一区亚洲| 神马国产精品三级电影在线观看| 国产色爽女视频免费观看| av视频在线观看入口| 亚洲av中文av极速乱| 日韩欧美国产在线观看| 亚洲七黄色美女视频| 欧美+亚洲+日韩+国产| 男女那种视频在线观看| 卡戴珊不雅视频在线播放| 日产精品乱码卡一卡2卡三| 国内精品久久久久精免费| 18禁在线无遮挡免费观看视频 | 精品久久久久久久久久久久久| 免费大片18禁| 精品人妻偷拍中文字幕| av在线老鸭窝| 色播亚洲综合网| 亚洲美女黄片视频| 美女被艹到高潮喷水动态| 在线观看午夜福利视频|