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

    沿駐點(diǎn)線化學(xué)非平衡程度評估方法

    2023-03-09 11:38:02楊天鵬嚴(yán)東升于江鵬劉文伶
    關(guān)鍵詞:化學(xué)特征模型

    楊天鵬,嚴(yán)東升,于江鵬,郭 陽,劉文伶

    沿駐點(diǎn)線化學(xué)非平衡程度評估方法

    楊天鵬,嚴(yán)東升,于江鵬,郭 陽,劉文伶

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

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

    數(shù)值模擬;高溫真實(shí)氣體效應(yīng);化學(xué)非平衡

    0 引 言

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

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

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

    1 沿駐點(diǎn)線化學(xué)反應(yīng)流動特征分析

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

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

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

    圖1 計(jì)算網(wǎng)格

    表1 化學(xué)反應(yīng)方程

    Tab.1 Chemical Reaction Equation

    序號化學(xué)反應(yīng)方程式 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種情況下沿駐點(diǎn)線化學(xué)反應(yīng)能量釋放情況。這里定義反應(yīng)的化學(xué)反應(yīng)功率密度E

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

    續(xù)圖2

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

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

    續(xù)圖3

    2 化學(xué)反應(yīng)特征時間評估

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

    圖4給出了取=15、高度30 km條件下的組元變化歷程,給出了以O(shè)2、N2和NO為參考,組元濃度分別達(dá)到80%、90%和99%平衡狀態(tài)情況下的時刻位置。可以看到將化學(xué)反應(yīng)特征時間定義為容器中的氧氣和氮?dú)膺_(dá)到80%平衡狀態(tài)時的反應(yīng)時間較為合適,之后化學(xué)反應(yīng)處于反應(yīng)緩慢的平臺區(qū),化學(xué)反應(yīng)導(dǎo)致的組元變化和能量變化主要集中于前80%部分,若選擇更高的非平衡標(biāo)準(zhǔn)會使得化學(xué)反應(yīng)特征時間成倍增加,化學(xué)非平衡評估模型的敏感性降低。

    圖4 Ma=15算例沿駐點(diǎn)線組元變化歷程

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

    圖5 化學(xué)反應(yīng)特征時間

    3 流動特征時間評估

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

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

    圖6 流動特征時間

    4 化學(xué)非平衡程度模型

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

    圖7 化學(xué)反應(yīng)特征時間與流動特征時間之比

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

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

    續(xù)圖8

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

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

    5 結(jié) 論

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

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

    [1] 歐陽水吾, 謝中強(qiáng). 高溫非平衡空氣繞流[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] 中國科學(xué)院. 中國學(xué)科發(fā)展戰(zhàn)略新型飛行器中的關(guān)鍵力學(xué)問題[M]. 北京: 科學(xué)出版社, 2018.

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

    [5] 王江峰, 等. 高超聲速復(fù)雜氣動問題數(shù)值方法研究進(jìn)展[J]. 航空學(xué)報, 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. 高超聲速和高溫氣體動力學(xué)(第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-),男,博士,工程師,主要研究方向?yàn)轱w行器氣動設(shè)計(jì)。

    嚴(yán)東升(1963-),男,研究員,主要研究方向?yàn)轱w行力學(xué)設(shè)計(jì)。

    于江鵬(1995-),男,博士,工程師,主要研究方向?yàn)轱w行器氣動設(shè)計(jì)。

    郭 陽(1991-),男,博士,工程師,主要研究方向?yàn)轱w行器氣動設(shè)計(jì)。

    劉文伶(1963-),女,研究員,主要研究方向?yàn)轱w行器總體設(shè)計(jì)。

    猜你喜歡
    化學(xué)特征模型
    一半模型
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    如何表達(dá)“特征”
    不忠誠的四個特征
    抓住特征巧觀察
    奇妙的化學(xué)
    奇妙的化學(xué)
    奇妙的化學(xué)
    3D打印中的模型分割與打包
    a级毛片免费高清观看在线播放| 亚洲不卡免费看| 少妇人妻一区二区三区视频| 中文字幕av成人在线电影| 亚洲av二区三区四区| 大香蕉久久网| 久久久欧美国产精品| av线在线观看网站| 国产黄色视频一区二区在线观看| 精品久久国产蜜桃| 国产亚洲91精品色在线| 精品一区二区免费观看| 少妇被粗大猛烈的视频| 一区在线观看完整版| 精华霜和精华液先用哪个| 黄片无遮挡物在线观看| 这个男人来自地球电影免费观看 | 中文天堂在线官网| 久久人妻熟女aⅴ| 99精国产麻豆久久婷婷| 日韩国内少妇激情av| 亚洲性久久影院| 色5月婷婷丁香| 国产毛片在线视频| 亚洲av在线观看美女高潮| 久久久久久人妻| 国产精品久久久久久精品古装| 精品少妇久久久久久888优播| 在线观看三级黄色| 亚洲精品自拍成人| 国产黄片视频在线免费观看| 久久久久久久久大av| 永久免费av网站大全| 中文字幕精品免费在线观看视频 | 久久精品久久久久久噜噜老黄| 免费观看性生交大片5| 免费大片18禁| 国产男女内射视频| 噜噜噜噜噜久久久久久91| 有码 亚洲区| 夫妻午夜视频| av一本久久久久| 国产视频首页在线观看| 国产精品国产三级国产av玫瑰| 五月伊人婷婷丁香| 观看av在线不卡| 亚洲精品亚洲一区二区| 纵有疾风起免费观看全集完整版| 人体艺术视频欧美日本| 国产成人精品久久久久久| 中文天堂在线官网| 日韩一本色道免费dvd| 中文字幕av成人在线电影| 边亲边吃奶的免费视频| 天天躁夜夜躁狠狠久久av| 我的老师免费观看完整版| 国内少妇人妻偷人精品xxx网站| 大陆偷拍与自拍| 欧美日韩综合久久久久久| 最近最新中文字幕大全电影3| 91久久精品国产一区二区成人| 中文天堂在线官网| 联通29元200g的流量卡| 午夜激情福利司机影院| 丰满乱子伦码专区| 久久久久精品久久久久真实原创| 亚洲,欧美,日韩| 如何舔出高潮| 一区二区三区乱码不卡18| av一本久久久久| 欧美日韩在线观看h| 欧美日韩综合久久久久久| 青春草视频在线免费观看| 热re99久久精品国产66热6| 久久人人爽人人片av| 欧美丝袜亚洲另类| 亚洲av综合色区一区| 欧美日韩亚洲高清精品| 久久精品国产亚洲av天美| 国产色婷婷99| av专区在线播放| 国产精品.久久久| 成人毛片60女人毛片免费| 亚洲精品久久久久久婷婷小说| 日本-黄色视频高清免费观看| 大陆偷拍与自拍| 99精国产麻豆久久婷婷| 国产爱豆传媒在线观看| 日本av免费视频播放| 亚洲伊人久久精品综合| 亚洲自偷自拍三级| 自拍欧美九色日韩亚洲蝌蚪91 | 国产男女内射视频| 色视频在线一区二区三区| 丰满迷人的少妇在线观看| 丰满迷人的少妇在线观看| 97热精品久久久久久| 国产精品一区www在线观看| 啦啦啦啦在线视频资源| 又爽又黄a免费视频| 欧美变态另类bdsm刘玥| 少妇高潮的动态图| 女人久久www免费人成看片| 亚洲欧美日韩东京热| 黄色配什么色好看| 99精国产麻豆久久婷婷| 国内少妇人妻偷人精品xxx网站| 亚洲精品日韩在线中文字幕| 免费观看a级毛片全部| 免费播放大片免费观看视频在线观看| 精品国产乱码久久久久久小说| 亚洲av男天堂| 国产精品爽爽va在线观看网站| 成人特级av手机在线观看| 欧美zozozo另类| 国产在视频线精品| 青春草国产在线视频| 老师上课跳d突然被开到最大视频| 极品少妇高潮喷水抽搐| 人妻一区二区av| 亚洲欧美精品自产自拍| 最黄视频免费看| 国产成人精品婷婷| 久久久久久久精品精品| 大香蕉97超碰在线| 99热全是精品| 欧美成人午夜免费资源| 在线观看免费视频网站a站| 亚洲国产日韩一区二区| 在线看a的网站| 草草在线视频免费看| 青青草视频在线视频观看| 女性被躁到高潮视频| 香蕉精品网在线| 久久久久久久亚洲中文字幕| 久久久久国产网址| 久久久久久久久久久免费av| 久久久午夜欧美精品| 亚洲国产欧美人成| 夜夜骑夜夜射夜夜干| 亚洲人成网站在线播| 国产伦在线观看视频一区| 国产av精品麻豆| 久久精品久久久久久噜噜老黄| 日韩大片免费观看网站| 一二三四中文在线观看免费高清| 97精品久久久久久久久久精品| 熟妇人妻不卡中文字幕| 大又大粗又爽又黄少妇毛片口| 色综合色国产| 欧美精品一区二区免费开放| 日韩av免费高清视频| 啦啦啦视频在线资源免费观看| 精品久久久精品久久久| 纯流量卡能插随身wifi吗| 日日摸夜夜添夜夜添av毛片| 中文字幕av成人在线电影| 免费播放大片免费观看视频在线观看| 九草在线视频观看| 91午夜精品亚洲一区二区三区| 男女免费视频国产| 精品久久国产蜜桃| 久久韩国三级中文字幕| 亚洲美女搞黄在线观看| 伦精品一区二区三区| 两个人的视频大全免费| 搡女人真爽免费视频火全软件| 国产精品精品国产色婷婷| 成人黄色视频免费在线看| 欧美日韩视频精品一区| 久久综合国产亚洲精品| 亚洲成人一二三区av| 国产成人精品福利久久| 夜夜看夜夜爽夜夜摸| 国产精品国产三级国产av玫瑰| 亚洲aⅴ乱码一区二区在线播放| 国产69精品久久久久777片| 国产精品麻豆人妻色哟哟久久| 国产精品一区二区在线观看99| 26uuu在线亚洲综合色| 视频中文字幕在线观看| 高清欧美精品videossex| 观看免费一级毛片| www.av在线官网国产| 中文欧美无线码| 九色成人免费人妻av| 五月天丁香电影| 老女人水多毛片| 欧美区成人在线视频| 国产精品一及| 亚洲欧美清纯卡通| 久久精品熟女亚洲av麻豆精品| 亚洲欧美日韩东京热| 久久热精品热| 自拍欧美九色日韩亚洲蝌蚪91 | 欧美日韩视频精品一区| 我要看日韩黄色一级片| 秋霞在线观看毛片| 人人妻人人爽人人添夜夜欢视频 | 看十八女毛片水多多多| 日日撸夜夜添| 国产老妇伦熟女老妇高清| 久久婷婷青草| 午夜日本视频在线| 天美传媒精品一区二区| 天堂俺去俺来也www色官网| 免费少妇av软件| 日本免费在线观看一区| 国产欧美另类精品又又久久亚洲欧美| 各种免费的搞黄视频| 免费播放大片免费观看视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 草草在线视频免费看| 黄片wwwwww| 日韩人妻高清精品专区| 男女下面进入的视频免费午夜| av免费观看日本| 中文欧美无线码| 岛国毛片在线播放| 在线天堂最新版资源| 身体一侧抽搐| 亚洲欧美精品专区久久| 亚洲av欧美aⅴ国产| 日韩中文字幕视频在线看片 | 男人添女人高潮全过程视频| 午夜福利高清视频| 97在线视频观看| 久久99热6这里只有精品| 国产成人一区二区在线| 国产真实伦视频高清在线观看| 国产视频首页在线观看| 少妇 在线观看| 国产av一区二区精品久久 | 亚州av有码| 国产精品久久久久成人av| 一个人免费看片子| 插逼视频在线观看| 又大又黄又爽视频免费| 亚洲婷婷狠狠爱综合网| 精品亚洲乱码少妇综合久久| h日本视频在线播放| www.色视频.com| 卡戴珊不雅视频在线播放| 国产69精品久久久久777片| 国产又色又爽无遮挡免| 久久热精品热| 最近最新中文字幕免费大全7| 久久久久久久精品精品| 免费黄网站久久成人精品| 久久国产精品大桥未久av | 久久国产乱子免费精品| 亚洲精品国产色婷婷电影| 在线免费观看不下载黄p国产| 在线观看免费视频网站a站| 亚洲精品,欧美精品| 久久97久久精品| 老师上课跳d突然被开到最大视频| 天堂中文最新版在线下载| 国产精品爽爽va在线观看网站| 国产精品国产三级国产专区5o| 国产午夜精品一二区理论片| 看非洲黑人一级黄片| av国产久精品久网站免费入址| 蜜桃亚洲精品一区二区三区| 久久久久久久久久成人| 亚洲怡红院男人天堂| 国产淫片久久久久久久久| 午夜免费男女啪啪视频观看| 婷婷色综合大香蕉| 日韩免费高清中文字幕av| 成人午夜精彩视频在线观看| 日韩中字成人| 欧美成人一区二区免费高清观看| 新久久久久国产一级毛片| 午夜日本视频在线| 91精品伊人久久大香线蕉| 伊人久久国产一区二区| 国产女主播在线喷水免费视频网站| 男女边吃奶边做爰视频| 亚洲av不卡在线观看| 国产熟女欧美一区二区| 高清日韩中文字幕在线| 亚洲欧美中文字幕日韩二区| 人妻一区二区av| 国内揄拍国产精品人妻在线| 五月天丁香电影| 久久亚洲国产成人精品v| 亚洲伊人久久精品综合| 国产精品一区www在线观看| 亚洲久久久国产精品| 亚洲综合色惰| 国产黄色视频一区二区在线观看| 一级毛片aaaaaa免费看小| 亚洲激情五月婷婷啪啪| 综合色丁香网| 国产欧美日韩一区二区三区在线 | 国产精品一区二区在线不卡| 韩国高清视频一区二区三区| 欧美97在线视频| 欧美亚洲 丝袜 人妻 在线| 一级毛片久久久久久久久女| 亚洲精品第二区| 伦理电影大哥的女人| 久久久久视频综合| 我要看黄色一级片免费的| 免费人成在线观看视频色| 国产 一区精品| 毛片女人毛片| 国产又色又爽无遮挡免| 国产精品秋霞免费鲁丝片| 久久精品夜色国产| 一本久久精品| 老师上课跳d突然被开到最大视频| 最近中文字幕2019免费版| 久久99精品国语久久久| 国产精品爽爽va在线观看网站| 青青草视频在线视频观看| 成人综合一区亚洲| 成人毛片a级毛片在线播放| 亚洲av.av天堂| 一级黄片播放器| 国产av码专区亚洲av| 日韩人妻高清精品专区| 久久 成人 亚洲| 亚洲人成网站在线观看播放| 国产成人91sexporn| 美女福利国产在线 | 国产高清三级在线| 国产精品国产三级专区第一集| 国产精品女同一区二区软件| 久久精品国产a三级三级三级| 亚洲av男天堂| 简卡轻食公司| 久久鲁丝午夜福利片| 亚洲一级一片aⅴ在线观看| 久久精品夜色国产| 国产黄片美女视频| 亚洲内射少妇av| 精品99又大又爽又粗少妇毛片| av福利片在线观看| 国产精品欧美亚洲77777| 熟妇人妻不卡中文字幕| 国内少妇人妻偷人精品xxx网站| 深爱激情五月婷婷| 18禁动态无遮挡网站| 国产黄片美女视频| 色哟哟·www| 久热久热在线精品观看| 寂寞人妻少妇视频99o| 天美传媒精品一区二区| 黄色欧美视频在线观看| av国产免费在线观看| 日韩亚洲欧美综合| 欧美激情极品国产一区二区三区 | 国产成人a∨麻豆精品| 人妻 亚洲 视频| 成人国产麻豆网| 51国产日韩欧美| 久久女婷五月综合色啪小说| 色哟哟·www| 亚洲精品乱码久久久v下载方式| 我的老师免费观看完整版| 国产高清有码在线观看视频| 免费人成在线观看视频色| 国内揄拍国产精品人妻在线| 狂野欧美激情性bbbbbb| 黄色日韩在线| 国产亚洲最大av| 国产亚洲av片在线观看秒播厂| 九九在线视频观看精品| 伦理电影免费视频| 97在线人人人人妻| av免费观看日本| 精品一区二区三区视频在线| 亚洲国产日韩一区二区| 亚洲av电影在线观看一区二区三区| 妹子高潮喷水视频| 国产av码专区亚洲av| 久久久久久久精品精品| 国产成人午夜福利电影在线观看| 久久热精品热| 王馨瑶露胸无遮挡在线观看| 亚洲无线观看免费| 亚洲国产欧美在线一区| 极品教师在线视频| 午夜福利高清视频| 18禁在线无遮挡免费观看视频| 日韩视频在线欧美| 国内揄拍国产精品人妻在线| 亚洲欧美一区二区三区国产| 美女脱内裤让男人舔精品视频| 国产视频首页在线观看| 国产一区有黄有色的免费视频| 婷婷色麻豆天堂久久| 中文字幕精品免费在线观看视频 | 亚洲精品,欧美精品| 色婷婷久久久亚洲欧美| h日本视频在线播放| 亚洲国产精品一区三区| 男人添女人高潮全过程视频| 日本wwww免费看| 麻豆精品久久久久久蜜桃| 寂寞人妻少妇视频99o| 男男h啪啪无遮挡| 亚洲精品一二三| xxx大片免费视频| 亚洲图色成人| 欧美日韩视频精品一区| 国产黄色视频一区二区在线观看| 人人妻人人爽人人添夜夜欢视频 | 精品国产乱码久久久久久小说| 久久久久精品性色| 最近中文字幕2019免费版| 91在线精品国自产拍蜜月| 午夜免费男女啪啪视频观看| 久久久久视频综合| 欧美高清性xxxxhd video| 在线观看免费日韩欧美大片 | 精品久久久久久电影网| 亚洲精品中文字幕在线视频 | 国产在线免费精品| 欧美性感艳星| 99久久精品一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 免费看av在线观看网站| 亚洲激情五月婷婷啪啪| 热99国产精品久久久久久7| 国产黄片美女视频| 午夜免费鲁丝| 另类亚洲欧美激情| 亚洲自偷自拍三级| 尤物成人国产欧美一区二区三区| 亚洲精品,欧美精品| 水蜜桃什么品种好| 欧美xxⅹ黑人| 午夜老司机福利剧场| 成人高潮视频无遮挡免费网站| 一边亲一边摸免费视频| 日日啪夜夜爽| 日韩不卡一区二区三区视频在线| 国产成人aa在线观看| 成年美女黄网站色视频大全免费 | 免费av不卡在线播放| www.av在线官网国产| 久久青草综合色| 我的老师免费观看完整版| 亚洲成人手机| 亚洲精华国产精华液的使用体验| 蜜臀久久99精品久久宅男| 精品一区二区三区视频在线| 亚洲第一区二区三区不卡| 丰满少妇做爰视频| 亚洲中文av在线| 天美传媒精品一区二区| 久久久欧美国产精品| 日韩成人av中文字幕在线观看| 成人高潮视频无遮挡免费网站| 日日摸夜夜添夜夜爱| 精品国产一区二区三区久久久樱花 | 一级毛片电影观看| 日韩,欧美,国产一区二区三区| 国产高清有码在线观看视频| 一级毛片黄色毛片免费观看视频| 在线 av 中文字幕| 国产精品福利在线免费观看| 久久久久久人妻| 精品午夜福利在线看| 久久久欧美国产精品| 女人十人毛片免费观看3o分钟| 人人妻人人澡人人爽人人夜夜| 欧美日韩国产mv在线观看视频 | 亚洲综合精品二区| 人人妻人人澡人人爽人人夜夜| 国产精品爽爽va在线观看网站| 中文在线观看免费www的网站| av播播在线观看一区| 国产精品久久久久久精品电影小说 | 久久精品人妻少妇| 99久国产av精品国产电影| 国产白丝娇喘喷水9色精品| 啦啦啦中文免费视频观看日本| 另类亚洲欧美激情| 国产亚洲午夜精品一区二区久久| 乱码一卡2卡4卡精品| 亚洲国产av新网站| 国产久久久一区二区三区| av国产免费在线观看| 嘟嘟电影网在线观看| 欧美精品一区二区免费开放| 热re99久久精品国产66热6| 只有这里有精品99| 99热这里只有是精品在线观看| 国产爱豆传媒在线观看| 国产片特级美女逼逼视频| 日本-黄色视频高清免费观看| av不卡在线播放| 亚洲av中文av极速乱| 国产爽快片一区二区三区| 中文精品一卡2卡3卡4更新| 欧美日本视频| 欧美成人一区二区免费高清观看| 午夜福利网站1000一区二区三区| 免费黄色在线免费观看| 国产精品av视频在线免费观看| 尤物成人国产欧美一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91 | 2022亚洲国产成人精品| av线在线观看网站| 夫妻午夜视频| 成人免费观看视频高清| 黄色视频在线播放观看不卡| 久久精品国产自在天天线| 成人美女网站在线观看视频| 熟女电影av网| 青春草国产在线视频| 少妇人妻久久综合中文| 18+在线观看网站| 男男h啪啪无遮挡| 久久久久久九九精品二区国产| 国产亚洲一区二区精品| 国产成人91sexporn| 人人妻人人看人人澡| 国产精品久久久久久av不卡| 国产乱人视频| 欧美成人午夜免费资源| 激情五月婷婷亚洲| .国产精品久久| 久久久久国产精品人妻一区二区| 精华霜和精华液先用哪个| 国产成人a区在线观看| 国产精品女同一区二区软件| 免费观看在线日韩| 免费看日本二区| 五月玫瑰六月丁香| 日本黄色日本黄色录像| 亚洲国产精品成人久久小说| 少妇丰满av| 秋霞伦理黄片| 六月丁香七月| 日产精品乱码卡一卡2卡三| 又爽又黄a免费视频| 中国国产av一级| 在线天堂最新版资源| 色视频www国产| 日韩三级伦理在线观看| 亚洲欧美日韩东京热| 亚洲国产精品成人久久小说| 少妇丰满av| 欧美变态另类bdsm刘玥| av国产精品久久久久影院| 国产女主播在线喷水免费视频网站| 亚洲成人一二三区av| 777米奇影视久久| 久久鲁丝午夜福利片| 下体分泌物呈黄色| 国产伦精品一区二区三区视频9| 中文字幕人妻熟人妻熟丝袜美| 麻豆国产97在线/欧美| 久久鲁丝午夜福利片| 人妻少妇偷人精品九色| 久久99热这里只有精品18| 日韩在线高清观看一区二区三区| 欧美成人a在线观看| 国产免费福利视频在线观看| 99热6这里只有精品| 五月伊人婷婷丁香| 又黄又爽又刺激的免费视频.| 国产成人精品婷婷| 久久精品国产a三级三级三级| 中文字幕久久专区| 大码成人一级视频| 免费黄频网站在线观看国产| 五月玫瑰六月丁香| 久热这里只有精品99| 国产精品欧美亚洲77777| 中国美白少妇内射xxxbb| 女人十人毛片免费观看3o分钟| 免费大片黄手机在线观看| 2021少妇久久久久久久久久久| 黄色配什么色好看| 亚洲成人手机| av天堂中文字幕网| 国产亚洲精品久久久com| 午夜福利高清视频| 蜜桃亚洲精品一区二区三区| 日本vs欧美在线观看视频 | 国产人妻一区二区三区在| av.在线天堂| 51国产日韩欧美| 国产黄频视频在线观看| 狠狠精品人妻久久久久久综合| 又大又黄又爽视频免费| 国产欧美日韩精品一区二区| 最近中文字幕高清免费大全6| 国产高清不卡午夜福利| 久久久久精品性色| 亚洲欧美精品专区久久| av线在线观看网站| 有码 亚洲区| 久久鲁丝午夜福利片| 成年人午夜在线观看视频| 日韩一本色道免费dvd| 久热久热在线精品观看| 亚洲av.av天堂| 日韩在线高清观看一区二区三区| 久久99热这里只频精品6学生| 伦理电影大哥的女人| 国产精品嫩草影院av在线观看| 汤姆久久久久久久影院中文字幕| 欧美日韩亚洲高清精品| 国产成人免费观看mmmm| 日日摸夜夜添夜夜添av毛片| 美女福利国产在线 | 又粗又硬又长又爽又黄的视频| 国产日韩欧美在线精品|