国人团队发文 | 这个体系也能跑MD模拟！

文献解读

搞学术不秃头 发文章更轻松

派言派语

液滴颗粒碰撞是许多技术应用中普遍存在的现象。在流化催化裂化(FCC)工业中，热的固体催化颗粒在FCC反应器中与雾化重质燃料液滴发生碰撞。此外，颗粒与液滴的碰撞在环境保护领域也备受关注。

迄今为止，人们对颗粒与液滴碰撞进行了广泛研究，主要集中在碰撞结果、碰撞动力学及其影响因素等方面。

实验研究通常采用液滴撞击静止颗粒的方式，通过高速摄像机捕捉液滴变形。由于长度和时间尺度较小，大多数实验采用大液滴-颗粒体系，影响因素有限。

数值模拟方法可弥补现有方法在捕捉碰撞细节方面的不足。但这些方法无法在分子尺度上准确计算气液界面与薄层。

近期，利用MD模拟对液滴行为进行了大量研究。这些研究具有无需事先假设、影响参数变化范围大等优点。

2023年3月16日，东南大学能源与环境学院杨林军教授课题组在《Powder Technology》上发表研究成果。

他们对液滴颗粒体系的正面碰撞进行了分子动力学模拟，使用的模拟软件是LAMMPS。通过改变广泛的影响参数，提出了不同的碰撞情况，并将其分为7组。

对薄膜厚度、涂覆面积和蒸发量等重要动力学参数的演变规律也进行了深入研究。MD模拟的结果也与之前宏观实验的结果进行了比较。

PDF原文后台回复“0329”获取。

关于LAMMPS软件，后台回复“lammps”，免费获取超详细安装教程。

分子动力学

分享关于Molecular Dynamics的一切。

模拟方法

1、原子间势

模拟中使用的液滴主要由纯水组成，因此其结果可以与使用纯水液滴的现有实验进行比较。使用LAMMPS对液滴颗粒碰撞进行了模拟。原子间的非键合相互作用遵循经典的Lennard-Jones 12-6势和库仑势。

2、初始构型

液滴颗粒碰撞的初始构型如图1所示。整个系统由液滴、颗粒和载气三部分组成。所有分子均置于模拟箱（300Å×300Å×600Å）中。

图1：初步配置液滴颗粒碰撞体系

在三个方向上都施加周期性边界条件。考虑到计算成本和代表性，液滴初始直径固定为8nm。

使用PACKMOL软件和LAMMPS中内置的建模命令创建所需的初始构型。最后，利用OVITO对初始构型与粒子轨迹进行了可视化。

3、模拟细节及后处理

液滴颗粒碰撞是一个复杂的物理过程，涉及惯性力、粘性力等之间的竞争。具体的模拟细节参数如表2所示。韦伯数的范围是精心选择的，因此在不同初始韦伯数下，液滴颗粒碰撞结果会有很大差异。

表2：液滴-颗粒碰撞的模拟细节

系统中原子的初始速度从高斯分布中随机分配，参考温度为293K。所制备的模拟箱首先在293K的NVT系综中平衡1ns以达到理想状态。

系统平衡后，弹簧力被移除，液滴中的原子增加了z方向的新速度。同时，环境氮气和颗粒温度升高至设定温度值。

环境温度采用Nosé-Hoover恒温器控制，颗粒温度采用郎之万方法控制。

图2：液膜厚度hf和等效扩散直径Ds示意图

研究总结

该研究采用分子动力学模拟方法对亚微米尺度下的液滴颗粒碰撞过程进行了研究。碰撞结果及其影响因素被赋予特别的重要性。

综合考虑到大范围内三个碰撞参数（0.5≤Ω≤1.5，45°≤θ≤145°，30≤We≤250）的影响，构建了一个描述碰撞结果的综合体系图。

主要结论如下：

1) 共发现7种碰撞结果，即沉积、带孔沉积、全涂层、韧带结构解体、锥形结构解体、反弹和飞溅。

2) 对液相平均速度vz*进行采样，分析整个碰撞过程中动能的耗散情况。对于反弹情况，液滴碰撞后会恢复约14%的初速度，大部分动能因粘性而耗散，而vz*在锥形结构解体中受影响较小。

3) 液滴颗粒碰撞过程中液滴蒸发速率与环境温度、初始韦伯数和表面润湿性有关。在这三个因素中，发现初始韦伯数是最重要的一个。

4) 在MDS中，沉积-分离制度的边界发生在比宏观研究中高得多的韦伯数上，临界韦伯数从75左右增加到200。重力的减弱以及粘性耗散和表面张力的增强导致了临界韦伯数的不同。

总而言之，本研究者对运动液滴与固定颗粒的正面碰撞进行了初步研究。今后将在斜碰撞的微观检验方面投入更多的精力。此外，还将考虑粒子运动的影响，包括平移与旋转。

号外号外：广告实在是无孔不入+防不胜防啊！所以不直接放群聊二维码了，想进分子动力学交流群的，请添加小编微信并备注昵称+专业，手动拉你们入群！

征稿启事

“分子动力学”自创号以来得到了广大科研人士的关注和支持。为更好地服务分子动力学研究和应用，本公众号现长期招聘供稿作者。

真诚欢迎MD相关科研爱好者加入“分子动力学”团队！可发送简历至likeapoem@163.com，邮件主题请注明“姓名+供稿作者”。

来源：分子动力学