[發明專利]一種基于粒子密度不均勻性的迭代蒙卡全局權窗參數生成方法

有效

• 申請號: 201610685275.2
• 申請日: 2016-08-18
• 公開（公告）號: CN106295213B
• 公開（公告）日: 2017-05-24
• 發明（設計）人: 吳宜燦;張澍;宋婧;郝麗娟;孫光耀
• 申請（專利權）人: 中國科學院合肥物質科學研究院
• 主分類號: G06F19/00
• 分類號: G06F19/00
• 代理公司: 合肥市上嘉專利代理事務所(普通合伙)34125
• 代理人: 王偉

1.一種基于粒子密度不均勻性的迭代蒙卡全局權窗參數生成方法，其特征在于包括以下步驟：

(1)生成初始權窗網格：

(11)建立覆蓋整個模型空間的網格；

(12)將所有網格的權窗參數都初始化為0；

(2)在粒子輸運計算過程中，統計進入每個網格的粒子數量N_i、總粒子權重W_i和每個網格對粒子密度均勻性的貢獻，以及所有網格的平均粒子數量Na，i指第i個網格；在每個粒子計算過程中需要：

(21)記錄粒子的徑跡信息，即粒子所穿過的所有網格的編號；

(22)當粒子第一次進入網格i后，如果粒子權重w大于當前所有網格權窗參數的最小值，將進入網格i的粒子數量N_i增加1，進入網格i的總粒子權重W_i增加當前粒子的權重w；

(23)將粒子每一步輸運中，終點所在的網格的粒子數量N_i和所有網格的平均粒子數量Na做比較，如果N_i小于Na，則該粒子徑跡上的所有網格的貢獻增加，增加的貢獻即為粒子進入一個網格i對粒子密度均勻性的貢獻；

(3)計算權窗參數

(31)計算每個網格的重要性，

${\mathrm{Iwp}}_i=\frac{C_i}{W_i}$

計算每個網格的重要性Imp_i，公式中C_i為所有計算的粒子進入網格i對粒子密度均勻性的總貢獻，W_i為進入網格i的總粒子權重；

(32)將每個網格的重要性根據源所在網格的平均重要性進行歸一化；

(33)將每個網格的權窗參數設為歸一化后的網格重要性的倒數的一半；

(4)迭代優化

(41)在權窗生成過程中，將粒子輸運分為若干個周期，根據前一個周期的網格重要性生成新的權窗參數，并應用于當前周期的輸運，同時每個周期生成的權窗參數都通過文本進行輸出；

(42)在每個周期粒子輸運開始時，將每個網格的貢獻和總粒子權重都乘以一個折扣因子f；

(43)迭代自動進行，在達到了用戶設定的時間上限后停止計算，當前的權窗參數即為最終生成的全局權窗參數。