高盐溶液蒸发过程中的分子动力学模拟与研究

随着我国工业的快速发展,如何减量化处理工业排放的高盐废水已经成为了亟待解决的环境问题.在处理高盐废水的各种工艺中,热蒸发技术因具备脱盐效果好,灵活性好等特点得到了广泛应用.本文从微观角度出发,采用分子动力学模拟方法研究了LiCl、KCl、CaCl₂三种溶液在400 K和500 K两种温度加热条件下的蒸发过程,分析了Li⁺、K⁺、Ca²⁺,Cl^-四种盐离子对蒸发速率、水分子取向、氢键、溶液结构等性质的影响.研究结果表明温度的提升会对蒸发速率产生极大影响,温度的提高不会改变配位水化层的位置但会明显减少离子的配位水分子数,有利于提高活跃水分子的占比和蒸发速率.     

高盐溶液蒸发过程中的分子动力学模拟与研究

随着我国工业的快速发展,如何减量化处理工业排放的高盐废水已经成为了亟待解决的环境问题.在处理高盐废水的各种工艺中,热蒸发技术因具备脱盐效果好,灵活性好等特点得到了广泛应用.本文从微观角度出发,采用分子动力学模拟方法研究了LiCl、KCl、CaCl_(2)三种溶液在400 K和500 K两种温度加热条件下的蒸发过程,分析了Li^(+)、K^(+)、Ca^(2+),Cl^(-)四种盐离子对蒸发速率、水分子取向、氢键、溶液结构等性质的影响.研究结果表明温度的提升会对蒸发速率产生极大影响,温度的提高不会改变配位水化层的位置但会明显减少离子的配位水分子数,有利于提高活跃水分子的占比和蒸发速率.     

纳米团簇熔化过程的分子动力学模拟

本文采用分子动力学结合嵌入原子多体势,模拟了不同半径的Ni纳米团簇的升温熔化过程,研究团簇尺寸对熔点和表面能的影响.模拟结果表明:团簇的熔点显著低于体材料的熔点.团簇熔化的过程首先是在团簇的表面出现预熔,然后向团簇内部扩展,直到整个团簇完全熔为液态.在模拟的纳米尺度范围内,团簇的熔点与团簇尺寸基本成线性关系.团簇的表面能随着团簇尺寸的增大而减小,而且表面能均高于体材料的表面能.     

纳米液滴撞击柱状固体表面动态行为的分子动力学模拟

液滴撞击固体表面是一种广泛存在于工农业生产中的现象.随着微纳技术的发展,纳米液滴撞击行为的定量描述有待完善.采用分子动力学模拟纳米水滴撞击柱状粗糙铜固体表面的动态行为.分别在液滴速度为2—15?/ps,五种方柱高度和六种固体表面特征能的情况下分析液滴的动态特征.结果表明,随着液滴初始速度V0的增加,其最终稳定状态先由C...     

多壁碳纳米管外壁高温蒸发的分子动力学模拟

采用经典分子动力学（MD）方法,使用EDIP（environment-dependent interatomic potential)势描述C纳米管内C原子之间相互作用,对多壁C纳米管由于Stone-Wales缺陷引起外层管高温剥落蒸发现象进行了计算模拟.研究结果表明,高温下多壁C纳米管外层管Stone-Wales缺陷处C原子剧烈振动导致C—C键断裂形成悬键,并逐渐向四周扩散导致外层管剥落蒸发.利用Lindemann指数作为判据,得出多壁C纳米管外层管出现剥落蒸发的温度为2290 K左右,与Huang Jianyu等实验中观测到多壁C纳米管外层管剥落蒸发现象产生的温度2000 ℃基本一致. 关键词： 多壁C纳米管 分子动力学 Stone-Wales缺陷 剥落蒸发     

多壁碳纳米管外壁高温蒸发的分子动力学模拟

采用经典分子动力学（MD）方法,使用EDIP（environment-dependent interatomic potential)势描述C纳米管内C原子之间相互作用,对多壁C纳米管由于Stone-Wales缺陷引起外层管高温剥落蒸发现象进行了计算模拟.研究结果表明,高温下多壁C纳米管外层管Stone-Wales缺陷处C原子剧烈振动导致C—C键断裂形成悬键,并逐渐向四周扩散导致外层管剥落蒸发.利用Lindemann指数作为判据,得出多壁C纳米管外层管出现剥落蒸发的温度为2290 K左右,与Huang Jianyu等实验中观测到多壁C纳米管外层管剥落蒸发现象产生的温度2000 ℃基本一致.     

