计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: 牧生
打印 Print 上一主题 Last thread 下一主题 Next thread

[GROMACS] 模拟水滴砸在金板上,低速下水滴会散开,且砸在金板上以后,水分子飞溅的太远

[复制链接 Copy URL]

395

帖子

8

威望

3908

eV
积分
4463

Level 6 (一方通行)

石墨

16#
发表于 Post on 2024-10-4 22:47:12 | 只看该作者 Only view this author
本帖最后由 Graphite 于 2024-10-4 22:49 编辑

用Lammps,fix setforce+fix gravity...原生就支持这种模拟,金板可以固定、可以锁定金板质心Z轴位置(但允许各个原子的振动),也可以像跳水板一样,一边固定,另外一边受冲击,等等。可以同时金板用多体势,可调硬度,水用MDPD,能溅射成若干个小水滴,等等。
自在飞花轻似梦,无边丝雨细如愁。

全自动反应动力学(ReaxFF、AIMD、NEP等)后处理工具网页版:http://cc-portal.xyz/reax_tools

21

帖子

0

威望

169

eV
积分
190

Level 3 能力者

17#
发表于 Post on 2024-10-8 18:39:31 | 只看该作者 Only view this author
牧生 发表于 2024-10-3 19:42
http://sobereva.com/153     根据这个博文的操作,建立一个金板,保存为pdb,如果金的参数抄这个博文,那 ...

感谢大佬,我在模拟中(100m/s,gmx2020),出现液滴保持形态下落,撞击过程也正常,但是液滴在撞击之后铺展碎裂,sob老师给我的意见是降低速度,我在阅读文献时看到其他人(lammps mw模型)在更大速度下也可以保持液滴完整性,以及液滴撞击回缩过程,这是什么原因呢?

21

帖子

0

威望

169

eV
积分
190

Level 3 能力者

18#
发表于 Post on 2024-10-8 18:48:31 | 只看该作者 Only view this author
Graphite 发表于 2024-10-4 22:47
用Lammps,fix setforce+fix gravity...原生就支持这种模拟,金板可以固定、可以锁定金板质心Z轴位置(但允 ...

好的,谢谢大佬

1478

帖子

0

威望

4537

eV
积分
6015

Level 6 (一方通行)

19#
 楼主 Author| 发表于 Post on 2024-10-8 20:14:41 | 只看该作者 Only view this author
本帖最后由 牧生 于 2024-10-8 21:07 编辑
Layman 发表于 2024-10-8 18:39
感谢大佬,我在模拟中(100m/s,gmx2020),出现液滴保持形态下落,撞击过程也正常,但是液滴在撞击之后 ...

液滴的速度过大的话,你用energy看一下液体的温度,这个温度是很高的。如果温度过高,在金属板上全部溅开,是正常的。

如果速度比较小,直接就滴落在金属板上,不溅开,温度也不高,符合直觉。

液滴在金属板上弹跳,我目前也没做出来,自我分析了一下,大概是因为水滴太小的缘故。
又菜又爱玩

21

帖子

0

威望

169

eV
积分
190

Level 3 能力者

20#
发表于 Post on 2024-10-8 21:13:21 | 只看该作者 Only view this author
牧生 发表于 2024-10-8 20:14
液滴的速度过大的话,你用energy看一下液体的温度,这个温度是很高的。如果温度过高,在金属板上全部溅开 ...

感谢大佬,液体温度的判断有没有一个区间范围呢?我看了下自己的模拟,感觉应该是液滴尺寸过小,速度过大,导致液体温度过高(600k左右),产生了飞溅现象

395

帖子

8

威望

3908

eV
积分
4463

Level 6 (一方通行)

石墨

21#
发表于 Post on 2024-10-9 15:12:23 | 只看该作者 Only view this author
本帖最后由 Graphite 于 2024-10-9 15:17 编辑
Layman 发表于 2024-10-8 18:39
感谢大佬,我在模拟中(100m/s,gmx2020),出现液滴保持形态下落,撞击过程也正常,但是液滴在撞击之后 ...

#以下为暴论
LAMMPS(或者自己开发任何基于动力学积分的软件),搞整活模拟就是在搞一种数值游戏,更粗暴地说,台球和弹簧的游戏,只要你会调、并且能合理化。
无非就是力(固定关系的键长键角二面角,不固定关系的相互作用势——LJ/coul/morse/eam/buckingham/soft/dpd...)、运动方程、外部条件。

如果能自洽、形成大的知识体系,解决不少问题,我们称之为一种范式,比如蛋白质的力场+1-2 fs的步长+平衡态NPT/NVT模拟——常见的蛋白质模拟套路。各类合金力场+切削/拉伸/扭曲/压痕/退火等一系列丝滑小连招——常见的材料的套路。
如果还不能自洽,那就得想办法从一定依据、佐证、实验现象出发,去调整和合理化这套东西。然后创造出一堆名词,XXMD之类的东西,有些经典永流传,有些纯属自嗨。

速度温度只是一种参数,可以理解为变动率。

自在飞花轻似梦,无边丝雨细如愁。

全自动反应动力学(ReaxFF、AIMD、NEP等)后处理工具网页版:http://cc-portal.xyz/reax_tools

21

帖子

0

威望

169

eV
积分
190

Level 3 能力者

22#
发表于 Post on 2024-10-9 16:55:31 | 只看该作者 Only view this author
Graphite 发表于 2024-10-9 15:12
#以下为暴论
LAMMPS(或者自己开发任何基于动力学积分的软件),搞整活模拟就是在搞一种数值游戏,更粗 ...

多谢大佬,我目前也是刚接触md,正处于摸索阶段,谢谢大佬的解答。

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-13 01:43 , Processed in 0.673275 second(s), 20 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list