24小时热门版块排行榜

>论坛更新日志 (3680)
>虫友互识 (1169)
>导师招生 (305)
>论文道贺祈福 (184)
>文献求助 (157)
>休闲灌水 (87)
>高分子 (80)
>考博 (69)
>硕博家园 (67)
>学术会议 (56)
>论文投稿 (44)
>招聘信息布告栏 (35)
>博后之家 (35)
>找工作 (34)
>考研 (28)
>教师之家 (18)

返回列表

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

新春愉快～

回复此楼