• <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à)值
    久久久国产一区二区| 午夜福利免费观看在线| 亚洲少妇的诱惑av| 亚洲精品国产区一区二| 国产精品久久久久久精品电影小说| 中文字幕高清在线视频| 成人漫画全彩无遮挡| 久久久久国产一级毛片高清牌| 精品久久久久久电影网| 欧美人与性动交α欧美精品济南到| av网站免费在线观看视频| 亚洲国产精品一区二区三区在线| 国产高清不卡午夜福利| 深夜精品福利| 熟女av电影| 国产av一区二区精品久久| 青春草视频在线免费观看| 国产成人系列免费观看| 巨乳人妻的诱惑在线观看| 十八禁高潮呻吟视频| 99久久人妻综合| 中文字幕亚洲精品专区| 国产深夜福利视频在线观看| 亚洲国产成人一精品久久久| 天堂俺去俺来也www色官网| 99久国产av精品国产电影| 天天躁日日躁夜夜躁夜夜| 高清不卡的av网站| 久久97久久精品| 啦啦啦视频在线资源免费观看| 日本av手机在线免费观看| 欧美日韩视频高清一区二区三区二| 午夜激情久久久久久久| 亚洲综合色网址| 99香蕉大伊视频| 中文欧美无线码| 久久99精品国语久久久| 男女边摸边吃奶| 这个男人来自地球电影免费观看 | 久久久亚洲精品成人影院| 交换朋友夫妻互换小说| 亚洲精品久久久久久婷婷小说| 国产极品天堂在线| 成人黄色视频免费在线看| 国产熟女午夜一区二区三区| 国产97色在线日韩免费| 一本色道久久久久久精品综合| 免费久久久久久久精品成人欧美视频| 一本色道久久久久久精品综合| 一本色道久久久久久精品综合| 爱豆传媒免费全集在线观看| 色吧在线观看| 777米奇影视久久| 国产精品三级大全| 在线观看三级黄色| 欧美少妇被猛烈插入视频| 黄片播放在线免费| 久久国产亚洲av麻豆专区| 夜夜骑夜夜射夜夜干| 最近2019中文字幕mv第一页| 在线 av 中文字幕| 在现免费观看毛片| 国产片特级美女逼逼视频| 狂野欧美激情性bbbbbb| 国产精品久久久久久精品电影小说| 97人妻天天添夜夜摸| 免费黄色在线免费观看| 亚洲成人一二三区av| 伊人久久大香线蕉亚洲五| 精品人妻一区二区三区麻豆| 欧美黑人精品巨大| 欧美老熟妇乱子伦牲交| 精品人妻一区二区三区麻豆| 一个人免费看片子| h视频一区二区三区| 午夜福利,免费看| 午夜福利在线免费观看网站| 亚洲人成77777在线视频| 欧美激情高清一区二区三区 | 国产精品一区二区在线不卡| 久久午夜综合久久蜜桃| 久久这里只有精品19| 日韩欧美一区视频在线观看| 熟女av电影| 中文字幕人妻丝袜一区二区 | 久久精品亚洲av国产电影网| 宅男免费午夜| 深夜精品福利| 成人免费观看视频高清| 亚洲国产精品一区二区三区在线| 这个男人来自地球电影免费观看 | 久久人人爽人人片av| 秋霞在线观看毛片| 久久人人97超碰香蕉20202| 久久婷婷青草| 久久精品国产亚洲av涩爱| 国产高清不卡午夜福利| 国产女主播在线喷水免费视频网站| 在线观看免费日韩欧美大片| 少妇的丰满在线观看| 久久午夜综合久久蜜桃| 国产成人a∨麻豆精品| 国产成人a∨麻豆精品| 最近的中文字幕免费完整| 国产成人欧美| videos熟女内射| 在线观看免费视频网站a站| 99re6热这里在线精品视频| 大陆偷拍与自拍| 国产精品 国内视频| 日本av免费视频播放| av片东京热男人的天堂| 成人三级做爰电影| 成人漫画全彩无遮挡| 欧美人与性动交α欧美软件| 精品少妇一区二区三区视频日本电影 | 最近最新中文字幕大全免费视频| 中国美女看黄片| 成年人黄色毛片网站| 非洲黑人性xxxx精品又粗又长| 女人精品久久久久毛片| 亚洲精品国产色婷婷电影| 一级毛片高清免费大全| 成人国产综合亚洲| 亚洲三区欧美一区| 神马国产精品三级电影在线观看 | 美女扒开内裤让男人捅视频| 人人澡人人妻人| 亚洲av片天天在线观看| 在线观看66精品国产| 韩国精品一区二区三区| 啦啦啦免费观看视频1| 日本精品一区二区三区蜜桃| 黑丝袜美女国产一区| 又大又爽又粗| 操美女的视频在线观看| 久久伊人香网站| 看片在线看免费视频| 9热在线视频观看99| 午夜视频精品福利| 美女午夜性视频免费| 免费高清视频大片| 18禁美女被吸乳视频| 如日韩欧美国产精品一区二区三区| 欧美黄色淫秽网站| 亚洲欧洲精品一区二区精品久久久| 黑人操中国人逼视频| 一区二区三区激情视频| 亚洲国产毛片av蜜桃av| 天堂动漫精品| 亚洲精品中文字幕一二三四区| 18禁黄网站禁片午夜丰满| 亚洲国产精品久久男人天堂| 国产成人精品久久二区二区免费| 国产精品1区2区在线观看.| 男人舔女人的私密视频| 欧美黑人精品巨大| 欧美成人性av电影在线观看| 久久久久久久久久久久大奶| 好男人在线观看高清免费视频 | 欧美色欧美亚洲另类二区 | 国产免费av片在线观看野外av| 亚洲精品粉嫩美女一区| a级毛片在线看网站| 后天国语完整版免费观看| 中文字幕色久视频| 国产成人精品久久二区二区免费| 一边摸一边做爽爽视频免费| 国产野战对白在线观看| 91麻豆精品激情在线观看国产| 两性午夜刺激爽爽歪歪视频在线观看 | 国产亚洲精品久久久久久毛片| 韩国精品一区二区三区| 免费av毛片视频| 香蕉丝袜av| 别揉我奶头~嗯~啊~动态视频| 国产精品永久免费网站| 精品久久久精品久久久| 每晚都被弄得嗷嗷叫到高潮| 一区福利在线观看| 亚洲人成77777在线视频| 久久久国产成人免费| 午夜福利成人在线免费观看| 国产在线精品亚洲第一网站| 国产精品亚洲美女久久久| 大型黄色视频在线免费观看| 宅男免费午夜| 亚洲七黄色美女视频| 一进一出抽搐动态| 国产日韩一区二区三区精品不卡| 亚洲最大成人中文| ponron亚洲| 他把我摸到了高潮在线观看| 日日干狠狠操夜夜爽| 亚洲人成伊人成综合网2020| 亚洲 欧美一区二区三区| 动漫黄色视频在线观看| 成人三级黄色视频| 亚洲无线在线观看| 国产精华一区二区三区| 精品久久久久久成人av| av电影中文网址| 国产高清视频在线播放一区| 国产xxxxx性猛交| 久久香蕉国产精品| 色精品久久人妻99蜜桃| 日韩一卡2卡3卡4卡2021年| 人妻丰满熟妇av一区二区三区| 国语自产精品视频在线第100页| 国产男靠女视频免费网站| 成人免费观看视频高清| 高清毛片免费观看视频网站| 久久精品国产99精品国产亚洲性色 | 欧美成人免费av一区二区三区| 丁香六月欧美| 精品一区二区三区视频在线观看免费| av欧美777| 91av网站免费观看| 中文字幕人妻熟女乱码| 亚洲欧洲精品一区二区精品久久久| 好看av亚洲va欧美ⅴa在| 国产午夜福利久久久久久| 男人舔女人下体高潮全视频| 午夜福利免费观看在线| 18禁裸乳无遮挡免费网站照片 | 美女高潮到喷水免费观看| 欧美在线一区亚洲| 久久久国产成人精品二区| 日日摸夜夜添夜夜添小说| 亚洲国产高清在线一区二区三 | 亚洲中文日韩欧美视频| 国产精品乱码一区二三区的特点 | 久久性视频一级片| 自拍欧美九色日韩亚洲蝌蚪91| 淫妇啪啪啪对白视频| 亚洲成人免费电影在线观看| 日本欧美视频一区| av有码第一页| tocl精华| 国产成人av激情在线播放| 在线观看66精品国产| 搡老岳熟女国产| 大陆偷拍与自拍| 99久久国产精品久久久| 精品第一国产精品| 变态另类成人亚洲欧美熟女 | 99久久精品国产亚洲精品| 国产精品野战在线观看| 亚洲欧美激情在线| 欧美黑人精品巨大| 欧美在线一区亚洲| 精品久久蜜臀av无| 黑人操中国人逼视频| 伊人久久大香线蕉亚洲五| 长腿黑丝高跟| 亚洲成av人片免费观看| av天堂久久9| 国语自产精品视频在线第100页| 国产精品久久久av美女十八| 亚洲狠狠婷婷综合久久图片| 日本一区二区免费在线视频| 精品第一国产精品| 国产亚洲av高清不卡| 在线观看一区二区三区| 99国产综合亚洲精品| 亚洲人成电影免费在线| 久久国产精品人妻蜜桃| 久久国产精品影院| 午夜福利,免费看| 国产av精品麻豆| av福利片在线| 人妻丰满熟妇av一区二区三区| 亚洲精品国产区一区二| 国产一区二区三区综合在线观看| 久久性视频一级片| 正在播放国产对白刺激| 亚洲专区国产一区二区| 日韩欧美三级三区| 午夜免费激情av| 国产日韩一区二区三区精品不卡| 亚洲欧美激情在线| 黄片小视频在线播放| www.www免费av| 正在播放国产对白刺激| 九色亚洲精品在线播放| 欧美日韩亚洲国产一区二区在线观看| 久久午夜综合久久蜜桃| 午夜两性在线视频| 在线永久观看黄色视频| 国产精品 国内视频| 伦理电影免费视频| 欧美日韩黄片免| 国产欧美日韩精品亚洲av| 麻豆国产av国片精品| av电影中文网址| 欧美av亚洲av综合av国产av| 两个人视频免费观看高清| 成人18禁在线播放| 真人做人爱边吃奶动态| 亚洲国产高清在线一区二区三 | 在线观看www视频免费| 欧美黑人精品巨大| 国产成人精品久久二区二区91| 日韩欧美在线二视频| 成年版毛片免费区| 两个人免费观看高清视频| a在线观看视频网站| 日韩大尺度精品在线看网址 | 日韩有码中文字幕| 成人精品一区二区免费| 亚洲伊人色综图| 亚洲九九香蕉| 久久草成人影院| 岛国在线观看网站| 男女下面进入的视频免费午夜 | 亚洲无线在线观看| www.熟女人妻精品国产| 好男人电影高清在线观看| 国产成人精品久久二区二区91| 免费观看精品视频网站| 国产成人系列免费观看| 香蕉国产在线看| 九色国产91popny在线| 亚洲三区欧美一区| 在线免费观看的www视频| 午夜精品久久久久久毛片777| 色播在线永久视频| 日韩成人在线观看一区二区三区| 日日干狠狠操夜夜爽| www.熟女人妻精品国产| 国产99白浆流出| 免费高清视频大片| 精品无人区乱码1区二区| 久久精品亚洲精品国产色婷小说| av欧美777| 成人手机av| 精品国内亚洲2022精品成人| 亚洲全国av大片| 精品国产乱码久久久久久男人| 亚洲男人的天堂狠狠| 在线观看www视频免费| 国产三级黄色录像| 日本撒尿小便嘘嘘汇集6| 美女国产高潮福利片在线看| 久久婷婷成人综合色麻豆| av天堂在线播放| 日韩 欧美 亚洲 中文字幕| 国产视频一区二区在线看| 中文字幕av电影在线播放| 欧美日韩一级在线毛片| 欧美激情高清一区二区三区| 精品欧美国产一区二区三| 多毛熟女@视频| 亚洲中文字幕一区二区三区有码在线看 | 国产欧美日韩综合在线一区二区| 国产成人av激情在线播放| 波多野结衣高清无吗| 欧美在线一区亚洲| 亚洲在线自拍视频| 国产男靠女视频免费网站| 亚洲国产精品成人综合色| 人人澡人人妻人| 亚洲精品一区av在线观看| 久久久久久久精品吃奶| 九色国产91popny在线| 亚洲专区中文字幕在线| 日韩视频一区二区在线观看| 亚洲av电影不卡..在线观看| 欧美亚洲日本最大视频资源| 欧美日韩中文字幕国产精品一区二区三区 | 精品熟女少妇八av免费久了| 国产精品亚洲美女久久久| 婷婷六月久久综合丁香| 久久九九热精品免费| 男男h啪啪无遮挡| 国产又爽黄色视频| 在线观看午夜福利视频| 50天的宝宝边吃奶边哭怎么回事| 乱人伦中国视频| 免费搜索国产男女视频| 免费观看精品视频网站| 亚洲精品一区av在线观看| 久久久久久久久久久久大奶| 亚洲欧美日韩高清在线视频| 午夜影院日韩av| 一本久久中文字幕| 黑人操中国人逼视频| 在线观看免费视频网站a站| 精品熟女少妇八av免费久了| 欧美日韩瑟瑟在线播放| 美女高潮喷水抽搐中文字幕| 性欧美人与动物交配| 一边摸一边做爽爽视频免费| 露出奶头的视频| 99久久综合精品五月天人人| 少妇粗大呻吟视频| 久久天躁狠狠躁夜夜2o2o| 麻豆一二三区av精品| 日韩精品免费视频一区二区三区| 在线观看免费日韩欧美大片| 女人爽到高潮嗷嗷叫在线视频| 国产蜜桃级精品一区二区三区| 可以在线观看的亚洲视频| 欧美激情久久久久久爽电影 | 国产精品 欧美亚洲| 久久精品91蜜桃| 叶爱在线成人免费视频播放| 校园春色视频在线观看| 亚洲一区二区三区色噜噜| 叶爱在线成人免费视频播放| 欧美大码av| 欧美av亚洲av综合av国产av| 国产成人影院久久av| 我的亚洲天堂| 啦啦啦免费观看视频1| 淫妇啪啪啪对白视频| 国产真人三级小视频在线观看| 精品欧美国产一区二区三| 757午夜福利合集在线观看| АⅤ资源中文在线天堂| 亚洲国产精品成人综合色| 成人欧美大片| av欧美777| x7x7x7水蜜桃| 精品少妇一区二区三区视频日本电影| 国产成年人精品一区二区| 久久久久国产一级毛片高清牌| 亚洲电影在线观看av| 免费观看人在逋| 人妻丰满熟妇av一区二区三区| 男女做爰动态图高潮gif福利片 | 九色国产91popny在线| 亚洲中文字幕一区二区三区有码在线看 | av中文乱码字幕在线| 亚洲精华国产精华精| 午夜免费鲁丝| 一a级毛片在线观看| 黄色片一级片一级黄色片| 国语自产精品视频在线第100页| 中文字幕精品免费在线观看视频| 久热爱精品视频在线9| 91成年电影在线观看| 两个人看的免费小视频| 午夜影院日韩av| 免费看十八禁软件| 午夜福利影视在线免费观看| 男女床上黄色一级片免费看| 夜夜夜夜夜久久久久| 亚洲电影在线观看av| 国产成人精品久久二区二区91| 国产1区2区3区精品| 国产私拍福利视频在线观看| 女警被强在线播放| 国产精品免费视频内射| 极品教师在线免费播放| 人人澡人人妻人| 国产又爽黄色视频| 免费在线观看黄色视频的| 又紧又爽又黄一区二区| 久久中文字幕一级| 啦啦啦韩国在线观看视频| 国产av一区二区精品久久| aaaaa片日本免费| 亚洲五月色婷婷综合| 免费女性裸体啪啪无遮挡网站| 99精品在免费线老司机午夜| 一区二区日韩欧美中文字幕| 国产麻豆成人av免费视频| 日本黄色视频三级网站网址| 国产成人精品无人区| 亚洲av片天天在线观看| 色综合亚洲欧美另类图片| 免费久久久久久久精品成人欧美视频| 此物有八面人人有两片| 亚洲五月婷婷丁香| 色老头精品视频在线观看| 制服丝袜大香蕉在线| 一级毛片高清免费大全| 国产91精品成人一区二区三区| 精品久久久久久成人av| 91成人精品电影| 可以在线观看的亚洲视频| 日本黄色视频三级网站网址| 亚洲国产欧美一区二区综合| 无限看片的www在线观看| 国产午夜精品久久久久久| 精品第一国产精品| 一区二区日韩欧美中文字幕| 久久久国产成人精品二区| 韩国精品一区二区三区| 欧美成狂野欧美在线观看| 又紧又爽又黄一区二区| 亚洲七黄色美女视频| 亚洲午夜理论影院| 久久精品影院6| 亚洲天堂国产精品一区在线| 国产精品香港三级国产av潘金莲| 97超级碰碰碰精品色视频在线观看| 色哟哟哟哟哟哟| 久久久久九九精品影院| 亚洲av五月六月丁香网| 国内毛片毛片毛片毛片毛片| 波多野结衣巨乳人妻| 十八禁网站免费在线| 国产人伦9x9x在线观看| avwww免费| 中文字幕色久视频| 男女午夜视频在线观看| 国产欧美日韩综合在线一区二区| 亚洲精品在线美女| 高清毛片免费观看视频网站| 精品福利观看| 国产成人精品无人区| 岛国在线观看网站| 免费女性裸体啪啪无遮挡网站| 亚洲av成人一区二区三| 久久精品影院6| 91老司机精品| 黄片小视频在线播放| 久久午夜综合久久蜜桃| 国产蜜桃级精品一区二区三区| 亚洲精品国产一区二区精华液| 巨乳人妻的诱惑在线观看| 久久久国产成人精品二区| 涩涩av久久男人的天堂| 97超级碰碰碰精品色视频在线观看| 黑人巨大精品欧美一区二区蜜桃| 18禁黄网站禁片午夜丰满| 每晚都被弄得嗷嗷叫到高潮| 制服诱惑二区| 999精品在线视频| 国产色视频综合| 999精品在线视频| 久久午夜亚洲精品久久| 亚洲国产高清在线一区二区三 | 亚洲激情在线av| 色尼玛亚洲综合影院| 婷婷六月久久综合丁香| 亚洲天堂国产精品一区在线| 一区二区日韩欧美中文字幕| 老司机午夜福利在线观看视频| 亚洲成人精品中文字幕电影| 女人精品久久久久毛片| 免费看a级黄色片| 亚洲免费av在线视频| 国产亚洲精品第一综合不卡| 夜夜夜夜夜久久久久| av中文乱码字幕在线| 免费在线观看亚洲国产| 伊人久久大香线蕉亚洲五| 欧美日韩福利视频一区二区| 免费在线观看影片大全网站| 真人一进一出gif抽搐免费| 免费在线观看影片大全网站| 国产熟女午夜一区二区三区| 女人被躁到高潮嗷嗷叫费观| 亚洲一卡2卡3卡4卡5卡精品中文| 99久久久亚洲精品蜜臀av| 波多野结衣一区麻豆| 99久久久亚洲精品蜜臀av| 亚洲人成电影免费在线| 18禁观看日本| 极品人妻少妇av视频| 日本五十路高清| 变态另类丝袜制服| 最近最新中文字幕大全电影3 | 精品国产超薄肉色丝袜足j| 两性夫妻黄色片| 色哟哟哟哟哟哟| 69精品国产乱码久久久| 亚洲国产欧美日韩在线播放| 国产精品美女特级片免费视频播放器 | 精品卡一卡二卡四卡免费| 国产精品自产拍在线观看55亚洲| 色综合婷婷激情| 国产精品精品国产色婷婷| 无限看片的www在线观看| 日日干狠狠操夜夜爽| 欧美不卡视频在线免费观看 | 亚洲成人久久性| 99久久国产精品久久久| 国产三级在线视频| 久久久久国产一级毛片高清牌| 午夜a级毛片| 91大片在线观看| 十八禁网站免费在线| 淫妇啪啪啪对白视频| 欧美激情极品国产一区二区三区| 男男h啪啪无遮挡| 91av网站免费观看| 一边摸一边抽搐一进一出视频| 在线观看午夜福利视频| 色哟哟哟哟哟哟| 女警被强在线播放| 亚洲成人精品中文字幕电影| 亚洲av电影不卡..在线观看| 午夜福利一区二区在线看| 国产aⅴ精品一区二区三区波| 99精品欧美一区二区三区四区| 亚洲色图综合在线观看| 一个人免费在线观看的高清视频| 女人爽到高潮嗷嗷叫在线视频| 亚洲七黄色美女视频| 97人妻天天添夜夜摸| 成年人黄色毛片网站| 18禁观看日本| 精品欧美一区二区三区在线| 国产主播在线观看一区二区| 最近最新中文字幕大全电影3 | 精品电影一区二区在线| 久久久久久亚洲精品国产蜜桃av| 日韩精品免费视频一区二区三区| 最新美女视频免费是黄的| 久久精品影院6| 亚洲国产日韩欧美精品在线观看 |