24小时热门版块排行榜    

查看: 3169  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +15 爱看书的可乐 2026-08-26 17/850 2026-08-27 07:31 by xxniao123
[基金申请] 国自然面上复盘~欢迎讨论 +6 晴天加油 2026-08-26 7/350 2026-08-27 07:30 by pinghe007
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 +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
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见