24小时热门版块排行榜

返回列表

yyfdemajia

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 2826.7
散金: 1100
红花: 2
沙发: 2
帖子: 858
在线: 113.7小时
虫号: 3295617
注册: 2014-06-27
专业: 胶体与界面化学

[求助] lammps 流體中移動的問題已有1人参与

我剛接觸lammps半個月，想要模擬孔在流體作用下剪切膨脹的過程，通過兩塊平行的固體平面的垂直移動?砟M孔的膨脹過程。我參考了lammps中couette flow 的例子，設置了移動的速度，問題是進行模擬時，上平面卻無法移動。我把lammps的代碼貼在下面，請問我是哪裏出現的問題？？？

boundary  p p p
dimension  3
units  lj

atom_style  atomic
lattice fcc 1.6
region box block 0 16 1 13.8 0 16.8

create_box 2 box

region upper block 10 11.6 1 13.8 2 14.8
region lower block 4.4 6 1 13.8 2 14.8

#create_atoms 1 random 80 1 box
create_atoms 1 region box
delete_atoms region upper
delete_atoms region lower
create_atoms 2 region upper
create_atoms 2 region lower

group  upper region upper
group  lower region lower
group  mobile region box
group  solid union lower upper
group  fluid subtract all solid

pair_style lj/cut 3.5
pair_coeff * * 1 1 3.5

mass 1 1.0
mass 2 1.0

set group upper type 2
set group lower type 2
set group fluid type 1

compute       mobile fluid temp
velocity    fluid create 1.0 482748 temp mobile
velocity    lower set 0.0 0.0 0.0
velocity    upper set 2.0 0.0 0.0

fix 1 fluid gcmc 1 10 11440 1 1234 1 -2.5 0.01
#fix 2 upper move linear NULL NULL NULL
#fix 3 lower move linear NULL NULL NULL
#fix 4 solid rigid/nvt group 2 upper lower temp 1 1 0.1

thermo 10
thermo_style multi temp etotal ethalpy
compute 11 upper stress/atom NULL

minimize 0.0 1.0e-8 1000 100000
dump box all xyz 1 sample.xyz

timestep 0.01
run 300

新春愉快～

回复此楼

» 收录本帖的淘帖专辑推荐

资源收集

» 猜你喜欢

售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急已经有6人回复
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急已经有5人回复
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急已经有5人回复
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急已经有5人回复
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急已经有3人回复
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急已经有5人回复
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有6人回复
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有5人回复
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急已经有6人回复
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急已经有5人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2015-02-19 14:29:04

已阅关注TA 给TA发消息送TA红花 TA的回帖

abinitio

版主 (著名写手)

模拟EPI: 2
应助: 221 (大学生)
贵宾: 0.093
金币: 7365.2
散金: 3917
红花: 57
帖子: 1317
在线: 485.2小时
虫号: 368545
注册: 2007-05-11
性别: GG
专业: 凝聚态物性I:结构、力学和
管辖: 分子模拟

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
月只蓝: 金币+10, 鼓励交流！春节快乐！ 2015-02-19 23:30:00
yyfdemajia(月只蓝代发): 金币+5, 感谢指导 2015-02-20 09:14:39

是不是boundary出了问题？把z方向上改为非周期性边界条件s，如何？新春快乐！

[ 发自小木虫客户端 ]

赞一下