Periodic Boundary Conditions for Long-time Nonequilibrium Molecular Dynamics Simulations of Incompressible Flows

Abstract

This work presents a generalization of the Kraynik-Reinelt (KR) boundary conditions for nonequilibrium molecular dynamics simulations. In the simulation of steady, homogeneous flows with periodic boundary conditions, the simulation box moves with the flow, and it is possible for particle replicas to become arbitrarily close, causing a breakdown in the simulation. The KR boundary conditions avoid this problem for planar elongational flow and general planar mixed flow [J. Chem. Phys 133, 14116 (2010)] through careful choice of the initial simulation box and by periodically remapping the simulation box in a way that conserves replica locations. In this work, the ideas are extended to a large class of three dimensional flows by using multiple remappings for the simulation box. The simulation box geometry is no longer time-periodic (which was shown to be impossible for uniaxial and biaxial stretching flows in the original work by Kraynik and Reinelt [Int. J. Multiphase Flow 18, 1045 (1992)]). The presented algorithm applies to all flows with nondefective flow matrices, and in particular, to uniaxial and biaxial flows.

11 pages, 3 figures

References in corpus (1)

• Periodic boundary conditions for the simulation of uniaxial extensional flow

Cited by in corpus (18)

• Relating Chain Conformations to Extensional Stress In Entangled Polymer Melts
• Topological Linking Drives Anomalous Thickening of Ring Polymers In Weak Extensional Flows
• Universal dynamics of dilute and semidilute solutions of flexible linear polymers
• Composite Entanglement Topology and Extensional Rheology of Symmetric Ring-Linear Polymer Blends
• A theoretical framework for steady-state rheometry in generic flow conditions
• Viscosity Overshoot in Biaxial Elongational Flow: Coarse-Grained Molecular Dynamics Simulation of Ring-Linear Polymer Mixtures
• Elongational viscosity of weakly entangled polymer melt via coarse-grained molecular dynamics simulation
• Simulation of semidilute polymer solutions in planar extensional via conformationally averaged Brownian noise
• Shear Is Not Always Simple: Rate-Dependent Effects of Flow Type on Granular Rheology
• Topological transition in multicyclic chains with structural symmetry inducing stress-overshoot phenomena in multicyclic/linear blends under biaxial elongational flow
• Universal behavior in fragmenting brittle, isotropic solids across material properties
• Coordinate transformation methodology for simulating quasi-static elastoplastic solids
• Comment on: "Relating Chain Conformations to Extensional Stress in Entangled Polymer Melts"
• On the importance of numerical integration details for homogeneous flow simulation
• Cell List Algorithms for Nonequilibrium Molecular Dynamics
• A multi-scale method for complex flows of non-Newtonian fluids
• Identifying the Threshold Chain Length for Stress Overshoot in Ring-Linear Polymer Blends under Uniaxial Elongation: The Role of Multiple Threading
• Convergence of Nonequilibrium Langevin Dynamics for Planar Flows