24小时热门版块排行榜    

查看: 3297  |  回复: 7

zhang_jaj

木虫 (正式写手)

[求助] 关于fix wall/region命令的问题

大家好,对于fix wall/region 命令里面有这样一句:Regions are defined using the [url=]region[/url] command.  Note thatthe region volume can be interior or exterior to the bounding surface,which will determine in which direction the surface interacts withparticles, i.e. the direction of the surface normal.

关于这句话大家怎么理解,是认为一个region的里侧和外侧均可当做作用面吗?

比如一般定义region wall-1 block 0 1 INF INF 0 1 和region wall-2 block 0 1 INF INF 0 1 side out 两个region,关于wall-2的region我用fix wall/region命令,这样得到的作用面是这个wall-2的里侧面呢还是wall-1的外侧面呢?如果我想得到wall-1这个region的外侧面的作用面时该怎么定义呢?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ddll_slsl

木虫 (正式写手)

【答案】应助回帖

★
感谢参与,应助指数 +1
jiaoyixiong: 金币+1, 赞,很有道理 2013-01-17 10:38:40
具体情况,你跑一下程序,看看轨迹结果,不就全明白了吗?
2楼2013-01-17 09:54:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhang_jaj

木虫 (正式写手)

引用回帖:
2楼: Originally posted by ddll_slsl at 2013-01-17 09:54:30
具体情况,你跑一下程序,看看轨迹结果,不就全明白了吗?

嗯,跑了,但是加的lj93 wall竟然对particles没有吸引作用,所以才看不出来,谢谢
3楼2013-01-18 18:16:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ddll_slsl

木虫 (正式写手)

【答案】应助回帖

★
zh1987hs: 金币+1, 谢谢 2013-01-19 09:16:37
引用回帖:
3楼: Originally posted by zhang_jaj at 2013-01-18 18:16:23
嗯,跑了,但是加的lj93 wall竟然对particles没有吸引作用,所以才看不出来,谢谢...

dump出原子的作用力,看看问题在哪里。 需要你耐心地解决问题。
4楼2013-01-18 21:36:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sdzjc1988

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
zhang_jaj: 金币+5, ★★★很有帮助 2013-01-19 08:57:25
zh1987hs: 金币+2, 谢谢 2013-01-19 09:16:46
这个命令我用过,region wall-1 block 0 1 INF INF 0 1就是里面的墙壁有作用力,region wall-2 block 0 1 INF INF 0 1 side out是外面的墙壁有作用力,你说没有吸引是不是因为原子不在cutoff里面。

» 本帖已获得的红花(最新10朵)

5楼2013-01-19 06:03:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhang_jaj

木虫 (正式写手)

送鲜花一朵
引用回帖:
5楼: Originally posted by sdzjc1988 at 2013-01-19 06:03:57
这个命令我用过,region wall-1 block 0 1 INF INF 0 1就是里面的墙壁有作用力,region wall-2 block 0 1 INF INF 0 1 side out是外面的墙壁有作用力,你说没有吸引是不是因为原子不在cutoff里面。

非常感谢,不可能不在cutoff里面。我设置的cutoff是13.0,有几个atom的距离到这个外壁都是3(lattice);关于这里我想问一下:units我用的是real,但是下面的box,region的units都用的是默认的lattice单位,而势函数的distance units应该是real units一致的吧(即13.0)。那这样单位岂不是不统一了吗?不用自己在转化吧,谢谢?
6楼2013-01-19 08:57:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhang_jaj

木虫 (正式写手)

引用回帖:
4楼: Originally posted by ddll_slsl at 2013-01-18 21:36:10
dump出原子的作用力,看看问题在哪里。 需要你耐心地解决问题。...

嗯,谢谢哦。纠结了很久了。
7楼2013-01-19 08:58:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhang_jaj

木虫 (正式写手)

送鲜花一朵
引用回帖:
5楼: Originally posted by sdzjc1988 at 2013-01-19 06:03:57
这个命令我用过,region wall-1 block 0 1 INF INF 0 1就是里面的墙壁有作用力,region wall-2 block 0 1 INF INF 0 1 side out是外面的墙壁有作用力,你说没有吸引是不是因为原子不在cutoff里面。

你好,对于region命令来说,wall-1和wall-2是两个单独的region。如果对这两个用fix wall/region的话,那么第一个就是wall-1的内壁作用,第二个就是wall-1的外壁作用而不是wall-2的内壁作用?
我这个in文件麻烦您帮忙看一下,为什么atom就是没有感觉到有吸引力,谢谢啊
units                        real
atom_style                atomic
boundary                p p ff

lattice                        sc 4.0
region                        box block 0 10 0 10 0 10
create_box                1 box
region                        atom-box block 2 8 2 8 2 8
create_atoms                1 region atom-box  

group                        ar type 1
mass                        1 39.948

velocity                all create 100.0 8765788
timestep                20.0

pair_style                lj/cut 13.0
pair_coeff                * * 0.2381 3.405

neighbor                2.0 bin
neigh_modify                every 10 delay 0 check yes

thermo                        500
fix                        zlowall all wall/lj93 zlo EDGE 0.5 3.405 13.0      
fix                        ref all wall/reflect zhi EDGE

fix                        1 all nvt temp 90.0 90.0 2.0
dump                        1 all cfg 500 ar.*.cfg id type xs ys zs vx vy vz
dump                        2 all xyz 1000 ar.xyz
run                        100000
8楼2013-01-21 14:28:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhang_jaj 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +3 wangyl_123 2026-10-10 6/300 2026-10-11 15:04 by Abin0406
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[教师之家] 课题组会帮你通过考核么? +3 好运来~9527 2026-10-10 5/250 2026-10-11 09:00 by 好运来~9527
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +5 13108017953 2026-10-08 5/250 2026-10-10 11:15 by wdwd123
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见