jianmu
/
lammps
forked from lijiext/lammps
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

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

This commit is contained in:
sjplimp 2011-08-27 20:09:27 +00:00
parent fe1132af69
commit 6ce9c2bb67
24 changed files with 1323 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
doc/Section_commands.html
View File

@ -354,8 +354,9 @@ of each style or click on the style itself for a full description:
package</A>.
</P>
<DIV ALIGN=center><TABLE BORDER=1 >
<TR ALIGN="center"><TD ><A HREF = "fix_atc.html">atc</A></TD><TD ><A HREF = "fix_imd.html">imd</A></TD><TD ><A HREF = "fix_langevin_eff.html">langevin/eff</A></TD><TD ><A HREF = "fix_nh_eff.html">nph/eff</A></TD><TD ><A HREF = "fix_nh_eff.html">npt/eff</A></TD><TD ><A HREF = "fix_nve_eff.html">nve/eff</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "fix_nh_eff.html">nvt/eff</A></TD><TD ><A HREF = "fix_nvt_sllod_eff.html">nvt/sllod/eff</A></TD><TD ><A HREF = "fix_qeq_reax.html">qeq/reax</A></TD><TD ><A HREF = "fix_smd.html">smd</A></TD><TD ><A HREF = "fix_temp_rescale_eff.html">temp/rescale/eff</A>
<TR ALIGN="center"><TD ><A HREF = "fix_atc.html">atc</A></TD><TD ><A HREF = "fix_imd.html">imd</A></TD><TD ><A HREF = "fix_langevin_eff.html">langevin/eff</A></TD><TD ><A HREF = "fix_meso.html">meso</A></TD><TD ><A HREF = "fix_meso_stationary.html">meso/stationary</A></TD><TD ><A HREF = "fix_nh_eff.html">nph/eff</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "fix_nh_eff.html">npt/eff</A></TD><TD ><A HREF = "fix_nve_eff.html">nve/eff</A></TD><TD ><A HREF = "fix_nh_eff.html">nvt/eff</A></TD><TD ><A HREF = "fix_nvt_sllod_eff.html">nvt/sllod/eff</A></TD><TD ><A HREF = "fix_qeq_reax.html">qeq/reax</A></TD><TD ><A HREF = "fix_smd.html">smd</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "fix_temp_rescale_eff.html">temp/rescale/eff</A>
</TD></TR></TABLE></DIV>
<P>These are accelerated fix styles, which can be used if LAMMPS is
@ -391,8 +392,8 @@ each style or click on the style itself for a full description:
package</A>.
</P>
<DIV ALIGN=center><TABLE BORDER=1 >
<TR ALIGN="center"><TD ><A HREF = "compute_ackland_atom.html">ackland/atom</A></TD><TD ><A HREF = "compute_ke_eff.html">ke/eff</A></TD><TD ><A HREF = "compute_ke_atom_eff.html">ke/atom/eff</A></TD><TD ><A HREF = "compute_temp_eff.html">temp/eff</A></TD><TD ><A HREF = "compute_temp_deform_eff.html">temp/deform/eff</A></TD><TD ><A HREF = "compute_temp_region_eff.html">temp/region/eff</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "compute_temp_rotate.html">temp/rotate</A>
<TR ALIGN="center"><TD ><A HREF = "compute_ackland_atom.html">ackland/atom</A></TD><TD ><A HREF = "compute_ke_eff.html">ke/eff</A></TD><TD ><A HREF = "compute_ke_atom_eff.html">ke/atom/eff</A></TD><TD ><A HREF = "compute_meso_e_atom.html">meso_e/atom</A></TD><TD ><A HREF = "compute_meso_rho_atom.html">meso_rho/atom</A></TD><TD ><A HREF = "compute_meso_t_atom.html">meso_t/atom</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "compute_temp_eff.html">temp/eff</A></TD><TD ><A HREF = "compute_temp_deform_eff.html">temp/deform/eff</A></TD><TD ><A HREF = "compute_temp_region_eff.html">temp/region/eff</A></TD><TD ><A HREF = "compute_temp_rotate.html">temp/rotate</A>
</TD></TR></TABLE></DIV>
<P>These are accelerated compute styles, which can be used if LAMMPS is
@ -436,7 +437,8 @@ package</A>.
<DIV ALIGN=center><TABLE BORDER=1 >
<TR ALIGN="center"><TD ><A HREF = "pair_awpmd.html">awpmd/cut</A></TD><TD ><A HREF = "pair_buck_coul.html">buck/coul</A></TD><TD ><A HREF = "pair_cmm.html">cg/cmm</A></TD><TD ><A HREF = "pair_cmm.html">cg/cmm/coul/cut</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_cmm.html">cg/cmm/coul/long</A></TD><TD ><A HREF = "pair_dipole.html">dipole/sf</A></TD><TD ><A HREF = "pair_eam.html">eam/cd</A></TD><TD ><A HREF = "pair_eff.html">eff/cut</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_lj_coul.html">lj/coul</A></TD><TD ><A HREF = "pair_lj_sf.html">lj/sf</A></TD><TD ><A HREF = "pair_reax_c.html">reax/c</A>
<TR ALIGN="center"><TD ><A HREF = "pair_lj_coul.html">lj/coul</A></TD><TD ><A HREF = "pair_lj_sf.html">lj/sf</A></TD><TD ><A HREF = "pair_reax_c.html">reax/c</A><A HREF = "pair_heatconduction.html">sph/heatconduction</A></TD><TD ><A HREF = "pair_idealgas.html">sph/idealgas</A></TD></TR>
<TR ALIGN="center"><TD ><A HREF = "pair_lj.html">sph/lj</A></TD><TD ><A HREF = "pair_rhosum.html">sph/rhosum</A></TD><TD ><A HREF = "pair_taitwater.html">sph/taitwater</A></TD><TD ><A HREF = "pair_taitwater_morris.html">sph/taitwater/morris</A>
</TD></TR></TABLE></DIV>
<P>These are accelerated pair styles, which can be used if LAMMPS is

