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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 冲击结果异常
121/2返回列表12下一页
查看: 2013  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 梦想 +5 myrtle 2026-08-26 8/400 2026-09-02 02:42 by Leogzhya
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +17 爱看书的可乐 2026-08-26 19/950 2026-09-02 00:19 by xiangy672
[基金申请] 为什么国自然不能直接公布 +5 bjdxyxy 2026-08-26 5/250 2026-09-01 17:54 by 小伟大博士
[基金申请] 学科评审组评审是指会评吗？ +4 瞬息宇宙 2026-08-31 4/200 2026-09-01 14:58 by jiaoxg
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 基金未中，这种答复是模板吗？ +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 能否申诉？ +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 为什么到现在没收到通知？ +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
关闭
请填处理意见
关闭