丁二酸-水纳米气溶胶液滴表面张力的分子动力学研究

表面张力在纳米气溶胶颗粒的吸湿生长研究中具有重要意义,然而现有实验方法不能对其准确测量.本文基于分子动力学方法模拟了丁二酸气溶胶颗粒吸湿生长形成稳定液滴的动力学过程,在此基础上,建立模型计算了液滴的表面张力,进而探究了温度、粒径和丁二酸浓度对纳米液滴表面张力的影响机制.结果表明,随着温度从260 K升高到320 K,液滴内分子间作用力的减弱导致了液滴表面张力的减小,且表面张力的减小程度随丁二酸浓度的增大而增大,究其主要原因在于液滴中丁二酸分子的径向分布随温度和丁二酸浓度变化的差异;随着粒径的增大,液滴表面张力先增大后趋于定值,且粒径对表面张力的显著影响区间随着丁二酸浓度的增大而缩短;研究还发现,丁二酸分子的表面活性导致液滴表面张力随着丁二酸浓度的增大而减小,且减小趋势符合对数函数形式,尤其是在粒径小于6.12 nm时,同时,基于Szyszkowski公式对液滴的表面张力进行了拟合.本文研究成果能为气溶胶颗粒的吸湿生长和相关动力学过程预测理论及模型的改进提供参数依据.     

纳米Cu颗粒等温晶化过程的分子动力学模拟研究

采用分子动力学方法和镶嵌原子势,模拟了500个Cu原子(简称Cu₅₀₀) 组成的纳米颗粒的等温晶化过程.利用修正的均方位移、键对分析技术和内在结构(IS) 等方法对该过程中的结构和动力学行为进行分析研究.结果显示:与块体金属不同的是, Cu₅₀₀纳米颗粒在某一温度保温时,其晶化时间并不是一个定值, 而是存在一个统计分布,并且保温温度越低其晶化时间的分布范围越广, 最长晶化时间越长.在低温晶化时, Cu₅₀₀经历了一系列中间构型的转变才达到晶态, 表现出多步晶化的特征.文章作者研究了颗粒的初始构型对晶化进程的影响, 发现颗粒的初始结构特征和能量状态对其随后的晶化过程有着重要的影响, 同一温度下,颗粒初始构型的IS能量越低其晶化时间越长,这一点在低温时尤其明显.     

正癸烷纳米液滴润湿特性的分子动力学研究

流体液滴在固体表面的浸润性对其润滑性能至关重要.本文利用分子动力学方法研究了正癸烷纳米液滴在铜表面上的润湿特性.结果表明：在平坦光滑表面上,壁面的厚度和分子数目对润湿效果影响不大.随着壁面能量势阱参数εs 增大,接触角线性减小.随着温度升高,液滴的接触角减小.在沟槽粗糙表面上,随着粗糙度因子增大,对于疏液表面,接触角增大到一定值后基本保持不变,符合Cassie理论;中性和亲液表面接触角则会减小,为Wenzel润湿模式.当表面分数增大时,疏液与亲液表面接触角整体呈减小的趋势,对中性表面影响不大.当温度升高时,粗糙疏液表面接触角会增大,润湿效果更差,而粗糙中性和亲液表面液滴润湿性会更好.     

二甲基醚液滴超临界蒸发的数值模拟

建立了二甲基醚(DME)在N2介质中超临界蒸发的数学模型.采用状态方程法计算了DME-N2体系的高压气液相平衡,在相关文献的基础上得出一套流体热物性和输运参数的计算方法.完成了DME超临界蒸发过程的数值模拟.从数值模拟的结果讨论了超临界条件下环境压力、环境温度和液滴运动等因素对液滴升温、滴径变化以及生存时间的影响.     

分子动力学模拟壁面附近流体特性及降温过程

