Reinitialization of the Level-Set Function in 3d Simulation of Moving Contact Lines

2016 ◽  
Vol 20 (5) ◽  
pp. 1163-1182 ◽  
Author(s):  
Shixin Xu ◽  
Weiqing Ren

AbstractThe level set method is one of the most successful methods for the simulation of multi-phase flows. To keep the level set function close the signed distance function, the level set function is constantly reinitialized by solving a Hamilton-Jacobi type of equation during the simulation. When the fluid interface intersects with a solid wall, a moving contact line forms and the reinitialization of the level set function requires a boundary condition in certain regions on the wall. In this work, we propose to use the dynamic contact angle, which is extended from the contact line, as the boundary condition for the reinitialization of the level set function. The reinitialization equation and the equation for the normal extension of the dynamic contact angle form a coupled system and are solved simultaneously. The extension equation is solved on the wall and it provides the boundary condition for the reinitialization equation; the level set function provides the directions along which the contact angle is extended from the contact line. The coupled system is solved using the 3rd order TVD Runge-Kutta method and the Godunov scheme. The Godunov scheme automatically identifies the regions where the angle condition needs to be imposed. The numerical method is illustrated by examples in three dimensions.

Author(s):  
Pablo Gómez ◽  
Adolfo Esteban ◽  
Claudio Zanzi ◽  
Joaquín López ◽  
Julio Hernández

Abstract We present a method based on a level set formulation to reproduce the behavior of the contact line on solid walls in the simulation of 3D unsteady interfacial flows characterized by large density ratios. The level set method poses a particular difficulty, related to the reinitialization procedure, when used in the simulation of interfacial flows in which the interface intersects a solid wall, due to the appearance of a blind zone where standard reinitialization procedures produce inconsistent results. The proposed method overcomes this difficulty by introducing a boundary condition for the level set function on the solid surface based on the normal extension of the contact angle from the interface along the solid wall. In order to reproduce the dynamics of the contact line we use a simplified model that imposes a boundary condition on the interface curvature based on the static contact angle, and define a thin slip zone at the solid wall around the contact line. To assess the accuracy and robustness of the proposed method, we conducted several preliminary numerical tests in three dimensions, whose results are compared with analytical solutions and other results available in the literature.


2020 ◽  
Vol 229 (10) ◽  
pp. 1945-1977 ◽  
Author(s):  
Yulii D. Shikhmurzaev

Abstract After a brief overview of the ‘moving contact-line problem’ as it emerged and evolved as a research topic, a ‘litmus test’ allowing one to assess adequacy of the mathematical models proposed as solutions to the problem is described. Its essence is in comparing the contact angle, an element inherent in every model, with what follows from a qualitative analysis of some simple flows. It is shown that, contrary to a widely held view, the dynamic contact angle is not a function of the contact-line speed as for different spontaneous spreading flows one has different paths in the contact angle-versus-speed plane. In particular, the dynamic contact angle can decrease as the contact-line speed increases. This completely undermines the search for the ‘right’ velocity-dependence of the dynamic contact angle, actual or apparent, as a direction of research. With a reference to an earlier publication, it is shown that, to date, the only mathematical model passing the ‘litmus test’ is the model of dynamic wetting as an interface formation process. The model, which was originated back in 1993, inscribes dynamic wetting into the general physical context as a particular case in a wide class of flows, which also includes coalescence, capillary breakup, free-surface cusping and some other flows, all sharing the same underlying physics. New challenges in the field of dynamic wetting are discussed.


2019 ◽  
Vol 873 ◽  
pp. 110-150
Author(s):  
Hsien-Hung Wei ◽  
Heng-Kwong Tsao ◽  
Kang-Ching Chu

