24小时热门版块排行榜    

查看: 2544  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 立项公示后,已经有几个老师及同学过来要本子,该怎么应对? +28 ty28ty28 2024-09-15 34/1700 2024-09-21 09:48 by guaiyaya2008
[考博] 浙江大学生物学博导推荐 +4 小小谢尔顿 2024-09-20 4/200 2024-09-21 09:30 by UpwindRose
[论文投稿] 论文投稿 5+3 love319 2024-09-20 5/250 2024-09-21 07:30 by nono2009
[考博] 帮帮我 +5 零差价 2024-09-19 5/250 2024-09-21 01:11 by ertwe1
[论文投稿] PCCP投稿 +4 wuxiaohuang1 2024-09-20 4/200 2024-09-21 01:09 by 褚师宣密
[硕博家园] 感觉在瑞典不适合做科研 +11 毕生所学 2024-09-19 14/700 2024-09-21 01:04 by 褚师宣密
[考博] 2025申博 +3 yqrxxxl 2024-09-18 4/200 2024-09-20 19:49 by ZHI1009
[硕博家园] 关于博士选择 +3 翰墨古风 2024-09-20 6/300 2024-09-20 17:25 by 大泡沫
[考博] 套磁信没有回复 +4 Wicoo 2024-09-18 5/250 2024-09-20 17:10 by No1Maxim
[基金申请] 山东省自然科学基金有查到专家意见么 +12 乱码哈哈哈 2024-09-18 16/800 2024-09-20 10:30 by andywei1028
[基金申请] QB 倒计时 +8 sztai 2024-09-17 10/500 2024-09-20 09:42 by adventure
[论文投稿] Acs nano投稿 +4 一毫米www 2024-09-14 6/300 2024-09-19 22:14 by nic111
[论文投稿] 顶级水刊,几乎不拒稿的SCI! +4 a433587770 2024-09-18 4/200 2024-09-19 21:50 by ndyjjcb
[论文投稿] 作者署名:nature大子刊共一还是共同通讯? 10+7 Eriknpeng 2024-09-16 18/900 2024-09-19 17:32 by Eriknpeng
[论文投稿] 有没有便宜一点的开源期刊 100+3 zhuzhuyaya 2024-09-19 4/200 2024-09-19 16:43 by 莱茵润色
[有机交流] 油状二胺 10+3 cui19236 2024-09-14 4/200 2024-09-19 14:18 by tianjunbo
[有机交流] 酰胺缩合反应 5+3 我努力三百 2024-09-14 3/150 2024-09-19 10:04 by czyzsu
[论文投稿] 期刊主题不符是真的主题不行吗还是内容不行 5+6 AlexaHu 2024-09-15 11/550 2024-09-19 09:40 by AlexaHu
[考博] 求助求助 15+3 恰的嘎嘣脆 2024-09-15 3/150 2024-09-19 09:20 by Kelaizhang
[论文投稿] proof 修改 20+3 力行成士 2024-09-14 9/450 2024-09-15 19:43 by 力行成士
信息提示
请填处理意见