ICODE-ADC/lammps
4
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

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

Browse Source
This commit is contained in:
sjplimp
2009-11-06 17:05:06 +00:00
parent 409f1e66c4
commit 486ab931c3
8 changed files with 299 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files

24
doc/Section_commands.html
View File

@ -371,18 +371,18 @@ potentials. Click on the style itself for a full description:
<TR ALIGN="center"><TD ><A HREF = "pair_none.html">none</A></TD><TD ><A HREF = "pair_hybrid.html">hybrid</A></TD><TD ><A HREF = "pair_hybrid.html">hybrid/overlay</A></TD><TD ><A HREF = "pair_airebo.html">airebo</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_born.html">born/coul/long</A></TD><TD ><A HREF = "pair_buck.html">buck</A></TD><TD ><A HREF = "pair_buck.html">buck/coul/cut</A></TD><TD ><A HREF = "pair_buck.html">buck/coul/long</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_colloid.html">colloid</A></TD><TD ><A HREF = "pair_coul.html">coul/cut</A></TD><TD ><A HREF = "pair_coul.html">coul/debye</A></TD><TD ><A HREF = "pair_coul.html">coul/long</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_dipole.html">dipole/cut</A></TD><TD ><A HREF = "pair_dpd.html">dpd</A></TD><TD ><A HREF = "pair_eam.html">eam</A></TD><TD ><A HREF = "pair_eam.html">eam/opt</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_eam.html">eam/alloy</A></TD><TD ><A HREF = "pair_eam.html">eam/alloy/opt</A></TD><TD ><A HREF = "pair_eam.html">eam/fs</A></TD><TD ><A HREF = "pair_eam.html">eam/fs/opt</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_gayberne.html">gayberne</A></TD><TD ><A HREF = "pair_gayberne.html">gayberne/gpu</A></TD><TD ><A HREF = "pair_gran.html">gran/hertz/history</A></TD><TD ><A HREF = "pair_gran.html">gran/hooke</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_gran.html">gran/hooke/history</A></TD><TD ><A HREF = "pair_charmm.html">lj/charmm/coul/charmm</A></TD><TD ><A HREF = "pair_charmm.html">lj/charmm/coul/charmm/implicit</A></TD><TD ><A HREF = "pair_charmm.html">lj/charmm/coul/long</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_charmm.html">lj/charmm/coul/long/opt</A></TD><TD ><A HREF = "pair_class2.html">lj/class2</A></TD><TD ><A HREF = "pair_class2.html">lj/class2/coul/cut</A></TD><TD ><A HREF = "pair_class2.html">lj/class2/coul/long</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_lj.html">lj/cut</A></TD><TD ><A HREF = "pair_lj.html">lj/cut/gpu</A></TD><TD ><A HREF = "pair_lj.html">lj/cut/opt</A></TD><TD ><A HREF = "pair_lj.html">lj/cut/coul/cut</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_lj.html">lj/cut/coul/debye</A></TD><TD ><A HREF = "pair_lj.html">lj/cut/coul/long</A></TD><TD ><A HREF = "pair_lj.html">lj/cut/coul/long/tip4p</A></TD><TD ><A HREF = "pair_lj_expand.html">lj/expand</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_gromacs.html">lj/gromacs</A></TD><TD ><A HREF = "pair_gromacs.html">lj/gromacs/coul/gromacs</A></TD><TD ><A HREF = "pair_lj_smooth.html">lj/smooth</A></TD><TD ><A HREF = "pair_lj96_cut.html">lj96/cut</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_lubricate.html">lubricate</A></TD><TD ><A HREF = "pair_meam.html">meam</A></TD><TD ><A HREF = "pair_morse.html">morse</A></TD><TD ><A HREF = "pair_morse.html">morse/opt</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_peri_pmb.html">peri/pmb</A></TD><TD ><A HREF = "pair_reax.html">reax</A></TD><TD ><A HREF = "pair_resquared.html">resquared</A></TD><TD ><A HREF = "pair_soft.html">soft</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_sw.html">sw</A></TD><TD ><A HREF = "pair_table.html">table</A></TD><TD ><A HREF = "pair_tersoff.html">tersoff</A></TD><TD ><A HREF = "pair_tersoff_zbl.html">tersoff/zbl</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_yukawa.html">yukawa</A>
<TR ALIGN="center"><TD ><A HREF = "pair_dipole.html">dipole/cut</A></TD><TD ><A HREF = "pair_dpd.html">dpd</A></TD><TD ><A HREF = "pair_dsmc.html">dsmc</A></TD><TD ><A HREF = "pair_eam.html">eam</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_eam.html">eam/opt</A></TD><TD ><A HREF = "pair_eam.html">eam/alloy</A></TD><TD ><A HREF = "pair_eam.html">eam/alloy/opt</A></TD><TD ><A HREF = "pair_eam.html">eam/fs</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_eam.html">eam/fs/opt</A></TD><TD ><A HREF = "pair_gayberne.html">gayberne</A></TD><TD ><A HREF = "pair_gayberne.html">gayberne/gpu</A></TD><TD ><A HREF = "pair_gran.html">gran/hertz/history</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_gran.html">gran/hooke</A></TD><TD ><A HREF = "pair_gran.html">gran/hooke/history</A></TD><TD ><A HREF = "pair_charmm.html">lj/charmm/coul/charmm</A></TD><TD ><A HREF = "pair_charmm.html">lj/charmm/coul/charmm/implicit</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_charmm.html">lj/charmm/coul/long</A></TD><TD ><A HREF = "pair_charmm.html">lj/charmm/coul/long/opt</A></TD><TD ><A HREF = "pair_class2.html">lj/class2</A></TD><TD ><A HREF = "pair_class2.html">lj/class2/coul/cut</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_class2.html">lj/class2/coul/long</A></TD><TD ><A HREF = "pair_lj.html">lj/cut</A></TD><TD ><A HREF = "pair_lj.html">lj/cut/gpu</A></TD><TD ><A HREF = "pair_lj.html">lj/cut/opt</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_lj.html">lj/cut/coul/cut</A></TD><TD ><A HREF = "pair_lj.html">lj/cut/coul/debye</A></TD><TD ><A HREF = "pair_lj.html">lj/cut/coul/long</A></TD><TD ><A HREF = "pair_lj.html">lj/cut/coul/long/tip4p</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_lj_expand.html">lj/expand</A></TD><TD ><A HREF = "pair_gromacs.html">lj/gromacs</A></TD><TD ><A HREF = "pair_gromacs.html">lj/gromacs/coul/gromacs</A></TD><TD ><A HREF = "pair_lj_smooth.html">lj/smooth</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_lj96_cut.html">lj96/cut</A></TD><TD ><A HREF = "pair_lubricate.html">lubricate</A></TD><TD ><A HREF = "pair_meam.html">meam</A></TD><TD ><A HREF = "pair_morse.html">morse</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_morse.html">morse/opt</A></TD><TD ><A HREF = "pair_peri_pmb.html">peri/pmb</A></TD><TD ><A HREF = "pair_reax.html">reax</A></TD><TD ><A HREF = "pair_resquared.html">resquared</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_soft.html">soft</A></TD><TD ><A HREF = "pair_sw.html">sw</A></TD><TD ><A HREF = "pair_table.html">table</A></TD><TD ><A HREF = "pair_tersoff.html">tersoff</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_tersoff_zbl.html">tersoff/zbl</A></TD><TD ><A HREF = "pair_yukawa.html">yukawa</A>
</TD></TR></TABLE></DIV>
<P>These are pair styles contributed by users, which can be used if