In the context of dynamic wetting, wall slip is often treated as a microscopic effect for removing viscous stress singularity at a moving contact line. In most drop spreading experiments, however, a considerable amount of slip may occur due to the use of polymer liquids such as silicone oils, which may cause significant deviations from the classical Tanner–de Gennes theory. Here we show that many classical results for complete wetting fluids may no longer hold due to wall slip, depending crucially on the extent of de Gennes’s slipping ‘foot’ to the relevant length scales at both the macroscopic and microscopic levels. At the macroscopic level, we find that for given liquid height $h$ and slip length $\unicode[STIX]{x1D706}$, the apparent dynamic contact angle $\unicode[STIX]{x1D703}_{d}$ can change from Tanner’s law $\unicode[STIX]{x1D703}_{d}\sim Ca^{1/3}$ for $h\gg \unicode[STIX]{x1D706}$ to the strong-slip law $\unicode[STIX]{x1D703}_{d}\sim Ca^{1/2}\,(L/\unicode[STIX]{x1D706})^{1/2}$ for $h\ll \unicode[STIX]{x1D706}$, where $Ca$ is the capillary number and $L$ is the macroscopic length scale. Such a no-slip-to-slip transition occurs at the critical capillary number $Ca^{\ast }\sim (\unicode[STIX]{x1D706}/L)^{3}$, accompanied by the switch of the ‘foot’ of size $\ell _{F}\sim \unicode[STIX]{x1D706}Ca^{-1/3}$ from the inner scale to the outer scale with respect to $L$. A more generalized dynamic contact angle relationship is also derived, capable of unifying Tanner’s law and the strong-slip law under $\unicode[STIX]{x1D706}\ll L/\unicode[STIX]{x1D703}_{d}$. We not only confirm the two distinct wetting laws using many-body dissipative particle dynamics simulations, but also provide a rational account for anomalous departures from Tanner’s law seen in experiments (Chen, J. Colloid Interface Sci., vol. 122, 1988, pp. 60–72; Albrecht et al., Phys. Rev. Lett., vol. 68, 1992, pp. 3192–3195). We also show that even for a common spreading drop with small macroscopic slip, slip effects can still be microscopically strong enough to change the microstructure of the contact line. The structure is identified to consist of a strongly slipping precursor film of length $\ell \sim (a\unicode[STIX]{x1D706})^{1/2}Ca^{-1/2}$ followed by a mesoscopic ‘foot’ of width $\ell _{F}\sim \unicode[STIX]{x1D706}Ca^{-1/3}$ ahead of the macroscopic wedge, where $a$ is the molecular length. It thus turns out that it is the ‘foot’, rather than the film, contributing to the microscopic length in Tanner’s law, in accordance with the experimental data reported by Kavehpour et al. (Phys. Rev. Lett., vol. 91, 2003, 196104) and Ueno et al. (Trans. ASME J. Heat Transfer, vol. 134, 2012, 051008). The advancement of the microscopic contact line is still led by the film whose length can grow as the $1/3$ power of time due to $\ell$, as supported by the experiments of Ueno et al. and Mate (Langmuir, vol. 28, 2012, pp. 16821–16827). The present work demonstrates that the behaviour of a moving contact line can be strongly influenced by wall slip. Such slip-mediated dynamic wetting might also provide an alternative means for probing slippery surfaces.


1997 ◽  
Vol 334 ◽  
pp. 211-249 ◽  
Author(s):  
YULII D. SHIKHMURZAEV

A general mathematical model which describes the motion of an interface between immiscible viscous fluids along a smooth homogeneous solid surface is examined in the case of small capillary and Reynolds numbers. The model stems from a conclusion that the Young equation, σ1 cos θ = σ2 − σ3, which expresses the balance of tangential projection of the forces acting on the three-phase contact line in terms of the surface tensions σi and the contact angle θ, together with the well-established experimental fact that the dynamic contact angle deviates from the static one, imply that the surface tensions of contacting interfaces in the immediate vicinity of the contact line deviate from their equilibrium values when the contact line is moving. The same conclusion also follows from the experimentally observed kinematics of the flow, which indicates that liquid particles belonging to interfaces traverse the three-phase interaction zone (i.e. the ‘contact line’) in a finite time and become elements of another interface – hence their surface properties have to relax to new equilibrium values giving rise to the surface tension gradients in the neighbourhood of the moving contact line. The kinematic picture of the flow also suggests that the contact-line motion is only a particular case of a more general phenomenon – the process of interface formation or disappearance – and the corresponding mathematical model should be derived from first principles for this general process and then applied to wetting as well as to other relevant flows. In the present paper, the simplest theory which uses this approach is formulated and applied to the moving contact-line problem. The model describes the true kinematics of the flow so that it allows for the ‘splitting’ of the free surface at the contact line, the appearance of the surface tension gradients near the contact line and their influence upon the contact angle and the flow field. An analytical expression for the dependence of the dynamic contact angle on the contact-line speed and parameters characterizing properties of contacting media is derived and examined. The role of a ‘thin’ microscopic residual film formed by adsorbed molecules of the receding fluid is considered. The flow field in the vicinity of the contact line is analysed. The results are compared with experimental data obtained for different fluid/liquid/solid systems.


2014 ◽  
Vol 746 ◽  
Author(s):  
E. Kirkinis ◽  
S. H. Davis

AbstractA recent hydrodynamic theory of liquid slippage on a solid substrate (Kirkinis & Davis, Phys. Rev. Lett., vol. 110, 2013, 234503) gives rise to a sequence of eddies (Moffatt vortices) that co-move with a moving contact line (CL) in a liquid wedge. The presence of these vortices is established through secular equations that depend on the dynamic contact angle $\alpha $ and capillary number Ca. The limiting case $\alpha \rightarrow 0$ is associated with the appearance of such vortices in a channel. The vortices are generated by the relative motion of the interfaces, which in turn is due to the motion of the CL. This effect has yet to be observed in experiment.


