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

    小型堆核主泵內(nèi)部流動(dòng)特性數(shù)值計(jì)算

    2024-06-05 00:00:00李天斌郭喜安龍?jiān)?/span>

    摘要: 為研究不同工況下小型堆核主泵內(nèi)部流動(dòng)情況,采用計(jì)算流體動(dòng)力學(xué)(CFD)數(shù)值模擬與試驗(yàn)相結(jié)合的方法,選取4種工況(0.6Qd,0.8Qd,1.0Qd與1.2Qd)進(jìn)行內(nèi)部流動(dòng)特性分析,并選取具有典型意義的出口中心截面,以三維速度流線(xiàn)、速度分布云圖、渦量分布云圖等形式,對(duì)比考察了不同流量工況條件下泵內(nèi)部流動(dòng)規(guī)律及其變化趨勢(shì).通過(guò)分析葉輪與導(dǎo)葉之間的通道回轉(zhuǎn)面壓力、速度分布云圖以及葉輪葉片與導(dǎo)葉葉片的葉片壓力載荷曲線(xiàn),解析了葉輪和導(dǎo)葉內(nèi)部的流動(dòng)分布和能量轉(zhuǎn)換機(jī)制,從而為小型堆核主泵的水力優(yōu)化設(shè)計(jì)提供直觀(guān)認(rèn)識(shí).研究結(jié)果表明:在設(shè)計(jì)流量工況1.0Qd下,小型堆核主泵內(nèi)部流線(xiàn)平順?lè)€(wěn)定,葉片工作面與背面壓力載荷較穩(wěn)定;在小流量工況0.6Qd和0.8Qd下,葉輪葉片上高壓區(qū)增大;在大流量工況1.2Qd運(yùn)行時(shí),泵內(nèi)壓力分布變化較大;試驗(yàn)結(jié)果與數(shù)值計(jì)算結(jié)果的一致性進(jìn)一步驗(yàn)證了計(jì)算模型的準(zhǔn)確性.研究結(jié)果不僅闡釋了小型堆核主泵內(nèi)部的流動(dòng)特性,而且為小型堆核主泵的設(shè)計(jì)提供了一定的理論依據(jù)和應(yīng)用指導(dǎo).

    關(guān)鍵詞: 小型堆核主泵;水動(dòng)力特性;內(nèi)部流動(dòng);數(shù)值計(jì)算;試驗(yàn)

    中圖分類(lèi)號(hào): S277.9;U661.3 文獻(xiàn)標(biāo)志碼: A 文章編號(hào):" 1674-8530(2024)05-0433-07

    DOI:10.3969/j.issn.1674-8530.22.0287

    李天斌,郭喜安,龍?jiān)?小型堆核主泵內(nèi)部流動(dòng)特性數(shù)值計(jì)算[J]. 排灌機(jī)械工程學(xué)報(bào),2024,42(5):433-439,455.

    LI Tianbin, GUO Xi′an, LONG Yun. Numerical calculation of internal flow characteristics of small modular reactors coolant pump[J]. Journal of drainage and irrigation machinery engineering(JDIME), 2024, 42(5): 433-439,455.(in Chinese)

    Numerical calculation of internal flow characteristics of

    small modular reactors coolant pump

    LI Tianbin1, GUO Xi′an2, LONG Yun2*

    (1. Shanghai Electric-KSB Nuclear Pumps and Valves Co.,Ltd., Shanghai 201306, China; 2. National Research Center of Pumps, Jiangsu University, Zhenjiang, Jiangsu 212013, China)

    Abstract:To study the internal flow of small modular reactor (SMR) coolant pumps under different working conditions, a combination method of computational fluid dynamics (CFD) numerical simulation and experimental validation was used. Four operating conditions(0.6Qd, 0.8Qd, 1.0Qd and 1.2Qd)were specifically selected for internal flow analysis. Three-dimensional streamlines were used to represent the flow inside the pump under different conditions. In addition, typical outlet center sections were selected to compare and investigate the internal flow structures and their variations under different flow rates through velocity streamlines, velocity distribution contour maps and vorticity distribution contour maps. Finally, the pressure and velocity distribution contour maps of the channel turning surface between the impeller and diffuser and the blade pressure load curves of the impeller and diffuser were analyzed to elucidate the flow distribution and energy conversion mechanisms within the impeller and diffuser. The results show that under the design flow rate, the internal streamline of the SMR coolant pumps is smooth and stable, and the pressure load on the pressure surface and suction surface of the blades is more stable. Under the partial flow rate of 0.6Qd and 0.8Qd, the high-pressure area on the impeller increases. Under the operation of the flow rate of 1.2Qd, the pressure distribution of the impeller increases significantly. The consistency between the experimental results and the numerical calculation results further confirms the accuracy of the computational model. This research comprehensively elucidates the flow characteristics inside SMR primary pumps, and provides important theoretical basis and practical guidance for their design and development.

    Key words: small modular reactors coolant pump;hydrodynamic characteristics;internal flow;numerical calculations;experiments

    小型核反應(yīng)堆廣泛應(yīng)用于城市熱電聯(lián)供、海水淡化、工業(yè)供氣、海上平臺(tái)能源以及為遠(yuǎn)航艦船提供動(dòng)力等領(lǐng)域[1].國(guó)際原子能機(jī)構(gòu)(IAEA)將小型核電機(jī)組定義為電功率小于300 MW的機(jī)組[2],如美國(guó)一體化堆Nuscale、俄羅斯緊湊型堆KLT-40S、中國(guó)中核集團(tuán)一體化堆ACP100等.先進(jìn)小型核反應(yīng)堆一般采用一體化設(shè)計(jì)及模塊化安裝,大多數(shù)部件在工廠(chǎng)完成組裝[3],對(duì)于其內(nèi)部的穩(wěn)定性具有極高的要求.相比大型核電站,小型核反應(yīng)堆具有較大的成本優(yōu)勢(shì)與極強(qiáng)的適應(yīng)性[4].考慮小型核反應(yīng)堆的一體化、緊湊化及多用途化特點(diǎn),其核主泵通常選用尺寸小和重量輕的高比轉(zhuǎn)數(shù)水泵型式[5].

    小型核反應(yīng)堆及其回路系統(tǒng)中蒸發(fā)器下腔室與核主泵泵殼空間較小,對(duì)葉輪進(jìn)口處的入流產(chǎn)生一定的影響.率志君等[6]采用FFT對(duì)泵蝸殼與葉輪流道內(nèi)的壓力脈動(dòng)進(jìn)行研究,發(fā)現(xiàn)蝸殼和葉輪流道內(nèi)的壓力脈動(dòng)主頻為葉輪葉頻及其倍頻.CHENG等[7]對(duì)AP1000核主泵與蒸汽發(fā)生器下腔室進(jìn)行聯(lián)合計(jì)算,并與均勻來(lái)流計(jì)算結(jié)果進(jìn)行對(duì)比,發(fā)現(xiàn)左邊的主泵揚(yáng)程提高1.10%,效率降低0.86%,而右邊的主泵揚(yáng)程降低2.20%,效率降低1.45%.這表明蒸汽發(fā)生器的下腔室對(duì)核主泵的外特性具有重要影響,甚至有可能改善核主泵的某些水力特性,但是缺乏相應(yīng)的試驗(yàn)驗(yàn)證.龍?jiān)疲?]對(duì)AP1000核主泵全流場(chǎng)特性進(jìn)行數(shù)值計(jì)算,并通過(guò)試驗(yàn)驗(yàn)證了優(yōu)化模型具有較高的水力效率.侯向陶[9]研究了熱管側(cè)蒸汽發(fā)生器下封頭流量分配對(duì)核主泵性能的影響,發(fā)現(xiàn)在下封頭的影響下,核主泵入口流動(dòng)產(chǎn)生畸變,入流畸變使葉輪的軸向力下降10.5%,使葉輪的徑向力增大8.47倍.倪丹等[10]研究了核主泵內(nèi)部非穩(wěn)態(tài)流動(dòng)結(jié)構(gòu)及其壓力脈動(dòng)特性,發(fā)現(xiàn)隨著流量增大,動(dòng)靜干涉作用在靠近出液管的導(dǎo)葉出口處逐漸增強(qiáng),而在小流量工況下,核主泵內(nèi)部流動(dòng)復(fù)雜,在類(lèi)球形蝸殼內(nèi)渦量強(qiáng)度較高,導(dǎo)葉流道堵塞嚴(yán)重.張興等[11]分別對(duì)有分流葉片和無(wú)分流葉片的離心泵模型進(jìn)行數(shù)值計(jì)算,發(fā)現(xiàn)增加分流葉片可提高泵抗空化性能,改善泵內(nèi)壓力脈動(dòng).NI等[12]采用數(shù)值計(jì)算和試驗(yàn)相結(jié)合的方法對(duì)核主泵瞬時(shí)流動(dòng)進(jìn)行研究,分析了渦結(jié)構(gòu)產(chǎn)生的主要位置以及壓力脈動(dòng)的瞬時(shí)特性.龍?jiān)频龋?3]對(duì)核主泵內(nèi)部流動(dòng)特性的研究現(xiàn)狀進(jìn)行綜述,認(rèn)為對(duì)核電站運(yùn)行中事故工況水動(dòng)力特性的研究尤為重要,現(xiàn)有的核主泵全特性曲線(xiàn)方程均是通過(guò)大量試驗(yàn)數(shù)據(jù)進(jìn)行擬合得到,成本較高且數(shù)學(xué)模型精度較低.黃濱等[14]通過(guò)試驗(yàn)對(duì)ACP100核主泵模型進(jìn)行全特性研究,為推廣和發(fā)展中國(guó)小型堆核主泵研制提供了有效數(shù)據(jù)支撐.葉竹等[15]使用APEOS對(duì)模塊化小型核反應(yīng)堆的瞬態(tài)工況進(jìn)行數(shù)值模擬,研究了降功率工況、升功率工況及故障工況的泵轉(zhuǎn)速等參數(shù),但并未對(duì)泵內(nèi)部流動(dòng)特性進(jìn)行研究.劉懿銳等[16]使用STAR-CCM+對(duì)小型反應(yīng)堆上腔室模型進(jìn)行研究,分析了小型反應(yīng)堆上腔室內(nèi)部的混合特性,但并未深入研究小型核反應(yīng)堆主泵內(nèi)部的特性.楊?lèi)?ài)民等[17]對(duì)海洋工況下小型核反應(yīng)堆的水力特性進(jìn)行研究,分析了小型核反應(yīng)堆6個(gè)自由度方向的作用力.

    綜上所述,現(xiàn)階段對(duì)小型堆核主泵內(nèi)部流動(dòng)特性的研究相對(duì)較少.文中采用數(shù)值計(jì)算與試驗(yàn)相結(jié)合的方法,對(duì)不同工況下小型堆核主泵內(nèi)部流動(dòng)特性進(jìn)行分析,揭示核主泵內(nèi)部流場(chǎng)分布和能量轉(zhuǎn)換機(jī)理,并對(duì)渦團(tuán)產(chǎn)生的主要位置進(jìn)行分析,從而為小型堆核主泵的水力設(shè)計(jì)和優(yōu)化提供一定依據(jù).

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

    1.1 計(jì)算模型

    以某一小型堆核主泵模型為研究對(duì)象,其設(shè)計(jì)性能參數(shù)分別為流量Qd=1 500 m3/h,轉(zhuǎn)速n=1 480 r/min,揚(yáng)程H=15 m,效率η=80%,比轉(zhuǎn)數(shù)ns=479.2.根據(jù)小型堆核主泵設(shè)計(jì)性能參數(shù)確定其水力部件基本幾何參數(shù),采用商用CFturbo軟件對(duì)小型堆核主泵的葉輪和導(dǎo)葉等過(guò)流部件進(jìn)行多參數(shù)建模,如圖1所示.小型堆核主泵計(jì)算域主要由直管入口水體、葉輪水體、導(dǎo)葉水體、蝸殼水體與出口段水體等組成.

    1.2 網(wǎng)格劃分

    對(duì)模型進(jìn)行網(wǎng)格劃分,其中葉輪和導(dǎo)葉采用結(jié)構(gòu)化網(wǎng)格,直管入口、蝸殼與出口段采用非結(jié)構(gòu)化網(wǎng)格.因各方案網(wǎng)格劃分均采用統(tǒng)一尺度網(wǎng)格單元,在同一尺度網(wǎng)格下對(duì)不同方案的計(jì)算結(jié)果影響不大[18-19],為了減小多方案計(jì)算周期并節(jié)約計(jì)算成本,文中選取葉輪水體、導(dǎo)葉水體、直管入口水體、蝸殼與出口段水體的網(wǎng)格數(shù)分別為181萬(wàn)、140萬(wàn)、69萬(wàn)、205萬(wàn)(總網(wǎng)格數(shù)約為595萬(wàn))進(jìn)行計(jì)算.

    1.3 數(shù)值計(jì)算方法及邊界條件

    應(yīng)用ANSYS CFX 2020軟件進(jìn)行小型堆核主泵內(nèi)部流動(dòng)計(jì)算,控制方程主要包括連續(xù)方程、動(dòng)量方程及能量方程等.采用SST k-ω湍流模型,該模型具有在壓力梯度大的邊界層內(nèi)計(jì)算較精確的優(yōu)點(diǎn)[20-21],能夠較好地適用于自由剪切流流動(dòng).

    流動(dòng)介質(zhì)為25 ℃水,其密度為997 kg/m3,動(dòng)力黏度為8.899×10-4 kg/(m·s).泵進(jìn)口采用質(zhì)量流量邊界條件,壁面為無(wú)滑移壁面.設(shè)置葉輪為轉(zhuǎn)動(dòng)域,導(dǎo)葉、進(jìn)口段和蝸殼為靜止域,轉(zhuǎn)動(dòng)部件和靜止部件交接面為Frozen Rotor,對(duì)流項(xiàng)采用UpWind格式.

    2 小型堆核主泵內(nèi)部流動(dòng)特性

    雖然核主泵長(zhǎng)期在設(shè)計(jì)工況下運(yùn)行,但工作過(guò)程中往往存在多工況的需求.為了分析小型堆核主泵的性能,揭示小型堆核主泵在小流量和大流量工況下內(nèi)部流動(dòng)的變化,文中選取設(shè)計(jì)工況1.0Qd及其附近的0.6Qd,0.8Qd和1.2Qd工況進(jìn)行計(jì)算.

    2.1 整泵內(nèi)部流動(dòng)

    圖2為不同工況下核主泵內(nèi)部流線(xiàn)分布.

    由圖2可以看出:在不同流量工況下,泵內(nèi)速度流線(xiàn)分布差異較大,隨著流量增大,泵內(nèi)流速增大;在小流量工況下,導(dǎo)葉出口與蝸殼內(nèi)部形成了明顯的低速渦團(tuán);在大流量工況下,導(dǎo)葉出口與蝸殼內(nèi)部流線(xiàn)更為平順;在設(shè)計(jì)流量工況下,由于小型堆核主泵未采用傳統(tǒng)的螺旋形蝸殼,而采用了對(duì)稱(chēng)的類(lèi)球形蝸殼,泵殼內(nèi)沿著周向截面的流速并不是均勻分布,靠近出口處A區(qū)域與B區(qū)域有著較大的速度差,導(dǎo)致泵殼內(nèi)速度流線(xiàn)極為混亂,受導(dǎo)葉內(nèi)出流影響,在泵殼內(nèi)形成大尺度渦團(tuán).這種渦團(tuán)一方面導(dǎo)致流道堵塞,另一方面反向影響導(dǎo)葉內(nèi)的流動(dòng),進(jìn)而對(duì)葉輪內(nèi)的流動(dòng)產(chǎn)生影響.

    2.2 葉輪出口速度流線(xiàn)云圖

    圖3為不同工況下葉輪出口中心截面速度流線(xiàn)分布,可以看出:在不同流量工況下,在導(dǎo)葉相鄰流道流出蝸殼的區(qū)域形成明顯的低速區(qū),低速區(qū)內(nèi)有較為明顯的旋渦,這種低速渦團(tuán)主要受導(dǎo)葉翼型結(jié)構(gòu)和葉輪內(nèi)高速射流共同影響,在小流量工況和大流量工況下均較為明顯,而在設(shè)計(jì)工況附近較弱;蝸殼出口處流線(xiàn)隨流量變化具有明顯的差異,小流量時(shí)流線(xiàn)明顯偏向A區(qū)域,尤其向泵殼壁面聚集,設(shè)計(jì)流量時(shí)流線(xiàn)基本處于分布均勻的狀態(tài),隨著流量增大,流線(xiàn)又向B區(qū)域偏移(見(jiàn)放大部分);在0.6Qd工況下,B區(qū)域存在顯著的回旋渦團(tuán),周?chē)牧鲃?dòng)軌跡受到渦結(jié)構(gòu)的影響而產(chǎn)生扭曲,流動(dòng)的順暢性受到較大的影響.

    2.3 蝸殼內(nèi)渦團(tuán)分布

    圖4為采用q準(zhǔn)則提取的導(dǎo)葉出口中心平面處渦量ω分布,可以看出:在不同流量工況下,各過(guò)流部件都有明顯的渦結(jié)構(gòu),葉輪中渦結(jié)構(gòu)主要分布在葉片前緣,導(dǎo)葉處渦結(jié)構(gòu)主要分布在導(dǎo)葉葉片間的通道內(nèi)與葉片尾緣處;在非設(shè)計(jì)工況下,各過(guò)流部件內(nèi)部渦結(jié)構(gòu)明顯增多,各尺度渦團(tuán)相互作用導(dǎo)致蝸殼內(nèi)部流動(dòng)紊亂;在1.2Qd工況下,蝸殼與出口段的渦結(jié)構(gòu)明顯增多,此時(shí)出口處的流動(dòng)更加紊亂,且B區(qū)域產(chǎn)生的大尺度渦團(tuán)隨流動(dòng)方向逐漸消散;在0.6Qd工況下,導(dǎo)葉內(nèi)部渦團(tuán)數(shù)量較多,嚴(yán)重堵塞流道;在1.0Qd工況下,渦團(tuán)快速消失.

    2.4 葉輪與導(dǎo)葉葉片通道回轉(zhuǎn)面壓力與速度云圖

    葉片通道回轉(zhuǎn)面為葉片通道內(nèi)壁面至外壁面某個(gè)位置展開(kāi)面,由span參數(shù)確定.圖5為span參數(shù)確定的葉片通道回轉(zhuǎn)面示意圖.

    圖6為不同流量工況下span=0.5處葉輪葉片通道回轉(zhuǎn)面壓力云圖,可以看出,隨著流量增大,葉輪出口區(qū)域的壓力逐漸減小.

    圖7為不同流量工況下span=0.5處葉輪葉片通道回轉(zhuǎn)面速度云圖,可以看出:在小流量工況下,葉片頭部存在明顯的低速區(qū),與壓力分布中的高壓區(qū)基本對(duì)應(yīng),在葉片頭部靠近吸力面存在局部高速區(qū);隨著流量增大,液流角度得到優(yōu)化,葉輪流道內(nèi)低速區(qū)明顯減小,流動(dòng)更為平順.

    圖8為不同流量工況下span=0.3處導(dǎo)葉葉片通道回轉(zhuǎn)面速度流線(xiàn)云圖,可以看出:在0.6Qd工況下,不僅在導(dǎo)葉流道末尾處,在導(dǎo)葉流道內(nèi)部也出現(xiàn)了渦結(jié)構(gòu),堵塞流道流動(dòng);隨著流量增大,導(dǎo)葉內(nèi)的低速渦團(tuán)逐漸減弱,流線(xiàn)變得平順.

    2.5 葉輪與導(dǎo)葉葉片壓力載荷分析

    整泵在運(yùn)行過(guò)程中主要承受液動(dòng)載荷與機(jī)械載荷作用,且2種載荷較為復(fù)雜.圖9為葉輪葉片壓力載荷分布曲線(xiàn),圖中橫坐標(biāo)S為沿葉片展向量綱一化位置(0~1.0表示從葉片前緣到尾緣).

    由圖9可以看出:在設(shè)計(jì)流量1.0Qd工況下,葉輪葉片的工作面與背面所受壓力差值較為穩(wěn)定,沒(méi)有較大的波動(dòng);span=0.3處,在小流量0.8Qd工況下S=0~0.2區(qū)間所受壓力差值明顯高于其他工況.span=0.9處,在0.8Qd下壓力差值增高區(qū)延展至S=0.4處;在偏大流量1.2Qd工況下,隨著span增大至0.9,靠近葉片頂端的壓力差值減小,甚至葉輪葉片的壓力載荷數(shù)值出現(xiàn)交匯.葉輪葉片工作面與背面壓力載荷差值越大,葉輪的性能越好,整體上,在設(shè)計(jì)流量工況下,葉輪性能較佳.

    由于葉片載荷與葉片兩側(cè)壓差的作用,工作面上的速度比對(duì)應(yīng)半徑上背面的速度低,而不穩(wěn)定的葉片載荷造成了導(dǎo)葉內(nèi)部流動(dòng)的紊亂,導(dǎo)葉葉片壓力載荷分布如圖10所示.

    由圖10可以看出,小型堆核主泵在非設(shè)計(jì)工況下運(yùn)行時(shí),在S為0.6~0.8位置內(nèi)導(dǎo)葉葉片壓力載荷出現(xiàn)明顯的突然增大的驟變,而在其他位置不同流量工況下的壓力載荷變化趨勢(shì)基本相似.導(dǎo)葉葉片出現(xiàn)驟變的區(qū)域同樣是出現(xiàn)低速渦團(tuán)的區(qū)域,在非設(shè)計(jì)工況下運(yùn)行時(shí),葉片正反面壓力載荷數(shù)值出現(xiàn)交匯的情況.

    3 小型堆核主泵水動(dòng)力性能試驗(yàn)

    試驗(yàn)標(biāo)準(zhǔn)參考GB/T 18149—2000,采用直管進(jìn)口,對(duì)核主泵在不同流量工況下進(jìn)行外特性試驗(yàn),并將試驗(yàn)數(shù)據(jù)與數(shù)值模擬結(jié)果進(jìn)行對(duì)比,如圖11所示.

    由圖11可以看出:在設(shè)計(jì)流量工況附近,數(shù)值計(jì)算和試驗(yàn)所得的外特性曲線(xiàn)變化趨勢(shì)相同,但在小流量工況附近數(shù)值計(jì)算和試驗(yàn)的數(shù)值偏較大;在設(shè)計(jì)流量工況下,數(shù)值計(jì)算揚(yáng)程H=16.0 m,試驗(yàn)揚(yáng)程H=15.5 m,試驗(yàn)結(jié)果滿(mǎn)足設(shè)計(jì)參數(shù)要求,且計(jì)算值高于試驗(yàn)值3.2%,誤差小于5%;在設(shè)計(jì)流量工況下,數(shù)值計(jì)算效率η=82.5%,試驗(yàn)效率η=81.2%,試驗(yàn)結(jié)果滿(mǎn)足設(shè)計(jì)參數(shù)要求,且數(shù)值計(jì)算效率高于試驗(yàn)效率1.3%.

    通過(guò)數(shù)值計(jì)算與試驗(yàn)對(duì)比,驗(yàn)證了數(shù)值計(jì)算的可靠性.同時(shí)分析兩者誤差產(chǎn)生的原因,主要是數(shù)值計(jì)算沒(méi)有考慮管路與泵內(nèi)壁的粗糙度,忽略了沿程水力損失與壁面效應(yīng).

    4 結(jié) 論

    通過(guò)對(duì)不同工況下核主泵內(nèi)部流場(chǎng)的數(shù)值模擬,獲得了不同工況下整泵內(nèi)部流動(dòng)情況,并進(jìn)行試驗(yàn)驗(yàn)證,得到結(jié)論如下:

    1) 通常核主泵蝸殼采用球形或類(lèi)球形設(shè)計(jì),該結(jié)構(gòu)承受應(yīng)力較為均勻,但這種結(jié)構(gòu)蝸殼使環(huán)向流速不均勻,將引起二次流,造成額外的流動(dòng)損失,同時(shí),蝸殼內(nèi)的流動(dòng)會(huì)反向干涉導(dǎo)葉和葉輪內(nèi)的流動(dòng),致使葉輪和導(dǎo)葉內(nèi)的流動(dòng)不穩(wěn)定,誘導(dǎo)更為嚴(yán)重和復(fù)雜的流致振動(dòng).

    2) 不同流量工況對(duì)小型堆核主泵內(nèi)部的流動(dòng)特性影響較大,在偏離設(shè)計(jì)工況下運(yùn)行所產(chǎn)生的影響往往是負(fù)面的.在設(shè)計(jì)流量工況下,小型堆核主泵內(nèi)部流線(xiàn)平順?lè)€(wěn)定,內(nèi)部渦結(jié)構(gòu)少且尺度小,葉片工作面與背面壓力載荷較穩(wěn)定.在小流量工況下,核主泵內(nèi)部渦結(jié)構(gòu)強(qiáng)度明顯增強(qiáng),在0.6Qd和0.8Qd小流量工況運(yùn)行時(shí),葉輪葉片上高壓區(qū)增大,導(dǎo)葉葉片間出現(xiàn)大量的渦結(jié)構(gòu),導(dǎo)致導(dǎo)葉流道堵塞,導(dǎo)葉葉片壓力由于渦結(jié)構(gòu)的影響出現(xiàn)驟變,蝸殼內(nèi)出現(xiàn)不同尺度的渦結(jié)構(gòu).在1.2Qd流量工況運(yùn)行時(shí),葉輪葉片上低壓區(qū)增大,壓力分布變化較大,導(dǎo)葉葉片尾緣處旋渦結(jié)構(gòu)呈現(xiàn)的形式與小流量相反,葉片上壓力同樣也會(huì)出現(xiàn)驟變.

    3) 對(duì)比核主泵外特性數(shù)值計(jì)算結(jié)果和試驗(yàn)數(shù)據(jù),在設(shè)計(jì)工況附近,兩者誤差較??;但在偏小流量工況,誤差較大.試驗(yàn)驗(yàn)證了數(shù)值計(jì)算結(jié)果的可靠性.

    參考文獻(xiàn)(References)

    [1] 張馨玉,郭慧芳,袁永龍. 全球小型模塊化反應(yīng)堆發(fā)展綜述[C]//中國(guó)核學(xué)會(huì)2021年學(xué)術(shù)年會(huì)論文集第8冊(cè)(核情報(bào)分卷),2021:240-247.

    [2] 孫小凱,張明,李林蔚. 先進(jìn)小型壓水堆的發(fā)展現(xiàn)狀及建議[C] //中國(guó)核學(xué)會(huì)2021年學(xué)術(shù)年會(huì)論文集第2冊(cè)(核能動(dòng)力分卷),2021:152-155.

    [3] 謝仕君,汪為慶,朱勃,等. 小型核反應(yīng)堆核主泵斜流泵水力部件.中國(guó):CN216642550U[P]. 2022-05-31.

    [4] 劉展,戚展飛,王國(guó)棟,等.國(guó)內(nèi)外小堆專(zhuān)設(shè)配置和CAP200專(zhuān)設(shè)論證[J].核科學(xué)與工程,2020,40(3):468-474.

    LlU Zhan,QI Zhanfei, WANG Guodong, et al. Small reactors engineered safety feature(ESF) and CAP200 passive ESF demonstration[J].Nuclear science and engineering, 2020,40(3):468-474. (in Chinese)

    [5] 張麗平,蔡龍.國(guó)內(nèi)小型壓水堆核電站核主泵技術(shù)[J].水泵技術(shù),2018 (5): 1-4.

    ZHANG Liping, CAI Long. Nuclear main pump techno-logy for a small-sized pressurized water reactor nuclear power station[J]. Pump technology, 2018 (5): 1-4. (in Chinese)

    [6] 率志君,張權(quán),陳春來(lái),等.多級(jí)離心泵整機(jī)流場(chǎng)三維非穩(wěn)態(tài)湍流壓力脈動(dòng)特性分析[J].哈爾濱工程大學(xué)學(xué)報(bào),2013, 34(3): 306-311.

    SHUAI Zhijun, ZHANG Quan, CHEN Chunlai, et al. Characteristic analysis of three-dimensional unsteady turbulent pressure fluctuation in whole fow field of multi-stage centrifugal pump[J]. Journal of Harbin Engineering University, 2013, 34(3): 306-311. (in Chinese)

    [7] CHENG H, LI H, YIN J, et al. Investigation of the distortion suction flow on the performance of the canned nuclear coolant pump[C]//Proceedings of the ISFMFE Conference, 2014.

    [8] 龍?jiān)? AP1000核主泵優(yōu)化設(shè)計(jì)及全特性數(shù)值模擬[D].鎮(zhèn)江:江蘇大學(xué),2014.

    [9] 侯向陶. 蒸汽發(fā)生器下封頭對(duì)核主泵性能影響研究[D].杭州: 浙江大學(xué), 2016.

    [10] 倪丹,楊敏官,高波,等.混流式核主泵內(nèi)流動(dòng)結(jié)構(gòu)與壓力脈動(dòng)特性關(guān)聯(lián)分析[J].工程熱物理學(xué)報(bào),2017,38(8):1676-1682.

    NI Dan, YANG Minguan, GAO Bo,et al. The internal correlations between unsteady flow and pressure pulsa-tions in a nuclear reactor coolant pump[J]. Journal of engineering thermophysics, 2017,38(8):1676-1682. (in Chinese)

    [11] 張興,賴(lài)喜德,廖姣.分流葉片對(duì)離心泵非定??栈鲃?dòng)特性影響分析[J]. 熱能動(dòng)力工程, 2017, 32(4): 42-43.

    ZHANG Xing, LAI Xide, LIAO Jiao. Influence of splitter blades on unsteady cavitating flow characteristics of centrifugal pump[J]. Journal of engineering for thermal energy and power, 2017, 32(4): 42-43.(in Chinese)

    [12] NI D, YANG M G, GAO B, et al. Experimental and numerical investigation on the pressure pulsation and instantaneous flow structure in a nuclear reactor coolant pump[J]. Nuclear engineering and design, 2018, 337:261-270.

    [13] 龍?jiān)疲瑝燮?,朱榮生,等.核主泵內(nèi)部流動(dòng)研究現(xiàn)狀與技術(shù)發(fā)展綜述[J].排灌機(jī)械工程學(xué)報(bào),2020,38(11):1081-1097.

    LONG Yun, YUAN Shouqi, ZHU Rongsheng, et al. Review on research status of internal flow and technolo-gical development of reactor coolant pump[J]. Journal of drainage and irrigation machinery engineering, 2020,38(11):1081-1097. (in Chinese)

    [14] 黃濱,武鵬,劉祥松,等.ACP100核主泵四象限全特性模型試驗(yàn)研究[J].流體機(jī)械,2020,48(2):8-11.

    HUANG Bin, WU Peng, LIU Xiangsong, et al. Four-quadrant full characteristic model test of the ACP100 reactor coolant pump[J]. Fluid machinery,2020,48(2):8-11. (in Chinese)

    [15] 葉竹,劉明皓,蔡志云,等.模塊化小型反應(yīng)堆給水系統(tǒng)仿真分析[J].能源與節(jié)能,2021(2):211-214.

    YE Zhu, LIU Minghao, CAI Zhiyun, et al. Simulation analysis of water supply system of modular small reactor[J]. Energy and energy conservation, 2021(2):211-214. (in Chinese)

    [16] 劉懿銳,董秀臣,張?chǎng)?,?小型反應(yīng)堆上腔室混合特性數(shù)值模擬[J].軟件,2021,42(10):111-113.

    LIU Yirui, DONG Xiuchen, ZHANG Xin, et al. Numerical simulation of mixing characteristics in upper chamber of small reactor[J]. Software, 2021,42(10):111-113. (in Chinese)

    [17] 楊?lèi)?ài)民,宋厚德,劉曉晶.海洋工況下海洋小堆熱工水力特性分析[J].核科學(xué)與工程,2022,42(1):34-43.

    YANG Aimin, SONG Houde,LIU Xiaojing. Analysis of thermal and hydraulic characteristics of a small marine reactor under ocean conditions[J]. Nuclear science and engineering, 2022,42(1):34-43. (in Chinese)

    [18] 池寶濤,朱玉麒,郭前建,等. 面向雙層插值邊界面法的非結(jié)構(gòu)自動(dòng)網(wǎng)格劃分[J]. 湖南大學(xué)學(xué)報(bào)(自然科學(xué)版),2023,50(4):76-86.

    CHI Baotao,ZHU Yuqi,GUO Qianjian,et al. Adaptive unstructured mesh generation in the dual interpolation boundary face method implementation[J]. Journal of Hunan University (natural sciences),2023,50(4):76-86. (in Chinese)

    [19] 王曉丁,奉策強(qiáng),郭正榮,等. 非能動(dòng)余熱排出系統(tǒng)空冷器設(shè)計(jì)研究[J].化工設(shè)備與管道,2023,60(4):48-54.

    WANG Xiaoding, FENG Ceqiang, GUO Zhengrong, et al. Design and research of air cooler for passive residual heat removal system[J]. Process equipment amp; piping,2023,60(4):48-54. (in Chinese)

    [20] 于群,章宇,王陽(yáng),等. 大流量工況下離心泵非定??栈鲃?dòng)特性分析[J].機(jī)電工程,2022,39(7): 935-941.

    YU Qun,ZHANG Yu,WANG Yang,et al. Analysis on unsteady cavitation flow characteristics of centrifugal pump at large flow rate[J]. Journal of mechanical & electrical engineering,2022,39(7): 935-941.(in Chinese)

    [21] 張寧,鄭凡坤,劉曉凱,等.離心泵偏工況下流動(dòng)結(jié)構(gòu)分布及演化特性研究[J].流體機(jī)械, 2023, 51 (8): 27-32.

    ZHANG Ning, ZHENG Fankun, LIU Xiaokai,et al. Distribution and evolution characteristics of complex flow field in a low specific speed centrifugal pump[J]. Fluid machinery, 2023, 51 (8): 27-32. (in Chinese)

    (責(zé)任編輯 陳建華)

    收稿日期: 2022-11-18; 修回日期: 2023-01-20; 網(wǎng)絡(luò)出版時(shí)間: 2024-04-25

    網(wǎng)絡(luò)出版地址: https://link.cnki.net/urlid/32.1814.TH.20240423.1113.016

    基金項(xiàng)目: 國(guó)家自然科學(xué)基金資助項(xiàng)目(51906085,U20A20292);中國(guó)博士后科學(xué)基金資助項(xiàng)目(2019M651734)

    第一作者簡(jiǎn)介: 李天斌(1979—),男,山東東阿人,高級(jí)工程師( Tianbin.li@sec-ksb.com),主要從事核主泵設(shè)計(jì)研究.

    通信作者簡(jiǎn)介: 龍?jiān)疲?988—),男,江蘇宿遷人,副研究員,博士(longyun@ujs.edu.cn),主要從事流體機(jī)械內(nèi)部流動(dòng)研究.

    国产亚洲欧美98| 午夜a级毛片| 88av欧美| 91成年电影在线观看| 日韩高清综合在线| 色尼玛亚洲综合影院| 中文在线观看免费www的网站 | 十八禁网站免费在线| 99热只有精品国产| 国模一区二区三区四区视频 | 久久久久免费精品人妻一区二区| 久久精品人妻少妇| 欧美性猛交黑人性爽| 香蕉国产在线看| 国产久久久一区二区三区| 欧美绝顶高潮抽搐喷水| 国产一区在线观看成人免费| 少妇裸体淫交视频免费看高清 | 亚洲中文字幕日韩| 亚洲精品久久成人aⅴ小说| 亚洲性夜色夜夜综合| 久99久视频精品免费| 精品第一国产精品| 亚洲人成网站在线播放欧美日韩| 黄色女人牲交| 久久精品91无色码中文字幕| 国产久久久一区二区三区| 欧美激情久久久久久爽电影| 亚洲国产精品sss在线观看| 最近视频中文字幕2019在线8| 亚洲人成电影免费在线| 真人做人爱边吃奶动态| 婷婷精品国产亚洲av在线| 88av欧美| 俺也久久电影网| 欧美成人免费av一区二区三区| 亚洲专区字幕在线| 美女午夜性视频免费| 精品国产乱码久久久久久男人| 亚洲人与动物交配视频| 国产亚洲精品综合一区在线观看 | 久久久国产成人精品二区| 在线国产一区二区在线| 大型黄色视频在线免费观看| 久久久久久久午夜电影| xxxwww97欧美| 亚洲欧美日韩东京热| 午夜福利18| 嫩草影院精品99| 国产精品美女特级片免费视频播放器 | 校园春色视频在线观看| 国产精品1区2区在线观看.| 国产不卡一卡二| 亚洲精品中文字幕一二三四区| 欧美丝袜亚洲另类 | 好男人电影高清在线观看| 哪里可以看免费的av片| 无人区码免费观看不卡| 一本综合久久免费| 天天躁狠狠躁夜夜躁狠狠躁| 香蕉丝袜av| 欧美色视频一区免费| 国内精品久久久久精免费| 国产一级毛片七仙女欲春2| 国产精品免费一区二区三区在线| 两性夫妻黄色片| 精品久久久久久久久久免费视频| 人妻夜夜爽99麻豆av| 国产91精品成人一区二区三区| 又大又爽又粗| 国产99久久九九免费精品| 麻豆久久精品国产亚洲av| 国产高清视频在线观看网站| 日本免费一区二区三区高清不卡| videosex国产| 两人在一起打扑克的视频| 国产亚洲精品综合一区在线观看 | √禁漫天堂资源中文www| 亚洲色图av天堂| av欧美777| 久久天躁狠狠躁夜夜2o2o| 国产久久久一区二区三区| 亚洲专区国产一区二区| 18禁国产床啪视频网站| 午夜免费激情av| 美女扒开内裤让男人捅视频| 色老头精品视频在线观看| 亚洲国产看品久久| 免费人成视频x8x8入口观看| 国产又黄又爽又无遮挡在线| 久久久久久人人人人人| 少妇裸体淫交视频免费看高清 | 又粗又爽又猛毛片免费看| 久久精品综合一区二区三区| 欧美黄色淫秽网站| 亚洲精品一卡2卡三卡4卡5卡| 国产av一区二区精品久久| 18美女黄网站色大片免费观看| 日韩欧美国产在线观看| 日本免费一区二区三区高清不卡| 露出奶头的视频| 在线看三级毛片| 无限看片的www在线观看| 亚洲国产精品999在线| 两个人视频免费观看高清| 成人三级做爰电影| 三级国产精品欧美在线观看 | 别揉我奶头~嗯~啊~动态视频| 国产精品电影一区二区三区| 国产男靠女视频免费网站| 亚洲九九香蕉| 国产高清videossex| 夜夜看夜夜爽夜夜摸| 免费在线观看完整版高清| 午夜精品在线福利| 脱女人内裤的视频| 国产爱豆传媒在线观看 | 欧美日韩国产亚洲二区| 麻豆久久精品国产亚洲av| 国产三级黄色录像| 熟女少妇亚洲综合色aaa.| 国产精品免费视频内射| av有码第一页| 色综合欧美亚洲国产小说| 麻豆成人av在线观看| 香蕉av资源在线| 制服诱惑二区| 超碰成人久久| 深夜精品福利| 亚洲人成电影免费在线| 白带黄色成豆腐渣| 精品乱码久久久久久99久播| 美女 人体艺术 gogo| 国产精品亚洲美女久久久| av欧美777| 国产亚洲av嫩草精品影院| 国产1区2区3区精品| 午夜福利18| 国产精品 国内视频| 黑人欧美特级aaaaaa片| 午夜福利在线在线| 免费看美女性在线毛片视频| 熟女少妇亚洲综合色aaa.| 在线a可以看的网站| 99在线人妻在线中文字幕| av片东京热男人的天堂| 日韩中文字幕欧美一区二区| 久久久久久国产a免费观看| 亚洲免费av在线视频| 国语自产精品视频在线第100页| 宅男免费午夜| 久久久久久亚洲精品国产蜜桃av| 变态另类成人亚洲欧美熟女| 久久精品综合一区二区三区| 国产v大片淫在线免费观看| 亚洲欧美精品综合一区二区三区| 老汉色av国产亚洲站长工具| 国产v大片淫在线免费观看| 女同久久另类99精品国产91| a级毛片a级免费在线| 国产单亲对白刺激| 国产一级毛片七仙女欲春2| 国产精品自产拍在线观看55亚洲| 黄色丝袜av网址大全| 我的老师免费观看完整版| 国产精品爽爽va在线观看网站| 久久草成人影院| 女同久久另类99精品国产91| 91成年电影在线观看| 一卡2卡三卡四卡精品乱码亚洲| 成人特级黄色片久久久久久久| 床上黄色一级片| 欧美丝袜亚洲另类 | 午夜a级毛片| 真人一进一出gif抽搐免费| 色综合婷婷激情| 色噜噜av男人的天堂激情| 欧美成狂野欧美在线观看| 国产亚洲精品av在线| 一a级毛片在线观看| 99国产精品一区二区蜜桃av| 日日夜夜操网爽| 97人妻精品一区二区三区麻豆| 国模一区二区三区四区视频 | 老熟妇仑乱视频hdxx| 欧美日韩乱码在线| 中文字幕久久专区| 午夜a级毛片| 99精品久久久久人妻精品| 一本大道久久a久久精品| 黄片小视频在线播放| 中文字幕精品亚洲无线码一区| 久久久久国产精品人妻aⅴ院| 91麻豆精品激情在线观看国产| 女警被强在线播放| 小说图片视频综合网站| 国产伦人伦偷精品视频| 成人欧美大片| 在线永久观看黄色视频| 国产97色在线日韩免费| 999精品在线视频| 亚洲精品中文字幕一二三四区| 国产在线精品亚洲第一网站| 国产一区在线观看成人免费| 久久香蕉激情| 国语自产精品视频在线第100页| 亚洲国产欧美人成| xxx96com| 午夜福利欧美成人| 狂野欧美白嫩少妇大欣赏| 亚洲国产精品合色在线| 舔av片在线| 法律面前人人平等表现在哪些方面| 久久婷婷成人综合色麻豆| 国产亚洲精品久久久久5区| 两个人看的免费小视频| 亚洲国产看品久久| 在线观看免费午夜福利视频| 久久国产精品影院| av视频在线观看入口| 日日夜夜操网爽| 国产精品香港三级国产av潘金莲| 在线看三级毛片| 中文亚洲av片在线观看爽| 极品教师在线免费播放| 国产探花在线观看一区二区| 一进一出好大好爽视频| 欧美极品一区二区三区四区| 午夜福利高清视频| 最近最新免费中文字幕在线| 亚洲av中文字字幕乱码综合| 亚洲成人国产一区在线观看| 91麻豆av在线| 日本免费一区二区三区高清不卡| 老鸭窝网址在线观看| 一区二区三区国产精品乱码| 亚洲av美国av| 国产精品久久电影中文字幕| 国产精品久久久久久亚洲av鲁大| 两性夫妻黄色片| 欧美日韩瑟瑟在线播放| 国产精品av久久久久免费| 两个人看的免费小视频| 国产精品久久久人人做人人爽| 欧美色欧美亚洲另类二区| 成人特级黄色片久久久久久久| 午夜亚洲福利在线播放| 午夜日韩欧美国产| www.999成人在线观看| 特大巨黑吊av在线直播| 精品人妻1区二区| a在线观看视频网站| 欧美日韩国产亚洲二区| 90打野战视频偷拍视频| 日本一本二区三区精品| 国产精品香港三级国产av潘金莲| 又爽又黄无遮挡网站| 欧美一区二区精品小视频在线| 亚洲五月天丁香| 好看av亚洲va欧美ⅴa在| 亚洲熟妇中文字幕五十中出| 中文字幕最新亚洲高清| 一区福利在线观看| 成人三级做爰电影| 亚洲国产精品999在线| 在线观看免费午夜福利视频| 老汉色∧v一级毛片| 欧美在线一区亚洲| 中文字幕av在线有码专区| 国产私拍福利视频在线观看| 国产成人啪精品午夜网站| 人人妻人人澡欧美一区二区| 日韩有码中文字幕| 变态另类成人亚洲欧美熟女| 免费人成视频x8x8入口观看| 特级一级黄色大片| 叶爱在线成人免费视频播放| 99久久99久久久精品蜜桃| 国产精品一区二区三区四区久久| 国产成人啪精品午夜网站| 亚洲乱码一区二区免费版| 在线观看日韩欧美| 国产精品久久久久久久电影 | 嫩草影院精品99| 怎么达到女性高潮| 国产精品一区二区三区四区久久| 欧美久久黑人一区二区| 久久中文看片网| 国产成人啪精品午夜网站| 久久人人精品亚洲av| av福利片在线| 国产精品久久久久久久电影 | 舔av片在线| 免费观看精品视频网站| 亚洲男人天堂网一区| 亚洲精品国产一区二区精华液| 人人妻,人人澡人人爽秒播| 亚洲成av人片免费观看| 黄色成人免费大全| 国产成人精品久久二区二区91| 色综合婷婷激情| 哪里可以看免费的av片| 怎么达到女性高潮| 久久精品影院6| 正在播放国产对白刺激| 精品国产美女av久久久久小说| 精品欧美一区二区三区在线| 国产精品久久久av美女十八| 可以免费在线观看a视频的电影网站| 国产三级中文精品| 午夜老司机福利片| 精品电影一区二区在线| av福利片在线| 天天躁夜夜躁狠狠躁躁| 日本一区二区免费在线视频| 麻豆久久精品国产亚洲av| 一本久久中文字幕| 国产精品亚洲一级av第二区| 久久精品成人免费网站| 三级男女做爰猛烈吃奶摸视频| 亚洲一码二码三码区别大吗| 色老头精品视频在线观看| 免费看美女性在线毛片视频| 成人永久免费在线观看视频| 黄频高清免费视频| 99国产精品一区二区三区| 亚洲狠狠婷婷综合久久图片| 看黄色毛片网站| 国产av不卡久久| 亚洲在线自拍视频| 长腿黑丝高跟| 国产av麻豆久久久久久久| 88av欧美| 免费观看精品视频网站| 长腿黑丝高跟| 黄色毛片三级朝国网站| 黄色视频不卡| 国产1区2区3区精品| 99在线人妻在线中文字幕| 国产成人av激情在线播放| 久久中文看片网| 正在播放国产对白刺激| 美女黄网站色视频| 两个人看的免费小视频| 中文字幕人妻丝袜一区二区| 精品久久久久久久久久免费视频| 亚洲av成人av| 日日爽夜夜爽网站| 90打野战视频偷拍视频| 国产熟女xx| 成在线人永久免费视频| 无限看片的www在线观看| av在线天堂中文字幕| 日韩欧美在线乱码| 亚洲九九香蕉| 国产精品乱码一区二三区的特点| 人成视频在线观看免费观看| 国产精品一区二区精品视频观看| 欧美成人一区二区免费高清观看 | 免费av毛片视频| svipshipincom国产片| av在线播放免费不卡| 久久草成人影院| 97人妻精品一区二区三区麻豆| 女人被狂操c到高潮| 亚洲成人精品中文字幕电影| 中出人妻视频一区二区| 美女午夜性视频免费| 婷婷丁香在线五月| 国产精品香港三级国产av潘金莲| 国内久久婷婷六月综合欲色啪| 男女下面进入的视频免费午夜| 精品欧美一区二区三区在线| 亚洲精品久久成人aⅴ小说| 99热这里只有精品一区 | 亚洲精品美女久久久久99蜜臀| 久久久久久久久免费视频了| 岛国视频午夜一区免费看| 久久热在线av| 国产黄a三级三级三级人| 啦啦啦韩国在线观看视频| 亚洲熟妇中文字幕五十中出| 麻豆一二三区av精品| 亚洲中文av在线| 亚洲中文日韩欧美视频| 日韩欧美免费精品| 欧美不卡视频在线免费观看 | 一进一出抽搐gif免费好疼| av福利片在线| 亚洲真实伦在线观看| 两人在一起打扑克的视频| 亚洲自拍偷在线| 免费看日本二区| 99re在线观看精品视频| 国产片内射在线| 国产精品,欧美在线| 久久热在线av| av有码第一页| 国产亚洲精品一区二区www| 91国产中文字幕| 亚洲 欧美 日韩 在线 免费| 欧美日韩福利视频一区二区| 日韩精品青青久久久久久| 在线看三级毛片| 无遮挡黄片免费观看| 亚洲第一欧美日韩一区二区三区| 午夜a级毛片| 国产片内射在线| 免费在线观看完整版高清| 欧美一级a爱片免费观看看 | 老司机在亚洲福利影院| 男人的好看免费观看在线视频 | 成人国语在线视频| 午夜福利欧美成人| 一本精品99久久精品77| 亚洲精品美女久久久久99蜜臀| 丁香六月欧美| 香蕉丝袜av| 又黄又爽又免费观看的视频| 麻豆av在线久日| 欧美成人午夜精品| 操出白浆在线播放| 亚洲欧美日韩无卡精品| 国产精品永久免费网站| 美女高潮喷水抽搐中文字幕| 精品高清国产在线一区| 亚洲欧美一区二区三区黑人| 久久久久国产精品人妻aⅴ院| 中亚洲国语对白在线视频| 成人特级黄色片久久久久久久| 久9热在线精品视频| 日韩免费av在线播放| 国产日本99.免费观看| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品永久免费网站| 久久精品影院6| 亚洲国产欧美网| 精品免费久久久久久久清纯| av免费在线观看网站| 青草久久国产| 久久久久亚洲av毛片大全| 亚洲免费av在线视频| 国产又色又爽无遮挡免费看| xxx96com| 一个人免费在线观看的高清视频| 欧美 亚洲 国产 日韩一| 亚洲国产欧美一区二区综合| 身体一侧抽搐| e午夜精品久久久久久久| 美女大奶头视频| 国产激情久久老熟女| 亚洲,欧美精品.| 国产99久久九九免费精品| 天天添夜夜摸| 老司机在亚洲福利影院| 国内少妇人妻偷人精品xxx网站 | 国产成人精品久久二区二区免费| 制服丝袜大香蕉在线| 色老头精品视频在线观看| 久久久久九九精品影院| 日韩免费av在线播放| 久久久国产成人精品二区| 两个人免费观看高清视频| 精品久久久久久久久久免费视频| 嫩草影院精品99| 日韩精品中文字幕看吧| 嫁个100分男人电影在线观看| 国产爱豆传媒在线观看 | 免费电影在线观看免费观看| 久久 成人 亚洲| 久久人妻福利社区极品人妻图片| 国产黄片美女视频| 别揉我奶头~嗯~啊~动态视频| 最好的美女福利视频网| 丰满人妻熟妇乱又伦精品不卡| 亚洲av成人av| 老司机深夜福利视频在线观看| 色精品久久人妻99蜜桃| 久久久久久九九精品二区国产 | 久久精品国产清高在天天线| 成人精品一区二区免费| av视频在线观看入口| 成人特级黄色片久久久久久久| 无人区码免费观看不卡| 日韩精品免费视频一区二区三区| 欧美久久黑人一区二区| 欧美性猛交╳xxx乱大交人| 激情在线观看视频在线高清| 亚洲欧美日韩高清专用| 亚洲国产精品999在线| 日韩欧美国产在线观看| 一个人免费在线观看的高清视频| 国产乱人伦免费视频| 最近在线观看免费完整版| 亚洲乱码一区二区免费版| 高清在线国产一区| 亚洲无线在线观看| 国产精品一区二区三区四区免费观看 | avwww免费| 亚洲欧美日韩高清在线视频| 久久人妻福利社区极品人妻图片| 国产视频一区二区在线看| 热99re8久久精品国产| 在线观看一区二区三区| 夜夜躁狠狠躁天天躁| 极品教师在线免费播放| 一边摸一边做爽爽视频免费| 免费在线观看成人毛片| 色综合亚洲欧美另类图片| 真人做人爱边吃奶动态| 亚洲黑人精品在线| 午夜福利高清视频| 精品久久久久久成人av| 黄色视频不卡| 婷婷亚洲欧美| 免费在线观看黄色视频的| 国产又色又爽无遮挡免费看| 亚洲 欧美一区二区三区| 制服丝袜大香蕉在线| 女同久久另类99精品国产91| aaaaa片日本免费| 精品国产美女av久久久久小说| 波多野结衣巨乳人妻| 69av精品久久久久久| 99久久久亚洲精品蜜臀av| 午夜福利成人在线免费观看| 在线免费观看的www视频| 国产成人精品久久二区二区免费| 夜夜爽天天搞| 又粗又爽又猛毛片免费看| 蜜桃久久精品国产亚洲av| 亚洲精品国产精品久久久不卡| 欧美性长视频在线观看| 最近在线观看免费完整版| 18禁黄网站禁片免费观看直播| 精品第一国产精品| 黄频高清免费视频| 黄色 视频免费看| 老司机福利观看| 日韩三级视频一区二区三区| 亚洲午夜精品一区,二区,三区| 别揉我奶头~嗯~啊~动态视频| 美女午夜性视频免费| 成人一区二区视频在线观看| 亚洲aⅴ乱码一区二区在线播放 | 俺也久久电影网| 亚洲精品一区av在线观看| 伦理电影免费视频| 国产精品av久久久久免费| 成人特级黄色片久久久久久久| 亚洲一区高清亚洲精品| 亚洲,欧美精品.| 一本大道久久a久久精品| 欧美成狂野欧美在线观看| 久久国产乱子伦精品免费另类| 精品一区二区三区av网在线观看| 免费观看人在逋| 色哟哟哟哟哟哟| 精品久久蜜臀av无| 长腿黑丝高跟| 国语自产精品视频在线第100页| 狂野欧美白嫩少妇大欣赏| 丰满的人妻完整版| 香蕉久久夜色| 久久性视频一级片| 成人三级黄色视频| 日日夜夜操网爽| av超薄肉色丝袜交足视频| 亚洲精品在线观看二区| 制服人妻中文乱码| 国产黄色小视频在线观看| 老汉色∧v一级毛片| 美女扒开内裤让男人捅视频| 国产成年人精品一区二区| 成人精品一区二区免费| 国产精品av久久久久免费| 免费在线观看视频国产中文字幕亚洲| 国产黄a三级三级三级人| 日日干狠狠操夜夜爽| 国产精品国产高清国产av| 俄罗斯特黄特色一大片| 欧美精品啪啪一区二区三区| a级毛片在线看网站| 又紧又爽又黄一区二区| 日本一二三区视频观看| 久热爱精品视频在线9| 亚洲成av人片免费观看| 亚洲全国av大片| 狂野欧美激情性xxxx| 日韩精品免费视频一区二区三区| 观看免费一级毛片| 久久久久久大精品| 亚洲最大成人中文| 十八禁人妻一区二区| 日韩欧美 国产精品| 久久久久久亚洲精品国产蜜桃av| 在线观看一区二区三区| 麻豆久久精品国产亚洲av| 久久久久久久午夜电影| 99国产精品一区二区三区| 国产精品久久视频播放| 黄色 视频免费看| 欧美色欧美亚洲另类二区| 国内揄拍国产精品人妻在线| 久久久久久久久免费视频了| 视频区欧美日本亚洲| 天天一区二区日本电影三级| 免费看十八禁软件| 日韩欧美 国产精品| av片东京热男人的天堂| 日本 av在线| 中文字幕久久专区| 久久中文字幕人妻熟女| 白带黄色成豆腐渣| 好男人在线观看高清免费视频|