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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » fix viscosity命令求助
51/1返回列表
查看: 1646  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

tsadaye

金虫 (初入文坛)

[求助] fix viscosity命令求助

刚刚开始学习lammps，想先试着用dpd方法计算一个单一珠子体系的粘度，其中粘度由fix viscosity命令计算。写的in文件如下：
dimension 3
units lj
atom_style atomic

read_data data.initial

neighbor 1.2 bin
neigh_modify delay 3

pair_style dpd 1.0 1.0 1234567
pair_coeff 1 1 25.0 4.5

communicate single vel yes

velocity all create 1.0 4928259 dist gaussian
fix 1 all nve
fix 2 all viscosity 50 x z 10

timestep 0.05
thermo 10

dump 1 all local 1000 initial.dump f_2
run 30000
write_restart initial.restart

问题是fix viscosity命令将计算一个粘度，但这个粘度应该由什么命令读取出来？我试着用dump local，但是貌似这个是个group attribute，读不出来。请高人解惑，谢谢。 第一次写in文件，难免很多错误，请轻喷
回复此楼

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

动力学或热力学计算

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-05-25 12:22:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

咕噜咕噜7631

新虫 (初入文坛)
这个问题楼主解决了吗？我也是刚入门，也想知道剪切率到底要怎么算呢
赞一下 回复此楼
5楼2016-11-02 11:24:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

WJ-MS

新虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
tsadaye: 金币+10, ★★★很有帮助 2013-06-19 22:34:10
楼主manual没有仔细看呀，这个命令还算不出粘度呢，只是得到体系的总的动量交换值，除以时间和模拟盒子的横截面积得到动量通量，再除以什么velocity profile的斜率就得到粘度了，这个命令的计算值可以用print f_2命令在屏幕上显示出来，然后在log文件里面也能找到了。我目前也刚弄出这些，对于velocity profile一筹莫展啊，希望能互相探讨探讨~~~
赞一下(1人) 回复此楼
加油！
2楼2013-05-28 16:27:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lixiaojun305

新虫 (初入文坛)
同求啊   楼主要是有进展  可否继续完善此贴？
赞一下 回复此楼
3楼2013-06-25 15:23:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

txcokokok

木虫 (小有名气)
直接看example/viscosity就好了
赞一下 回复此楼
bigbang
4楼2015-07-31 21:57:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 336化工调剂 +3 王大坦1 2026-03-23 4/200 2026-03-23 11:07 by laoshidan
[考研] 298求调剂 +8 上岸6666@ 2026-03-20 8/400 2026-03-23 11:02 by laoshidan
[考研] 085600材料与化工 +10 安全上岸！ 2026-03-16 10/500 2026-03-23 09:31 by Jy?
[考研] 310求调剂 +4 baibai1314 2026-03-16 4/200 2026-03-22 20:19 by edmund7
[考研] 一志愿上海交大生物与医药专硕324分，求调剂 +3 jiajunX 2026-03-22 3/150 2026-03-22 19:32 by brblmd
[考研] 寻找调剂 +4 倔强芒? 2026-03-21 4/200 2026-03-22 16:14 by 木托莫露露
[考研] 一志愿中南化学（0703）总分337求调剂 +9 niko- 2026-03-19 10/500 2026-03-22 16:08 by ColorlessPI
[考研] 289求调剂 +7 怀瑾握瑜l 2026-03-20 7/350 2026-03-22 15:57 by ColorlessPI
[考研] 求调剂院校信息 +6 CX 330 2026-03-21 6/300 2026-03-22 15:25 by 无懈可击111
[考研] 能源材料化学课题组招收硕士研究生8-10名 +5 脱颖而出 2026-03-16 17/850 2026-03-22 15:18 by 脱颖而出
[考研] 求调剂 +5 Zhangbod 2026-03-21 7/350 2026-03-22 13:13 by Zhangbod
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 299求调剂 +4 某某某某位 2026-03-21 4/200 2026-03-21 16:30 by barlinike
[基金申请] 学校已经提交到NSFC，还能修改吗？ 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 294求调剂材料与化工专硕 +15 陌の森林 2026-03-18 15/750 2026-03-20 23:28 by JourneyLucky
[考研] 一志愿 南京航空航天大学大学 ，080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[考研] 0856调剂，是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
[考研] 298-一志愿中国农业大学-求调剂 +9 手机用户 2026-03-17 9/450 2026-03-20 14:24 by 无懈可击111
[考研] 一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 一志愿南京大学，080500材料科学与工程，调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
信息提示
关闭
请填处理意见
关闭