阅读说明：本技术 一种辐照缺陷和纳米相协同演化的相场模拟方法 (Phase field simulation method for irradiation defect and nanophase co-evolution ) 是由 李永胜 严正伟 于 2020-06-28 设计创作，主要内容包括：本发明为一种辐照缺陷和纳米相协同演化的相场模拟方法。通过合金热力学自由能来表征系统的演变；建立辐照缺陷–合金体系的总自由能；创建辐照缺陷和纳米相共同演化的相场方程,实现辐照下合金缺陷和纳米相演化的模拟与表征；分析辐照温度和缺陷浓度变化条件下的缺陷和纳米相微观组织形貌的变化,得出辐照下缺陷和纳米析出相的演变规律；同时分析缺陷和纳米相的体积分数变化,得出不同辐照温度和缺陷浓度对缺陷和纳米相的形成及演化动力学的影响规律；本发明提供了一种辐照缺陷和纳米相协同演化的相场模拟方法,能够获得辐照下缺陷和相分解之间的相互作用关系,为辐照下合金材料的设计和服役安全提供指导。(The invention relates to a phase field simulation method for irradiation defect and nanophase co-evolution. The evolution of the system is characterized by the thermodynamic free energy of the alloy; establishing total free energy of an irradiation defect-alloy system; establishing a phase field equation of irradiation defect and nanophase co-evolution, and realizing simulation and characterization of alloy defect and nanophase evolution under irradiation; analyzing the defects under the condition of changing the irradiation temperature and the defect concentration and the change of the shape of the nanophase microstructure to obtain the evolution law of the defects and nanophase precipitated phases under irradiation; analyzing the volume fraction change of the defects and the nanophase simultaneously to obtain the influence rule of different irradiation temperatures and defect concentrations on the formation and evolution dynamics of the defects and the nanophase; the invention provides a phase field simulation method for irradiation defect and nanophase co-evolution, which can obtain the interaction relation between irradiation defect and phase decomposition and provide guidance for the design and service safety of an irradiated alloy material.)

一种辐照缺陷和纳米相协同演化的相场模拟方法

技术领域

本发明涉及金属材料辐照下微观组织控制技术领域，具体为一种辐照缺陷和纳米相协同演化的相场模拟方法。

背景技术

在核领域，核反应堆结构材料在高温和辐照下，合金中的原子将会离开原来晶格位置而形成空位和间隙原子，缺陷的存在促进了纳米相的析出，阻碍了位错运动引起合金脆硬性的增加，导致材料的使用寿命降低。因此研究辐照下缺陷和纳米相的演化对合金的安全使用尤为重要。

从目前的研究中可以看出辐照对合金相转变的影响因素主要有辐照温度和辐照剂量，辐照温度的升高使缺陷的迁移作用加强，原子扩散作用增强，继而促进了纳米相转变。辐照剂量的增加促使级联作用增强，缺陷数量的增加为原子扩散提供通道，加速原子扩散，使相分解加快。由此可知，辐照剂量和辐照温度都是通过加强缺陷作用来促进相分解，但是对于辐照下缺陷对纳米相演化的动力学影响规律确还不清楚。

目前较多的模拟研究都只针对辐照下缺陷或者纳米相单个的演变，对于辐照下缺陷和纳米相的协同演化的研究目前还很缺少，在相场法模拟中可以将缺陷简化成同等的元素变量，来模拟缺陷和纳米相的协同演变。同时利用相场法研究还可以弥补实验周期长，耗费成本大，研究对象单一和不易实时观测的不足。因此，利用相场模拟手段研究辐照下缺陷和纳米相协同演化是非常有优势且有必要的。

发明内容

本发明所解决的技术问题在于提供一种辐照缺陷和纳米相协同演化的相场模拟方法，解决现有研究方法无法同时研究缺陷和纳米相的问题。

实现本发明目的提供技术方案如下：

一种辐照缺陷和纳米相协同演化的相场模拟方法：

(1)利用合金热力学自由能进行缺陷和纳米相系统的热力学描述；

(2)根据热力学描述，求解辐照下缺陷和纳米相系统的化学自由能和弹性能；