13
doc/Section_commands.txt
View File

@ -487,6 +487,8 @@ package"_Section_start.html#start_3.
"atc"_fix_atc.html,
"imd"_fix_imd.html,
"langevin/eff"_fix_langevin_eff.html,
"meso"_fix_meso.html,
"meso/stationary"_fix_meso_stationary.html,
"nph/eff"_fix_nh_eff.html,
"npt/eff"_fix_nh_eff.html,
"nve/eff"_fix_nve_eff.html,
@ -578,6 +580,9 @@ package"_Section_start.html#start_3.
"ackland/atom"_compute_ackland_atom.html,
"ke/eff"_compute_ke_eff.html,
"ke/atom/eff"_compute_ke_atom_eff.html,
"meso_e/atom"_compute_meso_e_atom.html,
"meso_rho/atom"_compute_meso_rho_atom.html,
"meso_t/atom"_compute_meso_t_atom.html,
"temp/eff"_compute_temp_eff.html,
"temp/deform/eff"_compute_temp_deform_eff.html,
"temp/region/eff"_compute_temp_region_eff.html,
@ -675,7 +680,13 @@ package"_Section_start.html#start_3.
"eff/cut"_pair_eff.html,
"lj/coul"_pair_lj_coul.html,
"lj/sf"_pair_lj_sf.html,
"reax/c"_pair_reax_c.html :tb(c=4,ea=c)
"reax/c"_pair_reax_c.html
"sph/heatconduction"_pair_heatconduction.html,
"sph/idealgas"_pair_idealgas.html,
"sph/lj"_pair_lj.html,
"sph/rhosum"_pair_rhosum.html,
"sph/taitwater"_pair_taitwater.html,
"sph/taitwater/morris"_pair_taitwater_morris.html :tb(c=4,ea=c)
These are accelerated pair styles, which can be used if LAMMPS is
built with the "appropriate accelerated

54
doc/compute_meso_e_atom.html Normal file
View File

@ -0,0 +1,54 @@
<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>compute meso_e/atom command
</H3>
<P><B>Syntax:</B>
</P>
<PRE>compute ID group-ID meso_e/atom
</PRE>
<UL><LI>ID, group-ID are documented in <A HREF = "compute.html">compute</A> command
<LI>meso_e/atom = style name of this compute command
</UL>
<P><B>Examples:</B>
</P>
<PRE>compute 1 all meso_e/atom
</PRE>
<P><B>Description:</B>
</P>
<P>Define a computation that calculates the per-atom internal energy
for each atom in a group.
</P>
<P>The internal energy is the energy associated with the internal degrees
of freedom of a mesoscopic particles, e.g. a Smooth-Particle
Hydrodynamics particle.
</P>
<P>The value of the internal energy will be 0.0 for atoms not in the
specified compute group.
</P>
<P><B>Output info:</B>
</P>
<P>This compute calculates a per-atom vector, which can be accessed by
any command that uses per-atom values from a compute as input. See
<A HREF = "Section_howto.html#4_15">this section</A> for an overview of LAMMPS
output options.
</P>
<P>The per-atom vector values will be in energy <A HREF = "units.html">units</A>.
</P>
<P><B>Restrictions:</B> none
</P>
<P><B>Related commands:</B>
</P>
<P><A HREF = "dump.html">dump custom</A>
</P>
<P><B>Default:</B> none
</P>
</HTML>

49
doc/compute_meso_e_atom.txt Normal file
View File

@ -0,0 +1,49 @@
"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
compute meso_e/atom command :h3
[Syntax:]
compute ID group-ID meso_e/atom :pre
ID, group-ID are documented in "compute"_compute.html command
meso_e/atom = style name of this compute command :ul
[Examples:]
compute 1 all meso_e/atom :pre
[Description:]
Define a computation that calculates the per-atom internal energy
for each atom in a group.
The internal energy is the energy associated with the internal degrees
of freedom of a mesoscopic particles, e.g. a Smooth-Particle
Hydrodynamics particle.
The value of the internal energy will be 0.0 for atoms not in the
specified compute group.
[Output info:]
This compute calculates a per-atom vector, which can be accessed by
any command that uses per-atom values from a compute as input. See
"this section"_Section_howto.html#4_15 for an overview of LAMMPS
output options.
The per-atom vector values will be in energy "units"_units.html.
[Restrictions:] none
[Related commands:]
"dump custom"_dump.html
[Default:] none

54
doc/compute_meso_rho_atom.html Normal file
View File

@ -0,0 +1,54 @@
<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>compute meso_rho/atom command
</H3>
<P><B>Syntax:</B>
</P>
<PRE>compute ID group-ID meso_rho/atom
</PRE>
<UL><LI>ID, group-ID are documented in <A HREF = "compute.html">compute</A> command
<LI>meso_rho/atom = style name of this compute command
</UL>
<P><B>Examples:</B>
</P>
<PRE>compute 1 all meso_rho/atom
</PRE>
<P><B>Description:</B>
</P>
<P>Define a computation that calculates the per-atom mesoscopic density
for each atom in a group.
</P>
<P>The mesoscopic density is the mass density of a mesoscopic particle,
calculated by kernel function interpolation using
"pair style sph/rhosum".
</P>
<P>The value of the mesoscopic density will be 0.0 for atoms not in the
specified compute group.
</P>
<P><B>Output info:</B>
</P>
<P>This compute calculates a per-atom vector, which can be accessed by
any command that uses per-atom values from a compute as input. See
<A HREF = "Section_howto.html#4_15">this section</A> for an overview of LAMMPS
output options.
</P>
<P>The per-atom vector values will be in mass / volume <A HREF = "units.html">units</A>.
</P>
<P><B>Restrictions:</B> none
</P>
<P><B>Related commands:</B>
</P>
<P><A HREF = "dump.html">dump custom</A>
</P>
<P><B>Default:</B> none
</P>
</HTML>

