24小时热门版块排行榜    

查看: 3211  |  回复: 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的回帖

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的回帖
查看全部 8 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 8/400 2026-09-01 12:58 by ssyjh
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 15/750 2026-09-01 12:23 by icm639
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 投票:  有多少人是今天查系统知道结果的? +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
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +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
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 出来了 +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
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
信息提示
请填处理意见