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

    誘導(dǎo)輪超同步旋轉(zhuǎn)空化傳播機(jī)理

    2022-05-14 06:48:42李春樂許開富
    火箭推進(jìn) 2022年2期
    關(guān)鍵詞:液流空化前緣

    項(xiàng) 樂,李春樂,許開富,陳 暉

    (西安航天動(dòng)力研究所,陜西 西安 710100)

    0 引言

    誘導(dǎo)輪是安裝在渦輪泵主泵上游的一種軸流泵,通過適當(dāng)增壓提升來流壓力,從而提高主泵的空化性能。由于其安裝角小、稠度大等結(jié)構(gòu)特點(diǎn),能夠在一定程度的空化條件下正常工作,而不發(fā)生明顯的揚(yáng)程性能下降。然而即便在設(shè)計(jì)工況下,誘導(dǎo)輪內(nèi)也可能發(fā)生空化誘發(fā)的流動(dòng)不穩(wěn)定現(xiàn)象,例如我國新一代高壓補(bǔ)燃循環(huán)液氧煤油液體火箭發(fā)動(dòng)機(jī)研制過程中,超同步旋轉(zhuǎn)空化現(xiàn)象長(zhǎng)期以來是導(dǎo)致渦輪泵振動(dòng)量級(jí)過高的重要激振源之一,我國在研的500 tf級(jí)液氧煤油發(fā)動(dòng)機(jī)渦輪氧泵水試和熱試車中均出現(xiàn)了明顯的超同步旋轉(zhuǎn)空化現(xiàn)象,其特征頻率幅值顯著大于葉輪轉(zhuǎn)頻幅值,導(dǎo)致振動(dòng)幅值偏高,目前尚無有效技術(shù)手段徹底抑制,因此亟需對(duì)誘導(dǎo)輪超同步旋轉(zhuǎn)空化現(xiàn)象發(fā)生機(jī)理展開深入研究,在渦輪泵設(shè)計(jì)階段對(duì)該現(xiàn)象進(jìn)行抑制,從而提升發(fā)動(dòng)機(jī)的可靠性。

    超同步旋轉(zhuǎn)空化是旋轉(zhuǎn)空化現(xiàn)象的一種,是相對(duì)于同步旋轉(zhuǎn)空化和次同步旋轉(zhuǎn)空化而言的。旋轉(zhuǎn)空化現(xiàn)象是指誘導(dǎo)輪3個(gè)葉片上空化區(qū)波動(dòng)沿周向傳播的一種流動(dòng)不穩(wěn)定現(xiàn)象,按傳播的頻率不同可以分為3類:空化區(qū)波動(dòng)傳播頻率快于誘導(dǎo)輪轉(zhuǎn)速(或者傳播方向與誘導(dǎo)輪旋轉(zhuǎn)方向相同)則是超同步旋轉(zhuǎn)空化,對(duì)應(yīng)的頻率一般大于1;空化區(qū)波動(dòng)傳播頻率與誘導(dǎo)輪轉(zhuǎn)速一致,則是同步旋轉(zhuǎn)空化,空化區(qū)波動(dòng)相對(duì)于葉片是固定的,對(duì)應(yīng)的頻率為1;空化區(qū)波動(dòng)傳播頻率小于誘導(dǎo)輪轉(zhuǎn)速(或者傳播方向與誘導(dǎo)輪旋轉(zhuǎn)方向相反)則是次同步旋轉(zhuǎn)空化,對(duì)應(yīng)的頻率小于1,目前比較常見的是超同步旋轉(zhuǎn)空化現(xiàn)象。Kamijo等第一次在LE-7氧泵中發(fā)現(xiàn)超同步的軸振動(dòng),并通過一系列實(shí)驗(yàn)研究證實(shí)該振動(dòng)是由超同步旋轉(zhuǎn)空化引起。1998年日本HII火箭第8次發(fā)射失利,經(jīng)過分析確認(rèn)氫泵誘導(dǎo)輪旋轉(zhuǎn)空化誘發(fā)的壓力脈動(dòng)頻率與葉片的固有頻率發(fā)生了共振,使葉片發(fā)生了疲勞斷裂,導(dǎo)致發(fā)動(dòng)機(jī)停車,引起這次發(fā)射事故。該事件也使得超同步旋轉(zhuǎn)空化現(xiàn)象逐漸引起航天界的重視,Tsujimoto等對(duì)誘導(dǎo)輪內(nèi)的空化不穩(wěn)定現(xiàn)象開展了全面的研究,對(duì)旋轉(zhuǎn)空化的特征進(jìn)行詳細(xì)描述,并提出了完整的空化不穩(wěn)定發(fā)生判據(jù)。Choi等通過對(duì)入口殼體的結(jié)構(gòu)修正實(shí)現(xiàn)了旋轉(zhuǎn)空化較好的抑制效果,但是其機(jī)理并不清楚。Angelo等研究了多個(gè)誘導(dǎo)輪內(nèi)發(fā)生的空化不穩(wěn)定現(xiàn)象。陳暉等首次在我國某型號(hào)發(fā)動(dòng)機(jī)試車數(shù)據(jù)中識(shí)別出旋轉(zhuǎn)空化特征頻率。Li X等研究了葉頂間隙大小對(duì)同步旋轉(zhuǎn)空化的影響。Lettieri等基于可視化實(shí)驗(yàn)發(fā)現(xiàn)葉尖泄漏渦與旋轉(zhuǎn)空化的發(fā)生有密切聯(lián)系。Kim等利用PIV測(cè)量了誘導(dǎo)輪三葉片旋轉(zhuǎn)空化下葉尖的流場(chǎng),為揭示旋轉(zhuǎn)空化的機(jī)理提供了重要的實(shí)驗(yàn)數(shù)據(jù)支撐。Xiang L等結(jié)合高速攝像和動(dòng)態(tài)參數(shù)采集技術(shù)研究了某三葉片誘導(dǎo)輪中發(fā)生的空化不穩(wěn)定現(xiàn)象,清晰揭示了超同步旋轉(zhuǎn)空化工況下空化形態(tài)演變過程,填補(bǔ)了國內(nèi)該領(lǐng)域的可視化試驗(yàn)研究空白。

    上述實(shí)驗(yàn)研究為揭示旋轉(zhuǎn)空化的發(fā)生規(guī)律提供重要的支撐,關(guān)于其發(fā)生機(jī)理許多學(xué)者提出了不同的觀點(diǎn)。Iga等基于二維葉柵的仿真認(rèn)為旋轉(zhuǎn)空化是由潛在的“旋轉(zhuǎn)失速”元素在特定的工況下被激發(fā)引起,但目前為止該解釋尚未得到試驗(yàn)結(jié)果證實(shí)。Kimura等基于數(shù)值模擬發(fā)現(xiàn)泄漏渦與相鄰葉片的前緣相互作用與旋轉(zhuǎn)空化的發(fā)生有密切聯(lián)系。而Tani等數(shù)值研究工作證實(shí)泄漏渦與旋轉(zhuǎn)空化沒有直接關(guān)系,葉尖空化潰滅誘發(fā)的速度散度是引發(fā)旋轉(zhuǎn)空化的原因。

    綜上可知,目前關(guān)于旋轉(zhuǎn)空化的發(fā)生機(jī)理仍存在較大爭(zhēng)議,本文以某二維平板葉柵為研究對(duì)象,利用數(shù)值模擬研究了超同步旋轉(zhuǎn)空化現(xiàn)象,對(duì)超同步旋轉(zhuǎn)空化沿周向傳播的機(jī)理做出了明確解釋。

    1 數(shù)值計(jì)算方法

    1.1 控制方程

    采用均相流模型模擬空化流動(dòng),其假設(shè)液相和氣相為均勻混合的介質(zhì),兩相處于局部平衡狀態(tài),具有相同的速度、壓力等流場(chǎng)信息。利用氣相體積分?jǐn)?shù)或質(zhì)量分?jǐn)?shù)對(duì)兩相進(jìn)行區(qū)分,同時(shí)引入氣相體積分?jǐn)?shù)的輸運(yùn)方程來描述氣液兩相之間的質(zhì)量交換,方程為

    (1)

    (2)

    式中:混合密度=+(1-),、分別為液相、汽相密度;有效黏度=+,為湍流黏度,為混合黏度,=+(1-),、分別為液相、汽相黏度。由于本文研究的介質(zhì)為常溫水,空化過程熱效應(yīng)較弱,可以作等溫處理,因此這里不考慮能量方程。

    1.2 空化模型

    本文采用Singhal模型進(jìn)行空化流動(dòng)的仿真,其具體表達(dá)式為

    (3)

    (4)

    (5)

    1.3 數(shù)值計(jì)算與設(shè)置

    圖1 計(jì)算域Fig. 1 Computation domain

    數(shù)值計(jì)算是基于ANSYS CFX平臺(tái),利用CEL語言將上述Singhal模型嵌入軟件。常規(guī)二方程湍流模型由于高估了空化區(qū)尾緣的黏性,無法準(zhǔn)確預(yù)測(cè)空化脫落等非定?,F(xiàn)象。目前有研究者通過修正湍流模型中的渦黏系數(shù)表達(dá)式,實(shí)現(xiàn)對(duì)空化脫落過程的更準(zhǔn)確預(yù)測(cè)。根據(jù)Byungjin等的研究,誘導(dǎo)輪空化不穩(wěn)定是一種系統(tǒng)不穩(wěn)定現(xiàn)象,與空化脫落這種局部不穩(wěn)定現(xiàn)象沒有直接關(guān)系。這里為了簡(jiǎn)化流場(chǎng)結(jié)構(gòu),提煉經(jīng)典流動(dòng)模型,不考慮空化區(qū)的脫落,選擇SST-進(jìn)行仿真計(jì)算,其具體輸運(yùn)方程和相關(guān)表達(dá)式見文獻(xiàn)[20]。計(jì)算域網(wǎng)格數(shù)量約為2.2×10,具體網(wǎng)格細(xì)節(jié)如圖2所示,非定常計(jì)算時(shí)間步長(zhǎng)為0.000 1 s,根據(jù)文獻(xiàn)[21-23]計(jì)算的經(jīng)驗(yàn),這里采用的網(wǎng)格數(shù)和計(jì)算時(shí)間步長(zhǎng)滿足了無關(guān)性要求。計(jì)算結(jié)果顯示壁面值約20,滿足壁面函數(shù)對(duì)網(wǎng)格尺寸的要求。

    圖2 網(wǎng)格細(xì)節(jié)Fig.2 Mesh details

    2 結(jié)果與分析

    2.1 定常計(jì)算結(jié)果

    圖3為仿真預(yù)測(cè)的空化性能曲線與實(shí)驗(yàn)結(jié)果的對(duì)比,可以看到,仿真預(yù)測(cè)的斷裂點(diǎn)空化數(shù)與實(shí)驗(yàn)結(jié)果比較接近,但是揚(yáng)程變化趨勢(shì)有一定區(qū)別,這是由于簡(jiǎn)化的二維葉柵無法復(fù)現(xiàn)三維誘導(dǎo)輪內(nèi)的復(fù)雜流動(dòng)特性。

    圖3 仿真預(yù)測(cè)的空化性能與實(shí)驗(yàn)結(jié)果對(duì)比Fig.3 Comparison of predicted cavitation performance with experimental results

    圖中標(biāo)注了旋轉(zhuǎn)空化的發(fā)生范圍,仿真預(yù)測(cè)的旋轉(zhuǎn)空化范圍大于實(shí)驗(yàn)結(jié)果,但是二維葉柵的仿真結(jié)果可為揭示超同步旋轉(zhuǎn)空化的發(fā)生機(jī)理提供重要的流場(chǎng)信息。

    從氣相體積分?jǐn)?shù)分布來看,空化數(shù)比較高時(shí),空化區(qū)范圍比較小,而且3個(gè)葉片上空化區(qū)均勻?qū)ΨQ分布;隨空化數(shù)降低,空化區(qū)范圍逐漸擴(kuò)大;在旋轉(zhuǎn)空化發(fā)生范圍內(nèi),3個(gè)葉片上空化區(qū)分布明顯不均勻,下文將詳細(xì)分析;直至3個(gè)葉片上空化區(qū)均發(fā)展至喉部位置,造成流道較嚴(yán)重阻塞,開始發(fā)生空化性能斷裂,整體空化區(qū)發(fā)展規(guī)律與實(shí)驗(yàn)結(jié)果比較一致,因此整體來看,基于二維葉柵的仿真一定程度能夠反映典型的誘導(dǎo)輪內(nèi)空化流動(dòng)發(fā)展規(guī)律。

    2.2 非定常計(jì)算結(jié)果

    對(duì)葉片吸力面上氣相體積分?jǐn)?shù)分布進(jìn)行面積分可以得到空化區(qū)面積,圖4給出了=0.09時(shí)3個(gè)葉片上升力和空化區(qū)面積隨時(shí)間變化關(guān)系??梢钥吹剑S時(shí)間步推進(jìn),升力和空化區(qū)面積均由較小的穩(wěn)定值逐漸擴(kuò)大,最終形成穩(wěn)定的周期性波動(dòng),而且空化區(qū)波動(dòng)明顯沿葉片3-2-1方向傳播,而這正是誘導(dǎo)輪旋轉(zhuǎn)方向。

    圖4 升力和空化區(qū)面積隨時(shí)間步變化關(guān)系(σ=0.09)Fig.4 Lifting force and cavitation area versus timesteps(σ=0.09)

    圖5為實(shí)驗(yàn)結(jié)果捕捉的超同步旋轉(zhuǎn)空化現(xiàn)象,可以看到空化區(qū)波動(dòng)以一定的速度沿周向傳播,如圖5(a)所示,假如關(guān)注每一圈最短的空化區(qū),明顯能看到最短空化區(qū)沿著葉片3-2-1的方向傳播,大約在第6圈時(shí),最短空化區(qū)又回到葉片3,即對(duì)于單個(gè)葉片而言,空化區(qū)波動(dòng)的傳播周期約為5圈。圖5(b)更清晰地反映了該趨勢(shì),可以看到單個(gè)葉片表面空化長(zhǎng)度的波動(dòng)頻率為0.18倍轉(zhuǎn)頻,而且空化區(qū)波動(dòng)的傳播方向與葉輪旋轉(zhuǎn)方向一致,因此從絕對(duì)坐標(biāo)系來看,其特征頻率為1.18倍轉(zhuǎn)頻,詳細(xì)的分析見文獻(xiàn)[13]。

    圖5 超同步旋轉(zhuǎn)空化Fig.5 Super-synchronous rotating cavitation

    對(duì)比實(shí)驗(yàn)結(jié)果,可以得出仿真預(yù)測(cè)的正是超同步旋轉(zhuǎn)空化現(xiàn)象。同時(shí)注意到升力和空化區(qū)波動(dòng)呈明顯的相關(guān)性,即空化區(qū)面積越大,相應(yīng)的升力也越大。升力由葉片表面靜壓值積分所得,因此其絕對(duì)值取決于壓力面和吸力面的壓差,由于空化主要發(fā)生在葉片的吸力面,對(duì)于空化區(qū)面積較大的葉片,其吸力面低壓區(qū)范圍較大,壓差也更大,因此葉片升力較大,下文將重點(diǎn)分析葉片表面空化區(qū)的變化。

    圖6為不同空化數(shù)下3個(gè)葉片無量綱空化區(qū)面積(以葉片吸力面面積為基準(zhǔn))隨時(shí)間變化關(guān)系,可以看到,當(dāng)空化數(shù)較高或較低(=0.18、=0.02)時(shí),空化區(qū)均未出現(xiàn)波動(dòng);當(dāng)空化數(shù)較高時(shí)(=0.18),3個(gè)葉片上空化區(qū)較小,其發(fā)展都不受相鄰葉片的限制,呈現(xiàn)自由發(fā)展的特征;而當(dāng)空化數(shù)降低至一定程度時(shí)(=0.02),3個(gè)葉片上空化區(qū)范圍都較大,受相鄰葉片的限制,難以呈現(xiàn)出沿周向傳播的特征。而當(dāng)空化數(shù)處于一定范圍內(nèi),空化區(qū)面積呈明顯的周期性波動(dòng),對(duì)于單個(gè)葉片,空化區(qū)波動(dòng)頻率分別為28 Hz(=0.15)、26 Hz(=0.14)、22 Hz(=0.09)、9 Hz(=0.05),該結(jié)果也充分表明了空化區(qū)是否與相鄰葉片存在相互作用是旋轉(zhuǎn)空化現(xiàn)象發(fā)生與否的先決條件。

    圖6 不同工況的空化區(qū)波動(dòng)(仿真結(jié)果)Fig.6 Cavitation fluctuation atvarious conditions (simulated results)

    根據(jù)Iga等的研究,空化區(qū)沿周向傳播的頻率比可定義為PVR=(+),其中為空化區(qū)沿葉片旋轉(zhuǎn)方向的傳播速度。上述4個(gè)空化數(shù)下空化區(qū)的周向傳播頻率比PVR分別為1.336、1.312、1.264、1.108,這與旋轉(zhuǎn)空化的特征頻率范圍相符,而且PVR隨著空化數(shù)減小而不斷降低,這也是旋轉(zhuǎn)空化的典型特征,進(jìn)一步證實(shí)了數(shù)值仿真較好地預(yù)測(cè)了二維葉柵內(nèi)發(fā)生的超同步旋轉(zhuǎn)空化現(xiàn)象。

    進(jìn)一步以=0.09計(jì)算結(jié)果為例來揭示空化區(qū)變化規(guī)律,圖7(a)給出了連續(xù)一段時(shí)間內(nèi)氣相體積分?jǐn)?shù)和液流角(flow angle)的云圖,其中間隔時(shí)間Δ=0.048 s。首先從圖中能直觀地看出空化區(qū)波動(dòng)沿葉片3-2-1方向(即葉片旋轉(zhuǎn)方向)傳播。葉片前緣的液流角與空化區(qū)發(fā)展呈現(xiàn)出較好的一致性,即液流角越大,相應(yīng)的空化區(qū)也越大,表明空化區(qū)的增長(zhǎng)或縮短是液流角的變化引起的。圖7 (b)進(jìn)一步給出了空化區(qū)與液流角相互作用關(guān)系,其中液流角取自葉尖前緣的監(jiān)控點(diǎn),相同顏色曲線代表相同葉片,F(xiàn)_B1和V_B1分別表示葉片1前緣液流角和葉片1表面空化區(qū)面積,其他類同。

    圖7 空化區(qū)和液流角的變化(σ=0.09)Fig.7 Variation of cavitation area and flow angle(σ=0.09)

    可以看到,二者均呈現(xiàn)出非常規(guī)律的周期性,但是波動(dòng)曲線形狀有所不同;對(duì)于同一個(gè)葉片(如圖中藍(lán)色虛線框),空化區(qū)的發(fā)展滯后于液流角的變化,即只有當(dāng)液流角在最大值維持一定時(shí)間,空化區(qū)才會(huì)顯著增長(zhǎng);而液流角一旦開始減小,空化區(qū)面積隨即減小;如果關(guān)注兩個(gè)相鄰葉片(如圖中黑色虛線框),可以看到空化區(qū)和液流角的變化呈負(fù)相關(guān),葉片3上空化區(qū)面積從最大值開始減小時(shí),葉片1前緣液流角開始增加;而葉片3空化區(qū)從最小值開始增大時(shí),導(dǎo)致葉片1前緣液流角減小,同時(shí)注意到,只有當(dāng)葉片3空化區(qū)增大至一定程度時(shí)才會(huì)導(dǎo)致葉片1前緣液流角大幅減小,因此空化區(qū)發(fā)展不僅受液流角的影響,而且會(huì)反過來顯著影響相鄰葉片的液流角。

    總結(jié)空化區(qū)面積和液流角相互作用關(guān)系可以得出3個(gè)重要特征:

    1)葉片前緣液流角變化直接影響葉片表面空化區(qū)的發(fā)展,且兩者成正相關(guān);

    2)空化區(qū)變化會(huì)影響相鄰葉片的前緣液流角,且兩者呈負(fù)相關(guān);

    3)只有當(dāng)空化區(qū)增長(zhǎng)至一定程度才會(huì)影響到相鄰葉片的前緣液流角。

    2.3 超同步旋轉(zhuǎn)空化傳播機(jī)理

    為了更清晰揭示空化區(qū)變化對(duì)流場(chǎng)的擾動(dòng),圖8給出不同時(shí)刻瞬態(tài)計(jì)算結(jié)果的差值。

    圖8 不同時(shí)刻空化區(qū)差值和速度矢量差值分布(σ=0.09)Fig.8 Cavitation area difference and velocity vector difference at various moments(σ=0.09)

    例如圖8(a)為圖7中和兩個(gè)時(shí)刻瞬態(tài)計(jì)算結(jié)果的差值,其中負(fù)值代表空化區(qū)減小,正值代表空化區(qū)擴(kuò)大。從速度矢量差值來看,大部分流場(chǎng)速度幾乎不受影響,但是在空化區(qū)尾部,空化區(qū)的變化會(huì)對(duì)流場(chǎng)產(chǎn)生顯著的影響。

    ~時(shí)刻,葉片3上空化區(qū)擴(kuò)大,由于空化區(qū)內(nèi)氣泡的生長(zhǎng),在空化區(qū)尾緣誘發(fā)較大的順時(shí)針渦擾動(dòng)(如圖中白色箭頭),當(dāng)渦擾動(dòng)靠近喉部附近時(shí),引起正的流向速度分量,導(dǎo)致葉片1前緣液流角大幅減小,葉片1上空化區(qū)隨即減小,與此同時(shí)大量氣泡潰滅導(dǎo)致空化區(qū)尾緣形成較大的逆時(shí)針渦擾動(dòng),引起負(fù)的流向速度增量,導(dǎo)致葉片2前緣沖角顯著增大,因此葉片2上空化區(qū)開始增長(zhǎng);~時(shí)刻,由于前緣液流角較大,葉片2上空化區(qū)不斷擴(kuò)大,但是誘發(fā)的順時(shí)針渦擾動(dòng)尚未達(dá)到喉部位置,即不足以影響葉片3前緣的沖角,此時(shí)葉片3上空化區(qū)仍然在增長(zhǎng)直至最大,同時(shí)導(dǎo)致葉片1前緣沖角繼續(xù)減小,相應(yīng)的空化區(qū)也減小,進(jìn)一步導(dǎo)致葉片2前緣沖角增大,促進(jìn)了葉片上空化區(qū)的發(fā)展;~時(shí)刻,葉片2上空化區(qū)進(jìn)一步增長(zhǎng),誘發(fā)的順時(shí)針渦擾動(dòng)達(dá)到喉部位置,導(dǎo)致葉片3前緣沖角大幅減小,因此葉片3上空化區(qū)也顯著減小,并導(dǎo)致葉片1前緣沖角增大,由此葉片1上空化區(qū)開啟了新一輪生長(zhǎng)周期。可以看到,上述過程完成了最大空化區(qū)從葉片3向葉片2的傳播,進(jìn)一步重復(fù)該過程將會(huì)形成沿著葉片旋轉(zhuǎn)方向的周期性傳播。

    進(jìn)一步總結(jié)上述分析,可以得到旋轉(zhuǎn)空化的傳播機(jī)理,如圖9所示。

    圖9 超同步旋轉(zhuǎn)空化傳播機(jī)理Fig.9 Propagation mechanism of super-synchronous rotating cavitation

    假設(shè)在時(shí)刻,3個(gè)葉片表面空化區(qū)呈非對(duì)稱分布,其中葉片3上空化區(qū)最長(zhǎng),葉片1上空化區(qū)最短,葉片2上空化區(qū)剛好發(fā)展至足以影響葉片3的程度,結(jié)合文獻(xiàn)[6]的研究,可認(rèn)為當(dāng)空化區(qū)長(zhǎng)度達(dá)到65%的葉片間距,就會(huì)對(duì)相鄰葉片產(chǎn)生影響,開始出現(xiàn)各類空化不穩(wěn)定現(xiàn)象。葉片2上空化區(qū)持續(xù)增長(zhǎng),其空化區(qū)尾緣誘發(fā)的順時(shí)針渦擾動(dòng)到達(dá)喉部位置,使葉片3前緣沖角大幅減小,相應(yīng)空化區(qū)也大幅減小,同時(shí)在空化區(qū)尾緣誘發(fā)逆時(shí)針渦擾動(dòng),導(dǎo)致葉片1上空化區(qū)開始增大;在時(shí)刻,由于沖角增大,葉片1上空化區(qū)顯著擴(kuò)大,但是其長(zhǎng)度尚未達(dá)到65%的葉片間距,葉片2上空化區(qū)繼續(xù)增長(zhǎng),直至達(dá)到最大;在此過程中,受葉片2上空化區(qū)生長(zhǎng)的影響,葉片3上空化區(qū)持續(xù)減小至最??;在時(shí)刻,葉片1上空化區(qū)開始影響葉片2以后,繼續(xù)生長(zhǎng)至最大;受其影響葉片2上空化區(qū)減小至最小,并進(jìn)一步促進(jìn)了葉片3上空化區(qū)的發(fā)展,由此完成了空化區(qū)波動(dòng)沿葉片旋轉(zhuǎn)方向的傳播。

    上述討論清晰揭示了空化區(qū)與葉片沖角相互作用過程,解釋了空化區(qū)波動(dòng)沿周向傳播的機(jī)制。圖10為旋轉(zhuǎn)空化工況下兩個(gè)相鄰葉片上的空化區(qū)結(jié)構(gòu),可以看到,葉片3幾乎不受葉片2上空化區(qū)的影響,沖角較大(如圖中虛箭頭所示),空化區(qū)范圍也較大,幾乎發(fā)展至葉片1前緣,從葉片1上空化區(qū)的形態(tài)能看出,受上游葉片3空化區(qū)的影響,此時(shí)葉片1前緣的沖角顯著減小,因此實(shí)驗(yàn)結(jié)果一定程度證實(shí)了圖9中傳播機(jī)理的合理性。文獻(xiàn)[12]利用PIV測(cè)量了發(fā)生旋轉(zhuǎn)空化時(shí)誘導(dǎo)輪葉尖的液流角分布,首次通過實(shí)驗(yàn)證實(shí)了較大的空化區(qū)會(huì)導(dǎo)致其下游相鄰葉片沖角為負(fù)值,從而抑制泄漏空化的發(fā)展;而較小的空化區(qū)會(huì)增大下游相鄰葉片的沖角,促進(jìn)葉尖泄漏空化的發(fā)展,進(jìn)一步證實(shí)了圖9中的傳播機(jī)理。

    圖10 超同步旋轉(zhuǎn)工況下空化結(jié)構(gòu)Fig.10 Cavitation structure at SSRC condition

    3 結(jié)論

    本文利用數(shù)值模擬研究了二維葉柵中的超同步旋轉(zhuǎn)空化現(xiàn)象,其中葉柵取自真實(shí)誘導(dǎo)輪,根據(jù)實(shí)驗(yàn)結(jié)果中旋轉(zhuǎn)空化發(fā)生工況的轉(zhuǎn)速、流量換算成來流速度和液流角作為邊界條件,進(jìn)行非定常計(jì)算,得出結(jié)論如下:

    1)仿真預(yù)測(cè)的空化斷裂點(diǎn)與實(shí)驗(yàn)結(jié)果比較接近,且葉柵內(nèi)空化發(fā)展規(guī)律與實(shí)驗(yàn)結(jié)果一致,表明建立的二維葉柵計(jì)算模型能夠較可靠地反映真實(shí)誘導(dǎo)輪內(nèi)空化流動(dòng)特征;

    2)在一定空化數(shù)范圍內(nèi)3個(gè)葉片上空化區(qū)呈現(xiàn)規(guī)律的周期性波動(dòng),波動(dòng)頻率在1.1~1.4之間,且隨空化數(shù)減小而降低,空化區(qū)波動(dòng)傳播方向與誘導(dǎo)輪旋轉(zhuǎn)方向一致,證實(shí)仿真比較準(zhǔn)確地捕捉到超同步旋轉(zhuǎn)空化現(xiàn)象;

    3)計(jì)算了不同時(shí)刻瞬態(tài)流場(chǎng)差值,分析了空化區(qū)變化對(duì)流場(chǎng)的擾動(dòng)。發(fā)現(xiàn)空化區(qū)與葉片前緣沖角成正相關(guān),同時(shí)空化區(qū)擴(kuò)大會(huì)在尾緣誘發(fā)順時(shí)針渦擾動(dòng),空化區(qū)減小會(huì)在尾緣誘發(fā)逆時(shí)針渦擾動(dòng),當(dāng)擾動(dòng)達(dá)到葉片喉部位置時(shí),會(huì)對(duì)相鄰葉片產(chǎn)生影響,因此空化區(qū)擴(kuò)張會(huì)導(dǎo)致相鄰葉片沖角減小,相應(yīng)空化區(qū)也減小,空化區(qū)縮減會(huì)導(dǎo)致相鄰葉片前緣沖角增大,相應(yīng)空化區(qū)也擴(kuò)大,在空化區(qū)和葉片前緣沖角這種相互作用機(jī)制下,空化區(qū)波動(dòng)以一定速度沿周向傳播。

    猜你喜歡
    液流空化前緣
    功率超聲作用下鋼液中空化泡尺寸的演變特性
    鋼鐵釩鈦(2023年5期)2023-11-17 08:48:34
    一種飛機(jī)尾翼前緣除冰套安裝方式
    一種全釩液流電池并網(wǎng)控制系統(tǒng)的仿真設(shè)計(jì)
    三維扭曲水翼空化現(xiàn)象CFD模擬
    不同運(yùn)動(dòng)形式下水物相互作用空化數(shù)值模擬
    H2-Fe3+/Fe2+氧化還原液流電池初步研究
    深水沉積研究進(jìn)展及前緣問題
    前緣
    小說月刊(2015年2期)2015-04-23 08:49:47
    我國—項(xiàng)液流電池國際標(biāo)準(zhǔn)提案獲準(zhǔn)通過
    精品酒店卫生间| 美女福利国产在线| 日韩中文字幕视频在线看片| 美女主播在线视频| 天天躁夜夜躁狠狠久久av| 午夜福利,免费看| 老司机影院成人| 一级片'在线观看视频| 中文天堂在线官网| 久久久国产精品麻豆| 成年美女黄网站色视频大全免费| 黄片无遮挡物在线观看| xxx大片免费视频| 日本av手机在线免费观看| 夜夜骑夜夜射夜夜干| 麻豆乱淫一区二区| 国产日韩欧美在线精品| 黑丝袜美女国产一区| 超色免费av| 美女国产视频在线观看| 日本av手机在线免费观看| 国产亚洲一区二区精品| 亚洲欧美一区二区三区国产| 中文乱码字字幕精品一区二区三区| 亚洲婷婷狠狠爱综合网| 九色亚洲精品在线播放| 一个人免费看片子| 韩国高清视频一区二区三区| 欧美+日韩+精品| 国产精品不卡视频一区二区| 秋霞伦理黄片| 女人高潮潮喷娇喘18禁视频| 如日韩欧美国产精品一区二区三区| 成人亚洲精品一区在线观看| 亚洲av电影在线进入| 人成视频在线观看免费观看| 精品久久久久久电影网| 欧美日韩视频高清一区二区三区二| 热99国产精品久久久久久7| 亚洲av综合色区一区| 在线观看一区二区三区激情| 亚洲中文av在线| 欧美激情极品国产一区二区三区| 亚洲av电影在线进入| av有码第一页| 免费观看av网站的网址| 五月开心婷婷网| 亚洲视频免费观看视频| 欧美+日韩+精品| 免费观看av网站的网址| 精品国产一区二区久久| 9191精品国产免费久久| 欧美日韩亚洲高清精品| 国产又色又爽无遮挡免| 丝袜人妻中文字幕| 男女高潮啪啪啪动态图| 亚洲av免费高清在线观看| 国产亚洲精品第一综合不卡| 午夜福利在线免费观看网站| 中文欧美无线码| 亚洲国产精品一区三区| 99久久精品国产国产毛片| 老汉色∧v一级毛片| 深夜精品福利| 成人毛片a级毛片在线播放| 亚洲av福利一区| 亚洲,一卡二卡三卡| 亚洲国产欧美日韩在线播放| 一区福利在线观看| 中文字幕另类日韩欧美亚洲嫩草| 高清av免费在线| 午夜福利视频精品| 午夜影院在线不卡| 国产男女超爽视频在线观看| 视频区图区小说| 观看美女的网站| a 毛片基地| 久久精品夜色国产| 毛片一级片免费看久久久久| 色哟哟·www| 亚洲欧美成人综合另类久久久| 成人黄色视频免费在线看| 午夜av观看不卡| 国产极品天堂在线| 国产免费福利视频在线观看| 国产女主播在线喷水免费视频网站| 国产成人精品无人区| 欧美精品一区二区大全| 婷婷色麻豆天堂久久| 最近最新中文字幕免费大全7| 少妇的丰满在线观看| 欧美 亚洲 国产 日韩一| 欧美成人午夜免费资源| 久久久久久久久久久免费av| 中文天堂在线官网| 免费av中文字幕在线| 日韩欧美精品免费久久| 国产亚洲精品第一综合不卡| 男女无遮挡免费网站观看| 精品视频人人做人人爽| 久久精品国产亚洲av天美| 纯流量卡能插随身wifi吗| 久久久久久人人人人人| 九九爱精品视频在线观看| 国产又爽黄色视频| 久久久久久久亚洲中文字幕| 久久久久久人妻| 国产在线一区二区三区精| 男女午夜视频在线观看| 久久久国产欧美日韩av| a级片在线免费高清观看视频| 90打野战视频偷拍视频| 亚洲第一区二区三区不卡| 黑丝袜美女国产一区| 亚洲精品视频女| 精品一区二区三区四区五区乱码 | 久久ye,这里只有精品| 日本-黄色视频高清免费观看| 国产av一区二区精品久久| 99国产精品免费福利视频| 亚洲,欧美,日韩| 欧美激情高清一区二区三区 | 成人午夜精彩视频在线观看| 久久狼人影院| 18禁国产床啪视频网站| 捣出白浆h1v1| 久久99蜜桃精品久久| 桃花免费在线播放| 欧美少妇被猛烈插入视频| 国产1区2区3区精品| 国产片特级美女逼逼视频| 久久久久久久大尺度免费视频| 男人添女人高潮全过程视频| 成人二区视频| 精品亚洲成国产av| 国产精品不卡视频一区二区| 亚洲成色77777| 满18在线观看网站| av免费观看日本| 国产高清不卡午夜福利| 亚洲国产日韩一区二区| av.在线天堂| 日韩中字成人| 国产精品三级大全| 少妇猛男粗大的猛烈进出视频| 久久精品夜色国产| 亚洲伊人色综图| 蜜桃在线观看..| 叶爱在线成人免费视频播放| 亚洲av中文av极速乱| 日日撸夜夜添| 亚洲国产精品一区三区| 久久精品国产综合久久久| 中文乱码字字幕精品一区二区三区| 精品国产超薄肉色丝袜足j| 免费黄网站久久成人精品| 国产欧美日韩一区二区三区在线| 国产无遮挡羞羞视频在线观看| 欧美精品亚洲一区二区| 亚洲五月色婷婷综合| 亚洲精品中文字幕在线视频| 精品亚洲成a人片在线观看| 在线观看免费视频网站a站| 国产欧美亚洲国产| 亚洲av成人精品一二三区| 免费日韩欧美在线观看| 两个人免费观看高清视频| 电影成人av| 夫妻午夜视频| 人妻少妇偷人精品九色| 十八禁网站网址无遮挡| 日韩不卡一区二区三区视频在线| 成年人午夜在线观看视频| 日韩电影二区| 18在线观看网站| 久久久久人妻精品一区果冻| 中文字幕亚洲精品专区| 免费高清在线观看日韩| 一级片'在线观看视频| 一本久久精品| 亚洲国产日韩一区二区| 国产97色在线日韩免费| 九色亚洲精品在线播放| 美女脱内裤让男人舔精品视频| 99热网站在线观看| 最近手机中文字幕大全| 考比视频在线观看| 99久国产av精品国产电影| 999久久久国产精品视频| 国产成人精品福利久久| 只有这里有精品99| 久久午夜综合久久蜜桃| 欧美成人精品欧美一级黄| 亚洲,欧美精品.| 亚洲av男天堂| 国产 精品1| 热99久久久久精品小说推荐| 美女主播在线视频| 在现免费观看毛片| 18+在线观看网站| 亚洲av欧美aⅴ国产| 日韩成人av中文字幕在线观看| 另类精品久久| 夫妻性生交免费视频一级片| 午夜久久久在线观看| 久久久久久久大尺度免费视频| 免费大片黄手机在线观看| 欧美日韩一级在线毛片| 天天躁日日躁夜夜躁夜夜| 18禁国产床啪视频网站| 亚洲欧美成人综合另类久久久| 王馨瑶露胸无遮挡在线观看| 蜜桃在线观看..| 成人手机av| 我的亚洲天堂| 黄色 视频免费看| 久久亚洲国产成人精品v| 青草久久国产| 热re99久久国产66热| 999久久久国产精品视频| 老鸭窝网址在线观看| 激情五月婷婷亚洲| 国产精品蜜桃在线观看| 久久久久久久久久人人人人人人| 人人澡人人妻人| 国产日韩欧美亚洲二区| 两性夫妻黄色片| 色婷婷久久久亚洲欧美| 波野结衣二区三区在线| 我的亚洲天堂| 1024香蕉在线观看| 国产精品三级大全| 日韩一区二区三区影片| 狠狠婷婷综合久久久久久88av| 又粗又硬又长又爽又黄的视频| 亚洲欧美一区二区三区国产| 成人免费观看视频高清| 欧美精品国产亚洲| 老汉色∧v一级毛片| 久久99一区二区三区| 97在线视频观看| 国产亚洲一区二区精品| 99国产综合亚洲精品| 美女国产视频在线观看| a级毛片黄视频| 欧美成人午夜免费资源| 日本欧美视频一区| 亚洲国产精品成人久久小说| 一边摸一边做爽爽视频免费| 老熟女久久久| 国产一区二区激情短视频 | 午夜福利在线观看免费完整高清在| 亚洲精品久久久久久婷婷小说| 成人国产av品久久久| 色视频在线一区二区三区| 我要看黄色一级片免费的| 国产乱人偷精品视频| 久久久国产精品麻豆| 亚洲精品久久久久久婷婷小说| 有码 亚洲区| 人妻少妇偷人精品九色| 人妻系列 视频| 少妇精品久久久久久久| 国产精品一区二区在线观看99| av免费观看日本| 国产日韩一区二区三区精品不卡| 性色avwww在线观看| av在线app专区| 亚洲国产日韩一区二区| 欧美日韩成人在线一区二区| 大码成人一级视频| 王馨瑶露胸无遮挡在线观看| 妹子高潮喷水视频| 日本欧美视频一区| 国产精品二区激情视频| 男女下面插进去视频免费观看| 少妇人妻久久综合中文| 久久99热这里只频精品6学生| 韩国av在线不卡| 欧美精品高潮呻吟av久久| 亚洲av日韩在线播放| 亚洲av电影在线观看一区二区三区| 精品国产一区二区三区四区第35| 成人黄色视频免费在线看| 性色avwww在线观看| 日韩欧美一区视频在线观看| 亚洲av男天堂| 超碰97精品在线观看| 久久国产精品大桥未久av| 欧美日韩视频高清一区二区三区二| 国产在线一区二区三区精| 人成视频在线观看免费观看| 亚洲,一卡二卡三卡| 久久久a久久爽久久v久久| 欧美bdsm另类| 一级a爱视频在线免费观看| 黄色一级大片看看| 两个人看的免费小视频| 国产日韩欧美在线精品| 看非洲黑人一级黄片| 久久久精品区二区三区| 欧美日韩一区二区视频在线观看视频在线| 久久综合国产亚洲精品| 久热这里只有精品99| 最黄视频免费看| 最近手机中文字幕大全| 国产在线视频一区二区| √禁漫天堂资源中文www| 国产野战对白在线观看| 五月伊人婷婷丁香| 91精品伊人久久大香线蕉| 80岁老熟妇乱子伦牲交| 亚洲伊人久久精品综合| 五月天丁香电影| 国产探花极品一区二区| 国产成人精品婷婷| 久久狼人影院| 狠狠婷婷综合久久久久久88av| 99久久中文字幕三级久久日本| 精品少妇久久久久久888优播| 五月天丁香电影| 天天躁狠狠躁夜夜躁狠狠躁| 日本-黄色视频高清免费观看| 十分钟在线观看高清视频www| 女性被躁到高潮视频| 国产精品不卡视频一区二区| 搡老乐熟女国产| 极品少妇高潮喷水抽搐| 街头女战士在线观看网站| 国精品久久久久久国模美| 美女午夜性视频免费| 亚洲天堂av无毛| 我要看黄色一级片免费的| 大陆偷拍与自拍| 青春草亚洲视频在线观看| 亚洲伊人色综图| 午夜福利在线免费观看网站| 久久国产亚洲av麻豆专区| 亚洲国产看品久久| 亚洲国产av影院在线观看| 久久久久网色| 女性生殖器流出的白浆| 中文字幕最新亚洲高清| 高清在线视频一区二区三区| 亚洲色图 男人天堂 中文字幕| 2021少妇久久久久久久久久久| 精品人妻熟女毛片av久久网站| 亚洲精品国产一区二区精华液| 超碰97精品在线观看| 高清黄色对白视频在线免费看| 亚洲精品av麻豆狂野| 免费黄频网站在线观看国产| 亚洲综合色惰| 久久久久久久大尺度免费视频| 老鸭窝网址在线观看| 9热在线视频观看99| 久久精品人人爽人人爽视色| 久久综合国产亚洲精品| www.自偷自拍.com| 日韩不卡一区二区三区视频在线| 激情视频va一区二区三区| 精品国产乱码久久久久久男人| 亚洲国产欧美网| 亚洲精品美女久久久久99蜜臀 | 久久久久久久亚洲中文字幕| 久久精品亚洲av国产电影网| 久久久久精品性色| 在线天堂中文资源库| 亚洲色图综合在线观看| 少妇的逼水好多| 国产国语露脸激情在线看| 久久久久久人妻| 2022亚洲国产成人精品| 性色avwww在线观看| 精品国产露脸久久av麻豆| 人人妻人人添人人爽欧美一区卜| 在线观看免费视频网站a站| 有码 亚洲区| 精品亚洲乱码少妇综合久久| 成人黄色视频免费在线看| 欧美亚洲日本最大视频资源| 国产精品 国内视频| 亚洲国产欧美网| 赤兔流量卡办理| 欧美中文综合在线视频| 久久久久久久亚洲中文字幕| 少妇人妻 视频| 伊人久久大香线蕉亚洲五| 免费观看av网站的网址| 黄色毛片三级朝国网站| 成人亚洲精品一区在线观看| 一二三四中文在线观看免费高清| 18禁国产床啪视频网站| 搡老乐熟女国产| 成年女人在线观看亚洲视频| 亚洲av.av天堂| 97在线人人人人妻| 免费播放大片免费观看视频在线观看| 韩国av在线不卡| 两性夫妻黄色片| 国产成人aa在线观看| 国产精品香港三级国产av潘金莲 | 黑人巨大精品欧美一区二区蜜桃| 美女国产视频在线观看| 亚洲美女搞黄在线观看| 日韩中字成人| 亚洲伊人色综图| av免费观看日本| 国产精品久久久久久精品电影小说| 天天躁夜夜躁狠狠久久av| av电影中文网址| av有码第一页| 丰满少妇做爰视频| 五月开心婷婷网| 天天躁日日躁夜夜躁夜夜| 欧美激情高清一区二区三区 | 高清欧美精品videossex| 亚洲男人天堂网一区| 99久久中文字幕三级久久日本| 久久久久久久亚洲中文字幕| 久久这里只有精品19| 中文乱码字字幕精品一区二区三区| 久久久国产精品麻豆| videossex国产| 性高湖久久久久久久久免费观看| 天堂8中文在线网| 一级爰片在线观看| 亚洲av福利一区| 一个人免费看片子| 日韩一卡2卡3卡4卡2021年| 国产无遮挡羞羞视频在线观看| 中文字幕人妻丝袜制服| 美女国产高潮福利片在线看| 你懂的网址亚洲精品在线观看| 国产精品久久久久久av不卡| 黄色配什么色好看| 久久青草综合色| 少妇的逼水好多| 在线观看三级黄色| 久久久国产欧美日韩av| 菩萨蛮人人尽说江南好唐韦庄| 中文字幕制服av| 亚洲人成电影观看| 人妻少妇偷人精品九色| av.在线天堂| 久久人人爽av亚洲精品天堂| 成人午夜精彩视频在线观看| 少妇的丰满在线观看| 晚上一个人看的免费电影| 男女免费视频国产| 精品卡一卡二卡四卡免费| 伦理电影大哥的女人| 观看美女的网站| 一个人免费看片子| 十分钟在线观看高清视频www| 黑人巨大精品欧美一区二区蜜桃| 亚洲欧美成人综合另类久久久| 如何舔出高潮| 大香蕉久久成人网| 9热在线视频观看99| 寂寞人妻少妇视频99o| 热re99久久国产66热| 国产一区二区在线观看av| 9热在线视频观看99| 街头女战士在线观看网站| 久久久久久人妻| 中文字幕色久视频| 免费看不卡的av| 啦啦啦在线免费观看视频4| 一级毛片黄色毛片免费观看视频| 黄色视频在线播放观看不卡| 成人手机av| 亚洲精品一区蜜桃| 国产精品 欧美亚洲| 亚洲国产看品久久| 亚洲成国产人片在线观看| www.自偷自拍.com| 69精品国产乱码久久久| 在线观看www视频免费| 国产黄色免费在线视频| 久久久久国产精品人妻一区二区| 国产女主播在线喷水免费视频网站| 午夜老司机福利剧场| 大香蕉久久网| 亚洲,一卡二卡三卡| 99re6热这里在线精品视频| 制服人妻中文乱码| 午夜福利一区二区在线看| 亚洲av福利一区| av在线app专区| 日韩中字成人| 欧美 日韩 精品 国产| 只有这里有精品99| 免费黄色在线免费观看| a级毛片在线看网站| 日日摸夜夜添夜夜爱| 亚洲国产日韩一区二区| 久久 成人 亚洲| 在线天堂中文资源库| 久久久久久久精品精品| 黄片小视频在线播放| 91久久精品国产一区二区三区| 天天影视国产精品| 女人精品久久久久毛片| 久久久国产欧美日韩av| 国产97色在线日韩免费| 男女午夜视频在线观看| 精品久久久久久电影网| 亚洲图色成人| 免费播放大片免费观看视频在线观看| 999久久久国产精品视频| 国产熟女午夜一区二区三区| 老汉色∧v一级毛片| 成人黄色视频免费在线看| 如何舔出高潮| 日韩不卡一区二区三区视频在线| www.自偷自拍.com| 日日撸夜夜添| 天天躁狠狠躁夜夜躁狠狠躁| 精品午夜福利在线看| 亚洲第一av免费看| 大码成人一级视频| 欧美在线黄色| 一级片'在线观看视频| 天天躁狠狠躁夜夜躁狠狠躁| 中文字幕色久视频| 精品亚洲成a人片在线观看| av在线观看视频网站免费| 99热全是精品| 日本av免费视频播放| 在线亚洲精品国产二区图片欧美| 日韩伦理黄色片| 精品99又大又爽又粗少妇毛片| 国产精品无大码| 99久久人妻综合| 超色免费av| 免费黄色在线免费观看| freevideosex欧美| 在线观看www视频免费| 国产精品无大码| 另类精品久久| 久久久久精品久久久久真实原创| 老熟女久久久| 久久久久视频综合| 国产精品麻豆人妻色哟哟久久| 一级片'在线观看视频| 秋霞伦理黄片| 亚洲内射少妇av| 欧美在线黄色| 欧美国产精品一级二级三级| 久久av网站| 久久青草综合色| 波野结衣二区三区在线| 人人澡人人妻人| 亚洲精品久久成人aⅴ小说| 自拍欧美九色日韩亚洲蝌蚪91| 欧美在线黄色| av一本久久久久| 人人妻人人爽人人添夜夜欢视频| 在线观看国产h片| 欧美+日韩+精品| 久久国产亚洲av麻豆专区| 亚洲经典国产精华液单| 午夜日韩欧美国产| 国产有黄有色有爽视频| 国产老妇伦熟女老妇高清| 久久久久久久久久久久大奶| 亚洲五月色婷婷综合| 成人手机av| 久久久久久人人人人人| 亚洲国产av影院在线观看| 一级爰片在线观看| 男女高潮啪啪啪动态图| 十八禁网站网址无遮挡| 日韩三级伦理在线观看| 大香蕉久久网| 午夜91福利影院| 大陆偷拍与自拍| av不卡在线播放| 伊人亚洲综合成人网| 欧美 亚洲 国产 日韩一| 亚洲人成电影观看| 纯流量卡能插随身wifi吗| 国产97色在线日韩免费| 精品一区二区三卡| 亚洲精品乱久久久久久| 亚洲精品久久成人aⅴ小说| 激情五月婷婷亚洲| 日本黄色日本黄色录像| 午夜日韩欧美国产| 高清黄色对白视频在线免费看| 狂野欧美激情性bbbbbb| 91国产中文字幕| 三上悠亚av全集在线观看| 成人国语在线视频| 日韩一区二区视频免费看| 欧美日韩视频高清一区二区三区二| videosex国产| 国产精品嫩草影院av在线观看| 美女脱内裤让男人舔精品视频| 在线 av 中文字幕| 亚洲精品日韩在线中文字幕| 妹子高潮喷水视频| 国产精品国产三级专区第一集| 久久人人97超碰香蕉20202| 精品久久蜜臀av无| 久久久精品免费免费高清| 久久久久人妻精品一区果冻| 精品少妇黑人巨大在线播放| 啦啦啦在线免费观看视频4| 国产av码专区亚洲av| a 毛片基地| 丝袜人妻中文字幕| 久久久久久人人人人人| 亚洲综合色惰|