反应堆关键材料空位、间隙演化的速率理论并行模拟方法与流程

本发明涉及核材料辐照损伤计算机模拟技术，具体涉及一种反应堆关键材料空位、间隙演化的速率理论并行模拟方法。

背景技术：

速率理论(rt，ratetheory)是一种用于模拟材料辐照损伤长时间微观结构演化的重要方法，具有计算速度快、计算效率高且无时空尺度限制等优点，可模拟上百dpa下的缺陷演化过程，获得的缺陷速密度及尺寸分布可与实验直接进行对比，在辐照诱导微观结构演化研究中得到了广泛的应用，如材料辐照肿胀模拟、辐照生长计算、团簇析出研究等。近30年来，国际上采用速率理论成功解释/预测了辐照肿胀、辐照生长和辐照导致团簇析出等辐照效应。

速率理论是一种平均场方法，也称主方程或团簇动力学，它将材料微观结构视为空间均匀的有效介质，认为缺陷在空间中均匀分布。速率理论利用缺陷反应速率表示缺陷之间的相互作用，通过构建缺陷浓度(即缺陷数密度)随时间的演化方程来描述辐照诱导材料微观缺陷的演化过程，其模型被表述为一系列常微分方程(odes)，用来描述点缺陷和缺陷团簇的产生、消除，以及相应的微观结构演化。各种缺陷数密度cn随时间演化的一般性方程如下：

其中，cj表示尺寸为j的缺陷团簇的浓度(即数密度)；方程右边第一项表示级联碰撞过程产生的尺寸为j的缺陷团簇的速率；右边第二项表示其他尺寸的缺陷团簇反应生成尺寸为j的团簇的反应速率；右边第三项表示尺寸为j的缺陷团簇反应生成其他尺寸团簇的反应速率；右边第四项为材料内固有缺陷(如：位错、晶界等)对尺寸为j的缺陷团簇的吸收速率。

从rt模型上可以看出，一种缺陷对应一个常微分方程，当模拟一个含有n种缺陷的体系时，则需要求解n个刚性、病态的常微分方程来获得缺陷团簇的尺寸和数密度分布，速率方程数量通常为百万量级，且随着团簇类型和尺寸的增加，方程数呈指数级增长，当同时考虑空位和间隙的演化时，方程中还会增加间隙和空位的复合项。因此，rt模拟的难点在于常微分方程组的高效、精确求解上。

技术实现要素：

本发明的目的在于针对速率理论模拟方程刚性、大量、求解难的问题，提供一种基于指数时间差分(etd，exponentialtimedifferencing)算法的反应堆关键材料空位、间隙演化的速率理论并行模拟方法，从而加快方程的求解速度，实现速率方程的高效求解，以及空位、间隙的大规模快速模拟。

本发明的技术方案如下：一种反应堆关键材料空位、间隙演化的速率理论并行模拟方法，包括：

(1)将点缺陷及缺陷团簇的浓度随时间演化的方程，转化为带预估矫正的二阶指数时间差分格式，如下：

预估项：c^*＝e^δtlck+(e^δtl-i)f(ck，tk)/l

矫正项：ck+1＝e^δtlc^*+(e^δtl-i)(f(c^*，tk)-f(ck，tk))/2l

其中，ck为第k步迭代的缺陷x的浓度，δt为时间步长，tk为第k步的迭代时间，f(ck，tk)为缺陷x的速率方程的非线性部分，是一个2n×1矩阵，l是缺陷x的速率方程的线性部分，是一个2n阶稀疏矩阵，包含缺陷x的吸收速率和发射速率，i是2n阶单位矩阵；

(2)采用数据并行的方式，将指数时间差分格式的方程组进行划分，根据使用的进程数n，将总数为2n的方程组(含n个空位团簇演化方程、n个间隙团簇演化方程)划分为n组，每个进程求解个方程；

(3)由0号进程将点缺陷浓度c1v和c1i(v为空位、i为间隙原子、1为团簇尺寸)发送给其他进程，每个进程只求解自己所负责的m个方程；

(4)每一次迭代结束后，相邻的进程之间进行点对点通信，进程p将自己的最小尺寸团簇的浓度、发射速率和吸收速率发送给进程p-1，将最大尺寸团簇的浓度、发射速率和吸收速率发送给进程p+1，其中，p∈{0，...，n-1}。

进一步，如上所述的反应堆关键材料空位、间隙演化的速率理论并行模拟方法，步骤(1)中所述的点缺陷包括空位v、间隙原子i，缺陷团簇包括空位团簇、间隙团簇。

进一步，如上所述的反应堆关键材料空位、间隙演化的速率理论并行模拟方法，步骤(2)中将每个进程中原大小为m的数据块增加到m+2，以便于存储其它进程传来的数据。

进一步，如上所述的反应堆关键材料空位、间隙演化的速率理论并行模拟方法，步骤(4)中相邻的进程之间进行点对点通信时，引入虚拟进程实现0号进程、n－1号进程的通信。

