24小时热门版块排行榜    

查看: 1195  |  回复: 3

Bessel

木虫 (正式写手)

[求助] 请谁能解释一下Lammps里removing a "bias" from the atom velocities.

看了半天说明和其他文献,百思不得其解。
fix langevin command :

like other fixes that perform thermostatting, this fix can be used with compute commands that remove a "bias" from the atom velocities. E.g. removing the center-of-mass velocity from a group of atoms or removing the x-component of velocity from the calculation. This is not done by default, but only if the fix_modify command is used to assign a temperature compute to this fix that includes such a bias term. See the doc pages for individual compute commands to determine which ones include a bias. In this case, the thermostat works in the following manner: bias is removed from each atom, thermostatting is performed on the remaining thermal degrees of freedom, and the bias is added back in.

是不是说,比如我有一个沿X方向的couette flow,应为摩擦体系温度要升高,为了保持体系的温度,施加了langevin恒温器,但是在计算体系的温度时,要把x方向的动能刨除?不知道我理解的对不对?

关键问题是,“bias”到底是什么?我一直不理解。

谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sdzjc1988

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
zh1987hs: 金币+1, 谢谢 2012-10-03 08:17:25
Bessel: 金币+10, ★★★★★最佳答案 2012-10-04 15:25:29
算温度用的是能量均分公式:3/2kbT=1/2mv2,减去bias就是把v里面的宏观速度减去,否则的话系统宏观运动的v也会用来算温度。
2楼2012-10-03 05:38:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Bessel

木虫 (正式写手)

引用回帖:
2楼: Originally posted by sdzjc1988 at 2012-10-02 22:38:58
算温度用的是能量均分公式:3/2kbT=1/2mv2,减去bias就是把v里面的宏观速度减去,否则的话系统宏观运动的v也会用来算温度。

宏观速度 is the streaming velocity,  Am I right?

sorry, I cannot input Chinese.
3楼2012-10-03 14:00:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sdzjc1988

木虫 (小有名气)

【答案】应助回帖


zh1987hs: 金币+1, 谢谢 2012-10-03 22:49:18
引用回帖:
3楼: Originally posted by Bessel at 2012-10-03 14:00:41
宏观速度 is the streaming velocity,  Am I right?

sorry, I cannot input Chinese....

这要看你的steam veloctiy是不是系统的宏观速度了,这个应该比较好判断~
4楼2012-10-03 22:46:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Bessel 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 wsYGJnURYvm6 2026-08-31 3/150 2026-09-01 01:52 by o3EMxHUAQDtS
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 wsYGJnURYvm6 2026-08-31 3/150 2026-09-01 01:40 by o3EMxHUAQDtS
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 EZ67DG6CvcKL 2026-08-31 4/200 2026-09-01 01:30 by o3EMxHUAQDtS
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 为什么到现在没收到通知? +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
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国合里面能看到了 +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
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见