1
doc/Section_commands.txt
View File

@ -515,6 +515,7 @@ potentials. Click on the style itself for a full description:
"coul/long"_pair_coul.html,
"dipole/cut"_pair_dipole.html,
"dpd"_pair_dpd.html,
"dsmc"_pair_dsmc.html,
"eam"_pair_eam.html,
"eam/opt"_pair_eam.html,
"eam/alloy"_pair_eam.html,

1
doc/pair_coeff.html
View File

@ -98,6 +98,7 @@ the pair_style command, and coefficients specified by the associated
<LI><A HREF = "pair_coul.html">pair_style coul/long</A> - long-range Coulombic potential
<LI><A HREF = "pair_dipole.html">pair_style dipole/cut</A> - point dipoles with cutoff
<LI><A HREF = "pair_dpd.html">pair_style dpd</A> - dissipative particle dynamics (DPD)
<LI><A HREF = "pair_dsmc.html">pair_style dsmc</A> - Direct Simulation Monte Carlo (DSMC)
<LI><A HREF = "pair_eam.html">pair_style eam</A> - embedded atom method (EAM)
<LI><A HREF = "pair_eam.html">pair_style eam/opt</A> - optimized version of EAM
<LI><A HREF = "pair_eam.html">pair_style eam/alloy</A> - alloy EAM