49
doc/compute_meso_rho_atom.txt Normal file
View File

@ -0,0 +1,49 @@
"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
compute meso_rho/atom command :h3
[Syntax:]
compute ID group-ID meso_rho/atom :pre
ID, group-ID are documented in "compute"_compute.html command
meso_rho/atom = style name of this compute command :ul
[Examples:]
compute 1 all meso_rho/atom :pre
[Description:]
Define a computation that calculates the per-atom mesoscopic density
for each atom in a group.
The mesoscopic density is the mass density of a mesoscopic particle,
calculated by kernel function interpolation using
"pair style sph/rhosum".
The value of the mesoscopic density will be 0.0 for atoms not in the
specified compute group.
[Output info:]
This compute calculates a per-atom vector, which can be accessed by
any command that uses per-atom values from a compute as input. See
"this section"_Section_howto.html#4_15 for an overview of LAMMPS
output options.
The per-atom vector values will be in mass / volume "units"_units.html.
[Restrictions:] none
[Related commands:]
"dump custom"_dump.html
[Default:] none

56
doc/compute_meso_t_atom.html Normal file
View File

@ -0,0 +1,56 @@
<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>compute meso_t/atom command
</H3>
<P><B>Syntax:</B>
</P>
<PRE>compute ID group-ID meso_t/atom
</PRE>
<UL><LI>ID, group-ID are documented in <A HREF = "compute.html">compute</A> command
<LI>meso_t/atom = style name of this compute command
</UL>
<P><B>Examples:</B>
</P>
<PRE>compute 1 all meso_t/atom
</PRE>
<P><B>Description:</B>
</P>
<P>Define a computation that calculates the per-atom internal temperature
for each atom in a group.
</P>
<P>The internal temperature is the ratio of internal energy over the heat
capacity associated with the internal degrees of freedom of a mesoscopic
particles, e.g. a Smooth-Particle Hydrodynamics particle.
</P>
<P>T_<I>int</I> = E_<I>int</I> / C_<I>V, int</I>
</P>
<P>The value of the internal energy will be 0.0 for atoms not in the
specified compute group.
</P>
<P><B>Output info:</B>
</P>
<P>This compute calculates a per-atom vector, which can be accessed by
any command that uses per-atom values from a compute as input. See
<A HREF = "Section_howto.html#4_15">this section</A> for an overview of LAMMPS
output options.
</P>
<P>The per-atom vector values will be in temperature <A HREF = "units.html">units</A>.
</P>
<P><B>Restrictions:</B> none
</P>
<P><B>Related commands:</B>
</P>
<P><A HREF = "dump.html">dump custom</A>
</P>
<P><B>Default:</B> none
</P>
</HTML>

51
doc/compute_meso_t_atom.txt Normal file
View File

@ -0,0 +1,51 @@
"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
compute meso_t/atom command :h3
[Syntax:]
compute ID group-ID meso_t/atom :pre
ID, group-ID are documented in "compute"_compute.html command
meso_t/atom = style name of this compute command :ul
[Examples:]
compute 1 all meso_t/atom :pre
[Description:]
Define a computation that calculates the per-atom internal temperature
for each atom in a group.
The internal temperature is the ratio of internal energy over the heat
capacity associated with the internal degrees of freedom of a mesoscopic
particles, e.g. a Smooth-Particle Hydrodynamics particle.
T_{int} = E_{int} / C_{V, int}
The value of the internal energy will be 0.0 for atoms not in the
specified compute group.
[Output info:]
This compute calculates a per-atom vector, which can be accessed by
any command that uses per-atom values from a compute as input. See
"this section"_Section_howto.html#4_15 for an overview of LAMMPS
output options.
The per-atom vector values will be in temperature "units"_units.html.
[Restrictions:] none
[Related commands:]
"dump custom"_dump.html
[Default:] none

49
doc/fix_meso.html Normal file
View File

@ -0,0 +1,49 @@
<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>fix meso command
</H3>
<P><B>Syntax:</B>
</P>
<PRE>fix ID group-ID meso
</PRE>
<UL><LI>ID, group-ID are documented in <A HREF = "fix.html">fix</A> command
<LI>meso = style name of this fix command
</UL>
<P><B>Examples:</B>
</P>
<PRE>fix 1 all meso
</PRE>
<P><B>Description:</B>
</P>
<P>Perform time integration to update position, velocity, internal energy
and local density for atoms in the group each timestep. This fix is needed to
time-integrate mesoscopic systems where particles carry internal variables such as SPH or DPDE.
</P>
<P><B>Restart, fix_modify, output, run start/stop, minimize info:</B>
</P>
<P>No information about this fix is written to <A HREF = "restart.html">binary restart
files</A>. None of the <A HREF = "fix_modify.html">fix_modify</A> options
are relevant to this fix. No global or per-atom quantities are stored
by this fix for access by various <A HREF = "Section_howto.html#4_15">output
commands</A>. No parameter of this fix can be
used with the <I>start/stop</I> keywords of the <A HREF = "run.html">run</A> command.
This fix is not invoked during <A HREF = "minimize.html">energy minimization</A>.
</P>
<P><B>Restrictions:</B> none
</P>
<P><B>Related commands:</B>
</P>
<P>"fix meso/stationary"
</P>
<P><B>Default:</B> none
</P>
</HTML>

