"LIGGGHTS(R)-PUBLIC WWW Site"_liws - "LIGGGHTS(R)-PUBLIC Documentation"_ld - "LIGGGHTS(R)-PUBLIC Commands"_lc :c

:link(liws,http://www.cfdem.com)
:link(ld,Manual.html)
:link(lc,Section_commands.html#comm)

:line

boundary command :h3

[Syntax:]

boundary x y z :pre

x,y,z = {p} or {s} or {f} or {m}, one or two letters :ulb,l
  {p} is periodic
  {f} is non-periodic and fixed
  {s} is non-periodic and shrink-wrapped
  {m} is non-periodic and shrink-wrapped with a minimum value :pre
:ule

[Examples:]

boundary p p f
boundary p fs p
boundary s f fm :pre

[Description:]

Set the style of boundaries for the global simulation box in each
dimension.  A single letter assigns the same style to both the lower
and upper face of the box.  Two letters assigns the first style to the
lower face and the second style to the upper face.  The initial size
of the simulation box is set by the "read_data"_read_data.html,
"read_restart"_read_restart.html, or "create_box"_create_box.html
commands.

The style {p} means the box is periodic, so that particles interact
across the boundary, and they can exit one end of the box and re-enter
the other end.  A periodic dimension can change in size due to
constant pressure boundary conditions or box deformation (see the "fix
npt"_fix_nh.html and "fix deform"_fix_deform.html commands).  The {p}
style must be applied to both faces of a dimension.

The styles {f}, {s}, and {m} mean the box is non-periodic, so that
particles do not interact across the boundary and do not move from one
side of the box to the other.  For style {f}, the position of the face
is fixed.  If an atom moves outside the face it may be lost.  For
style {s}, the position of the face is set so as to encompass the
atoms in that dimension (shrink-wrapping), no matter how far they
move.  For style {m}, shrink-wrapping occurs, but is bounded by the
value specified in the data or restart file or set by the
"create_box"_create_box.html command.  For example, if the upper z
face has a value of 50.0 in the data file, the face will always be
positioned at 50.0 or above, even if the maximum z-extent of all the
atoms becomes less than 50.0.

For triclinic (non-orthogonal) simulation boxes, if the 2nd dimension
of a tilt factor (e.g. y for xy) is periodic, then the periodicity is
enforced with the tilt factor offset.  If the 1st dimension is
shrink-wrapped, then the shrink wrapping is applied to the tilted box
face, to encompass the atoms.  E.g. for a positive xy tilt, the xlo
and xhi faces of the box are planes tilting in the +y direction as y
increases.  These tilted planes are shrink-wrapped around the atoms to
determine the x extent of the box.

See "Section_howto 12"_Section_howto.html#howto_7 of the doc pages
for a geometric description of triclinic boxes, as defined by LIGGGHTS(R)-PUBLIC,
and how to transform these parameters to and from other commonly used
triclinic representations.

IMPORTANT NOTE: If mesh walls (e.g. "fix mesh/surface"_fix_mesh_surface.html)
are used, not only atom positions, but also the mesh nodes are used
for setting the boundaries.

[Restrictions:]

This command cannot be used after the simulation box is defined by a
"read_data"_read_data.html or "create_box"_create_box.html command or
"read_restart"_read_restart.html command.  See the
"change_box"_change_box.html command for how to change the simulation
box boundaries after it has been defined.

For 2d simulations, the z dimension must be periodic.

[Related commands:]

See the "thermo_modify"_thermo_modify.html command for a discussion
of lost atoms.

[Default:]

boundary p p p :pre
