| 查看: 2456 | 回复: 0 | |||
[交流]
分享一个可以用于高精度模拟液态碳氢燃料相变的LAMMPS程序
|
|
液态碳氢燃料比如柴油和航空煤油中含有大量的直链烷烃(n-c10h22, n-c12h26等）。分子动力学（md）模拟可以帮助理解微纳米燃料液滴在航空推进系统中相变过程。联合原子势（united-atom potential，简称 ua 势）是一种md中常用的模拟碳氢燃料的势能模型，它将碳原子（c）和与其直接相连的氢原子（h）合并，看作是一个单一的、虚拟的“联合原子”作用中心，从而与全原子模型相比大大提升了技术速度。 目前常用的ua模型有：sks potential [1], nerd potential [2], 和trappe potential [3]。这三种势能模型的共同特点是能够比较准确的复现碳氢燃料的饱和液体密度，但在复现燃料的饱和蒸汽压和饱和蒸汽密度时通常存在30%到40%的误差。而且如果在md模拟中，截断半径（rcut）以外的lennard-jones（lj）势被简单忽略，饱和蒸汽压和蒸汽密度的md模拟结果会超出实验结果2到3倍。这样大的偏差将严重影响md模拟碳氢燃料相变过程的可信度。 为了解决这个问题，我们近期对ua模型中的参数做了优化调试，把原有的ua模型改进成一个improved united-atom (i-ua) model。在i-ua模型中，我们把原来的lj势改为截断平移lj势并调整了参数，让它更适合md模拟并且能精确复现碳氢燃料的热性质。我们用i-ua模型模拟了七种不同的直链烷烃(n-c4h10 ~ n-c22h46）的九种热物理性质。如下表所示，md模拟结果与实验数据的平均绝对偏差都在4.5%以内（粘度除外）： 热物理性质 平均绝对偏差 ———————————————— 饱和液体密度 1.6% 饱和蒸汽密度 2.3% 饱和蒸汽压 2.8% 临界温度 0.3% 临界密度 2.7% 临界压强 4.3% 表面张力 4.4% 蒸发潜热 2.1% 液体粘度 14.3% 我们在最近的文章中 [4]，提供了i-ua模型的具体参数以及怎样通过md模拟获得这些热物理性质。在这片文章的supplementary material里，我们提供了lammps程序。这个程序是用来算n-dodecane (n-c12h26)的饱和热物理性质的。为方便大家使用，我们对每一段程序都做了注解。大家可以很容易把它改成计算其他直链烷烃的热物理性质的lammps程序。如果有任何关于文章和lammps程序的问题，或者不能从jcp下载lammps程序，可以给我们发邮件zlch5@mst.edu。邮件题目请不要写中文，否则会被学校邮箱自动放入垃圾信箱。 顺便做一个招聘广告，我们组招一个全奖（一年~$30k+免学费）博士研究生做分子模拟微观热流。学校：密苏里科技大学（原密苏里罗拉大学），机械系全美排名60。要求托福成绩>80，gre成绩有就行。27年春季或以后入学都可以。感兴趣的直接发邮件到zlch5@mst.edu并附上简历。 [1] b. smit, s. karaborni, and j. i. siepmann, “computer simulations of vapor╟liquid phase equilibria of n-alkanes,” j. chem. phys. 102, 2126 (1995); “erratum: ‘computer simulations of vapor-liquid phase equilibria of n-alkanes,’” 109, 352 (1998). [2] s. k. nath, f. a. escobedo, and j. j. de pablo, “on the simulation of vapor╟liquid equilibria for alkanes,” j. chem. phys. 108, 9905 (1998). [3] m. g. martin and j. i. siepmann, “transferable potentials for phase equilibria. 1. united-atom description of n-alkanes,” j. phys. chem. b 102, 2569 (1998). [4] w. tausif, j. hartfield, a. george, and z. liang, an improved united-atom potential for molecular dynamics simulation of saturated properties of n-alkanes, j. chem. phys. 164, 044902 (2026) https://doi.org/10.1063/5.0306295 |
回复此楼» 猜你喜欢
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有11人回复
-
极限拉扯，如何抉择
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有10人回复
» 本主题相关商家推荐: (我也要在这里推广)
