• <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.

    一区二区三区激情视频| 欧美成人一区二区免费高清观看 | 99久久无色码亚洲精品果冻| 99久久国产精品久久久| 精品一区二区三区av网在线观看| 国产精品99久久99久久久不卡| 男人舔女人的私密视频| 亚洲精品粉嫩美女一区| 在线观看日韩欧美| 99精品在免费线老司机午夜| 女同久久另类99精品国产91| 男人舔女人下体高潮全视频| 香蕉丝袜av| 一进一出抽搐动态| 亚洲中文日韩欧美视频| 大型黄色视频在线免费观看| 一二三四社区在线视频社区8| 色播亚洲综合网| 久久久国产成人精品二区| 国产激情欧美一区二区| 国产爱豆传媒在线观看 | 午夜福利视频1000在线观看| 精品午夜福利视频在线观看一区| 韩国av一区二区三区四区| 伊人久久大香线蕉亚洲五| 99在线人妻在线中文字幕| 一区二区三区精品91| 成人18禁在线播放| 黄网站色视频无遮挡免费观看| 久久精品aⅴ一区二区三区四区| 亚洲av第一区精品v没综合| 成年免费大片在线观看| 欧美国产精品va在线观看不卡| 自线自在国产av| 女生性感内裤真人,穿戴方法视频| 成年人黄色毛片网站| 亚洲天堂国产精品一区在线| 国产精品1区2区在线观看.| 美女午夜性视频免费| 国产黄色小视频在线观看| 亚洲成a人片在线一区二区| 9191精品国产免费久久| 在线视频色国产色| 欧美成人性av电影在线观看| 少妇的丰满在线观看| 香蕉久久夜色| 国产又爽黄色视频| 亚洲国产欧美一区二区综合| 熟女少妇亚洲综合色aaa.| 日韩国内少妇激情av| 2021天堂中文幕一二区在线观 | 激情在线观看视频在线高清| 国产真人三级小视频在线观看| 日本免费a在线| 国产单亲对白刺激| 手机成人av网站| 99re在线观看精品视频| 一级片免费观看大全| 亚洲av电影不卡..在线观看| 欧美日韩精品网址| 国产一卡二卡三卡精品| 在线观看免费午夜福利视频| 日日爽夜夜爽网站| 最好的美女福利视频网| 少妇裸体淫交视频免费看高清 | 欧美大码av| 麻豆成人av在线观看| 久久婷婷人人爽人人干人人爱| 亚洲国产精品成人综合色| 成人三级做爰电影| 又黄又爽又免费观看的视频| 国内揄拍国产精品人妻在线 | 国产熟女xx| 人人妻人人澡欧美一区二区| 国产黄色小视频在线观看| 精品午夜福利视频在线观看一区| 曰老女人黄片| 亚洲,欧美精品.| 两性夫妻黄色片| 久久青草综合色| 老熟妇仑乱视频hdxx| 亚洲精品粉嫩美女一区| 欧美中文日本在线观看视频| www日本在线高清视频| 国产一区在线观看成人免费| 母亲3免费完整高清在线观看| 欧洲精品卡2卡3卡4卡5卡区| 精品少妇一区二区三区视频日本电影| 国产免费av片在线观看野外av| 嫩草影院精品99| 在线免费观看的www视频| 国产成人欧美| 国产精品综合久久久久久久免费| 亚洲专区中文字幕在线| 色在线成人网| 热re99久久国产66热| 国产1区2区3区精品| 欧美 亚洲 国产 日韩一| 亚洲专区国产一区二区| 欧美一区二区精品小视频在线| 欧美日韩黄片免| 亚洲精品一区av在线观看| 亚洲精品久久成人aⅴ小说| videosex国产| 一个人观看的视频www高清免费观看 | 夜夜爽天天搞| 精品一区二区三区av网在线观看| 12—13女人毛片做爰片一| 国产精品九九99| 欧美色欧美亚洲另类二区| 国产成+人综合+亚洲专区| 成人三级黄色视频| 18禁裸乳无遮挡免费网站照片 | 黄色女人牲交| 别揉我奶头~嗯~啊~动态视频| 丁香六月欧美| 午夜免费成人在线视频| 亚洲一码二码三码区别大吗| 校园春色视频在线观看| 亚洲一区二区三区色噜噜| 在线天堂中文资源库| 日本撒尿小便嘘嘘汇集6| 免费在线观看视频国产中文字幕亚洲| 男女下面进入的视频免费午夜 | 日韩精品中文字幕看吧| 国产一卡二卡三卡精品| 国产视频内射| 中文亚洲av片在线观看爽| 国产精品一区二区三区四区久久 | 激情在线观看视频在线高清| 亚洲,欧美精品.| 亚洲一卡2卡3卡4卡5卡精品中文| 国产成人精品久久二区二区免费| 国产在线观看jvid| 给我免费播放毛片高清在线观看| 日日夜夜操网爽| 亚洲精品国产区一区二| 欧美日韩中文字幕国产精品一区二区三区| 色av中文字幕| 少妇粗大呻吟视频| 亚洲av熟女| 亚洲av电影不卡..在线观看| 在线观看66精品国产| 国产精品,欧美在线| 久9热在线精品视频| 久久精品国产综合久久久| 日韩精品免费视频一区二区三区| 嫁个100分男人电影在线观看| 午夜免费鲁丝| 国产精品久久久av美女十八| 欧美成人一区二区免费高清观看 | 天堂影院成人在线观看| 波多野结衣高清无吗| 国产黄色小视频在线观看| 亚洲欧洲精品一区二区精品久久久| 999久久久精品免费观看国产| 91成年电影在线观看| 性欧美人与动物交配| 在线观看66精品国产| 国产一区二区在线av高清观看| 国产精品亚洲av一区麻豆| 侵犯人妻中文字幕一二三四区| 亚洲国产欧洲综合997久久, | 男女做爰动态图高潮gif福利片| 大香蕉久久成人网| 久久精品亚洲精品国产色婷小说| 男女做爰动态图高潮gif福利片| 麻豆国产av国片精品| 美国免费a级毛片| 欧美黑人精品巨大| 午夜激情福利司机影院| 麻豆成人av在线观看| av有码第一页| 日韩欧美国产在线观看| 女人高潮潮喷娇喘18禁视频| 亚洲五月天丁香| 免费观看精品视频网站| av欧美777| 99在线视频只有这里精品首页| 国产伦一二天堂av在线观看| 国产精品久久久人人做人人爽| 久9热在线精品视频| 久久精品亚洲精品国产色婷小说| 91麻豆精品激情在线观看国产| 18禁黄网站禁片午夜丰满| 又黄又爽又免费观看的视频| 久久精品影院6| 成年女人毛片免费观看观看9| 婷婷精品国产亚洲av| 久久香蕉国产精品| 正在播放国产对白刺激| 波多野结衣高清无吗| 老鸭窝网址在线观看| 国产99久久九九免费精品| 国产成人系列免费观看| 亚洲精品久久国产高清桃花| 啦啦啦观看免费观看视频高清| 久久精品91蜜桃| 久久精品aⅴ一区二区三区四区| 久久中文字幕一级| 操出白浆在线播放| 成人一区二区视频在线观看| 欧美日韩黄片免| 色老头精品视频在线观看| 色老头精品视频在线观看| 一级a爱片免费观看的视频| 淫秽高清视频在线观看| 日本一区二区免费在线视频| 日日夜夜操网爽| 亚洲 欧美 日韩 在线 免费| 男女做爰动态图高潮gif福利片| 欧美黑人巨大hd| 法律面前人人平等表现在哪些方面| 老熟妇乱子伦视频在线观看| 男女下面进入的视频免费午夜 | 久久欧美精品欧美久久欧美| 国产精品av久久久久免费| 久久这里只有精品19| 窝窝影院91人妻| 欧美av亚洲av综合av国产av| 后天国语完整版免费观看| 色婷婷久久久亚洲欧美| 久久久国产精品麻豆| 国内久久婷婷六月综合欲色啪| 好男人电影高清在线观看| 亚洲国产欧美一区二区综合| 9191精品国产免费久久| 国产私拍福利视频在线观看| 非洲黑人性xxxx精品又粗又长| 在线看三级毛片| 欧美国产精品va在线观看不卡| 欧美黑人欧美精品刺激| 亚洲av电影不卡..在线观看| ponron亚洲| a级毛片在线看网站| 一级毛片高清免费大全| 国产精品二区激情视频| 美女扒开内裤让男人捅视频| 亚洲精品国产区一区二| 国产精品日韩av在线免费观看| 99热6这里只有精品| 国产亚洲精品久久久久久毛片| 欧美日韩亚洲综合一区二区三区_| 国产亚洲欧美精品永久| 久久欧美精品欧美久久欧美| 青草久久国产| 国产黄片美女视频| 嫩草影院精品99| 成人18禁高潮啪啪吃奶动态图| 成人三级黄色视频| 中文字幕人妻熟女乱码| 日本成人三级电影网站| 亚洲欧美一区二区三区黑人| 国产av在哪里看| 久久久久九九精品影院| 老熟妇乱子伦视频在线观看| 国产精品自产拍在线观看55亚洲| 久久热在线av| 久久99热这里只有精品18| 国产视频内射| 亚洲av成人不卡在线观看播放网| 俺也久久电影网| 国产伦在线观看视频一区| 巨乳人妻的诱惑在线观看| 久久午夜亚洲精品久久| 老熟妇乱子伦视频在线观看| 免费看美女性在线毛片视频| 国产精品野战在线观看| 99精品欧美一区二区三区四区| 久久久国产成人精品二区| 一进一出抽搐gif免费好疼| 亚洲av熟女| 亚洲av成人一区二区三| 99精品欧美一区二区三区四区| 熟女少妇亚洲综合色aaa.| 国产av在哪里看| 麻豆国产av国片精品| 婷婷精品国产亚洲av| 91av网站免费观看| 亚洲av成人一区二区三| videosex国产| 欧美一级a爱片免费观看看 | 亚洲aⅴ乱码一区二区在线播放 | 国产精品久久久av美女十八| 久久这里只有精品19| 国内毛片毛片毛片毛片毛片| 国产精品野战在线观看| 久久久水蜜桃国产精品网| 国产高清videossex| 亚洲欧美日韩无卡精品| 国产97色在线日韩免费| 久久久久久九九精品二区国产 | 精品一区二区三区av网在线观看| 黄色女人牲交| 国内久久婷婷六月综合欲色啪| 久久久国产欧美日韩av| 国产亚洲欧美98| 午夜福利高清视频| 一区二区三区精品91| 777久久人妻少妇嫩草av网站| 叶爱在线成人免费视频播放| 黄色毛片三级朝国网站| 国产精品九九99| 亚洲欧美日韩无卡精品| 久久久久久亚洲精品国产蜜桃av| www.自偷自拍.com| 亚洲人成电影免费在线| 99久久久亚洲精品蜜臀av| 视频在线观看一区二区三区| 亚洲色图av天堂| 听说在线观看完整版免费高清| 可以在线观看毛片的网站| 免费在线观看完整版高清| 丰满人妻熟妇乱又伦精品不卡| 亚洲第一av免费看| 特大巨黑吊av在线直播 | 亚洲欧美精品综合久久99| av中文乱码字幕在线| 久久香蕉激情| 国内少妇人妻偷人精品xxx网站 | 亚洲精品久久国产高清桃花| 亚洲av熟女| 一本大道久久a久久精品| 欧美日韩中文字幕国产精品一区二区三区| 看片在线看免费视频| 久久精品影院6| 国产精品一区二区免费欧美| 国产私拍福利视频在线观看| 精品久久久久久久久久久久久 | 在线av久久热| 国产av在哪里看| 国产成人一区二区三区免费视频网站| 亚洲一区二区三区色噜噜| 午夜福利在线观看吧| 欧美精品啪啪一区二区三区| 18禁美女被吸乳视频| 国产熟女xx| 91成年电影在线观看| 看免费av毛片| 99在线人妻在线中文字幕| 欧美日韩亚洲国产一区二区在线观看| 国产不卡一卡二| 特大巨黑吊av在线直播 | 亚洲国产中文字幕在线视频| 精品国产一区二区三区四区第35| 女人被狂操c到高潮| 日本五十路高清| 亚洲人成伊人成综合网2020| 亚洲精品在线观看二区| 欧美丝袜亚洲另类 | 看免费av毛片| 亚洲中文av在线| 日日摸夜夜添夜夜添小说| 天堂动漫精品| 久久青草综合色| 老汉色av国产亚洲站长工具| 99精品在免费线老司机午夜| 久久久久国产精品人妻aⅴ院| 欧美性猛交╳xxx乱大交人| 怎么达到女性高潮| 亚洲精品粉嫩美女一区| 色综合婷婷激情| a级毛片a级免费在线| 看黄色毛片网站| 色综合站精品国产| 欧美乱码精品一区二区三区| 欧美日韩亚洲综合一区二区三区_| 在线观看舔阴道视频| 嫩草影院精品99| 亚洲国产日韩欧美精品在线观看 | 激情在线观看视频在线高清| 啪啪无遮挡十八禁网站| 精品不卡国产一区二区三区| 丝袜人妻中文字幕| 黄网站色视频无遮挡免费观看| 午夜免费激情av| 精品不卡国产一区二区三区| 亚洲五月婷婷丁香| 亚洲成人久久性| 在线观看免费午夜福利视频| www.自偷自拍.com| 女人被狂操c到高潮| 亚洲国产精品sss在线观看| 老司机午夜福利在线观看视频| 丁香欧美五月| 中文字幕另类日韩欧美亚洲嫩草| 99久久99久久久精品蜜桃| 嫩草影院精品99| 国产熟女xx| 欧美日韩乱码在线| 老司机在亚洲福利影院| 黄片大片在线免费观看| 久久久久久人人人人人| 欧美色欧美亚洲另类二区| 亚洲精品一区av在线观看| 欧美黄色淫秽网站| 男女下面进入的视频免费午夜 | 亚洲国产毛片av蜜桃av| 少妇的丰满在线观看| 美女免费视频网站| 国产精品永久免费网站| 亚洲成av片中文字幕在线观看| 看免费av毛片| 嫩草影院精品99| 成人精品一区二区免费| 亚洲九九香蕉| 日本熟妇午夜| а√天堂www在线а√下载| 日韩精品中文字幕看吧| 18禁裸乳无遮挡免费网站照片 | 最近最新免费中文字幕在线| 亚洲人成网站在线播放欧美日韩| 国产精华一区二区三区| 99久久久亚洲精品蜜臀av| 人人妻人人澡人人看| 神马国产精品三级电影在线观看 | 欧美乱色亚洲激情| 99国产精品一区二区蜜桃av| 欧美人与性动交α欧美精品济南到| 老熟妇仑乱视频hdxx| 99热6这里只有精品| 无遮挡黄片免费观看| 欧美一级毛片孕妇| 在线观看免费午夜福利视频| 一本一本综合久久| 成人永久免费在线观看视频| 男女下面进入的视频免费午夜 | 91麻豆av在线| 婷婷精品国产亚洲av在线| 十八禁人妻一区二区| 国产免费男女视频| 啦啦啦 在线观看视频| 一本精品99久久精品77| 精品国产一区二区三区四区第35| 很黄的视频免费| 老汉色∧v一级毛片| 又大又爽又粗| 亚洲一区中文字幕在线| 国内毛片毛片毛片毛片毛片| 麻豆一二三区av精品| 亚洲成av片中文字幕在线观看| 在线播放国产精品三级| 国产在线观看jvid| x7x7x7水蜜桃| 国产精品综合久久久久久久免费| 黄片大片在线免费观看| 国产精品98久久久久久宅男小说| 一夜夜www| 黑人操中国人逼视频| 亚洲av成人av| 看片在线看免费视频| 制服丝袜大香蕉在线| 老司机在亚洲福利影院| 一本一本综合久久| 国产三级在线视频| 国产精品永久免费网站| 亚洲精品粉嫩美女一区| 两性午夜刺激爽爽歪歪视频在线观看 | 大型av网站在线播放| 变态另类成人亚洲欧美熟女| av福利片在线| 97人妻精品一区二区三区麻豆 | 精华霜和精华液先用哪个| www.熟女人妻精品国产| 长腿黑丝高跟| 国产1区2区3区精品| 国产黄色小视频在线观看| 久热爱精品视频在线9| 给我免费播放毛片高清在线观看| 午夜免费成人在线视频| svipshipincom国产片| 亚洲精品色激情综合| 在线国产一区二区在线| 搡老妇女老女人老熟妇| 夜夜看夜夜爽夜夜摸| 美女 人体艺术 gogo| 国产成人系列免费观看| 久久热在线av| 天天添夜夜摸| 免费无遮挡裸体视频| 亚洲成av人片免费观看| 99国产极品粉嫩在线观看| 成熟少妇高潮喷水视频| 99热只有精品国产| 国产成人精品久久二区二区免费| 国产精品久久久久久人妻精品电影| 激情在线观看视频在线高清| 91国产中文字幕| 老熟妇仑乱视频hdxx| 可以免费在线观看a视频的电影网站| 丰满的人妻完整版| 亚洲中文日韩欧美视频| 久久亚洲真实| 亚洲精品国产区一区二| 成人国语在线视频| 免费av毛片视频| 50天的宝宝边吃奶边哭怎么回事| 90打野战视频偷拍视频| 亚洲精品一卡2卡三卡4卡5卡| 后天国语完整版免费观看| 亚洲美女黄片视频| 黄色丝袜av网址大全| 女性生殖器流出的白浆| 婷婷亚洲欧美| 欧美色视频一区免费| 男女午夜视频在线观看| 亚洲天堂国产精品一区在线| 免费高清在线观看日韩| 侵犯人妻中文字幕一二三四区| 女警被强在线播放| 午夜福利欧美成人| 国产精品电影一区二区三区| 日本a在线网址| 动漫黄色视频在线观看| 国产视频内射| 欧美激情 高清一区二区三区| 精品卡一卡二卡四卡免费| 久久亚洲真实| 97碰自拍视频| 日本三级黄在线观看| 久久久水蜜桃国产精品网| 欧美精品啪啪一区二区三区| 久久久国产成人免费| 巨乳人妻的诱惑在线观看| 精品一区二区三区av网在线观看| 午夜亚洲福利在线播放| 国产一区二区三区视频了| www国产在线视频色| 男女床上黄色一级片免费看| 国产成+人综合+亚洲专区| 亚洲五月婷婷丁香| 日韩成人在线观看一区二区三区| 人人澡人人妻人| 叶爱在线成人免费视频播放| 日本黄色视频三级网站网址| 国产单亲对白刺激| 波多野结衣高清作品| 国产成+人综合+亚洲专区| 国产精华一区二区三区| 91成人精品电影| 成人三级做爰电影| 国产蜜桃级精品一区二区三区| 制服人妻中文乱码| 国产精品一区二区三区四区久久 | 日韩成人在线观看一区二区三区| 精品日产1卡2卡| 国产av一区二区精品久久| 中文字幕最新亚洲高清| 欧美激情 高清一区二区三区| 国产在线精品亚洲第一网站| 99精品欧美一区二区三区四区| 丰满人妻熟妇乱又伦精品不卡| 精华霜和精华液先用哪个| 久久久久久亚洲精品国产蜜桃av| 亚洲自拍偷在线| 亚洲国产精品久久男人天堂| 国产欧美日韩精品亚洲av| 成人一区二区视频在线观看| 亚洲国产欧洲综合997久久, | 草草在线视频免费看| 精品午夜福利视频在线观看一区| 国内揄拍国产精品人妻在线 | 亚洲午夜精品一区,二区,三区| 亚洲黑人精品在线| 美女国产高潮福利片在线看| 最近最新免费中文字幕在线| 精品一区二区三区视频在线观看免费| 久99久视频精品免费| 女人高潮潮喷娇喘18禁视频| 亚洲男人的天堂狠狠| 中文字幕人妻丝袜一区二区| 一级毛片精品| 国产精品爽爽va在线观看网站 | 日韩成人在线观看一区二区三区| 国产视频内射| 中文字幕人成人乱码亚洲影| 日本三级黄在线观看| 在线免费观看的www视频| 国产99久久九九免费精品| 亚洲片人在线观看| 色尼玛亚洲综合影院| 欧美黑人欧美精品刺激| 看片在线看免费视频| 免费观看人在逋| 777久久人妻少妇嫩草av网站| 禁无遮挡网站| www.自偷自拍.com| 99久久综合精品五月天人人| 久久精品aⅴ一区二区三区四区| 亚洲人成网站高清观看| 视频区欧美日本亚洲| 日本熟妇午夜| 俺也久久电影网| 精品电影一区二区在线| www国产在线视频色| 亚洲激情在线av| 久久 成人 亚洲| 国产黄色小视频在线观看| 一边摸一边做爽爽视频免费| 国内精品久久久久久久电影| 久久精品夜夜夜夜夜久久蜜豆 | xxxwww97欧美| 国产精品香港三级国产av潘金莲| 国产激情偷乱视频一区二区| 午夜福利在线观看吧| 天天躁狠狠躁夜夜躁狠狠躁| 欧美最黄视频在线播放免费| 视频在线观看一区二区三区| 97人妻精品一区二区三区麻豆 | 国产精品一区二区三区四区久久 | 一夜夜www| 好男人电影高清在线观看| 观看免费一级毛片| 欧美在线一区亚洲| 人人妻人人澡欧美一区二区|