本发明的有益效果如下：本发明通过将隐式方程组转化为显式格式，降低了方程求解的难度，同时使求解更加精确。通过将不同尺寸的缺陷演化方程进行并行求解，提高了求解速度，使模拟更加快速。在配备intel(r)xeon(r)cpue5-2620v3@2.40ghzcpu的集群上，求解1600个方程，使用16核的加速达5.7。

附图说明

图1是本发明的反应堆关键材料空位、间隙演化速率理论并行模拟流程图。其中，t表示当前的时间步，time表示模拟的总时长。δt为模拟的时间步长，“反应”包括：缺陷团簇吸收一个点缺陷、缺陷团簇发射一个点缺陷、空位与间隙的复合反应。

图2是本发明的数据划分示意图。其中p0～p3表示进程编号，c1(i，v)～cm(i，v)表示各进程负责的缺陷的尺寸范围。

图3是本发明的通信示意图。其中每一小格代表一个尺寸的团簇数据，虚线框存的是相邻尺寸的团簇数据。

图4是本发明方法的强可扩展性测试结果。

具体实施方式

下面结合附图和实施例对本发明进行详细的描述。

从速率理论rt模拟的一般性方程上可以看出，用rt方法模拟辐照缺陷的演化时，需要求解大量刚性的，且耦合的常微分方程组。传统串行算法虽然具有较高的准确性，但随着方程数量的增加，其模拟时间变得非常长，限制了所能模拟的缺陷数量。因此，本发明提供了一种基于指数时间差分的并行算法，该算法能够在保证准确性的基础上，同时利用并行计算技术，扩大rt模拟的缺陷数量，快速求解rt方程组。

针对速率方程刚性、大量、求解难的问题，本发明提供了一种基于指数时间差分(etd，exponentialtimedifferencing)算法的反应堆关键材料空位、间隙演化的速率理论并行模拟技术，通过将刚性微分方程转化为显式格式的方程，降低求解难度，采用数据并行的方式，将一组常微分方程(含n个空位团簇演化方程、n个间隙团簇演化方程)分为n组，每组m个方程(其中，)，分别由n个进程并行求解，从而加快方程的求解速度，实现速率方程的高效求解，以及空位、间隙的大规模快速模拟。

如图1所示，具体方法包括如下步骤：

(1)将点缺陷(空位v、间隙i)及缺陷团簇(空位团簇、间隙团簇)的浓度随时间演化的方程，转化为带预估矫正的二阶etd格式，如下：

预估项：c^*＝e^δtlck+(e^δtl-i)f(ck，tk)/l

矫正项：ck+1＝e^δtlc^*+(e^δtl-i)(f(c^*，tk)-f(ck，tk))/2l

其中，ck为第k步迭代的缺陷x的浓度，δt为时间步长，tk为第k步的迭代时间，f(ck，tk)为缺陷x的速率方程，是一个2n×1矩阵，l为缺陷x的吸收速率和发射速率之和，是一个2n阶矩阵。

(2)初始化模拟参数，包括模拟的缺陷最大规模2n、模拟的总时长time、时间步长δt等，其中时间步长设置为δt＝1×10^-6s，模拟温度、缺陷的扩散系数、扩散能量、结合能等参数。

(3)采用数据并行方式，将2n个方程划分到n个进程，这里假设需要计算20个不同尺寸的团簇浓度，进程数为4，则n＝10，n＝4，见附图2。主要步骤为：

a.获取总进程数；

b.每个进程负责5个方程的求解，数据块大小为5。

(4)判断当前模拟时间t是否小于总模拟时间time，

是：转步骤(5)；

否：结束模拟，输出结果。

(5)每个进程计算自己所负责缺陷的反应速率总和。若是其他进程，则将总反应速率发送给0号进程；若是0号进程，则对每个进程发来的总反应速率进行规约操作。

(6)各进程计算间隙/空位团簇的预估浓度c^*和矫正浓度ck+1(预估浓度对应etd格式中的预估项公式计算出的团簇浓度，不是每步迭代计算的最终团簇浓度；矫正浓度对应etd格式中的矫正项公式计算出的团簇浓度，矫正浓度是每步迭代要计算的最终团簇浓度)。附图2中数据块两边各增加一小块用来存储相邻数据块通信所得的数据，每个数据块大小增加至7，见附图3。其他进程接收上一时间步0号进程发来的点缺陷浓度，然后计算所负责缺陷团簇的浓度。相邻的进程之间使用mpi_sendrecv()函数进行点对点通信，实现捆绑发送/接收，避免死锁的产生，即进程p将自己的最小尺寸的团簇信息发送给进程p-1，将最大尺寸的团簇信息发送给进程p+1。通过建立虚拟进程实现0号进程、3号进程的通信，见附图3。

(7)0号进程将点缺陷浓度发送给其他进程。

(8)模拟时间前进δt，转步骤(4)。

显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若对本发明的这些修改和变型属于本发明权利要求及其同等技术的范围之内，则本发明也意图包含这些改动和变型在内。