【摘要】：在核反應(yīng)堆的真實(shí)物理世界中,會(huì)天然存在一定的不確定性,包括幾何不確定性、初始條件和邊界條件不確定性、物理混沌狀態(tài)不確定性等,這類不確定性會(huì)使根據(jù)確定論假設(shè)模擬結(jié)果做出的決策存在失誤風(fēng)險(xiǎn)。目前在核工程領(lǐng)域,采用“最佳估算程序+不確定性分析”方法取代傳統(tǒng)的保守安全分析手段被世界各國(guó)、各組織所推崇,并隨之產(chǎn)生了一系列針對(duì)系統(tǒng)程序的不確定性分析方法。隨著計(jì)算機(jī)技術(shù)的不斷發(fā)展和核安全分析對(duì)高保真度物理模型和有效精確計(jì)算程序的要求,諸如CFD數(shù)值模擬的高精度仿真技術(shù)被越來(lái)越多的應(yīng)用于核電廠設(shè)計(jì)和安全分析領(lǐng)域中,而仿真精度提高的代價(jià)便是大量計(jì)算資源的消耗,這使得很多已經(jīng)在系統(tǒng)程序應(yīng)用上取得成功的以多次程序運(yùn)行結(jié)果為量化基礎(chǔ)的不確定性分析方法變得難以實(shí)施,因此,開(kāi)展基于高精度仿真程序的堆芯熱工水力計(jì)算不確定性分析方法的應(yīng)用研究具有重要意義。為了減少不確定性分析所需計(jì)算量和充分利用成熟的商用CFD計(jì)算軟件,本文以非嵌入式多項(xiàng)式混沌理論為基礎(chǔ),給出了適合應(yīng)用于堆芯CFD計(jì)算的輸入不確定性傳播分析方案,并以核反應(yīng)堆堆芯5×5光棒區(qū)冷卻劑通道CFD模型為研究對(duì)象,開(kāi)展了由邊界條件引起的流場(chǎng)輸出參數(shù)不確定度量化和敏感性分析,同時(shí)對(duì)相同的模型和輸入?yún)?shù)使用傳統(tǒng)的蒙特卡洛法進(jìn)行了不確定性分析作為量化結(jié)果對(duì)比,驗(yàn)證了多項(xiàng)式混沌法具有足夠的計(jì)算精度和極佳的計(jì)算經(jīng)濟(jì)性,可以應(yīng)用于需要耗費(fèi)龐大計(jì)算時(shí)間的堆芯熱工水力計(jì)算中。之后在模型中引入燃料棒與冷卻劑的換熱模型,開(kāi)展了考慮傳熱過(guò)程的不同工況下的光棒區(qū)冷卻劑通道不確定性分析研究。除了光棒區(qū)之外,堆芯內(nèi)還存在著一些復(fù)雜幾何結(jié)構(gòu)區(qū)域,而這些特殊的幾何結(jié)構(gòu)有可能會(huì)對(duì)流場(chǎng)的輸入不確定性造成影響,基于這種考慮,本文以含定位格架、攪混翼的堆芯復(fù)雜結(jié)構(gòu)區(qū)冷卻劑通道CFD模型為研究對(duì)象,開(kāi)展了基于多項(xiàng)式混沌法的不確定性分析研究,得到了復(fù)雜結(jié)構(gòu)區(qū)的輸出參數(shù)不確定性特點(diǎn),同時(shí)給出了敏感性分析結(jié)果。最后,本文對(duì)比了使用高斯分布和均勻分布作為輸入?yún)?shù)分布類型的不確定性分析結(jié)果,探討了不同輸入?yún)?shù)概率分布類型對(duì)多項(xiàng)式混沌法不確定性分析結(jié)果的影響。
【學(xué)位授予單位】：哈爾濱工程大學(xué)
【學(xué)位級(jí)別】：碩士
【學(xué)位授予年份】：2019
【分類號(hào)】：TM623
【圖文】：

