CN111048158A - 基于spparks软件的kmc方法uo2势能力场 - Google Patents
基于spparks软件的kmc方法uo2势能力场 Download PDFInfo
- Publication number
- CN111048158A CN111048158A CN201911362612.4A CN201911362612A CN111048158A CN 111048158 A CN111048158 A CN 111048158A CN 201911362612 A CN201911362612 A CN 201911362612A CN 111048158 A CN111048158 A CN 111048158A
- Authority
- CN
- China
- Prior art keywords
- potential energy
- formula
- spparks
- pair
- field
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 16
- 238000005381 potential energy Methods 0.000 claims abstract description 43
- 238000004088 simulation Methods 0.000 claims abstract description 7
- 238000000329 molecular dynamics simulation Methods 0.000 claims abstract description 5
- 238000004364 calculation method Methods 0.000 claims description 11
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 9
- 229910052760 oxygen Inorganic materials 0.000 claims description 9
- 239000001301 oxygen Substances 0.000 claims description 9
- 229910052770 Uranium Inorganic materials 0.000 claims description 4
- 150000001224 Uranium Chemical group 0.000 claims 1
- 230000010354 integration Effects 0.000 claims 1
- FCTBKIHDJGHPPO-UHFFFAOYSA-N uranium dioxide Inorganic materials O=[U]=O FCTBKIHDJGHPPO-UHFFFAOYSA-N 0.000 abstract description 11
- 238000000342 Monte Carlo simulation Methods 0.000 abstract description 2
- OOAWCECZEHPMBX-UHFFFAOYSA-N oxygen(2-);uranium(4+) Chemical compound [O-2].[O-2].[U+4] OOAWCECZEHPMBX-UHFFFAOYSA-N 0.000 abstract description 2
- 239000011159 matrix material Substances 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical group [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 description 2
- WUKWITHWXAAZEY-UHFFFAOYSA-L calcium difluoride Chemical compound [F-].[F-].[Ca+2] WUKWITHWXAAZEY-UHFFFAOYSA-L 0.000 description 1
- 239000000919 ceramic Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 239000010436 fluorite Substances 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 229910001092 metal group alloy Inorganic materials 0.000 description 1
- 238000006386 neutralization reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
Landscapes
- Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本人公开了一种基于SPPARKS软件的UO2（二氧化铀）势能力场。该势能力场可用于SPPARKS软件中KMC（动态蒙特卡洛方法）的演化模拟的过程，进行能量的统计，能够进行比SPPARKS原有的Lennard‑Jones（LJ）势能力场更加精确，并进行了气泡穿过UO2的模拟演化过程进行了测试，结果能够与传统分子动力学MD方法吻合。
Description
技术领域
本人公开了一种基于SPPARKS软件的UO2（二氧化铀）势能力场。该势能力场可用于SPPARKS软件中KMC（动态蒙特卡洛方法）的演化模拟的过程，进行能量的统计，能够进行比SPPARKS原有的Lennard-Jones（LJ）势能力场更加精确，并进行了气泡穿过UO2的模拟演化过程进行了测试，结果能够与传统分子动力学MD方法吻合。
背景技术
在OKMC的模拟中，采用了SPPARKS软件，SPPARKS软件是圣地亚哥课题组开发的一款基于Monte Carlo的开源代码，主要具备两种算法拒绝式动态蒙特卡洛（rKMC）和大都会蒙特卡洛（MMC），能够进行点阵中和离散点阵计算，主要用以缺陷点阵和杂质迁移的演化过程模拟。在SPPARKS实现迁移和演化的过程中，采用Lennard-Jones（LJ）势能力场，但是这种势能力场对于金属合金来说，不能够精确的描述原子间的作用力和进行准确的能量统计，所以提出通过将势能力场修改为经过MD验证的Basak势能力场代替LJ力场实现更为精确的能量计算。其中，模拟过程采用的研究对象是UO2，它是典型陶瓷型氧化物材料，属于萤石型（CF2）结构，立方晶系，空间群Fm3m，以此为例进行新型势能力场的检验。
发明内容
本发明主要目的在于提供一种UO2在使用SPPARKS软甲中模拟时的势能力场，克服现有势能力场不够精确的问题。
为了达到上述目的，本文提供了以下的技术方案：
本发明实例提供了一种新的势能力场，用于UO2在SPPARKS软件中基于KMC的模拟。
通过C.B.Basak在用传统分子动力学进行UO2的热力学性能的预测中，提出了热力学上的Basak公式中热力学参数的具体值，通过对公式以及数值的整合，确定了三种原子对（别是U-U（铀原子之间）对，O-O（氧与氧之间）对，U-O（铀原子与氧之间）对）的势能力场公式，根据势能参数化简公式，可以得到U-U对公式为：
O-O对公式为：
U-O对公式为：
为了保证整个势能计算部分流畅运行，通过C语言将势能计算公式进行了代码编写，并运行了伪代码，结果显示正常，可以运行，伪代码图1。
在伪代码运行成功后，进行了SPPARKS软件中的LJ势能力场代码中修改势能写法，和计算公式，通过C语言进行势能力场编写，同时将主要参数进行了带入。主要计算部分可展示为图2所示。
势能力场建立之后，对新的势能力场进行了测试，测试的方法采用SPPARKS软件定义气泡穿过UO2之后，其晶粒表现，将主要的几个阶段进行了可视化，其过程如图3所示。
附图说明
图1 UO2势能计算公式伪代码。
图2 SPPARKS软件中的修改得到的势能力场代码。
图3势能力场检测的模型演化。
具体实施方式
通过C.B.Basak在用传统分子动力学进行UO2的热力学性能的预测中，提出了热力学上的Basak公式中热力学参数的具体值，通过对公式以及数值的整合，确定了三种原子对（别是U-U（铀原子之间）对，O-O（氧与氧之间）对，U-O（铀原子与氧之间）对）的势能力场公式，根据势能参数化简公式，可以得到U-U对公式为：
O-O对公式为：
U-O对公式为：
为了保证整个势能计算部分流畅运行，通过C语言将势能计算公式进行了代码编写，并运行了伪代码，结果显示正常，可以运行，伪代码图1。
在伪代码运行成功后，进行了SPPARKS软件中的LJ势能力场代码中修改势能写法，和计算公式，通过C语言进行势能力场编写，同时将主要参数进行了带入。主要计算部分可展示为图2所示。
势能力场建立之后，对新的势能力场进行了测试，测试的方法采用SPPARKS软件定义气泡穿过UO2之后，其晶粒表现，将主要的几个阶段进行了可视化，其过程如图3所示。
Claims (1)
1.基于SPPARKS软件的KMC方法UO2势能力场其特征在于，本发明实例提供了一种新的势能力场，用于UO2在SPPARKS软件中基于KMC的模拟，通过C.B.Basak在用传统分子动力学进行UO2的热力学性能的预测中，提出了热力学上的Basak公式中热力学参数的具体值，通过对公式以及数值的整合，确定了三种原子对（别是U-U（铀原子之间）对，O-O（氧与氧之间）对，U-O（铀原子与氧之间）对）的势能力场公式，根据势能参数化简公式，可以得到U-U对公式为：
O-O对公式为：
U-O对公式为：
为了保证整个势能计算部分流畅运行，通过C语言将势能计算公式进行了代码编写，并运行了伪代码。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911362612.4A CN111048158A (zh) | 2019-12-26 | 2019-12-26 | 基于spparks软件的kmc方法uo2势能力场 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911362612.4A CN111048158A (zh) | 2019-12-26 | 2019-12-26 | 基于spparks软件的kmc方法uo2势能力场 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111048158A true CN111048158A (zh) | 2020-04-21 |
Family
ID=70239987
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911362612.4A Pending CN111048158A (zh) | 2019-12-26 | 2019-12-26 | 基于spparks软件的kmc方法uo2势能力场 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111048158A (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102317764A (zh) * | 2008-12-24 | 2012-01-11 | 阿海珐核能公司 | 基本为球形的多层结构的非破坏性和无接触的表征方法及相关装置 |
| US20120014491A1 (en) * | 2009-07-13 | 2012-01-19 | Mike Deeth | Nuclear fusion power plant having a liquid reactor core of molten glass that is made laseractive and functions as a tritium breeding blanket which is capable of acousticly compressing/confining fuel so that it radiates and triggers outgoing laser cascades that will reflect from the blast chamber's spherical inside wall and return like photonic Tsunamis, crushing, heating, and causing thermonuclear ignition of the fuel so that heat engines and piezoelectric harvesters can convert the released energy into electricity |
| CN107423570A (zh) * | 2017-08-02 | 2017-12-01 | 南昌立德生物技术有限公司 | 快速准确计算蛋白酶与药物分子之间亲和自由能的算法 |
| CN109086497A (zh) * | 2018-07-16 | 2018-12-25 | 中国科学院宁波材料技术与工程研究所 | 基于粒子群算法的金属与合金势能力场开发方法 |
| CN110097927A (zh) * | 2019-05-10 | 2019-08-06 | 青岛理工大学 | 基于分子动力学测试电场作用下离子扩散系数的方法 |
-
2019
- 2019-12-26 CN CN201911362612.4A patent/CN111048158A/zh active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102317764A (zh) * | 2008-12-24 | 2012-01-11 | 阿海珐核能公司 | 基本为球形的多层结构的非破坏性和无接触的表征方法及相关装置 |
| US20120014491A1 (en) * | 2009-07-13 | 2012-01-19 | Mike Deeth | Nuclear fusion power plant having a liquid reactor core of molten glass that is made laseractive and functions as a tritium breeding blanket which is capable of acousticly compressing/confining fuel so that it radiates and triggers outgoing laser cascades that will reflect from the blast chamber's spherical inside wall and return like photonic Tsunamis, crushing, heating, and causing thermonuclear ignition of the fuel so that heat engines and piezoelectric harvesters can convert the released energy into electricity |
| CN107423570A (zh) * | 2017-08-02 | 2017-12-01 | 南昌立德生物技术有限公司 | 快速准确计算蛋白酶与药物分子之间亲和自由能的算法 |
| CN109086497A (zh) * | 2018-07-16 | 2018-12-25 | 中国科学院宁波材料技术与工程研究所 | 基于粒子群算法的金属与合金势能力场开发方法 |
| CN110097927A (zh) * | 2019-05-10 | 2019-08-06 | 青岛理工大学 | 基于分子动力学测试电场作用下离子扩散系数的方法 |
Non-Patent Citations (1)
| Title |
|---|
| 王艺博: "U3Si5氦损伤的计算与KMC方法下UO2势能力场的初步设计", 中国优秀硕士论文电子期刊网 基础科学辑, no. 02, pages 33 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110631908A (zh) | 一种岩石单轴压缩全过程蠕变损伤模型的构建方法 | |
| CN104615840B (zh) | 一种数字仿真模型的修正方法及系统 | |
| CN112785034B (zh) | 基于融合神经网络的台风路径预报方法、系统、介质及终端 | |
| CN103353866A (zh) | 一种支持xna技术的三维模型文件格式转换方法 | |
| CN119091085B (zh) | 一种面向三维飞行器气动系数预测的数据增广与模型训练方法、系统及介质 | |
| CN106777691A (zh) | 用于结构动力学仿真的橡胶o形圈有限元建模方法 | |
| CN114974465A (zh) | 基于迁移学习算法的材料流动应力预测方法及系统 | |
| CN117648829A (zh) | 非线性海工结构动力响应的拉普拉斯域计算方法 | |
| CN113742948B (zh) | 一种超高强度抽油杆p-s-n曲线拟合新模型及方法 | |
| CN111048158A (zh) | 基于spparks软件的kmc方法uo2势能力场 | |
| Liu et al. | Automatic generation of personalized Chinese handwriting characters | |
| CN117252063A (zh) | 基于机器学习的岩质高边坡开挖锚固变形预测方法及系统 | |
| CN106682450A (zh) | 一种基于状态转移模型的新一代测序拷贝数变异仿真方法 | |
| CN113705038A (zh) | 一种利用屈曲特征值判断不同模态转变的方法及设备 | |
| Zhu et al. | How robust is a large pre-trained language model for code generationƒ a case on attacking gpt2 | |
| CN106383946A (zh) | 一种加速退化模型参数估计方法 | |
| CN114491877A (zh) | 管道泄漏影响区域的确定方法、装置、系统、终端及介质 | |
| Huang et al. | Predicting brittle fracture surface shape from a versatile database | |
| CN120745353B (zh) | 基于原子尺度动力学获取含能材料爆速与爆压的方法 | |
| Broström et al. | Generalised Linear Models with Clustered Data: Fixed and random effects models with glmmML | |
| CN116486946A (zh) | 一种基于Matlab的力学本构模型构建方法 | |
| CN120181263B (zh) | 一种基于选择性干预的语言模型知识蒸馏方法及系统 | |
| CN113886784B (zh) | 基于语料库提升小训练集猜测效率的口令猜测方法 | |
| CN116631541A (zh) | 一种压缩条件下提高非晶合金塑性的模拟方法 | |
| CN115659498A (zh) | 一种Nastran计算文件的梁架尺寸缩比方法 |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20200421 |
|
| WD01 | Invention patent application deemed withdrawn after publication |