OCTP

A Tool for On-the-Fly Calculation of Transport Properties of Fluids with the Order- n Algorithm in LAMMPS

Journal article (2019)

Authors

S.H. Jamali Engineering Thermodynamics - Mechanical, Maritime and Materials Engineering
Ludger Wolff Rheinisch-Westfälische Technische Hochschule
T. Becker Engineering Thermodynamics - Mechanical, Maritime and Materials Engineering
M. de Groen Engineering Thermodynamics - Mechanical, Maritime and Materials Engineering
M. Ramdin Engineering Thermodynamics - Mechanical, Maritime and Materials Engineering
Remco Hartkamp Complex Fluid Processing - Mechanical, Maritime and Materials Engineering
André Bardow Rheinisch-Westfälische Technische Hochschule
T.J.H. Vlugt Engineering Thermodynamics - Mechanical, Maritime and Materials Engineering
O. Moultos Engineering Thermodynamics - Mechanical, Maritime and Materials Engineering
Research Group
Complex Fluid Processing (Mechanical, Maritime and Materials Engineering) (TU Delft)
Copyright: © 2019 S.H. Jamali, Ludger Wolff, T. Becker, M. de Groen, M. Ramdin, Remco Hartkamp, André Bardow, T.J.H. Vlugt, O. Moultos
DOI: https://doi.org/10.1021/acs.jcim.8b00939
More Info
expand_more

Published Date

2019

Language

English

Faculty

Mechanical, Maritime and Materials Engineering

Department

Process and Energy

Research Group

Complex Fluid Processing

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.