44
doc/fix_meso.txt Normal file
View File

@ -0,0 +1,44 @@
"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
fix meso command :h3
[Syntax:]
fix ID group-ID meso :pre
ID, group-ID are documented in "fix"_fix.html command
meso = style name of this fix command :ul
[Examples:]
fix 1 all meso :pre
[Description:]
Perform time integration to update position, velocity, internal energy
and local density for atoms in the group each timestep. This fix is needed to
time-integrate mesoscopic systems where particles carry internal variables such as SPH or DPDE.
[Restart, fix_modify, output, run start/stop, minimize info:]
No information about this fix is written to "binary restart
files"_restart.html. None of the "fix_modify"_fix_modify.html options
are relevant to this fix. No global or per-atom quantities are stored
by this fix for access by various "output
commands"_Section_howto.html#4_15. No parameter of this fix can be
used with the {start/stop} keywords of the "run"_run.html command.
This fix is not invoked during "energy minimization"_minimize.html.
[Restrictions:] none
[Related commands:]
"fix meso/stationary"
[Default:] none

50
doc/fix_meso_stationary.html Normal file
View File

@ -0,0 +1,50 @@
<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>fix meso/stationary command
</H3>
<P><B>Syntax:</B>
</P>
<PRE>fix ID group-ID meso/stationary
</PRE>
<UL><LI>ID, group-ID are documented in <A HREF = "fix.html">fix</A> command
<LI>meso = style name of this fix command
</UL>
<P><B>Examples:</B>
</P>
<PRE>fix 1 boundary meso/stationary
</PRE>
<P><B>Description:</B>
</P>
<P>Perform time integration to update internal energy and local density, but not
position or velocity for atoms in the group each timestep. This fix is needed
for SPH simulations to correctly time-integrate fixed boundary particles which
constrain a fluid to a given region in space.
</P>
<P><B>Restart, fix_modify, output, run start/stop, minimize info:</B>
</P>
<P>No information about this fix is written to <A HREF = "restart.html">binary restart
files</A>. None of the <A HREF = "fix_modify.html">fix_modify</A> options
are relevant to this fix. No global or per-atom quantities are stored
by this fix for access by various <A HREF = "Section_howto.html#4_15">output
commands</A>. No parameter of this fix can be
used with the <I>start/stop</I> keywords of the <A HREF = "run.html">run</A> command.
This fix is not invoked during <A HREF = "minimize.html">energy minimization</A>.
</P>
<P><B>Restrictions:</B> none
</P>
<P><B>Related commands:</B>
</P>
<P>"fix meso"
</P>
<P><B>Default:</B> none
</P>
</HTML>

45
doc/fix_meso_stationary.txt Normal file
View File

@ -0,0 +1,45 @@
"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
fix meso/stationary command :h3
[Syntax:]
fix ID group-ID meso/stationary :pre
ID, group-ID are documented in "fix"_fix.html command
meso = style name of this fix command :ul
[Examples:]
fix 1 boundary meso/stationary :pre
[Description:]
Perform time integration to update internal energy and local density, but not
position or velocity for atoms in the group each timestep. This fix is needed
for SPH simulations to correctly time-integrate fixed boundary particles which
constrain a fluid to a given region in space.
[Restart, fix_modify, output, run start/stop, minimize info:]
No information about this fix is written to "binary restart
files"_restart.html. None of the "fix_modify"_fix_modify.html options
are relevant to this fix. No global or per-atom quantities are stored
by this fix for access by various "output
commands"_Section_howto.html#4_15. No parameter of this fix can be
used with the {start/stop} keywords of the "run"_run.html command.
This fix is not invoked during "energy minimization"_minimize.html.
[Restrictions:] none
[Related commands:]
"fix meso"
[Default:] none

63
doc/pair_sph_heatconduction.html Normal file
View File

@ -0,0 +1,63 @@
<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 sph/heatconduction command
</H3>
<P><B>Syntax:</B>
</P>
<PRE>pair_style sph/heatconduction
</PRE>
<P><B>Examples:</B>
</P>
<P>pair_style sph/heatconduction
pair_coeff * * 1.0 2.4
</P>
<P><B>Description:</B>
</P>
<P>The sph/heatconduction style computes heat transport between SPH particles.
The transport model is the diffusion euqation for the internal energy.
</P>
<P>The following coefficients 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.
</P>
<P>D diffusion coefficient (length^2/time units)
h kernel function cutoff (distance units)
</P>
<HR>
<P><B>Mixing, shift, table, tail correction, restart, rRESPA info</B>:
</P>
<P>These pair styles do not support mixing. Thus, coefficients for all
I,J pairs must be specified explicitly.
</P>
<P>This style does not write information to <A HREF = "restart.html">binary restart
files</A>, so pair_style and pair_coeff commands do need
to be specified in an input script that reads a restart file.
</P>
<P>These styles can only be used via the <I>pair</I> keyword of the <A HREF = "run_style.html">run_style
respa</A> command. They do not support the <I>inner</I>,
<I>middle</I>, <I>outer</I> keywords.
</P>
<P><B>Restrictions:</B>
</P>
<P>The <I>pair sph/heatconduction</I> style is part of the "USER-SPH" package. It is
only enabled if LAMMPS was built with that package (which it is by
default). 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>, pair_sph/rhosum
</P>
<P><B>Default:</B> none
</P>
</HTML>