1
doc/pair_coeff.txt
View File

@ -95,6 +95,7 @@ the pair_style command, and coefficients specified by the associated
"pair_style coul/long"_pair_coul.html - long-range Coulombic potential
"pair_style dipole/cut"_pair_dipole.html - point dipoles with cutoff
"pair_style dpd"_pair_dpd.html - dissipative particle dynamics (DPD)
"pair_style dsmc"_pair_dsmc.html - Direct Simulation Monte Carlo (DSMC)
"pair_style eam"_pair_eam.html - embedded atom method (EAM)
"pair_style eam/opt"_pair_eam.html - optimized version of EAM
"pair_style eam/alloy"_pair_eam.html - alloy EAM

144
doc/pair_dsmc.html Normal file
View File

@ -0,0 +1,144 @@
<HTML>
<CENTER><A HREF = "http://lammps.sandia.gov">LAMMPS WWW Site</A> - <A HREF = "Manual.html">LAMMPS Documentation</A> - <A HREF = "Section_commands.html#comm">LAMMPS Commands</A>
</CENTER>
<HR>
<H3>pair_style dsmc command
</H3>
<P><B>Syntax:</B>
</P>
<PRE>pair_style dsmc max_cell_size seed weighting Tref Nrecompute Nsample
</PRE>
<UL><LI>max_cell_size = global maximum cell size for DSMC interactions (distance units)
<LI>seed = random # seed (positive integer)
<LI>weighting = macroparticle weighting
<LI>Tref = reference temperature (temperature units)
<LI>Nrecompute = recompute v*sigma_max every this many timesteps (timesteps)
<LI>Nsample = sample this many times in recomputing v*sigma_max
</UL>
<P><B>Examples:</B>
</P>
<PRE>pair_style dsmc 2.5 34387 10 1.0 100 20
pair_coeff * * 1.0
pair_coeff 1 1 1.0
</PRE>
<P><B>Description:</B>
</P>
<P>Style <I>dsmc</I> computes collisions between pairs of particles for a
direct simulation Monte Carlo (DSMC) model following the exposition in
<A HREF = "#Bird">(Bird)</A>. Each collision resets the velocities of the two
particles involved. The number of pairwise collisions for each pair
or particle types and the length scale within which they occur are
determined by the parameters of the pair_style and pair_coeff
commands.
</P>
<P>Stochastic collisions are performed using the variable hard sphere
(VHS) approach, with the user-defined <I>max_cell_size</I> value used as
the maximum DSMC cell size, and reference cross-sections for
collisions given using the pair_coeff command.
</P>
<P>There is no pairwise energy or virial contributions associated with
this pair style.
</P>
<P>The following coefficient must be defined for each pair of atoms types
via the <A HREF = "pair_coeff.html">pair_coeff</A> command as in the examples above,
or in the data file or restart files read by the
<A HREF = "read_data.html">read_data</A> or <A HREF = "read_restart.html">read_restart</A>
commands:
</P>
<UL><LI>sigma (area units, i.e. distance-squared)
</UL>
<P>The global DSMC <I>max_cell_size</I> determines the maximum cell length
used in the DSMC calculation. A structured mesh is overlayed on the
simulation box such that an integer number of cells are created in
each direction for each processor's sub-domain. Cell lengths are
adjusted up to the user-specified maximum cell size.
</P>
<HR>
<P>To perform a DSMC simulation with LAMMPS, several additional options
should be set in your input script, though LAMMPS does not check for
these settings.
</P>
<P>Since this pair style does not compute particle forces, you should use
the "fix nve/noforce" time integration fix for the DSMC particles,
e.g.
</P>
<PRE>fix 1 all nve/noforce
</PRE>
<P>This pair style assumes that all particles will communicated to
neighboring processors every timestep as they move. This makes it
possible to perform all collisions between pairs of particles that are
on the same processor. To ensure this occurs, you should use
these commands:
</P>
<PRE>neighbor 0.0 bin
neigh_modify every 1 delay 0 check no
communicate single cutoff 0.0
</PRE>
<P>These commands insure that LAMMPS communicates particles to
neighboring processors every timestep and that no ghost atoms are
created. The output statistics for a simulation run should indicate
there are no ghost particles or neighbors.
</P>
<HR>
<P><B>Mixing, shift, table, tail correction, restart, rRESPA info</B>:
</P>
<P>This pair style does not support mixing. Thus, coefficients for all
I,J pairs must be specified explicitly.
</P>
<P>This pair style does not support the <A HREF = "pair_modify.html">pair_modify</A>
shift option for the energy of the pair interaction.
</P>
<P>The <A HREF = "pair_modify.html">pair_modify</A> table option is not relevant
for this pair style.
</P>
<P>This pair style does not support the <A HREF = "pair_modify.html">pair_modify</A>
tail option for adding long-range tail corrections to energy and
pressure.
</P>
<P>This pair style writes its information to <A HREF = "restart.html">binary restart
files</A>, so pair_style and pair_coeff commands do not need
to be specified in an input script that reads a restart file. Note
that the user-specified random number seed is stored in the restart
file, so when a simulation is restarted, each processor will
re-initialize its random number generator the same way it did
initially. This means the random forces will be random, but will not
be the same as they would have been if the original simulation had
continued past the restart time.
</P>
<P>This pair style can only be used via the <I>pair</I> keyword of the
<A HREF = "run_style.html">run_style respa</A> command. It does not support the
<I>inner</I>, <I>middle</I>, <I>outer</I> keywords.
</P>
<HR>
<P><B>Restrictions:</B>
</P>
<P>This style is part of the "dsmc" package. It is only enabled if
LAMMPS was built with that package. See the <A HREF = "Section_start.html#2_3">Making
LAMMPS</A> section for more info.
</P>
<P><B>Related commands:</B>
</P>
<P><A HREF = "pair_coeff.html">pair_coeff</A>, <A HREF = "fix_nve_noforce.html">fix nve/noforce</A>,
<A HREF = "neigh_modify.html">neigh_modify</A>, <A HREF = "neighbor.html">neighbor</A>,
<A HREF = "communicate.html">communicate</A>
</P>
<P><B>Default:</B> none
</P>
<HR>
<A NAME = "Bird"></A>
<P><B>(Bird)</B> G. A. Bird, "Molecular Gas Dynamics and the Direct Simulation
of Gas Flows" (1994).
</P>
</HTML>

