分子动力学模拟超临界水微观结构及自扩散系数 microstructure and self-diffusion coefficient of supercritical water by molecular dynamics simulation.pdfVIP

{-。‘-v砷P 、i4毒。 分子动力学模拟超临界水微颜结构及自扩散系数 张乃强，徐鸿，白杨 (华北电力大学电站设备状态监测与控制教育部重点实验室，北京 102206) “瀛 摘要：采用分子动力学模拟方法．选用COMPASS力场，对超临界水的微观结构捐扩散系数进行了研究．： 结果表明，在超临界区水分子的氢键作用明屁减弱，但是仍有氢键存在。水分子的键}乇和键角存温度为 648～748 K时变化迅速，在温度为748～973 K时随 K时变化很小。径向分布函数特征峰值在温度为648～748 温度的升高而增大，在温度为748—973K时随温度升高而减小，表明水分子近程结构有序性随温度变化。 超临界水自扩散系数是常温常压水的2个数量级，随温度升高而增大．随压力的升高而减小。 关键词：分子动力学；超临界水；微观结构：自扩散系数 中图分类号：TK211 文献标志码：A 文章编号：1004．9649(2011)12-0047．04 组合而来的，未做严格考核，而COMPASS模型中非 0引言 键相互作用势雨数是从分子液体和品体中不断优化 出来的，更全面地考虑了分子内和分子间的相互作 超临界水作为清洁高效的反应介质得到了深入 用。力场中总能量的计算表达式为： 研究和广泛应用…．如超临界水氧化．超临界水气化 制氢，超临界机组发电等。超临界水微观结构和氢键 的变化引起很多宏观物性的奇异变化，2]．采用试验 该能量表达式由成键作用项和非键作用项2部 方法研究超临界水的微观结构和传递性质通常需要 分组成。成键项包括键、角、二面角和平面扭曲作用 很苛刻的试验条件和高昂的费用。随着计算机技术 的位能既，E。，瓦，E，分子内相互作用项有‰，，‰， 和理论化学的飞速发展．分子动力学模拟已经成为 ‰，风，，E眦，可以描述振动频率和结构变化。非键作 与理论和试验研究并行的一种研究方法【一J． 9—6位能 用项包括范德华相互作用力的“柔性”L—J 近年来。对水的分子动力学模拟进行了大龟的 项和计算静电作用力的Coulombic项． 研究17…。大多数研究都局限于648～748K的超临 氢键的计算实际上是通过综合考虑静电和范德 界水。 华力作用来实现的．文献『16]采用不同力场进行计 随着超超临界机组的投入运行和超高参数超临 算对比，得到COMPASS力场准确性最好，pcff力场 界机组的研制，超临界水的微观结构和扩散性质显 较合理，而cvff力场不能准确地反映氢键。因此． 得更加重要。因此，本文采用分子动力学的方法研究 本文采用CoMPASS力场模拟超临界水． 温度为648～973K。压力分别为25MPa和30MPa 的超I临界水微观结构和自扩散系数． 2模拟细节 1位能模型及力场