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

新春愉快～

回复此楼