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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 冲击结果异常
121/2返回列表12下一页
查看: 1855  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

lchunmei001

木虫 (正式写手)

[求助] 冲击结果异常

请教大家，我沿x方向给粒子2KM/s的速度冲击加载，具体命令如下，但得到的结果很奇怪，压强很低的情况下，温度就已达到几万K，显然很不合理，但是在不知道是哪里出问题了，测试了几种物质，均出现这种情况，应该是命令哪里没对，还请高手指点一二，感激不尽！
velocity        all set 20  NULL NULL  sum yes units box
fix             1 all nve
fix             2 all wall/reflect xlo xhi
回复此楼

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

condensed matter physics

» 猜你喜欢
1楼 2013-01-16 10:46:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ddll_slsl

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
jiaoyixiong: 金币+2, 鼓励交流 2013-01-16 14:29:01
把整个in文件，和相关输出的log贴出来，否则无法判断。
赞一下 回复此楼
2楼2013-01-16 11:11:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lchunmei001

木虫 (正式写手)
输入和结果如下，请予以指点，谢谢！
输入：
units                 metal
boundary        s p p
lattice         fcc 4.07  
atom_style        atomic
region                box block 0 100 0 10 0 10
create_box        1 box
create_atoms        1 box

pair_style eam
pair_coeff * * Pt_u3.eam

timestep        1.0e-3
neighbor        1.0 bin
neigh_modify    every 1 delay 10 check yes

thermo          100
thermo_style custom step press vol temp pxx pyy pzz lx ly lz

velocity        all set 20  NULL NULL  sum yes units box
fix             1 all nve
fix             2 all wall/reflect xlo xhi

dump                id all atom 10000   dump.melt
run             50000

结果log文件：

Step Press Volume Temp Pxx Pyy Pzz Lx Ly Lz
       0    19146.265    674326.27    31285.986    535234.08   -238897.64   -238897.64     407.0814         40.7         40.7
     100    23424.423    673363.58    31311.738     540451.5   -235089.11   -235089.11    406.50024         40.7         40.7
     200    27653.106    671960.25    31331.714    546177.73   -231609.21   -231609.21    405.65307         40.7         40.7
     300    31059.835    671084.06    31347.923    550838.79   -228829.64   -228829.64    405.12413         40.7         40.7
     400    36428.574    669579.59    31362.412    557393.91   -224054.09   -224054.09    404.21589         40.7         40.7
     500    39253.628    668750.35    31377.061    561471.39   -221855.25   -221855.25     403.7153         40.7         40.7
     600    44925.431     667277.7    31394.038    568443.51   -216833.61   -216833.61    402.82628         40.7         40.7
     700    47685.276    666310.63    31413.841     572717.7   -214830.93   -214830.93    402.24247         40.7         40.7
。。。。。。
赞一下 回复此楼
3楼2013-01-16 11:36:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ddll_slsl

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
jiaoyixiong: 金币+3, 赞！鼓励交流 2013-01-16 15:10:23
lchunmei001: 金币+5 2013-01-16 16:23:08
引用回帖:
3楼: Originally posted by lchunmei001 at 2013-01-16 11:36:01
输入和结果如下，请予以指点，谢谢！
输入：
units                 metal
boundary        s p p
lattice         fcc 4.07  
atom_style        atomic
region                box block 0 100 0 10 0 10
create_box        1 box
create_atoms        1 box

p ...

你的原子初速度很高，原子从边界反弹会与后面的原子频繁撞击。 X方向的压力很高， 有可能产生瞬间高温。 看能量，你的过程还算平稳。 有几点可以尝试：
1）你用dump命令输出各个原子的速度值，看看其最大值是多少？
2） 盒子的X尺寸在变小，是你所需要的吗？ 边界用p p p试试？
赞一下 回复此楼
4楼2013-01-16 14:55:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
lchunmei001: 金币+5 2013-01-16 16:23:13
zh1987hs: 金币+3, 谢谢 2013-01-19 14:04:12
楼主，我有点儿不懂，请教你哈：
1、你是要模拟一个什么样的冲击呢？我把你的in文件运行了一下，4万多个原子都一起向x方向运动吗？安装我的理解，所谓的“冲击”应该是一个物体以较高的速度撞击另一个物体，比如炮弹攻击掩体的时候。我看过使用lammps模拟纳米颗粒冲击石墨烯的。等等这些。我没弄明白你的模拟要干啥。
2、在你的in文件中，你能告诉我哪句命令是thermostat的？你设置的温度是多少？


PS 把dump命令改成  dump                id all atom 10000   dump.lammpstrj  吧，这样方便模拟完了使用VMD观看轨迹动画。
赞一下 回复此楼
5楼2013-01-16 15:20:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lchunmei001