圖 1.1 輸入?yún)?shù)不確定性傳播方法(1) 蒙特卡洛法(Monte-Carlo，MC)。蒙特卡洛法是最常用的確定參數(shù)樣本的方法，它的思想是在輸入?yún)?shù)的分布范圍內(nèi)對(duì)參數(shù)進(jìn)行簡(jiǎn)單隨機(jī)抽樣，通過(guò)對(duì)大樣本的重復(fù)計(jì)算，進(jìn)而得到輸出參數(shù)統(tǒng)計(jì)特性。該方法簡(jiǎn)便易行，只要抽樣規(guī)模足夠大，則最終計(jì)算出來(lái)的估計(jì)值就會(huì)趨于真實(shí)值，但該方法收斂速度很慢，需要大量的樣本才能取得可信的結(jié)果。法國(guó)的 CEA 將這一方法應(yīng)用于一個(gè)簡(jiǎn)單的“地板加熱”算例[16]，計(jì)算采用由 CEA 開(kāi)發(fā)的 CFD 計(jì)算程序 Trio_U，考慮存在對(duì)流和浮力效應(yīng)的 2 維穩(wěn)態(tài)流動(dòng)模型，在一開(kāi)始利用蒙特卡洛法抽取了 100 組輸入?yún)?shù)樣本，但是由于程序代碼漏洞，在實(shí)際運(yùn)行程序時(shí)，出現(xiàn)了很多次程序運(yùn)行失敗的情況，為了得到完整的統(tǒng)計(jì)信息，還需要重新抽樣補(bǔ)足樣本，最終運(yùn)行了三百多次程序才完成不確定度量化分析。從此例可以看出，抽樣方法需要的計(jì)算量非常大，這一缺點(diǎn)在需要占用大量 CPU 的 CFD 程序上顯得尤為突出。為了減少所需計(jì)算次數(shù)，數(shù)學(xué)領(lǐng)域有一些更為高效的抽樣方法可使樣本的標(biāo)準(zhǔn)差較[17][18]

圖 1.2 UMAE 方法分析流程上文對(duì) CFD 不確定性分析領(lǐng)域的現(xiàn)有工作進(jìn)行了回顧，一些學(xué)者正努力將已在系統(tǒng)程序應(yīng)用的輸入?yún)?shù)不確定性傳播方法和精度外推方法移植到 CFD 程序上，但不同于可以解決很多兩相問(wèn)題的系統(tǒng)程序模型，現(xiàn)在應(yīng)用于堆芯熱工水力安全分析的 CFD模型還是以單相為主，程序自身特性和物理過(guò)程建模的差異導(dǎo)致了相同的分析框架應(yīng)用到不同類型的程序上會(huì)出現(xiàn)很多需要改進(jìn)之處，從目前的幾個(gè)應(yīng)用研究中可以得到一些初步反饋：(1) CFD 程序不確定性來(lái)源多種多樣，包括：初始條件和邊界條件、材料物性、物理模型相關(guān)參數(shù)、無(wú)模型物理過(guò)程、數(shù)值解的誤差、幾何簡(jiǎn)化等；(2) 基于蒙特卡洛抽樣的不確定性傳播方法依然可以應(yīng)用于 CFD 領(lǐng)域，只是這種方法帶來(lái)的 CPU 代價(jià)過(guò)高；(3) 確定性抽樣法可以成為計(jì)算代價(jià)更低的一種選擇；(4) 在輸入?yún)?shù)個(gè)數(shù)較少時(shí)，多項(xiàng)式混沌法所需的計(jì)算次數(shù)也較少，當(dāng)使用一階多

