BROWSE

Related Researcher

Author's Photo

Baig, Chunggi
Theoretical and Computational study of Polymers & Nanomaterials Lab
Research Interests
  • Multiscale simulation, Polymer rheology, Nonequilibrium molecular dynamics/monte carlo

ITEM VIEW & DOWNLOAD

A proper approach for nonequilibrium molecular dynamics simulations of planar elongational flow

Cited 31 times inthomson ciCited 33 times inthomson ci
Title
A proper approach for nonequilibrium molecular dynamics simulations of planar elongational flow
Author
Baig, ChunggiEdwards, BJKeffer, DJCochran, HD
Issue Date
2005-03
Publisher
AMER INST PHYSICS
Citation
JOURNAL OF CHEMICAL PHYSICS, v.122, no.11, pp.114103
Abstract
We present nonequilibrium molecular dynamics simulations of planar elongational flow (PEF) by an algorithm proposed by Tuckerman et al. [J. Chem. Phys. 106, 5615 (1997)] and theoretically elaborated by Edwards and Dressier [J. Non-Newtonian, Fluid Mech. 96, 163 (2001)], which we shall call the proper-SLLOD algorithm, or p-SLLOD for short. [For background on names of algorithms see W. G. Hoover, D. J. Evans, R. B. Hickman, A. J. C. Ladd, W. T. Ashurst, and B. Moran, Phys. Rev. A 22, 1690 (1980) and D. J. Evans and G. P. Morriss, Phys. Rev. A 30, 1528 (1984).] We show that there are two sources for the exponential growth in PEF of the total linear momentum of the system in the contracting direction, which has been previously observed using the so-called SLLOD algorithm. The first comes from the SLLOD algorithm itself, and the second from the implementation of the Kraynik and Reinelt [Int. J. Multiphase Flow 18, 1045 (1992)] boundary conditions. Using the p-SLLOD algorithm (to eliminate the first source) implemented with our simulation strategy (to eliminate the second) in PEF simulations, we no longer observe the exponential growth. By analyzing the equations of motion, we also demonstrate that both the SLLOD and the DOLLS algorithms are intrinsically unsuitable for representing a nonequilibrium system with elongational flow. However, the p-SLLOD algorithm has a rigorously canonical structure in laboratory phase space, and thus can represent a nonequilibrium system not only for elongational flow but also for a general flow.
URI
https://scholarworks.unist.ac.kr/handle/201301/6964
URL
https://aip.scitation.org/doi/10.1063/1.1819869
DOI
10.1063/1.1819869
ISSN
0021-9606
Appears in Collections:
ECHE_Journal Papers
Files in This Item:
There are no files associated with this item.

find_unist can give you direct access to the published full text of this article. (UNISTARs only)

Show full item record

qrcode

  • mendeley

    citeulike

Items in DSpace are protected by copyright, with all rights reserved, unless otherwise indicated.

MENU