OCTP
A Tool for On-the-Fly Calculation of Transport Properties of Fluids with the Order- n Algorithm in LAMMPS
Seyed Hossein Jamali (TU Delft - Engineering Thermodynamics)
Ludger Wolff (RWTH Aachen University)
Tim Becker (TU Delft - Engineering Thermodynamics)
M De Groen (TU Delft - Engineering Thermodynamics)
M Ramdin (TU Delft - Engineering Thermodynamics)
Remco Hartkamp (TU Delft - Complex Fluid Processing)
Andre Bardow (RWTH Aachen University)
Thijs JH Vlugt (TU Delft - Engineering Thermodynamics)
Othon Moultos (TU Delft - Engineering Thermodynamics)
More Info
expand_more
Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.
Abstract
We present a new plugin for LAMMPS for on-the-fly computation of transport properties (OCTP) in equilibrium molecular dynamics. OCTP computes the self- and Maxwell-Stefan diffusivities, bulk and shear viscosities, and thermal conductivities of pure fluids and mixtures in a single simulation. OCTP is the first implementation in LAMMPS that uses the Einstein relations combined with the order-n algorithm for the efficient sampling of dynamic variables. OCTP has low computational requirements and is easy to use because it follows the native input file format of LAMMPS. A tool for calculating the radial distribution function (RDF) of the fluid beyond the cutoff radius, while taking into account the system size effects, is also part of the new plugin. The RDFs computed from OCTP are needed to obtain the thermodynamic factor, which relates Maxwell-Stefan and Fick diffusivities. To demonstrate the efficiency of the new plugin, the transport properties of an equimolar mixture of water-methanol were computed at 298 K and 1 bar.
help_outline