应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 关于lammps 计算温度时去除宏观速度的疑问
51/1返回列表
查看: 3923  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

7long

银虫 (小有名气)

[求助] 关于lammps 计算温度时去除宏观速度的疑问 已有3人参与

关于lammps 计算温度有时去除宏观速度，即subtract a bias?这在shear的例子里有体现，但是某一方向既有宏观速度又有热运动速度啊，
计算温度时简单去掉这个方向的速度分量会不会太粗暴了？况且热运动的速度量级似乎应在~100m/s，远高于宏观速度，这样宏观速度的
影响应该很小？
    求指教，谢谢！
回复此楼

» 收录本帖的淘帖专辑推荐

分子动力学模拟

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2016-08-23 09:03:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

7long

银虫 (小有名气)
引用回帖:
8楼: Originally posted by 一次性保温杯 at 2016-08-23 14:45:57
还有我们所研究的shear并不是真正的shear。也不会产生线性变化分布的速度，只是模拟剪切变化过程中的每个构象的动力学，通过构象-能量关系，计算出来的。
...

shear的例子in文件如下

# 3d metal shear simulation

units                metal
boundary        s s p

atom_style        atomic
lattice                fcc 3.52
region                box block 0 16.0 0 10.0 0 2.828427
create_box        3 box

lattice                fcc 3.52 orient        x 1 0 0 orient y 0 1 1 orient z 0 -1 1 &
                origin 0.5 0 0
create_atoms        1 box

pair_style        eam
pair_coeff        * * Ni_u3.eam

neighbor        0.3 bin
neigh_modify        delay 5

region                lower block INF INF INF 0.9 INF INF
region                upper block INF INF 6.1 INF INF INF
group                lower region lower
group                upper region upper
group                boundary union lower upper
group                mobile subtract all boundary

set                group lower type 2
set                group upper type 3

# void

#region                void cylinder z 8 5 2.5 INF INF
#delete_atoms        region void

# temp controllers

compute                new3d mobile temp
compute                new2d mobile temp/partial 0 1 1

# equilibrate

velocity        mobile create 300.0 5812775 temp new3d
fix                1 all nve
fix                2 boundary setforce 0.0 0.0 0.0

fix                3 mobile temp/rescale 10 300.0 300.0 10.0 1.0
fix_modify        3 temp new3d

thermo                25
thermo_modify        temp new3d

timestep        0.001
run                100

# shear

velocity        upper set 1.0 0 0
velocity        mobile ramp vx 0.0 1.0 y 1.4 8.6 sum yes

unfix                3
fix                3 mobile temp/rescale 10 300.0 300.0 10.0 1.0
fix_modify        3 temp new2d

dump                1 all atom 100 dump.shear

#dump                2 all image 100 image.*.jpg type type &
#                axes yes 0.8 0.02 view 0 0 zoom 1.5 up 0 1 0 adiam 2.0
#dump_modify        2 pad 4

#dump                3 all movie 100 movie.mpg type type &
#                axes yes 0.8 0.02 view 0 0 zoom 1.5 up 0 1 0 adiam 2.0
#dump_modify        3 pad 4

thermo                100
thermo_modify        temp new2d

reset_timestep        0
run                3000
赞一下 回复此楼
9楼2016-08-23 15:10:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

一次性保温杯

金虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
你所说的宏观速度是体系的质心移动速度吗，这个速度本来就是要消去的，要不然模拟的盒子会一直移动啊，就计算错误了
赞一下(1人) 回复此楼
2楼2016-08-23 09:24:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

7long

银虫 (小有名气)
引用回帖:
2楼: Originally posted by 一次性保温杯 at 2016-08-23 09:24:34
你所说的宏观速度是体系的质心移动速度吗，这个速度本来就是要消去的，要不然模拟的盒子会一直移动啊，就计算错误了

也不完全是质心的速度啊，感觉是外界加载的一个速度，比如说剪切，或者纳米压痕里的向下压，但加载方向也存在热速度啊，如何将这两种速度分离开呢？简单都去掉控温感觉不合理啊。
你看过shear的那个例子吗？话说怎么去掉质心的速度呢？
赞一下 回复此楼
3楼2016-08-23 10:53:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

一次性保温杯

金虫 (小有名气)

【答案】应助回帖

引用回帖:
3楼: Originally posted by 7long at 2016-08-22 18:53:23
也不完全是质心的速度啊，感觉是外界加载的一个速度，比如说剪切，或者纳米压痕里的向下压，但加载方向也存在热速度啊，如何将这两种速度分离开呢？简单都去掉控温感觉不合理啊。
你看过shear的那个例子吗？话说怎 ...

去掉质心速度是初始化速度时候做的，因为随机生成的速度不能保证质心平衡，所以会去掉质心速度，随后质心自己就不会在移动的（因为外界没有施加力）
赞一下 回复此楼
4楼2016-08-23 11:11:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