24小时热门版块排行榜    

查看: 10673  |  回复: 5

放逐孤岛

铜虫 (小有名气)

[求助] Non-numeric pressure - simulation unstable

我在使用lammps读入我刚刚写好的restart文件时,产生了这个错误:
ERROR: Non-numeric pressure - simulation unstable (../fix_nh.cpp:1038)
这会是什么原因呢,请大家帮我看一下
我的in文件:
read_restart            restart1.file

neighbor                2 bin
neigh_modify         delay 0 every 10 check yes

dump                     1 all atom 1000 dump.lammpstrj

thermo_style          custom step temp pe etotal press vol
thermo                   100

fix                           1 all npt temp 116.82 116.82 $(100.0*dt) iso 20 -400 1.0
timestep                 1
run                         1

fix                          1 all npt temp 116.82 116.82 $(100.0*dt) iso -400 -400 25.0
timestep                1
run                         100000

write_restart           "restart2.file"

链接: https://pan.baidu.com/s/1dyORn5KOWUOj2t5vugwI8g 提取码: ayxu
这个是源文件的地址
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

放逐孤岛

铜虫 (小有名气)

我找到了一部分原因,是我的数值发散导致了这个错误。我修改了参数后可以计算了,然后计算了部分数据,结果如下:
Per MPI rank memory allocation (min/avg/max) = 30.81 | 30.82 | 30.82 Mbytes
Step Temp PotEng TotEng Press Volume
   10001    116.05116   -451282.31   -278367.17   -71.398375     34852070
   10100    111.01823    -437462.8   -272046.68   -198.26378     36272261
   10200    108.08415   -426547.22   -265502.85   -198.13712     37045806
   10300    106.94912   -417128.38   -257775.19   -198.04581     37836073
   10400    106.46244   -407738.39   -249110.34   -197.62584     38697213
   10500    106.12565   -397667.18   -239540.94   -197.21792     39677059
   10600    106.07937    -387294.9   -229237.62   -196.68675     40766912
   10700     105.9978   -376011.89   -218076.15   -195.98214     42027638
   10800    105.86741   -363704.99   -205963.54    -194.9769     43503980
   10900    105.65782   -350145.11   -192715.94   -193.50238     45264696
   11000    105.04697   -334378.57   -177859.56    -191.1722     47464853
   11100     104.1841   -316208.34      -160975   -187.49958     50255070
   11200    103.04179   -294870.36   -141339.05   -181.61763     53928105
   11300    101.50664   -269045.43   -117801.47   -171.34534     59018740
   11400    99.424484   -236820.83   -88679.272   -153.10299     66568759
   11500    96.194298   -194718.13   -51389.526   -121.85042     79065126
   11600    88.997703    -133705.7   -1099.9542   -70.232476 1.0565983e+08
   11700    62.876195   -38285.259    55399.686   -2.5910112 2.2204799e+08
ERROR: Non-numeric pressure - simulation unstable (../fix_nh.cpp:1038)
Last command: run                     100000

明显数据发散的很严重,所以我猜测是数据发散引起的原因。
但是我现在的问题是出现这个问题后我该如何处理呢,请各位提供一些建议。
2楼2019-03-28 16:43:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独往烟火

铁虫 (初入文坛)

您好,请问您解决这个问题了吗?您修改了哪些参数呢?
3楼2019-11-05 20:18:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

TDTR

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 放逐孤岛 at 2019-03-28 16:43:12
我找到了一部分原因,是我的数值发散导致了这个错误。我修改了参数后可以计算了,然后计算了部分数据,结果如下:
Per MPI rank memory allocation (min/avg/max) = 30.81 | 30.82 | 30.82 Mbytes
Step Temp PotE ...

您好,请问您是如何处理的?您修改了哪些参数呢
4楼2021-07-26 18:34:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

搬砖的狗子

新虫 (初入文坛)

5楼2021-11-21 21:24:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

搬砖的狗子

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by TDTR at 2021-07-26 18:34:54
您好,请问您是如何处理的?您修改了哪些参数呢...

请问解决了吗?

发自小木虫IOS客户端
6楼2021-11-21 21:24:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 放逐孤岛 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 让我中一个面上吧! +12 大萍1987 2026-08-20 14/700 2026-08-23 10:30 by wrm
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 07:41 by OEbVnUOu01ol
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 10/500 2026-08-23 07:17 by OEbVnUOu01ol
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[基金申请] 范进中举一文的中心思想 +3 zjuhero 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见