24小时热门版块排行榜    

查看: 1802  |  回复: 7

追梦2016

铜虫 (小有名气)

[求助] 加热区域问题 已有1人参与

我想计算加热区域对距离热源不同位置的原子势能的影响程度,求相关lammps范例,谢谢。
回复此楼

» 猜你喜欢

只要肯学,没有学不会的
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alundilong

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
追梦2016: 金币+10, ★★★★★最佳答案 2016-07-28 09:20:27
Thermostatting Regions in LAMMPS
四行命令实现:
fix hot_rescale all temp/rescale 10 3.0 3.0 0.0001 1.0
region hot block 0 10 0 20 0 20 units box
compute hot_temp all temp/region hot
fix_modify hot_rescale temp hot_temp

详见:
https://niallpjackson.wordpress. ... -regions-in-lammps/
2楼2016-07-27 05:10:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

追梦2016

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by alundilong at 2016-07-27 05:10:31
Thermostatting Regions in LAMMPS
四行命令实现:
fix hot_rescale all temp/rescale 10 3.0 3.0 0.0001 1.0
region hot block 0 10 0 20 0 20 units box
compute hot_temp all temp/region hot
fix_modify h ...

您好,首先感谢您提供的这个例子,关于这个例子有个问题,就是下图中为什么用all呢,因为我想计算的是某一组原子中,加热某一区域原子,观察这个温度对周围原子的影响,如果用all的话,是不是对整个组进行加热了?
加热区域问题
QQ截图20160727093847.png

只要肯学,没有学不会的
3楼2016-07-27 09:47:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

追梦2016

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by alundilong at 2016-07-27 05:10:31
Thermostatting Regions in LAMMPS
四行命令实现:
fix hot_rescale all temp/rescale 10 3.0 3.0 0.0001 1.0
region hot block 0 10 0 20 0 20 units box
compute hot_temp all temp/region hot
fix_modify h ...

还有个问题 ,就是我用的下面的命令:
region           jiare block 0 30 0 5 0 30
......
fix              1 all nve
fix              hot_rescale all temp/rescale 10 300.0 300.0 0.01 300.0
compute          hot_temp all temp/region jiare
fix_modify       hot_rescale temp hot_temp
但运行的结果显示了:computed temperature for fix temp/rescale cannot be 0.0的错误,这个是什么原因造成的呢?谢谢
只要肯学,没有学不会的
4楼2016-07-27 11:02:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alundilong

木虫 (正式写手)

【答案】应助回帖

引用回帖:
3楼: Originally posted by 追梦2016 at 2016-07-26 19:47:30
您好,首先感谢您提供的这个例子,关于这个例子有个问题,就是下图中为什么用all呢,因为我想计算的是某一组原子中,加热某一区域原子,观察这个温度对周围原子的影响,如果用all的话,是不是对整个组进行加热了? ...

至于为什么用all:
楼主不必担心,fix_modify 会把这个fix中计算temp的方式进行修正的。欲知详情,可以看看code
5楼2016-07-28 01:56:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alundilong

木虫 (正式写手)

【答案】应助回帖

引用回帖:
4楼: Originally posted by 追梦2016 at 2016-07-26 21:02:27
还有个问题 ,就是我用的下面的命令:
region           jiare block 0 30 0 5 0 30
......
fix              1 all nve
fix              hot_rescale all temp/rescale 10 300.0 300.0 0.01 300.0
compute  ...

看看你的初始速度是否为0?又是值region内部原子的。
6楼2016-07-28 01:57:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

追梦2016

铜虫 (小有名气)

引用回帖:
6楼: Originally posted by alundilong at 2016-07-28 01:57:21
看看你的初始速度是否为0?又是值region内部原子的。...

嗯嗯。对了,还想请教您个问题,就是我能否一下子将温度加到某个值呢?我想模拟恒温下原子的变化情况,但每次都是温度从140多度慢慢加热到一定温度的。
只要肯学,没有学不会的
7楼2016-07-28 09:20:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alundilong

木虫 (正式写手)

引用回帖:
7楼: Originally posted by 追梦2016 at 2016-07-27 19:20:01
嗯嗯。对了,还想请教您个问题,就是我能否一下子将温度加到某个值呢?我想模拟恒温下原子的变化情况,但每次都是温度从140多度慢慢加热到一定温度的。...

问题1:可以
问题2:不知道你在说啥。。。

发自小木虫IOS客户端
8楼2016-07-28 11:35:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 追梦2016 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 280分求调剂 一志愿085802 +4 PUMPT 2026-03-22 7/350 2026-03-22 22:13 by 星空星月
[考研] 289材料与化工(085600)B区求调剂 +3 这么名字咋样 2026-03-22 4/200 2026-03-22 17:56 by 云民大李老师
[考研] 287求调剂 +8 晨昏线与星海 2026-03-19 9/450 2026-03-22 17:01 by i_cooler
[考研] 305分求调剂(食品工程) +4 Sxy112 2026-03-21 6/300 2026-03-22 15:26 by 无懈可击111
[考研] 260求调剂 +3 朱芷琳 2026-03-20 4/200 2026-03-22 15:12 by 朱芷琳
[考研] 085600材料与化工306 +4 z1z2z3879 2026-03-21 4/200 2026-03-21 23:44 by ms629
[考研] 280求调剂 +11 咕噜晓晓 2026-03-18 12/600 2026-03-21 22:40 by ACS Nano——
[考研] 求助 +5 梦里的无言 2026-03-21 6/300 2026-03-21 17:51 by 学员8dgXkO
[考研] 22 350 本科985求调剂,求老登收留 +3 李轶男003 2026-03-20 3/150 2026-03-21 13:28 by 搏击518
[考研] 求调剂 +6 Mqqqqqq 2026-03-19 6/300 2026-03-21 08:04 by JourneyLucky
[考研] 265求调剂 +3 Jack?k?y 2026-03-17 3/150 2026-03-21 03:17 by JourneyLucky
[考研] 083200学硕321分一志愿暨南大学求调剂 +3 innocenceF 2026-03-17 3/150 2026-03-21 02:35 by JourneyLucky
[考研] 332求调剂 +4 ydfyh 2026-03-17 4/200 2026-03-21 02:20 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +4 葵梓卫队 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 081700化工学硕调剂 +3 【1】 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
信息提示
请填处理意见