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

    熱浮力對(duì)球床反應(yīng)堆內(nèi)燃料球表面對(duì)流傳熱的影響

    2021-08-08 06:12:10姚強(qiáng)郭雪巖楊帆
    能源研究與信息 2021年4期

    姚強(qiáng) 郭雪巖 楊帆

    摘要:采用大渦模擬(LES )方法研究了低雷諾數(shù)(Re =1015)下熱浮力對(duì)球床面心立方(FCC )單元內(nèi)的局部流動(dòng)和傳熱的影響。為了精確求解燃料球接觸面附近的流場(chǎng),球間接觸點(diǎn)表達(dá)為面接觸,采用結(jié)構(gòu)化網(wǎng)格處理。研究結(jié)果表明:熱浮力的存在會(huì)抑制球床各層平均速度的波動(dòng);靠近燃料球表面處的相對(duì)時(shí)均速度和相對(duì)時(shí)均溫度受熱浮力影響變化較大;在中心流場(chǎng)區(qū)域,熱浮力的存在會(huì)減小流場(chǎng)中速度分布的不對(duì)稱性,使速度最大降低約10%,時(shí)均溫度至少升高約20%;球表面上尤其在頂部及接觸面附近受熱浮力影響明顯,努塞爾數(shù) Nu最大降低約6%。

    關(guān)鍵詞:熱浮力;球床反應(yīng)堆;大渦模擬;面心立方;面接觸

    中圖分類號(hào): TK11+2??? 文獻(xiàn)標(biāo)志碼: A

    Effect of thermal buoyancy on the convective heat transfer on the surface of fuel pebbles in a pebble bed reactor

    YAO Qiang,GUO Xueyan,YANG Fan

    (School of Energy and Power Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China)

    Abstract:In this paper, large eddy simulation (LES) method was adopted to study the influence of thermal buoyancy at low Reynolds number (Re =1015) on local flow and heat transfer in a face- centered cubic (FCC) of a pebble bed. The plane-plane contact among the pebbles and structured grid were adopted to accurately solve the flow field near the interface of fuel pebbles. It is found that the existence of thermal buoyancy could restrain the fluctuation of averaged velocity at each layer. The influence of thermal buoyancy on relative time-averaged velocity and relative time-

    averaged temperature near the surface of fuel pebbles was great. In the center of flow field, the existence of thermal buoyancy could reduce the asymmetry of velocity distribution. The velocity decreased by 10% at most and the time-averaged temperature increased by 20% at least. On the spherical ?surface,? especially? at? the? top? and? near? the? interface,? obvious? influence? of thermalbuoyancy was observed. And the Nusselt number Nu number decreased by 6% at most.

    Keywords:thermal buoyancy; pebble bed reactor; large eddy simulation; face-centered cubic; plane-plane contact

    球床式高溫氣冷堆是第四代核反應(yīng)堆類型之一。在球床反應(yīng)堆中,氦氣作為冷卻劑,核燃料嵌在石墨慢化劑中并可以維持很高的溫度。這種結(jié)構(gòu)可以保證燃料不會(huì)熔化,使反應(yīng)堆具有良好的固有安全性[1–4]。該堆型具有結(jié)構(gòu)簡(jiǎn)單、燃料元件適合批量生產(chǎn)以及燃料球裝卸方便等特點(diǎn)。

    目前,很多學(xué)者對(duì)球床中的流動(dòng)和換熱特點(diǎn)進(jìn)行了研究。孟現(xiàn)柯等[5]用由碳鋼球堆積成球床,以蒸餾水為工質(zhì),采用電磁感應(yīng)加熱方式研究了球床通道內(nèi)部的換熱特性; Jia 等[6]用黑色玻璃球堆積成球床,以水為工質(zhì)研究了不同床型對(duì)內(nèi)部工質(zhì)流動(dòng)情況的影響。考慮到實(shí)驗(yàn)成本和復(fù)雜性,多數(shù)實(shí)驗(yàn)結(jié)果很難為球床反應(yīng)堆提供可靠的參考依據(jù)。Lee 等[7]于2007年采用k 一"湍流模型比較了不同球間間隙對(duì)流場(chǎng)分布的影響;2014年,F(xiàn)erng等[8]用雷諾應(yīng)力湍流模型、非結(jié)構(gòu)網(wǎng)格對(duì)體心立方( BCC )人工間隙和點(diǎn)接觸兩種接觸方式的共軛模型的模擬結(jié)果進(jìn)行了分析比較;2012—2015年, Shams 等[9–11]采用 q?DNS 方法、多面體網(wǎng)格對(duì)面心立方(FCC )人工間隙模型進(jìn)行了分析,發(fā)現(xiàn)在兩球之間的狹窄區(qū)域存在射流現(xiàn)象及速度場(chǎng)、溫度場(chǎng)分布的不對(duì)稱性;2017年,蔣旭等[12]采用大渦模擬( LES )方法、結(jié)構(gòu)網(wǎng)格研究了 BCC堆積模型的流場(chǎng)和溫度場(chǎng),通過(guò)與k一"湍流模型計(jì)算結(jié)果對(duì)比發(fā)現(xiàn),采用 LES 方法能夠更好地捕捉流場(chǎng)中的渦結(jié)構(gòu),瞬時(shí)溫度變化更具參考價(jià)值。采用數(shù)值方法可以更容易獲得更多的流場(chǎng)、溫度場(chǎng)的信息,但多數(shù)學(xué)者采用數(shù)值方法時(shí)均只考慮球床在強(qiáng)制對(duì)流情況下的流動(dòng)傳熱特性,而忽略了熱浮力的影響。 Huang 等[13]對(duì)矩形腔內(nèi)湍流混合對(duì)流特性進(jìn)行了實(shí)驗(yàn)和數(shù)值研究,分析了浮力的影響,通過(guò)與實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比,確定了 LES 方法在混合對(duì)流計(jì)算中的精度;Guardo等[14–15]采用數(shù)值方法,對(duì)隨機(jī)球床內(nèi)超臨界二氧化碳在高壓下的混合對(duì)流傳質(zhì)進(jìn)行了分析,研究表明,對(duì)于層流狀態(tài)下的超臨界流體,可以直觀地看到流量和流動(dòng)方向?qū)髻|(zhì)量的影響,在助流狀態(tài)下,會(huì)得到更大的傳質(zhì)量。研究表明,混合對(duì)流中浮力的影響不可忽略,但對(duì)球床內(nèi)流動(dòng)和傳熱的研究還不夠充分。

    本文將采用 LES 方法、結(jié)構(gòu)網(wǎng)格研究熱浮力對(duì)球床反應(yīng)堆內(nèi) FCC 面接觸單元的燃料球表面及局部流動(dòng)與傳熱的影響。

    1幾何模型與數(shù)值方法

    球的規(guī)則排列方式主要有3種[16]:簡(jiǎn)單立方( SC ),體心立方( BCC ),面心立方(FCC )。球床內(nèi)部燃料球堆積結(jié)構(gòu)可以看作是3種方式的隨機(jī)組合。兩球的接觸點(diǎn)一般有3種處理方式:人工間隙、點(diǎn)接觸、面接觸。3種規(guī)則排列方式和接觸點(diǎn)處理方式如圖1所示。由于3種排列方式中 FCC 結(jié)構(gòu)空隙率最小,而實(shí)際球床中燃料球的無(wú)規(guī)則堆積必然使球床的空隙率最小化,因此本文選擇燃料球的 FCC 排列方式,又考慮到兩球在接觸位置會(huì)存在一定程度的變形,對(duì)接觸點(diǎn)的處理方式選擇了面接觸,由此得到了一個(gè) FCC 單元??紤]到計(jì)算資源有限及為避免受進(jìn)、出口效應(yīng)影響,本文中采用4個(gè) FCC 單元堆積模型,幾何模型如圖2所示。主要參數(shù)參考 Shams 等[11]中的值。燃料球球面材料為石墨,直徑為60 mm,接觸變形部分是燃料球直徑的1%,內(nèi)部冷卻工質(zhì)為氦氣。邊界條件及參數(shù)、物性參數(shù)分別如表1、2所示。

    考慮熱浮力時(shí),密度采用Boussinesq近似,計(jì)算 Gr/Re2= g*αv?tl/v2=0.024,其中: g*為本文中采用的重力加速度, g*= gαv( Tmax一Tmin),Tmax和Tmin分別為;αv為體積膨脹系數(shù);?t為進(jìn)、出口溫差;l為特征長(zhǎng)度;v為流體進(jìn)口速度; Gr 為格拉曉夫數(shù);Re 為雷諾數(shù)。復(fù)合對(duì)流中,采用Gr/Re2作為判斷自然對(duì)流影響程度的判據(jù)。一般認(rèn)為,當(dāng)Gr/Re2≥0.01時(shí),自然對(duì)流的影響不能忽略。本文中重力方向與流動(dòng)方向相同,熱浮力方向與流動(dòng)方向相反。

    為了能夠捕捉更多的瞬時(shí)流場(chǎng)信息,采用 LES 方法進(jìn)行研究。亞格子模型采用動(dòng)態(tài)Smagorinsky?Lily 模型,其優(yōu)點(diǎn)在于具有一定的自適應(yīng)性,得出的結(jié)果更可靠。

    2網(wǎng)格適用性驗(yàn)證

    為了提高計(jì)算結(jié)果的準(zhǔn)確性,對(duì)所采用的 FCC 面接觸模型進(jìn)行結(jié)構(gòu)網(wǎng)格創(chuàng)建,采用 ICEM 軟件將流場(chǎng)區(qū)域劃分成多個(gè)六面體幾何結(jié)構(gòu)并建立對(duì)應(yīng)的六面體塊,最后通過(guò)將幾何體與對(duì)應(yīng)的六面體塊進(jìn)行關(guān)聯(lián),成功完成了結(jié)構(gòu)網(wǎng)格的劃分,并對(duì)近壁面網(wǎng)格進(jìn)行了加密。 FCC 面接觸單元網(wǎng)格如圖3所示。

    為了驗(yàn)證近壁面網(wǎng)格加密方法的適用性,對(duì) Shams 等[12]采用的 FCC 人工間隙模型,采用相同的結(jié)構(gòu)網(wǎng)格劃分方法進(jìn)行全結(jié)構(gòu)化網(wǎng)格的劃分。在相同的邊界條件和進(jìn)、出口參數(shù)條件下,采用 LES 方法進(jìn)行計(jì)算,獲取了在垂直來(lái)流方向上近球面到中心流場(chǎng)區(qū)域的時(shí)均速度分布。圖4為近壁面加密網(wǎng)格適用性驗(yàn)證結(jié)果。4(a)為提取數(shù)據(jù)位置。圖4( b)為 LES 和 q?DNS 方法結(jié)果對(duì)比,其中,橫坐標(biāo)表示的位置與縱坐標(biāo)表示的時(shí)均速度都進(jìn)行了無(wú)量綱化,并與 Shams等[14]的 q?DNS 數(shù)據(jù)進(jìn)行對(duì)比;為時(shí)均速度, U為時(shí)均速度,Umax為最大時(shí)均速度;(Z 一Zmid)為在中心流場(chǎng)區(qū)域所提取線上的數(shù)據(jù)點(diǎn) Z 相對(duì)中間位置Zmid的坐標(biāo);D 為燃料球直徑。結(jié)果顯示,與近球面處計(jì)算結(jié)果高度一致,由此說(shuō)明近壁面的網(wǎng)格加密方法適用于該球床模型的數(shù)值計(jì)算。

    為驗(yàn)證網(wǎng)格無(wú)關(guān)性,分別對(duì)網(wǎng)格數(shù)為720萬(wàn)、503萬(wàn)、325萬(wàn)的模型進(jìn)行了計(jì)算,并整理對(duì)比了進(jìn)、出口速度,溫度差。以720萬(wàn)網(wǎng)格下的計(jì)算結(jié)果為標(biāo)準(zhǔn),計(jì)算出另外兩種網(wǎng)格數(shù)下得到的數(shù)據(jù)誤差。網(wǎng)格無(wú)關(guān)性驗(yàn)證結(jié)果如表3所示,結(jié)果顯示3種網(wǎng)格數(shù)的計(jì)算結(jié)果無(wú)明顯差異,這說(shuō)明325萬(wàn)網(wǎng)格已可滿足計(jì)算要求,

    網(wǎng)格數(shù)的增加對(duì)計(jì)算結(jié)果無(wú)明顯影響。由于采用了 LES 方法,為了能更好地捕捉內(nèi)部流場(chǎng)信息,本文中用于計(jì)算的網(wǎng)格數(shù)取為503萬(wàn)。

    3計(jì)算結(jié)果與分析

    在所研究的流場(chǎng)區(qū)域沿流動(dòng)方向取相同的截面,獲取截面上的平均速度,各層截面位置及平均速度、溫度變化如圖5所示。結(jié)果顯示,在進(jìn)、出口處速度變化較大,存在明顯的進(jìn)、出口效應(yīng),擾動(dòng)較為劇烈。進(jìn)口處受影響的程度相對(duì)較大,并持續(xù)到第3層截面之后。從第4層截面開始,流速呈現(xiàn)對(duì)稱性波動(dòng),表明流動(dòng)已基本穩(wěn)定。可以發(fā)現(xiàn),自第4層截面之后,熱浮力影響效果明顯,受熱浮力影響的流體流速相對(duì)較低,說(shuō)明熱浮力的存在對(duì)流速起到抑制作用。第2層截面之后,每層溫度均有升高,相鄰兩層溫差ΔT 為7 K。為了進(jìn)一步研究流場(chǎng)內(nèi)部的流動(dòng)和傳熱情況,通過(guò)以上的分析,選取第5層和第7層截面之間的 FCC 單元作為流場(chǎng)研究區(qū)域,選取3號(hào)中心球分析球面換熱情況。

    3.1熱浮力對(duì)速度和溫度場(chǎng)的影響

    在所研究的 FCC 單元內(nèi)選取截面1、2、3,流場(chǎng)區(qū)域各截面位置及截面上流線圖如圖6所示。獲取截面1、3上的時(shí)均速度和時(shí)均溫度分布,結(jié)果如圖7、8所示??梢园l(fā)現(xiàn)熱浮力對(duì)局部流場(chǎng)速度與溫度分布影響較為明顯。在截面1的速度場(chǎng)中,受熱浮力影響,低速區(qū)域增大,尤其在燃料球底部和頂部對(duì)應(yīng)的尾流區(qū)和滯止區(qū)位置,低速區(qū)域明顯增大。溫度場(chǎng)中,低速區(qū)域的增大使高溫區(qū)域隨之增大,且在中心區(qū)域,尾流區(qū)溫度較高。截面3與接觸面相切,對(duì)比速度分布,在尾流區(qū)和滯止區(qū)同樣可以發(fā)現(xiàn)與截面1中相同的現(xiàn)象。在接觸面處熱浮力存在時(shí),流體速度會(huì)在發(fā)生分離的位置略微增加,之后的速度都會(huì)相對(duì)減小。截面上的平均溫度相對(duì)較高,接觸面處較為明顯,高溫區(qū)域延伸較長(zhǎng)。這再次說(shuō)明熱浮力的存在對(duì)流速起到抑制作用。

    為了解速度和溫度在球面和中心流場(chǎng)區(qū)域之間的變化情況,在截面1上中心區(qū)域沿流動(dòng)方向和垂直流動(dòng)方向各取一條截線,如圖9所示。提取線上的時(shí)均速度和時(shí)均溫度分布,并與不考慮熱浮力的分布進(jìn)行了對(duì)比,結(jié)果如圖10所示,其中:Vb、Tb 分別為考慮熱浮力時(shí)的時(shí)均速度和時(shí)均溫度;Vn、Tn 分別為不考慮熱浮力時(shí)的時(shí)均速度和時(shí)均溫度。分析截線1上提取到的時(shí)均速度和溫度分布(相對(duì)于中心點(diǎn)):沿 X 正方向,在近球面位置熱浮力對(duì)速度和溫度的影響并不明顯;遠(yuǎn)離球面后,熱浮力的影響開始增強(qiáng)。在熱浮力影響下,相對(duì)中心點(diǎn)的對(duì)稱位置處較低的速度被提高,較高的速度被降低,最終減少了速度分布的不對(duì)稱性,也同樣減少了溫度分布的不對(duì)稱性,截線上的溫度有所提高,考慮到速度差的減小,降低了湍流強(qiáng)度,削弱了換熱。相對(duì)于不考慮熱浮力的情況,靠近球面位置,速度和溫度的相對(duì)變化較大(速度最大降低70%,相對(duì)層間溫差,溫度最大增加43%)。遠(yuǎn)離球面后速度和溫度相對(duì)變化趨于平緩,在中心流場(chǎng)區(qū)域速度相對(duì)變化范圍在10%以內(nèi),溫度平均增加約20%。

    分析截線2上提取到的時(shí)均速度分布,沿 Y 負(fù)方向(流動(dòng)方向),不考慮熱浮力時(shí),由上一個(gè)中心球底部到下一個(gè)中心球頂部,速度先快速增大,后緩慢增大,在中間偏下的位置達(dá)到峰值,后緩慢減少,至靠近球面位置急劇降低最后接近滯止,速度受滯止區(qū)和尾流區(qū)影響明顯。在考慮熱浮力后,由速度分布可以明顯發(fā)現(xiàn),相對(duì)中心點(diǎn)呈良好的對(duì)稱性,截線上速度相對(duì)減小,峰值位置更靠近中心點(diǎn)??拷蛎嫣幩俣葎∽兊那闆r有所緩和,從一定程度上減少了由于速度變化大引起的湍流強(qiáng)度。從時(shí)均溫度分布可以發(fā)現(xiàn),考慮熱浮力后,高低溫差減小,中心區(qū)域截線上平均溫度有所增加。靠近球面位置,時(shí)均速度和溫度的相對(duì)變化較大,速度幾乎滯止(相對(duì)層間溫差,溫度最大增加50%)??拷行牧鲌?chǎng)區(qū)域,速度和溫度的相對(duì)變化有所減小,溫度最少升高約20%。

    3.2熱浮力對(duì)燃料球表面對(duì)流換熱的影響

    圖11為3號(hào)中心球表面溫度分布。由圖中可見,熱浮力的存在使球面部分位置溫度相對(duì)升高,尤其在燃料球的頂部、側(cè)面的4個(gè)接觸面中間部分以及底部較為明顯。本文中定義努塞爾數(shù) Nu = hD/λ,h = q/(T 一 Tin),其中:h為對(duì)流換熱系數(shù);λ為氦氣熱導(dǎo)率;q為熱流密度;T為球面上不同位置溫度; Tin為氦氣進(jìn)口溫度。球面Nu 分布如圖12所示。對(duì)比發(fā)現(xiàn),熱浮力存在時(shí),燃料球頂部及接觸面附近的 Nu 明顯降低,高 Nu 區(qū)域明顯減小,且差異較明顯的位置分布在燃料球的上半球面。由時(shí)均速度分布情況,頂部附近流場(chǎng)速度相對(duì)減小,減弱了對(duì)流換熱,使局部 Nu 降低。由圖6截面1、2、3上流線圖可以發(fā)現(xiàn),流場(chǎng)中的渦多集中在接觸面后沿及燃料球底部的尾流區(qū)域,擾動(dòng)較為劇烈。這說(shuō)明熱浮力的存在削弱了換熱強(qiáng)度,而渦的存在增強(qiáng)了局部擾動(dòng),削弱了熱浮力的影響。

    為進(jìn)一步了解頂部和接觸面附近的換熱情況,在球面上3個(gè)特殊位置取3條截線,截線3、4、5,取線位置如圖13所示。球面3條截線上 Nu 分布如圖14所示。通過(guò)對(duì)比每條截線上的 Nu分布可以發(fā)現(xiàn),球面頂部換熱最好,接觸面附近由于接觸面的阻礙作用 Nu較小,換熱較差。在考慮熱浮力的情況下,球面相同位置的 Nu 會(huì)相對(duì)減少,最大差值約為球面平均 Nu 的6%,即熱浮力的存在對(duì)對(duì)流換熱起到抑制作用。沿流動(dòng)方向,位于球面不同位置,熱浮力影響的程度也有所不同。由圖中可以看出,燃料球的上半部分受熱浮力的影響較大,而下半部分由于流體擾動(dòng)較強(qiáng)受熱浮力影響較小,這與分析球面 Nu 分布的結(jié)果一致。

    4結(jié)論

    采用 LES 方法,在考慮熱浮力的情況下,對(duì)球床反應(yīng)堆 FCC 面接觸結(jié)構(gòu)內(nèi)的流動(dòng)和傳熱進(jìn)行了研究。通過(guò)對(duì)比分析整個(gè)球床各層的平均速度變化,局部流場(chǎng)中速度、溫度場(chǎng)的分布及中心燃料球表面的溫度、Nu 分布,得到熱浮力對(duì)球床反應(yīng)堆局部流動(dòng)和傳熱影響的幾點(diǎn)規(guī)律:

    (1)低雷諾數(shù)下熱浮力的影響不能忽略,熱浮力的存在對(duì)球床中各層平均速度變化起到阻礙作用,會(huì)抑制速度產(chǎn)生較大的波動(dòng)。

    (2)在局部流場(chǎng)區(qū)域,時(shí)均速度、時(shí)均溫度受熱浮力影響明顯??拷蛎嫖恢茫鄬?duì)時(shí)均速度和相對(duì)時(shí)均溫度變化較大。在中心區(qū)域,熱浮力的存在會(huì)減小流場(chǎng)中速度分布的不對(duì)稱性,速度最大降低約10%,時(shí)均溫度最少升高約20%。

    (3)相對(duì)來(lái)流方向,燃料球上半部分,尤其在燃料球頂部及接觸面附近,受熱浮力影響 Nu 減小,最大降低6%,削弱了換熱。流體流過(guò)接觸面接近球底部,渦結(jié)構(gòu)增加,擾動(dòng)劇烈,削弱了熱浮力影響。

    參考文獻(xiàn):

    [1] 雷鳴澤.高溫氣冷堆產(chǎn)業(yè)推廣及應(yīng)用前景[J].中國(guó)核電, 2018, 11(1):26-29.

    [2]吳宗鑫.?我國(guó)高溫氣冷堆的發(fā)展 [J]. 核動(dòng)力工程, 2000,?21(1):39 - 43,80

    [3] 符曉銘, 王捷.高溫氣冷堆在我國(guó)的發(fā)展綜述[J].現(xiàn)代電力, 2006, 23(5):70-75.

    [4] 吳宗鑫, 肖宏才.模塊式高溫氣冷堆的安全特性[J].高技術(shù)通訊, 1994(11):34-38.

    [5] 孟現(xiàn)珂, 孫中寧, 徐廣展,等.含內(nèi)熱源堆積球床對(duì)流換熱特性的實(shí)驗(yàn)研究[J].哈爾濱工程大學(xué)學(xué)報(bào), 2012, 33(9):1122-1126.

    [6] JIA X L, GUI N, YANG X T, et al. Experimental studyof flow field characteristics on bed configurations in thepebble bed reactor[J]. Annalsof Nuclear Energy, 2017, 102:1-10.

    [7] LEE? JJ,? PARK? G? C,? KIM? K? Y,? et? al. Numericaltreatment of pebble contact in the flow and heat transferanalysis? of? a? pebble ?bed? reactor? core[J]. Nuclear Engineering and Design, 2007, 237(22):2183-2196.

    [8] FERNG Y M, LIN K Y. CFD investigation of thermal-hydraulic characteristics in a PBR core using differentcontact? treatments? between? pebbles[J]. Annals? of Nuclear Energy, 2014, 72:156-165.

    [9] SHAMS? A,? ROELOFS? F,? KOMEN? E? M? J,? et? al.Optimization of a pebble bed configuration for quasi-direct numerical simulation[J]. Nuclear Engineeringand Design, 2012, 242:331-340.

    [10] SHAMS A, ROELOFS F, KOMEN E M J. Quasi-directnumerical simulation of a pebble bed configuration. PartⅠ: flow (velocity) field? analysis[J]. Nuclear EngineeringandDesign, 2013, 263:473-489.

    [11] SHAMS? A,? ROELOFS? F,? KOMEN? E? M? J,? et? al.Numerical? simulation? of? nuclear? pebble? bedconfigurations[J]. Nuclear Engineering and Design, 2015, 290:51-64.

    [12]蔣旭, 郭雪巖.球床反應(yīng)堆流動(dòng)與傳熱的 CFD 分析:燃料球尺度[J].能源工程, 2017(6):8-13,19.

    [13] HUANG? YY,? YANG? G,? WU? J? Y. Large? eddysimulation and experimental study of turbulent mixedconvection inside a cavity with large Rayleigh number:effect of buoyancy[J]. BuildingandEnvironment, 2019, 151:268-279.

    [14] GUARDO A, COUSSIRAT M, RECASENS F, et al.CFD study on particle-to-fluid heat transfer in fixed bedreactors: convective? heat? transfer? at? low? and? highpressure[J]. Chemical? Engineering? Science, 2006, 61(13):4341-4353.

    [15] GUARDO A, COUSSIRAT M, RECASENS F, et al.CFD studies on particle-to-fluid mass and heat transferin packed beds: free convection effects in supercriticalfluids[J]. Chemical Engineering Science, 2007, 62(18-20):5503-5511.

    [16] OOMS A. Pebble flow in a high temperature reactor[R].Physics? of Nuclear? Reactor. PNR-131-2008-003. TU Delft, The Netherlands, 2008.

    91精品国产国语对白视频| 少妇裸体淫交视频免费看高清 | 久久狼人影院| 18禁裸乳无遮挡动漫免费视频| 免费观看人在逋| 99久久人妻综合| 狂野欧美激情性xxxx| av在线老鸭窝| 熟女少妇亚洲综合色aaa.| 9191精品国产免费久久| bbb黄色大片| 夫妻性生交免费视频一级片| 中文字幕高清在线视频| 久久久亚洲精品成人影院| 亚洲午夜精品一区,二区,三区| 久久久久国产一级毛片高清牌| 亚洲视频免费观看视频| 亚洲精品成人av观看孕妇| 黑人欧美特级aaaaaa片| 国产色视频综合| 国产精品.久久久| 亚洲欧美日韩另类电影网站| 中文字幕av电影在线播放| 操美女的视频在线观看| 又大又爽又粗| 日本黄色日本黄色录像| 久久久精品免费免费高清| 亚洲成人手机| 日韩制服骚丝袜av| 日韩熟女老妇一区二区性免费视频| 精品亚洲成国产av| 国产片特级美女逼逼视频| 国产精品欧美亚洲77777| 在线观看免费日韩欧美大片| 亚洲av日韩精品久久久久久密 | 一区二区三区激情视频| 99国产精品99久久久久| 女警被强在线播放| 免费高清在线观看视频在线观看| 国产伦人伦偷精品视频| 国产伦理片在线播放av一区| 男男h啪啪无遮挡| 国产精品免费大片| 少妇人妻久久综合中文| 性高湖久久久久久久久免费观看| 欧美+亚洲+日韩+国产| 超碰成人久久| 日日夜夜操网爽| 日本午夜av视频| 男女免费视频国产| 97在线人人人人妻| 女人爽到高潮嗷嗷叫在线视频| 国产精品一区二区在线不卡| 亚洲色图 男人天堂 中文字幕| 成年人免费黄色播放视频| 免费黄频网站在线观看国产| 精品少妇一区二区三区视频日本电影| 在线观看免费高清a一片| 精品熟女少妇八av免费久了| 男女免费视频国产| 亚洲男人天堂网一区| 丝袜在线中文字幕| 色播在线永久视频| 欧美大码av| 男女国产视频网站| 美女脱内裤让男人舔精品视频| 午夜福利一区二区在线看| 国产爽快片一区二区三区| 2021少妇久久久久久久久久久| 久久久久久久久免费视频了| 国产成人av教育| 一级片免费观看大全| 久久久久精品人妻al黑| 亚洲专区中文字幕在线| 女性被躁到高潮视频| 亚洲av欧美aⅴ国产| 天天添夜夜摸| 老司机影院毛片| 亚洲人成网站在线观看播放| 三上悠亚av全集在线观看| 欧美性长视频在线观看| 中国国产av一级| 多毛熟女@视频| 久久性视频一级片| 久久国产精品大桥未久av| 国产精品一区二区在线观看99| 波多野结衣av一区二区av| 激情视频va一区二区三区| 老司机在亚洲福利影院| 久久青草综合色| xxx大片免费视频| 两性夫妻黄色片| 国产熟女欧美一区二区| 免费看不卡的av| 亚洲国产毛片av蜜桃av| 在线亚洲精品国产二区图片欧美| 2021少妇久久久久久久久久久| 青春草亚洲视频在线观看| 搡老乐熟女国产| 国产精品 欧美亚洲| 黄色毛片三级朝国网站| 女人精品久久久久毛片| 亚洲一区二区三区欧美精品| 天天添夜夜摸| 国产伦理片在线播放av一区| 国产又色又爽无遮挡免| 国产av一区二区精品久久| 少妇精品久久久久久久| 久久人人97超碰香蕉20202| 国产日韩欧美视频二区| 99国产精品一区二区三区| 日日夜夜操网爽| 国产精品.久久久| 国产91精品成人一区二区三区 | 岛国毛片在线播放| 午夜福利乱码中文字幕| 国产爽快片一区二区三区| 麻豆乱淫一区二区| 午夜91福利影院| 啦啦啦中文免费视频观看日本| 丰满迷人的少妇在线观看| 国产精品亚洲av一区麻豆| 制服诱惑二区| 母亲3免费完整高清在线观看| 国产成人精品久久久久久| 午夜免费男女啪啪视频观看| 天天影视国产精品| 晚上一个人看的免费电影| 亚洲欧美激情在线| 亚洲精品日本国产第一区| 50天的宝宝边吃奶边哭怎么回事| 91麻豆av在线| 性少妇av在线| 熟女少妇亚洲综合色aaa.| 老司机午夜十八禁免费视频| 狠狠婷婷综合久久久久久88av| 91精品三级在线观看| 久久毛片免费看一区二区三区| 亚洲av成人不卡在线观看播放网 | 久久青草综合色| 国产精品一区二区免费欧美 | e午夜精品久久久久久久| 操出白浆在线播放| 久9热在线精品视频| 男女边摸边吃奶| 精品免费久久久久久久清纯 | 91精品三级在线观看| 美女午夜性视频免费| 亚洲图色成人| 嫩草影视91久久| 亚洲 国产 在线| 飞空精品影院首页| 一边摸一边做爽爽视频免费| 欧美精品亚洲一区二区| 五月天丁香电影| 久久精品国产综合久久久| 午夜日韩欧美国产| 一级毛片我不卡| 老司机在亚洲福利影院| 一边摸一边做爽爽视频免费| 妹子高潮喷水视频| 老司机靠b影院| 亚洲国产欧美日韩在线播放| 久久久久久久大尺度免费视频| 老汉色∧v一级毛片| 中文字幕制服av| 欧美亚洲日本最大视频资源| 午夜视频精品福利| 十八禁高潮呻吟视频| 自线自在国产av| 99久久人妻综合| 久久人人97超碰香蕉20202| 成年av动漫网址| 亚洲 国产 在线| 亚洲国产精品国产精品| 蜜桃国产av成人99| 又紧又爽又黄一区二区| 老司机靠b影院| 成年人午夜在线观看视频| 日韩免费高清中文字幕av| 色婷婷久久久亚洲欧美| 国产一级毛片在线| av国产久精品久网站免费入址| 9色porny在线观看| 久9热在线精品视频| 别揉我奶头~嗯~啊~动态视频 | 夫妻性生交免费视频一级片| 国产在线免费精品| 亚洲av日韩精品久久久久久密 | 久久性视频一级片| 亚洲熟女精品中文字幕| 青草久久国产| 观看av在线不卡| 亚洲国产av影院在线观看| 桃花免费在线播放| 首页视频小说图片口味搜索 | 老熟女久久久| avwww免费| 国产精品国产av在线观看| 国产精品二区激情视频| 美女高潮到喷水免费观看| 成人免费观看视频高清| 一本一本久久a久久精品综合妖精| 国产又色又爽无遮挡免| 赤兔流量卡办理| 久久精品国产亚洲av涩爱| 色94色欧美一区二区| cao死你这个sao货| 日韩av不卡免费在线播放| 亚洲精品自拍成人| 欧美人与性动交α欧美精品济南到| 后天国语完整版免费观看| 日韩伦理黄色片| 国产精品秋霞免费鲁丝片| 久久久久久亚洲精品国产蜜桃av| 国产97色在线日韩免费| 色综合欧美亚洲国产小说| 熟女少妇亚洲综合色aaa.| 搡老乐熟女国产| 免费一级毛片在线播放高清视频 | 黄片播放在线免费| 午夜福利视频精品| 国产精品久久久久久精品电影小说| 一级黄色大片毛片| 99久久人妻综合| 国产欧美日韩一区二区三 | 丝袜在线中文字幕| 国产xxxxx性猛交| 久久久久精品人妻al黑| 国产女主播在线喷水免费视频网站| 久久精品亚洲av国产电影网| 性少妇av在线| www.精华液| 啦啦啦中文免费视频观看日本| 中文字幕av电影在线播放| 久久国产亚洲av麻豆专区| 精品国产国语对白av| 亚洲五月色婷婷综合| 日韩中文字幕欧美一区二区 | 欧美激情高清一区二区三区| 永久免费av网站大全| 久久免费观看电影| 大香蕉久久成人网| 国产成人av教育| 多毛熟女@视频| 精品一区二区三区av网在线观看 | 亚洲精品国产av成人精品| 九色亚洲精品在线播放| 亚洲精品国产区一区二| av天堂在线播放| 午夜福利视频在线观看免费| 日韩电影二区| 亚洲色图综合在线观看| 日本欧美国产在线视频| 9热在线视频观看99| 国产日韩欧美视频二区| 免费日韩欧美在线观看| 国产精品免费视频内射| 久久人人97超碰香蕉20202| 亚洲精品一区蜜桃| 国产又爽黄色视频| 免费在线观看视频国产中文字幕亚洲 | 日韩一区二区三区影片| 在线观看www视频免费| 亚洲三区欧美一区| 亚洲国产欧美日韩在线播放| 丰满饥渴人妻一区二区三| 嫩草影视91久久| 18禁观看日本| 纯流量卡能插随身wifi吗| 国产爽快片一区二区三区| 亚洲精品一二三| 日本欧美视频一区| 日韩电影二区| 亚洲国产看品久久| 一区二区日韩欧美中文字幕| 91麻豆精品激情在线观看国产 | 国产精品.久久久| 日本黄色日本黄色录像| 看十八女毛片水多多多| 亚洲成色77777| 欧美日韩亚洲综合一区二区三区_| 大话2 男鬼变身卡| 久久99一区二区三区| 肉色欧美久久久久久久蜜桃| 日韩av在线免费看完整版不卡| 久久性视频一级片| 极品人妻少妇av视频| 国产精品欧美亚洲77777| 制服人妻中文乱码| 脱女人内裤的视频| 侵犯人妻中文字幕一二三四区| 麻豆乱淫一区二区| 伊人久久大香线蕉亚洲五| 成人午夜精彩视频在线观看| 午夜影院在线不卡| 久久久国产精品麻豆| 最近中文字幕2019免费版| 色94色欧美一区二区| 亚洲一码二码三码区别大吗| 只有这里有精品99| 免费不卡黄色视频| 黄色a级毛片大全视频| 国产成人免费无遮挡视频| 高清不卡的av网站| av电影中文网址| 国产精品熟女久久久久浪| 日韩免费高清中文字幕av| 免费在线观看影片大全网站 | 亚洲中文字幕日韩| 好男人电影高清在线观看| 日本91视频免费播放| 国产成人系列免费观看| 国产又色又爽无遮挡免| 在线精品无人区一区二区三| 国产一区二区三区综合在线观看| 日韩,欧美,国产一区二区三区| 国产欧美亚洲国产| 国产精品久久久久久人妻精品电影 | 久久精品国产亚洲av涩爱| 99精国产麻豆久久婷婷| 亚洲精品日韩在线中文字幕| 精品国产乱码久久久久久小说| 国产欧美日韩综合在线一区二区| 国产麻豆69| 欧美黄色淫秽网站| 熟女少妇亚洲综合色aaa.| 久久久精品94久久精品| 久久 成人 亚洲| 久久精品亚洲熟妇少妇任你| 一级黄片播放器| 欧美日韩黄片免| 男人添女人高潮全过程视频| 啦啦啦 在线观看视频| 男女下面插进去视频免费观看| 制服诱惑二区| 中文字幕色久视频| 男女边吃奶边做爰视频| 人人澡人人妻人| 咕卡用的链子| 少妇 在线观看| 亚洲成人免费av在线播放| 久久毛片免费看一区二区三区| 黄色视频不卡| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品久久久久久婷婷小说| 黑人猛操日本美女一级片| 热re99久久国产66热| 中文乱码字字幕精品一区二区三区| 亚洲专区国产一区二区| 好男人电影高清在线观看| 国产成人系列免费观看| 尾随美女入室| 日韩av在线免费看完整版不卡| 麻豆国产av国片精品| 国产男女内射视频| 日日夜夜操网爽| 亚洲国产精品成人久久小说| 精品亚洲成a人片在线观看| 亚洲国产精品999| 成年美女黄网站色视频大全免费| 视频在线观看一区二区三区| 亚洲欧美成人综合另类久久久| 久久午夜综合久久蜜桃| 一区二区三区激情视频| 欧美在线黄色| 午夜福利免费观看在线| 国产女主播在线喷水免费视频网站| 最新的欧美精品一区二区| 亚洲欧洲日产国产| 国产片特级美女逼逼视频| 亚洲 国产 在线| 91麻豆精品激情在线观看国产 | 欧美xxⅹ黑人| 狠狠婷婷综合久久久久久88av| 欧美日韩福利视频一区二区| 美女扒开内裤让男人捅视频| 亚洲,一卡二卡三卡| 老汉色∧v一级毛片| tube8黄色片| 久久久欧美国产精品| a级毛片在线看网站| 99热网站在线观看| 50天的宝宝边吃奶边哭怎么回事| 亚洲欧美日韩另类电影网站| 国产精品久久久av美女十八| 男女床上黄色一级片免费看| 久久99一区二区三区| av又黄又爽大尺度在线免费看| 91精品国产国语对白视频| 亚洲五月色婷婷综合| 国产精品久久久人人做人人爽| 久久人妻熟女aⅴ| 女警被强在线播放| 男男h啪啪无遮挡| 丰满饥渴人妻一区二区三| 人人妻,人人澡人人爽秒播 | 人体艺术视频欧美日本| 电影成人av| 亚洲国产中文字幕在线视频| 国产精品熟女久久久久浪| 亚洲第一青青草原| 男女之事视频高清在线观看 | 国产精品偷伦视频观看了| 国产精品一区二区在线不卡| 99国产精品一区二区三区| 在线观看www视频免费| 人人妻人人爽人人添夜夜欢视频| 欧美人与性动交α欧美精品济南到| 国产精品一区二区免费欧美 | 天天躁日日躁夜夜躁夜夜| 欧美+亚洲+日韩+国产| a级片在线免费高清观看视频| 波野结衣二区三区在线| 91国产中文字幕| 黑人欧美特级aaaaaa片| 欧美日韩av久久| 日本91视频免费播放| 国产精品久久久久久人妻精品电影 | 久久99一区二区三区| 亚洲国产av影院在线观看| 午夜日韩欧美国产| 国产精品久久久久久精品古装| 久久精品熟女亚洲av麻豆精品| 国产三级黄色录像| 日韩av不卡免费在线播放| 性少妇av在线| 午夜福利视频精品| 亚洲欧美日韩另类电影网站| 亚洲国产日韩一区二区| 在线观看免费午夜福利视频| 青青草视频在线视频观看| 亚洲欧美日韩高清在线视频 | 天天影视国产精品| 女性生殖器流出的白浆| 最新在线观看一区二区三区 | 美女扒开内裤让男人捅视频| 丰满少妇做爰视频| 免费少妇av软件| 国产欧美日韩一区二区三 | 一边摸一边抽搐一进一出视频| 伊人久久大香线蕉亚洲五| 国产精品二区激情视频| 最黄视频免费看| 亚洲综合色网址| 黄色毛片三级朝国网站| 久久久久久人人人人人| 人人妻人人澡人人爽人人夜夜| bbb黄色大片| 久久 成人 亚洲| 久久人人97超碰香蕉20202| 国产成人一区二区在线| videos熟女内射| 一级a爱视频在线免费观看| 色播在线永久视频| 老司机靠b影院| 午夜免费男女啪啪视频观看| 亚洲精品自拍成人| 丰满人妻熟妇乱又伦精品不卡| 久热爱精品视频在线9| 亚洲第一青青草原| 免费观看a级毛片全部| 精品久久久精品久久久| 涩涩av久久男人的天堂| kizo精华| 久久九九热精品免费| 久久 成人 亚洲| 欧美日韩黄片免| 欧美xxⅹ黑人| 美女扒开内裤让男人捅视频| 色婷婷久久久亚洲欧美| a级片在线免费高清观看视频| 欧美成狂野欧美在线观看| 国产1区2区3区精品| 久久久久久久国产电影| 亚洲成色77777| 男女高潮啪啪啪动态图| 一级黄片播放器| 91麻豆精品激情在线观看国产 | 日本av免费视频播放| 成人影院久久| 少妇 在线观看| 午夜福利免费观看在线| 99热网站在线观看| 啦啦啦视频在线资源免费观看| 亚洲精品国产av蜜桃| 国产成人精品在线电影| 精品一区在线观看国产| 欧美日韩黄片免| 国产精品九九99| 成人18禁高潮啪啪吃奶动态图| 三上悠亚av全集在线观看| 国产成人精品久久久久久| 黄色毛片三级朝国网站| 国产日韩欧美亚洲二区| 成在线人永久免费视频| 黄色a级毛片大全视频| 在线观看人妻少妇| 国产成人精品久久久久久| 狂野欧美激情性xxxx| 欧美成人精品欧美一级黄| 国产在线免费精品| 亚洲伊人久久精品综合| 一区二区三区激情视频| 成人免费观看视频高清| 久久热在线av| 亚洲 国产 在线| 少妇裸体淫交视频免费看高清 | 国产成人欧美| 考比视频在线观看| e午夜精品久久久久久久| 手机成人av网站| 久久久国产欧美日韩av| 国产熟女午夜一区二区三区| 咕卡用的链子| 久久av网站| 亚洲欧洲日产国产| 国产日韩欧美亚洲二区| e午夜精品久久久久久久| 狠狠精品人妻久久久久久综合| 视频在线观看一区二区三区| 亚洲 欧美一区二区三区| 精品人妻一区二区三区麻豆| 欧美成狂野欧美在线观看| netflix在线观看网站| 丁香六月天网| 久9热在线精品视频| 咕卡用的链子| 国产成人一区二区在线| 日本vs欧美在线观看视频| 手机成人av网站| 最近中文字幕2019免费版| 国产精品免费大片| 一边摸一边做爽爽视频免费| 美女脱内裤让男人舔精品视频| www.熟女人妻精品国产| 欧美日韩亚洲高清精品| 麻豆av在线久日| 色播在线永久视频| 久久天堂一区二区三区四区| 国产成人91sexporn| 九色亚洲精品在线播放| 欧美日韩综合久久久久久| 多毛熟女@视频| 91字幕亚洲| 亚洲欧美精品自产自拍| 亚洲欧洲精品一区二区精品久久久| 精品亚洲成国产av| 日韩大码丰满熟妇| 王馨瑶露胸无遮挡在线观看| 久久精品熟女亚洲av麻豆精品| 青草久久国产| 777久久人妻少妇嫩草av网站| 亚洲成av片中文字幕在线观看| 日韩制服骚丝袜av| 国产欧美日韩一区二区三区在线| 免费观看av网站的网址| 日本猛色少妇xxxxx猛交久久| 国产精品麻豆人妻色哟哟久久| cao死你这个sao货| 看免费av毛片| 另类亚洲欧美激情| 亚洲av成人精品一二三区| 亚洲欧美中文字幕日韩二区| 中文字幕另类日韩欧美亚洲嫩草| 精品第一国产精品| 中文字幕亚洲精品专区| 久久性视频一级片| 我要看黄色一级片免费的| 男女床上黄色一级片免费看| 久久影院123| 国产一区有黄有色的免费视频| 精品熟女少妇八av免费久了| 亚洲欧美中文字幕日韩二区| 51午夜福利影视在线观看| 91精品国产国语对白视频| videosex国产| 青春草亚洲视频在线观看| 国产日韩欧美视频二区| 伊人亚洲综合成人网| 美女福利国产在线| 一区福利在线观看| 2018国产大陆天天弄谢| 狠狠精品人妻久久久久久综合| 亚洲av男天堂| 国产深夜福利视频在线观看| 国产熟女午夜一区二区三区| 99热全是精品| 叶爱在线成人免费视频播放| 丰满人妻熟妇乱又伦精品不卡| 久久精品国产a三级三级三级| 国产成人av教育| 国产亚洲一区二区精品| 人体艺术视频欧美日本| 咕卡用的链子| 99国产精品免费福利视频| 午夜激情久久久久久久| 精品亚洲成国产av| 欧美精品一区二区免费开放| 日韩免费高清中文字幕av| 丝袜美腿诱惑在线| 国产99久久九九免费精品| 精品视频人人做人人爽| 一边亲一边摸免费视频| 亚洲 国产 在线| 久久99精品国语久久久| 久热爱精品视频在线9| 香蕉丝袜av| 99国产精品一区二区三区| 看十八女毛片水多多多| av天堂久久9| 亚洲,欧美精品.| 欧美国产精品一级二级三级| 精品人妻在线不人妻| 一本色道久久久久久精品综合|