24小时热门版块排行榜    

查看: 674  |  回复: 3
【悬赏金币】回答本帖问题,作者ypzzu将赠送您 200 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ypzzu

新虫 (著名写手)

[求助] lammps中compute viscosity/cos command命令如何使用已有1人参与

求助大神!!!小弟初学lammps,想模拟液体粘度,看到“compute viscosity/cos command”命令可以做相关计算,但是不清楚具体使用方法,请各位大神指点迷津。
下面是in文件
units                                           real                        #单位
dimension                                3                                #维度
boundary                        p p p                #边界类型

pair_style                                lj/cut 10.0
atom_style                      full                                                                #原子类型
bond_style                                 harmonic                                #键类型-弹簧势
angle_style                         harmonic                                #角类型
dihedral_style                 harmonic                                #二面角类型

read_data                                216-PEG200-NVE1.data  #读入data文件
timestep        0.1                                                #时间步长1fs
neighbor       4.0 bin                                #近邻边界2埃
neigh_modify     every 10 delay 0 check yes                #每10步更新一次
thermo                100 #
fix                                                        1 all nve

run              20000

compute  cos all viscosity/cos
variable V equal c_cos[7]
variable A equal 0.02E-5
variable density equal density
variable lz equal lz
variable reciprocalViscosity equal v_V/${A}/v_density*39.4784/v_lz/v_lz*100

fix 2 all accelerate/cos 0.02e-5
run  5000
print            "reciprocalViscosity : ${reciprocalViscosity }"
执行以后提示:
compute  cos all viscosity/cos
ERROR: Unrecognized compute style 'viscosity/cos' (../modify.cpp:1206)
Last command: compute  cos all viscosity/cos
已经安装了官方说明书里提示的“USER-MISC package”
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ypzzu

新虫 (著名写手)

自己顶一下,拜求各位!
3楼2021-06-18 18:59:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ypzzu

新虫 (著名写手)

不要沉
4楼2021-06-20 20:35:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ypzzu 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 请问这审稿意见准确地说是啥意思 100+3 枯禅 2024-05-11 4/200 2024-05-12 11:17 by xli1984
[教师之家] 学生家长私下联系老师修改成绩不成,唆使19名学生联名要求复核成绩 +12 sjtu2012 2024-05-11 14/700 2024-05-12 10:43 by tfang
[基金申请] 特聘教授,算高级职称吗?青年基金还没结题的情况下还能报面上吗? +22 aasahr 2024-05-05 31/1550 2024-05-12 09:00 by feng6531
[教师之家] 相比过去,现在高校导师水平越来越高(可见招聘条件),研究生越来越差 +7 苏东坡二世 2024-05-11 10/500 2024-05-11 21:12 by feng6531
[论文投稿] 投稿竟然倒在了Scientific Reports上。 +7 lizhengke06 2024-05-10 7/350 2024-05-11 20:53 by feng6531
[找工作] 985博士毕业,前几天通过了河南省科学院下面一个所的面试,值得去吗? +12 wendao_2016 2024-05-07 16/800 2024-05-11 20:23 by 无知!无畏?
[硕博家园] 矿大本-中科大硕(推免)-24年科研助理求职-方向不限 +8 Kelaizhang 2024-05-08 16/800 2024-05-11 19:58 by Kelaizhang
[考博] 一篇一作AM,寻求西安985博导 +6 LYWwrz 2024-05-06 7/350 2024-05-11 19:37 by 东昌莘
[基金申请] E05送审了么? +5 wly78528 2024-05-08 5/250 2024-05-11 18:54 by 小博shine
[基金申请] E06送审了吗? +4 edge099 2024-05-08 4/200 2024-05-11 18:52 by pangmingju
[有机交流] 常见的固体有机酸由哪些呢? 40+4 whw818 2024-05-09 8/400 2024-05-11 18:45 by mapenggao
[基金申请] 科研之友老是给我发消息 +4 问君611 2024-05-11 5/250 2024-05-11 16:42 by 问君611
[基金申请] NSFC函评几月几日截止提交评审意见? +7 瞬息宇宙 2024-05-11 7/350 2024-05-11 16:17 by 南海猪
[药学] 本科药学,未来怎么选择 17+3 Roywjq777 2024-05-07 7/350 2024-05-11 15:13 by jy1230910
[基金申请] 是这样的? +17 xiangsu121 2024-05-10 19/950 2024-05-11 11:43 by zqwcr
[硕博家园] 还能保住学位证书吗 +9 屠师! 2024-05-09 14/700 2024-05-11 11:32 by 屠师!
[硕博家园] 盲审 +13 暮色恋伊人 2024-05-07 15/750 2024-05-10 21:28 by 暮色恋伊人
[基金申请] 函评这一关首先干掉一大半,今年卷上天了 (金币+5) +26 mycaas 2024-05-07 36/1800 2024-05-10 17:22 by 降温器
[基金申请] 解某透露,生命学部的项目质量在评审中明显高于医学部 +12 babu2015 2024-05-06 12/600 2024-05-09 10:00 by 星火12
[考博] 24-25申博 +5 muscling 2024-05-06 10/500 2024-05-08 21:33 by lizhengke06
信息提示
请填处理意见