原子動(dòng)力學(xué)蒙特卡洛程序OpenKMC在反應(yīng)堆壓力容器鋼缺陷損傷研究中的優(yōu)化與應(yīng)用
【文章頁(yè)數(shù)】:12 頁(yè)
【部分圖文】:
圖1 bcc晶格體心原子的三重近鄰
KMC模擬方法是一種通過構(gòu)造隨機(jī)過程來模擬體系長(zhǎng)時(shí)間演化的方法。經(jīng)典的原子動(dòng)力學(xué)蒙特卡洛AKMC(AtomicKineticMonteCarlo)方法[10]利用晶格點(diǎn)陣來映射原子與格點(diǎn)的位置關(guān)系,應(yīng)用于提供控制微觀系統(tǒng)演化的精確解。α-Fe的晶格為體心立方(bcc)結(jié)構(gòu),....
圖2 Pair截?cái)喟霃酵茖?dǎo)立體圖
截?cái)喟霃降耐茖?dǎo)過程在三維坐標(biāo)系中完成。假設(shè)發(fā)生躍遷交換的原子分別為i和j,并在坐標(biāo)系中固定i和j的位置坐標(biāo)。確定Pair勢(shì)能模型下最遠(yuǎn)處需要被更新的空位位置可以分為下面3個(gè)步驟:(1)找到j(luò)原子附近2層近鄰內(nèi)的全部原子并記錄坐標(biāo);(2)分別計(jì)算這些原子與i原子的距離,找到其中一個(gè)....
圖3 Pair截?cái)喟霃酵茖?dǎo)三視圖
圖2Pair截?cái)喟霃酵茖?dǎo)立體圖若記i原子坐標(biāo)為i(1,1,1),則有k(4,4,4),故i與k的距離值dPair=33,并將這個(gè)長(zhǎng)度作為Pair勢(shì)能模型下空位選取的截?cái)喟霃介L(zhǎng)度。同時(shí)本文給出原子間近鄰層數(shù)與距離的關(guān)系如表1所示。表1中數(shù)據(jù)是體系中晶格的邊長(zhǎng)為2的情況下,一....
圖4 EAM截?cái)喟霃酵茖?dǎo)立體圖
通過相同的思路,可以確定EAM勢(shì)能模型下空位更新的截?cái)喟霃。與Pair勢(shì)能模型的區(qū)別在于,EAM勢(shì)能模型的原子間相互作用考慮到第3近鄰位置,因此在尋找空位k步驟的第1步里要找到j(luò)原子附近3層近鄰內(nèi)的全部原子并記錄坐標(biāo),接著執(zhí)行相同的第2步和第3步。計(jì)算出的EAM勢(shì)能模型下的截?cái)喟?...
本文編號(hào):4050040
本文鏈接:http://www.wukwdryxk.cn/projectlw/hkxlw/4050040.html