24小时热门版块排行榜    

查看: 1974  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

奈落之夜

木虫 (正式写手)

[求助] 【求助】lammps关于fix ave/spatial命令应用问题

想利用fix ave/spatial的分层和求平均来做一个热传导,传导后希望能够得到每一层的温度能量来画图,如何能输出每一层的能量啊 。可有例子。用compute 1 all ke/atom 或者计算temp 添加在后面 能否输出每一层的能量和温度

[ Last edited by 奈落之夜 on 2014-2-18 at 17:50 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhang53201

铜虫 (小有名气)


月只蓝: 金币+1, 鼓励交流! 2014-02-28 08:14:48
引用回帖:
3楼: Originally posted by 奈落之夜 at 2014-02-26 20:23:53
多谢 我会注意原子数目的。 有几个问题请问一下,
问题1:compute设置成特定的那块层,是说要在前面设置compute的时候就设定吗,这个命令会不会自己输出特定层的。也就是说我再前面设定compute ke/atom,例如下面的 ...

这么多问题啊。我说下问题一吧,第一种方法是,compute的时候,不要写all,写group1,group2,等。
第二种,就是写个脚本,按照你几何形状的左边位置,对输出的原子提取一下。

问题二是,什么norm。。?
4楼2014-02-27 20:57:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

zhang53201

铜虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★
月只蓝: 金币+1, 鼓励交流! 2014-02-25 21:45:56
奈落之夜: 金币+5, ★★★很有帮助 2014-02-26 20:25:57
温度是统计平均的值,只有原子数目多了,才有意义。不过你可以对region设置很多小的片层,用compute 计算每个片层的温度,但这个厚度不能太少。compute 那里就写成特定的那块层。
2楼2014-02-25 20:38:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

奈落之夜

木虫 (正式写手)

★ ★
月只蓝: 金币+2, 鼓励交流! 2014-02-26 21:59:04
引用回帖:
2楼: Originally posted by zhang53201 at 2014-02-25 20:38:49
温度是统计平均的值,只有原子数目多了,才有意义。不过你可以对region设置很多小的片层,用compute 计算每个片层的温度,但这个厚度不能太少。compute 那里就写成特定的那块层。

多谢 我会注意原子数目的。 有几个问题请问一下,
问题1:compute设置成特定的那块层,是说要在前面设置compute的时候就设定吗,这个命令会不会自己输出特定层的。也就是说我再前面设定compute ke/atom,例如下面的例子compute           ke all ke/atom
compute           ke all ke/atom
variable          temp atom c_ke/1.5/8.617343e-5
variable          x  atom  x
variable          y  atom  y
variable          z  atom  z
fix               temp_profile_graphene   mobile     ave/spatial  1  100 1000  z  lower  ${delta_CNT}  c_ke v_z  v_temp  file  temp_profile_graphene_hss2.dat   units  reduced
定义的时候只是定义了计算每个原子的动能之类。输出的时候能不能对应于每个层输出相应层的数据.

问题2:我建的就是石墨烯 不太规则所以lattice不知道如何定义。如果用reduced应该不用设置lattice。假设在x方向分25层是不是就应该定义为fix               temp_profile_graphene   mobile     ave/spatial  1  100 1000  x  lower  0.04  c_ke v_z  v_temp  file  temp_profile_graphene_hss2.dat   units  reduced。
norm这个命令用不用加上,
问题3:有解释说norm关键字影响到每Nfreq步的输出结果如何平均出来。对于all选项,层属性对所有Nrepeat个采样点的所有原子进行平均,输出值是总量和/总原子数。是不是加上norm all之后算的值才是整个层所有原子平均的温度或者能量的。
问题:4:里面的ave 说设为one为每个 输出步 都是独立的输出结果。而runing则输出步的计算还要与以前输出的计算值进行平均。这是如何与以前的计算进行平均。例如上面设置的 每1000步 把前面100步的结果进行计算。这个如果设置了ave runing或者one会有什么不同呢。
感觉是自己知识不够,看说明书,看别人的翻译看不太懂。希望能帮我解答一下
3楼2014-02-26 20:23:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

奈落之夜

木虫 (正式写手)

引用回帖:
4楼: Originally posted by zhang53201 at 2014-02-27 20:57:06
这么多问题啊。我说下问题一吧,第一种方法是,compute的时候,不要写all,写group1,group2,等。
第二种,就是写个脚本,按照你几何形状的左边位置,对输出的原子提取一下。

问题二是,什么norm。。?...

晕,不好意思这么久才回复你。第一个问题就是group只能设置21组 所以才不能用grouop的  这个命令感觉解释像能直接输出每个层而且每个层都可以输出局部属性才这样的。第二个问题和第三个是一起看的我弄错了。不过这边不让再次修改了。这个fix   ave/spatial   后面有个小参数就是norm,我后来看了下 这个应该设置成all 这样就是所有的原子的平均。一般人都不写其实默认的也是all。就是第一个问题我不太确定是否能够像我想要的那样输出每个层的数据。
5楼2014-03-04 08:19:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 广西大学材料导师推荐 +3 夏夏夏小正 2026-03-17 5/250 2026-03-21 22:20 by 金昊ML
[考研] 333求调剂 +5 87639 2026-03-21 7/350 2026-03-21 19:31 by ColorlessPI
[考研] 考研化学学硕调剂,一志愿985 +5 张vvvv 2026-03-15 7/350 2026-03-21 19:23 by ColorlessPI
[考研] 0703化学调剂 +11 妮妮ninicgb 2026-03-15 15/750 2026-03-21 19:15 by ColorlessPI
[基金申请] 山东省面上项目限额评审 +3 石瑞0426 2026-03-19 3/150 2026-03-21 18:59 by RRRKKK
[考研] 277材料科学与工程080500求调剂 +6 自由煎饼果子 2026-03-16 6/300 2026-03-21 17:21 by 学员8dgXkO
[基金申请] 学校已经提交到NSFC,还能修改吗? 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 0856材料专硕353求调剂 +3 NIFFFfff 2026-03-20 3/150 2026-03-21 10:23 by luoyongfeng
[考研] 324分 085600材料化工求调剂 +4 llllkkkhh 2026-03-18 4/200 2026-03-21 01:24 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-18 3/150 2026-03-21 00:46 by JourneyLucky
[考研] 材料专业求调剂 +6 hanamiko 2026-03-18 6/300 2026-03-21 00:24 by JourneyLucky
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +4 葵梓卫队 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[考研] 考研调剂求学校推荐 +3 伯乐29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
信息提示
请填处理意见