TY - JOUR
T1 - Stability and bifurcation of dynamic contact lines in two dimensions
AU - Keeler, Jack S.
AU - Lockerby, Duncan A.
AU - Kumar, Satish
AU - Sprittles, James E.
N1 - Data availability statement: The data that support the findings of this study are openly available in figshare at https://doi.org/10.6084/m9.figshare.17277812.
PY - 2022/7/27
Y1 - 2022/7/27
N2 - The moving contact line between a fluid, liquid and solid is a ubiquitous phenomenon, and determining the maximum speed at which a liquid can wet/dewet a solid is a practically important problem. Using continuum models, previous studies have shown that the maximum speed of wetting/dewetting can be found by calculating steady solutions of the governing equations and locating the critical capillary number, Cacrit , above which no steady-state solution can be found. Below Cacrit , both stable and unstable steady-state solutions exist and if some appropriate measure of these solutions is plotted against Ca , a fold bifurcation appears where the stable and unstable branches meet. Interestingly, the significance of this bifurcation structure to the transient dynamics has yet to be explored. This article develops a computational model and uses ideas from dynamical systems theory to show the profound importance of the unstable solutions on the transient behaviour. By perturbing the stable state by the eigenmodes calculated from a linear stability analysis it is shown that the unstable branch is an ‘edge’ state that is responsible for the eventual dynamical outcomes and that the system can become transient when Ca<Cacrit due to finite-amplitude perturbations. Furthermore, when Ca>Cacrit , we show that the trajectories in phase space closely follow the unstable branch.
AB - The moving contact line between a fluid, liquid and solid is a ubiquitous phenomenon, and determining the maximum speed at which a liquid can wet/dewet a solid is a practically important problem. Using continuum models, previous studies have shown that the maximum speed of wetting/dewetting can be found by calculating steady solutions of the governing equations and locating the critical capillary number, Cacrit , above which no steady-state solution can be found. Below Cacrit , both stable and unstable steady-state solutions exist and if some appropriate measure of these solutions is plotted against Ca , a fold bifurcation appears where the stable and unstable branches meet. Interestingly, the significance of this bifurcation structure to the transient dynamics has yet to be explored. This article develops a computational model and uses ideas from dynamical systems theory to show the profound importance of the unstable solutions on the transient behaviour. By perturbing the stable state by the eigenmodes calculated from a linear stability analysis it is shown that the unstable branch is an ‘edge’ state that is responsible for the eventual dynamical outcomes and that the system can become transient when Ca<Cacrit due to finite-amplitude perturbations. Furthermore, when Ca>Cacrit , we show that the trajectories in phase space closely follow the unstable branch.
U2 - 10.1017/jfm.2022.526
DO - 10.1017/jfm.2022.526
M3 - Article
VL - 945
JO - Journal of Fluid Mechanics
JF - Journal of Fluid Mechanics
SN - 0022-1120
M1 - A34
ER -