圖 3.2 5×5 光棒區(qū)冷卻劑通道計(jì)算域幾何結(jié)構(gòu)與參數(shù)示意圖.1.2 網(wǎng)格劃分方案網(wǎng)格劃分在 CFD 分析中是人員參與度最高、所需工作時(shí)間最長(zhǎng)的一個(gè)研究環(huán)節(jié)格質(zhì)量的高低會(huì)對(duì)計(jì)算速度、收斂性和結(jié)果準(zhǔn)確性產(chǎn)生直接影響，而在后續(xù)的不確分析中，很難對(duì)網(wǎng)格優(yōu)劣導(dǎo)致的不確定度進(jìn)行量化，因此必須在初始建模階段尋找理的網(wǎng)格劃分方案，既能滿足計(jì)算時(shí)間成本可接受的條件，又能使仿真模型的計(jì)算達(dá)到要求，因此網(wǎng)格劃分是 CFD 不確定性分析中的重要基礎(chǔ)性工作之一。如圖 3.2 所示，光棒區(qū)冷卻劑通道在沿軸向的不同高度上有著相同的橫截面幾何，這種幾何特征適合選用棱柱形的網(wǎng)格進(jìn)行建模計(jì)算，棱柱形網(wǎng)格不僅網(wǎng)格質(zhì)量較且還可以在很大程度上節(jié)省建模所需網(wǎng)格量，這對(duì)節(jié)約計(jì)算資源有著重要意義。棱網(wǎng)格又可以分為三棱柱網(wǎng)格(Tetrahedral Prism，TP)、四棱柱網(wǎng)格(Hexahedral Prism 多棱柱網(wǎng)格(Polyhedral Prism，PP)，各類網(wǎng)格形狀如圖 3.3 所示，圖中 HP2[59]和 HP兩種應(yīng)用較為廣泛的結(jié)構(gòu)化 HP 網(wǎng)格。哈爾濱工程大學(xué)的陳廣亮[61]等人針對(duì)各類棱

【相似文獻(xiàn)】

相關(guān)期刊論文 前10條

1 郭文元 ,羅璋琳;核電系統(tǒng)故障概率不確定性分析的解析方法及其計(jì)算程序[J];核動(dòng)力工程;1988年02期

2 石林;;淺談Execl在不確定性分析中的應(yīng)用[J];金融經(jīng)濟(jì);2006年10期

3 侯普;房地產(chǎn)項(xiàng)目不確定性分析方法研究[J];建筑經(jīng)濟(jì);1995年06期

4 冉旭;吳丹;陳炳德;張渝;;最佳估算加不確定性分析方法及其應(yīng)用研究[J];核動(dòng)力工程;2013年03期

5 李錦;何平;;基于不確定性分析的高校科技創(chuàng)新與競(jìng)爭(zhēng)力研究[J];中國(guó)高�？萍寂c產(chǎn)業(yè)化(學(xué)術(shù)版);2006年S2期

6 王寶森,鄭丕諤,李秋英;在投資項(xiàng)目不確定性分析中盲數(shù)法與概率法的比較[J];天津大學(xué)學(xué)報(bào);2003年05期

7 吳雪春;工程經(jīng)濟(jì)評(píng)價(jià)的不確定性分析[J];安徽電力職工大學(xué)學(xué)報(bào);2002年03期

8 孟德彪;黃洪鐘;許煥衛(wèi);張小玲;張旭東;;一種多學(xué)科系統(tǒng)不確定性分析方法——協(xié)同不確定性分析法的改進(jìn)[J];機(jī)械工程學(xué)報(bào);2011年19期

9 王玲;;關(guān)于公路工程經(jīng)濟(jì)評(píng)價(jià)中不確定性分析與應(yīng)用探討[J];經(jīng)濟(jì)師;2018年07期

10 亢春萍;;淺析不確定性分析在水利工程項(xiàng)目經(jīng)濟(jì)評(píng)價(jià)中的應(yīng)用[J];財(cái)經(jīng)界(學(xué)術(shù)版);2013年16期

相關(guān)會(huì)議論文 前10條

1 靖劍平;張春明;孫微;賈斌;莊少欣;;最佳估算程序不確定性分析方法研究[A];第六屆反應(yīng)堆物理與核材料學(xué)術(shù)研討會(huì)第三屆核能軟件自主化研討會(huì)會(huì)議摘要集[C];2013年

2 王可;蔡杰進(jìn);任志豪;朱元兵;王婷;厲井鋼;;基于不確定性分析方法的程序驗(yàn)證研究[A];第十四屆全國(guó)反應(yīng)堆熱工流體學(xué)術(shù)會(huì)議暨中核核反應(yīng)堆熱工水力技術(shù)重點(diǎn)實(shí)驗(yàn)室2015年度學(xué)術(shù)年會(huì)論文集[C];2015年

3 王迪;王智峰;董勝;;海上風(fēng)機(jī)設(shè)計(jì)風(fēng)速的不確定性分析[A];第十九屆中國(guó)海洋（岸）工程學(xué)術(shù)討論會(huì)論文集（上）[C];2019年

4 鹿力成;馬力;;地聲反演不確定性分析[A];中國(guó)聲學(xué)學(xué)會(huì)2009年青年學(xué)術(shù)會(huì)議[CYCA’09]論文集[C];2009年

5 郭濤;;珠寶評(píng)估的不確定性分析[A];2013中國(guó)珠寶首飾學(xué)術(shù)交流會(huì)論文集[C];2013年

6 李爽;楊斌;齊飛;;基于人工神經(jīng)網(wǎng)絡(luò)的燃燒反應(yīng)動(dòng)力學(xué)模型的全局靈敏性和不確定性分析[A];中國(guó)化學(xué)會(huì)第29屆學(xué)術(shù)年會(huì)摘要集——第41分會(huì)：燃料與燃燒化學(xué)[C];2014年

7 林超;李廷秋;;黏性自由表面船舶繞流問(wèn)題的不確定性分析[A];第二十七屆全國(guó)水動(dòng)力學(xué)研討會(huì)文集（下冊(cè)）[C];2015年

8 石石;;基于概率法的油氣儲(chǔ)量不確定性分析[A];第八屆中國(guó)不確定系統(tǒng)年會(huì)論文集[C];2010年

9 殷建;宋松柏;;基于貝葉斯理論的P-Ⅲ型分布參數(shù)不確定性分析[A];流域水循環(huán)與水安全——第十一屆中國(guó)水論壇論文集[C];2013年

10 劉信恩;肖世富;莫軍;;用于不確定性分析的高斯過(guò)程響應(yīng)面模型設(shè)計(jì)點(diǎn)選擇方法研究[A];中國(guó)計(jì)算力學(xué)大會(huì)'2010（CCCM2010）暨第八屆南方計(jì)算力學(xué)學(xué)術(shù)會(huì)議（SCCM8）論文集[C];2010年

相關(guān)博士學(xué)位論文 前10條

1 馬遠(yuǎn)卓;不確定性分析和結(jié)構(gòu)優(yōu)化新方法及其在熱防護(hù)系統(tǒng)的應(yīng)用[D];南京航空航天大學(xué);2018年

2 許鑫;矩獨(dú)立和鄰域重要性理論及應(yīng)用研究[D];西北工業(yè)大學(xué);2015年

3 夏慧;小型無(wú)人直升機(jī)的模型辨識(shí)與魯棒控制技術(shù)研究[D];南京理工大學(xué);2015年

4 張江江;地下水污染源解析的貝葉斯監(jiān)測(cè)設(shè)計(jì)與參數(shù)反演方法[D];浙江大學(xué);2017年

5 易俐娜;面向?qū)ο筮b感影像分類不確定性分析[D];武漢大學(xué);2011年