59
doc/pair_sph_heatconduction.txt Normal file
View File

@ -0,0 +1,59 @@
"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 sph/heatconduction command :h3
[Syntax:]
pair_style sph/heatconduction :pre
[Examples:]
pair_style sph/heatconduction
pair_coeff * * 1.0 2.4
[Description:]
The sph/heatconduction style computes heat transport between SPH particles.
The transport model is the diffusion euqation for the internal energy.
The following coefficients must be defined for each pair of atoms
types via the "pair_coeff"_pair_coeff.html command as in the examples
above.
D diffusion coefficient (length^2/time units)
h kernel function cutoff (distance units)
:line
[Mixing, shift, table, tail correction, restart, rRESPA info]:
These pair styles do not support mixing. Thus, coefficients for all
I,J pairs must be specified explicitly.
This style does not write information to "binary restart
files"_restart.html, so pair_style and pair_coeff commands do need
to be specified in an input script that reads a restart file.
These styles can only be used via the {pair} keyword of the "run_style
respa"_run_style.html command. They do not support the {inner},
{middle}, {outer} keywords.
[Restrictions:]
The {pair sph/heatconduction} style is part of the "USER-SPH" package. It is
only enabled if LAMMPS was built with that package (which it is by
default). See the "Making LAMMPS"_Section_start.html#2_3 section for
more info.
[Related commands:]
"pair_coeff"_pair_coeff.html, pair_sph/rhosum
[Default:] none

68
doc/pair_sph_idealgas.html Normal file
View File

@ -0,0 +1,68 @@
<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 sph/idealgas command
</H3>
<P><B>Syntax:</B>
</P>
<PRE>pair_style sph/idealgas
</PRE>
<P><B>Examples:</B>
</P>
<P>pair_style sph/idealgas
pair_coeff * * 1.0 2.4
</P>
<P><B>Description:</B>
</P>
<P>The sph/idealgas style computes pressure forces between particles according to the ideal gas equation of state,
\begin<I>equation</I>
p = (\gamma - 1) \rho e,
\end<I>equation</I>
where $\gamma=1.4$ is the heat capacity ratio, $\rho$ is the local density, and $e$ is the internal energy per unit mass. This pair style also computes Monaghan's artificial viscosity to prevent particles from interpentrating (Monaghan1983).
</P>
<P>The following coefficients 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.
</P>
<P>nu artificial viscosity (no units)
h kernel function cutoff (distance units)
</P>
<HR>
<P><B>Mixing, shift, table, tail correction, restart, rRESPA info</B>:
</P>
<P>These pair styles do not support mixing. Thus, coefficients for all
I,J pairs must be specified explicitly.
</P>
<P>This style does not write information to <A HREF = "restart.html">binary restart
files</A>, so pair_style and pair_coeff commands do need
to be specified in an input script that reads a restart file.
</P>
<P>These styles can only be used via the <I>pair</I> keyword of the <A HREF = "run_style.html">run_style
respa</A> command. They do not support the <I>inner</I>,
<I>middle</I>, <I>outer</I> keywords.
</P>
<P><B>Restrictions:</B>
</P>
<P>The <I>pair sph/idealgas</I> style is part of the "USER-SPH" package. It is
only enabled if LAMMPS was built with that package (which it is by
default). 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>, pair_sph/rhosum
</P>
<P><B>Default:</B> none
</P>
<P><B>(Monaghan1983)</B> J. J. Monaghan and R. A. Gingold, “Shock simulation by the particle method SPH,” Journal of Computational Physics, vol. 52, no. 2, pp. 374-389, Nov. 1983.
</P>
</HTML>

65
doc/pair_sph_idealgas.txt Normal file
View File

@ -0,0 +1,65 @@
"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 sph/idealgas command :h3
[Syntax:]
pair_style sph/idealgas :pre
[Examples:]
pair_style sph/idealgas
pair_coeff * * 1.0 2.4
[Description:]
The sph/idealgas style computes pressure forces between particles according to the ideal gas equation of state,
\begin{equation}
p = (\gamma - 1) \rho e,
\end{equation}
where $\gamma=1.4$ is the heat capacity ratio, $\rho$ is the local density, and $e$ is the internal energy per unit mass. This pair style also computes Monaghan's artificial viscosity to prevent particles from interpentrating (Monaghan1983).
The following coefficients must be defined for each pair of atoms
types via the "pair_coeff"_pair_coeff.html command as in the examples
above.
nu artificial viscosity (no units)
h kernel function cutoff (distance units)
:line
[Mixing, shift, table, tail correction, restart, rRESPA info]:
These pair styles do not support mixing. Thus, coefficients for all
I,J pairs must be specified explicitly.
This style does not write information to "binary restart
files"_restart.html, so pair_style and pair_coeff commands do need
to be specified in an input script that reads a restart file.
These styles can only be used via the {pair} keyword of the "run_style
respa"_run_style.html command. They do not support the {inner},
{middle}, {outer} keywords.
[Restrictions:]
The {pair sph/idealgas} style is part of the "USER-SPH" package. It is
only enabled if LAMMPS was built with that package (which it is by
default). See the "Making LAMMPS"_Section_start.html#2_3 section for
more info.
[Related commands:]
"pair_coeff"_pair_coeff.html, pair_sph/rhosum
[Default:] none
[(Monaghan1983)] J. J. Monaghan and R. A. Gingold, “Shock simulation by the particle method SPH,” Journal of Computational Physics, vol. 52, no. 2, pp. 374-389, Nov. 1983.

72
doc/pair_sph_lj.html Normal file
View File