2018 ◽  
Vol 841 ◽  
pp. 767-783 ◽  
Author(s):  
Yi Xia ◽  
Paul H. Steen

Contact-line mobility characterizes how fast a liquid can wet or unwet a solid support by relating the contact angle $\unicode[STIX]{x0394}\unicode[STIX]{x1D6FC}$ to the contact-line speed $U_{CL}$. The contact angle changes dynamically with contact-line speeds during rapid movement of liquid across a solid. Speeds beyond the region of stick–slip are the focus of this experimental paper. For these speeds, liquid inertia and surface tension compete while damping is weak. The mobility parameter $M$ is defined empirically as the proportionality, when it exists, between $\unicode[STIX]{x0394}\unicode[STIX]{x1D6FC}$ and $U_{CL}$, $M\unicode[STIX]{x0394}\unicode[STIX]{x1D6FC}=U_{CL}$. We discover that $M$ exists and measure it. The experimental approach is to drive the contact line of a sessile drop by a plane-normal oscillation of the drop’s support. Contact angles, displacements and speeds of the contact line are measured. To unmask the mobility away from stick–slip, the diagram of $\unicode[STIX]{x0394}\unicode[STIX]{x1D6FC}$ against $U_{CL}$, the traditional diagram, is remapped to a new diagram by rescaling with displacement. This new diagram reveals a regime where $\unicode[STIX]{x0394}\unicode[STIX]{x1D6FC}$ is proportional to $U_{CL}$ and the slope yields the mobility $M$. The experimental approach reported introduces the cyclically dynamic contact angle goniometer. The concept and method of the goniometer are illustrated with data mappings for water on a low-hysteresis non-wetting substrate.


Author(s):  
Ndivhuwo Musehane ◽  
Rhameez Herbst

Computational Fluid Dynamics (CFD) is used to study the spreading process of a water droplet with a radius of 0.00275mm impacting a wax surface at a velocity of 1.18ms−1 . This type of flow is considered to be Multiphase, incompressible, laminar, surface tension dominated and is governed by the Navier stokes and continuity equations. To accurately model the spreading process 3 different contact angle models are investigated, two of which take into account the moving contact line. The governing equations are solved using the open source C++ library OpenFOAM, which uses a Finite Volume Method (FVM) of discretization and a Volume Of Fluid (VOF) interface capturing method. The VOF method is known to produce unphysical velocities when high pressure gradients exist between the two phases, thus a numerical improvement is implemented to reduce the magnitudes of the unphysical velocities. The improvement reduces the magnitudes of the unphysical velocities and as shown in literature their magnitudes increase with an increase in surface tension dominance. The improvement is implemented together with different contact angle models and results obtained show that contact angle models that take into account the moving contact line gives a good correlation of the spreading diameter obtained numerically with the one obtained experimentally.


1994 ◽  
Vol 116 (4) ◽  
pp. 938-945 ◽  
Author(s):  
P. C. Wayner

A new physical model for the spreading dynamics of fluids with an apparent finite contact angle on solid substrates is presented. The model is based on the premise that both interfacial intermolecular forces and temperature control change-of-phase heat transfer and (therefore) motion in the moving contact line region. Classical change-of-phase kinetics and interfacial concepts like the Kelvin–Clapeyron, Young–Dupre, and augmented Young–Laplace equations are used to compare the effects of stress (change in apparent dynamic contact angle) and temperature (superheat). Explicit equations are obtained for the velocity, heat flux, and superheat in the contact line region as a function of the change in the apparent contact angle. Comparisons with experimental data demonstrate that the resulting interfacial model of evaporation/condensation not only describes the “apparently isothermal” contact line movement in these systems at 20°C but also describes the substrate superheat at the critical heat flux.


1997 ◽  
Vol 342 ◽  
pp. 87-96 ◽  
Author(s):  
E. RAMÉ

A good approximation to modelling the shape of a liquid–air meniscus advancing or receding in a capillary tube of radius a can be constructed by balancing the curvature of the interface with the sum of a viscous stress valid near the contact line and a constant static pressure. This model has unique solutions for each value of the boundary condition, i.e. the dynamic contact angle. When the meniscus recedes at very small capillary numbers, the model predicts a critical receding velocity beyond which a liquid layer of the receding fluid (a liquid tail) develops along the solid (see figure 4). The length of the layer increases as the receding speed and the contact angle decrease. This layer regime is characterized by menisci whose macroscopic curvature is >1/a.


Sign in / Sign up

Export Citation Format

Share Document