Cell-level canonical sampling by velocity scaling for multiparticle collision dynamics simulations

  • Authors:

  • C. C. Huang;A. Chatterji;G. Sutmann;G. Gompper;R. G. Winkler

  • Affiliations:

  • Institut für Festkörperforschung, Forschungszentrum Jülich, 52425 Jülich, Germany;Institut für Festkörperforschung, Forschungszentrum Jülich, 52425 Jülich, Germany and Indian Institute of Science Education and Research (IISER), Pune 411021, India;Jülich Supercomputing Centre, Institute for Advanced Simulation, Forschungszentrum Jülich, 52425 Jülich, Germany;Institut für Festkörperforschung, Forschungszentrum Jülich, 52425 Jülich, Germany and Theoretical Soft-Matter and Biophysics, Institute for Advanced Simulation, Forschungszentr ...;Theoretical Soft-Matter and Biophysics, Institute for Advanced Simulation, Forschungszentrum Jülich, 2425 Jülich, Germany

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2010

Quantified Score

Hi-index 31.45

Visualization

Abstract

A local Maxwellian thermostat for the multiparticle collision dynamics algorithm is proposed. The algorithm is based on a scaling of the relative velocities of the fluid particles within a collision cell. The scaling factor is determined from the distribution of the kinetic energy within such a cell. Thereby the algorithm ensures that the distribution of the relative velocities is given by the Maxwell-Boltzmann distribution. The algorithm is particularly useful for non-equilibrium systems, where temperature has to be controlled locally. We perform various non-equilibrium simulations for fluids in shear and pressure-driven flow, which confirm the validity of the proposed simulation scheme. In addition, we determine the dynamic structure factors for fluids with and without thermostat, which exhibit significant differences due to suppression of the diffusive part of the energy transport of the isothermal system.