138
doc/pair_dsmc.txt Normal file
View File

@ -0,0 +1,138 @@
"LAMMPS WWW Site"_lws - "LAMMPS Documentation"_ld - "LAMMPS Commands"_lc :c
:link(lws,http://lammps.sandia.gov)
:link(ld,Manual.html)
:link(lc,Section_commands.html#comm)
:line
pair_style dsmc command :h3
[Syntax:]
pair_style dsmc max_cell_size seed weighting Tref Nrecompute Nsample :pre
max_cell_size = global maximum cell size for DSMC interactions (distance units)
seed = random # seed (positive integer)
weighting = macroparticle weighting
Tref = reference temperature (temperature units)
Nrecompute = recompute v*sigma_max every this many timesteps (timesteps)
Nsample = sample this many times in recomputing v*sigma_max :ul
[Examples:]
pair_style dsmc 2.5 34387 10 1.0 100 20
pair_coeff * * 1.0
pair_coeff 1 1 1.0 :pre
[Description:]
Style {dsmc} computes collisions between pairs of particles for a
direct simulation Monte Carlo (DSMC) model following the exposition in
"(Bird)"_#Bird. Each collision resets the velocities of the two
particles involved. The number of pairwise collisions for each pair
or particle types and the length scale within which they occur are
determined by the parameters of the pair_style and pair_coeff
commands.
Stochastic collisions are performed using the variable hard sphere
(VHS) approach, with the user-defined {max_cell_size} value used as
the maximum DSMC cell size, and reference cross-sections for
collisions given using the pair_coeff command.
There is no pairwise energy or virial contributions associated with
this pair style.
The following coefficient must be defined for each pair of atoms types
via the "pair_coeff"_pair_coeff.html command as in the examples above,
or in the data file or restart files read by the
"read_data"_read_data.html or "read_restart"_read_restart.html
commands:
sigma (area units, i.e. distance-squared) :ul
The global DSMC {max_cell_size} determines the maximum cell length
used in the DSMC calculation. A structured mesh is overlayed on the
simulation box such that an integer number of cells are created in
each direction for each processor's sub-domain. Cell lengths are
adjusted up to the user-specified maximum cell size.
:line
To perform a DSMC simulation with LAMMPS, several additional options
should be set in your input script, though LAMMPS does not check for
these settings.
Since this pair style does not compute particle forces, you should use
the "fix nve/noforce" time integration fix for the DSMC particles,
e.g.
fix 1 all nve/noforce :pre
This pair style assumes that all particles will communicated to
neighboring processors every timestep as they move. This makes it
possible to perform all collisions between pairs of particles that are
on the same processor. To ensure this occurs, you should use
these commands:
neighbor 0.0 bin
neigh_modify every 1 delay 0 check no
communicate single cutoff 0.0 :pre
These commands insure that LAMMPS communicates particles to
neighboring processors every timestep and that no ghost atoms are
created. The output statistics for a simulation run should indicate
there are no ghost particles or neighbors.
:line
[Mixing, shift, table, tail correction, restart, rRESPA info]:
This pair style does not support mixing. Thus, coefficients for all
I,J pairs must be specified explicitly.
This pair style does not support the "pair_modify"_pair_modify.html
shift option for the energy of the pair interaction.
The "pair_modify"_pair_modify.html table option is not relevant
for this pair style.
This pair style does not support the "pair_modify"_pair_modify.html
tail option for adding long-range tail corrections to energy and
pressure.
This pair style writes its information to "binary restart
files"_restart.html, so pair_style and pair_coeff commands do not need
to be specified in an input script that reads a restart file. Note
that the user-specified random number seed is stored in the restart
file, so when a simulation is restarted, each processor will
re-initialize its random number generator the same way it did
initially. This means the random forces will be random, but will not
be the same as they would have been if the original simulation had
continued past the restart time.
This pair style can only be used via the {pair} keyword of the
"run_style respa"_run_style.html command. It does not support the
{inner}, {middle}, {outer} keywords.
:line
[Restrictions:]
This style is part of the "dsmc" package. It is only enabled if
LAMMPS was built with that package. See the "Making
LAMMPS"_Section_start.html#2_3 section for more info.
[Related commands:]
"pair_coeff"_pair_coeff.html, "fix nve/noforce"_fix_nve_noforce.html,
"neigh_modify"_neigh_modify.html, "neighbor"_neighbor.html,
"communicate"_communicate.html
[Default:] none
:line
:link(Bird)
[(Bird)] G. A. Bird, "Molecular Gas Dynamics and the Direct Simulation
of Gas Flows" (1994).

1
doc/pair_style.html
View File

@ -100,6 +100,7 @@ the pair_style command, and coefficients specified by the associated
<LI><A HREF = "pair_coul.html">pair_style coul/long</A> - long-range Coulombic potential
<LI><A HREF = "pair_dipole.html">pair_style dipole/cut</A> - point dipoles with cutoff
<LI><A HREF = "pair_dpd.html">pair_style dpd</A> - dissipative particle dynamics (DPD)
<LI><A HREF = "pair_dsmc.html">pair_style dsmc</A> - Direct Simulation Monte Carlo (DSMC)
<LI><A HREF = "pair_eam.html">pair_style eam</A> - embedded atom method (EAM)
<LI><A HREF = "pair_eam.html">pair_style eam/opt</A> - optimized version of EAM
<LI><A HREF = "pair_eam.html">pair_style eam/alloy</A> - alloy EAM

1
doc/pair_style.txt
View File

@ -97,6 +97,7 @@ the pair_style command, and coefficients specified by the associated
"pair_style coul/long"_pair_coul.html - long-range Coulombic potential
"pair_style dipole/cut"_pair_dipole.html - point dipoles with cutoff
"pair_style dpd"_pair_dpd.html - dissipative particle dynamics (DPD)
"pair_style dsmc"_pair_dsmc.html - Direct Simulation Monte Carlo (DSMC)
"pair_style eam"_pair_eam.html - embedded atom method (EAM)
"pair_style eam/opt"_pair_eam.html - optimized version of EAM
"pair_style eam/alloy"_pair_eam.html - alloy EAM