← Back to branch summary

~ubuntu-branches/debian/sid/lammps/sid

« back to all changes in this revision

Viewing changes to examples/USER/lb/dragforce/in.setgamma_drag

  • Committer: Package Import Robot
  • Author(s): Anton Gladky
  • Date: 2013-11-20 22:41:36 UTC
  • mfrom: (1.2.2)
  • Revision ID: package-import@ubuntu.com-20131120224136-tzx7leh606fqnckm
Tags: 0~20131119.git7162cf0-1
* [e65b919] Imported Upstream version 0~20131119.git7162cf0
* [f7bddd4] Fix some problems, introduced by upstream recently.
* [3616dfc] Use wrap-and-sort script.
* [7e92030] Ignore quilt dir

Show diffs side-by-side

added added

removed removed

Lines of Context:
examples/USER/lb/dragforce/in.setgamma_drag
 
1
#===========================================================================#
 
2
# Drag force on a single sphere.                                            #
 
3
#                                                                           #
 
4
# Here, gamma (used in the calculation of the particle-fluid interaction    #
 
5
#   force) is set by the user (gamma = 13.0 for this simulation.)  This     #
 
6
#   type of simulation is used to calibrate the value for gamma which will  #
 
7
#   give the desired Stokes drag force.                                     #
 
8
#                                                                           #
 
9
# Sample output from this run can be found in the file:                     # 
 
10
#   'setgamma13d0_drag.out'                                                 #
 
11
#===========================================================================# 
 
12
 
 
13
units          micro
 
14
dimension      3
 
15
boundary       p p f
 
16
atom_style     atomic
 
17
 
 
18
#----------------------------------------------------------------------------
 
19
# Need a neighbor bin size smaller than the lattice-Boltzmann grid spacing
 
20
#   to ensure that the particles belonging to a given processor remain inside
 
21
#   that processors lattice-Boltzmann grid.
 
22
# The arguments for neigh_modify have been set to "delay 0 every 1", again
 
23
#   to ensure that the particles belonging to a given processor remain inside
 
24
#   that processors lattice-Boltzmann grid.  However, these values can likely
 
25
#   be somewhat increased without issue.  If a problem does arise (a particle
 
26
#   is outside of its processors LB grid) an error message is printed and 
 
27
#   the simulation is terminated.  
 
28
#----------------------------------------------------------------------------
 
29
neighbor       1.0 bin
 
30
neigh_modify   delay 0 every 1
 
31
 
 
32
read_data      data.one_radius16d2
 
33
 
 
34
#----------------------------------------------------------------------------
 
35
# None of the particles comprising the spherical colloidal object should
 
36
#   interact with one another.
 
37
#----------------------------------------------------------------------------
 
38
pair_style      lj/cut 2.45
 
39
pair_coeff      * * 0.0 0.0 2.45
 
40
neigh_modify    exclude type 1 1 
 
41
 
 
42
mass * 1.0
 
43
timestep 4.0
 
44
 
 
45
velocity all set 0.0 0.0001 0.0 units box
 
46
 
 
47
#---------------------------------------------------------------------------
 
48
# Create a lattice-Boltzmann fluid covering the simulation domain.
 
49
# All of the particles in the simulation apply a force to the fluid.
 
50
# Use the LB integration scheme of Ollila et. al. (for stability reasons,
 
51
#   this integration scheme should be used when a large user set value for 
 
52
#   gamma is specified), a fluid density = 1.0, fluid viscosity = 1.0, value 
 
53
#   for gamma=13.0, lattice spacing dx=4.0, and mass unit, dm=10.0.
 
54
# Print the force and torque acting on the particle to the screen at each 
 
55
#   timestep.
 
56
#----------------------------------------------------------------------------
 
57
fix     1 all lb/fluid 1 2 1.0 1.0 setGamma 13.0 dx 4.0 dm 10.0 calcforce 10 all 
 
58
 
 
59
#---------------------------------------------------------------------------
 
60
# For this simulation the colloidal particle moves at a constant velocity
 
61
#   through the fluid.  As such, we do not wish to apply the force from
 
62
#   the fluid back onto the object.  Therefore, we do not use any of the
 
63
#   viscous_lb, rigid_pc_sphere, or pc fixes, and simply integrate the 
 
64
#   particle motion using one of the built-in LAMMPS integrators.
 
65
#---------------------------------------------------------------------------
 
66
fix     2 all nve
 
67
 
 
68
run     100000