Poiseuille flow of liquid methane in nanoscopic graphite channels by molecular dynamics simulation

Martin Horsch,Jadran Vrabec,Martin Bernreuther,Hans Hasse
DOI: https://doi.org/10.48550/arXiv.0911.5481
2009-11-29
Abstract:MD simulations of methane confined between graphite walls with up to 4,800,000 interaction sites, i.e., carbon atoms and methane molecules, are conducted, where the channel width is varied to include both the boundary-dominated regime and the transition to the continuum regime. This proves that MD can be used today to cover the entire range of characteristic lengths for which continuum methods fail.
Mesoscale and Nanoscale Physics,Soft Condensed Matter
What problem does this paper attempt to address?
This paper aims to solve the problem of how to study the flow behavior of liquid methane in graphite nano - channels when the continuum methods (such as Navier - Stokes equations) become invalid at the nanoscale. Specifically, the author uses molecular dynamics (MD) simulations to explore the Poiseuille flow characteristics of liquid methane at the nanoscale. ### Main problems 1. **Limitations of continuum methods**: At the nanoscale, traditional continuum methods (such as Navier - Stokes equations) are no longer applicable, so research needs to be carried out at the molecular level. This requires the use of molecular dynamics simulations to understand the transport processes at the nanoscale. 2. **Flow behavior in nanostructures**: Graphite and carbon nanotubes have become the focus of research because of their anisotropy. This paper pays special attention to the flow behavior of liquid methane between graphite walls, especially Poiseuille flow, that is, laminar flow driven by external forces. 3. **Simulation methods and scalability**: Due to the limitation of computing resources, early molecular dynamics simulations could only handle small systems containing several thousand particles. With the improvement of computing power, larger - scale systems can be simulated now, but the simulation methods need to be optimized to achieve efficient parallel processing. 4. **Boundary slip phenomenon**: At the nanoscale, boundary slip has a significant impact on the flow behavior. Especially when the channel width is less than 2 nanometers, the slip velocity may reach more than 99% of the characteristic velocity. This paper studies the boundary slip phenomenon under different channel widths and its impact on the flow characteristics. ### Specific objectives - **Verify Darcy's law**: Study whether the flow of liquid methane conforms to Darcy's law under different channel widths. - **Explore the transition region**: Determine the transition range from the boundary - dominated region to the continuum region, that is, the flow characteristics when the channel width is between 20 and 50 nanometers. - **Optimize the simulation program**: Evaluate and optimize the LS1/Mardyn program for large - scale parallel molecular dynamics simulations to ensure its efficiency and scalability when dealing with complex systems. Through these studies, the author hopes to demonstrate the application potential of molecular dynamics simulations in the study of nanoscale flow and provide valuable references for future research.