The influence of boundary conditions and interfacial slip on the time taken to achieve a nonequilibrium steady-state for boundary-driven flows

Carmelo Riccardo Civello, Luca Maffioli, Joseph Johnson, Edward R. Smith, James P. Ewen, Peter J. Daivis, Daniele Dini, B. D. Todd

Research output: Contribution to journalArticlepeer-review

Abstract

In this study, we investigate the equilibration time to attain steady-state for a system of liquid molecules under boundary-driven planar Couette flow via nonequilibrium molecular dynamics (NEMD) simulation. In particular, we examine the equilibration time for the two common types of boundary-driven flow: one in which both walls slide with equal and opposite velocity (±û/2), and the other in which one wall is fixed and the other moves with twice the velocity (û). Both flows give identical steady-state strain rates and, hence, flow properties, but the transient behavior is completely different. We find that in the case of no-slip boundary conditions, the equilibration times for the counter-sliding walls flow are exactly four times faster than those of the single-sliding wall system, and this is independent of the atomistic nature of the fluid, i.e., it is an entirely hydrodynamic feature. We also find that systems that exhibit slip have longer equilibration times in general, and the ratio of equilibration times for the two types of boundary-driven flow is even more pronounced. We analyze the problem by decomposing a generic planar Couette flow into a linear sum of purely symmetric and antisymmetric flows. We find that the no-slip equilibration time is dominated by the slowest decaying eigenvalue of the solution to the Navier–Stokes equation. In the case of slip, the longest relaxation time is now dominated by the transient slip velocity response, which is longer than the no-slip response time. In the case of a high-slip system of water confined to graphene channels, the enhancement is over two orders of magnitude. We propose a simple universal relation that predicts the enhanced equilibration time, which agrees well with our NEMD results for simple Lennard-Jones fluids and the water–graphene system. The implications of this significant speed-up in attaining steady-state, which is especially pronounced in the presence of slip, are discussed in general.

Original languageEnglish
Article number122124
Number of pages14
JournalPhysics of Fluids
Volume37
Issue number12
Early online date12 Dec 2025
DOIs
Publication statusPublished - 12 Dec 2025

Bibliographical note

Publisher Copyright:
© 2025 Author(s).

Data Availability Statement

The data that support the findings of this study are available within the article and its supplementary material.

Acknowledgements

The authors acknowledge the Swinburne OzSTAR Supercomputing facility and the Imperial College London Research Computing Service (DOI:10.14469/hpc/223) for providing computational resources for
this work.

Funding

The authors thank the Australian Research Council for a grant obtained through the Discovery Projects Scheme (Grant No. DP200100422) and the Royal Society for support via International Exchanges (Grant No. IES/R3/170/233). J.P.E. was supported by the Royal Academy of Engineering (RAEng) through their Research Fellowships scheme. D.D. was supported through a Shell/RAEng Research Chair in Complex Engineering Interfaces.

Fingerprint

Dive into the research topics of 'The influence of boundary conditions and interfacial slip on the time taken to achieve a nonequilibrium steady-state for boundary-driven flows'. Together they form a unique fingerprint.

Cite this