From fbfd916c6998c131132ddf95069361900898b336 Mon Sep 17 00:00:00 2001
From: sjplimp langevin/eff lb/fluid lb/momentum lb/pc lb/rigid/pc/sphere lb/viscous meso meso/stationary nph/eff npt/eff nve/eff nvt/eff
-nvt/sllod/eff phonon pimd qbmsst qeq/reax qmmm
-qtb reax/c/bonds reax/c/species saed/vtk smd temp/rescale/eff ti/rs ti/spring ttm/mod
+
+qtb reax/c/bonds reax/c/species saed/vtk smd smd/adjust/dt
+smd/integrate/tlsph smd/integrate/ulsph smd/move/triangulated/surface smd/setvel smd/tlsph/reference/configuration smd/wall/surface temp/rescale/eff ti/rs ti/spring ttm/mod
@@ -474,8 +475,11 @@ package.
This package implements smoothed Mach dynamics (SMD) in +LAMMPS. Currently, the package has the following features: +
+??? +
+See the file doc/USER/smd/SPH_user_guide.pdf to get started. +
+There are example scripts for using this package in examples/USER/smd. +
+The person who created this package is Georg Ganzenmuller at the +Fraunhofer-Institute for High-Speed Dynamics, Ernst Mach Institute in +Germany (georg.ganzenmueller at emi.fhg.de). Contact him directly if +you have questions. +
This package implements smoothed particle hydrodynamics (SPH) in diff --git a/doc/doc2/atom_style.html b/doc/doc2/atom_style.html index 6c690b135f..a8a9cc45d8 100644 --- a/doc/doc2/atom_style.html +++ b/doc/doc2/atom_style.html @@ -15,7 +15,7 @@
atom_style style args-
args = none for any style except body and hybrid body args = bstyle bstyle-args @@ -80,6 +80,7 @@ quantities.meso rho, e, cv SPH particles molecular bonds, angles, dihedrals, impropers uncharged molecules + peri mass, volume mesocopic Peridynamic models smd volume, kernel diameter, contact radius, mass solid and fluid SPH particles sphere diameter, mass, angular velocity granular models template template index, template atom small molecules with fixed topology @@ -128,6 +129,13 @@ per-particle mass and volume. particles which store a density (rho), energy (e), and heat capacity (cv). + tri corner points, angular momentum rigid bodies The smd style is for a general formulation of Smooth Particle +Hydrodynamics. Both fluids and solids can be modeled. Particles +store the mass and volume of an integration point, a kernel diameter +used for calculating the field variables (e.g. stress and deformation) +and a contact radius for calculating repulsive forces which prevent +individual physical bodies from penetretating each other. +
The wavepacket style is similar to electron, but the electrons may consist of several Gaussian wave packets, summed up with coefficients cs= (cs_re,cs_im). Each of the wave packets is treated as a separate diff --git a/doc/doc2/compute_smd_contact_radius.html b/doc/doc2/compute_smd_contact_radius.html new file mode 100644 index 0000000000..e2c2881821 --- /dev/null +++ b/doc/doc2/compute_smd_contact_radius.html @@ -0,0 +1,58 @@ + +
LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Commands + + + + + + + +
+ +compute smd/contact_radius command +
+Syntax: +
+compute ID group-ID smd/contact_radius ++
Examples: +
+compute 1 all smd/contact_radius ++
Description: +
+Define a computation which outputs the contact radius, i.e., the radius used to prevent particles from penetrating each other. +The contact radius is used only to prevent particles belonging to different physical bodies from penetrating each other. It is used by the contact pair styles, +e.g., smd/hertz and smd/tri_surface. +
+See this PDF guide to using Smooth Mach Dynamics in +LAMMPS. +
+The value of the contact radius will be 0.0 for particles not in the +specified compute group. +
+Output info: +
+This compute calculates a per-particle vector, which can be accessed by +any command that uses per-particle values from a compute as input. See +Section_howto 15 for an overview of +LAMMPS output options. +
+The per-particle vector values will be in distance units. +
+Restrictions: +
+This compute is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making +LAMMPS section for more info. +
+Related commands: +
+dump custom smd/hertz smd/tri_surface +
+Default: none +
+ diff --git a/doc/doc2/compute_smd_damage.html b/doc/doc2/compute_smd_damage.html new file mode 100644 index 0000000000..f4b10c50b2 --- /dev/null +++ b/doc/doc2/compute_smd_damage.html @@ -0,0 +1,53 @@ + +Syntax: +
+compute ID group-ID smd/damage ++
Examples: +
+compute 1 all smd/damage ++
Description: +
+Define a computation that calculates the damage status of SPH particles +according to the damage model which is defined via the SMD SPH pair styles, e.g., the maximum plastic strain failure criterion. +
+See this PDF guide to use Smooth Mach Dynamics in LAMMPS. +
+Output Info: +
+This compute calculates a per-particle vector, which can be accessed by +any command that uses per-particle values from a compute as input. See +How-to discussions, section 6.15 +for an overview of LAMMPS output options. +
+The per-particle values are dimensionless an in the range of zero to one. +
+Restrictions: +
+This compute is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. +
+Related commands: +
+smd/plastic_strain, smd/tlsph_stress +
+Default: none +
+ diff --git a/doc/doc2/compute_smd_hourglass_error.html b/doc/doc2/compute_smd_hourglass_error.html new file mode 100644 index 0000000000..f77e3d40b7 --- /dev/null +++ b/doc/doc2/compute_smd_hourglass_error.html @@ -0,0 +1,61 @@ + +Syntax: +
+compute ID group-ID smd/hourglass_error ++
Examples: +
+compute 1 all smd/hourglass_error ++
Description: +
+Define a computation which outputs the error of the approximated +relative separation with respect to the actual relative separation of +the particles i and j. Ideally, if the deformation gradient is exact, +and there exists a unique mapping between all particles' positions +within the neighborhood of the central node and the deformation gradient, +the approximated relative separation will coincide with the actual relative +separation of the particles i and j in the deformed configuration. +This compute is only really useful for debugging the hourglass control mechanim which is part of the Total-Lagrangian SPH pair style. +
+See this PDF guide to use Smooth Mach Dynamics in LAMMPS. +
+Output Info: +
+This compute calculates a per-particle vector, which can be accessed by +any command that uses per-particle values from a compute as input. See +How-to discussions, section 6.15 +for an overview of LAMMPS output options. +
+The per-particle vector values will are dimensionless. See units. +
+Restrictions: +
+This compute is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. +
+This quantity will be computed only for particles which interact with tlsph pair style. +
+Related Commands: +
+ +Default: +
+ diff --git a/doc/doc2/compute_smd_internal_energy.html b/doc/doc2/compute_smd_internal_energy.html new file mode 100644 index 0000000000..f7852089e5 --- /dev/null +++ b/doc/doc2/compute_smd_internal_energy.html @@ -0,0 +1,51 @@ + +Syntax: +
+compute ID group-ID smd/internal_energy ++
Examples: +
+compute 1 all smd/internal_energy ++
Description: +
+Define a computation which outputs the per-particle enthalpy, i.e., the sum of potential energy and heat. +
+See this PDF guide to use Smooth Mach Dynamics in LAMMPS. +
+Output Info: +
+This compute calculates a per-particle vector, which can be accessed by +any command that uses per-particle values from a compute as input. See +How-to discussions, section 6.15 +for an overview of LAMMPS output options. +
+The per-particle vector values will be given in units of energy. +
+Restrictions: +
+This compute is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. This compute can only be used for particles which interact via the +updated Lagrangian or total Lagrangian SPH pair styles. +
+Related Commands: +
+Default: +
+ diff --git a/doc/doc2/compute_smd_plastic_strain.html b/doc/doc2/compute_smd_plastic_strain.html new file mode 100644 index 0000000000..0a776ee033 --- /dev/null +++ b/doc/doc2/compute_smd_plastic_strain.html @@ -0,0 +1,54 @@ + +Syntax: +
+compute ID group-ID smd/plastic_strain ++
Examples: +
+compute 1 all smd/plastic_strain ++
Description: +
+Define a computation that outputs the equivalent plastic strain per particle. +This command is only meaningful if a material model with plasticity is defined. +
+See this PDF guide to use Smooth Mach Dynamics in LAMMPS. +
+Output Info: +
+This compute calculates a per-particle vector, which can be accessed by +any command that uses per-particle values from a compute as input. See +How-to discussions, section 6.15 +for an overview of LAMMPS output options. +
+The per-particle values will be given dimensionless. See units. +
+Restrictions: +
+This compute is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. This compute can only be used for particles which interact via the +updated Lagrangian or total Lagrangian SPH pair styles. +
+Related commands: +
+smd/plastic_strain_rate, smd/tlsph_strain_rate, +smd/tlsph_strain +
+Default: none
+ diff --git a/doc/doc2/compute_smd_plastic_strain_rate.html b/doc/doc2/compute_smd_plastic_strain_rate.html new file mode 100644 index 0000000000..243c5d7d2c --- /dev/null +++ b/doc/doc2/compute_smd_plastic_strain_rate.html @@ -0,0 +1,54 @@ + +Syntax: +
+compute ID group-ID smd/plastic_strain_rate ++
Examples: +
+compute 1 all smd/plastic_strain_rate ++
Description: +
+Define a computation that outputs the time rate of the equivalent plastic strain. +This command is only meaningful if a material model with plasticity is defined. +
+See this PDF guide to use Smooth Mach Dynamics in LAMMPS. +
+Output Info: +
+This compute calculates a per-particle vector, which can be accessed by +any command that uses per-particle values from a compute as input. See +How-to discussions, section 6.15 +for an overview of LAMMPS output options. +
+The per-particle values will be given in units of one over time. +
+Restrictions: +
+This compute is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. This compute can only be used for particles which interact via the +updated Lagrangian or total Lagrangian SPH pair styles. +
+Related commands: +
+smd/plastic_strain, smd/tlsph_strain_rate, +smd/tlsph_strain +
+Default: none
+ diff --git a/doc/doc2/compute_smd_rho.html b/doc/doc2/compute_smd_rho.html new file mode 100644 index 0000000000..281153a69a --- /dev/null +++ b/doc/doc2/compute_smd_rho.html @@ -0,0 +1,54 @@ + +Syntax: +
+compute ID group-ID smd/rho ++
Examples: +
+compute 1 all smd/rho ++
Description: +
+Define a computation that calculates the per-particle mass density. +The mass density is the mass of a particle which is constant during +the course of a simulation, divided by its volume, which can change due to mechanical deformation. +
+See this PDF guide to use Smooth Mach Dynamics in LAMMPS. +
+Output info: +
+This compute calculates a per-particle vector, which can be accessed by +any command that uses per-particle values from a compute as input. See +How-to discussions, section 6.15 +for an overview of LAMMPS output options. +
+The per-particle values will be in units of mass over volume. +
+Restrictions: +
+This compute is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. +
+Related commands: +
+smd/vol +
+Default: none +
+ diff --git a/doc/doc2/compute_smd_tlsph_defgrad.html b/doc/doc2/compute_smd_tlsph_defgrad.html new file mode 100644 index 0000000000..56c52169db --- /dev/null +++ b/doc/doc2/compute_smd_tlsph_defgrad.html @@ -0,0 +1,57 @@ + +Syntax: +
+compute ID group-ID smd/tlsph_defgrad ++
Examples: +
+compute 1 all smd/tlsph_defgrad ++
Description: +
+Define a computation that calculates the deformation gradient. +It is only meaningful for particles which interact according to the Total-Lagrangian SPH pair style. +
+See this PDF guide to use Smooth Mach Dynamics in +LAMMPS. +
+Output info: +
+This compute outputss a per-particle vector of vectors (tensors), which can be +accessed by any command that uses per-particle values from a compute as input. See +How-to discussions, section 6.15 +for an overview of LAMMPS output options. +
+The per-particle vector values will be given dimensionless. See units. +The per-particle vector has 10 entries. The first nine entries correspond to the xx, xy, xz, yx, yy, yz, zx, zy, zz components +of the asymmetric deformation gradient tensor. The tenth entry is the determinant of the deformation gradient. +
+Restrictions: +
+This compute is part of the USER-SMD package. It is only enabled if LAMMPS was +built with that package. See the Making LAMMPS +section for more info. TThis compute can only be used for particles which interact via the +total Lagrangian SPH pair style. +
+Related commands: +
+ +Default: none +
+ diff --git a/doc/doc2/compute_smd_tlsph_dt.html b/doc/doc2/compute_smd_tlsph_dt.html new file mode 100644 index 0000000000..bec2e7ca8a --- /dev/null +++ b/doc/doc2/compute_smd_tlsph_dt.html @@ -0,0 +1,59 @@ + +Syntax: +
+compute ID group-ID smd/tlsph_dt ++
Examples: +
+compute 1 all smd/tlsph_dt ++
Description: +
+Define a computation that outputs the CFL-stable time increment per particle. +This time increment is essentially given by the speed of sound, divided by the SPH smoothing length. +Because both the speed of sound and the smoothing length typically change during the course of a simulation, +the stable time increment needs to be recomputed every time step. +This calculation is performed automatically in the relevant SPH pair styles and this compute only serves to make the +stable time increment accessible for output purposes. +
+See this PDF guide to using Smooth Mach Dynamics in LAMMPS. +
+Output info: +
+This compute calculates a per-particle vector, which can be accessed by +any command that uses per-particle values from a compute as input. See +How-to discussions, section 6.15 +for an overview of LAMMPS output options. +
+The per-particle values will be given in units of time. +
+Restrictions: +
+This compute is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. +
+This compute can only be used for particles interacting with the Total-Lagrangian SPH pair style. +
+Related commands: +
+ +Default: none +
+ diff --git a/doc/doc2/compute_smd_tlsph_num_neighs.html b/doc/doc2/compute_smd_tlsph_num_neighs.html new file mode 100644 index 0000000000..15a1ccdcae --- /dev/null +++ b/doc/doc2/compute_smd_tlsph_num_neighs.html @@ -0,0 +1,55 @@ + +Syntax: +
+compute ID group-ID smd/tlsph_num_neighs ++
Examples: +
+compute 1 all smd/tlsph_num_neighs ++
Description: +
+Define a computation that calculates the number of particles +inside of the smoothing kernel radius for particles interacting via the Total-Lagrangian SPH pair style. +
+See this PDF guide to using Smooth Mach Dynamics in LAMMPS. +
+Output info: +
+This compute calculates a per-particle vector, which can be accessed by +any command that uses per-particle values from a compute as input. See +How-to discussions, section 6.15 +for an overview of LAMMPS output options. +
+The per-particle values are dimensionless. See units. +
+Restrictions: +
+This compute is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. +
+This quantity will be computed only for particles which interact with the Total-Lagrangian pair style. +
+Related commands: +
+ +Default: none +
+ diff --git a/doc/doc2/compute_smd_tlsph_shape.html b/doc/doc2/compute_smd_tlsph_shape.html new file mode 100644 index 0000000000..ff4557350e --- /dev/null +++ b/doc/doc2/compute_smd_tlsph_shape.html @@ -0,0 +1,57 @@ + +Syntax: +
+compute ID group-ID smd/tlsph_shape ++
Examples: +
+compute 1 all smd/tlsph_shape ++
Description: +
+Define a computation that outputs the current shape of the volume associated with a particle as a rotated ellipsoid. +It is only meaningful for particles which interact according to the Total-Lagrangian SPH pair style. +
+See this PDF guide to use Smooth Mach Dynamics in LAMMPS. +
+Output info: +
+This compute calculates a per-particle vector of vectors, which can be +accessed by any command that uses per-particle values from a compute as input. See +How-to discussions, section 6.15 +for an overview of LAMMPS output options. +
+The per-particle vector has 7 entries. The first three entries correspond to the lengths of the ellipsoid's axes and have units of length. +These axis valus are computed as the contact radius times the xx, yy, or zz components of the Green-Lagrange strain tensor associated with the particle. +The next 4 values are quaternions (order: q, x, y, z) which describe the spatial rotation of the particle relative to its initial state. +
+Restrictions: +
+This compute is part of the USER-SMD package. It is only enabled if LAMMPS was +built with that package. See the Making LAMMPS +section for more info. +
+This quantity will be computed only for particles which interact with the Total-Lagrangian SPH pair style. +
+Related commands: +
+ +Default: none +
+ diff --git a/doc/doc2/compute_smd_tlsph_strain.html b/doc/doc2/compute_smd_tlsph_strain.html new file mode 100644 index 0000000000..d1a44dadc1 --- /dev/null +++ b/doc/doc2/compute_smd_tlsph_strain.html @@ -0,0 +1,58 @@ + +Syntax: +
+compute ID group-ID smd/tlsph_strain ++
Examples: +
+compute 1 all smd/tlsph_strain ++
Description: +
+Define a computation that calculates the Green-Lagrange strain tensor for particles interacting via the Total-Lagrangian SPH pair style. +
+See this PDF guide to using Smooth Mach Dynamics in LAMMPS. +
+Output info: +
+This compute calculates a per-particle vector of vectors (tensors), +which can be accessed by any command that uses per-particle values +from a compute as input. See +How-to discussions, section 6.15 + for an overview of LAMMPS output options. +
+The per-particle tensor values will be given dimensionless. See units. +
+The per-particle vector has 6 entries, corresponding to the xx, yy, zz, xy, xz, yz components +of the symmetric strain tensor. +
+Restrictions: +
+This compute is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. +
+This quantity will be computed only for particles which interact with the Total-Lagrangian SPH pair style. +
+Related commands: +
+smd/tlsph_strain_rate, smd/tlsph_stress +
+Default: none +
+ diff --git a/doc/doc2/compute_smd_tlsph_strain_rate.html b/doc/doc2/compute_smd_tlsph_strain_rate.html new file mode 100644 index 0000000000..7705697208 --- /dev/null +++ b/doc/doc2/compute_smd_tlsph_strain_rate.html @@ -0,0 +1,57 @@ + +Syntax: +
+compute ID group-ID smd/tlsph_strain_rate ++
Examples: +
+compute 1 all smd/tlsph_strain_rate ++
Description: +
+Define a computation that calculates the rate of the strain tensor for particles interacting via the Total-Lagrangian SPH pair style. +
+See this PDF guide to using Smooth Mach Dynamics in LAMMPS. +
+Output info: +
+This compute calculates a per-particle vector of vectors (tensors), which can be +accessed by any command that uses per-particle values from a compute as input. See +How-to discussions, section 6.15 +for an overview of LAMMPS output options. +
+The values will be given in units of one over time. +
+The per-particle vector has 6 entries, corresponding to the xx, yy, zz, xy, xz, yz components +of the symmetric strain rate tensor. +
+Restrictions: +
+This compute is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. +
+This quantity will be computed only for particles which interact with Total-Lagrangian SPH pair style. +
+Related commands: +
+smd/tlsph_strain, smd/tlsph_stress +
+Default: none +
+ diff --git a/doc/doc2/compute_smd_tlsph_stress.html b/doc/doc2/compute_smd_tlsph_stress.html new file mode 100644 index 0000000000..3f22be60cd --- /dev/null +++ b/doc/doc2/compute_smd_tlsph_stress.html @@ -0,0 +1,58 @@ + +Syntax: +
+compute ID group-ID smd/tlsph_stress ++
Examples: +
+compute 1 all smd/tlsph_stress ++
Description: +
+Define a computation that outputs the Cauchy stress tensor for particles interacting via the Total-Lagrangian SPH pair style. +
+See this PDF guide to using Smooth Mach Dynamics in LAMMPS. +
+Output info: +
+This compute calculates a per-particle vector of vectors (tensors), which can be +accessed by any command that uses per-particle values from a compute as input. See +How-to discussions, section 6.15 +for an overview of LAMMPS output options. +
+The values will be given in units of pressure. +
+The per-particle vector has 7 entries. The first six entries correspond to the xx, yy, zz, xy, xz and yz components +of the symmetric Cauchy stress tensor. The seventh entry is the second invariant of the stress tensor, i.e., +the von Mises equivalent stress. +
+Restrictions: +
+This compute is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. +
+This quantity will be computed only for particles which interact with the Total-Lagrangian SPH pair style. +
+Related commands: +
+smd/tlsph_strain, smd/tlsph_strain_rate +
+Default: none +
+ diff --git a/doc/doc2/compute_smd_triangle_mesh_vertices.html b/doc/doc2/compute_smd_triangle_mesh_vertices.html new file mode 100644 index 0000000000..a4fc72d701 --- /dev/null +++ b/doc/doc2/compute_smd_triangle_mesh_vertices.html @@ -0,0 +1,64 @@ + +Syntax: +
+compute ID group-ID smd/triangle_mesh_vertices ++
Examples: +
+compute 1 all smd/triangle_mesh_vertices ++
Description: +
+Define a computation that returns the coordinates of the vertices +corresponding to the triangle-elements of a mesh created by the fix smd/wall_surface. +
+See this PDF guide to using Smooth Mach Dynamics in LAMMPS. +
+Output info: +
+This compute returns a per-particle vector of vectors, which can be +accessed by any command that uses per-particle values from a compute as +input. See How-to discussions, section 6.15 +for an overview of LAMMPS output options. +
+The per-particle vector has nine entries, (x1/y1/z1), (x2/y2/z2), and (x3/y3/z3) corresponding +to the first, second, and third vertex of each triangle. +
+It is only meaningful to use this compute for a group of particles which is created via the +fix smd/wall_surface command. +
+The output of this compute can be used with the dump2vtk_tris tool to generate a VTK representation of the +smd/wall_surace mesh for visualization purposes. +
+The values will be given in units of distance. +
+Restrictions: +
+This compute is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. +
+Related commands: +
+smd/move_tri_surf, smd/wall_surface +
+Default: none +
+ diff --git a/doc/doc2/compute_smd_ulsph_num_neighs.html b/doc/doc2/compute_smd_ulsph_num_neighs.html new file mode 100644 index 0000000000..e3e1a43e85 --- /dev/null +++ b/doc/doc2/compute_smd_ulsph_num_neighs.html @@ -0,0 +1,53 @@ + +Syntax: +
+compute ID group-ID smd/ulsph_num_neighs ++
Examples: +
+compute 1 all smd/ulsph_num_neighs ++
Description: +
+Define a computation that returns the number of neighbor particles +inside of the smoothing kernel radius for particles interacting via the updated Lagrangian SPH pair style. +
+See this PDF guide to using Smooth Mach Dynamics in LAMMPS. +
+Output info: +
+This compute returns a per-particle vector, which can be accessed by +any command that uses per-particle values from a compute as input. See +Section_howto 15 for an overview of +LAMMPS output options. +
+The per-particle values will be given dimentionless, see units. +
+Restrictions: +
+This compute is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. This compute can only be used for particles which interact with the updated Lagrangian SPH pair style. +
+Related commands: +
+smd/tlsph_num_neighs +
+Default: none +
+ diff --git a/doc/doc2/compute_smd_ulsph_strain.html b/doc/doc2/compute_smd_ulsph_strain.html new file mode 100644 index 0000000000..8aca12a2a0 --- /dev/null +++ b/doc/doc2/compute_smd_ulsph_strain.html @@ -0,0 +1,56 @@ + +Syntax: +
+compute ID group-ID smd/ulsph_strain ++
Examples: +
+compute 1 all smd/ulsph_strain ++
Description: +
+Define a computation that outputs the logarithmic strain tensor. +for particles interacting via the updated Lagrangian SPH pair style. +
+See this PDF guide to using Smooth Mach Dynamics in LAMMPS. +
+Output info: +
+This compute calculates a per-particle tensor, which can be accessed by +any command that uses per-particle values from a compute as input. See +Section_howto 15 for an overview of +LAMMPS output options. +
+The per-particle vector has 6 entries, corresponding to the xx, yy, zz, xy, xz, yz components +of the symmetric strain rate tensor. +
+The per-particle tensor values will be given dimensionless, see units. +
+Restrictions: +
+This compute is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. This compute can only be used for particles which interact with the updated Lagrangian SPH pair style. +
+Related commands: +
+smd/tlsph_strain +
+Default: none +
+ diff --git a/doc/doc2/compute_smd_ulsph_strain_rate.html b/doc/doc2/compute_smd_ulsph_strain_rate.html new file mode 100644 index 0000000000..85015d69f5 --- /dev/null +++ b/doc/doc2/compute_smd_ulsph_strain_rate.html @@ -0,0 +1,56 @@ + +Syntax: +
+compute ID group-ID smd/ulsph_strain_rate ++
Examples: +
+compute 1 all smd/ulsph_strain_rate ++
Description: +
+Define a computation that outputs the rate of the logarithmic strain tensor +for particles interacting via the updated Lagrangian SPH pair style. +
+See this PDF guide to using Smooth Mach Dynamics in LAMMPS. +
+Output info: +
+This compute calculates a per-particle vector of vectors (tensors), which +can be accessed by any command that uses per-particle values from a compute +as input. See Section_howto 15 for an overview of +LAMMPS output options. +
+The values will be given in units of one over time. +
+The per-particle vector has 6 entries, corresponding to the xx, yy, zz, xy, xz, yz components +of the symmetric strain rate tensor. +
+Restrictions: +
+This compute is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. This compute can only be used for particles which interact with the updated Lagrangian SPH pair style. +
+Related commands: +
+smd/utsph_strain_rate +
+Default: none +
+ diff --git a/doc/doc2/compute_smd_ulsph_stress.html b/doc/doc2/compute_smd_ulsph_stress.html new file mode 100644 index 0000000000..0af4d67087 --- /dev/null +++ b/doc/doc2/compute_smd_ulsph_stress.html @@ -0,0 +1,57 @@ + +Syntax: +
+compute ID group-ID smd/ulsph_stress ++
Examples: +
+compute 1 all smd/ulsph_stress ++
Description: +
+Define a computation that outputs the Cauchy stress tensor. +
+See this PDF guide to using Smooth Mach Dynamics in LAMMPS. +
+Output info: +
+This compute calculates a per-particle vector of vectors (tensors), +which can be accessed by any command that uses per-particle values +from a compute as input. See Section_howto 15 +for an overview of LAMMPS output options. +
+The values will be given in units of pressure. +
+The per-particle vector has 7 entries. The first six entries correspond to the xx, yy, zz, xy, xz, yz components +of the symmetric Cauchy stress tensor. The seventh entry is the second invariant of the stress tensor, i.e., +the von Mises equivalent stress. +
+Restrictions: +
+This compute is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. This compute can only be used for particles which interact with the updated Lagrangian SPH pair style. +
+Related commands: +
+smd/ulsph_strain, smd/ulsph_strain_rate +smd/tlsph_stress ++
Default: none +
+ diff --git a/doc/doc2/compute_smd_vol.html b/doc/doc2/compute_smd_vol.html new file mode 100644 index 0000000000..d0a1c70911 --- /dev/null +++ b/doc/doc2/compute_smd_vol.html @@ -0,0 +1,57 @@ + +Syntax: +
+compute ID group-ID smd/vol ++
Examples: +
+compute 1 all smd/vol ++
Description: +
+Define a computation that provides the per-particle volume and the sum of the per-particle volumes of the +group for which the fix is defined. +
+See this PDF guide to using Smooth Mach Dynamics in +LAMMPS. +
+Output info: +
+This compute calculates a per-particle vector, which can be accessed by +any command that uses per-particle values from a compute as input. See +How-to discussions, section 6.15 +for an overview of LAMMPS output options. +
+The per-particle vector values will be given in units of volume. +
+Additionally, the compute returns a scalar, which is the sum of the per-particle volumes of the +group for which the fix is defined. +
+Restrictions: +
+This compute is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS + section for more info. +
+Related commands: +
+smd/rho +
+Default: none +
+ diff --git a/doc/doc2/fix_gcmc.html b/doc/doc2/fix_gcmc.html index 7a02e5b4b1..e01d624782 100644 --- a/doc/doc2/fix_gcmc.html +++ b/doc/doc2/fix_gcmc.html @@ -148,13 +148,9 @@ coordinates are chosen at a random position within the current simulation cell or region, and new atom velocities are randomly chosen from the specified temperature distribution given by T. Relative coordinates for atoms in a molecule are taken from the template -molecule provided by the user, with the origin of the relative -coordinates coinciding with the chosen insertion point. This means -that if the origin of the template molecule coordinate system -lies far from the center of the template molecule, -the inserted molecule will lie far from the insertion point. -A random initial rotation is used in -the case of molecule insertions. +molecule provided by the user. The center of mass of the molecule +is placed at the insertion point. The orientation of the molecule +is chosen at random by rotating about this point.Individual atoms are inserted, unless the mol keyword is used. It specifies a template-ID previously defined using the @@ -182,11 +178,9 @@ the angle in degrees. Rotations are performed by generating a random point on the unit sphere and a random rotation angle on the range [0,maxangle). The molecule is then rotated by that angle about an axis passing through the molecule center of mass. The axis is parallel -to the unit vector defined by the point on the unit sphere. The -same procedure is used for randomly rotating molecules when they -are inserted, except that the rotation axis passes through whatever -origin is used for the molecule template, and the maximum angle is -360 degrees. +to the unit vector defined by the point on the unit sphere. +The same procedure is used for randomly rotating molecules when they +are inserted, except that the maximum angle is 360 degrees.
Note that fix GCMC does not use configurational bias MC or any other kind of sampling of intramolecular degrees of freedom. diff --git a/doc/doc2/fix_smd_adjust_dt.html b/doc/doc2/fix_smd_adjust_dt.html new file mode 100644 index 0000000000..e7fab431b1 --- /dev/null +++ b/doc/doc2/fix_smd_adjust_dt.html @@ -0,0 +1,63 @@ + +
Syntax: +
+fix ID group-ID smd/adjust_dt arg ++
s_fact = safety factor ++ +
Examples: +
+fix 1 all smd/adjust_dt 0.1 ++
Description: +
+The fix calculates a new stable time increment for use with the SMD time integrators. +
+The stable time increment is based on multiple conditions. For the SPH pair styles, a +CFL criterion (Courant, Friedrichs & Lewy, 1928) is evaluated, which determines the the speed of +sound cannot propagate further than a typical spacing between particles within a single time step to ensure +no information is lost. For the contact pair styles, a linear analysis of the pair potential determines a +stable maximum time step. +
+This fix inquires the minimum stable time increment across all particles contained in the group for which this +fix is defined. An additional safety factor s_fact is applied to the time increment. +
+See this PDF guide to use Smooth Mach Dynamics in LAMMPS. +
+Restart, fix_modify, output, run start/stop, minimize info: +
+Currently, no part of USER-SMD supports restarting nor minimization. +
+Restrictions: +
+This fix is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. +
+Related commands: +
+ +Default: none +
+ diff --git a/doc/doc2/fix_smd_integrate_tlsph.html b/doc/doc2/fix_smd_integrate_tlsph.html new file mode 100644 index 0000000000..f6a89384b4 --- /dev/null +++ b/doc/doc2/fix_smd_integrate_tlsph.html @@ -0,0 +1,59 @@ + +Syntax: +
+fix ID group-ID smd/integrate_tlsph keyword values ++
limit_velocity value = max_vel + max_vel = maximum allowed velocity ++ +
Examples: +
+fix 1 all smd/integrate_tlsph +fix 1 all smd/integrate_tlsph limit_velocity 1000 ++
Description: +
+The fix performs explicit time integration for particles which interact according with the Total-Lagrangian SPH pair style. +
+See this PDF guide to using Smooth Mach Dynamics in LAMMPS. +
+The limit_velocity keyword will control the velocity, scaling the norm of +the velocity vector to max_vel in case it exceeds this velocity limit. +
+Restart, fix_modify, output, run start/stop, minimize info: +
+Currently, no part of USER-SMD supports restarting nor minimization. This fix has no outputs. +
+Restrictions: +
+This fix is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. +
+Related commands: +
+ +Default: none +
+ diff --git a/doc/doc2/fix_smd_integrate_ulsph.html b/doc/doc2/fix_smd_integrate_ulsph.html new file mode 100644 index 0000000000..0ea578c2fc --- /dev/null +++ b/doc/doc2/fix_smd_integrate_ulsph.html @@ -0,0 +1,63 @@ + +Syntax: +
+fix ID group-ID smd/integrate_ulsph keyword ++
keyword = adjust_radius or limit_velocity +
+adjust_radius values = adjust_radius_factor min_nn max_nn + adjust_radius_factor = factor which scale the smooth/kernel radius + min_nn = minimum number of neighbors + max_nn = maximum number of neighbors +limit_velocity values = max_velocity + max_velocity = maximum allowed velocity. +
+Examples: +
+fix 1 all smd/integrate_ulsph adjust_radius 1.02 25 50 ++
fix 1 all smd/integrate_ulsph limit_velocity 1000 ++
Description: +
+The fix performs explicit time integration for particles which interact with the updated Lagrangian SPH pair style. +See this PDF guide to using Smooth Mach Dynamics in LAMMPS. +
+The adjust_radius keyword activates dynamic adjustment of the per-particle SPH smoothing kernel radius such that the number of neighbors per particles remains +within the interval min_nn to max_nn. The parameter adjust_radius_factor determines the amount of adjustment per timestep. Typical values are +adjust_radius_factor=1.02, min_nn=15, and max_nn=20. +
+The limit_velocity keyword will control the velocity, scaling the norm of +the velocity vector to max_vel in case it exceeds this velocity limit. +
+Restart, fix_modify, output, run start/stop, minimize info: +
+Currently, no part of USER-SMD supports restarting nor minimization. This fix has no outputs. +
+Restrictions: +
+This fix is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. +
+Related commands: +
+Default: none +
+ diff --git a/doc/doc2/fix_smd_move_triangulated_surface.html b/doc/doc2/fix_smd_move_triangulated_surface.html new file mode 100644 index 0000000000..e3d1cea464 --- /dev/null +++ b/doc/doc2/fix_smd_move_triangulated_surface.html @@ -0,0 +1,81 @@ + +Syntax: +
+fix ID group-ID smd/move_tri_surf keyword ++
*LINEAR args = Vx Vy Vz + Vx,Vy,Vz = components of velocity vector (velocity units), any component can be specified as NULL + *WIGGLE args = Vx Vy Vz max_travel + vx,vy,vz = components of velocity vector (velocity units), any component can be specified as NULL + max_travel = wiggle amplitude + *ROTATE args = Px Py Pz Rx Ry Rz period + Px,Py,Pz = origin point of axis of rotation (distance units) + Rx,Ry,Rz = axis of rotation vector + period = period of rotation (time units) ++ +
Examples: +
+fix 1 tool smd/move_tri_surf *LINEAR 20 20 10 +fix 2 tool smd/move_tri_surf *WIGGLE 20 20 10 +fix 2 tool smd/move_tri_surf *ROTATE 0 0 0 5 2 1 ++
Description: +
+This fix applies only to rigid surfaces read from .STL files via fix smd/wall_surface . +It updates position and velocity for the particles in the group each timestep without regard to forces on the particles. +The rigid surfaces can thus be moved along simple trajectories during the simulation. +
+The *LINEAR style moves particles with the specified constant velocity +vector V = (Vx,Vy,Vz). This style also sets the velocity of each particle +to V = (Vx,Vy,Vz). +
+The *WIGGLE style moves particles in an oscillatory fashion. +Particles are moved along (vx, vy, vz) with constant velocity until a +displacement of max_travel is reached. Then, the velocity vector is +reversed. This process is repeated. +
+The *ROTATE style rotates particles around a rotation axis R = (Rx,Ry,Rz) that +goes through a point P = (Px,Py,Pz). The period of the rotation is also +specified. This style also sets the velocity of each particle to (omega cross +Rperp) where omega is its angular velocity around the rotation axis and +Rperp is a perpendicular vector from the rotation axis to the particle. +
+See this PDF guide to using Smooth Mach Dynamics in LAMMPS. +
+Restart, fix_modify, output, run start/stop, minimize info: +
+Currently, no part of USER-SMD supports restarting nor minimization. This fix has no outputs. +
+Restrictions: +
+This fix is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. +
+Related commands: +
+smd/triangle_mesh_vertices, smd/wall_surface +
+Default: none +
+ diff --git a/doc/doc2/fix_smd_setvel.html b/doc/doc2/fix_smd_setvel.html new file mode 100644 index 0000000000..02124aa05a --- /dev/null +++ b/doc/doc2/fix_smd_setvel.html @@ -0,0 +1,96 @@ + +Syntax: +
+fix ID group-ID smd/setvel vx vy vz keyword value ... ++
region value = region-ID + region-ID = ID of region particles must be in to have their velocities set ++ +
Examples: +
+fix top_velocity top_group setvel 1.0 0.0 0.0 ++
Description: +
+Set each component of velocity on each particle in the group to the specified +values vx,vy,vz, regardless of the forces acting on the particle. This command can +be used to impose velocity boundary conditions. +
+Any of the vx,vy,vz values can be specified as NULL which means do not +alter the velocity component in that dimension. +
+This fix is indented to be used together with a time integration fix. +
+Any of the 3 quantities defining the velocity components can be specified +as an equal-style or atom-style variable, namely vx, +vy, vz. If the value is a variable, it should be specified as +v_name, where name is the variable name. In this case, the variable +will be evaluated each timestep, and its value used to determine the +force component. +
+Equal-style variables can specify formulas with various mathematical +functions, and include thermo_style command +keywords for the simulation box parameters and timestep and elapsed +time. Thus it is easy to specify a time-dependent velocity field. +
+Atom-style variables can specify the same formulas as equal-style +variables but can also include per-atom values, such as atom +coordinates. Thus it is easy to specify a spatially-dependent velocity +field with optional time-dependence as well. +
+If the region keyword is used, the particle must also be in the +specified geometric region in order to have its velocity set by this command. +
+Restart, fix_modify, output, run start/stop, minimize info: +
+Currently, no part of USER-SMD supports restarting nor minimization +None of the fix_modify options +are relevant to this fix. +
+This fix computes a global 3-vector of forces, which can be accessed +by various output commands. This is the +total force on the group of atoms. The vector values calculated by this +fix are "extensive". +
+No parameter of this fix can be used with the start/stop keywords of +the run command. +
+Restrictions: +
+This fix is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. +
+Related commands: none +
+Default: none +
+ diff --git a/doc/doc2/fix_smd_tlsph_reference_configuration.html b/doc/doc2/fix_smd_tlsph_reference_configuration.html new file mode 100644 index 0000000000..a0fa4f0f02 --- /dev/null +++ b/doc/doc2/fix_smd_tlsph_reference_configuration.html @@ -0,0 +1,2 @@ + + diff --git a/doc/doc2/fix_smd_wall_surface.html b/doc/doc2/fix_smd_wall_surface.html new file mode 100644 index 0000000000..8ddb59372c --- /dev/null +++ b/doc/doc2/fix_smd_wall_surface.html @@ -0,0 +1,68 @@ + +Syntax: +
+fix ID group-ID smd/wall_surface arg type mol-ID ++
file = file name of a triangular mesh in stl format ++
Examples: +
+fix stl_surf all smd/wall_surface tool.stl 2 65535 ++
Description: +
+This fix creates reads a traingulated surface from a file in .STL format. +For each triangle, a new particle is created which stores the barycenter of the triangle and the vertex positions. +The radius of the new particle is that of the minimum circle which encompasses the triangle vertices. +
+The triangulated surface can be used as a complex rigid wall via the smd/tri_surface pair style. +It is possible to move the triangulated surface via the smd/move_tri_surf fix style. +
+Immediately after a .STL file has been read, the simulation needs to be run for 0 timesteps in order to properly register the new particles +in the system. See the "funnel_flow" example in the USER-SMD examples directory. +
+See this PDF guide to use Smooth Mach Dynamics in LAMMPS. +
+Restart, fix_modify, output, run start/stop, minimize info: +
+Currently, no part of USER-SMD supports restarting nor minimization. This fix has no outputs. +
+Restrictions: +
+This fix is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. The molecule ID given to the particles created by this fix have to be equal to or larger than 65535. +
+Within each .STL file, only a single triangulated object must be present, even though the STL format allows for the possibility of multiple objects in one file. +
+Related commands: +
+smd/triangle_mesh_vertices, smd/move_tri_surf, smd/tri_surface +
+Default: none +
+ diff --git a/doc/doc2/pair_smd_hertz.html b/doc/doc2/pair_smd_hertz.html new file mode 100644 index 0000000000..8a7f29107c --- /dev/null +++ b/doc/doc2/pair_smd_hertz.html @@ -0,0 +1,60 @@ + +Syntax: +
+pair_style smd/hertz scale_factor ++
Examples: +
+pair_style smd/hertz 1.0
+pair_coeff 1 1
Description: +
+The smd/hertz style calculates contact forces between SPH particles belonging to different physical bodies. +
+The contact forces are calculated using a Hertz potential, which evaluates the overlap between two particles
+(whose spatial extents are defined via its contact radius).
+The effect is that a particles cannot penetrate into each other.
+The parameter
The parameter scale_factor can be used to scale the particles' contact radii. This can be useful to control how close +particles can approach each other. Usually, scale_factor=1.0. +
+Mixing, shift, table, tail correction, restart, rRESPA info: +
+No mixing is performed automatically. +Currently, no part of USER-SMD supports restarting nor minimization. +rRESPA does not apply to this pair style. +
+Restrictions: +
+This fix is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. +
+Related commands: +
+ +Default: none +
+Syntax: +
+pair_style smd/tlsph args ++
Examples: +
+pair_style smd/tlsph +
+Description: +
+The smd/tlsph style computes particle interactions according to continuum mechanics constitutive laws and a Total-Lagrangian Smooth-Particle Hydrodynamics algorithm. +
+This pair style is invoked with the following command: +
+pair_style smd/tlsph +pair_coeff i j *COMMON rho0 E nu Q1 Q2 hg Cp & + *END ++
Here, i and j denote the LAMMPS particle types for which this pair style is +defined. Note that i and j must be equal, i.e., no tlsph cross interactions +between different particle types are allowed. +In contrast to the usual LAMMPS pair coeff definitions, which are given solely a +number of floats and integers, the tlsph pair coeff definition is organised using +keywords. These keywords mark the beginning of different sets of parameters for particle properties, +material constitutive models, and damage models. The pair coeff line must be terminated with +the *END keyword. The use the line continuation operator & is recommended. A typical +invocation of the tlsph for a solid body would consist of an equation of state for computing +the pressure (the diagonal components of the stress tensor), and a material model to compute shear +stresses (the off-diagonal components of the stress tensor). Damage and failure models can also be added. +
+Please see the SMD user guide for a complete listing of the possible keywords and material models. +
+Mixing, shift, table, tail correction, restart, rRESPA info: +
+No mixing is performed automatically. +Currently, no part of USER-SMD supports restarting nor minimization. +rRESPA does not apply to this pair style. +
+Restrictions: +
+This fix is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. +
+Related commands: +
+ +Default: none +
+Syntax: +
+pair_style smd/tri_surface scale_factor ++
Examples: +
+pair_style smd/tri_surface 1.0
+pair_coeff 1 1
Description: +
+The smd/tri_surface style calculates contact forces between SPH particles and a rigid wall boundary defined via the +smd/wall_surface fix. +
+The contact forces are calculated using a Hertz potential, which evaluates the overlap between a particle
+(whose spatial extents are defined via its contact radius) and the triangle.
+The effect is that a particle cannot penetrate into the triangular surface.
+The parameter
The parameter scale_factor can be used to scale the particles' contact radii. This can be useful to control how close +particles can approach the triangulated surface. Usually, scale_factor=1.0. +
+Mixing, shift, table, tail correction, restart, rRESPA info: +
+No mixing is performed automatically. +Currently, no part of USER-SMD supports restarting nor minimization. +rRESPA does not apply to this pair style. +
+Restrictions: +
+This fix is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. +
+Related commands: +
+ +Default: none +
+Syntax: +
+pair_style smd/ulsph args ++
keyword = *DENSITY_SUMMATION or *DENSITY_CONTINUITY and *VELOCITY_GRADIENT or *NO_VELOCITY_GRADIENT and *GRADIENT_CORRECTION or *NO_GRADIENT_CORRECTION +
+Examples: +
+pair_style smd/ulsph *DENSITY_CONTINUITY *VELOCITY_GRADIENT *NO_GRADIENT_CORRECTION +
+Description: +
+The smd/ulsph style computes particle interactions according to continuum mechanics constitutive laws and an updated Lagrangian Smooth-Particle Hydrodynamics algorithm. +
+This pair style is invoked similar to the following command: +
+pair_style smd/ulsph *DENSITY_CONTINUITY *VELOCITY_GRADIENT *NO_GRADIENT_CORRECTION +pair_coeff i j *COMMON rho0 c0 Q1 Cp hg & + *END ++
Here, i and j denote the LAMMPS particle types for which this pair style is +defined. Note that i and j can be different, i.e., ulsph cross interactions +between different particle types are allowed. However, i--i respectively j--j pair_coeff lines have to preceed a cross interaction. +In contrast to the usual LAMMPS pair coeff definitions, which are given solely a +number of floats and integers, the ulsph pair coeff definition is organised using +keywords. These keywords mark the beginning of different sets of parameters for particle properties, +material constitutive models, and damage models. The pair coeff line must be terminated with +the *END keyword. The use the line continuation operator & is recommended. A typical +invocation of the ulsph for a solid body would consist of an equation of state for computing +the pressure (the diagonal components of the stress tensor), and a material model to compute shear +stresses (the off-diagonal components of the stress tensor). +
+Note that the use of *GRADIENT_CORRECTION can lead to severe numerical instabilities. For a general fluid simulation, *NO_GRADIENT_CORRECTION is recommended. +
+Please see the SMD user guide for a complete listing of the possible keywords and material models. +
+Mixing, shift, table, tail correction, restart, rRESPA info: +
+No mixing is performed automatically. +Currently, no part of USER-SMD supports restarting nor minimization. +rRESPA does not apply to this pair style. +
+Restrictions: +
+This fix is part of the USER-SMD package. It is only enabled if +LAMMPS was built with that package. See the Making LAMMPS +section for more info. +
+Related commands: +
+ +Default: none +
+Default:
The option defaults are min = 0.1 0.1 40 50, no temp setting, vel = -geom gaussian, and time = step. +geom gaussian, and time = steps.
The units for these quantities depend on the unit style; see the units command for details. diff --git a/doc/doc2/set.html b/doc/doc2/set.html index e314b1f62c..1d5aa86924 100644 --- a/doc/doc2/set.html +++ b/doc/doc2/set.html @@ -23,7 +23,7 @@
type value = atom type value can be an atom-style variable (see below) @@ -83,10 +83,12 @@ value can be an atom-style variable (see below) meso_rho value = density of SPH particles (need units) value can be an atom-style variable (see below) - i_name value = value for custom integer vector with name - value can be an atom-style variable (see below) - d_name value = value for custom floating-point vector with name + smd_contact_radius = radius for short range interactions, i.e. contact and friction value can be an atom-style variable (see below) + smd_mass_density = set particle mass based on volume by providing a mass density + value can be an atom-style variable (see below) + i_name value = value for custom integer vector with name + d_name value = value for custom floating-point vector with name@@ -316,6 +318,20 @@ assumed to be in mass/distance^2 units). If none of these cases are valid, then the mass is set to the density value directly (the input density is assumed to be in mass units). +
Keyword smd_mass_density also sets the mass of all selected +particles, but it is only applicable to the Smooth Mach Dynamics +package USER-SMD. This command assumes that the particle volume has +already been correctly set and calculates particle mass from the +provided mass density value. +
+Keyword smd_contact_radius only applies to simulations with the +Smooth Mach Dynamics package USER-SMD. This command sets an +interaction radius for computing short-range interactions, +e.g. repulsive forces to prevent different individual physical bodies +from penetrating each other. Note that the SPH smoothing kernel +diameter used for computing long range, nonlocal interactions, is set +using the diameter keyword. +
Keyword volume sets the volume of all selected particles. Currently, only the atom_style peri command defines particles with a volume attribute. Note that this command does not @@ -331,8 +347,7 @@ the simulation. The flags can be output with atom snapshots via the nx,ny,nz, then the current image value for that dimension is unchanged. For non-periodic dimensions only a value of 0 can be specified. This keyword does not allow use of atom-style variables. -
-This command can be useful after a system has been equilibrated and +This command can be useful after a system has been equilibrated and atoms have diffused one or more box lengths in various directions. This command can then reset the image values for atoms so that they are effectively inside the simulation box, e.g if a diffusion