6 龔偉;基于信息熵和互信息的流域水文模型不確定性分析[D];清華大學(xué);2012年

7 譚君位;作物模型參數(shù)敏感性和不確定性分析方法研究[D];武漢大學(xué);2017年

8 劉廣怡;基于貝葉斯理論的通信網(wǎng)絡(luò)數(shù)據(jù)處理技術(shù)研究[D];解放軍信息工程大學(xué);2015年

9 龍剛;深井復(fù)雜地層井筒套管強(qiáng)度可靠性評(píng)價(jià)[D];中國(guó)石油大學(xué)(華東);2013年

10 劉冀;徑流分類組合預(yù)報(bào)方法及其應(yīng)用研究[D];大連理工大學(xué);2008年

相關(guān)碩士學(xué)位論文 前10條

1 葉潛;小型堆小破口失水事故最佳估算及不確定性分析[D];華中科技大學(xué);2019年

2 張梟羽;多項(xiàng)式混沌法在堆芯計(jì)算不確定性分析中的應(yīng)用研究[D];哈爾濱工程大學(xué);2019年

3 石濤;多種工況下的中子活化計(jì)算不確定性分析[D];華北電力大學(xué)(北京);2019年

4 阮曉行;考慮認(rèn)知參數(shù)或多峰隨機(jī)參數(shù)的不確定性分析方法[D];湖南大學(xué);2018年

5 徐佳意;考慮不同核素干涉效應(yīng)的截面不確定性分析[D];華北電力大學(xué)(北京);2018年

6 黨英杰;運(yùn)移轉(zhuǎn)換模型的不確定性分析方法及其應(yīng)用[D];長(zhǎng)安大學(xué);2014年

7 張惠雁;基于微擾理論的計(jì)算不確定性分析方法研究[D];哈爾濱工程大學(xué);2016年

8 鄒濤;基于整體敏感性分析與整體不確定性分析的海洋及海岸結(jié)構(gòu)風(fēng)險(xiǎn)評(píng)估新方法[D];中國(guó)海洋大學(xué);2012年

9 李霞;印尼R區(qū)塊復(fù)雜斷塊油藏描述及建模不確定性分析[D];長(zhǎng)安大學(xué);2014年

10 石磊;堆芯物理計(jì)算不確定性分析方法研究[D];哈爾濱工程大學(xué);2014年

本文編號(hào)：2801723