木虫 (正式写手)
引用回帖:
4楼: Originally posted by ddll_slsl at 2013-01-16 14:55:28
你的原子初速度很高，原子从边界反弹会与后面的原子频繁撞击。 X方向的压力很高， 有可能产生瞬间高温。 看能量，你的过程还算平稳。 有几点可以尝试：
1）你用dump命令输出各个原子的速度值，看看其最大值是多少 ...

各个原子的速度差距不大的，1.8几-2.1几之间波动，wall/reflect只针对非周期性边界条件，x方向的尺寸肯定会变
赞一下 回复此楼
6楼2013-01-16 15:30:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lchunmei001

木虫 (正式写手)
引用回帖:
5楼: Originally posted by jiaoyixiong at 2013-01-16 15:20:55
楼主，我有点儿不懂，请教你哈：
1、你是要模拟一个什么样的冲击呢？我把你的in文件运行了一下，4万多个原子都一起向x方向运动吗？安装我的理解，所谓的“冲击”应该是一个物体以较高的速度撞击另一个物体，比如炮 ...

谢谢您的建议。另外这个是采用动量镜像法，让这样一个体系撞向无穷大的墙，所以所有的原子都一起向一个方向运动
赞一下 回复此楼
7楼2013-01-16 15:34:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ddll_slsl

木虫 (正式写手)

【答案】应助回帖

★
jiaoyixiong: 金币+1, 鼓励交流 2013-01-16 21:23:52
引用回帖:
6楼: Originally posted by lchunmei001 at 2013-01-16 15:30:24
各个原子的速度差距不大的，1.8几-2.1几之间波动，wall/reflect只针对非周期性边界条件，x方向的尺寸肯定会变...

如果这样的话，温度似乎可以不那么在意。 看看动态轨迹是否合理呢？
赞一下 回复此楼
8楼2013-01-16 15:56:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lchunmei001

木虫 (正式写手)
我需要温度值啊!而且这个结果确实是很不合理的。我想知道是不是我的输入哪里有问题？
赞一下 回复此楼
9楼2013-01-16 16:12:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ldy85367996
10楼2013-01-16 16:41:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lchunmei001 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 4/200 2026-05-15 17:43 by x0mp7owy2b
[教师之家] 上海大学实验技术岗位非升即走 +3 嘻嘻哈哈乐呵呵 2026-05-15 3/150 2026-05-15 15:40 by 187385hgh
[论文投稿] 有带发论文的吗 +4 山楂之术 2026-05-09 4/200 2026-05-15 15:40 by 妹子不好惹
[基金申请] 青C资助名额大幅增加！ +11 西葫芦炒鸡蛋 2026-05-13 15/750 2026-05-15 14:36 by ambravo
[公派出国] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 3/150 2026-05-15 12:03 by onwj4wpxp2
[基金申请] 精华III:  评审感受-评审感受-评审感受 +14 ferrarichen 2026-05-11 18/900 2026-05-15 11:12 by cmhchen
[基金申请] 这年头没有找到涵评专家，还有中面上的可能吗 +9 dd921ww 2026-05-12 10/500 2026-05-15 10:41 by muyiliuhui
[考博] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 3/150 2026-05-15 09:23 by onwj4wpxp2
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 cjf4bx70cj 2026-05-14 4/200 2026-05-15 09:03 by gagyerk94e
[教师之家] 教学课件你会给同学吗 +8 硕士研究生吗 2026-05-13 8/400 2026-05-14 22:23 by 常规沥青
[考博] 26应届毕业生考博求助 +3 wo一定上岸 2026-05-13 3/150 2026-05-14 21:47 by 明海天涯
[基金申请] 重磅！青年科学基金项目（C类）资助增幅预计超过50% +5 水和泥不是水泥 2026-05-13 7/350 2026-05-14 20:57 by 水和泥不是水泥
[有机交流] 求助2，4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 5/250 2026-05-14 20:15 by 一切都是空工
[高分子] 本人最近太闲了，谁有问题可以提，每天会统一回复 +8 一切都是空工 2026-05-12 19/950 2026-05-14 20:03 by 一切都是空工
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
[基金申请] 请问大佬b0816评完了吗 +3 市民华南虎 2026-05-12 7/350 2026-05-14 07:41 by 市民华南虎
[硕博家园] 导师各种操作恶心咋办 +11 苍白的小青天 2026-05-09 13/650 2026-05-13 17:11 by 六两废铜
[考博] 西南大学考核制博士 +3 lijunjie84 2026-05-11 6/300 2026-05-12 18:09 by lijunjie84
[文学芳草园] 窗边初夏的小雨 +7 阿美_Lml888 2026-05-09 10/500 2026-05-12 15:27 by 阿美_Lml888
[考博] 现在不知道怎么办，感觉很痛苦 +4 qweww 2026-05-11 5/250 2026-05-11 20:23 by Oversize
信息提示
关闭
请填处理意见
关闭