On dual-grid level-set method for contact line modeling during impact of a droplet on hydrophobic and superhydrophobic surfaces
Graphical abstract
Introduction
Droplet impact dynamics on hydrophobic and superhydrophobic surfaces is both scientifically exciting as well as relevant to many practical applications. The potential applications are in the design and manufacturing of effective spray coating (Bergeron et al., 2000), self-cleaning surfaces (Blossey, 2003), drag reduction for fluid flow (Truesdell et al., 2006), anti-snow adhesion surfaces in aviation applications (Cao et al., 2009) and hot surface cooling (Betz et al., 2010, Patil and Bhardwaj, 2014). The impact of droplet is a highly transient process which depends on forces due to surface tension, inertia, viscosity, gravity and solid surface wettability. The interplay of these forces determines spreading, bouncing, non-bouncing or splashing of the droplet (see review by Yarin, 2006, Marengo et al., 2011). The wetting at the solid–liquid–gas contact line plays an important role during the impact dynamics. The dynamics at the wetting line is mainly governed by the equilibrium, advancing and receding contact angles of the droplet on a surface. The equilibrium contact angle value characterizes the surface property, i.e., hydrophilic or hydrophobic. The equilibrium contact angle varies between 0 and less than 90°, for hydrophilic; 90 and less than 150°, for hydrophobic; and 150 and 180° for superhydrophobic surfaces. The hydrophobic/superhydrophobic surfaces exhibit larger equilibrium contact angle and lower contact angle hysteresis as compared to the hydrophilic surfaces; thus, such surfaces when inclined leads to rolling of the water droplets – taking away the dirt particles and clean up the surfaces (Blossey, 2003).
In the literature, several numerical investigations for the droplet impact on various solid surfaces are found. These simulations were performed on two types of grid: fixed (Renardy et al., 2001, Pasandideh-Fard et al., 2002, Gunjal et al., 2005, Yokoi et al., 2009) and moving (Fukai et al., 1995, Bhardwaj and Attinger, 2008). The former one corresponds to marker and cell (MAC) (Harlow and Shannon, 1967, Unverdi and Tryggvason, 1992), volume of fluid (VOF) (Gunjal et al., 2005, Dupont and Legendre, 2010, Malgarinos et al., 2014), level-set method (LSM) (Ding and Spelt, 2007, Caviezel et al., 2008, Lee and Son, 2011), and combined LS-VOF (Yokoi et al., 2009) methods, and the latter corresponds to the finite element method (FEM) (Fukai et al., 1995, Bhardwaj and Attinger, 2008). In recent times, the Lattice–Boltzmann method (Mukherjee and Abraham, 2007, Tanaka et al., 2011) is also used for the simulation of the droplet impact problem. VOF and LSM based formulations can handle severe interface deformations easily, and have been comprehensively employed for the droplet dynamics simulations.
There are numerous works on the improvement of the biggest disadvantage of LSM – mass-error. One such method involves improvement in the accuracy of capturing interface, by solving level-set equations on a finer grid resolution as compared to the grid used for solving Navier–Stokes (NS) equations for the flow, proposed by Gómez et al. (2005) and called as Dual Resolution LSM (DR-LSM). The LS/advection equations take much less computational time as compared to the nonlinear coupled NS/flow equations, leading to a slight increase in the total computational time for DR-LSM as compared to LSM; however, there is a substantial improvement of the accuracy/mass-error of the results.
The initial work on DR-LSM – by Gómez et al. (2005) for structured/Cartesian and by Herrmann (2008) for unstructured grid – was on locally-refined adaptive grid; the grid-refinement was restricted to a narrow band near the interface. However, to avoid dynamic adaption of the locally refined grid, a globally-refined DR-LSM was proposed by Gada and Sharma (2011) – on a uniform Cartesian grid. They called it a Dual-grid level-set method (DGLSM), with a built-in grid refinement everywhere in the domain; one additional grid for the level-set function in-between two uniformly spaced grid points for velocity/pressure/temperature is used. The partially refined DGLSM takes slightly more computational time (as compared to completely coarse-grid based LSM), to achieve a computational accuracy slightly less than that obtained by a completely refined-grid based LSM. This was demonstrated – qualitatively as well quantitatively – for various two-phase flow problems in Cartesian coordinates by Gada and Sharma (2011), and recently for two fluid electro-dynamic axisymmetric flow by Lakdawala et al. (2015a). The globally-refined DGLSM is much easier to implement as compared to the locally-refined adaptive-grid DR-LSM. The extension of the globally-refined dual-grid strategy for the volume-of-fluid method – DR-VOFM – was proposed recently by Ding and Yuan (2014). A review on LSM based development, application and analysis for multi-phase flow was reported recently by Sharma (2014).
The challenging task in the various numerical models used for the droplet impact simulations is the implementation/modeling of the moving three phase (solid-liquid-gas) wetting/contact line which characterizes the droplet fate. The difficulty arises due to the presence of hydrodynamic stress singularity at the moving contact line because of no-slip boundary conditions at the solid surface (Huh and Scriven, 1971, Shikhmurzaev, 2006). Thus, at the contact point, a slip BC is required to move the contact line (Shikhmurzaev, 2006). Several numerical studies can be found for the modeling of moving contact line (see review by Bonn et al., 2009 and Sui et al., 2014). The correct implementation of dynamic contact angle models is prominent factor for the accuracy of the results. At the contact line, dynamic contact angle (θd) and contact line velocity (UCL) are interdependent parameters, which are responsible for the droplet spreading as well as receding/rebounding (Ngan and Dussan, 1989). Till date, there is a scarcity of universal contact line boundary condition at the contact line. Therefore, the experimentally measured equilibrium (θeq), advancing (θadv) and receding (θrec) contact angles are used in the contact line model.
In the literature, two types of contact line models are found namely, hydrodynamic (Dussan and Davis, 1974, Cox, 1986) and molecular-kinetic (Blake and De Coninck, 2002) models. In the former (used in the present paper), the growth of an apparent dynamic contact angle is accounted by the hydrodynamic deformation of the interface at the contact line (Cox, 1986). Such approach can be found in Fukai et al., 1995, Bussmann et al., 2000, Spelt, 2005, Liu et al., 2005, Ding and Spelt, 2007, Shin and Juric, 2009 and Lee and Son (2011). They treated contact line boundary condition mainly by prescribing the dynamic contact angle as a function of θadv and θrec based on the sign of the contact line velocity; for example, one of the first such model – by Fukai et al. (1995) – is for UCL ≥ 0 and for UCL < 0. The droplet impact simulations by Ngan and Dussan, 1989, Spelt, 2005, Mukherjee and Abraham, 2007, Shikhmurzaev, 2006, Yokoi et al., 2009 and Griebel and Klitz (2013) calculated the dynamic contact angle as a function of several parameters such as equilibrium contact angle, dynamic advancing and receding contact angle, Capillary number (Ca), Weber number (We) or substrate material related constants (kadv and krec). Earlier, Pasandideh-Fard et al. (1996) implemented contact angle effects by measuring its temporal variation from the experiments and using these data as boundary conditions for the numerical simulations. Furthermore, Renardy et al. (2001) and Caviezel et al. (2008) used a constant equilibrium contact angle value throughout the droplet impact simulations and also found good agreement with measurements. In the molecular-kinetic approach, the molecular movement at the wetting line is accounted. The model is developed by Blake and De Coninck (2002). They explained that in the vicinity of a moving contact point the interfacial inertial and surface tension forces are present which favor the motion of molecules in a advancing and receding direction. This model can be found in the works of Ren and Weinan, 2007, Bhardwaj and Attinger, 2008 and Bhardwaj et al. (2010). They treated contact line boundary condition by prescribing the contact line velocity dependent on the value of dynamic contact angle.
From the literature survey, it is found that most of the numerical studies are reported for non-bouncing cases on hydrophilic surfaces; there is scarcity of numerical investigation (except Mukherjee and Abraham, 2007, Caviezel et al., 2008, Sprittles and Shikhmurzaev, 2012) for the impact of a droplet on hydrophobic and superhydrophobic surfaces which results in bouncing or non-bouncing fate. Thus, with a long term objective of a numerical impact-dynamics study for the hydrophobic and superhydrophobic surfaces, the present paper was initiated with three objectives as follows:
- i.
To incorporate the modeling of contact line motion (as an extension to a series of work, Gada and Sharma, 2011, Gada and Sharma, 2012, Lakdawala et al., 2014, Lakdawala et al., 2015a, Lakdawala et al., 2015b; in our research group) in the DGLSM based in-house code in the 2D axisymmetric cylindrical coordinate system developed by Gada (2012).
- ii.
To perform the experiments on the impact of a droplet on hydrophobic and superhydrophobic surfaces, at various droplet impact velocities, and use them for the validation and a performance study of our numerical results.
- iii.
To demonstrate an effectiveness of the DGLSM as compared to the traditional LSM for the present problem.
In the present paper, a quasi-dynamic contact angle model of Fukai et al. (1995) is implemented and experimentally recorded advancing and receding contact angles are used as inputs in the model.
Section snippets
Physical and computational description of the problem
In the present study on the level-set method for the contact line modeling, the problem considered is a gravity induced impact of a spherical droplet on a flat or micropillared surface, shown in Fig. 1. The figure shows a water droplet (of initial diameter d0) impacting the surface with a velocity v0 in air. For the present non-dimensional study, d0 is the characteristic length and v0 is the velocity scale, and water (fluid 1) is considered as the reference fluid. A time wise variation of
Governing equations
The Navier–Stokes equations are coupled with the LS method (Sussman et al., 1994) by invoking the single field formulation, wherein a single velocity and pressure field is described for both the fluids. The material properties of the liquid droplet and the air are used, assuming that fluids are incompressible. The surface tension force at the interface is modeled as a volumetric source term in the momentum equation; this source term is non-zero only at the interface. Gada and Sharma (2009)
Experimental details
In this section, an experimental details performed in the present paper are described in brief. The micropillared surfaces were fabricated using an ultraviolet lithography technique (Milionis et al., 2013), after depositing SU-8(2025) epoxy polymer on a 2 inch silicon wafer. The surfaces are coated with 10 nm platinum layer. The SEM images of the fabricated micropillared surface, with 47 µm pitch, are shown in Fig. 4a. The measured width and height of the square pillars are 20 ± 2 µm and 27 ± 2 µm,
Results and discussion
This section is presented mainly in two parts. In section “Code validation study of DGLSM”, the code validation study of the DGLSM based simulations is shown with the present as well as published experimental results. In section “Performance study of the dual-grid as compared to the traditional LSM”, the performance evaluation of the DGLSM over the traditional LSM is presented. For the simulations, four cases are selected as shown in Table 1.
During the droplet dynamics, three stages are
Conclusions
A detailed numerical methodology for the modeling of contact line motion (quasi-dynamic contact angle model of Fukai et al., 1995), during the impact of a droplet on a solid surface, is presented for a dual-grid level-set method (DGLSM). Experiments are done for the validation of the DGLSM based code and also for demonstrating the better performance of the DGLSM as compared to the traditional LSM. The validation and performance study – with our experimental results – is presented for the
Acknowledgments
Authors would like to acknowledge Mr. Absar Lakdawala as well as Mr. Javed Shaikh for continuous support in the development of code and sharing ideas of implementation. The first author was supported by fellowship awarded by Industrial Research and Consultancy Centre, IIT Bombay.
References (66)
- et al.
Inertia dominated flow and heat transfer in liquid drop spreading on a hot substrate
Int. J. Heat Fluid Flow
(2011) - et al.
Non-isothermal wetting during impact of millimeter-size water drop on a flat substrate: numerical investigation and comparison with high-speed visualization experiments
Int. J. Heat Fluid Flow
(2008) - et al.
Interfacial temperature measurements, high- speed visualization and finite-element simulations of droplet impact and evaporation on a solid surface
Int. J. Heat Mass Transf.
(2010) - et al.
The influence of solid-liquid interactions on dynamic wetting
Adv. Coll. Interface Sci.
(2002) - et al.
On the diffuse interface method using a dual-resolution cartesian grid
J. Comp. Phys.
(2014) - et al.
Numerical simulation of static and sliding drop with contact angle hysteresis
J. Comp. Phys.
(2010) - et al.
Analytical and level-set method based numerical study on oil- water smooth/wavy stratified-flow in an inclined plane-channel
Int. J. Multiph. Flow
(2012) - et al.
Hydrodynamic model of steady movement of a solid/liquid/fluid contact line
J. Coll. Interface Sci.
(1971) - et al.
A dual grid level set method based study of interface-dynamics for a liquid jet injected upwards into another liquid
Int. J. Multiph. Flow
(2014) - et al.
DGLSM based study of temporal instability and formation of satellite drop in a capillary jet breakup
Chem. Eng. Sci.
(2015)
Numerical study of droplet impact and coalescence in a microline patterning process
Comput. Fluids
Sharp interface cartesian grid method ii: a technique for simulating droplet interactions with surfaces of arbitrary shape
J. Comp. Phys.
VOF simulations of the contact angle dynamics during the drop spreading: standard models and a new wetting force model
Adv. Coll. Interface Sci.
Drop collisions with simple and complex surfaces
Curr. Opin. Coll. Interface Sci.
Investigations of drop impact on dry walls with a lattice- boltzmann model
J. Coll. Interface Sci.
A three-dimensional model of droplet impact and solidification
Int. J. Heat Mass Transf.
Numerical simulation of moving contact line problems using a volume-of-fluid method
J. Comp. Physics
Singularities at the moving contact line: mathematical, physical and computational aspects
Phys. D: Nonlinear Phenom.
A level-set approach for simulations of flows with multiple moving contact lines with hysteresis
J. Comp. Phys.
A level set approach for computing solutions to incompressible two-phase flow
J. Comp. Phys.
Numerical simulation of dynamic behavior of droplet on solid surface by the two-phase lattice boltzmann method
Comput. Fluids
A front-tracking method for viscous, incompressible, multi-fluid flows
J. Comp. Phys.
A coupled quadrilateral grid level set projection method applied to ink jet simulation
J. Comp. Phys.
Phase-field simulations of interfacial dynamics in viscoelastic fluids using finite elements with adaptive meshing
J. Comp. Phys.
Controlling droplet deposition with polymer additives
Nature
Do surfaces with mixed hydrophilic and hydrophobic areas enhance pool boiling?
Appl. Phys. Lett.
Micro-, nano- and hierarchical structures for superhydrophobicity, self-cleaning and low adhesion
Philos. Trans. R. Soc. A
Self-cleaning surfaces- virtual realities
Nature Mater.
Wetting and spreading
Rev. Mod. Phys.
Modeling the splash of a droplet impacting a solid surface
Phys. Fluids
Anti-icing superhydrophobic coatings
Langmuir
Adherence and bouncing of liquid droplets impacting on dry surfaces
Microfluid. Nanofluid.
Maximal deformation of an impacting drop
J. Fluid Mech.
Cited by (34)
Droplet collision with hydrophobic and superhydrophobic surfaces: Experimental studies and numerical modeling
2024, Surfaces and InterfacesDesign principle of ridge-textured superhydrophobic surfaces for inducing pancake bouncing
2022, International Communications in Heat and Mass TransferDroplet spreading dynamics on hydrophobic textured surfaces: A lattice Boltzmann study
2021, Computers and FluidsCitation Excerpt :This phenomenon depends on two parameters, the contact time and the contact area. Most researches have focused on inducing a more rapid detachment of the impacting droplet from the surface to reduce the contact time [11-13]. Some investigations in this category indicate that the surface morphology in sub-millimeter scales plays a vital role in reducing the contact time of impact droplets by breaking the symmetry mode of the spreading process [14-16].
Anti-freezing characteristics of water droplet impinging the superhydrophobic surface: An experimental and predictive study
2021, Applied Surface ScienceCitation Excerpt :Hence, it is of great value to investigate the dynamic behaviors of droplet impacting the supercooled superhydrophobic surfaces. Over the past two decades, droplet impact has been extensively investigated, mainly concentrating on the influence of droplet size [8–10], initial impact velocity [11–15] and surface morphology and structure [16–21] on the droplet dynamics. Tran et al. [8] performed an experimental study of droplet impinging a superheated surface and reported that the contact time of droplet strongly depends on the droplet size and the maximum spreading factor Γ displays a universal scaling behavior Γ ~ We3/10.
Role of anisotropic pinning and liquid properties during partial rebound of droplets on unidirectionally structured hydrophobic surfaces
2021, Chemical Engineering ScienceCitation Excerpt :Surface wettability can be altered either by chemical means (Almohammadi and Amirfazli, 2017; Antonini et al., 2012) or by varying topography (Bandyopadhyay and Sharma, 2010; Parihar et al., 2018a; Tsai et al., 2009). Recently topographically dependent wettability of surfaces has gained importance for various applications (Kumar et al., 2019; Liu et al., 2015; Patil et al., 2016a; Yeong et al., 2014). Microstructures on flat substrates increase the surface roughness, thereby modifying the contact angle in the process and make the surfaces hydrophobic/super-hydrophobic.
Numerical investigation on impingement dynamics and freezing performance of micrometer-sized water droplet on dry flat surface in supercooled environment
2019, International Journal of Multiphase FlowCitation Excerpt :After the droplets transformed to two-phase coexistence state at higher temperature than the initial supercooling temperature, the temporal evolutions of the volume fraction of the water phase, the velocity distribution in the flow field, the wetting diameter of the impinging droplet and the surface heat exchange were tracked and discussed in this study. Since volume of fluid (VOF) method (Cimpeanu and Papageorgiou, 2018) has the advantage of good fluid conservation, which has been proved practicable in tracking the interface (Wang et al., 2018), while level-set method (Patil et al., 2016) is coupled to improve the accuracy of interfacial forces, the dynamics of the droplet were simulated through a coupled VOF and level-set method (Sussman and Puckett, 2000). As the strong compatibility of the coupled VOF and level-set method, the dynamic contact angle model and the phase-change model were both incorporated.