(3)基于Chan-Hilliard扩散方程，建立辐照下缺陷和纳米相演变的相场动力学方程；

(4)根据输入变量来模拟不同辐照条件下缺陷和纳米相的组织形貌以及成分演变，分析不同辐照条件下缺陷对相分解动力学规律的影响。

具体内容如下：

根据合金热力学特征对(1)中的热力学进行描述，主要包括：系统的化学自由能、弹性能、梯度能，具体如下：

其中F为系统总的自由能，f_chem为合金的化学自由能，V_m为合金的摩尔体积，i代表合金元素，c_i为元素i的原子分数，c_V和c_I分为空位和间隙原子的原子分数，k_j(j＝i,V,I)为梯度能系数，f_elas为弹性能密度。

在本发明中将缺陷的形成能及其浓度耦合到系统的化学自由能中，将(2)中的化学自由能表示为：

其中和分别为空位和间隙原子形成能，为元素i的吉布斯自由能，k代表合金元素，c_k为元素k的原子分数，L_ik为原子间相互作用参数。

在本发明中将缺陷和应变所引起成分的不均匀耦合到本征应变中，(2)中的弹性能密度可表示如下：

其中C_ijkl为弹性模量，为本征应变具体公式为：

其中ε₀，ε_V0，和ε_I0分别为成分不均匀以及空位和间隙原子引入产生的晶格膨胀系数，弹性应变ε_kl与位移有关，其表达式为：

通过将缺陷有关量耦合到化学自由能和弹性能中，便于利用相场法来描述缺陷的演变及其对相分解的影响。

对于本发明的步骤(3)中所描述的辐照缺陷和纳米相协同演变的动力学方程可表示为：

其中M为合金的化学迁移率，M_j(j＝V,I)为空位和间隙原子迁移率，ξ_j(j＝i,V,I)为热力学起伏，P_j(j＝V,I)和S_j(j＝V,I)分别为缺陷的产生速率和湮灭速率，R_VI为缺陷的合并速率。在相场模拟计算中利用半隐式傅里叶谱算法对动力学演变方程进行求解，得出成分与时间和空间位置关系，通过定量的分析来完成(4)中不同辐照条件下缺陷和纳米相组织形貌特征，同时探寻缺陷对相分解动力学的影响规律。

对比于现有技术，本发明的显著优点如下：

(1)该方法通过将缺陷的形成能和浓度耦合到系统的自由能中，使得辐照点缺陷和纳米相能够同时被捕获，弥补了现有模拟方法只能研究单一因素的局限性，同时为合金辐照效应的研究提供了具体有效的方法；

(2)该方法考虑纳米相转变的同时结合了晶体缺陷的本质，通过Chan-Hilliard扩散方程，建立了空位和间隙原子缺陷和纳米相协同演化的动力学方程，得到一种辐照缺陷和纳米相共同演化的相场模拟方法，能更好的解析不同辐照条件对缺陷和纳米相演变的影响，以及缺陷和纳米相间的相互影响。

附图说明

图1为Fe-35at.％Cr合金中α′相、空位和间隙原子的微观组织图，辐照温度为750K，缺陷初始浓度为10^-5，辐照时间t*分别为44、52和482。

图2为不同辐照温度下Fe-35at.％Cr合金中α′相、空位和间隙原子的微观组织图，辐照温度分别为710K，730K和760K，缺陷初始浓度为10^-5，辐照时间t*为382。

图3为不同初始缺陷浓度下Fe-35at.％Cr中α′相体积分数随时间的变化，缺陷初始浓度分别为0、10^-7、10^-5和10^-4，辐照温度为750K。

具体实施方式

为了更加清楚的说明本发明的具体实施方式，通过附图具体的实施例来进一步加以说明。

本发明为一种辐照缺陷和纳米相协同演化的相场模拟方法，其步骤在于：

(1)利用合金热力学自由能进行缺陷和纳米相系统的热力学描述；

(2)根据热力学描述，求解辐照下缺陷和纳米相系统的化学自由能和弹性能；

(3)基于Chan-Hilliard扩散方程，建立辐照下缺陷和纳米相演变的相场动力学方程；

