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

    基于熱聲解耦方法的燃燒不穩(wěn)定性預(yù)測

    2022-03-22 03:24:56高宇鵬李俊杰劉巍李敬軒
    實驗流體力學(xué) 2022年1期

    高宇鵬,李俊杰 ,劉巍,李敬軒

    1.北京動力機(jī)械研究所,北京 100074 2.北京航空航天大學(xué) 宇航學(xué)院,北京 100191

    0 引 言

    燃燒不穩(wěn)定現(xiàn)象廣泛存在于火箭發(fā)動機(jī)、航空發(fā)動機(jī)以及大型空氣加熱器等裝置中,會使燃燒器內(nèi)產(chǎn)生不可控的巨大壓力脈動,造成發(fā)動機(jī)回火、熄火,燃燒器壁面過熱,導(dǎo)致發(fā)動機(jī)壽命縮短甚至損毀等一系列問題。20世紀(jì)60年代,在對F-1火箭發(fā)動機(jī)初始方案的44次試車中,超過半數(shù)的試車都發(fā)生了燃燒不穩(wěn)定現(xiàn)象;其振幅甚至達(dá)到了燃燒室平均壓力的100%。因此,為盡可能減少不穩(wěn)定燃燒帶來的巨大危害,通常采用實驗和數(shù)值仿真等方法對燃燒不穩(wěn)定性機(jī)理進(jìn)行研究。

    王迪等對單噴嘴模型發(fā)動機(jī)進(jìn)行了實驗研究,測量了不同燃燒室和噴嘴縮進(jìn)長度對燃燒穩(wěn)定性的影響,結(jié)果發(fā)現(xiàn)增加縮進(jìn)長度會對縱向高頻燃燒不穩(wěn)定性產(chǎn)生阻尼作用,燃燒室長度對縱向高頻燃燒不穩(wěn)定性的影響比縮進(jìn)長度更為明顯。Richecoeur等對液氧/氣態(tài)甲烷多噴嘴模型燃燒室進(jìn)行了試驗研究,結(jié)果發(fā)現(xiàn)橫向壓力震蕩與燃燒和火焰之間存在強(qiáng)烈的耦合作用,可增大火焰膨脹率以及火焰向壁面的傳熱率。Bibik等拍攝了燃燒室內(nèi)部氣流的旋轉(zhuǎn)過程,發(fā)現(xiàn)氣流旋轉(zhuǎn)方向與中心噴嘴旋轉(zhuǎn)相同,認(rèn)為中心噴嘴是激發(fā)切向不穩(wěn)定的重要因素。Méry等實驗研究了液氧/甲烷燃燒室在存在擾動的條件下,聲學(xué)擾動與火焰之間的相互作用。結(jié)果表明:在中等擾動條件下,液體射流存在橫向周期性振蕩;大幅度擾動時,火焰隨聲學(xué)擾動而周期性振蕩。Sliphorst等通過實驗研究發(fā)現(xiàn),燃燒不穩(wěn)定激勵主要來自于火焰與聲學(xué)振蕩的耦合作用。

    Urbano等通過大渦模擬(Large Eddy Simulation, LES)對功率為80 MW的全尺寸燃燒器進(jìn)行仿真,結(jié)果發(fā)現(xiàn)一階切向和一階縱向不穩(wěn)定對燃燒室的不穩(wěn)定性起重要作用,并且與混合燃料的噴注過程強(qiáng)烈耦合。Shimizu等等對液體火箭發(fā)動機(jī)(Liquid Rocket Engine, LRE)進(jìn)行了切向數(shù)值仿真,認(rèn)為平均流動分布控制一階切向不穩(wěn)定的發(fā)生。Staffelbach等對自激振蕩的環(huán)形燃燒室進(jìn)行了大渦模擬,成功捕捉了燃燒室自激振蕩頻率,發(fā)現(xiàn)其特征是由兩種不同振幅的旋轉(zhuǎn)模式疊加而成。王方、張志浩和夏朝陽等均采用LES對發(fā)動機(jī)內(nèi)燃燒不穩(wěn)定性進(jìn)行詳細(xì)的研究。Li等對ORACLESL長火焰燃燒室進(jìn)行了預(yù)測。采用大渦模擬獲得系統(tǒng)的火焰?zhèn)鬟f函數(shù),并將此傳遞函數(shù)帶入到Li和Morgans等開發(fā)的燃燒不穩(wěn)定預(yù)測和模擬程序—OSCILOS中,最終得到了精準(zhǔn)的熱–聲振蕩頻率和幅值的預(yù)測值。

    綜上所述,前人已通過實驗和數(shù)值的方法對熱聲耦合不穩(wěn)定燃燒進(jìn)行了研究,但是大尺寸燃燒器燃燒不穩(wěn)定性的預(yù)測一直存在計算量太大的問題。本文將通過熱–聲解耦的方法,對某型空氣加熱器燃燒不穩(wěn)定性進(jìn)行預(yù)測,為加熱器熱試提供技術(shù)支撐。

    1 燃燒-聲學(xué)解耦方法

    本文基于Li等的方法,將燃燒對擾動的響應(yīng)計算與聲學(xué)計算解耦,即分別計算火焰對聲學(xué)擾動的響應(yīng)(通常表征為火焰?zhèn)鬟f函數(shù))和聲學(xué)系統(tǒng)在熱源擾動下的響應(yīng)。因此,對于燃燒的數(shù)值模擬只需考慮燃燒室內(nèi)的火焰對來流擾動的響應(yīng),不需要考慮包括集氣環(huán)和噴管在內(nèi)的整個發(fā)動機(jī)系統(tǒng),顯著降低了計算量,提高了計算速度。

    1.1 噴管聲學(xué)邊界條件計算

    在聲學(xué)邊界條件計算中,除了噴管外,都可處理為壅塞或聲學(xué)閉口。由于噴管復(fù)雜的頻率特性,需要考慮其聲學(xué)邊界條件隨頻率的變化,通常表征為傳遞函數(shù)的形式。一般拉瓦爾噴管的聲學(xué)系統(tǒng)可以表示成如圖1所示模型。當(dāng)噴管上游來流的物性參數(shù)發(fā)生擾動時,下游氣流會產(chǎn)生壓力與速度的振蕩,產(chǎn)

    圖1 噴管的聲學(xué)模型Fig.1 The acoustic model of nozzle

    對于亞聲速流噴管(圖1),將向下游傳播的聲波施加在入口處,向上游傳播的聲波施加在出口處。假設(shè)噴管中的氣體是理想氣體,建立一維線性化歐拉方程如下:

    式中,u為氣體速度,p為氣體壓力,為比熱比,s為氣體熵,為 氣體定壓比熱容,為當(dāng)?shù)芈曀?,上?biāo)“”為擾動量。

    將上述控制方程進(jìn)行無量綱化處理,當(dāng)馬赫數(shù)不為1時,進(jìn)行傅里葉變換可以得到:

    式中,=/L為 噴管長度,x為距噴管入口的軸向距離;矩陣是由不變量構(gòu)成的矢量。矩陣和的形式分別如下:

    利用Magnus展開法可以求得式(2)的解,并求得最終的聲學(xué)反射系數(shù)。

    1.2 數(shù)值計算方法

    1.2.1 燃燒室三維建模

    加熱器整體模型如圖2所示,包含噴注面板、燃燒室、噴管等部件。其中,噴注面板上包含48個燃料噴嘴和36個冷卻空氣噴嘴。燃燒室直徑為350 mm,長度為800 mm。在對比了考慮和不考慮集氣環(huán)的燃燒室聲學(xué)特性后,發(fā)現(xiàn)集氣環(huán)對燃燒室不穩(wěn)定模態(tài)的穩(wěn)定性影響基本可以忽略。且由于各種連接件的壅塞特性,因此,在燃燒不穩(wěn)定性計算中忽略了集氣環(huán)這一部件。同時,由于噴注面板上噴嘴分布的高度對稱性,將采用1/6模型進(jìn)行計算,如圖3、4所示。

    圖2 加熱器各部件示意圖Fig.2 Schematic diagram of the heater components

    圖3 燃燒室模型示意圖Fig.3 Schematic diagram of combustion chamber model

    為確保火焰相關(guān)參數(shù)計算準(zhǔn)確,對噴注面板下游200 mm范圍內(nèi)的網(wǎng)格進(jìn)行加密,其各向尺寸小于1 mm。經(jīng)驗證計算,在網(wǎng)格量達(dá)到265萬時,滿足網(wǎng)格無關(guān)性要求;瞬態(tài)計算時間步長為10s時,滿足時間無關(guān)性要求,此時全局最大Courant數(shù)小于5。

    1.2.2 模型選擇和參數(shù)設(shè)置

    由于燃燒室內(nèi)總溫為1 743 K,總壓達(dá)8.2 MPa,當(dāng)預(yù)混燃料采用乙醇/液氧時,二者均處于超臨界態(tài)。超臨界流體流出噴嘴后,其氣液界面較為模糊,密度為液態(tài)的1/3,而黏度只有液態(tài)的1/12~1/4,因此其流動行為更偏向于氣體。在數(shù)值仿真中,燃料按連續(xù)相進(jìn)行處理。

    圖4 燃燒室網(wǎng)格劃分示意圖Fig.4 Divided mesh of combustion chamber model

    采用ANSYS Fluent對燃燒室內(nèi)的燃燒過程進(jìn)行數(shù)值仿真,其中內(nèi)部湍流流動采用Realizable模型求解;燃燒采用Eddy-Disspation模型求解。該燃燒模型假定燃燒速率遠(yuǎn)大于燃料混合速率,即燃料一經(jīng)混合便完全燃燒,忽略復(fù)雜的化學(xué)反應(yīng)速率,采用乙醇一步總包反應(yīng)機(jī)理:

    為確保在高溫高壓的燃燒室環(huán)境下,燃料物性參數(shù)(特別是對當(dāng)?shù)芈曀儆绊憳O大的密度)數(shù)值準(zhǔn)確,采用Peng-Robinson狀態(tài)方程計算。燃料入口參數(shù)如表1所示,各噴嘴入口邊界條件采用質(zhì)量流量入口,燃燒室出口邊界條件則采用壓力出口。

    表1 預(yù)混燃料入口工況Table 1 The working conditions of premixed fuel at the inlet

    得到燃燒室穩(wěn)態(tài)場后,在燃料入口處施加不同頻率、振幅的正弦擾動(式(6)),即可獲得各噴嘴的熱釋放率隨時間的變化關(guān)系,進(jìn)而得到各噴嘴的火焰?zhèn)鬟f函數(shù):

    1.2.3 火焰?zhèn)鬟f函數(shù)的定義

    通過對每個噴嘴出口的火焰熱釋放率進(jìn)行處理,提取激振頻率上的擾動幅值和相位,并帶入火焰?zhèn)鬟f函數(shù)的定義式(式(7)),即可得到各火焰在不同頻率和不同擾動幅值下的火焰?zhèn)鬟f函數(shù)。

    1.3 聲學(xué)模型

    1.3.1 聲學(xué)模擬模型

    對整個燃燒室建立Helmholtz方程并耦合火焰?zhèn)鬟f函數(shù),采用噴管聲學(xué)邊界條件,即可求得整個系統(tǒng)的聲學(xué)模態(tài)。

    對于一個有熱源的區(qū)域,流體的密度隨壓力和 熵的變化而變化,即:

    結(jié)合火焰?zhèn)鬟f函數(shù)的定義式(7),最終可推導(dǎo)得到單極域源的表達(dá)式為:

    式中,下標(biāo)b和u分別代表已燃?xì)怏w和未燃?xì)怏w,S為噴嘴截面積,V為熱源體積。

    1.3.2 聲學(xué)邊界阻抗模型

    在實際工況中,聲波在出口處存在一定的反射,即為聲阻抗R。如圖5所示,設(shè)入射波強(qiáng)度為A,反射波強(qiáng)度為A。

    圖5 聲波反射與聲學(xué)邊界阻抗示意圖Fig.5 Schematic diagram of sound reflection and acoustic boundary impedance

    當(dāng)集氣環(huán)與噴注器和燃燒室解耦后,噴注器進(jìn)口可認(rèn)為是壓力節(jié)點(diǎn),空氣燃料混合噴嘴入口可看作是硬聲場邊界,空氣噴嘴入口為軟聲場邊界。燃燒室出口與噴管相連,其聲學(xué)邊界條件采用1.1小節(jié)中推導(dǎo)得到的聲學(xué)反射系數(shù)模型。

    聲阻抗定義為:

    2 燃燒室不穩(wěn)定模態(tài)估計

    由于空氣加熱器燃燒室的噴嘴面積總和與燃燒室橫截面積之比遠(yuǎn)遠(yuǎn)小于1,且噴嘴內(nèi)與燃燒室的溫度相差較大(約1 500 K),因此噴嘴前的聲學(xué)系統(tǒng)可與燃燒室內(nèi)的聲學(xué)系統(tǒng)解耦。燃燒室出口為一壅塞的噴管,燃燒室內(nèi)的聲學(xué)系統(tǒng)也可以與噴管下游的聲學(xué)系統(tǒng)解耦,故只需計算燃燒室內(nèi)的聲學(xué)模態(tài)。在粗估燃燒室聲學(xué)模態(tài)的頻率時,燃燒室進(jìn)口可認(rèn)為是聲學(xué)閉口條件,出口由于噴管前馬赫數(shù)小于0.1,也可近似認(rèn)為是聲學(xué)閉口。

    燃燒室的各個聲學(xué)模態(tài)的頻率f可以根據(jù)經(jīng)典公式進(jìn)行粗估:

    式中,、、1、、分別為縱向模態(tài)數(shù)、切向模態(tài)數(shù)、徑向模態(tài)數(shù)、燃燒室聲速、徑向波數(shù)、縱向波數(shù)、垂直于縱向方向的波數(shù)和燃燒室的直徑。為 貝塞爾函數(shù) dJ(k)/d=0的根,可通過查表2得到。

    表2 不同模態(tài)下 k m,j 的值Table 2 The value of k m,j in different modes

    結(jié)合加熱器燃燒室結(jié)構(gòu)參數(shù)以及燃燒室內(nèi)部聲速(約為837 m/s),可以粗估燃燒室的各聲學(xué)不穩(wěn)定模態(tài)頻率,如表3所示。

    表3 燃燒室各不穩(wěn)定模態(tài)頻率估計值Table 3 Estimated frequency of each unstable mode of the combustion chamber

    考慮到實際加熱器燃燒室最易出現(xiàn)(n=1,m=1)的不穩(wěn)定模態(tài),因此最可能出現(xiàn)的不穩(wěn)定頻率應(yīng)為(n=1,m=1,j=1)模態(tài)的頻率1 406 Hz。這一粗估頻率將作為后續(xù)計算的初始參考值。

    3 結(jié)果與討論

    3.1 聲學(xué)反射系數(shù)計算結(jié)果

    對1.1小節(jié)所得的數(shù)學(xué)方法進(jìn)行求解,可得0~3 000 Hz頻率范圍內(nèi)噴管聲學(xué)反射系數(shù)的幅頻、相頻特性(圖6藍(lán)色圈)。對計算結(jié)果進(jìn)行五階多項式擬合,可求得噴管聲學(xué)反射系數(shù)的近似解析表達(dá)式,其頻響特性曲線如圖6中紅色實線所示。此結(jié)果可作為燃燒室聲學(xué)計算的出口邊界條件,用于計算燃燒室的各聲學(xué)模態(tài)。

    圖6 聲學(xué)反射系數(shù)的幅頻、相頻特性曲線Fig.6 Amplitude frequency and phase frequency characteristic curve of acoustic reflection coefficient

    3.2 燃燒室內(nèi)數(shù)值計算結(jié)果

    3.2.1 溫度場分析

    由于粗估燃燒室最不穩(wěn)定頻率為1 406 Hz,因此取擾動頻率f =1 400 Hz、擾動振幅=04的工況進(jìn)行分析。圖7為t=0.021 161 s時刻,燃燒室距噴注面板50、100、200、400 mm處截面的總溫云圖。

    圖7 不同截面處的總溫云圖Fig.7 Total temperature contour map at different cross-sections

    從圖7中可以看到,燃料噴嘴出口近場截面的溫度分布并不均勻,中心噴嘴受周圍冷卻氣流的影響更為明顯,存在明顯的低溫環(huán)(圖7中黑色區(qū)域);隨著軸向距離的增加,冷熱空氣逐步混合,在距噴注面板約400 mm處幾乎混合均勻,燃燒室內(nèi)部溫度分布趨于穩(wěn)定。

    3.2.2 熱釋放率分析

    圖8是t=0.021 161 s、f=1 400 Hz、3種不同擾動振幅下,燃燒室中心面上的火焰熱釋放率鋒面與流線圖。其中,熱釋放率鋒面可近似認(rèn)為是火焰鋒面。在施加擾動后,燃燒室的中心與邊緣靠近壁面處會分別產(chǎn)生兩個較為強(qiáng)烈的渦旋。處于中心位置的渦旋對火焰鋒面的拉伸作用明顯。在擾動振幅相對較小時,渦旋強(qiáng)度較低(流線相對稀疏平整),火焰鋒面平整,形狀較為舒展;隨擾動振幅增加,中心渦旋強(qiáng)度增加(流線更為緊湊卷曲),對火焰鋒面的影響愈發(fā)明顯,火焰產(chǎn)生明顯褶皺,且長度變短。邊緣渦旋對于流量脈動幅值的響應(yīng)并不是很明顯。在高頻擾動下,中心渦旋造成的速度擾動直接對火焰產(chǎn)生擾動,造成火焰面周期性卷曲和拉伸,從而產(chǎn)生火焰?zhèn)鞑ニ俣葦_動,引起火焰熱釋放率的波動。

    本文分別檢測了每一個噴嘴的熱釋放率隨時間的變化,以圖8中標(biāo)注的噴嘴為例,特別對受到冷空氣擠壓的火焰進(jìn)行討論。圖9(a)、(b)分別為該噴嘴在擾動振幅=04、不同擾動頻率f下和擾動頻率f1 400 Hz、不同振幅下,火焰熱釋放率隨時間變化的曲線。

    圖8 不同振幅下的燃燒室內(nèi)熱釋放率云圖及流線圖Fig.8 Contour of heat release rate and streamline inside the combustion chamber under different amplitudes

    從圖9(a)可以看到,當(dāng)振幅相同時,該噴嘴出口火焰熱釋放率對頻率的響應(yīng)不同—不同擾動頻率對應(yīng)不同的熱釋放率脈動幅值。理論上,當(dāng)對入口流量給定相同的擾動振幅時,擾動頻率越高,熱釋放脈動幅值應(yīng)更低。但對于該噴嘴,在擾動頻率f=1 400 Hz時,其幅值相對f=1 250 Hz并未出現(xiàn)明顯的下降,即燃燒室內(nèi)部的不穩(wěn)定頻率可能出現(xiàn)在f=1 400 Hz附近,這一結(jié)果也證明了第2節(jié)的粗估結(jié)果是可靠的。

    從圖9(b)可以看到,當(dāng)擾動頻率f=1 400 Hz時,該噴嘴出口火焰對于不同振幅的響應(yīng)同樣存在差異。隨擾動振幅的增加,熱釋放脈動幅值增加;但熱釋放率脈動幅值并不隨流量擾動幅值的增加而持續(xù)增加。當(dāng)=06 時,熱釋放脈動幅值與=04 時并無明顯差異,這也說明了該噴嘴出口火焰對于大流量脈動的響應(yīng)并不敏感。

    圖9 火焰熱釋放率隨時間變化曲線Fig.9 Dimensionless flame heat release rate vs.time curve

    3.3 燃燒室內(nèi)聲學(xué)模態(tài)預(yù)測結(jié)果

    圖10為示例噴嘴(圖8)對應(yīng)的火焰?zhèn)鬟f函數(shù)幅值和相位圖。從圖中可以看到,在任意擾動頻率和擾動振幅下,該噴嘴的火焰?zhèn)鬟f函數(shù)幅值均小于1,即火焰熱釋放率相對于入口流量擾動的響應(yīng)呈抑制作用。并且隨著擾動流量的增加,傳遞函數(shù)幅值在各頻率下均有明顯的下降趨勢。

    圖10 火焰?zhèn)鬟f函數(shù)的幅值和相位Fig.10 The amplitude and phase of the flame transfer function

    將各噴嘴的火焰?zhèn)鬟f函數(shù)結(jié)果代入1.3小節(jié)的聲學(xué)模型中,即可對存在流量擾動的燃燒室的穩(wěn)定性進(jìn)行預(yù)測,預(yù)測結(jié)果如圖11所示。該圖為燃燒室無量綱(n=1,m=1,j=1)模態(tài)的壓力分布。壓力分布模式主要存在2種,分布在沿燃燒室直徑方向上互相垂直。經(jīng)1.3小節(jié)中的模型計算,當(dāng)擾動振幅分別為 02、04、06時,(n=1,m=1,j=1)模態(tài)的特征值分別為1 389.9+5.975 i、1 389.9+6.060 3 i和1 389.9+6.142 4 i。該特征值為復(fù)數(shù),其實部為模態(tài)頻率,虛部為增長率。當(dāng)虛部大于0時,代表燃燒室在該聲學(xué)模態(tài)上是穩(wěn)定的;當(dāng)虛部小于0時,代表燃燒室在該聲學(xué)模態(tài)上是不穩(wěn)定的。

    圖11 不同擾動振幅下,燃燒室最不穩(wěn)定模態(tài)的總聲壓場Fig.11 The total sound pressure field of the most unstable mode in the combustion chamber

    結(jié)果表明,在不同擾動幅值下,燃燒室的(n=1,m=1,j=1)模態(tài)均穩(wěn)定。這與該型加熱器的試驗結(jié)果一致。該型加熱器在后期試驗中并未出現(xiàn)燃燒不穩(wěn)定現(xiàn)象。另一型加熱器(燃燒室直徑為460 mm,長度為800 mm)在試驗過程中出現(xiàn)燃燒不穩(wěn)定現(xiàn)象,振蕩頻率為1 100 Hz。基于第2節(jié)中的方法對該燃燒室的不穩(wěn)定模態(tài)進(jìn)行預(yù)測,得到振蕩頻率為1 135 Hz,與試驗測量頻率符合很好,進(jìn)一步驗證了本文方法良好的預(yù)測能力。

    4 結(jié) 論

    1)采用Magnus展開法獲得了噴管的聲學(xué)邊界條件。通過五階多項式擬合獲得了噴管聲反射系數(shù)的近似解析表達(dá)式。

    2)采用ANSYS Fluent對全尺寸燃燒器進(jìn)行了非穩(wěn)態(tài)數(shù)值仿真,得到了燃燒室內(nèi)部存在擾動時的物理場信息,發(fā)現(xiàn)位于噴注面板中心處的噴嘴受到周圍冷卻空氣的影響較大。提取了每一個噴嘴的熱釋放率信息,獲得了噴嘴熱釋放率脈動與速度擾動之間的火焰?zhèn)鬟f函數(shù)。

    3)將所得聲學(xué)邊界條件與火焰?zhèn)鬟f函數(shù)代入燃燒室聲學(xué)模型中,成功預(yù)測該燃燒室最易發(fā)生不穩(wěn)定燃燒的頻率為1 389.9 Hz,且不同擾動振幅下,燃燒室的(n=1,m=1,j=1)模態(tài)均穩(wěn)定,不會發(fā)生不穩(wěn)定燃燒現(xiàn)象,與試驗結(jié)果一致。

    a级片在线免费高清观看视频| 青春草视频在线免费观看| 777久久人妻少妇嫩草av网站| 国产日韩一区二区三区精品不卡| 亚洲精品国产一区二区精华液| 青草久久国产| 欧美精品一区二区大全| 午夜免费观看性视频| av女优亚洲男人天堂| 观看美女的网站| 操出白浆在线播放| 1024视频免费在线观看| 精品人妻熟女毛片av久久网站| 黄网站色视频无遮挡免费观看| 一级毛片电影观看| 麻豆av在线久日| 捣出白浆h1v1| 天天影视国产精品| 国产精品国产三级国产专区5o| 99热网站在线观看| 亚洲第一区二区三区不卡| 欧美黑人精品巨大| 久久久精品区二区三区| 亚洲精品国产av成人精品| 日本一区二区免费在线视频| 肉色欧美久久久久久久蜜桃| 黄片小视频在线播放| 天堂俺去俺来也www色官网| 人人澡人人妻人| 最新在线观看一区二区三区 | 成人国语在线视频| 天堂中文最新版在线下载| 尾随美女入室| 午夜精品国产一区二区电影| 国产免费视频播放在线视频| 爱豆传媒免费全集在线观看| 亚洲人成电影观看| 天天躁夜夜躁狠狠久久av| 久久鲁丝午夜福利片| 日日啪夜夜爽| 亚洲熟女毛片儿| 午夜福利视频精品| 久久久久视频综合| 哪个播放器可以免费观看大片| 久久精品熟女亚洲av麻豆精品| 我要看黄色一级片免费的| 亚洲精品乱久久久久久| 91aial.com中文字幕在线观看| 欧美日韩综合久久久久久| 精品人妻熟女毛片av久久网站| 一区福利在线观看| 国产亚洲av片在线观看秒播厂| 亚洲五月色婷婷综合| 国产精品久久久久成人av| av电影中文网址| 热re99久久国产66热| 国产1区2区3区精品| 中文字幕最新亚洲高清| 日本猛色少妇xxxxx猛交久久| 香蕉国产在线看| 男女边摸边吃奶| 亚洲欧美精品自产自拍| 色网站视频免费| 国产在线免费精品| 国产又爽黄色视频| 亚洲精品一区蜜桃| 人人妻,人人澡人人爽秒播 | 99香蕉大伊视频| 国产一区二区三区av在线| 久久天躁狠狠躁夜夜2o2o | 欧美日本中文国产一区发布| 久久久久久人妻| 国产无遮挡羞羞视频在线观看| 啦啦啦中文免费视频观看日本| 国产欧美日韩综合在线一区二区| 欧美激情高清一区二区三区 | 亚洲国产av影院在线观看| 99香蕉大伊视频| 国产黄频视频在线观看| 在线看a的网站| 国产男女超爽视频在线观看| 中文字幕人妻丝袜一区二区 | 晚上一个人看的免费电影| 国产精品99久久99久久久不卡 | 老司机影院成人| 在线免费观看不下载黄p国产| 色婷婷久久久亚洲欧美| 又黄又粗又硬又大视频| 99热全是精品| 老司机深夜福利视频在线观看 | 人妻一区二区av| 黄片小视频在线播放| 熟妇人妻不卡中文字幕| 精品少妇内射三级| 卡戴珊不雅视频在线播放| 美女福利国产在线| 午夜福利视频精品| 国产99久久九九免费精品| 2018国产大陆天天弄谢| 高清不卡的av网站| 日韩欧美精品免费久久| 亚洲熟女毛片儿| 日本欧美视频一区| 精品国产乱码久久久久久男人| 91aial.com中文字幕在线观看| 亚洲av成人精品一二三区| 成人国产麻豆网| 别揉我奶头~嗯~啊~动态视频 | 看免费av毛片| 男男h啪啪无遮挡| 91精品伊人久久大香线蕉| 毛片一级片免费看久久久久| 精品少妇内射三级| 一边摸一边抽搐一进一出视频| 中国国产av一级| www日本在线高清视频| avwww免费| 亚洲国产欧美一区二区综合| 欧美日韩视频高清一区二区三区二| 午夜激情av网站| 宅男免费午夜| 国产一区二区在线观看av| 9191精品国产免费久久| 亚洲精品久久成人aⅴ小说| 午夜福利网站1000一区二区三区| 亚洲av福利一区| 夜夜骑夜夜射夜夜干| 91成人精品电影| 最近中文字幕2019免费版| 日日撸夜夜添| 99热全是精品| 欧美人与善性xxx| 欧美亚洲 丝袜 人妻 在线| 亚洲情色 制服丝袜| 久久久久久久国产电影| 男女之事视频高清在线观看 | 久久久精品国产亚洲av高清涩受| 久久精品久久精品一区二区三区| 精品国产超薄肉色丝袜足j| 精品少妇黑人巨大在线播放| 亚洲国产av影院在线观看| 女人爽到高潮嗷嗷叫在线视频| 自线自在国产av| 色网站视频免费| av免费观看日本| 国产在线视频一区二区| 一级黄片播放器| www日本在线高清视频| 亚洲精品在线美女| 乱人伦中国视频| h视频一区二区三区| 国产日韩欧美在线精品| 国产 一区精品| av免费观看日本| 亚洲伊人久久精品综合| 中文字幕色久视频| 大片免费播放器 马上看| 妹子高潮喷水视频| 中文字幕精品免费在线观看视频| 久久精品国产综合久久久| 在线观看www视频免费| 免费少妇av软件| 国产一区二区三区综合在线观看| 精品酒店卫生间| 国产av码专区亚洲av| 51午夜福利影视在线观看| 亚洲美女黄色视频免费看| 亚洲国产欧美在线一区| 99久久99久久久精品蜜桃| 婷婷成人精品国产| 亚洲精品久久成人aⅴ小说| 精品免费久久久久久久清纯 | 狠狠精品人妻久久久久久综合| 极品少妇高潮喷水抽搐| 久久久久久久久久久免费av| 国产欧美日韩综合在线一区二区| 国产高清国产精品国产三级| a 毛片基地| 国产一卡二卡三卡精品 | 午夜激情久久久久久久| 免费观看性生交大片5| 美国免费a级毛片| 欧美日韩亚洲高清精品| 免费少妇av软件| 久久精品aⅴ一区二区三区四区| 亚洲在久久综合| 男女免费视频国产| 极品人妻少妇av视频| 日本av手机在线免费观看| 99热全是精品| 99re6热这里在线精品视频| 免费不卡黄色视频| 少妇 在线观看| 亚洲,欧美精品.| 操出白浆在线播放| 一级a爱视频在线免费观看| 国产不卡av网站在线观看| a级毛片在线看网站| 成年人免费黄色播放视频| 一边摸一边做爽爽视频免费| kizo精华| 成年人免费黄色播放视频| 亚洲人成电影观看| 欧美 日韩 精品 国产| 热99久久久久精品小说推荐| 国产野战对白在线观看| 亚洲精品一区蜜桃| 日韩精品免费视频一区二区三区| 在线观看一区二区三区激情| 夫妻性生交免费视频一级片| av在线观看视频网站免费| 国产精品三级大全| 久久久久网色| 天天躁狠狠躁夜夜躁狠狠躁| 伦理电影大哥的女人| 免费在线观看黄色视频的| av国产精品久久久久影院| 热re99久久精品国产66热6| 久久久欧美国产精品| 天天操日日干夜夜撸| 国产男人的电影天堂91| 亚洲欧美中文字幕日韩二区| 日日啪夜夜爽| 亚洲精品自拍成人| 毛片一级片免费看久久久久| 午夜福利,免费看| 亚洲伊人久久精品综合| 中文字幕制服av| 黄频高清免费视频| 日本爱情动作片www.在线观看| 日韩制服丝袜自拍偷拍| 在线观看免费日韩欧美大片| 在现免费观看毛片| 女人久久www免费人成看片| 啦啦啦在线观看免费高清www| 亚洲欧美激情在线| 欧美激情高清一区二区三区 | 伊人久久大香线蕉亚洲五| av网站免费在线观看视频| 日日摸夜夜添夜夜爱| 一区福利在线观看| 国产国语露脸激情在线看| 97人妻天天添夜夜摸| 国产精品.久久久| 欧美中文综合在线视频| 丁香六月天网| 最黄视频免费看| 午夜福利影视在线免费观看| 国产成人精品福利久久| 成人影院久久| 亚洲视频免费观看视频| 男女床上黄色一级片免费看| 亚洲美女黄色视频免费看| 色吧在线观看| 色吧在线观看| 最新在线观看一区二区三区 | 国产亚洲一区二区精品| 国产成人系列免费观看| 亚洲av日韩在线播放| 国产男女内射视频| 国产 精品1| 亚洲精品久久成人aⅴ小说| 国产精品 国内视频| 免费观看av网站的网址| 日韩熟女老妇一区二区性免费视频| 街头女战士在线观看网站| 9191精品国产免费久久| 亚洲国产日韩一区二区| 国产片内射在线| 亚洲国产毛片av蜜桃av| 国产一区二区在线观看av| 日韩中文字幕欧美一区二区 | 99热国产这里只有精品6| 最近的中文字幕免费完整| 国产精品av久久久久免费| 天堂8中文在线网| 欧美97在线视频| 日本一区二区免费在线视频| 亚洲av综合色区一区| 老鸭窝网址在线观看| 久久精品国产综合久久久| 国产精品二区激情视频| 中文字幕最新亚洲高清| 赤兔流量卡办理| 菩萨蛮人人尽说江南好唐韦庄| 在线 av 中文字幕| 最近中文字幕2019免费版| 最近中文字幕高清免费大全6| 搡老岳熟女国产| 免费在线观看黄色视频的| 亚洲欧美日韩另类电影网站| 亚洲av日韩在线播放| 国产免费视频播放在线视频| 亚洲在久久综合| av一本久久久久| 丝袜喷水一区| 精品一区二区三区四区五区乱码 | 一级a爱视频在线免费观看| 99精国产麻豆久久婷婷| 男男h啪啪无遮挡| 亚洲男人天堂网一区| 日韩av免费高清视频| 在线观看人妻少妇| 性少妇av在线| 亚洲av成人不卡在线观看播放网 | 在线天堂最新版资源| 久久99热这里只频精品6学生| 日本欧美视频一区| 一本色道久久久久久精品综合| 97精品久久久久久久久久精品| 捣出白浆h1v1| 看免费av毛片| 99国产精品免费福利视频| 免费久久久久久久精品成人欧美视频| 国产乱人偷精品视频| 亚洲色图 男人天堂 中文字幕| 我要看黄色一级片免费的| 国产在线免费精品| av国产精品久久久久影院| 亚洲精华国产精华液的使用体验| 性少妇av在线| 乱人伦中国视频| 欧美日韩一级在线毛片| 超碰97精品在线观看| 色精品久久人妻99蜜桃| 日韩 亚洲 欧美在线| 国产高清不卡午夜福利| 国产1区2区3区精品| 99国产精品免费福利视频| 激情五月婷婷亚洲| 国产淫语在线视频| 国产成人啪精品午夜网站| 少妇精品久久久久久久| 亚洲av成人精品一二三区| 如何舔出高潮| 美女主播在线视频| 晚上一个人看的免费电影| 黄片无遮挡物在线观看| 精品国产超薄肉色丝袜足j| 电影成人av| 国产片特级美女逼逼视频| 国产熟女欧美一区二区| 国语对白做爰xxxⅹ性视频网站| 国产精品国产av在线观看| 国产男人的电影天堂91| 亚洲成人免费av在线播放| 少妇人妻 视频| 一边摸一边做爽爽视频免费| 亚洲成色77777| 亚洲免费av在线视频| 国产黄色视频一区二区在线观看| 国产又爽黄色视频| 亚洲精品,欧美精品| 美女福利国产在线| 久久婷婷青草| 日韩精品有码人妻一区| 2021少妇久久久久久久久久久| 亚洲国产毛片av蜜桃av| 国产一区二区三区av在线| 亚洲成人手机| 午夜影院在线不卡| 日韩av不卡免费在线播放| 男女国产视频网站| 飞空精品影院首页| 免费女性裸体啪啪无遮挡网站| 国产极品粉嫩免费观看在线| 免费在线观看视频国产中文字幕亚洲 | 2021少妇久久久久久久久久久| 黄网站色视频无遮挡免费观看| 亚洲男人天堂网一区| 日本一区二区免费在线视频| 久久 成人 亚洲| 91成人精品电影| 电影成人av| 成人影院久久| 国产精品蜜桃在线观看| 午夜福利免费观看在线| 爱豆传媒免费全集在线观看| 欧美精品亚洲一区二区| 久久毛片免费看一区二区三区| 老鸭窝网址在线观看| 精品第一国产精品| 悠悠久久av| 丝袜美腿诱惑在线| 日本猛色少妇xxxxx猛交久久| 国产精品久久久久久精品电影小说| 一区福利在线观看| 亚洲av中文av极速乱| 最近中文字幕高清免费大全6| 飞空精品影院首页| 美女脱内裤让男人舔精品视频| 久久性视频一级片| 一二三四中文在线观看免费高清| 亚洲成人一二三区av| 黑丝袜美女国产一区| 自拍欧美九色日韩亚洲蝌蚪91| 不卡av一区二区三区| 18禁国产床啪视频网站| 日本一区二区免费在线视频| av电影中文网址| 亚洲精品久久成人aⅴ小说| 国产一区二区三区av在线| av卡一久久| 久久婷婷青草| av免费观看日本| 老司机亚洲免费影院| 1024香蕉在线观看| 操出白浆在线播放| 五月天丁香电影| 亚洲国产精品999| 久久久久精品久久久久真实原创| 午夜福利视频精品| 欧美最新免费一区二区三区| 久久天躁狠狠躁夜夜2o2o | 国产黄频视频在线观看| 高清不卡的av网站| 久久国产精品男人的天堂亚洲| 成人毛片60女人毛片免费| 91老司机精品| 精品亚洲成a人片在线观看| 欧美精品一区二区大全| 最近中文字幕高清免费大全6| 男人舔女人的私密视频| 国产福利在线免费观看视频| 人人妻人人添人人爽欧美一区卜| 伊人久久大香线蕉亚洲五| 国产精品久久久久久人妻精品电影 | 欧美日韩一区二区视频在线观看视频在线| 日韩大片免费观看网站| 精品少妇一区二区三区视频日本电影 | 十八禁高潮呻吟视频| a级毛片在线看网站| 男男h啪啪无遮挡| 午夜福利乱码中文字幕| 中文字幕人妻熟女乱码| 天天躁日日躁夜夜躁夜夜| 午夜免费观看性视频| 免费观看a级毛片全部| 夜夜骑夜夜射夜夜干| 另类亚洲欧美激情| 精品一区二区免费观看| 在线免费观看不下载黄p国产| 国产不卡av网站在线观看| 国产极品粉嫩免费观看在线| 免费日韩欧美在线观看| 深夜精品福利| 国产成人91sexporn| 九色亚洲精品在线播放| 女性生殖器流出的白浆| 多毛熟女@视频| 男女午夜视频在线观看| 十八禁人妻一区二区| 老司机深夜福利视频在线观看 | 欧美精品一区二区大全| 各种免费的搞黄视频| 国产乱来视频区| 免费女性裸体啪啪无遮挡网站| www.精华液| 乱人伦中国视频| 亚洲人成电影观看| 不卡av一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 日韩精品有码人妻一区| 亚洲四区av| 国产xxxxx性猛交| 国产精品免费大片| 久久精品久久精品一区二区三区| 一级黄片播放器| 亚洲精品aⅴ在线观看| 中文字幕av电影在线播放| 久久久欧美国产精品| 一本—道久久a久久精品蜜桃钙片| 国产 精品1| 国产男女超爽视频在线观看| 国产日韩一区二区三区精品不卡| 欧美 日韩 精品 国产| 永久免费av网站大全| 在线观看免费高清a一片| 免费看不卡的av| 99久久精品国产亚洲精品| 日韩中文字幕欧美一区二区 | 国产精品国产av在线观看| 成人国产麻豆网| 亚洲五月色婷婷综合| 免费观看av网站的网址| 午夜免费鲁丝| 两个人免费观看高清视频| av片东京热男人的天堂| 亚洲精品国产区一区二| 亚洲国产欧美一区二区综合| 天天躁日日躁夜夜躁夜夜| 中文字幕人妻丝袜一区二区 | 热re99久久精品国产66热6| 肉色欧美久久久久久久蜜桃| 啦啦啦中文免费视频观看日本| 一个人免费看片子| bbb黄色大片| 啦啦啦啦在线视频资源| 久久人妻熟女aⅴ| 热re99久久国产66热| 高清视频免费观看一区二区| 国产精品免费大片| 97在线人人人人妻| 只有这里有精品99| 一本大道久久a久久精品| 免费在线观看完整版高清| 亚洲一级一片aⅴ在线观看| 最近2019中文字幕mv第一页| 精品人妻熟女毛片av久久网站| 街头女战士在线观看网站| 欧美人与性动交α欧美软件| 亚洲人成网站在线观看播放| 最近最新中文字幕免费大全7| 激情视频va一区二区三区| 黑丝袜美女国产一区| av在线播放精品| 成年动漫av网址| 国产成人免费观看mmmm| 美女视频免费永久观看网站| 母亲3免费完整高清在线观看| xxxhd国产人妻xxx| 一区二区日韩欧美中文字幕| 各种免费的搞黄视频| 欧美亚洲日本最大视频资源| 交换朋友夫妻互换小说| 一级毛片我不卡| 岛国毛片在线播放| av天堂久久9| h视频一区二区三区| 少妇人妻 视频| 妹子高潮喷水视频| 一级,二级,三级黄色视频| 欧美 亚洲 国产 日韩一| 国产精品.久久久| xxxhd国产人妻xxx| 精品人妻在线不人妻| 亚洲免费av在线视频| 国产亚洲最大av| 美女午夜性视频免费| 极品人妻少妇av视频| 久久人人97超碰香蕉20202| 国产女主播在线喷水免费视频网站| 成人漫画全彩无遮挡| 搡老岳熟女国产| 精品一区二区免费观看| 女人被躁到高潮嗷嗷叫费观| 丝袜人妻中文字幕| 免费日韩欧美在线观看| 波野结衣二区三区在线| 看免费av毛片| 亚洲,欧美精品.| 男女边吃奶边做爰视频| 亚洲国产毛片av蜜桃av| 99久久精品国产亚洲精品| 亚洲精品国产av蜜桃| 观看av在线不卡| 亚洲五月色婷婷综合| 亚洲精品aⅴ在线观看| 精品久久久精品久久久| 亚洲中文av在线| 天美传媒精品一区二区| 国产淫语在线视频| 一本色道久久久久久精品综合| 久久精品久久久久久久性| 欧美在线一区亚洲| 不卡av一区二区三区| 久久久国产精品麻豆| 国精品久久久久久国模美| 人人妻人人爽人人添夜夜欢视频| 中文字幕色久视频| 国产亚洲av片在线观看秒播厂| 成人黄色视频免费在线看| 老司机在亚洲福利影院| 亚洲欧美精品自产自拍| 一区二区三区乱码不卡18| 成人亚洲欧美一区二区av| 少妇猛男粗大的猛烈进出视频| 一区二区日韩欧美中文字幕| 精品午夜福利在线看| 精品亚洲成国产av| 日本一区二区免费在线视频| 啦啦啦啦在线视频资源| 18在线观看网站| 国产黄色视频一区二区在线观看| 性高湖久久久久久久久免费观看| 色综合欧美亚洲国产小说| tube8黄色片| 操美女的视频在线观看| 极品少妇高潮喷水抽搐| 国产熟女午夜一区二区三区| 亚洲av综合色区一区| 黄片无遮挡物在线观看| 欧美精品亚洲一区二区| 亚洲av中文av极速乱| 丰满饥渴人妻一区二区三| 91aial.com中文字幕在线观看| 国产片特级美女逼逼视频| 国产高清不卡午夜福利| av国产精品久久久久影院| 19禁男女啪啪无遮挡网站| 日韩成人av中文字幕在线观看| 少妇精品久久久久久久| 99热网站在线观看| 欧美亚洲日本最大视频资源| 麻豆av在线久日| a级毛片黄视频| 久久精品国产亚洲av高清一级| 丁香六月欧美| 在线天堂中文资源库| 成人18禁高潮啪啪吃奶动态图| 国产一区有黄有色的免费视频| 国产 一区精品| 久久久久久免费高清国产稀缺| 视频在线观看一区二区三区| 超碰成人久久|