@ -0,0 +1,72 @@
<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 sph/lj command
</H3>
<P><B>Syntax:</B>
</P>
<PRE>pair_style sph/lj
</PRE>
<P><B>Examples:</B>
</P>
<P>pair_style sph/lj
pair_coeff * * 1.0 2.4
</P>
<P><B>Description:</B>
</P>
<P>The sph/lj style computes pressure forces between particles according to
the Lennard-Jones equation of state, which is computed according to
Ree's 1980 polynomial fit (Ree1980). The Lennard-Jones parameters
$\epsilon$ and $\sigma$ are set to unity. This pair style also computes
Monaghan's artificial viscosity to prevent particles from
interpentrating (Monaghan1983).
</P>
<P>The following coefficients 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.
</P>
<P>nu artificial viscosity (no units)
h kernel function cutoff (distance units)
</P>
<HR>
<P><B>Mixing, shift, table, tail correction, restart, rRESPA info</B>:
</P>
<P>These pair styles do not support mixing. Thus, coefficients for all
I,J pairs must be specified explicitly.
</P>
<P>This style does not write information to <A HREF = "restart.html">binary restart
files</A>, so pair_style and pair_coeff commands do need
to be specified in an input script that reads a restart file.
</P>
<P>These styles can only be used via the <I>pair</I> keyword of the <A HREF = "run_style.html">run_style
respa</A> command. They do not support the <I>inner</I>,
<I>middle</I>, <I>outer</I> keywords.
</P>
<P><B>Restrictions:</B>
The Lennard-Jones parameters $\epsilon$ and $\sigma$ are set to unity.
The <I>pair sph/lj</I> style is part of the "USER-SPH" package. It is
only enabled if LAMMPS was built with that package (which it is by
default). 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>, pair_sph/rhosum
</P>
<P><B>Default:</B> none
</P>
<P><B>(Ree1980)</B> F. H. Ree. Analytic representation of thermodynamic data for the Lennard-Jones fluid.
The Journal of Chemical Physics, 73(10):5401, 1980.
</P>
<P><B>(Monaghan1983)</B> J. J. Monaghan and R. A. Gingold, “Shock simulation by the particle method SPH,” Journal of Computational Physics, vol. 52, no. 2, pp. 374-389, Nov. 1983.
</P>
</HTML>

70
doc/pair_sph_lj.txt Normal file
View File

@ -0,0 +1,70 @@
"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 sph/lj command :h3
[Syntax:]
pair_style sph/lj :pre
[Examples:]
pair_style sph/lj
pair_coeff * * 1.0 2.4
[Description:]
The sph/lj style computes pressure forces between particles according to
the Lennard-Jones equation of state, which is computed according to
Ree's 1980 polynomial fit (Ree1980). The Lennard-Jones parameters
$\epsilon$ and $\sigma$ are set to unity. This pair style also computes
Monaghan's artificial viscosity to prevent particles from
interpentrating (Monaghan1983).
The following coefficients must be defined for each pair of atoms
types via the "pair_coeff"_pair_coeff.html command as in the examples
above.
nu artificial viscosity (no units)
h kernel function cutoff (distance units)
:line
[Mixing, shift, table, tail correction, restart, rRESPA info]:
These pair styles do not support mixing. Thus, coefficients for all
I,J pairs must be specified explicitly.
This style does not write information to "binary restart
files"_restart.html, so pair_style and pair_coeff commands do need
to be specified in an input script that reads a restart file.
These styles can only be used via the {pair} keyword of the "run_style
respa"_run_style.html command. They do not support the {inner},
{middle}, {outer} keywords.
[Restrictions:]
The Lennard-Jones parameters $\epsilon$ and $\sigma$ are set to unity.
The {pair sph/lj} style is part of the "USER-SPH" package. It is
only enabled if LAMMPS was built with that package (which it is by
default). See the "Making LAMMPS"_Section_start.html#2_3 section for
more info.
[Related commands:]
"pair_coeff"_pair_coeff.html, pair_sph/rhosum
[Default:] none
[(Ree1980)] F. H. Ree. Analytic representation of thermodynamic data for the Lennard-Jones fluid.
The Journal of Chemical Physics, 73(10):5401, 1980.
[(Monaghan1983)] J. J. Monaghan and R. A. Gingold, “Shock simulation by the particle method SPH,” Journal of Computational Physics, vol. 52, no. 2, pp. 374-389, Nov. 1983.

64
doc/pair_sph_rhosum.html Normal file
View File

@ -0,0 +1,64 @@
<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 sph/rhosum command
</H3>
<P><B>Syntax:</B>
</P>
<PRE>pair_style sph/rhosum nstep
</PRE>
<P>nstep = timestep interval
</P>
<P><B>Examples:</B>
</P>
<P>pair_style sph/rhosum 10
pair_coeff * * 2.4
</P>
<P><B>Description:</B>
</P>
<P>The sph/rhosum style computes the local particle mass density rho by kernel function interpolation,
see the SPH-USER documentation.
</P>
<P>The following coefficients 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.
</P>
<P>h (distance units)
</P>
<HR>
<P><B>Mixing, shift, table, tail correction, restart, rRESPA info</B>:
</P>
<P>These pair styles do not support mixing. Thus, coefficients for all
I,J pairs must be specified explicitly.
</P>
<P>This style does not write information to <A HREF = "restart.html">binary restart
files</A>, so pair_style and pair_coeff commands do need
to be specified in an input script that reads a restart file.
</P>
<P>These styles can only be used via the <I>pair</I> keyword of the <A HREF = "run_style.html">run_style
respa</A> command. They do not support the <I>inner</I>,
<I>middle</I>, <I>outer</I> keywords.
</P>
<P><B>Restrictions:</B>
</P>
<P>The <I>pair sph/rhosum</I> style is part of the "USER-SPH" package. It is
only enabled if LAMMPS was built with that package (which it is by
default). 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>, pair_sph/taitwater
</P>
<P><B>Default:</B> none
</P>
</HTML>

