![核工程與核技術(shù)畢業(yè)設(shè)計(jì)(論文)閃爍探測器探測效率與γ射線能量關(guān)系的數(shù)值模擬_第1頁](http://file2.renrendoc.com/fileroot_temp3/2021-10/4/f6d6db73-ef5e-4d5a-9983-de33b9ae6e60/f6d6db73-ef5e-4d5a-9983-de33b9ae6e601.gif)
![核工程與核技術(shù)畢業(yè)設(shè)計(jì)(論文)閃爍探測器探測效率與γ射線能量關(guān)系的數(shù)值模擬_第2頁](http://file2.renrendoc.com/fileroot_temp3/2021-10/4/f6d6db73-ef5e-4d5a-9983-de33b9ae6e60/f6d6db73-ef5e-4d5a-9983-de33b9ae6e602.gif)
![核工程與核技術(shù)畢業(yè)設(shè)計(jì)(論文)閃爍探測器探測效率與γ射線能量關(guān)系的數(shù)值模擬_第3頁](http://file2.renrendoc.com/fileroot_temp3/2021-10/4/f6d6db73-ef5e-4d5a-9983-de33b9ae6e60/f6d6db73-ef5e-4d5a-9983-de33b9ae6e603.gif)
![核工程與核技術(shù)畢業(yè)設(shè)計(jì)(論文)閃爍探測器探測效率與γ射線能量關(guān)系的數(shù)值模擬_第4頁](http://file2.renrendoc.com/fileroot_temp3/2021-10/4/f6d6db73-ef5e-4d5a-9983-de33b9ae6e60/f6d6db73-ef5e-4d5a-9983-de33b9ae6e604.gif)
![核工程與核技術(shù)畢業(yè)設(shè)計(jì)(論文)閃爍探測器探測效率與γ射線能量關(guān)系的數(shù)值模擬_第5頁](http://file2.renrendoc.com/fileroot_temp3/2021-10/4/f6d6db73-ef5e-4d5a-9983-de33b9ae6e60/f6d6db73-ef5e-4d5a-9983-de33b9ae6e605.gif)
版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡介
1、 畢業(yè)設(shè)計(jì)(論文) 題 目 閃爍探測器探測效率與射線能量關(guān)系的數(shù)值模擬英文題the numerical simulation relationship of scintillation detector detection efficiency and -ray energy學(xué)生姓名: 專 業(yè): 核工程與核技術(shù)班 級(jí): 指導(dǎo)教師: 二零一一年六月摘 要閃爍能譜儀的探測效率與射線能量關(guān)系之間的關(guān)系,通過使用蒙特卡羅的數(shù)值模擬方法進(jìn)行研究。nai閃爍探測器由于其探測效率高、在常溫下可以使用的特點(diǎn),越來越受到人們的重視。因此研究其探測效率對(duì)進(jìn)一步的開發(fā)和研究閃爍探測器是很有必要的,特別是閃爍探測器與
2、能量之間的關(guān)系決定了今后的研究方向。射線主要分為兩大類:一類是,自然射線;另一類是,人工射線,用nai探測器了對(duì)射線進(jìn)行探測。然后用mcnp4c軟件進(jìn)行模擬實(shí)驗(yàn),得到了不同射線能量下的能譜。通過能譜圖我們可以根據(jù)全能峰內(nèi)的相對(duì)計(jì)數(shù)率,得出相對(duì)探測效率。再得出相對(duì)探測效率以后,根據(jù)相對(duì)應(yīng)的射線能量,擬合出圖形,得到探測效率與射線能量關(guān)系的數(shù)值模擬。通過對(duì)上述內(nèi)容的研究,本課題基本上達(dá)到了預(yù)期研究目的,取得較好的階段性成果,對(duì)閃爍探測器探測效率與能量關(guān)系的研究起到一定的科學(xué)指導(dǎo)作用,為實(shí)驗(yàn)提供了一定的依據(jù)。關(guān)鍵詞:探測效率; 射線; 蒙特卡羅模擬 abstractthe relationship
3、between detection efficiency of the -ray scintillation spectrometer and the -ray energy, we are using monte carlo simulation method to study. because of its high detection efficiency, can be used at room temperature characteristics, nai scintillation detector more attention has been paid. 朗讀顯示對(duì)應(yīng)的拉丁字
4、符的拼音therefore, it is very necessary to study the detection efficiency for further development and research, especially the relationship between scintillation detectors and -ray energy direction of our future research. -ray can be divided into two categories: one is the natural -ray; the other is art
5、ificial -ray, with the nai detector for -ray detection. then simulation experiments by mcnp4c software. obtained under different -ray energy spectrum of -ray. we obtained the relative detection efficiency by the energy spectra we can according to all-powerful within the relative peak count rate. and
6、 then come after the relative detection efficiency, according to the corresponding -ray energy, fitting out graphics, and -ray detection efficiency has been the relationship between energy simulation. through the above research, this research project is basically to achieve the desired objectives an
7、d achieve good initial results, the efficiency of the scintillation detector and -energy relationship of scientific research play a guiding role for the experiment provided some basis.key words: detection efficiency; -rays; monte carlo simulation 目 錄緒 論11.1 選題意義及目的11.2 實(shí)驗(yàn)背景12 模擬實(shí)驗(yàn)基礎(chǔ)32.1 射線與物質(zhì)的相互作用32
8、.1.1 光電效應(yīng)32.1.2 康普頓效應(yīng)32.1.3 光電子對(duì)效應(yīng)42.2 閃爍能譜儀52.2.1 閃爍探測器的原理62.2.2 閃爍探測效率73 構(gòu)建模型83.1 蒙特卡羅方法及mcnp程序簡介83.1.1 蒙特卡羅計(jì)算方法及其應(yīng)用范圍83.1.2 mcnp程序簡介83.2 建立數(shù)學(xué)模型113.2.1 點(diǎn)源模擬113.2.2 探測器柵元卡片填寫113.2.3 模擬模型143.2.4 源定義143.2.5 輸入卡片154 計(jì)算結(jié)果與處理174.1 模擬輸出結(jié)果174.1.1 自然射線的模擬輸出結(jié)果174.1.2 人工射線的模擬輸出結(jié)果204.2 數(shù)據(jù)處理與結(jié)果23結(jié)論26致 謝27參考文獻(xiàn)2
9、830東華理工大學(xué)畢業(yè)設(shè)計(jì)(論文) 緒論緒 論1.1 選題意義及目的隨著我國經(jīng)濟(jì)的發(fā)展,國民生活水平不斷提高,國內(nèi)需求不斷增長。核輻射探測器也逐漸被使用到各行各業(yè),在各類探測器中,閃爍體探測器由于其對(duì)射線的探測分辨時(shí)間短、探測效率高、能測量射線的能量的優(yōu)點(diǎn),所以是目前應(yīng)用的最廣的射線探測器。因此研究閃爍探測器探測效率與射線能量關(guān)系對(duì)我們提高閃爍探測器探測效率有著重要的意義。本課題主要是應(yīng)用蒙特卡羅的方法對(duì)探測器探測效率與射線能量關(guān)系進(jìn)行數(shù)值模擬,通過這個(gè)數(shù)值關(guān)系,我們可以找到最佳的射線能量對(duì)應(yīng)著最高的探測效率。這樣就對(duì)我們以后的實(shí)驗(yàn)、研究提供了一定的依據(jù)。1.2 實(shí)驗(yàn)背景蒙特卡羅方法又名隨機(jī)模
10、擬法或統(tǒng)計(jì)實(shí)驗(yàn)法,它是一種具有獨(dú)特風(fēng)格的數(shù)值計(jì)算方法,它既能求解定性的數(shù)學(xué)問題,也能求解隨即問題。半個(gè)多世紀(jì)以來,由于科學(xué)技術(shù)的發(fā)展和電子計(jì)算機(jī)的發(fā)明,這種方法作為一種獨(dú)立的方法被提出來,但與一般數(shù)值計(jì)算方法有很大區(qū)別,它是以概率統(tǒng)計(jì)理論為基礎(chǔ)的一種方法。由于蒙特卡羅方法能夠比較逼真地描述事物的特點(diǎn)及物理實(shí)驗(yàn)過程,解決一些數(shù)值方法難以解決的問題,因而該方法的應(yīng)用領(lǐng)域日趨廣泛1。蒙特卡洛方法由于其簡單性、靈活性和普遍性在獲得廣泛應(yīng)用。尤其在核物理上更是得到很大應(yīng)用,20世紀(jì)80年代,schlumberger、lanl和美國北卡羅來納州立大學(xué)開始將蒙特卡羅方法用于核測井的模擬,取得很大的成功2。以
11、后的發(fā)展里,隨著計(jì)算機(jī)技術(shù)的不斷進(jìn)步,現(xiàn)代計(jì)算機(jī)技術(shù)使隨即采樣成為解決復(fù)雜問題的實(shí)用方法。到了20世紀(jì)90年代,蒙特卡羅模擬已經(jīng)成為核測井科研和生產(chǎn)中的重要工具。近一些年來,一些大型的蒙特卡羅方法應(yīng)用軟件先后生產(chǎn)并發(fā)展,例如,美國的橡樹嶺國家實(shí)驗(yàn)室的輻射屏蔽信息中心rsic、阿貢國家實(shí)驗(yàn)的軟件中心nesc和計(jì)算機(jī)物理庫cpc等所收集的大量應(yīng)用軟件,就有相當(dāng)數(shù)量的蒙特卡洛方法應(yīng)用軟件。具體的用在核技術(shù)領(lǐng)域的軟件有美國橡樹嶺國家實(shí)驗(yàn)室的mose程序、美國stanford linear accelerator center的egs程序、美國生地國家實(shí)驗(yàn)室的sandyl和tiger程序、美國洛斯阿拉莫
12、斯國家實(shí)驗(yàn)室的mcnp程序、歐洲核子研究中心的gernt程序、英國的mcbend程序以及法國的tripoli程序等。這些程序在粒子輸運(yùn)模擬中發(fā)揮了重要作用,但這些程序一般為通用程序,并且操作復(fù)雜不便計(jì)算和探測問題。因而,研究核探測問題的蒙特卡羅計(jì)算方法及軟件成為當(dāng)今輻射輸運(yùn)領(lǐng)域的熱門課題2,3。在運(yùn)用蒙特卡羅方法對(duì)射線的響應(yīng)函數(shù)進(jìn)行計(jì)算以及探測效率刻度方面,國內(nèi)外重20世紀(jì)60年代就已經(jīng)有很多人開始研究工作,而從20世紀(jì)70年代以來,國內(nèi)外有關(guān)這方面的研究工作可分為:對(duì)裸晶體,隨著光子能量越來越高,考慮的反應(yīng)機(jī)制日趨完善;對(duì)低能光子,在反應(yīng)機(jī)制較簡單的情況下,研究了體源、反射層及源與探測器的幾
13、何關(guān)系對(duì)響應(yīng)函數(shù)及探測效率的影響2,3。進(jìn)入20世紀(jì)90年代以后,國內(nèi)外持續(xù)有學(xué)者對(duì)蒙特卡洛方法計(jì)算響應(yīng)函數(shù)和探測效率的方法進(jìn)行改進(jìn)完善,像h.vnguyen等利用ngrc程序、fayezh.h等利用egs4程序模擬了nai(tl)閃爍探測器對(duì)射線的響應(yīng)函數(shù),清華大學(xué)的hu-xia-shi等人利用mcnp模擬了nai(tl) 閃爍探測器對(duì)射線的響應(yīng)函數(shù)3,4??v觀前人的研究成果,可以發(fā)現(xiàn)蒙特卡洛方法在計(jì)算射線在探測器中與探測器效率刻度方面,為我們研究提供了極大的便利,并可以減少大量的實(shí)驗(yàn)工作,但是在大多停留在實(shí)用國外的通用軟件上,而國內(nèi)還沒有比較成功的自主研發(fā)陳鼓工我們方便使用,并且在探測效率
14、刻度大多是正對(duì)源與探測器的空間幾何位置來進(jìn)行,很少有人對(duì)射線能量與探測效率之間的關(guān)系進(jìn)行研究3,4。東華理工大學(xué)畢業(yè)設(shè)計(jì)(論文) 模擬實(shí)驗(yàn)基礎(chǔ)2 模擬實(shí)驗(yàn)基礎(chǔ)2.1 射線與物質(zhì)的相互作用2.1.1 光電效應(yīng)當(dāng)射線與原子殼層軌道上的內(nèi)層電子碰撞時(shí),將所有能量交給殼層的一個(gè)內(nèi)層電子,電子克服了電離能,脫離原子而運(yùn)動(dòng),產(chǎn)生光電子,而射線被完全吸收,這種作用稱為光電效應(yīng)。為使光電效應(yīng)保持動(dòng)量守恒,必須有第三者,即反沖核。由于自由電子不能吸收射線而形成光電子,因此光電效應(yīng)產(chǎn)生的幾率將隨著電子結(jié)合能的增加而增加,靠近k層、l層產(chǎn)生光電效應(yīng)的幾率最大1,2。圖21 光電效應(yīng)示意圖當(dāng)入射射線的能量大于原子的電
15、離能,才能產(chǎn)生該物質(zhì)的光電效應(yīng)。光電子的動(dòng)能是入射射線能量與電離能之差。一般認(rèn)為,光電效應(yīng)視整個(gè)原子看成一個(gè)統(tǒng)一的體系與入射光子發(fā)生作用,吸收能量,發(fā)射光電子。沒有足夠的能量克服軌道電子的結(jié)合能,電子就不可能脫離原子的束縛而成為自由電子1,2。2.1.2 康普頓效應(yīng) 當(dāng)入射的射線與一個(gè)自由電子發(fā)生一次碰撞時(shí),射線將部分能量傳給電子,使它與射線的初始運(yùn)動(dòng)方向呈角射出,而射線則與初始運(yùn)動(dòng)方向呈角散射。這種現(xiàn)象稱為康普頓吳有訓(xùn)效應(yīng),簡稱康普頓效應(yīng)1,2。將入射光子看成是具有一定能量和動(dòng)量的粒子與電子發(fā)生剛性碰撞。在這一碰撞過程中,入射光子將一部分能量傳遞給電子。而光子發(fā)生能量和方向的改變,形成能量為
16、,動(dòng)量為的光量子。散射角為。電子在獲得一定能量之后,克服結(jié)合能,沿方向出射。稱為反沖電子或康普頓電子。其動(dòng)能取決與電子獲得能量與電子結(jié)合能之差;并具有相應(yīng)的動(dòng)量,作用過程見圖22所示。根據(jù)能量和動(dòng)量守恒定律,可以得到沿方向散射的射線能量: (21)沿方向出射的反沖電子能量: (22) 式中:電子靜止質(zhì)量; 光速。圖22 康普頓效應(yīng)示意圖2.1.3 光電子對(duì)效應(yīng)隨著入射光子能量的增高,光電效應(yīng)的吸收作用很寬減弱,康普頓效應(yīng)也逐漸減弱。當(dāng)光子能量大于1.02mev時(shí),就存在形成電子對(duì)效應(yīng)的幾率,即光子完全被吸收而產(chǎn)生一正、負(fù)電子對(duì)。電子對(duì)的動(dòng)能為 (2-3)式中,入射光仔能量; 正、負(fù)電子的靜止質(zhì)
17、量; 光速。正電子與負(fù)電子不同,它很不穩(wěn)定,其壽命大約為。正電子在物理中由于電力而損失能量后,將和電子結(jié)合,通常轉(zhuǎn)化為兩個(gè)能來那個(gè)為0.51mev且方向相反的光子,這種現(xiàn)象即正電子湮沒1,2。當(dāng)能量大于1.02mev的光子通過原子序數(shù)較大的物質(zhì)時(shí),將容易測到這種射線。圖23 電子對(duì)效應(yīng)示意圖2.2 閃爍能譜儀閃爍體是透明的絕緣體。射線與閃爍體相互作用所產(chǎn)生的次級(jí)電子可以再和它作用,使其原子分子電離和激發(fā)。顯然,產(chǎn)生的電子離子是無法從物質(zhì)中引出來的,但其過程中產(chǎn)生的閃光可從物質(zhì)中透射出來,被收集轉(zhuǎn)換成點(diǎn)脈沖放大后,可供分析記錄。因此,閃爍體可作射線探測器的探測介質(zhì)。閃爍體可分為有機(jī)和無機(jī)兩大類,
18、有機(jī)閃爍體又分為液體閃爍體、有機(jī)晶體和塑料閃爍體三種。這些閃爍體對(duì)射線有中等大小的吸收本領(lǐng)。有機(jī)閃爍體的探測器有比較高的探測效率。盡管其能量分辨率不好,但響應(yīng)快、時(shí)間分辨好。因此,常用這類探測器用作快符合裝置上的探測器。無機(jī)閃爍體的典型代表是鉈激發(fā)的碘化鈉晶體和碘化銫晶體。他們原子序數(shù)高、密度大,對(duì)射線有較高的吸收本領(lǐng)。其時(shí)間響應(yīng)快,能量分辨率也好,是20世紀(jì)50年代以來探測射線的主要探測器2,3。2.2.1 閃爍探測器的原理閃爍探測器主要由閃爍體、光電倍增管和相應(yīng)的電子儀器贊歌主要部分組成。圖2-4是閃爍能譜的結(jié)構(gòu)框圖,探頭部分主要包括閃爍體和光電倍增管和前置放大器;主機(jī)部分包括主放大器,多
19、道脈沖幅度分析器,接口電路、電源電路以及計(jì)算機(jī)和相應(yīng)的應(yīng)用軟件等。圖2-4 閃爍能譜儀原理框圖目前,閃爍晶體主要使用nai(tl)晶體,它的密度大(=3.67g/cm3),而且高原子系數(shù)的碘(z=53)占重要的85%,所以對(duì)射線探測效率特別高,相對(duì)發(fā)光效率大,它的光譜最強(qiáng)波長為415nm左右,能與光電倍增管的光譜響應(yīng)較好的一種。但是nai(tl)晶體容易潮解,吸收空氣水分會(huì)變質(zhì)失效,所以一般nai(tl)晶體使用時(shí)都是裝在密封的金屬盒中。射線與nai(tl)晶體相互作用,會(huì)使其電離、激發(fā)而發(fā)射熒光,從閃爍體出來的光子通過光導(dǎo)射向光電倍增管的光陰極,由于光電效應(yīng),在光陰極上打出電子。光電子經(jīng)過各
20、級(jí)之間的電壓加速,聚焦后射向打拿極。每個(gè)光子在打拿極上擊出幾個(gè)電子,這些電子射向后續(xù)的打拿極。所以,最后射向陽極的電子數(shù)目是很多的,陰極把這些電子收集起來,轉(zhuǎn)變成電信號(hào)并輸出。經(jīng)過nai(tl)探測器和光電倍增管輸出的信號(hào)往往較弱,需要加前置放大器。其主要作用是:提高系統(tǒng)的信噪比;減小信號(hào)經(jīng)電纜傳送時(shí)外界干擾的影響;阻抗匹配,功率放大。探頭的信號(hào)經(jīng)前置放大器輸出,但前置放大器輸出脈沖幅度和波形并不適合后面的分析測量設(shè)備的要求,所以對(duì)信號(hào)還需要進(jìn)一步放大和成形,在放大和成形的過程中必須嚴(yán)格保持探測器輸出的有用信號(hào)(如射線的能量和時(shí)間信息),可能減少它們的失真。這樣的放大和成形就由放大器來完成。放
21、大器輸出的信號(hào)要適應(yīng)分析測量設(shè)備的要求,必須解決兩個(gè)問題。一個(gè)是把小信號(hào)放大到需要的幅度;另一個(gè)就是改造信號(hào)形狀,通常稱為濾波成形。目的是放大有用的信號(hào),降低造成,提高信噪比,適合于后續(xù)電路的測量。在這個(gè)過程中盡可能不損失有用的信息。主放大器輸出的信號(hào)就可以給多道脈沖幅度分析器(mca)進(jìn)行脈沖的幅度分析。多道譜儀是按照脈沖的幅度對(duì)輻射脈沖信號(hào)進(jìn)行分類計(jì)數(shù),以及對(duì)后續(xù)輻射信息保存、處理以及獲取有用信息的儀器。它的核心部件就是多道脈沖幅度分析器,其基本功能就是按輸入脈沖的幅度分類計(jì)數(shù)。多道脈沖幅度分析器的任務(wù)是將被測量的脈沖幅度范圍平均分成2n個(gè)幅度間隔,然后測量脈沖幅度在每一個(gè)“幅度間隔”內(nèi)的
22、輸入脈沖個(gè)數(shù),最后得到輸入信號(hào)的脈沖幅度分布曲線。其測量方法是用計(jì)算機(jī)技術(shù)中的a/d轉(zhuǎn)換以及數(shù)據(jù)存儲(chǔ)技術(shù)。它將能夠分析的脈沖幅度范圍分成多個(gè)幅度間隔,這些幅度間隔的個(gè)數(shù)就是脈沖幅度分析器的道數(shù);而幅度間隔的寬度就是道寬。道數(shù)越多,幅度分布分析得越精細(xì),各個(gè)道的計(jì)數(shù)也相應(yīng)減少,則需要測量的時(shí)間也要加長,當(dāng)然硬件電路也就越復(fù)雜,所以不應(yīng)盲目地追求道數(shù)。通常,要求在幅度峰的半寬度范圍內(nèi)應(yīng)有510道,對(duì)于采用nai(tl)探測器的多道能譜儀,由于它的能量分辨率比較差,128道至512道就能滿足測量要求2,3。2.2.2 閃爍探測效率要知道探測效率與射線能量之間的數(shù)值關(guān)系,我們首先要探測效率的定義和如何
23、計(jì)算探測效率。目前,作為探測效率的定義是多樣的,例如集合效率、本征效率、峰-總比效率、本征全能效率、絕對(duì)全能峰效率、絕對(duì)雙逃逸峰效率和絕對(duì)單逃逸峰效率等3。在本文中,我們是用絕對(duì)全能峰效率來進(jìn)行對(duì)探測器探測效率的計(jì)算。如果只測量全能峰內(nèi)的計(jì)數(shù)率,則探測裝置的噪聲和干擾輻射對(duì)計(jì)數(shù)率的影響都可減少。全能峰很容易辨認(rèn),從測得的能譜中求得全能峰內(nèi)的凈計(jì)數(shù)率都是比較容易的。根據(jù)全能峰內(nèi)的凈計(jì)數(shù)率與全能峰探測效率,求樣品射線發(fā)射率的測量方法稱為全能峰法。利用全能峰法我們就可以測出全能峰探測效率。例如,在一定條件,測量發(fā)射率為n0的標(biāo)準(zhǔn)源,得到全能峰內(nèi)的凈計(jì)數(shù)率為np則全能峰探測效率為5: (2-1) 東華
24、理工大學(xué)畢業(yè)設(shè)計(jì)(論文) 構(gòu)建模型3 構(gòu)建模型3.1 蒙特卡羅方法及mcnp程序簡介蒙特卡羅方法又稱隨機(jī)抽樣技巧或統(tǒng)計(jì)試驗(yàn)方法。半個(gè)多世紀(jì)以來,由于科學(xué)技術(shù)的發(fā)展和電子計(jì)算機(jī)的出現(xiàn)與發(fā)展,這種方法作為一種獨(dú)立的方法被提出來,并首先在核武器的試驗(yàn)與研制中得到應(yīng)用。3.1.1 蒙特卡羅計(jì)算方法及其應(yīng)用范圍蒙特卡羅方法是一種計(jì)算方法,但與一般數(shù)值計(jì)算方法有很大區(qū)別。它是以概率統(tǒng)計(jì)理論為基礎(chǔ)的一種方法。由于蒙特卡羅方法能夠比較逼真地描述事物的特點(diǎn)及物理實(shí)驗(yàn)過程,解決一些數(shù)值方法難以解決的問題,因而該方法的應(yīng)用領(lǐng)域日趨廣泛4。蒙特卡羅方法的基本思想,當(dāng)所求問題的解是某個(gè)事件的概率,或者是某個(gè)隨機(jī)變量的數(shù)
25、學(xué)期望,或者是與概率、數(shù)學(xué)期望有關(guān)的量時(shí),通過某種試驗(yàn)的方法,得出該事件發(fā)生的頻率,或者該隨機(jī)變量若干個(gè)具體觀察值的算術(shù)平均值,通過它得到問題的解4,5。蒙特卡羅方法的優(yōu)點(diǎn),能夠比較逼真地描述具有隨機(jī)性質(zhì)的事物的特點(diǎn)及物理實(shí)驗(yàn)過程;受幾何條件限制?。皇諗克俣扰c問題的維數(shù)無關(guān);具有同時(shí)計(jì)算多個(gè)方案與多個(gè)未知量的能力;誤差容易確定;程序結(jié)構(gòu)簡單,易于實(shí)現(xiàn)。同時(shí)它的缺點(diǎn),收斂速度慢;誤差具有概率性;在粒子輸運(yùn)問題中,計(jì)算結(jié)果與系統(tǒng)大小有關(guān)。故使用蒙特卡羅方法需“揚(yáng)長避短”,發(fā)揮蒙特卡羅的特長,使其應(yīng)用范圍更加廣泛4,5,6。蒙特卡羅方法應(yīng)用范圍包括:粒子輸運(yùn)問題,統(tǒng)計(jì)物理,典型數(shù)學(xué)問題,真空技術(shù),激
26、光技術(shù)以及醫(yī)學(xué),生物,探礦等方面。其在粒子輸運(yùn)問題中的應(yīng)用范圍主要包括:實(shí)驗(yàn)核物理,反應(yīng)堆物理,高能物理等方面。在實(shí)驗(yàn)核物理中的應(yīng)用范圍主要包括:通量及反應(yīng)率,中子探測效率,光子探測效率,光子能量沉積譜及響應(yīng)函數(shù),氣體正比計(jì)數(shù)管反沖質(zhì)子譜,多次散射與通量衰減修正等方面6。3.1.2 mcnp程序簡介mcnp程序全名為monte carlo neutron and photo transport code。它是美國los alamos國家實(shí)驗(yàn)室開發(fā)的大型多功能通用蒙特卡羅程序,可以計(jì)算中子、光子和電子的聯(lián)合輸運(yùn)問題以及臨界問題,中子能量范圍從10-11mev至20mev,光子和電子的能量范圍從1
27、kev至1000mev。程序采用獨(dú)特的曲面組合幾何結(jié)構(gòu),使用點(diǎn)截面數(shù)據(jù),程序通用性較強(qiáng),與其它程序相比,mcnp程序中的減方差技巧是比較多而全的4,5,6。mcnp程序主要特點(diǎn),程序中的幾何是三維任意組態(tài);可以使用精細(xì)的點(diǎn)截面數(shù)據(jù);該程序功能齊全;在減小方差技巧方面,內(nèi)容十分豐富;具有很強(qiáng)的通用性7。mcnp的輸入包括幾個(gè)文件,但主要的一個(gè)是由用戶編寫的inp文件,該文件包括描述問題所必須的全部輸入信息。文件采用卡片結(jié)構(gòu),每行代表一張卡片,文件由一系列卡片組成,對(duì)于任一特定的問題,只需用到inp全部輸入卡片的一小部分。例如對(duì)不同的源,只須修改源的描述卡;不同的計(jì)數(shù)類型,有對(duì)應(yīng)的計(jì)數(shù)描述卡8。m
28、cnp輸入數(shù)據(jù)中,物理量的單位見表3-1:表3-1 mcnp輸入文件中物理量的單位長度厘米能量mev時(shí)間10-8秒溫度mev(kt)原子密度1024個(gè)原子/厘米3截面10-24厘米2原子量中子質(zhì)量的1.008664967質(zhì)量密度克/厘米3阿伏伽德羅常數(shù)6.0231023mcnp程序輸入文件的基本形式,使用這種形式建立一個(gè)蒙特卡羅計(jì)算問題,對(duì)問題的幾何結(jié)構(gòu)、材料、記數(shù)要求等等給以描述,如果需要,可以直接運(yùn)行。該文件的格式如下:表3-2 mcnp初始運(yùn)行的輸入文件信息塊空行分隔符選擇項(xiàng)標(biāo)題卡僅一行,占用地1-80列。作為輸入標(biāo)題。柵元卡空行分隔符定義構(gòu)成整個(gè)系統(tǒng)的各個(gè)基本戒指單元以及相應(yīng)的物理信息
29、。曲面卡空行分隔符定義組成柵元的曲面信息。數(shù)據(jù)卡空行分隔符其它數(shù)據(jù),包括問題類型、源描述、材料描述、計(jì)數(shù)描述,問題截?cái)鄺l件等。其它選擇項(xiàng)信息塊空行分隔符選擇項(xiàng)continue寫在1-8列數(shù)據(jù)卡空行分隔符只允許部分?jǐn)?shù)據(jù)卡。(fq,dd,nps,ctme,idum,rdum,prdmp,lost,dbcn,print,kcode,mplot,za,zb和zc)其它選擇項(xiàng)卡片填寫格式注意問題,inp輸入文件的每一行(稱之為一張卡片)都限于使用第180列并構(gòu)成卡片映象。大部分輸入卡片按行填寫;然而,對(duì)數(shù)據(jù)卡允許按列填寫。$符號(hào)為它所在那行數(shù)據(jù)的結(jié)束符,在$符號(hào)后面的內(nèi)容作為注釋,它可從$符號(hào)后面的任一
30、列開始。標(biāo)題卡只占一行,整行都可填入用戶需要的信息,也可以是空行。但要注意在其它地方使用空行是作為結(jié)束符或者分隔符。輸入文件中,在標(biāo)題卡之后及最后的空行結(jié)束卡之前的任何地方都可插入注釋卡。注釋卡必須是字母“c”寫在15列中的任意位置,且至少用一個(gè)空格隔開后面的注釋內(nèi)容9。其中行輸入格式,柵元卡、曲面卡和數(shù)據(jù)卡的書寫格式是相同的。必須從15列開始填寫這些卡片相應(yīng)的名字(或編號(hào))和粒子標(biāo)識(shí)符,后面填寫用空格分隔的數(shù)據(jù)項(xiàng)。如果15列為空,則表示它是前一張卡片的繼續(xù)卡。如果在一行的末尾有一個(gè)用空格隔開的符號(hào)“&”,則表示下一行是該行的繼續(xù)卡,數(shù)據(jù)可填寫在180列。一個(gè)數(shù)據(jù)項(xiàng)必須在一張卡片上寫完,不得跨
31、到下一張卡片上。完全空白的一行則為兩組卡片的分隔符。對(duì)任何給定的帶有粒子標(biāo)識(shí)符的類型卡只能有一張。需要整數(shù)的數(shù)據(jù)項(xiàng)必須填寫整數(shù),其它數(shù)據(jù)可填寫為整數(shù)或浮點(diǎn)數(shù)以及mcnp能讀的數(shù)據(jù)10。列輸入格式對(duì)輸入柵元參數(shù)及源描述是非常有用的。對(duì)于柵元重要性及體積等參數(shù),輸入按行排列時(shí)其可讀性差,且在增加或刪除一些柵元時(shí)容易出錯(cuò)。用列輸入格式,一個(gè)柵元的所有柵元參數(shù)是放在標(biāo)有該柵元名字的那行上。如果刪掉一個(gè)柵元,用戶只需刪除該柵元參數(shù)行,而不需要在每一個(gè)柵元參數(shù)卡上尋找該柵元所對(duì)應(yīng)的數(shù)據(jù)項(xiàng)。對(duì)于源描述也有類似的情況。用列格式,卡片名字逐個(gè)放在一個(gè)輸入行上,并且在這些卡片名字下面按列列出數(shù)據(jù)項(xiàng)。后續(xù)各行為各個(gè)
32、柵元的數(shù)據(jù)。如果填寫某個(gè)柵元名字,則必須填寫全部柵元名字,且柵元的順序可以任意排列;如果沒有指定柵元名字,則按柵元卡描述的順序排列。在一個(gè)輸入文件中,允許有多個(gè)列數(shù)據(jù)塊11。3.2 建立數(shù)學(xué)模型3.2.1 點(diǎn)源模擬本實(shí)驗(yàn),模擬探測器使用nai(tl)多道能譜儀(圖3-1)替代閃爍探測器。為了得到探測器探測效率和射線能量之間的數(shù)值關(guān)系,我們使用點(diǎn)源進(jìn)行模擬,所選用用的點(diǎn)源包括自然射線和人工射線兩大類,最低能量值為0.352mev,最大的能量值2.62mev(詳見表3-5) 7,8。 圖3-1 nai(tl)多道能譜儀及其探頭部分3.2.2 探測器柵元卡片填寫模擬nai(tl)探測器內(nèi)部結(jié)構(gòu)示意見
33、圖3-2,探測器分6個(gè)柵元,使用15個(gè)界面組成。下圖中阿拉伯?dāng)?shù)字表示該曲面的曲面號(hào),圓圈數(shù)字表示各柵元號(hào)。表3-3為探測器各部分的尺寸大小8。 圖3-2 nai(tl)探測器內(nèi)部結(jié)構(gòu)示意圖表3-3 探測器各部分尺寸大小單位(cm)內(nèi)徑外徑厚度高nai-3.8100-7.062mgo3.81003.99500.18507.850sio2-3.9950-0.300fe-3.9950-0.300air4.04504.36460.3196-fe4.36464.44060.07609.705柵元卡填寫:1 1 -3.67 1 -2 -3 imp:p=1 imp:e=12 2 -2.00 (5 -6 -3
34、)#1 imp:p=1 imp:e=13 3 -2.65 3 -6 -4 imp:p=1 imp:e=14 4 -7.93 (7 -8 -9)#1#2#3 imp:p=1 imp:e=15 5 -0.001293 (10 -11 -12)#1#2#3#4 imp:p=1 imp:e=16 6 -7.93 (13 -14 -15)#1#2#3#4#5 imp:p=1 imp:e=1第一列表示柵元號(hào)(參考圖 3-2);第二列表示各柵元對(duì)應(yīng)的材料號(hào);第三列表示該材料的密度;第四列表示該柵元組成的曲面號(hào)(參考圖 3-2),imp:p=1 表示該柵元中光子的重要性為 1,同樣 imp:e=1 表示該柵元
35、中電子的重要性為 1。對(duì)應(yīng)各個(gè)曲面的幾何位置,第一列為曲面號(hào);第二列中 pz 表示垂直 z 軸的平面,cz 表示軸心在 z 軸的圓柱;第三列表示對(duì)應(yīng)第二列方程數(shù)據(jù)項(xiàng)。1 pz 1.8112 cz 3.813 pz 9.4314 pz 9.7315 pz 1.6266 cz 3.9957 pz 1.5768 cz 4.0459 pz 9.78110 pz 0.07611 cz 4.364612 pz 11.28113 pz 014 cz 4.440615 pz 11.357探測器材料組成見表3-3表3-4 探測器組成材料及其各元素含量-物質(zhì)密度(g/cm3)成分比探測器nai&tli3.67na
36、:-0.1532i:-0.8462tl:-0.0006mgo2.00mg:0.5o:0.5-sio22.65si:0.3333o:0.6667-fe&c7.93fe:-0.99995c:-0.00005-空氣air0.001293n:-0.7562o:-0.2438-注:成分比部分帶負(fù)值表示物質(zhì)中元素的質(zhì)量比,正值表示物質(zhì)中元素的原子成分比材料卡描述:m1 11023 -0.15318 53127 -0.846201 81204 -0.000616m2 12000 0.5 8016 0.5m3 14000 0.33333 8016 0.66667m4 26000 -0.99995 6012 -
37、0.00005m5 8016 -0.24377 7014 -0.75623第一列表示表示材料號(hào);第二列表示對(duì)應(yīng)各材料的物質(zhì)組成比例;以上材料分別對(duì)應(yīng) nai(tl)晶體、反射體 mgo、有機(jī)玻璃 sio2、不銹鋼罐以及空氣12。3.2.3 模擬模型建立點(diǎn)源模擬系統(tǒng),見圖 3-3,建立三維坐標(biāo)系,將探測器置于坐標(biāo)原點(diǎn),原點(diǎn)正下方 10cm 處放置點(diǎn)源,現(xiàn)以cs137源為例。以原點(diǎn)為中心建立一個(gè)半徑為 150cm 的球體屏蔽系統(tǒng),以便快速結(jié)束粒子追蹤,球體填充空氣,密度大小為 0.001293 g/cm3,球體外為真空8。圖3-3 點(diǎn)源模擬系統(tǒng)示意剖面圖3.2.4 源定義我們選取多種射線能量的進(jìn)行
38、模擬,本文中我們選取自然射線和人工放射性兩大類射線。各個(gè)能量和元素見下表8。表3-5 各個(gè)核素與能量自然射線核素能量(mev)pb-2140.352bi-2140.609bi-2141.12k-401.46bi-2141.76ti-2082.62人工射線ba-1310.356ba-1310.389cs-1370.662co-601.173co-601.332na-241.369下文中我們用137cs源進(jìn)行模擬示范。137cs源為各向同性點(diǎn)源,能量為 662kev,位置在原點(diǎn)正下方 10cm 處。發(fā)射粒子權(quán)重為 1,發(fā)射源的粒子類型是光子,同時(shí)考慮光子作用。描述如下12:mode p esdef
39、 erg=0.662 pos=0 0 -10 par=2 wgt=13.2.5 輸入卡片# nai探測器點(diǎn)源137cs模擬81 1 -3.67 1 -2 -3 imp:p=1 imp:e=1 $nai晶體2 2 -2.00 (5 -6 -3)#1 imp:p=1 imp:e=1 $mgo3 3 -2.65 3 -6 -4 imp:p=1 imp:e=1 $sio24 4 -7.93 (7 -8 -9)#1#2#3 imp:p=1 imp:e=1 $不銹鋼殼5 5 -0.001293 (10 -11 -12)#1#2#3#4 imp:p=1 imp:e=1 $空氣6 6 -7.93 (13 -1
40、4 -15)#1#2#3#4#5 imp:p=1 imp:e=1 $不銹鋼殼7 7 -0.001293 -16#1#2#3#4#5#6 imp:p=1 imp:e=1 $空氣8 0 16 imp:p=0 imp:e=0 $真空1 pz 1.811 $曲面定義2 cz 3.813 pz 9.4314 pz 9.7315 pz 1.6266 cz 3.9957 pz 1.5768 cz 4.0459 pz 9.78110 pz 0.07611 cz 4.364612 pz 11.28113 pz 014 cz 4.440615 pz 11.35716 so 150mode p e $源定義sdef
41、 erg=0.662 pos=0 0 -10 par=2 wgt=1m1 11023 -0.15318 53127 -0.846201 81204 -0.000616 $材料定義卡m2 12000 0.5 8016 0.5m3 14000 0.33333 8016 0.66667m4 26000 -0.99995 6012 -0.00005m5 8016 -0.24377 7014 -0.75623m6 26000 -0.99995 6012 -0.00005m7 8016 -0.24377 7014 -0.75623f8:p 1 $記數(shù)卡e8 0 1e-5 0.01 400i 1.0 $記數(shù)
42、能量卡ft8 geb 0 0.0252 $記數(shù)特殊處理卡nps 50000000 $計(jì)算粒子數(shù)截?cái)嗫╬rint $打印輸出結(jié)果東華理工大學(xué)畢業(yè)設(shè)計(jì)(論文) 計(jì)算結(jié)果與處理4 計(jì)算結(jié)果與處理4.1 模擬輸出結(jié)果4.1.1 自然射線的模擬輸出結(jié)果本次模擬粒子數(shù)為50000000pb-214能量為0.352mev的能譜圖如下:圖4-1214pb0.352mev能譜圖bi-214能量為0.609mev的能譜圖如下:圖4-2215bi0.609mev能譜圖 bi-214能量為1.12mev的能譜圖如下:圖4-3 214bi1.12mev能譜圖k-40能量為1.46mev的能譜圖如下:圖4-4 40k1.
43、46mev能譜圖bi-214能量為1.76mev的能譜圖如下:圖4-5 214bi1.76mev能譜圖ti-208能量為2.26mev的能譜圖如下:圖4-6 208ti2.26mev能譜圖對(duì)于能量大于1.02mev的射線與探測介質(zhì)相互作用就有可能產(chǎn)生電子對(duì)事件,在此事件中,正電子湮滅產(chǎn)生兩個(gè)0.511mev的光子,一個(gè)或兩個(gè)均逃離探測介質(zhì),探測介質(zhì)吸收余下的能量,則產(chǎn)生單逃逸峰和雙逃逸峰。在實(shí)驗(yàn)結(jié)果中這一現(xiàn)象都表現(xiàn)出來了,特別是在1.5mev以上,譜形上又出現(xiàn)單逃逸峰和雙逃逸峰13,見圖4-5和圖4-6中。圖4-1到4-6是所選的自然射線能譜圖,根據(jù)測量各個(gè)能譜上全能峰上的相對(duì)計(jì)數(shù)率,與對(duì)應(yīng)的
44、能量,列出下表:表4-1 各個(gè)元素能量與對(duì)應(yīng)的計(jì)數(shù)率元素能量(mev)相對(duì)計(jì)數(shù)率pb-2140.3520.0132bi-2140.6090.00854bi-2141.120.00513k-401.460.00409bi-2141.760.00350ti-2082.620.002424.1.2 人工射線的模擬輸出結(jié)果以下是所選取的人工射線的能譜圖:ba-131能量為0.356mev的能譜圖如下: 圖4-7 131ba0.356mev能譜圖ba-131能量為0.389mev的能譜圖如下: 圖4-8 131ba0.389mev能譜圖cs-137能量為0.662mev的能譜圖如下: 圖4-9 137c
45、s0.662mev能譜圖co-60能量為1.173mev的能譜圖如下: 圖4-10 60co1.173mev能譜圖co-60能量為1.332mev的能譜圖如下: 圖4-11 60co1.332mev能譜圖na-24能量為1.369mev的能譜圖如下: 圖4-12 24na1.369mev能譜圖圖4-7到4-12是所選的人工射線能譜圖,根據(jù)測量各個(gè)能譜上全能峰上的相對(duì)計(jì)數(shù)率,與對(duì)應(yīng)的能量,列出下表:表4-2 各個(gè)元素能量與對(duì)應(yīng)的計(jì)數(shù)率元素能量(mev)相對(duì)計(jì)數(shù)率ba-1310.3560.0131ba-1310.3890.0123cs-1370.6620.00796co-601.1730.00491co-601.3320.00442na-241.3690004314.2 數(shù)據(jù)處理與結(jié)果根據(jù)表4-1與4-2可分別模擬出人工射線與自然射線的射線能量和閃爍探測器的探測效率之間的數(shù)值模擬圖和擬合關(guān)系式。自然射線能量值和相對(duì)計(jì)數(shù)率的模擬圖如下: 圖4-13 自然射線能量值和相對(duì)計(jì)數(shù)率的模擬圖則根據(jù)圖4-13擬合出指數(shù)衰減函數(shù)下面關(guān)系:p=0.02052*exp(-e/0.51684)+0.00263 (4-1)人工射線能量值和相對(duì)計(jì)數(shù)率的模擬圖如下:圖4-14 人工射線能量值和相對(duì)計(jì)數(shù)率的模擬圖則根據(jù)圖4-14擬合出指數(shù)衰減函數(shù)下面關(guān)系:p=0.0233*exp(-e/0
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 路邊廣告位轉(zhuǎn)讓合同
- 美國自費(fèi)出國留學(xué)咨詢服務(wù)合同年
- 居間合同傭金承諾書
- 事故車買賣合同協(xié)議
- 連車帶人租賃合同
- 荒山承包合同范本
- 叉車租賃合同協(xié)議書范本大全
- 工地材料運(yùn)輸合同
- 借款合同答辯狀范本范本
- 個(gè)人工作總結(jié)范文20篇
- 2024年廣東省公務(wù)員錄用考試《行測》真題及解析
- 高中英語必背3500單詞表(完整版)
- 禁止送禮的協(xié)議書
- 2024年版《輸變電工程標(biāo)準(zhǔn)工藝應(yīng)用圖冊》
- 2024年高考數(shù)學(xué)試卷(北京)(空白卷)
- 2024從洞見到生意:阿里健康特色人群消費(fèi)趨勢報(bào)告-阿里健康x一財(cái)商學(xué)院
- 人教版2024年新教材七年級(jí)上冊英語starter unit 1 -unit7重點(diǎn)短語句型清單
- 護(hù)理服務(wù)在產(chǎn)科中的應(yīng)用課件
- 2024年小升初語文入學(xué)分班測試卷四(統(tǒng)編版)
- 流行文化對(duì)青少年價(jià)值觀的影響研究
- 中國保險(xiǎn)行業(yè)協(xié)會(huì)官方-2023年度商業(yè)健康保險(xiǎn)經(jīng)營數(shù)據(jù)分析報(bào)告-2024年3月
評(píng)論
0/150
提交評(píng)論