采用三相模型计算结果表明:近壁区内流体密度呈现出"层状化"的特点,切向应力张量呈"振荡"分布;固液界面张力随温度增加而降低;进一步增强壁面浸润性,固液界面张力变化不大;同温壁面对近壁区流体特性的影响较大,近壁区之外的流体特性与两相模型吻合较好.降低壁面温度的模拟表明,当温差减小时,传热量减小,界面接触热阻占主导地位.     

分子动力学模拟纳米晶体铜的结构与性能

通过分子动力学方法模拟纳米面心立方晶体铜的结构，并对模拟的结果进行了不同晶粒尺寸的纳米晶体的密度、能量分布以及弛豫前后的X射线衍射、径向分布函数等计算.结果表明，大体积分数的晶界和畸变的晶粒是纳米晶体有别于传统的粗晶粒晶体材料结构的重要方面，由此导致纳米晶体一系列不同性能. 关键词：     

扩散系数的分子动力学模拟

本文选用平衡分子动力学模拟方法,应用不同的势能模型对强极性分子水的扩散系数进行了模拟计算．结果表明,用分子动力学模拟方法计算水的扩散系数时,模拟结果对模型的选取极为敏感。在目前应用较为广泛的几种模型中以 SPCE模型较为适用,所得的扩散系数与实验值较为接近。     

水导热系数的分子动力学模拟

本文选用平衡分子动力学模拟方法,应用不同的势能模型对强极性分子水的导热系数进行了模拟计算.结果表明,用分子动力学模拟方法计算水的导热系数时,模拟结果对模型的选取极为敏感.在目前应用较为广泛的几种模型中以TIP4P模型较为适用,所得的导热系数与实验值较为接近.     

液态Au和Ag冷却过程的分子动力学模拟

采用反映原子多体相互作用的FS势模型,通过分子动力学方法研究了贵金属Au、Ag在温度2000~300 K之间的冷却过程,考察其结构变化特点,给出了冷却过程微观结构转变的重要信息.     

多形性相变的分子动力学模拟

 本文利用分子动力学方法研究了KCl晶体在ρ＝ρ₀时的温度相变。面心立方（fcc）和体心立方（bcc）两种结构的径向分布函数随温度的变化的情况说明，在高压下，发生着bcc结构相fcc结构的转变，bcc结构是不稳定的。为了选取合理的势参数，利用了分子动力学程序在T＝0时的性质，计算了NaCl和KCl晶体的零温状态方程，研究了它们在压力作用下发生的多形性相变。计算表明，NaCl和KCl晶体将分别于18.8 GPa和5.9 GPa发生从fcc到bcc的多形性相变，这些值相当接近实验结果。着说明本文选用的势参数势有一定精度的。     

液氩多体作用势研究及其Hugoniot曲线的分子动力学模拟

应用Hartree－Fock方法,计算了氩体系的多体相互作用.对Parson等提出的MSVⅢ势进行了 多体修正,并将修正后的MSVⅢ势用于液氩物态方程的分子动力学模拟.计算表明,修正后 的MSVⅢ势能较好地反映液氩在所研究压力范围的情况,计算的Hugoniot曲线与实验结果在 冲击压力高达40?GPa密度区的符合程度得到改善.     

分子动力学模拟铜薄膜的热导率

采用分子动力学(MD)方法模拟铜薄膜的热导率,给出了厚度在100~400nm、温度在100~600K范围内铜薄膜热导率对尺寸及温度的依赖关系.     

多核环境下的分子动力学模拟

本文在多核环境下,使用OpenMP实现了经典分子动力学模拟程序的并行;同时对分子动力学模拟进行了两项主要的优化:分子排序及运用SIMD指令运算.在4核下获得了4.13倍的计算性能提升,将经典分子动力学模拟的模拟规模提高至4000分子×10~7模拟总步数.     

纤维素热解过程的分子动力学模拟

本文利用 PM3 方法对生物质主要组分纤维素进行了结构优化并得到了一系列的结构参数.对不同力场下,聚合度为 9 的纤维素单链热分解进行了分子动力学方法模拟研究,得到不同力场下的模拟过程中参数.通过与相关生物质热解实验结果对比,利用 Amber 力场模拟得到的结果与实验值吻合较好.基于 Amber 力场并结合量子力学对纤维素单元热解过程进行了研究,模拟得到纤维素单元分子链在加热过程中的主要分解温度范围、断键顺序以及一次热解的基团,并对一次产物进行了分析.