碳纳米管对氢气的吸收模拟.ppt

CCMS课题设计报告单壁碳纳米管吸附氢气的计算模拟 目录 小组成员 课题介绍 模拟设计与结果分析 课程设计小结 程序演示 V_r拟合 化学位与压力的关系 一维模型 小组成员 （E3-EST005） 组 长： 龚红波 刘文庆 其他成员： 杜宝贵 杨 梅 李 茜 雷光华 林 杰 黄 盛 肖 敏 朱红伟 （1）参考图 * back 课题介绍 1.“望氢兴叹” 为了能达到和现有燃油汽车相当的存储量,氢气存 储的重量密度需达到6.2%(氢气重量/系统重量),体积 密度需达到65%(氢气质量/系统体积). 氢气主要通过压缩、液化、转变为金属氢化物和物理吸附等四种方式进行存储.但目前尚无一种方式能达到上述存储目标. 课题介绍 图1 锯齿管结构 2.碳纳米管： 纳米级的同轴管状炭纤维； 1991年5月,日本ＮＥＣ公司 的饭岛澄男(Ｓ.lijima)博士首次 发现； 具有很大的表面积,长径比可 达到100～1000； 在氢气存储方面很大的潜力。 课题介绍 3.计算机分子模拟（优点） 大大降低实验成本； 有可重复性和能模拟真实实验很难达到的超高压和超低温的情况； 旨在通过计算机模拟氢气在碳纳米管中的吸附, 揭示影响吸附的内部机理,如单壁碳纳米管半径、管间距等，进行理论指导。 课题介绍 三个基本步骤： （1）粒子坐标的变化,即move,在该过程中接受概率ζmov为min(1,exp(-ΔＵ/ＫBＴ))。从所示GCMC方法的角度来看,首先假设一个具体粒子的位置发生变化,相应整个系统的能量改变,满足给定条件则移动成功,做大量的模拟后得到粒子在管中的运动轨迹。 4.GCMC方法 巨正则系综的蒙特卡罗方法 （2）粒子的产生 被称为create,接受概率ζcre为 min(1,[aV /(Ｎ+1)]exp(-ΔＵ/ＫBＴ)) （3）粒子的消灭 被称为destory,接受概率ζdes为 min(1,(aV/Ｎ)exp(-ΔＵ/ＫBＴ)) 课题介绍 氢气分子与碳原子之间的相互作用采用L-J模型.其 中,将氢气对整个碳纳米管中所有碳原子的作用进行加 和,以获得氢气分子与单根给定半径的碳纳米管的相互 作用模型. 课题介绍 back 系统的化学位：模拟中，设为所有氢分子的势能和。 L-J势能模型及势能曲线的拟合 氢气分子分布的模拟 模拟设计与结果分析 1.L-J势能模型及势能曲线的拟合 模拟设计与结果分析 拟合 (2)参考公式:为了方便编程，将参考图的结果拟合成 多项式形式: (3)公式中的疑问： 变量在哪？如何体现其“势能随距离变化而变化”？ 所谓参数A、B、C、D、E 是多少？ 模拟设计与结果分析 (4)由相关资料分析，猜测：上式中的R即为参考图中 的 r/R ； (5)以变量 r/R 代替拟合公式中的 R，采用LabVIEW 编程，适当调整A、B、C、D 四个参数直到模拟的图 形与参考图形相吻合为止。 模拟设计与结果分析 (6)图形曲线 其中红线对应b,表示管外氢分子的势能随r/R变化情 况，白线对应a，表示管内氢分子的势能的变化情况。 模拟设计与结果分析 (6)模拟公式及参数 a) 管内情况 A= 0.000285 B=-1.1650 C= 9.5900 D=-24.8000 Z= 0.3179 模拟设计与结果分析 *为方便编程，用Z表示σCH（0.3179 nm） b) 管外情况 A= 0.000264 B=-1.0980 C= 8.9800 D= -18.8000 Z= 0.3179 模拟设计与结果分析 模拟设计与结果分析 (7)图形分析 通过模拟与图形对比知： 所作猜测基本正确，即 R 应为 r/R ; 所尝试参数比较合理，所的结果可作为势能的产生公式； 管内、外氢分子所受纳米管作用随距离的变化有相似 的变化趋势，碳氢分子相距很近时，斥力其主要作用， 随距离增大而减小；距离较大时引力其主要作用，一 定范围内随距离增大而增大。引力与斥力平衡处作用 势应最小，平衡点约为0.058R. 2.氢气分子分布的模拟 模拟设计与结果分析 （1）氢在纳米管中的轴向分布 轴向分布 径向分布 由图可见，氢在管中的分布与 管子的直