Implicit particle wall boundary condition in molecular dynamics

Thin film and nano-tube manufacturing, micro-channel cooling, and many other similar interesting techniques demand the prediction of heat transfer characteristics at the nanometre scale. In this respect, the transport properties at gas—solid and liquid—solid interfaces are very important. The proces...

Full description

Saved in:
Bibliographic Details
Published in:Proceedings of the Institution of Mechanical Engineers. Part C, Journal of mechanical engineering science Journal of mechanical engineering science, 2008-05, Vol.222 (5), p.855-864
Main Authors: Spijker, P, ten Eikelder, H M M, Markvoort, A J, Nedea, S V, Hilbers, P A J
Format: Article
Language:English
Subjects:
Algorithms
Boundary conditions
Channels
Computer simulation
Cooling
Crystal lattices
Demand
Heat flux
Heat transfer
Mathematical analysis
Molecular dynamics
Nanostructure
Nanotubes
Parallel plates
Phase transitions
Simulation
Temperature effects
Transport properties
Velocity
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Thin film and nano-tube manufacturing, micro-channel cooling, and many other similar interesting techniques demand the prediction of heat transfer characteristics at the nanometre scale. In this respect, the transport properties at gas—solid and liquid—solid interfaces are very important. The processes at these interfaces can be studied in detail with molecular dynamics (MD) simulations. However, the computational cost involved in simulating the solid wall currently restrains the size of channels, which can be simulated. Therefore, the solid wall is sometimes replaced by boundary conditions, which often compromise on macroscopic quantities, such as density, temperature, pressure, and heat flux. In the current paper, a new particle wall boundary condition is presented, which is in good agreement with existing boundary conditions, but allows for the pressure calculation. This new boundary condition is based on averaging the contributions of an explicit solid wall and is derived using knowledge on common practices in MD algorithms, such as truncation and shifting. Moreover, it allows for different crystal lattices to be included in the new potential. The applicability of the new method is demonstrated by MD simulations of a gas between two parallel plates at different temperatures and densities. Furthermore, these simulations are compared with explicit wall simulations and existing boundary conditions.
ISSN:0954-4062
2041-2983
DOI:10.1243/09544062JMES713