首页/文章/ 详情

基于LAMMPS实现页岩纳米孔多组分气体和页岩油的竞争吸附

11月前浏览982

关键词:分子动力学,吸附,竞争吸附,CH4/CO2/H2,页岩油,LAMMPS,GROMACS

摘要:MS软件可视化程度高,操作简便,但是速度慢,效率低。为此,我们使用LAMMPS/GROMACS软件进行气体和页岩油吸附,竞争吸附的研究。我们基于LAMMPS/GROMACS研究不同页岩孔隙(石英,蒙脱石,伊利石,高岭石,方解石,石墨烯,真实干酪根)中气体(CO2,CH4,H2,N2等)和 页岩油(NC6,NC8,NC12等)的吸附和竞争吸附行为,涉及到不同壁面类型,不同宽度,不同油气类型,不同摩尔比例,不同温度压力的设置,可以分析模拟快照,密度分布,扩散系数,相互作用能等。

整个分子动力学过程通过Lammps实现,其中部分壁面需要在MS中进行调整。全部流程如下:

1,建立壁面模型(干酪根,石墨烯,二氧化硅,蒙脱石,高岭石,伊利石,方解石等);建立原油组分分子结构;建立注入气(CO2, CH4, N2)的分子结构;

2,赋予干酪根CVFF力场,粘土矿物ClayFF力场,CO2, N2分别用fix rigid设为刚体,CH4用联合原子/OPLS力场,原油组分用OPLS-AA力场。

3,利用in文件将壁面,油,CO2气体等组合成data文件。

4,进行能量最小化,设置平板移动速度进行分子动力学模拟。

可以提供代码和一对一教学,可以进行答疑,可以针对特定需求来建模和代码。

可以在论文中分析的图和数据有:

1 初始模型


2 平衡模型(不同摩尔比)


3 密度分布


4 气体吸附量


5 相互作用力(静电力与范德华力)


6 不同含水量的吸附图像


7 不同含水的密度分布


8 不同含水下的气体吸附量


9 气体损失率


10 扩散能力


11 封存效率


12 力场验证


13 吸附气/自由气判定


14 模型随机性


来源:320科技工作室
电力油气GID分子动力学LAMMPSGROMACS
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-10-19
最近编辑:11月前
320科技工作室
硕士 | 结构工程师 微❤️ CAE320
获赞 243粉丝 447文章 503课程 0
点赞
收藏
作者推荐

基于COMSOL计算扭曲光子晶体中偏振可调的BIC

连续体中的束缚态(BIC)由于其非辐射的特点已被证明是相关的在动量空间中具有拓扑电荷和涡旋极化奇点。对于常规对称的光子晶体板中,BIC被线极化远场所包围,不利于高容量和多功能集成光学应用。动量空间中如何调控其周围极化偏振是一个有趣的问题。利用COMSOL来复现一篇国产小子刊,题为“Arbitrarily polarized bound states in the continuum with twisted photonic crystal slabs”。本文采用扭曲光子晶体结构实现了BIC附近线偏振往圆偏振的转化,BIC附近的偏振态具有相同的椭偏率。椭偏率由扭转角度直接决定。如下图所示:图1:扭转光子晶体采用本征求解器,首先采用挖孔结构完成建模,材料采用介质硅,设置一定的空气层高度并上下添加完美匹配层。前后左右的边界条件采用周期性边界条件,kx和ky代表x和y方向的波矢,如下图所示。图2:建模图3:周期性边界条件设置扫描波矢kx,将ky设为0,我们计算得到方向的能带,如下图所示。此处纵坐标为波长,横坐标为波矢kx。插图为原文三维能带,结果一致。品质因子在点处趋于无穷大,证明了BIC的存在。图4:能带计算图5:品质因子计算我们给予一定的扭转角度,设为45°,同时扫描波矢kx和ky。在上表面添加积分算子来计算斯托克斯参数,具体计算公式参考[ Physical Review Letters, 124(15), 153904.]的补充材料。最后绘制动量空间中的极化分布,此时BIC全部为圆偏振。图6:椭偏率计算图7:圆偏振分布来源:320科技工作室

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