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

    沿駐點線化學非平衡程度評估方法

    2023-03-09 11:38:02楊天鵬嚴東升于江鵬劉文伶
    導彈與航天運載技術(shù) 2023年1期
    關(guān)鍵詞:化學特征模型

    楊天鵬,嚴東升,于江鵬,郭 陽,劉文伶

    沿駐點線化學非平衡程度評估方法

    楊天鵬,嚴東升,于江鵬,郭 陽,劉文伶

    (北京航天長征飛行器研究所,北京,100076)

    為了研究高溫真實氣體效應中的化學非平衡的程度,針對球頭的沿駐點線流動,通過高溫高壓氣體的化學反應歷程估算化學反應特征時間,使用正激波流動關(guān)系來估算流動特征時間,建立了化學非平衡評估模型,給出了化學非平衡程度隨馬赫數(shù)和高度的變化關(guān)系。計算結(jié)果表明,馬赫數(shù)越大、高度越低,化學反應特征時間比流動特征時間越小,定義的化學反應特征時間與流動特征時間量級相當時為化學非平衡流動,時間比大于10則接近化學凍結(jié)流,時間比小于0.1則接近化學平衡流。在設定的球頭模型下,化學反應平衡流在高度0km、10km、20km、30km、40km分別需要馬赫數(shù)大于8.0、9.5、10.8、12.7、20.0,化學反應凍結(jié)流在高度30km、40km、50km、60km分別需要馬赫數(shù)小于8.4、10.2、12.9、17.5,中間來流狀態(tài)對應化學反應非平衡狀態(tài)。

    數(shù)值模擬;高溫真實氣體效應;化學非平衡

    0 引 言

    飛行器在大氣中高速運動時,其前緣附近產(chǎn)生的強激波會使得激波層內(nèi)的空氣處于高溫高壓狀態(tài),此時空氣會發(fā)生一系列復雜的化學反應和能量變化,這種現(xiàn)象被稱為高溫真實氣體效應。對于高溫真實氣體繞流,一般可將流動分為3種情況[1]:第1種是凍結(jié)流,此時流場中的化學反應速率為零,氣體組元的組成不隨空間和時間變化,但粘性流動會使得氣體組元由于擴散而改變;第2種是平衡流,此時流場中的氣體處于熱力學平衡狀態(tài),且流體處于化學平衡狀態(tài),此時流場中的化學反應速率為無窮大;第3種是非平衡流,流場中氣體的熱力學性質(zhì)和化學反應狀態(tài)隨時間空間發(fā)生變化,此時化學反應或能量變化的特征時間與流動特征時間處于相近量級。

    對于高溫真實氣體效應中涉及到的化學平衡態(tài)問題,在多個文獻中都給出了大致的參考[2,3],但是這些文獻并沒有對平衡態(tài)問題進行定量分析。文獻[4]認為可以通過達姆科勒數(shù)am,用于描述同一系統(tǒng)中化學反應相比其他現(xiàn)象的相對時間尺度,使用粒子碰撞特征時間與流場特征時間尺度之比來定義流動的化學狀態(tài)。但是該定義是一種寬泛的定性認識,一般認為特征時間之比高于4個量級后才能被劃分到另一流態(tài),不便于工程應用。由于化學反應與組元構(gòu)成、能量輸運、化學反應機理等等密切相關(guān),并且流動中還需要關(guān)注湍流、邊界層、壁面催化、燒蝕等一系列與化學反應相耦合的問題,其表現(xiàn)出強烈的非線性,建模難度非常高,難以給出明確的公式對其進行計算。

    為了給出高溫真實氣體效應中對化學非平衡程度的評估,本文從化學反應的宏觀表現(xiàn)出發(fā),使用高溫高壓空氣的化學反應歷程時間與流動特征時間作為比較,建立化學非平衡評估模型,給出化學非平衡程度隨高度和馬赫數(shù)的變化規(guī)律,為流場化學反應模型選取、非平衡流場計算結(jié)果不確定度等研究提供參考,提升高速飛行器氣動特性預示能力。

    1 沿駐點線化學反應流動特征分析

    為了定量地建立空氣化學反應非平衡程度評估模型,需要首先對相應流場的物理特征進行分析。本文選取化學非平衡程度最為嚴重的駐點線區(qū)域作為評估對象,使用數(shù)值模擬方法對化學反應效應和流動特征進行分析。

    本文使用多組元可壓縮雷諾平均Navier-Stokes方程作為控制方程,其積分形式[5]為

    計算模型選用半徑50 mm的球頭,網(wǎng)格量為100×100,計算網(wǎng)格如圖1所示。邊界設置上來流為遠場邊界,物面使用無滑移、絕熱壁和非催化壁邊界條件,排除物面模型對流動非平衡程度的影響。來流組元為空氣,使用Gupta給出的5組元六反應模型及其化學反應速率系數(shù)[6],化學反應模型如表1所示。

    圖1 計算網(wǎng)格

    表1 化學反應方程

    Tab.1 Chemical Reaction Equation

    序號化學反應方程式 1O2+M=2O+M 2N2+M=2N+M 3NO+M=N+O+M 4N2+O=NO+N 5NO+O=O2+N 6N2+N=2N+N

    下面給出高度30 km條件下,分別為10、15、20的3種情況下沿駐點線化學反應能量釋放情況。這里定義反應的化學反應功率密度E

    式中 fi為化學反應i的凈反應速率;ΔHi為化學反應i的焓變。通過該表達式可以得知在空間中某化學反應所造成的產(chǎn)熱功率,其反映了某化學反應對流場內(nèi)能的影響程度,正值為吸熱,負值為放熱。從圖2中可以看出,不同來流馬赫數(shù)下,化學反應流場表現(xiàn)出不同的能量釋放模式,在Ma=10的情況下,化學反應吸熱主要以O2分解為主;在Ma=15的情況下,開始出現(xiàn)NO分解;在Ma=20的情況下,NO分解與O2分解并重,且N2分解反應也表現(xiàn)出較高的吸熱量。這說明在評估化學非平衡時,需要同時將O2、N2和NO 3個組元選為評估指標。

    續(xù)圖2

    圖3給出了=15情況下沿駐點線壓強、溫度和速度分布??梢钥吹剑伛v點線的壓強在激波后略有升高,基本上保持不變。激波后溫度由于激波層較窄,首先呈現(xiàn)出化學凍結(jié)的高溫,之后由于空氣分解吸熱導致溫度逐漸降低,最后形成一個變化較小的溫度平臺。而激波后的速度則呈現(xiàn)出接近直線的分布。激波后流動參數(shù)分布的規(guī)律性使得后文中建立的化學反應非平衡評估模型更加接近于物理實際情況。

    圖3 Ma=15算例沿駐點線流動參數(shù)分布

    續(xù)圖3

    2 化學反應特征時間評估

    根據(jù)上文分析得到的沿駐點線流動特征,本文使用密閉容器中的化學反應進程來評估化學反應特征時間。容器內(nèi)的起始溫度和壓強取來流總溫和總壓,隨著反應進程的進行,容器內(nèi)的溫度逐漸降低,該模型與圖3給出的沿駐點線上流動參數(shù)變化過程較為接近。

    圖4給出了取=15、高度30 km條件下的組元變化歷程,給出了以O2、N2和NO為參考,組元濃度分別達到80%、90%和99%平衡狀態(tài)情況下的時刻位置??梢钥吹綄⒒瘜W反應特征時間定義為容器中的氧氣和氮氣達到80%平衡狀態(tài)時的反應時間較為合適,之后化學反應處于反應緩慢的平臺區(qū),化學反應導致的組元變化和能量變化主要集中于前80%部分,若選擇更高的非平衡標準會使得化學反應特征時間成倍增加,化學非平衡評估模型的敏感性降低。

    圖4 Ma=15算例沿駐點線組元變化歷程

    圖5給出了化學反應特征時間隨馬赫數(shù)和高度的變化關(guān)系,可以看到,馬赫數(shù)越大、高度越低,化學反應特征時間越小,這主要是由于壓強和溫度增加有助于分子碰撞,從而有利于化學反應進行。但是壓強增加并不總有利于空氣化學反應進行,可以看到=12附近的曲線會出現(xiàn)轉(zhuǎn)折,在轉(zhuǎn)折兩側(cè)的曲線斜率發(fā)生了明顯變化。這是由于壓強增加會增大空氣化學反應中各組元的濃度,使得分解反應難以進行,出現(xiàn)局部溫度壓強增大并不促進化學反應發(fā)生的現(xiàn)象。

    圖5 化學反應特征時間

    3 流動特征時間評估

    本文以流體微團沿駐點線從正激波后運動到物面的時間為流動特征時間。首先對激波脫體距離進行估算,球頭激波脫體距離近似公式[7]如下:

    圖6給出了不同高度下流動特征時間隨馬赫數(shù)變化情況,由圖6可知,在本文所使用的模型中,馬赫數(shù)對流動特征時間起主導作用,不同高度的流動特征時間變化不大,其整體上表現(xiàn)出隨著馬赫數(shù)增大而減小的趨勢。

    圖6 流動特征時間

    4 化學非平衡程度模型

    本文將化學非平衡定義為化學反應特征時間與流動特征時間處于相近量級,圖7給出了化學反應特征時間與流動特征時間之比隨馬赫數(shù)的變化??梢钥吹?,由于不同高度下的流動特征時間較為接近,因此化學非平衡程度整體變化情況與化學特征時間模型基本一致。對于不同高度,高度越高,流動越接近于化學凍結(jié)流,這是由于低壓使得化學反應更難發(fā)生,化學反應時間更長;而對于低高度,流動越接近于化學平衡流,這是由于高壓促進分子碰撞,使得化學反應時間更短。

    圖7 化學反應特征時間與流動特征時間之比

    下面進一步對比實際情況與建模之間的不同。為方便起見,選擇高度40 km、馬赫數(shù)為10、13、20的來流參數(shù)進行計算,其化學反應特征時間與流動特征時間之比分別為12.3、1.13和0.094。圖8給出了沿駐點線不同組元摩爾分數(shù)分布的對比,其中0D代表化學反應特征時間的評估結(jié)果,其零時間對應于激波面,右側(cè)終點為達到80%平衡的時間。從圖8中可以看到,對于=10的情況,流動接近于化學凍結(jié)流,組元整體變化程度較小;對于=13的情況,流動處于明顯的化學平衡態(tài),且在接近物面處組元的分布恰好接近于平衡態(tài),可以認為此時化學反應時間與流動時間基本相當;對于=20的情況,可以看到沿駐點線組元很快達到了平衡態(tài),流動接近于化學平衡流。由于本文采用了正激波模型,因此數(shù)值模擬得到的結(jié)果相比于建模結(jié)果,其激波后溫度較低,化學反應速率較低。

    圖8 沿駐點線組元變化歷程與非平衡建模結(jié)果對比

    續(xù)圖8

    以上3個算例對比說明本文所建立的化學非平衡程度評估模型是可行的,可以用于指導化學非平衡數(shù)值模擬。對于流動中化學非平衡的判據(jù)則需要根據(jù)數(shù)值模擬的需要進行,從結(jié)果來看,化學反應特征時間比流動特征時間的值取為0.1~10之間可以認為是化學非平衡,大于10則接近化學平衡流,小于0.1則接近化學凍結(jié)流。需要注意的是,在低馬赫數(shù)下化學非平衡程度呈現(xiàn)出向化學平衡流移動的趨勢,但是此時流場中的溫度不夠高,高溫真實氣體效應不明顯,即使是化學非平衡流也可使用化學凍結(jié)流模型進行計算。

    由于前緣半徑與激波脫體距離近似成正比關(guān)系,前緣半徑越大沿駐點線的流動特征時間越長,但化學反應特征時間是不具有幾何相似性的,因此前緣半徑越大,化學反應特征時間比流動特征時間越小。實際飛行器前緣半徑的取值與本文所選數(shù)值不會有量級上的差異,相比于化學反應特征時間約3~4個量級的跨度,前緣半徑改變帶來的時間變化相對較小。

    5 結(jié) 論

    本文針對球頭繞流中的駐點線流動,通過高溫高壓氣體的化學反應歷程給出了化學反應特征時間,使用正激波流動關(guān)系來估算流動特征時間,建立了化學非平衡評估模型,通過化學反應特征時間與流動特征時間之比給出了流場中的化學非平衡程度隨馬赫數(shù)和高度的變化結(jié)果,并通過數(shù)值模擬結(jié)果驗證了非平衡評估結(jié)果的正確性。可用于指導高速流動中的化學反應模型選取。

    計算結(jié)果表明,化學反應特征時間隨馬赫數(shù)增大而減小,隨高度增加而增加;流動特征時間隨高度變化不大,隨馬赫數(shù)增大而減?。换瘜W反應特征時間與流動特征時間之比與化學反應特征時間規(guī)律相近,其隨馬赫數(shù)增大而減小,隨高度增加而增加。與數(shù)值模擬結(jié)果對比表明,時間比處于0.1~10之間可以認為是化學非平衡,大于10則接近化學凍結(jié)流,小于0.1則接近化學平衡流。在本文設定的球頭模型下,化學反應平衡流在高度0 km、10 km、20 km、30 km、40 km分別需要馬赫數(shù)大于8、9.5、10.8、12.7、20.0,化學反應凍結(jié)流在高度30 km、40 km、50 km、60 km分別需要馬赫數(shù)小于8.4、10.2、12.9、17.5,中間來流狀態(tài)對應化學反應非平衡狀態(tài)。

    [1] 歐陽水吾, 謝中強. 高溫非平衡空氣繞流[M]. 北京: 國防工業(yè)出版社, 2001.

    Ouyang Shuiwu, Xie Zhongqiang. High temperature non-equilibrium flow[M]. Beijing: National Defence Industry Press, 2001.

    [2] Miller J H. Computational aerothermodynamic datasets for hypersonic heat transfer on reentry vehicles[C]. San Francisco: AIAA Atmospheric Flight Mechanics Conference and Exhibit, 2005.

    [3] Tauber M E, Meneses G P. Aerothermodynamics of transatmospheric vehicles[C]. Boston: 4th Thermophysics and Heat Transfer Conference, 1986.

    [4] 中國科學院. 中國學科發(fā)展戰(zhàn)略新型飛行器中的關(guān)鍵力學問題[M]. 北京: 科學出版社, 2018.

    Chinese Academy of Sciences. Chinese discipline development strategy key mechanical problems in new aircraft[M]. Beijing: Science Press, 2018.

    [5] 王江峰, 等. 高超聲速復雜氣動問題數(shù)值方法研究進展[J]. 航空學報, 2015, 36(1): 159-175.

    Wang Jiangfeng, et al. Progress in numerical simulation techniques of hypersonic aerodynamic problems[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(1): 159-175.

    [6] Gupta R N, et al. A review of reaction rates and thermodynamic and transport properties for an 11-species air model for chemical and thermal nonequilibrium calculations to 30000K[R]. NASA-RP-1232, 1989.

    [7] Anderson J D. 高超聲速和高溫氣體動力學(第2版)[M]. 北京: 航空工業(yè)出版社, 2013.

    Anderson J D. Hypersonic and high-temperature gas dynamics(Second Edition)[M]. Beijing: Aviation Industry Press, 2013.

    Evaluation Method of Chemical Non-equilibrium along Stagnation Line

    Yang Tian-peng, Yan Dong-sheng, Yu Jiang-peng, Guo Yang, Liu Wen-ling

    (Beijng Institute of Space Long March Vehicle, Beijing, 100076)

    In order to study the degree of chemical non-equilibrium in high-temperature real gas effect, this paper begins with the flow of the ball head along the stagnation line to establish a chemical non-equilibrium evaluation model. In this model, the chemical reaction characteristic time is calculated by using the chemical reaction process, and the flow characteristic time is evaluated by using the relationship before and after the normal shock wave. The degree of chemical non-equilibrium variation law with Mach number and height is given. Calculation results show that the larger the Mach number and the lower the height makes the ratio of chemical reaction characteristic time and flow characteristic time smaller. The flow is chemical non-equilibrium when the chemical reaction characteristic time and flow characteristic time is in the near order of magnitude. If the time ratio is greater than 10, the flow is close to the chemical equilibrium flow, and if the time ratio is less than 0.1, the flow is close to the chemical frozen flow. Under this ball head model, the chemical equilibrium flow requires Mach numbers greater than 8.0, 9.5, 10.8, 12.7, and 20.0 at heights of 0 km, 10 km, 20 km, 30 km, and 40 km respectively. The chemical frozen flow requires Mach numbers less than 8.4, 10.2, 12.9, and 17.5 at heights of 30 km, 40 km, 50 km, and 60 km respectively. The intermediate inflow state is chemical non-equilibrium flow.

    numerical simulation; high-temperature real gas effect; chemical non-equilibrium

    2097-1974(2023)01-0122-05

    10.7654/j.issn.2097-1974.20230124

    V211.3

    A

    2022-11-09;

    2022-12-04

    楊天鵬(1993-),男,博士,工程師,主要研究方向為飛行器氣動設計。

    嚴東升(1963-),男,研究員,主要研究方向為飛行力學設計。

    于江鵬(1995-),男,博士,工程師,主要研究方向為飛行器氣動設計。

    郭 陽(1991-),男,博士,工程師,主要研究方向為飛行器氣動設計。

    劉文伶(1963-),女,研究員,主要研究方向為飛行器總體設計。

    猜你喜歡
    化學特征模型
    一半模型
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    如何表達“特征”
    不忠誠的四個特征
    當代陜西(2019年10期)2019-06-03 10:12:04
    抓住特征巧觀察
    奇妙的化學
    奇妙的化學
    奇妙的化學
    3D打印中的模型分割與打包
    18禁在线无遮挡免费观看视频| 亚洲精品乱码久久久久久按摩| 欧美人与善性xxx| 丝袜在线中文字幕| 国产色婷婷99| av在线app专区| 午夜免费男女啪啪视频观看| 99re6热这里在线精品视频| 亚洲精品一二三| 国产在线免费精品| 国产永久视频网站| 2022亚洲国产成人精品| 亚洲精品视频女| 最新的欧美精品一区二区| 久久国产乱子免费精品| xxx大片免费视频| 男女国产视频网站| 亚洲综合精品二区| 精品久久久精品久久久| 亚洲精品日本国产第一区| 免费看日本二区| 国产在线一区二区三区精| 美女中出高潮动态图| 99热这里只有精品一区| 久久午夜福利片| 精品一区二区三区视频在线| 精品99又大又爽又粗少妇毛片| 丝袜在线中文字幕| 观看免费一级毛片| 嘟嘟电影网在线观看| 精品99又大又爽又粗少妇毛片| 欧美一级a爱片免费观看看| www.av在线官网国产| 国内少妇人妻偷人精品xxx网站| 美女主播在线视频| 最近最新中文字幕免费大全7| 日韩不卡一区二区三区视频在线| 97在线视频观看| 国产成人a∨麻豆精品| 国产淫片久久久久久久久| 亚洲国产成人一精品久久久| 日本欧美视频一区| 久久影院123| 伊人亚洲综合成人网| 亚洲国产精品一区二区三区在线| 高清午夜精品一区二区三区| 高清午夜精品一区二区三区| 成人毛片a级毛片在线播放| 久久久久久久久久久久大奶| 国产亚洲5aaaaa淫片| 日韩大片免费观看网站| 国产亚洲欧美精品永久| 一级a做视频免费观看| 99热这里只有是精品在线观看| 简卡轻食公司| 久久久久精品性色| 美女cb高潮喷水在线观看| 色视频www国产| 亚洲图色成人| 人人妻人人看人人澡| 亚洲经典国产精华液单| 肉色欧美久久久久久久蜜桃| 男的添女的下面高潮视频| 婷婷色麻豆天堂久久| 欧美日韩在线观看h| 国产熟女午夜一区二区三区 | 黄色一级大片看看| 免费看光身美女| 国产黄片美女视频| 国产成人精品无人区| 久久热精品热| 少妇被粗大猛烈的视频| 丰满迷人的少妇在线观看| 女性被躁到高潮视频| 国产精品蜜桃在线观看| 免费人成在线观看视频色| 亚洲丝袜综合中文字幕| 国产色爽女视频免费观看| 观看免费一级毛片| 亚洲国产欧美日韩在线播放 | 欧美+日韩+精品| 午夜免费观看性视频| 亚洲国产成人一精品久久久| 国产精品蜜桃在线观看| 三上悠亚av全集在线观看 | 国产成人精品福利久久| 在线观看免费视频网站a站| 成人午夜精彩视频在线观看| 涩涩av久久男人的天堂| 亚洲av成人精品一区久久| 观看免费一级毛片| 久久久a久久爽久久v久久| 亚洲av综合色区一区| 99热这里只有是精品50| 国产淫语在线视频| 赤兔流量卡办理| 91精品一卡2卡3卡4卡| 国产成人freesex在线| 高清视频免费观看一区二区| 91精品一卡2卡3卡4卡| 十八禁网站网址无遮挡 | 国产男女内射视频| 免费大片18禁| 亚洲四区av| 久久久久国产精品人妻一区二区| 久久精品国产亚洲网站| 人妻少妇偷人精品九色| 国产精品久久久久成人av| 亚洲av电影在线观看一区二区三区| 日韩大片免费观看网站| 深夜a级毛片| 亚洲精品,欧美精品| 亚洲av国产av综合av卡| 免费人妻精品一区二区三区视频| 国产 一区精品| 日韩制服骚丝袜av| av免费在线看不卡| 成人漫画全彩无遮挡| 男男h啪啪无遮挡| 91在线精品国自产拍蜜月| 国产av国产精品国产| 最后的刺客免费高清国语| 国产片特级美女逼逼视频| 午夜91福利影院| 精品少妇内射三级| 国产精品欧美亚洲77777| 女人精品久久久久毛片| 精品久久久久久电影网| 街头女战士在线观看网站| 99视频精品全部免费 在线| 国产伦在线观看视频一区| 丁香六月天网| 如何舔出高潮| 亚洲精品第二区| 久久 成人 亚洲| 精品国产一区二区久久| 国产女主播在线喷水免费视频网站| 伦理电影大哥的女人| 精品少妇黑人巨大在线播放| 美女中出高潮动态图| 少妇精品久久久久久久| 久久精品国产鲁丝片午夜精品| 一级av片app| 亚洲人与动物交配视频| 国产免费一区二区三区四区乱码| 啦啦啦视频在线资源免费观看| 亚洲av电影在线观看一区二区三区| 亚洲av.av天堂| 亚洲怡红院男人天堂| 一区在线观看完整版| 最近手机中文字幕大全| 最近中文字幕高清免费大全6| 9色porny在线观看| 欧美另类一区| 桃花免费在线播放| 女人精品久久久久毛片| 亚洲一级一片aⅴ在线观看| 人人妻人人澡人人看| 日本av免费视频播放| 国产精品福利在线免费观看| 91aial.com中文字幕在线观看| 亚洲欧美精品专区久久| 免费观看a级毛片全部| 亚洲av.av天堂| 成人国产av品久久久| 欧美高清成人免费视频www| 一级爰片在线观看| videos熟女内射| 另类亚洲欧美激情| 国产精品国产av在线观看| 成人综合一区亚洲| 秋霞伦理黄片| 久久毛片免费看一区二区三区| 国产精品人妻久久久久久| 日韩,欧美,国产一区二区三区| 波野结衣二区三区在线| 街头女战士在线观看网站| 国产精品秋霞免费鲁丝片| 一边亲一边摸免费视频| 精品少妇久久久久久888优播| 日本欧美国产在线视频| 欧美xxxx性猛交bbbb| 欧美日韩亚洲高清精品| 午夜视频国产福利| 国产一区二区三区综合在线观看 | 免费大片18禁| 国产成人a∨麻豆精品| 国产精品久久久久成人av| 免费在线观看成人毛片| 中文字幕亚洲精品专区| 久久久久久久亚洲中文字幕| 久久久欧美国产精品| 国产探花极品一区二区| 亚洲成色77777| 欧美精品一区二区大全| 久久免费观看电影| 日本色播在线视频| 爱豆传媒免费全集在线观看| 男女国产视频网站| 热re99久久国产66热| 下体分泌物呈黄色| 国产精品伦人一区二区| 一个人免费看片子| 日韩成人伦理影院| 免费看光身美女| 国产av码专区亚洲av| 22中文网久久字幕| 69精品国产乱码久久久| 色吧在线观看| 色视频在线一区二区三区| 久久午夜福利片| 韩国av在线不卡| 极品教师在线视频| 国产精品麻豆人妻色哟哟久久| 午夜视频国产福利| 女人久久www免费人成看片| 成年美女黄网站色视频大全免费 | 亚洲av日韩在线播放| 亚洲精品日韩av片在线观看| 中国三级夫妇交换| 一二三四中文在线观看免费高清| 一级毛片我不卡| 久久精品国产亚洲网站| 亚洲国产精品国产精品| 精品久久久噜噜| 午夜福利,免费看| 两个人的视频大全免费| 亚洲欧美成人综合另类久久久| 久久综合国产亚洲精品| 夜夜骑夜夜射夜夜干| 成人二区视频| 欧美老熟妇乱子伦牲交| 亚洲精品aⅴ在线观看| 亚洲精品乱码久久久久久按摩| 久久久久精品久久久久真实原创| 欧美国产精品一级二级三级 | 精品视频人人做人人爽| 人人妻人人爽人人添夜夜欢视频 | 2022亚洲国产成人精品| 成人黄色视频免费在线看| 国产伦精品一区二区三区视频9| 国产亚洲午夜精品一区二区久久| 欧美性感艳星| 日本爱情动作片www.在线观看| 菩萨蛮人人尽说江南好唐韦庄| 欧美精品国产亚洲| 美女福利国产在线| 简卡轻食公司| 久久精品国产a三级三级三级| 亚洲性久久影院| 男的添女的下面高潮视频| 色婷婷av一区二区三区视频| 18禁在线播放成人免费| 亚洲欧美一区二区三区国产| 黑人巨大精品欧美一区二区蜜桃 | 美女主播在线视频| 国产一区亚洲一区在线观看| 视频中文字幕在线观看| av不卡在线播放| 女的被弄到高潮叫床怎么办| 国产伦精品一区二区三区四那| 久久女婷五月综合色啪小说| 这个男人来自地球电影免费观看 | 曰老女人黄片| 九九爱精品视频在线观看| 久久这里有精品视频免费| 国产精品一区二区在线观看99| 亚洲精品亚洲一区二区| 日韩精品有码人妻一区| 精品一区在线观看国产| 一区二区三区四区激情视频| 日本欧美国产在线视频| 特大巨黑吊av在线直播| 日韩电影二区| a级片在线免费高清观看视频| 国产黄片美女视频| 伊人亚洲综合成人网| 亚洲国产av新网站| 国产视频首页在线观看| 能在线免费看毛片的网站| 精品久久久噜噜| 女人久久www免费人成看片| a级毛片在线看网站| 日韩中文字幕视频在线看片| 国精品久久久久久国模美| 一级av片app| 蜜桃在线观看..| 91午夜精品亚洲一区二区三区| 一本色道久久久久久精品综合| 18禁在线无遮挡免费观看视频| 欧美日韩亚洲高清精品| 国产探花极品一区二区| 如日韩欧美国产精品一区二区三区 | 丰满迷人的少妇在线观看| 国内精品宾馆在线| 亚洲精品一区蜜桃| .国产精品久久| 日本vs欧美在线观看视频 | 赤兔流量卡办理| 人妻夜夜爽99麻豆av| 国产成人精品婷婷| av不卡在线播放| 精品99又大又爽又粗少妇毛片| a级片在线免费高清观看视频| 交换朋友夫妻互换小说| 一本色道久久久久久精品综合| 国产免费视频播放在线视频| 青春草国产在线视频| 免费人成在线观看视频色| 亚洲av二区三区四区| 成人国产麻豆网| 免费观看无遮挡的男女| 国产精品99久久99久久久不卡 | 成人毛片60女人毛片免费| 大码成人一级视频| 欧美精品亚洲一区二区| 国产伦理片在线播放av一区| 婷婷色麻豆天堂久久| 色哟哟·www| 久久久久精品性色| 成人综合一区亚洲| 中文字幕人妻丝袜制服| 一区二区av电影网| 精品少妇久久久久久888优播| 亚洲欧美日韩另类电影网站| 国产成人精品一,二区| 精品国产一区二区三区久久久樱花| 美女主播在线视频| 欧美性感艳星| 99热全是精品| 久久毛片免费看一区二区三区| 天堂中文最新版在线下载| 成人国产av品久久久| 日本色播在线视频| 日本欧美国产在线视频| 男人狂女人下面高潮的视频| 成人午夜精彩视频在线观看| 日本黄大片高清| 日韩三级伦理在线观看| 人妻人人澡人人爽人人| 久久久久久久国产电影| av专区在线播放| 男女免费视频国产| 一本—道久久a久久精品蜜桃钙片| 国产成人a∨麻豆精品| 亚洲精品国产av蜜桃| 国产欧美亚洲国产| 亚洲天堂av无毛| 观看美女的网站| 成人二区视频| 不卡视频在线观看欧美| 日本av免费视频播放| 少妇丰满av| 国产高清有码在线观看视频| 国产无遮挡羞羞视频在线观看| 尾随美女入室| 免费人妻精品一区二区三区视频| 91久久精品国产一区二区三区| 在线 av 中文字幕| 在线观看人妻少妇| 免费观看av网站的网址| 成人国产av品久久久| 亚洲一级一片aⅴ在线观看| 国产成人a∨麻豆精品| 亚洲图色成人| 国产精品一区二区在线观看99| 国产精品秋霞免费鲁丝片| 国精品久久久久久国模美| 国产一级毛片在线| 欧美人与善性xxx| 欧美精品一区二区免费开放| 在线观看美女被高潮喷水网站| 免费看av在线观看网站| 国产成人freesex在线| 3wmmmm亚洲av在线观看| 黄色毛片三级朝国网站 | a级毛片在线看网站| 国产视频首页在线观看| 麻豆乱淫一区二区| 欧美人与善性xxx| av免费观看日本| 91精品伊人久久大香线蕉| 国产深夜福利视频在线观看| 久久国产精品大桥未久av | 国产精品久久久久久精品古装| 欧美激情国产日韩精品一区| 69精品国产乱码久久久| 啦啦啦视频在线资源免费观看| 中文字幕免费在线视频6| 女性生殖器流出的白浆| 亚洲美女视频黄频| 亚洲国产精品成人久久小说| 亚洲精华国产精华液的使用体验| 亚洲欧美精品专区久久| 国产综合精华液| 国产91av在线免费观看| 久久久久久久久久人人人人人人| 七月丁香在线播放| 国产老妇伦熟女老妇高清| h视频一区二区三区| 亚洲精品,欧美精品| 久久97久久精品| 老熟女久久久| 天堂中文最新版在线下载| 2018国产大陆天天弄谢| 国产中年淑女户外野战色| 午夜福利,免费看| 最黄视频免费看| 精品久久久精品久久久| 亚洲三级黄色毛片| 国产极品天堂在线| 免费观看性生交大片5| 国产在线一区二区三区精| 国产av码专区亚洲av| 高清av免费在线| 色视频www国产| 五月开心婷婷网| 丰满人妻一区二区三区视频av| 久久 成人 亚洲| 国产一区有黄有色的免费视频| 伊人亚洲综合成人网| 综合色丁香网| 高清在线视频一区二区三区| 国产成人精品无人区| 99九九线精品视频在线观看视频| 特大巨黑吊av在线直播| 少妇高潮的动态图| 十分钟在线观看高清视频www | 国产无遮挡羞羞视频在线观看| 人人妻人人爽人人添夜夜欢视频 | 七月丁香在线播放| 最近2019中文字幕mv第一页| 人人妻人人看人人澡| 久久精品国产自在天天线| 久久午夜福利片| 久久精品国产a三级三级三级| 日韩人妻高清精品专区| 啦啦啦中文免费视频观看日本| 精品亚洲成a人片在线观看| 国产永久视频网站| 日本av手机在线免费观看| 极品少妇高潮喷水抽搐| 午夜福利网站1000一区二区三区| 成人18禁高潮啪啪吃奶动态图 | 午夜久久久在线观看| 日韩一区二区三区影片| 成人亚洲精品一区在线观看| xxx大片免费视频| 国产69精品久久久久777片| 欧美日韩综合久久久久久| 狂野欧美激情性xxxx在线观看| h日本视频在线播放| 午夜福利在线观看免费完整高清在| 一级毛片电影观看| 亚洲不卡免费看| 久久热精品热| 亚洲婷婷狠狠爱综合网| 少妇的逼水好多| 国产伦精品一区二区三区视频9| 国产 精品1| 国产成人精品福利久久| 少妇的逼好多水| 免费黄网站久久成人精品| 国产永久视频网站| av国产精品久久久久影院| 免费在线观看成人毛片| 国产在线视频一区二区| 国产伦精品一区二区三区视频9| 天美传媒精品一区二区| 久久久a久久爽久久v久久| 只有这里有精品99| 国产精品久久久久成人av| 精品少妇黑人巨大在线播放| 国产亚洲欧美精品永久| 自线自在国产av| 久久99精品国语久久久| 亚洲国产最新在线播放| 人人妻人人看人人澡| 丰满人妻一区二区三区视频av| 在线观看www视频免费| 国产免费一级a男人的天堂| 亚洲精品久久午夜乱码| 国产淫语在线视频| 日本wwww免费看| 国产精品一二三区在线看| 全区人妻精品视频| 黄色怎么调成土黄色| 大又大粗又爽又黄少妇毛片口| 日本午夜av视频| 国产精品一区www在线观看| 亚洲av二区三区四区| 成人漫画全彩无遮挡| 在线观看一区二区三区激情| 午夜av观看不卡| 熟女人妻精品中文字幕| 成人亚洲欧美一区二区av| 亚洲精品国产色婷婷电影| 国产亚洲av片在线观看秒播厂| 国产精品无大码| 老司机影院成人| 三级经典国产精品| 成人特级av手机在线观看| 高清在线视频一区二区三区| 深夜a级毛片| 麻豆成人午夜福利视频| 国产成人freesex在线| 亚洲欧洲精品一区二区精品久久久 | 美女cb高潮喷水在线观看| 十八禁网站网址无遮挡 | 国产亚洲av片在线观看秒播厂| av福利片在线| 最近的中文字幕免费完整| 免费观看在线日韩| 最近中文字幕高清免费大全6| 精品久久久噜噜| 又大又黄又爽视频免费| 亚洲精品第二区| 我的老师免费观看完整版| 亚洲精品国产av成人精品| 国产有黄有色有爽视频| 能在线免费看毛片的网站| 不卡视频在线观看欧美| 久久久精品94久久精品| 国产一区有黄有色的免费视频| 看非洲黑人一级黄片| 秋霞在线观看毛片| 永久免费av网站大全| 美女cb高潮喷水在线观看| 黄色日韩在线| 少妇裸体淫交视频免费看高清| av视频免费观看在线观看| 各种免费的搞黄视频| 久久久久久久国产电影| 黄片无遮挡物在线观看| 少妇 在线观看| 亚洲欧美一区二区三区黑人 | 日韩av免费高清视频| 91成人精品电影| 亚洲av成人精品一二三区| 国产精品久久久久久精品古装| 91午夜精品亚洲一区二区三区| 性色avwww在线观看| 精品酒店卫生间| 天美传媒精品一区二区| 亚洲人成网站在线观看播放| 九色成人免费人妻av| 夜夜爽夜夜爽视频| 亚洲精品色激情综合| 亚洲不卡免费看| 精华霜和精华液先用哪个| 国产乱来视频区| 日韩av在线免费看完整版不卡| 男人爽女人下面视频在线观看| 制服丝袜香蕉在线| 99久久精品国产国产毛片| 看十八女毛片水多多多| 久久6这里有精品| 麻豆成人av视频| 多毛熟女@视频| 大又大粗又爽又黄少妇毛片口| 最近最新中文字幕免费大全7| 精品亚洲乱码少妇综合久久| av免费观看日本| 啦啦啦啦在线视频资源| 日产精品乱码卡一卡2卡三| 熟女人妻精品中文字幕| 欧美精品一区二区大全| 成年人午夜在线观看视频| 在线看a的网站| 男人和女人高潮做爰伦理| 插逼视频在线观看| 日韩欧美精品免费久久| 色94色欧美一区二区| 91精品伊人久久大香线蕉| 大香蕉久久网| 免费观看性生交大片5| 国产成人一区二区在线| 18禁在线播放成人免费| 日韩在线高清观看一区二区三区| 一级a做视频免费观看| 男人添女人高潮全过程视频| 97超视频在线观看视频| 丁香六月天网| 丰满人妻一区二区三区视频av| 国产色爽女视频免费观看| 18禁裸乳无遮挡动漫免费视频| av在线老鸭窝| 国产男人的电影天堂91| 久久精品熟女亚洲av麻豆精品| 高清av免费在线| 久热这里只有精品99| av一本久久久久| 午夜福利,免费看| 国产黄片视频在线免费观看| 亚洲国产日韩一区二区| 一二三四中文在线观看免费高清| 黄色日韩在线| 最近手机中文字幕大全| 男女免费视频国产| 欧美精品亚洲一区二区| 黄片无遮挡物在线观看| 国产精品蜜桃在线观看| 日韩制服骚丝袜av| 在线观看免费视频网站a站| 精品少妇久久久久久888优播| 在线天堂最新版资源| 少妇人妻精品综合一区二区| 亚洲天堂av无毛| 国产深夜福利视频在线观看| 97超视频在线观看视频| 日韩精品免费视频一区二区三区 | 国产免费一区二区三区四区乱码| 久久国产精品大桥未久av | 桃花免费在线播放| 精品一区二区三卡|