59
doc/pair_sph_rhosum.txt Normal file
View File

@ -0,0 +1,59 @@
"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 sph/rhosum command :h3
[Syntax:]
pair_style sph/rhosum nstep :pre
nstep = timestep interval
[Examples:]
pair_style sph/rhosum 10
pair_coeff * * 2.4
[Description:]
The sph/rhosum style computes the local particle mass density rho by kernel function interpolation,
see the SPH-USER documentation.
The following coefficients must be defined for each pair of atoms
types via the "pair_coeff"_pair_coeff.html command as in the examples
above.
h (distance units)
:line
[Mixing, shift, table, tail correction, restart, rRESPA info]:
These pair styles do not support mixing. Thus, coefficients for all
I,J pairs must be specified explicitly.
This style does not write information to "binary restart
files"_restart.html, so pair_style and pair_coeff commands do need
to be specified in an input script that reads a restart file.
These styles can only be used via the {pair} keyword of the "run_style
respa"_run_style.html command. They do not support the {inner},
{middle}, {outer} keywords.
[Restrictions:]
The {pair sph/rhosum} style is part of the "USER-SPH" package. It is
only enabled if LAMMPS was built with that package (which it is by
default). See the "Making LAMMPS"_Section_start.html#2_3 section for
more info.
[Related commands:]
"pair_coeff"_pair_coeff.html, pair_sph/taitwater
[Default:] none

72
doc/pair_sph_taitwater.html Normal file
View File

@ -0,0 +1,72 @@
<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 sph/taitwater command
</H3>
<P><B>Syntax:</B>
</P>
<PRE>pair_style sph/taitwater
</PRE>
<P><B>Examples:</B>
</P>
<P>pair_style sph/taitwater
pair_coeff * * 1000.0 1430.0 1.0 2.4
</P>
<P><B>Description:</B>
</P>
<P>The sph/taitwater style computes pressure forces between particles according to Tait's equation of state,
\begin<I>equation</I>
p = B <B>(\frac<I>\rho</I><I>\rho_0</I>)^<I>\gamma</I> - 1</B>,
\end<I>equation</I>
where $\gamma=7$ and $B=c_0^2 \rho_0 / \gamma$, with $\rho_0$ being the reference density and $c_0$ the reference speed of sound.
This pair style also computes Monaghan's artificial viscosity to prevent particles from interpentrating (Monaghan1983).
</P>
<P>The following coefficients 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.
</P>
<P>rho0 reference density (mass/volume units)
c0 reference soundspeed (distance/time units)
nu artificial viscosity (no units)
h kernel function cutoff (distance units)
</P>
<HR>
<P><B>Mixing, shift, table, tail correction, restart, rRESPA info</B>:
</P>
<P>These pair styles do not support mixing. Thus, coefficients for all
I,J pairs must be specified explicitly. Additionally, because rho0 and c0 are defined on a per-type basis,
an I,I as well a J,J coefficients line must be specified prior to specifying an I,J interaction.
</P>
<P>This style does not write information to <A HREF = "restart.html">binary restart
files</A>, so pair_style and pair_coeff commands do need
to be specified in an input script that reads a restart file.
</P>
<P>These styles can only be used via the <I>pair</I> keyword of the <A HREF = "run_style.html">run_style
respa</A> command. They do not support the <I>inner</I>,
<I>middle</I>, <I>outer</I> keywords.
</P>
<P><B>Restrictions:</B>
</P>
<P>The <I>pair sph/taitwater</I> style is part of the "USER-SPH" package. It is
only enabled if LAMMPS was built with that package (which it is by
default). 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>, pair_sph/rhosum
</P>
<P><B>Default:</B> none
</P>
<P><B>(Monaghan1983)</B> J. J. Monaghan and R. A. Gingold, “Shock simulation by the particle method SPH,” Journal of Computational Physics, vol. 52, no. 2, pp. 374-389, Nov. 1983.
</P>
</HTML>

69
doc/pair_sph_taitwater.txt Normal file
View File

@ -0,0 +1,69 @@
"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 sph/taitwater command :h3
[Syntax:]
pair_style sph/taitwater :pre
[Examples:]
pair_style sph/taitwater
pair_coeff * * 1000.0 1430.0 1.0 2.4
[Description:]
The sph/taitwater style computes pressure forces between particles according to Tait's equation of state,
\begin{equation}
p = B [(\frac{\rho}{\rho_0})^{\gamma} - 1],
\end{equation}
where $\gamma=7$ and $B=c_0^2 \rho_0 / \gamma$, with $\rho_0$ being the reference density and $c_0$ the reference speed of sound.
This pair style also computes Monaghan's artificial viscosity to prevent particles from interpentrating (Monaghan1983).
The following coefficients must be defined for each pair of atoms
types via the "pair_coeff"_pair_coeff.html command as in the examples
above.
rho0 reference density (mass/volume units)
c0 reference soundspeed (distance/time units)
nu artificial viscosity (no units)
h kernel function cutoff (distance units)
:line
[Mixing, shift, table, tail correction, restart, rRESPA info]:
These pair styles do not support mixing. Thus, coefficients for all
I,J pairs must be specified explicitly. Additionally, because rho0 and c0 are defined on a per-type basis,
an I,I as well a J,J coefficients line must be specified prior to specifying an I,J interaction.
This style does not write information to "binary restart
files"_restart.html, so pair_style and pair_coeff commands do need
to be specified in an input script that reads a restart file.
These styles can only be used via the {pair} keyword of the "run_style
respa"_run_style.html command. They do not support the {inner},
{middle}, {outer} keywords.
[Restrictions:]
The {pair sph/taitwater} style is part of the "USER-SPH" package. It is
only enabled if LAMMPS was built with that package (which it is by
default). See the "Making LAMMPS"_Section_start.html#2_3 section for
more info.
[Related commands:]
"pair_coeff"_pair_coeff.html, pair_sph/rhosum
[Default:] none
[(Monaghan1983)] J. J. Monaghan and R. A. Gingold, “Shock simulation by the particle method SPH,” Journal of Computational Physics, vol. 52, no. 2, pp. 374-389, Nov. 1983.

