EAST面对等离子体材料的分子动力学模拟与实验研究课件.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《EAST面对等离子体材料的分子动力学模拟与实验研究课件.ppt》由用户(晟晟文业)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- EAST 面对 等离子体 材料 分子 动力学 模拟 实验 研究 课件
- 资源描述:
-
1、 ASIPP报告人:杨钟时报告人:杨钟时导导 师:罗广南师:罗广南 研究员研究员钨面对等离子体材料的分子动力学模拟中期报告(中期报告(2009-12-09)ASIPP2022-11-301主要内容主要内容u研究背景研究背景u研究方法研究方法u研究目标研究目标u论文进展情况论文进展情况u存在问题与改进措施存在问题与改进措施 ASIPP2022-11-302主要内容主要内容u研究背景研究背景u研究方法研究方法u研究目标研究目标u论文进展情况论文进展情况u存在问题与改进措施存在问题与改进措施 ASIPP2022-11-303人类能源危机与聚变能研究进展人类能源危机与聚变能研究进展能源需求日益增大能源
2、需求日益增大u传统能源:传统能源:化石燃料化石燃料,存量有限,并带来严峻环境问题,存量有限,并带来严峻环境问题u替代能源替代能源:太阳能、生物能、风能太阳能、生物能、风能等产量与效率很低等产量与效率很低 裂变能裂变能:铀储量有限,核废料造成的环境和安全问题:铀储量有限,核废料造成的环境和安全问题 聚变能聚变能:受控核聚变将为人类提供终极理想清洁能源:受控核聚变将为人类提供终极理想清洁能源ITER cutaway 磁约束等离子体与托卡马克热核聚变装置磁约束等离子体与托卡马克热核聚变装置 ASIPP2022-11-304聚变能领域聚变能领域W W材料的应用材料的应用钨(钨(W W)材料优良特性)材
3、料优良特性高熔点高熔点低物理溅射率和高溅射阈能低物理溅射率和高溅射阈能不与氢发生化学刻蚀以及低的不与氢发生化学刻蚀以及低的H H滞留滞留方便用喷涂方法进行大规模生产方便用喷涂方法进行大规模生产钨(钨(W W)材料在聚变装置中的应用)材料在聚变装置中的应用ITERITER:第一阶段:第一阶段W W用作偏滤器非靶板的用作偏滤器非靶板的PFMPFM,最后实现全钨第一壁。,最后实现全钨第一壁。ITERITER后的堆型设计:全后的堆型设计:全W W的的PFMPFM成为共识成为共识EASTEAST:在:在3-53-5年实现主动冷却年实现主动冷却VPS-W/Cu-VPS-W/Cu-偏滤器偏滤器德国德国ASD
4、EX-UASDEX-U:已实现全:已实现全W W第一壁第一壁英国英国JETJET托卡马克:托卡马克:ITER-like WallITER-like Wall,W W块材作为偏滤器的块材作为偏滤器的PFMPFM日本九州大学日本九州大学QUESTQUEST球马克:将球马克:将PFMPFM升级为升级为W W涂层材料涂层材料 ASIPP2022-11-305主要内容主要内容u研究背景研究背景u研究方法研究方法u研究目标研究目标u论文进展情况论文进展情况u存在问题与改进存在问题与改进 ASIPP2022-11-306分子动力学在聚变材料中的应用分子动力学在聚变材料中的应用材料计算模拟:材料计算模拟:与理
5、论和实验相并列的一个新的分支与理论和实验相并列的一个新的分支分子动力学分子动力学(Molecular Dynamics-MD):l原理:根据量子力学等基本物理理论,建立物理模型,构造势函数,在原理:根据量子力学等基本物理理论,建立物理模型,构造势函数,在微观的分子和原子层面上对材料进行研究,预测材料服役行为微观的分子和原子层面上对材料进行研究,预测材料服役行为l方法:建立一组粒子的运动方程,通过数值计算得到每个粒子在各时刻方法:建立一组粒子的运动方程,通过数值计算得到每个粒子在各时刻的坐标与动量,求得相空间的运动轨迹,并用统计方法得到系统的静态的坐标与动量,求得相空间的运动轨迹,并用统计方法得
6、到系统的静态和动态特性,从而得到系统的宏观特性和动态特性,从而得到系统的宏观特性核聚变研究的极端条件与特殊性:核聚变研究的极端条件与特殊性:超高温、强辐照、强磁场等极端条件超高温、强辐照、强磁场等极端条件下,等离子体与下,等离子体与PFMPFM的相互作用涉及到微观尺度和很短的时间尺度的相互作用涉及到微观尺度和很短的时间尺度分子动力学的应用:分子动力学的应用:材料表面溅射、高能粒子的射程分布、材料的内部材料表面溅射、高能粒子的射程分布、材料的内部辐照损伤,材料的熔化以及相变等辐照损伤,材料的熔化以及相变等 ASIPP2022-11-307BOP势函数势函数(Bond Order Potentia
7、ls)BOP势:势:表述具有不同化学键和表述具有不同化学键和不同结构材料中的原子间相互作不同结构材料中的原子间相互作用进行用进行Brenner势:势:用于不同的用于不同的C C或者或者CHCH系统,可描述纯系统,可描述纯C C结构材料和结构材料和CHCH分子以及键合和键断等动力学分子以及键合和键断等动力学效应效应ZjattrepirbVrqVE1)()(21Z:最近邻数最近邻数q:每个原子的价电子数每个原子的价电子数b:健级健级(bond order)Vrep(r):紧邻原子间排斥力紧邻原子间排斥力Vatt(r):引力函数引力函数 修订修订Brenner势:势:JuslinJuslin等对的等
8、对的CHCH势进行改造,并重新计算拟合了势进行改造,并重新计算拟合了W-WW-W,W-CW-C和和W-HW-H之间的参数设置,有效模拟之间的参数设置,有效模拟H H和和CHCH与金属与金属 W W之间的相互作之间的相互作用用 用于包含多元材料的不同性质,用于包含多元材料的不同性质,如表面,缺陷,融化性质等如表面,缺陷,融化性质等 ASIPP2022-11-308主要内容主要内容u研究背景研究背景u研究方法研究方法u研究目标研究目标u论文进展情况论文进展情况u存在问题与改进措施存在问题与改进措施 ASIPP论文的计划目标论文的计划目标主要研究内容主要研究内容:l 采用采用MDMD方法对等离子体与
9、方法对等离子体与W-PFMW-PFM相互作用过程进行计算模拟研究相互作用过程进行计算模拟研究 低能低能H H粒子在粒子在W W表面的反射和吸附,表面的反射和吸附,H H在在W W体内的俘获等体内的俘获等 C C粒子在粒子在W W表面的反射和吸附,体内的射程分布、沟道效应以及扩散行为等表面的反射和吸附,体内的射程分布、沟道效应以及扩散行为等达到的目标达到的目标:运用运用MDMD方法,研究方法,研究H H、C C和和W W材料相互作用的物理机制材料相互作用的物理机制 W W在等离子体辐照下材料本身结构与性能的变化规律及其机制在等离子体辐照下材料本身结构与性能的变化规律及其机制 了解了解W W材料在
10、材料在EASTEAST和和ITERITER中应用的物理基础,为聚变中应用的物理基础,为聚变PFMPFM的设计和应用提供的设计和应用提供参考依据参考依据2022-11-309 ASIPP2022-11-3010主要内容主要内容u研究背景研究背景u研究方法研究方法u研究目标研究目标u论文进展情况论文进展情况u存在问题与改进措施存在问题与改进措施 ASIPP非晶非晶W的的MD模拟模拟2022-11-301124681005101520012345 g(r)r()W-W(crystalline W)(b)g(r)W-W(amorphous W)(a)晶态和非晶态晶态和非晶态W中的中的 W-W 对分布函
11、数对分布函数g(r)非晶胞模拟过程:初始结构:初始结构:bccbcc晶态晶态W W元胞,密度为元胞,密度为19.25 19.25 g/cmg/cm3 3,三个方向施加周期性边界条件三个方向施加周期性边界条件 系统在系统在300K300K室温下平衡室温下平衡1010个个psps后,升温至后,升温至4000K4000K 在在4000K4000K融化状态下保持融化状态下保持200ps200ps 将元胞冷却至将元胞冷却至300K300K,冷却速率为,冷却速率为40K.ps40K.ps-1-1 在在300K300K下平衡下平衡50ps50psl 第一紧邻位置和数值相似揭示非晶第一紧邻位置和数值相似揭示非
12、晶W W的的的短程有序性。的短程有序性。l 第二紧邻位置偏移,更高级的近邻峰值第二紧邻位置偏移,更高级的近邻峰值拓展或消失揭示长程无序。拓展或消失揭示长程无序。ASIPP模拟元胞初始化模拟元胞初始化2022-11-3012 表面模拟:X X和和Y Y方向施加周期性边界条方向施加周期性边界条件。最底三层原子保持固定模拟为衬底。件。最底三层原子保持固定模拟为衬底。温控:衬底以上已经四个侧面的三个原衬底以上已经四个侧面的三个原子层内施加温控条件。子层内施加温控条件。初始位置:入射原子在元胞外势函数截入射原子在元胞外势函数截断距离之外,非温度控制区域以上作为初断距离之外,非温度控制区域以上作为初始位置
13、。随机改变原子的入射位置,进行始位置。随机改变原子的入射位置,进行统计。统计。初始动能:改变原子的入射动能,研究改变原子的入射动能,研究能量依赖关系。入射粒子的能量分为两种,能量依赖关系。入射粒子的能量分为两种,一种是固定动能,另外一种就是载能粒子一种是固定动能,另外一种就是载能粒子的能量服从麦克斯韦分。的能量服从麦克斯韦分。初始角度:改变原子入射角度,研究角改变原子入射角度,研究角度依赖关系。度依赖关系。元胞温度:改变元胞温度,研究衬底温改变元胞温度,研究衬底温度依赖关系。选择模拟的步长。度依赖关系。选择模拟的步长。两种辐照模拟:非累计型模拟和累计型非累计型模拟和累计型模拟。模拟。ASIPP
14、表面损伤表面损伤2022-11-30135050个入射个入射C C原子在表面的堆积情原子在表面的堆积情况与入射能量的依赖关系况与入射能量的依赖关系(a)1 eV(b)5 eV (c)10 eV(d)50 eV (e)100 eV(f)200 eV*red balls represent the C atoms *black ones represent W atoms.l表面损伤:损伤程度随着入损伤程度随着入射能量的升高而增加。其表射能量的升高而增加。其表征参数为对函数分布函数。征参数为对函数分布函数。l表面损伤的机理:动能释放动能释放与晶格碰撞。与晶格碰撞。ASIPP表面反射表面反射2022
15、-11-3014l10eV时,能量和粒子反时,能量和粒子反射在两种表面上均达到最大射在两种表面上均达到最大值。值。l大于大于10 eV,反射系数单反射系数单调递减,源于较高的能量,调递减,源于较高的能量,粒子更容易进入元胞内部。粒子更容易进入元胞内部。l低于低于10eV时,反射系数时,反射系数单调递增,源于较低的能量,单调递增,源于较低的能量,粒子更容易吸附到表面。粒子更容易吸附到表面。lsbe(surface binding energy):):8.68ev(W););esb(projectile binding energy):):4.00eV(C););1.00eV(H)(a)(a)粒子
16、和粒子和(b)(b)能量反射系数与入射能量依赖关系能量反射系数与入射能量依赖关系(Eckstein(Eckstein*用用TRIM.SP TRIM.SP 计算的结果也一并列出计算的结果也一并列出 )*W.Eckstein,Calculated Sputtering,Reflection and Range Values,IPP 9/132,2002 1101000.00.10.20.30.40.50.60.70.80.91.0 crystalline-W amorphous-W EcksteinEnergy Reflection CoefficinetIncident Energy(eV)(b
17、)1101000.00.10.20.30.40.50.60.70.80.91.0 Particle Reflection Coefficinet crystalline-W amorphous-W Eckstein(a)1100.30.40.50.60.70.80.91.0 MD EcksteinParticle Reflection Coefficient(a)1100.30.40.50.60.70.80.91.0 MD EcksteinEnergy Reflection CoefficientIncident Energy(eV)(b)C WH W ASIPP射程分布射程分布2022-11
18、-3015-2024681 01 21 4024681 01 21 41 e V Fraction of counts(%)D e p th ()(a)-2024681 01 21 4024681 01 21 4 Fraction of counts(%)D e p t h ()5 e V(b)-2024681 01 21 4024681 01 21 4 Fraction of counts(%)D e p th ()1 0 e V(c)-2024681 01 21 4024681 01 21 4 Fraction of counts(%)D e p th ()5 0 e V(d)-20246
19、81 01 21 4024681 01 21 4 Fraction of counts(%)D e p t h ()1 0 0 e V(e)11010005101520 Average depth()Incident Energy(eV)crystalline-W amorphous-W Eckstein l 低于10 eV,未被散射的原子吸附在表面的几个原子层不能穿透的体胞内部l 在10 eV,入射原子被散射回真空的几率最大l 高于10 eV,荷能粒子就有更大的几率进入的W的体内,同时平均射程分布以及射程摆动都会相应增加l在50eV附近,粒子在晶态内 平均射程超越了非晶表面,其原因在于沟道效
展开阅读全文