git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@23 f3b2605a-c512-4ea7-a41b-209d697bcdaa

examples/obstacle/in.obstacle

@@ -0,0 +1,74 @@
+# 2d LJ obstacle flow
+
+dimension	2
+boundary	p s p
+
+atom_style	atomic
+neighbor	0.3 bin
+neigh_modify	delay 5
+
+# create geometry
+
+lattice		hex 0.7
+region		box block 0 40 0 10 -0.25 0.25
+create_box	3 box
+create_atoms	1
+
+mass		1 1.0
+mass		2 1.0
+mass		3 1.0
+
+# LJ potentials
+
+pair_style	lj/cut 1.12246
+pair_coeff	* * 1.0 1.0 1.12246
+
+# define groups
+
+region	        1 block INF INF INF 1.25 INF INF
+group		lower region 1
+region		2 block INF INF 8.75 INF INF INF
+group		upper region 2
+group		boundary union lower upper
+group		flow subtract all boundary
+
+set		lower atom 2
+set		upper atom 3
+
+# initial velocities
+
+temperature     mobile flow full
+velocity	flow create 1.0 482748 temp mobile
+fix		1 all nve
+fix		2 flow temp/rescale 200 1.0 1.0 0.02 1.0
+fix_modify	2 temp mobile
+
+# Poiselle flow
+
+velocity	boundary set 0.0 0.0 0.0
+fix		3 lower setforce 0.0 0.0 0.0
+fix		4 upper setforce 0.0 NULL 0.0
+fix		5 upper aveforce 0.0 -0.5 0.0
+fix		6 flow addforce 1.0 0.0 0.0
+
+# 2 obstacles
+
+region		void1 sphere 10 4 0 3
+delete_atoms	region void1
+region		void2 sphere 20 7 0 3
+delete_atoms	region void2
+
+fix		7 flow indent 100 sphere 10 4 0 4
+fix		8 flow indent 100 sphere 20 7 0 4
+fix		9 all enforce2d
+
+# Run
+
+timestep	0.003
+thermo		1000
+thermo_modify	temp mobile
+
+dump		1 all atom 100 dump.obstacle
+#dump		1 all custom 100 dump.custom tag type x y z vx vy 
+
+run		25000