72
doc/pair_sph_taitwater_morris.html Normal file
View File

@ -0,0 +1,72 @@
<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 sph/taitwater/morris command
</H3>
<P><B>Syntax:</B>
</P>
<PRE>pair_style sph/taitwater/morris
</PRE>
<P><B>Examples:</B>
</P>
<P>pair_style sph/taitwater/morris
pair_coeff * * 1000.0 1430.0 1.0 2.4
</P>
<P><B>Description:</B>
</P>
<P>The sph/taitwater/morris style computes pressure forces between particles according to Tait's equation of state,
\begin<I>equation</I>
p = B <B>(\frac<I>\rho</I><I>\rho_0</I>)^<I>\gamma</I> - 1</B>,
\end<I>equation</I>
where $\gamma=7$ and $B=c_0^2 \rho_0 / \gamma$, with $\rho_0$ being the reference density and $c_0$ the reference speed of sound.
This pair style also computes laminar viscosity (Morris1997).
</P>
<P>The following coefficients 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.
</P>
<P>rho0 reference density (mass/volume units)
c0 reference soundspeed (distance/time units)
nu dynamic viscosity (mass*distance/time units)
h kernel function cutoff (distance units)
</P>
<HR>
<P><B>Mixing, shift, table, tail correction, restart, rRESPA info</B>:
</P>
<P>These pair styles do not support mixing. Thus, coefficients for all
I,J pairs must be specified explicitly. Additionally, because rho0 and c0 are defined on a per-type basis,
an I,I as well a J,J coefficients line must be specified prior to specifying an I,J interaction.
</P>
<P>This style does not write information to <A HREF = "restart.html">binary restart
files</A>, so pair_style and pair_coeff commands do need
to be specified in an input script that reads a restart file.
</P>
<P>These styles can only be used via the <I>pair</I> keyword of the <A HREF = "run_style.html">run_style
respa</A> command. They do not support the <I>inner</I>,
<I>middle</I>, <I>outer</I> keywords.
</P>
<P><B>Restrictions:</B>
</P>
<P>The <I>pair sph/taitwater/morris</I> style is part of the "USER-SPH" package. It is
only enabled if LAMMPS was built with that package (which it is by
default). 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>, pair_sph/rhosum
</P>
<P><B>Default:</B> none
</P>
<P><B>(Morris1997)</B> J.P. Morris, P.J. Fox, and Y. Zhu. Modeling Low Reynolds Number Incompressible Flows Using SPH. Journal of Computational Physics, 136:214–226, 1997.
</P>
</HTML>

70
doc/pair_sph_taitwater_morris.txt Normal file
View File

@ -0,0 +1,70 @@
"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 sph/taitwater/morris command :h3
[Syntax:]
pair_style sph/taitwater/morris :pre
[Examples:]
pair_style sph/taitwater/morris
pair_coeff * * 1000.0 1430.0 1.0 2.4
[Description:]
The sph/taitwater/morris style computes pressure forces between particles according to Tait's equation of state,
\begin{equation}
p = B [(\frac{\rho}{\rho_0})^{\gamma} - 1],
\end{equation}
where $\gamma=7$ and $B=c_0^2 \rho_0 / \gamma$, with $\rho_0$ being the reference density and $c_0$ the reference speed of sound.
This pair style also computes laminar viscosity (Morris1997).
The following coefficients must be defined for each pair of atoms
types via the "pair_coeff"_pair_coeff.html command as in the examples
above.
rho0 reference density (mass/volume units)
c0 reference soundspeed (distance/time units)
nu dynamic viscosity (mass*distance/time units)
h kernel function cutoff (distance units)
:line
[Mixing, shift, table, tail correction, restart, rRESPA info]:
These pair styles do not support mixing. Thus, coefficients for all
I,J pairs must be specified explicitly. Additionally, because rho0 and c0 are defined on a per-type basis,
an I,I as well a J,J coefficients line must be specified prior to specifying an I,J interaction.
This style does not write information to "binary restart
files"_restart.html, so pair_style and pair_coeff commands do need
to be specified in an input script that reads a restart file.
These styles can only be used via the {pair} keyword of the "run_style
respa"_run_style.html command. They do not support the {inner},
{middle}, {outer} keywords.
[Restrictions:]
The {pair sph/taitwater/morris} style is part of the "USER-SPH" package. It is
only enabled if LAMMPS was built with that package (which it is by
default). See the "Making LAMMPS"_Section_start.html#2_3 section for
more info.
[Related commands:]
"pair_coeff"_pair_coeff.html, pair_sph/rhosum
[Default:] none
[(Morris1997)] J.P. Morris, P.J. Fox, and Y. Zhu. Modeling Low Reynolds Number Incompressible Flows Using SPH. Journal of Computational Physics, 136:214–226, 1997.