CN108052745B - 一种模拟辐照缺陷与晶界协同演化的方法 - Google Patents
一种模拟辐照缺陷与晶界协同演化的方法 Download PDFInfo
- Publication number
- CN108052745B CN108052745B CN201711352182.9A CN201711352182A CN108052745B CN 108052745 B CN108052745 B CN 108052745B CN 201711352182 A CN201711352182 A CN 201711352182A CN 108052745 B CN108052745 B CN 108052745B
- Authority
- CN
- China
- Prior art keywords
- grain boundary
- transition
- irradiation
- grain
- defect
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Testing Or Measuring Of Semiconductors Or The Like (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种模拟辐照缺陷与晶界协同演化的方法,包括以下步骤:获得晶界能量与平移向量关系,获得基态对应的平移向量集合与基态间跃迁能垒,分析最易跃迁模式,考察点缺陷对晶界跃迁过程的影响,考察晶界附近空位团簇对晶界跃迁过程的影响,考察辐照缺陷对VTk的影响,建立包含晶界运动的速率表格,辐照缺陷与晶界协同演化的OKMC模拟。本发明可以处理辐照缺陷对晶界运动以及辐照缺陷偏聚对晶界运动动力学的影响、运动的晶界与晶界附近活性低的空位团簇或空洞的作用、辐照缺陷与静态的晶界作用、辐照缺陷微结构演化过程中可能出现的晶界运动及其对微结构演化的影响;尤其适用于模拟复杂缺陷‑晶界交互作用下纳米结构辐照损伤微结构演化过程。
Description
技术领域
本发明涉及核材料辐照损伤模拟领域,尤其涉及一种模拟辐照缺陷与晶界协同演化的方法。
背景技术
材料遭受高能粒子辐照时,会生成扩散活性不同的点缺陷及其团簇,如自间隙原子、空位和这些缺陷的团簇。晶界常常作为缺陷阱,吸收辐照缺陷,影响辐照缺陷微结构的演化和材料最终的服役性能。研究辐照缺陷与晶界相互作用具有科学和工程意义。
在研究辐照缺陷与界面跨尺度相互作用机理时,人们通常将晶界看作稳定的辐照缺陷演化环境,关注晶界对点缺陷行为的影响。这是因为一方面点缺陷活性较高,在辐照缺陷微结构演化中可能扮演较为重要的角色,而且其扩散过程较为简单;另一方面缺陷团簇或者晶界活性较低。然而,随着辐照缺陷的演化,如扩散、偏聚,不仅可能改变了晶界结构,进而改变了其热力学和动力学性质,诱导晶界滑移或迁移,反过来,受辐照缺陷修饰的晶界还可能以一种不同的方式与晶界附近活性低的缺陷团簇发生作用,即辐照缺陷与晶界可能发生协同演化。由于晶界本身自由度很多,其势能面较为复杂,目前还没有很好的模拟技术直接考察这种涉及到晶界运动或相变的缺陷-晶界协同演化过程。
因此,需要开发一种能计算晶界滑移/迁移的热力学和动力学性质的模拟技术。
发明内容
本发明所要解决的技术问题在于提供一种模拟辐照缺陷与晶界协同演化的方法,该方法同时将晶界运动事件纳入到OKMC模拟缺陷演化框架中,以考察晶界运动或相变参与的缺陷-晶界多尺度相互作用,尤其是缺陷团簇与晶界的交互作用。
本发明所要解决的技术问题采用以下技术方案来实现:
一种模拟辐照缺陷与晶界协同演化的方法,包括以下步骤
S1、获得晶界能量与平移向量关系:建立初始晶界模型,设其在平行晶界两个方向的最小周期分别为px和py,将一个晶粒相对于另一个晶粒平移时,相应的平移步长分别为dx、dy,平移向量VT=(i*dx,j*dy),其每个方向的值在0-px、0-py之间变化,每一次平移之后,对整个晶界模型进行充分原子弛豫,完成所有的平移-弛豫之后,得到晶界能量E与平移向量VT的关系E=E(i*dx,j*dy);
S2、获得基态对应的平移向量集合与基态间跃迁能垒:对比每次平移-弛豫后的体系能量,得到能量最低的相应平移向量集合VT1,同时得到平行弛豫后相应的体系构型R集合,对R集合中任意两个构型Ri和Rj,比较其距离,若其距离大于跃迁值时,认为两个态之间发生了跃迁,并调用NEB(Nudged Elastic Band)计算跃迁能垒;
S3、分析最易跃迁模式:分析晶界跃迁能垒与跃迁距离关系,获得跃迁能垒最低的跃迁及其对应的平移向量VTk;
S4、考察点缺陷对晶界跃迁过程的影响;
S5、考察晶界附近空位团簇对晶界跃迁过程的影响;
S6、考察辐照缺陷对VTk的影响;
S7、建立包含晶界运动的速率表格;
S8、进行辐照缺陷与晶界协同演化的OKMC(实体动力学蒙特卡洛)模拟。
进一步改进在于,步骤S1中所述建立初始晶界模型,其包括设置模型尺寸、体系中缺陷类型和缺陷数、模型边界条件、模拟体系温度、空位和自间隙的运动向量。
进一步改进在于,步骤S2中所述的跃迁值为2.0埃。
进一步改进在于,步骤S4具体为:在晶界内部放置一定浓度的空位或自间隙原子,按照平移向量VTk平移晶界模型,采用NEB计算当前条件下的跃迁能垒。
进一步改进在于,步骤S5的具体操作为:将模型体系(一般选则较小的模型体系)周期性复制,在晶界附近典型位置放置典型尺寸的空位团簇,然后以平移向量VTk平移其中一个晶粒,采用NEB方法计算此时的跃迁能垒。
进一步改进在于,步骤S6的具体操作为:当晶界含有偏聚的辐照缺陷或晶界附近存在活性低(通过观察比较其活性)的空位团簇或者空洞时,采用步骤S1-S3中的方法重新确定晶界容易跃迁的向量和能垒。
进一步改进在于,步骤S8的具体操作为:将晶界运动看作一个OKMC事件,经过步骤S1-S3,得到纯晶界以及含有辐照缺陷的晶界跃迁向量及能垒,进而得到相应的速率,然后在OKMC循环中,调用相应的速率,选择、执行事件,并对可能出现的新的缺陷-晶界作用重新参数化相关过程。
本发明的有益效果是:
本发明通过晶界滑移诱导的迁移,得到晶界迁移的方向向量和能垒,进而考察辐照缺陷对晶界运动的影响。同现有的模拟方法相比,能够处理辐照缺陷偏聚对晶界运动动力学的影响,同时可处理运动的晶界与晶界附近活性低的空位团簇或空洞的作用。
本发明既能处理辐照缺陷与静态的晶界作用,又能处理辐照缺陷微结构演化过程中可能出现的晶界运动及其对微结构演化的影响。本发明尤其适用于模拟复杂的缺陷-晶界交互作用下纳米结构辐照损伤微结构演化过程。
附图说明
图1为本发明中计算晶界能量与晶粒间相对滑移向量关系示意图;
图2为本发明中晶界最易跃迁模式计算示意图;
图3为本发明的主要步骤;
图4为本发明算法流程图。
具体实施方式
为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体图示和实施例,进一步阐述本发明。
结合图1至图3所示,本发明提出的一种包含晶界运动的模拟辐照辐照缺陷与晶界作用的新方法包括以下步骤:
S1:获得晶界能量与平移向量关系。对于建好的初始晶界模型,设其在平行晶界两个方向的最小周期分别为px和py,将一个晶粒相对于另一个晶粒平移时,相应的平移步长分别为dx、dy。平移向量VT=(i*dx,j*dy),其每个方向的值在0-px、0-py之间变化。每一次平移之后,对整个晶界模型进行充分原子弛豫。完成所有的平移-弛豫之后,得到晶界能量E与平移向量VT的关系E=E(i*dx,j*dy)。
S2:获得基态对应的平移向量集合与基态间跃迁能垒。对比每次平移-弛豫后的体系能量,得到能量最低的相应平移向量集合VT1。同时得到平行弛豫后相应的体系构型R集合。对R集合中任意两个构型Ri和Rj,比较其距离,若排除原子位置周期性调整带来的影响后,其距离大于跃迁值(2.0埃)时,认为两个态之间发生了跃迁,调用NEB计算跃迁能垒。
S3:分析最易跃迁模式。分析晶界跃迁能垒与跃迁距离关系,获得跃迁能垒最低的跃迁及其对应的平移向量VTk。
S4:考察点缺陷对晶界跃迁过程的影响。此步骤中,对于给定的晶界,在晶界内部放置一定浓度的空位或自间隙原子。按照平移向量VTk平移晶界,采用NEB计算此时的跃迁能垒。
S5:考察晶界附近空位团簇对晶界跃迁过程的影响。将较小的模型体系周期性复制。在晶界附近典型位置放置典型尺寸的空位团簇,然后以平移向量VTk平移其中一个晶粒,采用NEB方法计算此时的跃迁能垒。
S6:考察辐照缺陷对VTk的影响。辐照缺陷不仅仅影响晶界跃迁能垒,还可能影响其最易跃迁路径。当晶界含有偏聚的辐照缺陷或晶界附近存在活性低的空位团簇或者空洞时,采用步骤S1-S3中类似方法重新确定晶界容易跃迁的向量和能垒。
S7:建立包含晶界运动的速率表格。
S8:辐照缺陷与晶界协同演化的的OKMC模拟方法。将晶界运动看作一个OKMC事件。经过上述步骤,得到纯晶界以及含有辐照缺陷的晶界跃迁向量及能垒,进而得到相应的速率。然后在OKMC循环中,调用相应的速率,选择、执行事件,并对可能出现的新的缺陷-晶界作用重新参数化相关过程。
参照图4,本发明算法流程如下:
(1)初始化模型参数,包括模型尺寸、体系中缺陷类型及缺陷数。设置模型边界条件,如周期性边界,自由边界条件等。设置模拟体系温度,读入缺陷构型文件。空位、自间隙的运动向量。
(2)读入点缺陷及其团簇扩散、复合、溶解能垒数据,以及作用范围数据,计算相应过程速率。
(3)计算晶界能量与平移向量关系。
(4)计算晶界平移最易跃迁向量和跃迁能垒。
(5)计算当前缺陷缺陷体系跃迁速率。根据事件类型,调用不同的速率。
(6)采用轮盘赌的方式筛选要执行的事件:包括Vn/SIAn扩散、Vn-Vn和SIAn-SIAn结合、晶界运动、Vn-SIAn复合、Vn/SIAn溶解、SIAn旋转等。
(7)执行选择的事件。对于辐照缺陷,按照其运动向量移动空位或自间隙。对于晶界运动,调整晶界位置,同时整体移动缺陷位置。
(7)更新计时器和缺陷位置,输出跃迁时间,当前体系构型等量。
(8)根据缺陷分布判断是否需要对晶界运动参数重参数化。是则转至步骤(3)。
(9)判断体系演化时间或运行步数是否超过最大设定值,是则退出,否则继续进行上述步骤(5)。
以上显示和描述了本发明的基本原理、主要特征和优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (3)
1.一种模拟辐照缺陷与晶界协同演化的方法,其特征在于:包括以下步骤
S1、获得晶界能量与平移向量关系:建立初始晶界模型,设其在平行晶界两个方向的最小周期分别为px和py,将一个晶粒相对于另一个晶粒平移时,相应的平移步长分别为dx、dy,平移向量VT=(i*dx,j*dy),其每个方向的值在0-px、0-py之间变化,每一次平移之后,对整个晶界模型进行充分原子弛豫,完成所有的平移-弛豫之后,得到晶界能量E与平移向量VT的关系E=E(i*dx,j*dy);
S2、获得基态对应的平移向量集合与基态间跃迁能垒:对比每次平移-弛豫后的体系能量,得到能量最低的相应平移向量集合VT1,同时得到平行弛豫后相应的体系构型R集合,对R集合中任意两个构型Ri和Rj,比较其距离,若其距离大于跃迁值时,认为两个态之间发生了跃迁,并调用微动弹性带NEB计算跃迁能垒;
S3、分析最易跃迁模式:分析晶界跃迁能垒与跃迁距离关系,获得跃迁能垒最低的跃迁及其对应的平移向量VTk;
S4、考察点缺陷对晶界跃迁过程的影响:在晶界内部放置一定浓度的空位或自间隙原子,按照平移向量VTk平移晶界模型,采用微动弹性带NEB计算当前条件下的跃迁能垒;
S5、考察晶界附近空位团簇对晶界跃迁过程的影响:将模型体系周期性复制,在晶界附近典型位置放置典型尺寸的空位团簇,然后以平移向量VTk平移其中一个晶粒,采用微动弹性带NEB方法计算当前条件下的跃迁能垒;
S6、考察辐照缺陷对VTk的影响:当晶界含有偏聚的辐照缺陷或晶界附近存在活性低的空位团簇或者空洞时,采用步骤S1-S3中的方法重新确定晶界容易跃迁的向量和能垒;
S7、建立包含晶界运动的速率表格;
S8、进行辐照缺陷与晶界协同演化的实体动力学蒙特卡洛OKMC模拟:将晶界运动看作一个实体动力学蒙特卡洛OKMC事件,经过步骤S1-S3,得到纯晶界以及含有辐照缺陷的晶界跃迁向量及能垒,进而得到相应的速率,然后在实体动力学蒙特卡洛OKMC循环中,调用相应的速率,选择、执行事件,并对可能出现的新的缺陷-晶界作用重新参数化相关过程。
2.根据权利要求1所述的一种模拟辐照缺陷与晶界协同演化的方法,其特征在于:步骤S1中所述建立初始晶界模型,其包括设置模型尺寸、体系中缺陷类型和缺陷数、模型边界条件、模拟体系温度、空位和自间隙的运动向量。
3.根据权利要求1所述的一种模拟辐照缺陷与晶界协同演化的方法,其特征在于:步骤S2中所述的跃迁值为2.0埃。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711352182.9A CN108052745B (zh) | 2017-12-15 | 2017-12-15 | 一种模拟辐照缺陷与晶界协同演化的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711352182.9A CN108052745B (zh) | 2017-12-15 | 2017-12-15 | 一种模拟辐照缺陷与晶界协同演化的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108052745A CN108052745A (zh) | 2018-05-18 |
CN108052745B true CN108052745B (zh) | 2021-04-06 |
Family
ID=62133250
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711352182.9A Active CN108052745B (zh) | 2017-12-15 | 2017-12-15 | 一种模拟辐照缺陷与晶界协同演化的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108052745B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108959709B (zh) * | 2018-06-04 | 2022-05-03 | 中国科学院合肥物质科学研究院 | 基于缺陷性质和多尺度模拟的晶界结构搜索方法 |
CN109727647B (zh) * | 2018-12-29 | 2020-11-24 | 中国原子能科学研究院 | 一种辐照导致材料结构演化的模拟方法和系统 |
CN109920487B (zh) * | 2019-01-17 | 2021-04-02 | 中国科学院合肥物质科学研究院 | 一种辐照环境下钨中氢滞留与脱附的多尺度模拟方法 |
CN110120248B (zh) * | 2019-04-08 | 2020-12-25 | 中国科学院合肥物质科学研究院 | 模拟纳米晶金属累积离位损伤的方法 |
CN110660453B (zh) * | 2019-10-09 | 2023-03-07 | 中国原子能科学研究院 | 基于指数时间差分格式求解速率理论方程的并行计算方法 |
CN111814315B (zh) * | 2020-06-17 | 2024-02-20 | 中国科学院合肥物质科学研究院 | 一种金属材料中缺陷团簇动力学性质的计算方法 |
CN113223641B (zh) * | 2021-04-28 | 2024-02-09 | 中国科学院合肥物质科学研究院 | 一种空位团簇跳出超势能谷的模拟方法 |
CN114112892B (zh) * | 2021-11-05 | 2024-05-24 | 南通爱尔思轻合金精密成型有限公司 | 一种模拟晶界迁移的实验方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101685530B (zh) * | 2008-09-23 | 2011-12-14 | 中国科学院过程工程研究所 | 一种在gpu上利用多体作用模型进行粒子计算的方法 |
CN105787227B (zh) * | 2016-05-11 | 2018-10-09 | 中国科学院近代物理研究所 | 结构材料辐照损伤的多gpu分子动力学模拟方法 |
CN107195345B (zh) * | 2017-06-09 | 2018-12-18 | 中国科学院合肥物质科学研究院 | 全能谱中子辐照损伤精确模拟系统及其算法 |
-
2017
- 2017-12-15 CN CN201711352182.9A patent/CN108052745B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108052745A (zh) | 2018-05-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108052745B (zh) | 一种模拟辐照缺陷与晶界协同演化的方法 | |
Brax et al. | Nonlinear structure formation with the environmentally dependent dilaton | |
Huber et al. | Ab initio modelling of solute segregation energies to a general grain boundary | |
Mattsson | Dark energy as a mirage | |
CN108090279B (zh) | 基于杂化lkmc和okmc模拟纳米结构核材料辐照损伤的方法 | |
Zhang et al. | Relative motion across the eastern Tibetan plateau: Contributions from faulting, internal strain and rotation rates | |
Niu et al. | Radial basis function approach in nuclear mass predictions | |
Hayward et al. | Energetics of small hydrogen–vacancy clusters in bcc iron | |
Govers et al. | Molecular dynamics simulation of helium and oxygen diffusion in UO2±x | |
CN111814315B (zh) | 一种金属材料中缺陷团簇动力学性质的计算方法 | |
Goudarzi et al. | Effects of three-body forces on the maximum mass of neutron stars in the lowest-order constrained variational formalism | |
Cusentino et al. | A comparison of interatomic potentials for modeling tungsten–hydrogen–helium plasma–surface interactions | |
Liang et al. | Three-dimensional phase-field simulations of intragranular gas bubble evolution in irradiated U-Mo fuel | |
Ikeno et al. | η-nucleus interaction from the d+ d reaction around the η production threshold | |
An et al. | Phase transition during inflation and the gravitational wave signal at pulsar timing arrays | |
Ashraf et al. | Exploring the effect of compiler optimizations on the reliability of HPC applications | |
Wolf et al. | Collapse and revival oscillations as a probe for the tunneling amplitude in an ultracold Bose gas | |
Bianchini et al. | Enabling QM-accurate simulation of dislocation motion in γ-Ni and α-Fe using a hybrid multiscale approach | |
CN110120248B (zh) | 模拟纳米晶金属累积离位损伤的方法 | |
Mareček et al. | Faster and lower-dose X-ray reflectivity measurements enabled by physics-informed modeling and artificial intelligence co-refinement | |
Aldinucci et al. | StochKit-FF: Efficient systems biology on multicore architectures | |
Becquart et al. | Effect of displacement cascade structure and defect mobility on the growth of point defect clusters under irradiation | |
Mikelsons et al. | Relationship between Hirsch-Fye and weak-coupling diagrammatic quantum Monte Carlo methods | |
Sanders et al. | Self-consistent triaxial models | |
Germann et al. | Tflop/s multibillion-atom molecular dynamics simulations and visualization/analysis on BlueGene/L |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |