CN111899819B - 一种辐照缺陷和纳米相协同演化的相场模拟方法 - Google Patents
一种辐照缺陷和纳米相协同演化的相场模拟方法 Download PDFInfo
- Publication number
- CN111899819B CN111899819B CN202010599137.9A CN202010599137A CN111899819B CN 111899819 B CN111899819 B CN 111899819B CN 202010599137 A CN202010599137 A CN 202010599137A CN 111899819 B CN111899819 B CN 111899819B
- Authority
- CN
- China
- Prior art keywords
- defect
- nanophase
- irradiation
- evolution
- defects
- 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
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Computing Systems (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Geometry (AREA)
- Operations Research (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明为一种辐照缺陷和纳米相协同演化的相场模拟方法。通过合金热力学自由能来表征系统的演变;建立辐照缺陷–合金体系的总自由能;创建辐照缺陷和纳米相共同演化的相场方程,实现辐照下合金缺陷和纳米相演化的模拟与表征;分析辐照温度和缺陷浓度变化条件下的缺陷和纳米相微观组织形貌的变化,得出辐照下缺陷和纳米析出相的演变规律;同时分析缺陷和纳米相的体积分数变化,得出不同辐照温度和缺陷浓度对缺陷和纳米相的形成及演化动力学的影响规律;本发明提供了一种辐照缺陷和纳米相协同演化的相场模拟方法,能够获得辐照下缺陷和相分解之间的相互作用关系,为辐照下合金材料的设计和服役安全提供指导。
Description
技术领域
本发明涉及金属材料辐照下微观组织控制技术领域,具体为一种辐照缺陷和纳米相协同演化的相场模拟方法。
背景技术
在核领域,核反应堆结构材料在高温和辐照下,合金中的原子将会离开原来晶格位置而形成空位和间隙原子,缺陷的存在促进了纳米相的析出,阻碍了位错运动引起合金脆硬性的增加,导致材料的使用寿命降低。因此研究辐照下缺陷和纳米相的演化对合金的安全使用尤为重要。
从目前的研究中可以看出辐照对合金相转变的影响因素主要有辐照温度和辐照剂量,辐照温度的升高使缺陷的迁移作用加强,原子扩散作用增强,继而促进了纳米相转变。辐照剂量的增加促使级联作用增强,缺陷数量的增加为原子扩散提供通道,加速原子扩散,使相分解加快。由此可知,辐照剂量和辐照温度都是通过加强缺陷作用来促进相分解,但是对于辐照下缺陷对纳米相演化的动力学影响规律确还不清楚。
目前较多的模拟研究都只针对辐照下缺陷或者纳米相单个的演变,对于辐照下缺陷和纳米相的协同演化的研究目前还很缺少,在相场法模拟中可以将缺陷简化成同等的元素变量,来模拟缺陷和纳米相的协同演变。同时利用相场法研究还可以弥补实验周期长,耗费成本大,研究对象单一和不易实时观测的不足。因此,利用相场模拟手段研究辐照下缺陷和纳米相协同演化是非常有优势且有必要的。
发明内容
本发明所解决的技术问题在于提供一种辐照缺陷和纳米相协同演化的相场模拟方法,解决现有研究方法无法同时研究缺陷和纳米相的问题。
实现本发明目的提供技术方案如下:
一种辐照缺陷和纳米相协同演化的相场模拟方法:
(1)利用合金热力学自由能进行缺陷和纳米相系统的热力学描述;
(2)根据热力学描述,求解辐照下缺陷和纳米相系统的化学自由能和弹性能;
(3)基于Chan-Hilliard扩散方程,建立辐照下缺陷和纳米相演变的相场动力学方程;
(4)根据输入变量来模拟不同辐照条件下缺陷和纳米相的组织形貌以及成分演变,分析不同辐照条件下缺陷对相分解动力学规律的影响。
具体内容如下:
根据合金热力学特征对(1)中的热力学进行描述,主要包括:系统的化学自由能、弹性能、梯度能,具体如下:
其中F为系统总的自由能,fchem为合金的化学自由能,Vm为合金的摩尔体积,i代表合金元素,ci为元素i的原子分数,cV和cI分为空位和间隙原子的原子分数,kj(j=i,V,I)为梯度能系数,felas为弹性能密度。
在本发明中将缺陷的形成能及其浓度耦合到系统的化学自由能中,将(2)中的化学自由能表示为:
在本发明中将缺陷和应变所引起成分的不均匀耦合到本征应变中,(2)中的弹性能密度可表示如下:
其中ε0,εV0,和εI0分别为成分不均匀以及空位和间隙原子引入产生的晶格膨胀系数,弹性应变εkl与位移有关,其表达式为:
通过将缺陷有关量耦合到化学自由能和弹性能中,便于利用相场法来描述缺陷的演变及其对相分解的影响。
对于本发明的步骤(3)中所描述的辐照缺陷和纳米相协同演变的动力学方程可表示为:
其中M为合金的化学迁移率,Mj(j=V,I)为空位和间隙原子迁移率,ξj(j=i,V,I)为热力学起伏,Pj(j=V,I)和Sj(j=V,I)分别为缺陷的产生速率和湮灭速率,RVI为缺陷的合并速率。在相场模拟计算中利用半隐式傅里叶谱算法对动力学演变方程进行求解,得出成分与时间和空间位置关系,通过定量的分析来完成(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为系统总的自由能,fchem为合金的化学自由能,Vm为合金的摩尔体积,i代表合金元素,ci为元素i的原子分数,cV和cI分为空位和间隙原子的原子分数,kj(j=i,V,I)为梯度能系数,felas为弹性能密度。
将缺陷的形成能及其浓度耦合到系统的化学自由能中,步骤(2)中Fe-Cr合金的化学自由能可表示为:
(i=Fe,Cr;k=Fe,Cr;i≠k)
将缺陷和应变所引起成分的不均匀耦合到本征应变中,步骤(2)中Fe-Cr合金的弹性能密度可表示如下:
其中ε0,εV0,和εI0分别为成分不均匀以及空位和间隙原子引入产生的晶格膨胀系数,弹性应变εkl与位移有关,其表达式为:
步骤(3)中,Fe-Cr合金辐照缺陷和相协同演变的动力学方程可表示为:
其中M为Fe-Cr合金的化学迁移率可表示为M=[cCrMFe+(1-cCr)MCr]cCr(1-cCr),Mj(j=V,I)为空位和间隙原子迁移率,从动力学方程中可知缺陷和原子可通过迁移率来建立内在联系,而空位扩散机制表明原子通过与空位交换来完成其迁移,因此该相场方程在理论上是符合的。ξj(j=i,V,I)为热力学起伏,Pj(j=V,I)和Sj(j=V,I)分别为缺陷的产生速率和湮灭速率,RVI为缺陷的合并速率,速率的引入使得相场法更加符合实际。在模拟计算过程中利用半隐式傅里叶谱算法对动力学演变方程进行求解。
步骤(4),设置Fe-Cr合金初始成分为cCr=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;从图中可以看出,随着缺陷浓度的升高,α′相体积分数开始变化的时间提前,相分解时间提前,相变受到促进。
总的来说,辐照温度和缺陷浓度都可以促进相分解,并且温度越高作用越强。这可作为辐照下核反应堆结构材料设计应用时的理论参考,并以此优化合金的组织性能。
以上所述是对本发明的优选实施案列及原理进行了详细的说明,对于本邻域技术人员依据本发明所提供的思想,在具体实施上的差别都将视为本发明的保护范畴。
Claims (3)
1.一种辐照缺陷和纳米相协同演化的相场模拟方法,其特征在于,包括以下步骤:
(1)利用合金热力学自由能进行缺陷和纳米相系统的热力学描述;
(2)根据热力学描述,求解辐照下缺陷和纳米相系统的化学自由能和弹性能;
(3)基于Chan-Hilliard扩散方程,建立辐照下缺陷和纳米相演变的相场动力学方程;
(4)根据输入变量来模拟不同辐照条件下缺陷和纳米相的组织形貌以及成分演变,分析不同辐照条件下缺陷对相分解动力学规律的影响;步骤(1)中的热力学描述包括:系统的化学自由能、弹性能、梯度能,具体如下:
其中F为系统总的自由能,fchem为合金的化学自由能,Vm为合金的摩尔体积,i代表合金元素,ci为元素i的原子分数,cV和cI分为空位和间隙原子的原子分数,kj,其中j=i,V,I为梯度能系数,felas为弹性能密度;将缺陷的形成能及其浓度耦合到系统的化学自由能中,步骤(2)中的化学自由能可表示为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010599137.9A CN111899819B (zh) | 2020-06-28 | 2020-06-28 | 一种辐照缺陷和纳米相协同演化的相场模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010599137.9A CN111899819B (zh) | 2020-06-28 | 2020-06-28 | 一种辐照缺陷和纳米相协同演化的相场模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111899819A CN111899819A (zh) | 2020-11-06 |
CN111899819B true CN111899819B (zh) | 2022-09-27 |
Family
ID=73207223
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010599137.9A Active CN111899819B (zh) | 2020-06-28 | 2020-06-28 | 一种辐照缺陷和纳米相协同演化的相场模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111899819B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112926205B (zh) * | 2021-02-24 | 2022-05-13 | 中国核动力研究设计院 | 基于团簇动力学的锆基合金辐照损伤模拟方法及模型系统 |
CN113255136B (zh) * | 2021-05-28 | 2024-05-10 | 大连理工大学 | 一种预测辐照气泡演化的相场模拟方法及系统 |
CN114330053B (zh) * | 2021-12-16 | 2022-12-02 | 北京应用物理与计算数学研究所 | 模拟金属表面氢化点腐蚀形貌演化的方法及装置 |
CN114970106B (zh) * | 2022-04-28 | 2023-03-21 | 天津大学 | 一种基于微观组织的辐照硬化预测方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108052794A (zh) * | 2018-01-16 | 2018-05-18 | 桂林电子科技大学 | 不同尺寸的纳米Ag-Cu合金体系的共晶温度的预测方法 |
CN110379470A (zh) * | 2019-04-10 | 2019-10-25 | 西北工业大学 | 单晶高温合金的组织和缺陷的模拟方法 |
CN111324969A (zh) * | 2020-03-06 | 2020-06-23 | 南京理工大学 | 一种温度梯度调控层状高温合金微观组织的相场模拟方法 |
-
2020
- 2020-06-28 CN CN202010599137.9A patent/CN111899819B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108052794A (zh) * | 2018-01-16 | 2018-05-18 | 桂林电子科技大学 | 不同尺寸的纳米Ag-Cu合金体系的共晶温度的预测方法 |
CN110379470A (zh) * | 2019-04-10 | 2019-10-25 | 西北工业大学 | 单晶高温合金的组织和缺陷的模拟方法 |
CN111324969A (zh) * | 2020-03-06 | 2020-06-23 | 南京理工大学 | 一种温度梯度调控层状高温合金微观组织的相场模拟方法 |
Non-Patent Citations (1)
Title |
---|
辐照及温度作用下二氧化铀中空洞的动力学相场模拟;彭玉艺;《中国优秀硕士论文全文数据库》;20190115;正文第12-30页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111899819A (zh) | 2020-11-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111899819B (zh) | 一种辐照缺陷和纳米相协同演化的相场模拟方法 | |
Pourbabak et al. | Microscopic investigation of as built and hot isostatic pressed Hastelloy X processed by Selective Laser Melting | |
Odette | Recent progress in developing and qualifying nanostructured ferritic alloys for advanced fission and fusion applications | |
Swenson et al. | Nanocluster irradiation evolution in Fe-9% Cr ODS and ferritic-martensitic alloys | |
Hémery et al. | Microplasticity at room temperature in α/β titanium alloys | |
Caballero et al. | New experimental evidence of the diffusionless transformation nature of bainite | |
Li et al. | Experimental verification of segregation of carbon and precipitation of carbides due to deep cryogenic treatment for tool steel by internal friction method | |
Wakai et al. | Effect of helium production on swelling of F82H irradiated in HFIR | |
Haghdadi et al. | Five-parameter characterization of intervariant boundaries in additively manufactured Ti-6Al-4V | |
Banerjee et al. | Kinetics of Order–Disorder Transformation in Alloys under Electron Irradiation | |
Gorsse et al. | Modeling the precipitation processes and the formation of hierarchical microstructures in a single crystal high entropy superalloy | |
Chai et al. | Microstructural characteristics of a commercially pure Zr treated by pulsed laser at different powers | |
CN103789530A (zh) | 一种用离子预辐照提高低活化材料抗辐照性能的方法 | |
Guo et al. | Electroshock treatment dependent microstructural evolution and mechanical properties of near-β titanium alloy manufactured by directed energy deposition | |
Wu et al. | High temperature creep mechanisms of a single crystal superalloy: A phase-field simulation and microstructure characterization | |
Watanabe et al. | Hardening and microstructural evolution of A533b steels irradiated with Fe ions and electrons | |
Sun et al. | Multi-component phase-field simulation of microstructural evolution and elemental distribution in Fe–Cu–Mn–Ni–Al alloy | |
Losertová et al. | Microstructure evolution of heat treated NiTi alloys | |
Wang et al. | Effect of Nb addition on the microstructure of Ni–12Mo–7Cr based superalloy | |
Kamboj et al. | The effect of phosphorus on precipitation in irradiated reactor pressure vessel (RPV) steels | |
Perrier et al. | Characterization and modeling of precipitation kinetics in a Fe-Si-Ti alloy | |
Penisten | The Mechanism of Radiation-Induced Segregation in Ferritic-Martensitic Steels. | |
Zhu et al. | Irradiation behaviour of α2 and γ phases in He ion implanted titanium aluminide alloy | |
Ye et al. | Positron annihilation study of open volume defects and Cr segregation in deformed CoCrFeMnNi alloy | |
Lee et al. | Observation and rate theory modeling of grain boundary segregation in Σ3 twin boundaries in ion-irradiated stainless steel 316 |
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 |