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

    ACP100集成式堆頂結(jié)構(gòu)可壓縮氣體對(duì)流傳熱數(shù)值模擬

    2016-11-03 07:54:43何培峰馬梓淇孫善文周進(jìn)雄
    核技術(shù) 2016年10期
    關(guān)鍵詞:控制棒壁面流場(chǎng)

    何培峰 許 斌 羅 英 余 豪 馬梓淇 孫善文 周進(jìn)雄

    1(中國(guó)核動(dòng)力研究設(shè)計(jì)院 核反應(yīng)堆系統(tǒng)設(shè)計(jì)技術(shù)重點(diǎn)實(shí)驗(yàn)室 成都 610200)2(西安交通大學(xué) 航天航空學(xué)院 機(jī)械結(jié)構(gòu)強(qiáng)度與振動(dòng)國(guó)家重點(diǎn)實(shí)驗(yàn)室 西安 710049)

    ACP100集成式堆頂結(jié)構(gòu)可壓縮氣體對(duì)流傳熱數(shù)值模擬

    何培峰1許斌1羅英1余豪2馬梓淇2孫善文2周進(jìn)雄2

    1(中國(guó)核動(dòng)力研究設(shè)計(jì)院核反應(yīng)堆系統(tǒng)設(shè)計(jì)技術(shù)重點(diǎn)實(shí)驗(yàn)室成都610200)2(西安交通大學(xué)航天航空學(xué)院機(jī)械結(jié)構(gòu)強(qiáng)度與振動(dòng)國(guó)家重點(diǎn)實(shí)驗(yàn)室西安710049)

    反應(yīng)堆集成式堆頂結(jié)構(gòu)的功能之一是通過(guò)冷卻氣體的對(duì)流傳熱為控制棒驅(qū)動(dòng)機(jī)構(gòu)提供冷卻。針對(duì)ACP100反應(yīng)堆集成式堆頂結(jié)構(gòu)建立完整的模型并劃分流場(chǎng)網(wǎng)格,基于ANSYS/CFX軟件分別對(duì)可壓縮氣體和不可壓縮氣體進(jìn)行計(jì)算,并嚴(yán)格對(duì)比分析了其結(jié)果,重點(diǎn)研究氣體可壓縮性對(duì)對(duì)流傳熱計(jì)算結(jié)果的影響。計(jì)算結(jié)果表明,氣體可壓縮性對(duì)速度、溫度分布和壓降有較大影響,忽略氣體可壓縮性所預(yù)測(cè)的控制棒驅(qū)動(dòng)機(jī)構(gòu)表面最高溫度偏低,控制棒驅(qū)動(dòng)機(jī)構(gòu)間隙氣體最大流速和壓降也偏低。

    可壓縮氣體,集成式堆頂結(jié)構(gòu),計(jì)算流體動(dòng)力學(xué),對(duì)流傳熱

    隨著數(shù)值仿真計(jì)算技術(shù)的進(jìn)步,計(jì)算流體動(dòng)力學(xué)(Computational fluid dynamic, CFD)方法已經(jīng)被廣泛應(yīng)用于各種流場(chǎng)問(wèn)題的數(shù)值仿真計(jì)算[1?3],反應(yīng)堆的流場(chǎng)分析問(wèn)題就是其中之一。Martinez等[4]基于STAR-CD對(duì)反應(yīng)堆一回路熱工水力進(jìn)行了數(shù)值模擬。Kao等[5]對(duì)反應(yīng)堆壓力容器上腔室流場(chǎng)進(jìn)行了分析。Baliga等[6]討論了集成式堆頂結(jié)構(gòu)在商用壓水反應(yīng)堆結(jié)構(gòu)上的應(yīng)用,論述了堆頂結(jié)構(gòu)內(nèi)部控制棒驅(qū)動(dòng)機(jī)構(gòu)(Control rod drive mechanism, CRDM)冷卻系統(tǒng)的設(shè)計(jì)和基本參數(shù)。余志偉等[7]對(duì)M310堆頂結(jié)構(gòu)通風(fēng)冷卻系統(tǒng)進(jìn)行了不可壓縮氣體流場(chǎng)和溫度場(chǎng)數(shù)值分析。于浩等[8]基于可壓縮理想氣體對(duì)堆頂結(jié)構(gòu)風(fēng)冷系統(tǒng)進(jìn)行了數(shù)值模擬。

    在對(duì)反應(yīng)堆堆頂結(jié)構(gòu)進(jìn)行對(duì)流傳熱數(shù)值模擬時(shí),氣體的可壓縮性是一個(gè)基本問(wèn)題。對(duì)于氣體可壓縮性問(wèn)題,已有相關(guān)學(xué)者進(jìn)行了研究[9?10]。通常情況下認(rèn)為當(dāng)氣體流動(dòng)速度不高、壓力變化較小時(shí),可以近似為不可壓縮氣體;但是當(dāng)氣體流動(dòng)馬赫數(shù)大于0.3 Ma時(shí),必然會(huì)引起很大的壓力變化,從而導(dǎo)致密度也會(huì)發(fā)生顯著的變化,這時(shí)就需要考慮氣體可壓縮性的影響。由于堆頂結(jié)構(gòu)復(fù)雜,局部區(qū)域氣流馬赫數(shù)要高于0.3 Ma,因此在對(duì)堆頂結(jié)構(gòu)進(jìn)行流場(chǎng)分析時(shí),氣體的可壓縮性會(huì)對(duì)堆頂結(jié)構(gòu)流場(chǎng)分析產(chǎn)生較大影響?,F(xiàn)有關(guān)于堆頂對(duì)流傳熱模擬的文獻(xiàn)中,既有采用可壓縮假設(shè)的,也有采用不可壓縮假設(shè)的,更缺乏對(duì)兩種不同假設(shè)計(jì)算結(jié)果的定量對(duì)比和分析。基于此,本文以ACP100反應(yīng)堆集成式堆頂結(jié)構(gòu)為例,深入研究了氣體可壓縮性對(duì)集成式堆頂結(jié)構(gòu)對(duì)流傳熱特性的影響。首先利用HyperMesh軟件對(duì)堆頂結(jié)構(gòu)進(jìn)行高質(zhì)量網(wǎng)格劃分;然后基于ANSYS/CFX15.0軟件,分別分析在氣體可壓縮與不可壓縮條件下堆頂結(jié)構(gòu)的速度場(chǎng)、溫度場(chǎng)以及壓力場(chǎng)。

    1 計(jì)算模型與求解設(shè)置

    1.1幾何模型

    ACP100集成式堆頂結(jié)構(gòu)通風(fēng)冷卻系統(tǒng)主要由保溫層、風(fēng)板、控制棒驅(qū)動(dòng)機(jī)構(gòu)、圍筒、出口、支撐板、抗震板、固定保溫層等組成。其工作原理是冷卻空氣首先通過(guò)頂部的進(jìn)風(fēng)口進(jìn)入,流經(jīng)控制棒驅(qū)動(dòng)機(jī)構(gòu)間隙,之后通過(guò)圍筒下部和圍板,最后經(jīng)出風(fēng)口流出,從而達(dá)到為控制棒驅(qū)動(dòng)機(jī)構(gòu)降溫的目的,其幾何模型如圖1所示。

    圖1 ACP100堆頂結(jié)構(gòu)幾何模型Fig.1 ACP100 3D model.

    1.2網(wǎng)格劃分

    由于堆頂結(jié)構(gòu)比較復(fù)雜,使用四面體網(wǎng)格自動(dòng)劃分方法雖然可以快速地對(duì)堆頂結(jié)構(gòu)進(jìn)行網(wǎng)格劃分,但是這種方法劃分的網(wǎng)格質(zhì)量不高,而且需要的網(wǎng)格數(shù)量較多。本文基于HyperMesh有限元?jiǎng)澐周浖?duì)堆頂結(jié)構(gòu)進(jìn)行了高質(zhì)量的六面體網(wǎng)格劃分,相比于四面體網(wǎng)格而言,六面體網(wǎng)格在同樣的網(wǎng)格尺寸下數(shù)量要少很多并且其網(wǎng)格方向更能迎合流場(chǎng)方向(例如在邊界層附近區(qū)域),因此采用六面體網(wǎng)格在計(jì)算分析時(shí)更容易收斂,結(jié)果更加精確。由于網(wǎng)格數(shù)量在很大程度上會(huì)影響計(jì)算結(jié)果,因此我們對(duì)堆頂結(jié)構(gòu)分別劃分了6.48×106、9.42×106、1.319×107三種不同數(shù)量網(wǎng)格。表1給出了三種不同數(shù)量網(wǎng)格所計(jì)算的堆頂結(jié)構(gòu)進(jìn)出口壓差以及出口速度對(duì)比,從表1中結(jié)果可以看出,當(dāng)網(wǎng)格數(shù)量超過(guò)9.42×106之后,繼續(xù)增加網(wǎng)格數(shù)量計(jì)算結(jié)果基本保持不變。由此可知網(wǎng)格數(shù)量為9.42×106滿足網(wǎng)格無(wú)關(guān)性要求,因此對(duì)堆頂結(jié)構(gòu)的計(jì)算分析采用網(wǎng)格數(shù)量為9.42×106的網(wǎng)格。堆頂結(jié)構(gòu)整體流場(chǎng)以及頂部流場(chǎng)網(wǎng)格如圖2所示。

    表1 不同網(wǎng)格數(shù)量分析Table 1 Analysis for different amounts of mesh.

    圖2 ACP100堆頂結(jié)構(gòu)流場(chǎng)網(wǎng)格 (a) 整體流場(chǎng)網(wǎng)格,(b) 頂部流場(chǎng)網(wǎng)格Fig.2 ACP100 mesh model. (a) Overall flow field mesh, (b) Top flow field mesh

    1.3邊界條件與求解設(shè)置

    ACP100集成式堆頂結(jié)構(gòu)的工作環(huán)境為反應(yīng)堆安全殼內(nèi)部,流體介質(zhì)為空氣,通常在運(yùn)行工況下,堆頂結(jié)構(gòu)入口附近的氣體壓力為一個(gè)大氣壓,入口溫度在5?80 oC之間,堆頂結(jié)構(gòu)出口冷卻風(fēng)量的質(zhì)量流量在9.5458 kg·s?1左右。我們基于ANSYS/CFX軟件分別對(duì)不同工況下堆頂結(jié)構(gòu)流場(chǎng)在氣體可壓縮與不可壓縮條件下進(jìn)行分析計(jì)算,邊界條件與求解設(shè)置如下:

    1) 流體域設(shè)置:在可壓縮條件下,采用可壓縮的理想氣體模型作為實(shí)際空氣的近似,流體介質(zhì)選擇理想氣體,熱量傳輸模型選擇全熱模型。參考?jí)毫υO(shè)置為1.01×105Pa,湍流模型選擇K-Epsilon湍流模型。在不可壓縮條件下,采用不可壓縮的氣體模型作為實(shí)際空氣的近似,流體介質(zhì)選擇25 oC空氣,熱量傳輸模型選擇熱焓模型,其他設(shè)置與可壓縮條件相同。

    2) 入口及出口邊界條件設(shè)置:入口邊界類型選擇開(kāi)放式邊界,質(zhì)量與動(dòng)量參量設(shè)置為開(kāi)放壓強(qiáng)與方向,相對(duì)壓力選擇0 Pa,方向選擇與邊界垂直,湍流參量選擇中等湍流密度,熱量傳輸設(shè)定為開(kāi)放式溫度,環(huán)境溫度根據(jù)工況的不同分別設(shè)為5 oC、20 oC、55 oC、80 oC。堆頂結(jié)構(gòu)的出口邊界類型設(shè)置為出口,出口流體的質(zhì)量與動(dòng)量設(shè)置為質(zhì)量流量速率,其根據(jù)工況的不同設(shè)為9.5458 kg·s?1,以及在此基礎(chǔ)上分別增加和減少10%。

    3) 壁面邊界條件設(shè)置:控制棒驅(qū)動(dòng)機(jī)構(gòu)壁面設(shè)置為不可滑移的光滑壁面,熱流量為3592.1 W·m?2,絕熱壁面設(shè)置為不可滑移的光滑絕熱壁面,保溫層恒溫壁面設(shè)置為不可滑移的光滑壁面,溫度設(shè)置為60 oC。

    4) 求解步長(zhǎng)設(shè)置:由于堆頂結(jié)構(gòu)網(wǎng)格數(shù)量巨大,所以合理的步長(zhǎng)選擇對(duì)于求解精度以及收斂性至關(guān)重要。步長(zhǎng)過(guò)小,收斂速度慢,但是比較穩(wěn)定;步長(zhǎng)過(guò)大,收斂速度快,但是容易造成震蕩甚至發(fā)散。經(jīng)過(guò)不斷試算,最終得到了合適的求解步長(zhǎng),如表2所示,在該求解步長(zhǎng)的計(jì)算下,質(zhì)量與動(dòng)量、熱傳導(dǎo)、湍流耗散與湍流動(dòng)能均收斂到5×10?5左右。

    表2 求解步長(zhǎng)設(shè)定Table 2 Analysis step length.

    2 可壓縮與不可壓縮氣體計(jì)算結(jié)果對(duì)比

    2.1環(huán)境溫度為20 oC、冷卻風(fēng)量為9.5458 kg·s-1時(shí)結(jié)果分析

    圖3分別給出了在氣體可壓與不可壓條件下堆頂結(jié)構(gòu)整體流線圖,從圖3中可以看出,冷卻空氣從堆頂入口開(kāi)始進(jìn)入,在經(jīng)過(guò)控制棒驅(qū)動(dòng)機(jī)構(gòu)間隙時(shí),由于流動(dòng)區(qū)域變小,氣流流速增大。氣流在通過(guò)控制棒驅(qū)動(dòng)機(jī)構(gòu)間隙之后由于流動(dòng)空間變大,氣流流速減小,在通過(guò)圍板之后氣流流速增大,最終通過(guò)出口流出。從可壓與不可壓條件下的流線圖均可以看出氣流在圍筒底部區(qū)域以及出口與圍板的連接位置附近存在較多漩渦,這會(huì)導(dǎo)致流場(chǎng)在這些位置產(chǎn)生較大壓力損失。在可壓條件下,空氣最大流速為151.7 m·s?1,發(fā)生在圍板靠近出口區(qū)域。在不可壓條件下,空氣最大流速為134.1 m·s?1,相比可壓條件減小較大,但出現(xiàn)位置基本相同。

    圖3 堆頂結(jié)構(gòu)整體流線圖?。╝) 可壓縮,(b) 不可壓縮Fig.3 ACP100 streamline. (a) Compressible, (b) Incompressible

    圖4分別給出了在氣體可壓與不可壓條件下控制棒驅(qū)動(dòng)機(jī)構(gòu)壁面溫度云圖。從圖4(a)可以看出,在可壓條件下,控制棒驅(qū)動(dòng)機(jī)構(gòu)壁面最高溫度為669.8K,其中外圍控制棒驅(qū)動(dòng)機(jī)構(gòu)壁面溫度較低,靠近中心區(qū)域的溫度較高。氣體在流經(jīng)磁軛線圈時(shí)由于間隙較小,流速較大,因此磁軛線圈壁面溫度較低,在兩層磁軛線圈中間區(qū)域,由于流動(dòng)空間突然變大,會(huì)在該區(qū)域產(chǎn)生漩渦,且流速不高,因此在磁軛線圈中間區(qū)域壁面溫度較高。在流過(guò)下層磁軛線圈后,流動(dòng)空間增大會(huì)導(dǎo)致氣體流速的降低,所以下層磁軛線圈下部壁面溫度較高。在不可壓條件下,控制棒驅(qū)動(dòng)機(jī)構(gòu)壁面最高溫度為536.9 K,相比可壓條件下最高溫度有所降低,同時(shí)整體溫度分布也更加平均。

    圖4 控制棒驅(qū)動(dòng)機(jī)構(gòu)壁面溫度云圖?。╝) 可壓縮,(b)不可壓縮Fig.4 CRDM wall temperature. (a) Compressible, (b) Incompressible

    圖5給出了Y=0 m的XZ平面位置示意圖以及在該平面上的速度以及壓力分布云圖。從圖5(b)、(c)可以看出,在可壓條件下,從上層磁軛線圈間隙到下層磁軛線圈間隙,氣流流速逐漸增高,大部分速度高于20 m·s?1,在下層磁軛線圈間隙氣流速度最大,達(dá)到30.96 m·s?1。氣流在流過(guò)下層磁軛線圈間隙之后,由于流動(dòng)空間增大導(dǎo)致速度逐漸降低,受底部外側(cè)負(fù)壓影響,氣流流向偏向控制棒驅(qū)動(dòng)機(jī)構(gòu)外側(cè)。在不可壓條件下,其速度分布規(guī)律與可壓條件基本相同,氣體最高流速為30.06 m·s?1,低于可壓條件下的30.96 m·s?1,在下層磁軛線圈底部,氣流流向偏向控制棒驅(qū)動(dòng)機(jī)構(gòu)外側(cè)趨勢(shì)沒(méi)有可壓縮氣體明顯。從圖5(d)、(e)可以看出,在可壓條件下,沿著控制棒驅(qū)動(dòng)機(jī)構(gòu)向下,其壓力逐漸降低,由上部的最大值?6.6 Pa減小到底部最小值?710 Pa,壓降為703.4 Pa。在兩層磁軛線圈中間區(qū)域,由于流動(dòng)空間的突然增大會(huì)出現(xiàn)漩渦,因此在這些區(qū)域會(huì)導(dǎo)致很大的壓力損失。在不可壓條件下,其壓力相比可壓條件有所增大,在該平面上壓力最大值為?5.5 Pa,最小值為?681.7 Pa,壓降為676.2 Pa,低于可壓條件下的703.4 Pa。

    圖5 XZ平面速度及壓力分布云圖(a) XZ平面位置示意圖,(b) 可壓縮速度矢量圖,(c) 不可壓縮速度矢量圖,(d) 可壓縮壓力云圖,(e) 不可壓縮壓力云圖Fig.5 Velocity and pressure distribution in XZ plane. (a) XZ plane position, (b) Velocity of compressible, (c) Velocity of incompressible, (d) Pressure of compressible, (e) Pressure of incompressible

    圖6給出了控制棒驅(qū)動(dòng)機(jī)構(gòu)三個(gè)線圈平面位置示意圖及壓力分布云圖。從圖6中可以看出,在可壓條件下,上層線圈平均壓力在?330 Pa左右,壓力較大位置出現(xiàn)在靠近圍板附近區(qū)域;中層線圈平均壓力在?445 Pa左右,壓力較大位置也出現(xiàn)在靠近圍板附近區(qū)域;下層線圈平均壓力在?615 Pa左右,壓力較大區(qū)域出現(xiàn)在線圈平面中間位置。在不可壓條件下,不同線圈平面平均壓力相比可壓條件偏小,其中上、中、下三層線圈平均壓力分別在?317Pa、?463 Pa、?594 Pa左右,上層線圈平面內(nèi)最大壓力同樣出現(xiàn)在靠近圍板附近區(qū)域;下層線圈壓力分布規(guī)律則與可壓條件不同,壓力最大位置沒(méi)有出現(xiàn)在線圈平面中心區(qū)域。

    圖6 控制棒驅(qū)動(dòng)機(jī)構(gòu)線圈平面壓力云圖?。╝) 線圈平面位置示意圖,(b) 可壓縮提升線圈平面,(c) 不可壓縮提升線圈平面,(d) 可壓縮傳遞線圈平面,(e) 不可壓縮傳遞線圈,(f) 可壓縮夾持線圈平面,(g) 不可壓縮夾持線圈平面Fig.6 Pressure distribution in coil plane. (a) Coil plane position, (b) Lift of compressible, (c) Lift of incompressible, (d) Movable of compressible, (e) Movable of incompressible, (f) Stationary of compressible, (g) Stationary of incompressible

    2.2不同環(huán)境溫度與不同冷卻風(fēng)量結(jié)果分析

    通常情況下,堆頂結(jié)構(gòu)工作的環(huán)境溫度為5?80oC之間,因此我們分別討論了環(huán)境溫度為5oC、20 oC、55 oC、80 oC下氣體可壓縮性對(duì)堆頂結(jié)構(gòu)對(duì)流傳熱特性的影響。圖7給出了在可壓與不可壓條件下控制棒驅(qū)動(dòng)機(jī)構(gòu)壁面平均溫度以及間隙平均速度,由圖7(a)可知,隨著環(huán)境溫度的升高,控制棒驅(qū)動(dòng)機(jī)構(gòu)壁面平均溫度也隨之升高且氣體的可壓縮性對(duì)控制棒驅(qū)動(dòng)機(jī)構(gòu)壁面平均溫度幾乎沒(méi)有影響。由圖 7(b)可知,在20 oC、55 oC、80 oC三種環(huán)境溫度下,可壓條件下的控制棒驅(qū)動(dòng)機(jī)構(gòu)間隙平均速度相比不可壓條件大,而5 oC環(huán)境溫度下,可壓條件下的平均速度相比不可壓小。

    圖7 控制棒驅(qū)動(dòng)機(jī)構(gòu)壁面平均溫度(a)和間隙平均速度(b)Fig.7 Average temperature in CRDM wall (a) and average velocity in CRDM gap (b).

    在原方案的基礎(chǔ)上減小10%冷卻風(fēng)量以及增加10%冷卻風(fēng)量,其他條件不變,分別在氣體可壓與不可壓的條件下進(jìn)行計(jì)算分析,以此來(lái)考察氣體的可壓縮性對(duì)堆頂結(jié)構(gòu)對(duì)流傳熱特性的影響。圖8給出了在可壓與不可壓條件下中層線圈平面速度與出口平均速度。從圖8(a)中可以看出,隨著冷卻風(fēng)量的增多,中層線圈平面的平均速度也隨之增大,氣體的可壓縮性對(duì)中層線圈平面平均速度影響不大。從圖8(b)可以看出,隨著冷卻風(fēng)量的增多,出口平均速度隨之增大,可壓條件下出口平均速度大于不可壓。

    圖8 傳遞線圈平面(a)和出口(b)平均速度Fig.8 Average velocity in movable gripper coil plane (a) and outlet (b).

    3 結(jié)語(yǔ)

    本文以ACP100集成式堆頂結(jié)構(gòu)為例,基于ANSYS/CFX軟件,對(duì)氣體可壓縮性對(duì)堆頂結(jié)構(gòu)對(duì)流傳熱特性的影響進(jìn)行了深入研究。分別計(jì)算分析了可壓縮與不可壓縮條件下堆頂結(jié)構(gòu)的速度場(chǎng)、溫度場(chǎng)、壓力場(chǎng)分布特性,進(jìn)而對(duì)比了不同環(huán)境溫度和不同冷卻風(fēng)量下堆頂結(jié)構(gòu)的流場(chǎng)分布特性。計(jì)算結(jié)果表明,氣體可壓縮性對(duì)速度、溫度分布和壓降有較大影響,忽略氣體可壓縮性所預(yù)測(cè)的控制棒驅(qū)動(dòng)結(jié)構(gòu)表面最高溫度偏低,控制棒驅(qū)動(dòng)機(jī)構(gòu)間隙氣體最大流速和壓降也偏低。因此在進(jìn)行堆頂結(jié)構(gòu)流場(chǎng)分析時(shí)需要將氣體的可壓縮性考慮在內(nèi),這樣得到的仿真結(jié)果更加接近真實(shí)情況。

    1 Beukelmann D, Guo W, Holzer W, et al. Safety assessment of reactor pressure vessel integrity for loss of coolant accident conditions[J]. Journal of Pressure Vessel Technology, 2012, 134(1): 345?355. DOI: 10.1115/ 1.4004799

    2 葛志浩, 彭勇升, 呂逸君, 等. 液態(tài)金屬在堆芯子通道內(nèi)的湍流換熱[J]. 核技術(shù), 2015, 38(9): 090603. DOI: 10.11889/j.0253-3219.2015.hjs.38.090603

    GE Zhihao, PENG Yongsheng, LYU Yijun, et al. Turbulent heat transfer of liquid metal inside the sub-channels of reactor core[J]. Nuclear Techniques, 2015, 38(9): 090603. DOI: 10.11889/j.0253-3219.2015. hjs.38.090603

    3 曹寅, 吳燕華, 林超, 等. 渦腔倒角結(jié)構(gòu)對(duì)渦流二極管性能影響的數(shù)值模擬[J]. 核技術(shù), 2015, 38(1): 010602. DOI: 10.11889/j.0253-3219.2015.hjs.38.010602

    CAO Yin, WU Yanhua, LIN Chao, et al. Numerical simulation of the performance of vortex diodes with chamfered vortex chambers[J]. Nuclear Techniques, 2015, 38(1): 010602. DOI: 10.11889/j.0253-3219.2015.hjs. 38.010602

    4 Martinez P, Galpin J. CFD modeling of the EPR primary circuit[J]. Nuclear Engineering and Design, 2014, 278: 529?541. DOI: 10.1016/j.nucengdes.2014.08.013

    5 Kao M T, Wu C Y, Chieng C C, et al. CFD analysis of PWR core top and reactor vessel upper plenum internal subdomain models[J]. Nuclear Engineering and Design, 2011, 241(10): 4181?4193. DOI: 10.1016/j.nucengdes. 2011.08.007

    6 Baliga R, Watts T N, Kamath H. Application of an integrated head assembly concept at pressurized water reactor commercial nuclear plants[C]. 2014 22nd International Conference on Nuclear Engineering, American Society of Mechanical Engineers, Prague, Czech Republic, 2014. DOI: 10.1115/ ICONE22-30916

    7 余志偉, 何培峰, 李燕, 等. M310堆頂冷卻結(jié)構(gòu)流場(chǎng)和溫度場(chǎng)數(shù)值仿真[J]. CAD/CAM與制造業(yè)信息化, 2014, 7(7): 43?48. DOI: 10.3969/j.issn.1671-8186.2014.07.039

    YU Zhiwei, HE Peifeng, LI Yan, et al. Numerical simulation on the flow and temperature field in M310 integrated head package[J]. Digital Manufacturing Industry, 2014, 7(7): 43?48. DOI: 10.3969/j.issn. 1671-8186.2014.07.039

    8 于浩, 張明, 馮少東, 等. CAP1000一體化堆頂組件風(fēng)冷系統(tǒng)流場(chǎng)分析[J]. 核技術(shù), 2013, 36(4): 040624. DOI: 10.11889/j.0253-3219.2013.hjs.36.040624

    YU Hao, ZHANG Ming, FENG Shaodong, et al. CAP1000 integrated head package airflow system fluid field analysis[J]. Nuclear Techniques, 2013, 36(4): 040624. DOI: 10.11889/j.0253-3219.2013.hjs.36.040624

    9 任重, 黃興元, 柳和生, 等. 基于可壓縮氣輔的聚合物擠出成型非等溫黏彈數(shù)值分析[J]. 化工學(xué)報(bào), 2015, 66(4): 1615?1623. DOI: 10.11949/j.issn.0438-1157. 20141674

    REN Zhong, HUANG Xingyuan, LIU Hesheng, et al. Non-isothermal viscoelastic numerical analysis of compressible gas-assisted polymer extrusion molding[J]. Journal of Chemical Industry and Engineering, 2015, 66(4): 1615?1623. DOI: 10.11949/j.issn.0438-1157. 20141674

    10 張錫文, 李亨, 姚朝暉. 多孔介質(zhì)/純流體耦合區(qū)域內(nèi)可壓縮氣體的流動(dòng)[J]. 化工學(xué)報(bào), 2003, 54(9): 1209?1214

    ZHANG Xiwen, LI Heng, YAO Zhaohui. Compressible gas flow in porous media/fluid coupled areas[J]. Journal of Chemical Industry and Engineering, 2003, 54(9): 1209?1214

    Numerical simulation on the convective heat transfer of compressive gas in ACP100 integrated head package

    HE Peifeng1XU Bin1LUO Ying1YU Hao2MA Ziqi2SUN Shanwen2ZHOU Jinxiong2

    1(Key Laboratory of Nuclear Reactor System Design Technology, Nuclear Power Institute of China, Chengdu 610200, China)
    2(State Key Laboratory for Strength and Vibration of Mechanical Structures, School of Aerospace, Xi’an Jiaotong University, Xi’an 710049, China)

    Background: One of the functions of the integrated head package is to cool down control rod drive mechanism (CRDM), which is realized through convective flow and heat transfer of cooling gas. Purpose: Concerning the integrated head package of ACP100 nuclear reactor, we rigorously compared the thermal-fluid computation results. The emphasis of numerical simulation and comparison was laid on the effect of gas compressibility on heat transfer of cooling gas. Methods: The complete mesh model is built by HyperMesh and the thermal-fluid computation results are simulated with commercial ANSYS/CFX software and treating fluid media as both compressible and incompressible gases. Results: The simulation results show that the gas compressibility has great impact on the distribution of temperature, velocity and pressure fields. Conclusion: Ignoring the compressibility of gas would give an underestimation of maximum temperature on CRDM surfaces, lower maximum velocity and smaller pressure drop on CRDM gaps.

    Compressible gas, Integrated head package, Computational fluid dynamic (CFD), Convective heat transfer

    ZHOU Jinxiong, E-mail: jxzhouxx@mail.xjtu.edu.cn

    TL45

    10.11889/j.0253-3219.2016.hjs.39.100601

    中國(guó)核動(dòng)力研究設(shè)計(jì)院核反應(yīng)堆系統(tǒng)設(shè)計(jì)技術(shù)重點(diǎn)實(shí)驗(yàn)室課題(No.HT-A100-02-2015002)資助

    何培峰,男,1981年出生,2006年于浙江大學(xué)獲碩士學(xué)位,研究領(lǐng)域?yàn)榉磻?yīng)堆結(jié)構(gòu)設(shè)計(jì)

    周進(jìn)雄,E-mail: jxzhouxx@mail.xjtu.edu.cn

    Supported by Key Laboratory of Nuclear Reactor System Design Technology Project (No.HT-A100-02-2015002)First author: HE Peifeng, male, born in 1981, graduated from Zhejiang University with a master’s degree in 2006, focusing on reactor structural design

    2016-05-03,

    2016-07-28

    猜你喜歡
    控制棒壁面流場(chǎng)
    二維有限長(zhǎng)度柔性壁面上T-S波演化的數(shù)值研究
    大型空冷汽輪發(fā)電機(jī)轉(zhuǎn)子三維流場(chǎng)計(jì)算
    CARR壽期對(duì)控制棒價(jià)值的影響研究
    耐事故燃料控制棒備選材料的燃耗分析
    轉(zhuǎn)杯紡排雜區(qū)流場(chǎng)與排雜性能
    基于HYCOM的斯里蘭卡南部海域溫、鹽、流場(chǎng)統(tǒng)計(jì)分析
    壁面溫度對(duì)微型內(nèi)燃機(jī)燃燒特性的影響
    控制棒驅(qū)動(dòng)桿解鎖工具探討
    基于瞬態(tài)流場(chǎng)計(jì)算的滑動(dòng)軸承靜平衡位置求解
    改進(jìn)的源倍增方法測(cè)量控制棒價(jià)值
    少妇的逼水好多| 欧美日韩综合久久久久久| 夜夜骑夜夜射夜夜干| 久久久久精品久久久久真实原创| 寂寞人妻少妇视频99o| 国产在线视频一区二区| 高清在线视频一区二区三区| 99热全是精品| 午夜福利视频精品| 久久久精品94久久精品| 一级片'在线观看视频| 成人免费观看视频高清| 熟妇人妻不卡中文字幕| 另类亚洲欧美激情| 亚洲国产最新在线播放| 美女国产视频在线观看| 国产女主播在线喷水免费视频网站| 免费人成在线观看视频色| 国产欧美日韩一区二区三区在线 | 少妇熟女欧美另类| 九色成人免费人妻av| 超色免费av| 国产精品久久久久久久电影| 亚洲av二区三区四区| 国产成人精品久久久久久| 肉色欧美久久久久久久蜜桃| www.av在线官网国产| 99视频精品全部免费 在线| 天天操日日干夜夜撸| 汤姆久久久久久久影院中文字幕| 最黄视频免费看| 最近最新中文字幕免费大全7| 日韩精品有码人妻一区| 成人毛片60女人毛片免费| 亚洲人成网站在线观看播放| 国产一级毛片在线| 精品99又大又爽又粗少妇毛片| 国产亚洲一区二区精品| av卡一久久| 一边亲一边摸免费视频| 22中文网久久字幕| 一个人看视频在线观看www免费| 一级黄片播放器| 亚洲国产欧美在线一区| 超色免费av| 免费观看在线日韩| 日韩在线高清观看一区二区三区| 免费少妇av软件| 欧美成人精品欧美一级黄| 久久青草综合色| 九色成人免费人妻av| 超色免费av| 九九在线视频观看精品| 免费看不卡的av| 99久久人妻综合| 狂野欧美白嫩少妇大欣赏| 少妇猛男粗大的猛烈进出视频| 国产男女内射视频| 亚洲美女黄色视频免费看| 久久人妻熟女aⅴ| 人妻 亚洲 视频| 久热久热在线精品观看| 老司机亚洲免费影院| 黄色怎么调成土黄色| 91在线精品国自产拍蜜月| 国产极品粉嫩免费观看在线 | 久久99蜜桃精品久久| 久久精品人人爽人人爽视色| 80岁老熟妇乱子伦牲交| 亚洲精品一区蜜桃| 韩国av在线不卡| 免费观看的影片在线观看| 久久ye,这里只有精品| 菩萨蛮人人尽说江南好唐韦庄| 这个男人来自地球电影免费观看 | 一区二区日韩欧美中文字幕 | 欧美少妇被猛烈插入视频| 国产极品粉嫩免费观看在线 | 我的女老师完整版在线观看| 国产精品偷伦视频观看了| 你懂的网址亚洲精品在线观看| 欧美激情国产日韩精品一区| 免费黄频网站在线观看国产| 国产一级毛片在线| 少妇人妻精品综合一区二区| 美女cb高潮喷水在线观看| 久久影院123| 少妇丰满av| 91久久精品电影网| 免费av不卡在线播放| 欧美激情 高清一区二区三区| 久久久久久久亚洲中文字幕| 久久99热6这里只有精品| 18禁在线无遮挡免费观看视频| a级毛片免费高清观看在线播放| 国产精品一国产av| 三级国产精品欧美在线观看| 亚洲,一卡二卡三卡| 欧美少妇被猛烈插入视频| 国产欧美另类精品又又久久亚洲欧美| 日韩不卡一区二区三区视频在线| 亚洲精品自拍成人| 成年美女黄网站色视频大全免费 | 丰满乱子伦码专区| 免费不卡的大黄色大毛片视频在线观看| 91国产中文字幕| 亚洲精品乱码久久久久久按摩| 看十八女毛片水多多多| 午夜免费观看性视频| 国产日韩欧美亚洲二区| 久热久热在线精品观看| 国产成人aa在线观看| 国产精品熟女久久久久浪| 嫩草影院入口| 777米奇影视久久| 亚洲国产av影院在线观看| 91国产中文字幕| 午夜福利网站1000一区二区三区| 亚洲国产精品999| 日韩一区二区视频免费看| 丝袜在线中文字幕| 老司机亚洲免费影院| 午夜91福利影院| 9色porny在线观看| 亚洲国产欧美在线一区| 亚洲av.av天堂| 亚洲综合色惰| 亚洲色图 男人天堂 中文字幕 | 三级国产精品片| 亚洲国产精品成人久久小说| 一区二区三区免费毛片| 人体艺术视频欧美日本| 18+在线观看网站| 亚洲三级黄色毛片| 久久毛片免费看一区二区三区| 两个人免费观看高清视频| 男女无遮挡免费网站观看| 亚洲精品乱码久久久v下载方式| 日日啪夜夜爽| 丁香六月天网| tube8黄色片| 日韩视频在线欧美| av在线老鸭窝| 亚洲国产精品国产精品| 在线亚洲精品国产二区图片欧美 | 美女大奶头黄色视频| 国产极品天堂在线| 一级,二级,三级黄色视频| 精品熟女少妇av免费看| 精品国产一区二区三区久久久樱花| 色94色欧美一区二区| 久久精品人人爽人人爽视色| 色94色欧美一区二区| 国产成人精品福利久久| 看十八女毛片水多多多| 黑丝袜美女国产一区| 亚洲中文av在线| 国产一区二区三区综合在线观看 | 国产深夜福利视频在线观看| 97精品久久久久久久久久精品| 最近的中文字幕免费完整| 欧美3d第一页| 久久精品国产亚洲av天美| 乱人伦中国视频| 欧美精品高潮呻吟av久久| 国产精品一二三区在线看| 国产日韩欧美视频二区| 国产69精品久久久久777片| 另类精品久久| 国产精品蜜桃在线观看| 亚洲精品中文字幕在线视频| 中文字幕制服av| 亚洲av不卡在线观看| 自线自在国产av| 亚洲精品一区蜜桃| 亚洲欧美色中文字幕在线| 亚洲丝袜综合中文字幕| 在线播放无遮挡| 肉色欧美久久久久久久蜜桃| 天天躁夜夜躁狠狠久久av| 大香蕉久久成人网| 日韩制服骚丝袜av| 菩萨蛮人人尽说江南好唐韦庄| 久久久久视频综合| 午夜日本视频在线| 日韩av在线免费看完整版不卡| 国产成人一区二区在线| 国产探花极品一区二区| 国产成人一区二区在线| 亚洲欧洲国产日韩| 极品人妻少妇av视频| 日韩伦理黄色片| 搡女人真爽免费视频火全软件| 菩萨蛮人人尽说江南好唐韦庄| 免费不卡的大黄色大毛片视频在线观看| 午夜激情久久久久久久| 色哟哟·www| 成人18禁高潮啪啪吃奶动态图 | 亚洲欧美中文字幕日韩二区| 少妇人妻精品综合一区二区| 免费黄网站久久成人精品| 精品一区二区三卡| 国产精品麻豆人妻色哟哟久久| 午夜福利视频在线观看免费| 国产成人午夜福利电影在线观看| 欧美人与善性xxx| 国产精品一区二区在线观看99| 日本黄色日本黄色录像| 国产免费一级a男人的天堂| 亚洲美女视频黄频| 精品久久蜜臀av无| 男女高潮啪啪啪动态图| 国产高清有码在线观看视频| 能在线免费看毛片的网站| 嫩草影院入口| 大香蕉97超碰在线| 一边摸一边做爽爽视频免费| 亚洲人成网站在线播| 十八禁高潮呻吟视频| 久久毛片免费看一区二区三区| 一区二区日韩欧美中文字幕 | 少妇熟女欧美另类| 亚洲国产欧美日韩在线播放| 国产色爽女视频免费观看| 欧美变态另类bdsm刘玥| 久久久久久久久久成人| 国产av码专区亚洲av| 高清黄色对白视频在线免费看| 亚洲精品日韩av片在线观看| 欧美日韩在线观看h| 成人毛片a级毛片在线播放| 最新的欧美精品一区二区| 久久精品国产自在天天线| 三级国产精品欧美在线观看| 成年美女黄网站色视频大全免费 | 国产一区亚洲一区在线观看| 亚洲精品国产色婷婷电影| 久久午夜福利片| 老司机亚洲免费影院| 一级爰片在线观看| 十八禁网站网址无遮挡| 亚洲av二区三区四区| 男女啪啪激烈高潮av片| 少妇丰满av| 国产精品三级大全| 国产av国产精品国产| 婷婷色综合www| 亚洲五月色婷婷综合| 十八禁高潮呻吟视频| 另类亚洲欧美激情| 国产成人精品在线电影| 中文天堂在线官网| 国产一区二区三区av在线| 成人国产麻豆网| 中国国产av一级| 丰满少妇做爰视频| 国产日韩欧美视频二区| 午夜福利视频在线观看免费| 精品人妻一区二区三区麻豆| 欧美变态另类bdsm刘玥| 五月开心婷婷网| 成年美女黄网站色视频大全免费 | 春色校园在线视频观看| 九色亚洲精品在线播放| 熟妇人妻不卡中文字幕| 日日摸夜夜添夜夜添av毛片| 99国产精品免费福利视频| 久久鲁丝午夜福利片| 欧美三级亚洲精品| 午夜日本视频在线| 人人妻人人澡人人看| 午夜久久久在线观看| 亚洲,一卡二卡三卡| 成年人午夜在线观看视频| 超碰97精品在线观看| 亚洲成人一二三区av| 国产乱人偷精品视频| 精品亚洲乱码少妇综合久久| 免费少妇av软件| 十八禁网站网址无遮挡| 边亲边吃奶的免费视频| 午夜av观看不卡| 一级片'在线观看视频| 插阴视频在线观看视频| 亚洲欧美日韩另类电影网站| 老女人水多毛片| 黄色一级大片看看| 狠狠婷婷综合久久久久久88av| 国产午夜精品久久久久久一区二区三区| 啦啦啦在线观看免费高清www| 国产精品免费大片| av在线播放精品| 一边摸一边做爽爽视频免费| 亚洲三级黄色毛片| 免费高清在线观看日韩| 国产一区二区在线观看日韩| 久久久精品94久久精品| 一本久久精品| 狠狠精品人妻久久久久久综合| 国产精品不卡视频一区二区| 亚洲精品亚洲一区二区| 777米奇影视久久| 免费观看的影片在线观看| 免费黄网站久久成人精品| 亚洲性久久影院| 一级,二级,三级黄色视频| 爱豆传媒免费全集在线观看| 啦啦啦中文免费视频观看日本| 午夜激情福利司机影院| 一个人免费看片子| 国产成人午夜福利电影在线观看| 国产免费现黄频在线看| av在线观看视频网站免费| 中文精品一卡2卡3卡4更新| 国产成人免费无遮挡视频| 久久人人爽人人片av| tube8黄色片| 亚洲经典国产精华液单| 亚洲av中文av极速乱| 一本—道久久a久久精品蜜桃钙片| 国产黄片视频在线免费观看| 另类精品久久| 一个人看视频在线观看www免费| 亚洲人成77777在线视频| 如日韩欧美国产精品一区二区三区 | 精品一区二区三卡| 九九在线视频观看精品| 日本色播在线视频| a级毛片免费高清观看在线播放| 日本欧美国产在线视频| 精品亚洲乱码少妇综合久久| 欧美老熟妇乱子伦牲交| 高清黄色对白视频在线免费看| 五月伊人婷婷丁香| 国产男女内射视频| 免费大片18禁| 国产 精品1| 日本欧美视频一区| 久久精品人人爽人人爽视色| 制服诱惑二区| 欧美另类一区| 亚洲五月色婷婷综合| 黄色一级大片看看| 久久99精品国语久久久| 少妇的逼好多水| 黄色毛片三级朝国网站| 三级国产精品片| 国产又色又爽无遮挡免| 三级国产精品欧美在线观看| 色吧在线观看| 王馨瑶露胸无遮挡在线观看| 蜜桃国产av成人99| 亚洲性久久影院| 91国产中文字幕| 中国国产av一级| 国产日韩欧美视频二区| 国产午夜精品一二区理论片| 两个人的视频大全免费| 在线观看三级黄色| 插阴视频在线观看视频| 99久久人妻综合| 精品少妇黑人巨大在线播放| 久久热精品热| 亚洲国产欧美在线一区| 热re99久久国产66热| 免费看av在线观看网站| 欧美变态另类bdsm刘玥| a级毛片在线看网站| 大片电影免费在线观看免费| 一本久久精品| 涩涩av久久男人的天堂| 成人亚洲欧美一区二区av| 少妇人妻久久综合中文| 日韩免费高清中文字幕av| av在线app专区| 久久久国产精品麻豆| 3wmmmm亚洲av在线观看| 天堂8中文在线网| 久久99热6这里只有精品| 国产精品人妻久久久久久| 两个人的视频大全免费| 2018国产大陆天天弄谢| 欧美日韩av久久| 精品久久久噜噜| 亚洲av欧美aⅴ国产| 日韩强制内射视频| 欧美 日韩 精品 国产| 全区人妻精品视频| 国产精品免费大片| 国产日韩欧美视频二区| 免费看光身美女| 99久久中文字幕三级久久日本| 热99国产精品久久久久久7| av视频免费观看在线观看| 免费观看a级毛片全部| 日韩欧美精品免费久久| 赤兔流量卡办理| 最黄视频免费看| 女性被躁到高潮视频| 十八禁网站网址无遮挡| 十分钟在线观看高清视频www| 99热全是精品| 美女脱内裤让男人舔精品视频| 大片电影免费在线观看免费| a级毛色黄片| 免费高清在线观看日韩| 国产精品国产三级国产av玫瑰| 国产精品99久久99久久久不卡 | 18+在线观看网站| 亚洲天堂av无毛| av不卡在线播放| a级片在线免费高清观看视频| 人人妻人人澡人人看| 色婷婷久久久亚洲欧美| 久久久久久久久久成人| 亚洲av.av天堂| 少妇人妻精品综合一区二区| 日韩熟女老妇一区二区性免费视频| 各种免费的搞黄视频| 在线观看www视频免费| 国产高清不卡午夜福利| 黑人猛操日本美女一级片| www.av在线官网国产| 日韩制服骚丝袜av| 大香蕉97超碰在线| 人妻一区二区av| 久久久久精品久久久久真实原创| 五月天丁香电影| 久久久久国产网址| 少妇的逼好多水| 中文字幕亚洲精品专区| a级毛片在线看网站| 欧美日韩国产mv在线观看视频| 国产不卡av网站在线观看| 我的女老师完整版在线观看| 久久人人爽人人片av| 国产成人精品在线电影| 国产精品一区二区在线不卡| 老女人水多毛片| 一区在线观看完整版| 免费黄色在线免费观看| 夫妻性生交免费视频一级片| 欧美性感艳星| 久久精品熟女亚洲av麻豆精品| 男人操女人黄网站| 99久久精品国产国产毛片| 国产精品不卡视频一区二区| www.色视频.com| 青春草国产在线视频| 国产熟女欧美一区二区| 黑人猛操日本美女一级片| 女的被弄到高潮叫床怎么办| 亚洲av日韩在线播放| 一区二区三区精品91| 成人漫画全彩无遮挡| 午夜福利视频在线观看免费| 日本色播在线视频| 亚洲国产精品成人久久小说| 日本黄色日本黄色录像| 精品久久国产蜜桃| 国产爽快片一区二区三区| 制服丝袜香蕉在线| 免费大片18禁| 免费人成在线观看视频色| 亚洲欧洲国产日韩| 在线免费观看不下载黄p国产| 国产伦理片在线播放av一区| 一本—道久久a久久精品蜜桃钙片| 日韩一区二区三区影片| av网站免费在线观看视频| 18禁在线无遮挡免费观看视频| 国产成人午夜福利电影在线观看| 亚洲国产av影院在线观看| 少妇熟女欧美另类| 精品国产一区二区三区久久久樱花| 精品一区在线观看国产| 久久久a久久爽久久v久久| 少妇被粗大猛烈的视频| 亚洲欧美一区二区三区国产| 欧美xxⅹ黑人| 男男h啪啪无遮挡| 尾随美女入室| 国产男人的电影天堂91| 18禁动态无遮挡网站| 天天躁夜夜躁狠狠久久av| 欧美日韩国产mv在线观看视频| 日日摸夜夜添夜夜添av毛片| 777米奇影视久久| 18禁观看日本| 日韩一本色道免费dvd| 国模一区二区三区四区视频| 热99久久久久精品小说推荐| 大片电影免费在线观看免费| 久久精品久久久久久久性| 制服丝袜香蕉在线| 丝袜美足系列| 亚洲国产日韩一区二区| 丰满乱子伦码专区| 精品久久久噜噜| 色婷婷av一区二区三区视频| 国产69精品久久久久777片| 男的添女的下面高潮视频| 一本大道久久a久久精品| 日本黄色日本黄色录像| 一级毛片电影观看| 国产免费又黄又爽又色| 亚洲伊人久久精品综合| 国产精品蜜桃在线观看| 国产黄频视频在线观看| 天堂俺去俺来也www色官网| 天天操日日干夜夜撸| 亚洲综合精品二区| 国产在线免费精品| 制服人妻中文乱码| 男人添女人高潮全过程视频| 日本vs欧美在线观看视频| 天堂俺去俺来也www色官网| 满18在线观看网站| 午夜老司机福利剧场| 国产精品一区二区在线观看99| 天美传媒精品一区二区| 热99久久久久精品小说推荐| 777米奇影视久久| 一个人免费看片子| 久久久久久久久大av| 国产爽快片一区二区三区| 欧美日韩视频精品一区| 水蜜桃什么品种好| 国产成人精品无人区| 国模一区二区三区四区视频| 女性生殖器流出的白浆| 免费观看性生交大片5| 久久av网站| 激情五月婷婷亚洲| 久久婷婷青草| 国产国拍精品亚洲av在线观看| 亚洲少妇的诱惑av| 日本wwww免费看| 亚洲精品国产av蜜桃| 超色免费av| 久久人人爽人人片av| 久久久久久久久久久久大奶| 我要看黄色一级片免费的| 欧美日韩综合久久久久久| 美女福利国产在线| 欧美激情极品国产一区二区三区 | 日本vs欧美在线观看视频| 亚洲美女搞黄在线观看| 精品国产乱码久久久久久小说| 蜜桃国产av成人99| 黄片无遮挡物在线观看| 欧美精品亚洲一区二区| 亚洲第一区二区三区不卡| 在现免费观看毛片| 亚洲欧美精品自产自拍| 国产在视频线精品| 99国产综合亚洲精品| 亚洲精品国产av成人精品| 成年人免费黄色播放视频| 日日摸夜夜添夜夜爱| 日韩,欧美,国产一区二区三区| 久久人人爽av亚洲精品天堂| 熟女av电影| 亚洲情色 制服丝袜| 十八禁网站网址无遮挡| 免费大片18禁| 亚洲国产精品999| 99九九线精品视频在线观看视频| 国产午夜精品一二区理论片| 国产亚洲午夜精品一区二区久久| 狠狠精品人妻久久久久久综合| 美女主播在线视频| 美女视频免费永久观看网站| 日本午夜av视频| 国产亚洲精品久久久com| 美女脱内裤让男人舔精品视频| av卡一久久| 亚洲人成77777在线视频| 最后的刺客免费高清国语| 高清午夜精品一区二区三区| 亚洲精品国产av成人精品| 免费观看在线日韩| 色婷婷久久久亚洲欧美| 边亲边吃奶的免费视频| 丝袜脚勾引网站| 亚洲精品视频女| 18禁动态无遮挡网站| 26uuu在线亚洲综合色| 丝瓜视频免费看黄片| 国产成人免费观看mmmm| av在线app专区| 大片电影免费在线观看免费| 国产又色又爽无遮挡免| 亚洲精品久久午夜乱码| 久久精品人人爽人人爽视色| 大片免费播放器 马上看| 肉色欧美久久久久久久蜜桃| 亚洲欧美一区二区三区国产| 男人操女人黄网站| 久久精品国产亚洲av天美| 日日撸夜夜添| 精品人妻在线不人妻| 日韩av免费高清视频| 欧美变态另类bdsm刘玥| 色视频在线一区二区三区| 亚洲av二区三区四区| av国产久精品久网站免费入址| 欧美日韩一区二区视频在线观看视频在线| 婷婷色综合www| 亚洲国产精品一区二区三区在线| 国产片特级美女逼逼视频| 男女边摸边吃奶| 晚上一个人看的免费电影| 十分钟在线观看高清视频www| 男女无遮挡免费网站观看|