(4)根据输入变量来模拟不同辐照条件下缺陷和纳米相的组织形貌以及成分演变，分析不同辐照条件下缺陷对相分解动力学规律的影响。

以Fe-Cr合金为例。

步骤(1)中，Fe-Cr合金的热力学描述包括，系统的化学自由能、弹性能、梯度能，具体公式如下：

其中F为系统总的自由能，f_chem为合金的化学自由能，V_m为合金的摩尔体积，i代表合金元素，c_i为元素i的原子分数，c_V和c_I分为空位和间隙原子的原子分数，k_j(j＝i,V,I)为梯度能系数，f_elas为弹性能密度。

将缺陷的形成能及其浓度耦合到系统的化学自由能中，步骤(2)中Fe-Cr合金的化学自由能可表示为：

(i＝Fe,Cr；k＝Fe,Cr；i≠k)

其中

和

分别为空位和间隙原子形成能，

为元素i的吉布斯自由能，k代表合金元素，c_k为元素k的原子分数，L_ik为原子间相互作用参数。

将缺陷和应变所引起成分的不均匀耦合到本征应变中，步骤(2)中Fe-Cr合金的弹性能密度可表示如下：

其中C_ijkl为弹性模量，为本征应变具体公式为：

其中ε₀，ε_V0，和ε_I0分别为成分不均匀以及空位和间隙原子引入产生的晶格膨胀系数，弹性应变ε_kl与位移有关，其表达式为：

步骤(3)中，Fe-Cr合金辐照缺陷和相协同演变的动力学方程可表示为：

其中M为Fe-Cr合金的化学迁移率可表示为M＝[c_CrM_Fe+(1-c_Cr)M_Cr]c_Cr(1-c_Cr)，M_j(j＝V,I)为空位和间隙原子迁移率，从动力学方程中可知缺陷和原子可通过迁移率来建立内在联系，而空位扩散机制表明原子通过与空位交换来完成其迁移，因此该相场方程在理论上是符合的。ξ_j(j＝i,V,I)为热力学起伏，P_j(j＝V,I)和S_j(j＝V,I)分别为缺陷的产生速率和湮灭速率，R_VI为缺陷的合并速率，速率的引入使得相场法更加符合实际。在模拟计算过程中利用半隐式傅里叶谱算法对动力学演变方程进行求解。

步骤(4)，设置Fe-Cr合金初始成分为c_Cr＝0.35，空位和间隙原子初始浓度范围为10^-7～10^-4，辐照温度范围为710～760K，计算网格大小为128△x^*×128△y^*，根据模拟计算得到的数值绘制Fe-Cr合金在不同辐照条件下的微观组织演化图。

图1为Fe-35at.％Cr合金中α′相、空位和间隙原子的微观组织图，辐照温度为750K，缺陷初始浓度为10^-5，辐照时间t*分别为44、52和482。从图中可以知，随着辐照时间的延长，α′相在高缺陷浓度处开始聚集，同时缺陷逐渐从中间向相边界聚集，在两相间形成缺陷环。这是因为在空位扩散机制中，原子易于和空位进行交换继而完成其迁移，所以辐照产生的大量空位能够促进相分解。

图2为不同辐照温度下Fe-35at.％Cr合金中α′相、空位和间隙原子的微观组织图，辐照温度分别为710K，730K和760K，缺陷初始浓度为10^-5，辐照时间t*为382。从图中可知，随着辐照温度的增加，α′相长大粗化加剧，颗粒增大，同时缺陷也逐渐从初始位置向边界聚集，并且形成环。

图3为不同初始缺陷浓度下Fe-35at.％Cr中α′相体积分数随时间的变化，缺陷初始浓度分别为0、10^-7、10^-5和10^-4，辐照温度为750K；从图中可以看出，随着缺陷浓度的升高，α′相体积分数开始变化的时间提前，相分解时间提前，相变受到促进。

总的来说，辐照温度和缺陷浓度都可以促进相分解，并且温度越高作用越强。这可作为辐照下核反应堆结构材料设计应用时的理论参考，并以此优化合金的组织性能。

以上所述是对本发明的优选实施案列及原理进行了详细的说明，对于本邻域技术人员依据本发明所提供的思想，在具体实施上的差别都将视为本发明的保护范畴。