lammps教程:薄膜渗透模拟(3)–不同孔隙率对过滤效果的影响

lammps教程:薄膜渗透模拟(3)–不同孔隙率对过滤效果的影响本文是薄膜渗透过滤的最后一篇文章:不同孔隙率薄膜建模。孔隙或空位缺陷的建模原理比较简单:删除一定数量的原子就可以。lammps自带delete_atoms可以随机删除一定比例的原子,如果对孔隙或空位的形状、尺寸等有特殊需求,需要用编程的方法删除原子。delete_atomsporosity命令可随时产生设定比例的原子,如删除50%的原子:delete_atomsporositymembrane0.5482793membrane为原子组0.5为删除原子的比例482793为随机数种子

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE使用 1年只要46元 售后保障 童叟无欺

本文是薄膜渗透过滤的最后一篇文章:不同孔隙率薄膜建模。
孔隙或空位缺陷的建模原理比较简单:删除一定数量的原子就可以。
lammps自带delete_atoms可以随机删除一定比例的原子,如果对孔隙或空位的形状、尺寸等有特殊需求,需要用编程的方法删除原子。
delete_atoms porosity命令可随时产生设定比例的原子,如删除50%的原子:

delete_atoms porosity membrane 0.5 482793
  • membrane为原子组
    0.5为删除原子的比例 482793为随机数种子,取任一正整数即可。

初始模型和删除50%原子后的薄膜对比:
在这里插入图片描述
在这里插入图片描述

模拟结果为:
在这里插入图片描述

模拟in文件为:

#编写:Simon Gravelle
boundary	s p p
pair_style	lj/cut 2.5
neigh_modify 	every 1 delay 5 check yes
#读取模型文件
read_data	data.lammps
region		right block 0 INF INF INF INF INF
#原子分组
group		solid type 1 2 3
group		piston_right type 3
group		membrane type 2
group		piston_left type 1
group		fluid type 4 5
group		solvent type 4
group		solute type 5
#排除墙原子之间的作用力
neigh_modify 	exclude group solid solid
#按比例删除薄膜原子
region        membrane block -0.25 0.25 INF INF INF INF
delete_atoms     porosity membrane 0.5 482793
#温度初始化
velocity 	fluid create 1.0 4928459 mom yes rot yes dist gaussian
fix        mynve all nve
compute        temperature_fluid fluid temp
fix        mylgv fluid langevin 1.0 1.0 0.1 1530917 zero yes
fix_modify    mylgv temp temperature_fluid
thermo_modify    temp temperature_fluid
#固定薄膜
fix        mysf1 membrane setforce 0 0 0
fix        mysf2 piston_left setforce NULL 0 0
fix        mysf3 piston_right setforce NULL 0 0
#两侧施加压力
variable    F equal 0.01
fix        myaf1 piston_left aveforce ${ 
   F} NULL NULL
fix        myaf2 piston_right aveforce -${ 
   F} NULL NULL
#统计过滤后的溶质和溶剂原子数量
variable    solvent_right equal count(solvent,right)
variable    solute_right equal count(solute,right)
variable    position_piston_left equal xcm(piston_left,x)
variable    position_piston_right equal xcm(piston_right,x)
#输出统计记过
fix        myat1 all ave/time 10000 1 10000 v_solvent_right file solvent_right.dat
fix        myat2 all ave/time 10000 1 10000 v_solute_right file solute_right.dat
fix        myat3 all ave/time 10000 1 10000 v_position_piston_left file position_piston_left.dat
fix        myat4 all ave/time 10000 1 10000 v_position_piston_right file position_piston_right.dat
fix        myat5 all ave/time 10 1000 10000 f_mysf1[1] file force_membrane.dat
#输出轨迹文件
dump        mydmp all atom 10000 dump.xyz
thermo        10000
run        5000000

公众号:lammps加油站

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/197630.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)
blank

相关推荐

  • IntelliJ IDEA全局内容搜索和替换

    IntelliJ IDEA全局内容搜索和替换在做项目时,有时会在整个项目里或指定文件夹下进行全局搜索和替换,这是一个很方便功能。使用方法如下:一、全局搜索1、使用快捷键Ctrl+Shift+F打开搜索窗口,或者通过点击Edit–>Find–>Findinpath打开搜索窗口,如下图:2、搜索界面如下,主要分为上中下三部分,上部主要为搜索条件(要搜索的内容、范围、方式等),中间部分为包含搜索内容的文件列表,点击单个文件可以在下面部分

  • navicat连接mysql时报错1251怎么办

    navicat连接mysql时报错1251怎么办1、新安装的mysql8,使用激活成功教程版的navicat连接的时候一直报错,如图所示:2、网上查找原因发现是mysql8之前的版本中加密规则是mysql_native_password,而在mysql8之后,加密规则是caching_sha2_password。解决问题方法有两种,一种是升级navicat驱动;一种是把mysql用户登录密码加密规则还原成mysql_native_password。由于用的是激活成功教程版的navicat,所以只能用第二种方法解决了。3、首先win.

    2022年10月10日
  • linux修改密码长度限制_linux文件名长度限制修改

    linux修改密码长度限制_linux文件名长度限制修改修改密码长度:设置为不少于8位的。修改最短密码长度需要编辑login.defs文件(vi/etc/login.defs),把下面这行PASS_MIN_LEN5(默认的情况)改为PASS_MIN_LEN8(修改后的情况)login.defs文件是login程序的配置文件 转载于:https://blog.51ct…

  • fmp helper下载_sftp下载

    fmp helper下载_sftp下载1.下载安装https://teleport-pro.en.softonic.com/或者其他激活成功教程网站2.运行扒站打开软件后File->NewProjectWizard里面可以选择许多TP能做的事情,比如创建一个网站的副本,根据关键词查询网站等。下载网站选择第一个。输入想要下载的网址,和对应的下载深度。选择仅下载文本或是全部。其次如果该网站需要登陆的话需要输入登陆网…

  • 哈佛大学《幸福课》笔记

    哈佛大学《幸福课》笔记前述整个九月份以来一直在找工作,工作签好后有各种杂事需要做。当然,更主要的原因是自己懒,所以博客有一个多月没有更新了。今天将自己在大二暑假跟着视频学过的哈佛大学的《幸福课》整理成一篇文章,希望自己以后能够克服困难,多阅读,多学习,多写作,每天能够过的开心充实有意义。《幸福课》笔记(Tal的幸福秘籍)1、你问的问题经常会决定你所做的探索,问题会带来探索。2、信念创造现实,如何理解

  • web前端开发面试中常见的算法题(JS)

    web前端开发面试中常见的算法题(JS)前言最近在准备秋招,做过了大大小小的公司的面试题,发现除了基础知识外,算法还是挺重要的。特意整理了一些常见的算法题,添加了自己的理解并实现。除此之外,建议大家还可以刷刷《剑指offer》(但我还没刷完?,任重道远呐)。此外,左神在牛客网上也有算法课程,听了基础班的感觉还不错,起码让我这个算法小白也能快速地理解了很多问题,知识付费的时代,这个真的是良心课程了。就我个人而言的话,平时为了解决一…

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号