Modelling Capillary Rise with a Slip Boundary Condition: Well-posedness and Long-time Dynamics of Solutions to Washburn's Equation
Isidora Rapajić, Srboljub Simić, Endre Süli
TL;DR
The paper derives a Washburn-type capillary-rise equation that incorporates a slip boundary condition via a slip parameter $\beta$, starting from first principles and a cylindrical geometry. Through nondimensionalization, it obtains the dimensionless model $\omega(HH')' + \beta HH' + H = 1$ with equilibrium height $H_e=1$, and provides a rigorous analysis establishing global existence, uniqueness, and continuous dependence of solutions for initial heights $H(0)=\alpha$ with $\alpha\in[0,3/2]$. It further analyzes long-time behavior by recasting the problem as a 2D autonomous system, identifying a critical threshold $\omega_*=\beta^2/4$ that distinguishes monotone and oscillatory approaches, and proving convergence to the equilibrium via a Lyapunov function and LaSalle’s invariance principle. The results show that while the slip parameter shapes transient dynamics, the unique equilibrium height and global convergence are preserved for all positive $\beta$ and $\omega$, providing a solid mathematical foundation for capillary rise with slip and guiding future numerical and physical investigations.
Abstract
The aim of this paper is to extend Washburn's capillary rise equation by incorporating a slip condition at the pipe wall. The governing equation is derived using fundamental principles from continuum mechanics. A new scaling is introduced, allowing for a systematic analysis of different flow regimes. We prove the global-in-time existence and uniqueness of a bounded positive solution to Washburn's equation that includes the slip parameter, as well as the continuous dependence of the solution in the maximum norm on the initial data. Thus, the initial-value problem for Washburn's equation is shown to be well-posed in the sense of Hadamard. Additionally, we show that the unique equilibrium solution may be reached either monotonically or in an oscillatory fashion, similarly to the no-slip case. Finally, we determine the basin of attraction for the system, ensuring that the equilibrium state will be reached from the initial data we impose. These results hold for any positive value of the nondimensional slip parameter in the model, and for all values of the ratio $h_0/h_e$ in the range $[0,3/2]$, where $h_0$ is the initial height of the fluid column and $h_e$ is its equilibrium height.
