Contribution of wedge and bulk viscous forces in droplets moving on inclined surfaces

Francisco Bodziony,Xiaomei Li,Mariana Yin,Rüdiger Berger,Hans-Jürgen Butt,Holger Marschall
DOI: https://doi.org/10.1007/s00162-024-00692-3
IF: 2.892
2024-04-19
Theoretical and Computational Fluid Dynamics
Abstract:Employing direct numerical simulations, we investigate water and water-glycerol (85 wt%) droplets ( 25 μL) moving on smooth surfaces, with contact angles of around 90 , at varying inclinations. Our focus is on elucidating the relative contribution of local viscous forces in the wedge and bulk regions in droplets to the total viscous force. We observe that, for fast-moving droplets, both regions contribute comparably, while the contribution of the wedge region dominates in slow-moving cases. Comparisons with existing estimates reveal the inadequacy of previous predictions in capturing the contributions of wedge and bulk viscous forces in fast-moving droplets. Furthermore, we demonstrate that droplets with identical velocities can exhibit disparate viscous forces due to variations in internal fluid dynamics.
mechanics,physics, fluids & plasmas
What problem does this paper attempt to address?
This paper aims to explore the relative contributions of the viscous forces in the wedge region and the main body region to the total viscous force when a droplet moves on an inclined surface. Specifically, through direct numerical simulation (DNS), the researchers studied the movement of water and water - glycerol mixture droplets (about 25 microliters) on a smooth surface with a contact angle of about 90 degrees at different inclination angles. They focused on how the local viscous forces in the wedge region and the main body region work together in a fast - moving droplet, and how the contribution of the wedge region dominates in a slow - moving droplet. In addition, the study also found that even droplets with the same speed may exhibit different viscous forces due to different internal fluid dynamics. These findings challenge the deficiencies in previous predictions of the contributions of the viscous forces in the wedge region and the main body region in fast - moving droplets. ### Problems the paper attempts to solve 1. **Quantify the contributions of the viscous forces in the wedge region and the main body region**: - Through numerical simulation, the researchers analyzed in detail the relative contributions of the viscous forces in the wedge region and the main body region to the total viscous force when a droplet moves on an inclined surface. In particular, how the contributions of these two forces change in a fast - moving droplet. 2. **Verify the deficiencies of existing theories**: - Through the simulation of droplets at different speeds, the researchers found that the existing theoretical predictions are insufficient in capturing the contributions of the viscous forces in the wedge region and the main body region in fast - moving droplets. This indicates that more accurate methods are needed to describe the actions of these forces. 3. **Explore the influence of internal fluid dynamics**: - The study also revealed that even droplets with the same speed may exhibit different viscous forces due to different internal fluid dynamics. This finding emphasizes the importance of internal fluid dynamics in droplet movement. ### Formulas and methods - **Navier - Stokes equation**: \[ \nabla \cdot \mathbf{u} = 0 \] \[ \frac{\partial (\rho \mathbf{u})}{\partial t} + \nabla \cdot (\rho \mathbf{u} \otimes \mathbf{u}) = - \nabla \tilde{p} + \nabla \cdot [\mu (\nabla \mathbf{u} + (\nabla \mathbf{u})^T)] + \rho \mathbf{g} - \nabla \cdot (\mathbf{u} \otimes \mathbf{J}) + \mathbf{f}_s \] where \(\mathbf{u}\) is the flow velocity, \(\mathbf{g}\) is the gravity vector, \(\tilde{p}\) is the modified pressure, \(\mu\) is the dynamic viscosity, \(\mathbf{J}\) is the diffusive phase - field mass flux, and \(\mathbf{f}_s\) is the source term. - **Cahn - Hilliard equation**: \[ \frac{\partial C_p}{\partial t} + \nabla \cdot (C_p \mathbf{u}) = - \nabla \cdot \mathbf{J}_p \] \[ \mathbf{J}_p = \sum_{q = 1}^N - M_{p,q} \nabla \Phi_q \] where \(C_p\) is the phase - field order parameter, \(\mathbf{J}_p\) is the phase - field flux, \(\Phi_q\) is the chemical potential, and \(M_{p,q}\) is the inter - phase mobility. - **Chemical potential**: \[ \Phi_p = \lambda_{p,q} \left( \frac{1}{\epsilon^2} \Psi'(C_p) + \nabla^2 C_p \right) \] where \(\lambda_{p,q}\) is the mixing energy, \(\epsilon\) is the capillary width, and \(\Psi(C_p)\)