Page MenuHomec4science
Files F76923141

pair_dsmc.html
No OneTemporary
Actions

File Metadata

Created
Sun, Aug 11, 05:28

pair_dsmc.html
View Options

<!DOCTYPE html>
<!--[if IE 8]><html class="no-js lt-ie9" lang="en" > <![endif]-->
<!--[if gt IE 8]><!--> <html class="no-js" lang="en" > <!--<![endif]-->
<head>
<meta charset="utf-8">
<meta name="viewport" content="width=device-width, initial-scale=1.0">
<title>pair_style dsmc command &mdash; LAMMPS documentation</title>
<link rel="stylesheet" href="_static/css/theme.css" type="text/css" />
<link rel="stylesheet" href="_static/sphinxcontrib-images/LightBox2/lightbox2/css/lightbox.css" type="text/css" />
<link rel="top" title="LAMMPS documentation" href="index.html"/>
<script src="_static/js/modernizr.min.js"></script>
</head>
<body class="wy-body-for-nav" role="document">
<div class="wy-grid-for-nav">
<nav data-toggle="wy-nav-shift" class="wy-nav-side">
<div class="wy-side-nav-search">
<a href="Manual.html" class="icon icon-home"> LAMMPS
</a>
<div role="search">
<form id="rtd-search-form" class="wy-form" action="search.html" method="get">
<input type="text" name="q" placeholder="Search docs" />
<input type="hidden" name="check_keywords" value="yes" />
<input type="hidden" name="area" value="default" />
</form>
</div>
</div>
<div class="wy-menu wy-menu-vertical" data-spy="affix" role="navigation" aria-label="main navigation">
<ul>
<li class="toctree-l1"><a class="reference internal" href="Section_intro.html">1. Introduction</a></li>
<li class="toctree-l1"><a class="reference internal" href="Section_start.html">2. Getting Started</a></li>
<li class="toctree-l1"><a class="reference internal" href="Section_commands.html">3. Commands</a></li>
<li class="toctree-l1"><a class="reference internal" href="Section_packages.html">4. Packages</a></li>
<li class="toctree-l1"><a class="reference internal" href="Section_accelerate.html">5. Accelerating LAMMPS performance</a></li>
<li class="toctree-l1"><a class="reference internal" href="Section_howto.html">6. How-to discussions</a></li>
<li class="toctree-l1"><a class="reference internal" href="Section_example.html">7. Example problems</a></li>
<li class="toctree-l1"><a class="reference internal" href="Section_perf.html">8. Performance &amp; scalability</a></li>
<li class="toctree-l1"><a class="reference internal" href="Section_tools.html">9. Additional tools</a></li>
<li class="toctree-l1"><a class="reference internal" href="Section_modify.html">10. Modifying &amp; extending LAMMPS</a></li>
<li class="toctree-l1"><a class="reference internal" href="Section_python.html">11. Python interface to LAMMPS</a></li>
<li class="toctree-l1"><a class="reference internal" href="Section_errors.html">12. Errors</a></li>
<li class="toctree-l1"><a class="reference internal" href="Section_history.html">13. Future and history</a></li>
</ul>
</div>
&nbsp;
</nav>
<section data-toggle="wy-nav-shift" class="wy-nav-content-wrap">
<nav class="wy-nav-top" role="navigation" aria-label="top navigation">
<i data-toggle="wy-nav-top" class="fa fa-bars"></i>
<a href="Manual.html">LAMMPS</a>
</nav>
<div class="wy-nav-content">
<div class="rst-content">
<div role="navigation" aria-label="breadcrumbs navigation">
<ul class="wy-breadcrumbs">
<li><a href="Manual.html">Docs</a> &raquo;</li>
<li>pair_style dsmc command</li>
<li class="wy-breadcrumbs-aside">
<a href="http://lammps.sandia.gov">Website</a>
<a href="Section_commands.html#comm">Commands</a>
</li>
</ul>
<hr/>
</div>
<div role="main" class="document" itemscope="itemscope" itemtype="http://schema.org/Article">
<div itemprop="articleBody">
<div class="section" id="pair-style-dsmc-command">
<span id="index-0"></span><h1>pair_style dsmc command</h1>
<div class="section" id="syntax">
<h2>Syntax</h2>
<div class="highlight-default"><div class="highlight"><pre><span></span><span class="n">pair_style</span> <span class="n">dsmc</span> <span class="n">max_cell_size</span> <span class="n">seed</span> <span class="n">weighting</span> <span class="n">Tref</span> <span class="n">Nrecompute</span> <span class="n">Nsample</span>
</pre></div>
</div>
<ul class="simple">
<li>max_cell_size = global maximum cell size for DSMC interactions (distance units)</li>
<li>seed = random # seed (positive integer)</li>
<li>weighting = macroparticle weighting</li>
<li>Tref = reference temperature (temperature units)</li>
<li>Nrecompute = recompute v*sigma_max every this many timesteps (timesteps)</li>
<li>Nsample = sample this many times in recomputing v*sigma_max</li>
</ul>
</div>
<div class="section" id="examples">
<h2>Examples</h2>
<div class="highlight-default"><div class="highlight"><pre><span></span><span class="n">pair_style</span> <span class="n">dsmc</span> <span class="mf">2.5</span> <span class="mi">34387</span> <span class="mi">10</span> <span class="mf">1.0</span> <span class="mi">100</span> <span class="mi">20</span>
<span class="n">pair_coeff</span> <span class="o">*</span> <span class="o">*</span> <span class="mf">1.0</span>
<span class="n">pair_coeff</span> <span class="mi">1</span> <span class="mi">1</span> <span class="mf">1.0</span>
</pre></div>
</div>
</div>
<div class="section" id="description">
<h2>Description</h2>
<p>Style <em>dsmc</em> computes collisions between pairs of particles for a
direct simulation Monte Carlo (DSMC) model following the exposition in
<a class="reference internal" href="#bird"><span class="std std-ref">(Bird)</span></a>. Each collision resets the velocities of the two
particles involved. The number of pairwise collisions for each pair
or particle types and the length scale within which they occur are
determined by the parameters of the pair_style and pair_coeff
commands.</p>
<p>Stochastic collisions are performed using the variable hard sphere
(VHS) approach, with the user-defined <em>max_cell_size</em> value used as
the maximum DSMC cell size, and reference cross-sections for
collisions given using the pair_coeff command.</p>
<p>There is no pairwise energy or virial contributions associated with
this pair style.</p>
<p>The following coefficient must be defined for each pair of atoms types
via the <a class="reference internal" href="pair_coeff.html"><span class="doc">pair_coeff</span></a> command as in the examples above,
or in the data file or restart files read by the
<a class="reference internal" href="read_data.html"><span class="doc">read_data</span></a> or <a class="reference internal" href="read_restart.html"><span class="doc">read_restart</span></a>
commands:</p>
<ul class="simple">
<li>sigma (area units, i.e. distance-squared)</li>
</ul>
<p>The global DSMC <em>max_cell_size</em> determines the maximum cell length
used in the DSMC calculation. A structured mesh is overlayed on the
simulation box such that an integer number of cells are created in
each direction for each processor&#8217;s sub-domain. Cell lengths are
adjusted up to the user-specified maximum cell size.</p>
<hr class="docutils" />
<p>To perform a DSMC simulation with LAMMPS, several additional options
should be set in your input script, though LAMMPS does not check for
these settings.</p>
<p>Since this pair style does not compute particle forces, you should use
the &#8220;fix nve/noforce&#8221; time integration fix for the DSMC particles,
e.g.</p>
<div class="highlight-default"><div class="highlight"><pre><span></span><span class="n">fix</span> <span class="mi">1</span> <span class="nb">all</span> <span class="n">nve</span><span class="o">/</span><span class="n">noforce</span>
</pre></div>
</div>
<p>This pair style assumes that all particles will communicated to
neighboring processors every timestep as they move. This makes it
possible to perform all collisions between pairs of particles that are
on the same processor. To ensure this occurs, you should use
these commands:</p>
<div class="highlight-default"><div class="highlight"><pre><span></span><span class="n">neighbor</span> <span class="mf">0.0</span> <span class="nb">bin</span>
<span class="n">neigh_modify</span> <span class="n">every</span> <span class="mi">1</span> <span class="n">delay</span> <span class="mi">0</span> <span class="n">check</span> <span class="n">no</span>
<span class="n">atom_modify</span> <span class="n">sort</span> <span class="mi">0</span> <span class="mf">0.0</span>
<span class="n">communicate</span> <span class="n">single</span> <span class="n">cutoff</span> <span class="mf">0.0</span>
</pre></div>
</div>
<p>These commands ensure that LAMMPS communicates particles to
neighboring processors every timestep and that no ghost atoms are
created. The output statistics for a simulation run should indicate
there are no ghost particles or neighbors.</p>
<p>In order to get correct DSMC collision statistics, users should
specify a Gaussian velocity distribution when populating the
simulation domain. Note that the default velocity distribution is
uniform, which will not give good DSMC collision rates. Specify
&#8220;dist gaussian&#8221; when using the <a class="reference internal" href="velocity.html"><span class="doc">velocity</span></a> command
as in the following:</p>
<div class="highlight-default"><div class="highlight"><pre><span></span><span class="n">velocity</span> <span class="nb">all</span> <span class="n">create</span> <span class="mf">594.6</span> <span class="mi">87287</span> <span class="n">loop</span> <span class="n">geom</span> <span class="n">dist</span> <span class="n">gaussian</span>
</pre></div>
</div>
<hr class="docutils" />
<p><strong>Mixing, shift, table, tail correction, restart, rRESPA info</strong>:</p>
<p>This pair style does not support mixing. Thus, coefficients for all
I,J pairs must be specified explicitly.</p>
<p>This pair style does not support the <a class="reference internal" href="pair_modify.html"><span class="doc">pair_modify</span></a>
shift option for the energy of the pair interaction.</p>
<p>The <a class="reference internal" href="pair_modify.html"><span class="doc">pair_modify</span></a> table option is not relevant
for this pair style.</p>
<p>This pair style does not support the <a class="reference internal" href="pair_modify.html"><span class="doc">pair_modify</span></a>
tail option for adding long-range tail corrections to energy and
pressure.</p>
<p>This pair style writes its information to <a class="reference internal" href="restart.html"><span class="doc">binary restart files</span></a>, so pair_style and pair_coeff commands do not need
to be specified in an input script that reads a restart file. Note
that the user-specified random number seed is stored in the restart
file, so when a simulation is restarted, each processor will
re-initialize its random number generator the same way it did
initially. This means the random forces will be random, but will not
be the same as they would have been if the original simulation had
continued past the restart time.</p>
<p>This pair style can only be used via the <em>pair</em> keyword of the
<a class="reference internal" href="run_style.html"><span class="doc">run_style respa</span></a> command. It does not support the
<em>inner</em>, <em>middle</em>, <em>outer</em> keywords.</p>
</div>
<hr class="docutils" />
<div class="section" id="restrictions">
<h2>Restrictions</h2>
<p>This style is part of the MC package. It is only enabled if LAMMPS
was built with that package. See the <a class="reference internal" href="Section_start.html#start-3"><span class="std std-ref">Making LAMMPS</span></a> section for more info.</p>
</div>
<div class="section" id="related-commands">
<h2>Related commands</h2>
<p><a class="reference internal" href="pair_coeff.html"><span class="doc">pair_coeff</span></a>, <a class="reference internal" href="fix_nve_noforce.html"><span class="doc">fix nve/noforce</span></a>,
<a class="reference internal" href="neigh_modify.html"><span class="doc">neigh_modify</span></a>, <a class="reference internal" href="neighbor.html"><span class="doc">neighbor</span></a>,
<a class="reference internal" href="comm_modify.html"><span class="doc">comm_modify</span></a></p>
<p><strong>Default:</strong> none</p>
<hr class="docutils" />
<p id="bird"><strong>(Bird)</strong> G. A. Bird, &#8220;Molecular Gas Dynamics and the Direct Simulation
of Gas Flows&#8221; (1994).</p>
</div>
</div>
</div>
</div>
<footer>
<hr/>
<div role="contentinfo">
<p>
&copy; Copyright 2013 Sandia Corporation.
</p>
</div>
Built with <a href="http://sphinx-doc.org/">Sphinx</a> using a <a href="https://github.com/snide/sphinx_rtd_theme">theme</a> provided by <a href="https://readthedocs.org">Read the Docs</a>.
</footer>
</div>
</div>
</section>
</div>
<script type="text/javascript">
var DOCUMENTATION_OPTIONS = {
URL_ROOT:'./',
VERSION:'',
COLLAPSE_INDEX:false,
FILE_SUFFIX:'.html',
HAS_SOURCE: true
};
</script>
<script type="text/javascript" src="_static/jquery.js"></script>
<script type="text/javascript" src="_static/underscore.js"></script>
<script type="text/javascript" src="_static/doctools.js"></script>
<script type="text/javascript" src="https://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML"></script>
<script type="text/javascript" src="_static/sphinxcontrib-images/LightBox2/lightbox2/js/jquery-1.11.0.min.js"></script>
<script type="text/javascript" src="_static/sphinxcontrib-images/LightBox2/lightbox2/js/lightbox.min.js"></script>
<script type="text/javascript" src="_static/sphinxcontrib-images/LightBox2/lightbox2-customize/jquery-noconflict.js"></script>
<script type="text/javascript" src="_static/js/theme.js"></script>
<script type="text/javascript">
jQuery(function () {
SphinxRtdTheme.StickyNav.enable();
});
</script>
</body>
</html>

Event Timeline

c4science · Help