Next Article in Journal
Facile Synthesis of Nitrogen-Doped Carbon-Supported Rhodium–Cobalt Alloy Electrocatalyst for Oxygen Reduction Reaction
Next Article in Special Issue
Mixed Seeds of Oat and Vetch Based on DEM-Fluent Coupling Motion Simulation in a Venturi Tube
Previous Article in Journal
The Liquid Phase Oxidation of Light Hydrocarbons for Thermo-Gas-Chemical Enhanced Oil Recovery Method
Previous Article in Special Issue
Experimental and Numerical Studies of Fine Quartz Single-Particle Sedimentation Based on Particle Morphology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Lattice Boltzmann Modeling of a Sessile and a Body Force-Driven Sliding Droplet over a Grooved Surface

by
Assetbek Ashirbekov
1,
Nursultan Zhumatay
1,
Alibek Kuljabekov
1,
Bagdagul Kabdenova
1,
Ernesto Monaco
2,
Lei Wang
3 and
Luis R. Rojas-Solórzano
1,*
1
School of Engineering and Digital Sciences, Department of Mechanical and Aerospace Engineering, Nazarbayev University, Kabanbay batyr Avenue 53, Astana 010000, Kazakhstan
2
Engineering Software Steyr (ESS), Berggasse 35, 4400 Steyr, Austria
3
State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, College of Energy, Chengdu University of Technology, Chengdu 610059, China
*
Author to whom correspondence should be addressed.
Processes 2022, 10(11), 2356; https://doi.org/10.3390/pr10112356
Submission received: 29 August 2022 / Revised: 11 October 2022 / Accepted: 27 October 2022 / Published: 11 November 2022
(This article belongs to the Special Issue Computational Modeling of Multiphase Flow (II))

Abstract

:
This work presents the numerical modeling of a droplet’s sessile and dynamic behavior on a grooved surface. A droplet is placed on horizontal and vertical sliding conditions to observe its behavior under wettable and non-wettable conditions. The numerical analysis uses the multicomponent multiphase Shan-Chen Lattice Boltzmann Model (SC-LBM). The Cassie–Baxter and Wenzel states are reproduced for the sessile condition, and the enhancement of the contact angle is appreciated under the action of the grooved-ridged horizontal surface. The sliding droplet is analyzed through the Bond number by varying the ratio between the body force and the surface tension number. For Cassie–Baxter and Wenzel wettability conditions, a critical Bond number was discovered above which the sliding droplet will continue to deform indefinitely. The numerical model proved its suitability to predict the gradual deformation of a droplet over a grooved vertical surface subject to a tangential body force until the droplet eventually reaches a sessile condition or a breakup.

1. Introduction

Modeling the dynamics of a sliding droplet on a vertical surface excited by external forces such as gravity, pressure gradient, and the attractive/repulsive wetting molecular forces at the fluid–solid contact point can help to describe the processes of many practical applications such as inkjet printing [1], enhanced oil recovery [2], and film coating [3].
However, most substrate surfaces are not perfectly smooth and present irregularities that can dramatically change the dynamic behavior of the sliding droplet. Furthermore, some applications require engineered surfaces corrugated with periodic groove ridges to enhance their wettability properties [4]. Therefore, it is essential to understand how the wettability of such grooved-ridged surfaces is affected by their primary dimensions and spacing. Hence, this study examines the droplet movement and deformation on a grooved-ridged vertical surface under varying wettability conditions.
Two primary states, namely Cassie–Baxter [5] and Wenzel [6], have been historically identified on grooved-ridged surfaces. Both are produced by minimizing the system’s free energy under various conditions. The droplet contact line advances over the surface in the Wenzel state, filling the substrate groove voids between the ridges. The surface roughness enlarges the effective surface area of the wall, and the contact angle measured on a flat surface is adjusted according to [5]:
c o s θ W e n z e l   = r   c o s   θ f l a t
where the surface roughness is accounted for by r, which is generally more than 1. Thus, Equation (1) illustrates how the surface roughness affects the apparent wettability of a given surface under the presence of a sessile droplet in the Wenzel state. On the other hand, Cassie and Baxter [5] suggested a correlation between the apparent contact angle and surface roughness when the droplet does not fill the cavities. Their proposal was based on the hypothesis that the droplet constantly rests on top of the roughness. The apparent contact angle is determined by Equation (2) [6]:
c o s θ C a s s i e B a x t e r   = χ c o s   θ f l a t ( 1 f )
where χ is the fraction of the area of the surface moistened with a drop.
Droplets in the Cassie–Baxter and Wenzel states are ubiquitous in nature and engineering applications. For example, Cassie–Baxter droplets are frequently spotted on hydrophobic surfaces such as butterfly wings and lotus leaves [7]. Surfaces with the Cassie–Baxter state have numerous desirable qualities, including low adhesion, water repellency, and high droplet mobility, which may be utilized to minimize drag in microfluidic devices, fuel cell chips, and self-cleaning surfaces [8]. On the other hand, irregularities on hydrophilic surfaces increase their capacity to entrain droplets in the Wenzel state, which is ideal for applications requiring good solid–liquid contacts, such as heat control microfluidics and pesticide retention on plants [9]. Therefore, it is necessary to study the behavior of drops in the Cassie–Baxter, and Wenzel states to maximize the benefits of rough surfaces in engineering applications. Numerical simulations can help us design microfluidic grooved-ridged surfaces by controlling the different wettability conditions under horizontal and vertical configurations, provided the numerical model can capture the rich physics behind the interface motion and deformation of the sliding droplet.
The Lattice Boltzmann method (LBM) has emerged as a suitable numerical tool for modeling fluid flows. The basic idea behind the Lattice Boltzmann technique is to build simpler kinetic models that include the critical physics of mesoscopic or microscopic processes so that the macroscopic averaged attributes also reproduce macroscopic continuum equations. In addition, the straightforward implementation of boundary conditions, programming simplicity, and the ability to incorporate microscopic interactions are advantages that the LBM provides [10]. The LBM has been used in previous studies of multiphase and multicomponent flows in porous media and single domains [11,12,13,14,15]. For example, Suo, Liu, and Gan [16] studied the transport and spreading of droplets through porous media. Wang and Wang [17] successfully investigated the movement of a droplet along a filament under the effect of gravity force. Additionally, Chen et al. [18] studied the spreading of a 3D droplet on an inclined wall using the LBM with a single-component multiphase pseudopotential model. Bhardwaj and Dalal [19] did a mesoscopic analysis of 3D droplet displacement on wetted grooved wall microchannel under the gravity influence. Hyväluoma et al. [20] analyzed a 3D droplet sliding on a micro-grooved wall under various angles, utilizing SC-LBM. However, the variation of geometric parameters of the wall and modeling the detachment were not considered. Zhou et al. [21] performed a similar 2D study featuring a flow model in a channel lined with grooves. Vrancken et al. [22] studied the reversible transition of a droplet between Cassie–Baxter and Wenzel states using the free-energy LBM model with experimental validation, but the droplet detachment was not considered. Droplet detachment was considered instead by Kang, Zhang and Chen [23] in a similar study featuring 3D SC-LBM, analyzing the effect of contact angle, but only with smooth walls. Wang, Deng and Zhang [24] developed a model of droplet detachment from a grooved surface, but in a downward-horizontal orientation, with a droplet dripping without controlling the Bond number. Despite the 3D approach being a more realistic depiction of a physical system, 2D models have proven to be a valid alternative in preliminary analyses successfully describing essential features of droplet dynamics interacting with solid surfaces.
Thus, in the present study, we used a 2D SC-LBM, following our earlier work [25], to study droplet spreading on a horizontal grooved surface and, for the first time to our knowledge, the droplet detachment from a vertical grooved surface as a case of practical interest. Cassie–Baxter and Wenzel’s states were explored, primarily focusing on observing the contact angle dynamics under the interaction of solid–fluid adhesion forces, driving body force intensity, and fluid–fluid surface tension.

2. Methodology

The LBM formulation considers the fluid as a set of fictitious particle packages interacting and moving according to probability distribution functions. The particles collide instantaneously and stream in a predetermined pattern and speed to neighboring nodes every time step in a Cartesian grid (i.e., a lattice). A single probability distribution function describes each fluid particle’s physical condition at a particular time step in a given node. After the particles collide, the system transits to reach its local equilibrium condition in a relaxation time before streaming. As a result, every probability distribution function in the model is associated with a fluid component and a given streaming velocity vector, governed by the following Lattice Boltzmann Equation (3):
f α , i ( x + e i Δ t , t + Δ t ) f α , i ( x , t ) =   1 τ α   [ f α , i e q ( x , t ) f α , i ( x , t ) ]
where f α , i   represents the particle’s probability distribution function for the fluid component “α” in the “ith” direction,   f α , i e q ( x , t )   is the equilibrium distribution function in the corresponding lattice velocity direction and fluid component, e i   is the microscopic velocity in the ith direction, Δ t   is the time-step, and τ is the relaxation time.
The collision operator on the right side of Equation (3) is based on the Bhatnagar–Gross–Krook (BGK) approximation [10], and the streaming step is on the left side. Because the collision is only local to a single node, the LBM is suitable for massively parallel calculations on many CPU or GP-GPU systems. On the other hand, streaming provides data transfer from one node to the following [26]. The dimensionless relaxation time of component “α”, τ α , is connected with its kinematic viscosity by υ α   = c s 2   ( τ α −0.5 Δ t ).
In the present work, the droplet analysis will be undertaken as a fluid immersed in another fluid phase. The proposed two-component multiphase model is the Shan-Chen (SC) LBM with a D2Q9 lattice scheme, i.e., two-dimensional nine-velocity directions.
For the current suggested model, the equilibrium distribution function is:
f α , i e q ( x , t ) = ω i ρ α   [ 1 + 3 ( e i u α e q ) + 9 2 ( e i u α e q ) 2 3 2 ( u α e q ) 2 ]
where u and ρ represent the macroscopic velocity and density, and ω i designates the weighting factor. The weighting factor ω i for i = 0 is 4/9, while i = 1, 2, 3, and 4 is 1/9, and i = 5, 6, 7, and 8 is 1/36. The moments of the probability distribution functions are shown in Equations (5) and (6) to calculate the density   ρ α and momentum density   ρ α u α .
ρ α = i f α , i  
ρ α u α = i f α , i e i
The cohesive force between the two components is given by the Equation (7) [27]:
  F α ,   C ( x , t ) = G α α ˜ , C ρ α ( x , t ) i ω i ρ α ( x + e i Δ t , t ) e i
where G C represents the inter-particle interaction strength, either positive or negative depending on the nature of the resulting force, i.e., attraction or repulsion, respectively. If the threshold point of interaction strength is surpassed, phase separation occurs spontaneously. The behavior of a non-ideal fluid is represented using an interaction parameter termed effective mass, which is meant to replicate the function of a Lennard-Jones potential [28]. The density plays this role in Equation (7).
Moreover, Equation (8) defines the adhesive force between the solid substrate and the fluid [29]:
F α ,   a d s ( x , t ) = - G α , a d s   ρ α   ( x , t ) i ω i s ( x + e i Δ t , t ) e i
where s = 0 denotes a fluid node, s = 1 denotes a solid substrate, and G a d s is the wettability factor, also known as the solid–fluid interaction parameter, which indicates the strength of the interaction between the solid and fluid. A non-wetting fluid–surface contact has a positive G a d s , while a wetting fluid–surface interaction has a negative G a d s [26]. A body force is used to create a flow in vertical cases:
F α ,   b   ( x , t ) = g ( a f a ( x ) )  
where g is body force per unit mass. Throughout every spatial point on the grid nodes, the forces stated above operate on the respective fluid component.
The effect of forces is accounted for by calculating the velocity under the influence of forces u α e q , as shown in Equation (10), a component of Equation (4).
u α e q   = u α + τ α ρ α   (   F α ,   C + F α ,   a d s + F α ,   b   ) Δ t

3. Verification and Validation of the LBM Model

The present study continues our previous work [25], where the earlier described LBM model was validated on a 2D sessile droplet and a two-fluid stratified Poiseuille flow. The present study further analyzes the droplet spreading on a roughened surface represented by horizontal and vertical groove-ridge patterned surfaces. Firstly, the static contact angle on the groove-ridge surface was measured and compared with the theoretical contact angle determined by Equations (1) and (2) to prove the validity of the LBM code used in the study. The equivalent surface roughness r was defined by   r = [ ( w + 2 s ) 2 + 4 w h ]   ( w + 2 s ) 2 ; where w and h are the ridge width and height, and s denotes the spacing between the ridges, as shown in Figure 1. The static contact angle is prescribed for the flat surface, as in our previous work [25]. The solid fraction χ in Equation (2) was determined by χ = w 2 ( w + 2 s ) 2 . As explained in the next section, changing the groove-ridge dimensions triggered the Cassie–Baxter and Wenzel wetting states.
The horizontal grooved-ridged surface is firstly studied without including gravitational effects. As indicated in our previous work [25], a 2D liquid droplet is placed at the chamber’s center, surrounded by an immiscible background fluid. The domain geometry, initial, and boundary conditions are shown in Figure 2, with the weightless droplet in its initial state. Hence, a droplet with a radius of 20 lattice units (lu) was placed in the middle between two contiguous ridges. The no-slip boundary condition (BC) is imposed on the top and bottom walls, whereas periodic BCs are prescribed on the left and right sides of the domain. The no-slip boundary condition is applied by reversing the distribution function on each wall node x w , f i ( x w , t ) = f j ( x w , t ) , where j is the direction opposite to the initial direction i, so e i   = e j . The periodic boundary condition is applied by translating the values of the probability distribution function at the boundary node x p to the opposing side node x o , preserving the directions, f i ( x p , t ) = f i ( x o , t ) .
The size of the initial computational domain is 151 × 101 lu2. This mesh configuration was chosen as an initial set-up because it was used in our previous work [25] for correctly reproducing static contact angles.
The height (h) and width (w) of the ridges, and the spacing (s) between the two ridges, were varied to observe the droplet spreading over the ribbed surface to depict either the Wenzel and Cassie–Baxter states. Geometrical considerations were used to determine the droplet contact angle once it reached equilibrium, as reported in [25]. The SC-LBM model showed its excellent capability in recreating these classical hydrophilic or hydrophobic states.
The first set of geometrical parameters used was w = h = s = 5 lu, with a roughness of r = 1.44. Wenzel’s wettability state was attained from this initial simulation. For the initial mesh, Wenzel’s state contact angle obtained from the numerical simulation is 69.3°, close to the theoretical value of 69.6° calculated via Equation (1).
The required verification of the domain discretization was conducted before validating the numerical findings using a progressive grid refinement method. The grid highlighted on top of the grooved wall geometry is shown in Figure 3. The domain height (y-direction) was discretized with a fixed medium resolution of 101 lu-grid because the vertical grid size was not expected to affect Wenzel’s contact angle, as it greatly depends on the x-discretization. As a result, the horizontal grid size (x-direction) was the focus of the grid verification. The domain’s x-grid was refined by 33% twice (see Table 1). The contact angle at the Wenzel state between the first and second grid had a relative error of 0.14% and a much smaller error between the second and third grids; hence, the computational domain with a size of 151 × 101 lu2 was chosen for subsequent static wetting regime simulations.
The grid verification and validation of the LBM platform for the droplet contact angle at the Wenzel state were completed when the final contact angle in Table 1 (69.4°) differed by 0.29% from the theoretical value (69.6°).
Table 2 shows how the surface roughness with the Wenzel state enhances the hydrophilic/hydrophobic nature of the surface.
The Cassie–Baxter state was obtained at geometrical parameters of w = h = 10 lu, s = 5 lu, and a solid fraction χ of 0.25, as shown in Table 3.
The assessment of the numerical results revealed that the SC-LBM could accurately represent the wettability effects for grooved horizontal surfaces over a wide range of contact angles. Table 2 and Table 3 showed that all cases’ errors between numerical and analytical contact angles for Cassie–Baxter and Wenzel wettability states are lower than 1%.

4. Transition from Cassie–Baxter to Wenzel State

As indicated above, the SC-LBM model is capable of recreating the classical Cassie–Baxter and Wenzel states. However, it is crucial to investigate which of the three geometrical parameters (width, height, spacing) directly impacts the contact angle at the Cassie–Baxter and Wenzel wettability states. Thus, the next step is to determine under what condition the transition between the two states occurs. According to Zhumatay et al. [25], the Wenzel wettability state occurs when the sum of geometrical parameters is lower or equal to the droplet radius. Additionally, according to Yagub et al. [30], the wettability model is defined by the aspect ratio R, representing the ratio of the ridge width and the spacing between two consecutive ridges. Zheng et al. [31] revealed that a high aspect ratio promotes the Cassie–Baxter wettability condition. In contrast, a decrease in the aspect ratio hastens the transition from the Cassie–Baxter state to the Wenzel state.
In the present study, the aspect ratio featured in the initial set-up was equal to 1 (w = 5 lu, s = 5 lu), and the Wenzel state was obtained. This outcome was expected since the sum of groove-ridge geometrical parameters perfectly matches the theoretical expectation [25]. Therefore, the width and spacing were increased to 10 lu and 6 lu, respectively, leading to an aspect ratio of 1.67. The new groove-ridge dimensions added up more than the droplet radius (10 + 6 + 5 = 21 lu > 20 lu), causing a metastable Cassie–Baxter state (see Figure 4) after oscillations between the Cassie–Baxter and Wenzel states for several thousands of time-steps. Further, the width was increased to 15 lu, and the spacing decreased to 5 lu, thus, resulting in an aspect ratio of 3 with unchanged height. The metastable Cassie–Baxter state settled up again, meaning that the occurrence of this state does not depend solely on the aspect ratio but also on the size of the ridge. Therefore, the height was increased to 10 lu, and the width was decreased to 10 lu; although spacing has not changed, the droplet ends up in a classical Cassie–Baxter state. The height was reduced to 9 lu while keeping the other parameters (s = 5 lu, w = 10 lu) to determine the critical value of the ridge height that can produce a Cassie–Baxter state. However, the metastable Cassie–Baxter state was obtained. Based on these findings, it is reasonable to deduce that the following three conditions must be met to produce the Cassie–Baxter state for a droplet immersed in a background fluid with unit density and viscosity ratios:
  • The sum of geometrical parameters should be higher than the droplet radius;
  • The aspect ratio should be higher or equal to 1;
  • The height of the ridge should be at least half of the droplet radius.
Naturally, these results only apply to the current model with the rectangular groove-ridge surface roughness.

5. Dynamic Regime: Droplet Sliding over a Grooved-Ridged Surface

The SC-LBM mimics the dynamic regime of a sliding droplet by putting a 2D droplet under the influence of a body force. This force models the effect of gravity and pressure gradient onto a vertical wall with rectangular groove ridges. It lets the droplet slide, starting from a weightless initial static condition, subject to different wetting regimes. Once the weightless droplet settles down on the vertical wall, the body force enters into action in the vertical direction. Fluid 1 designates the liquid droplet throughout the analysis, while the surrounding background fluid is prescribed as Fluid 2. The density and viscosity ratios between fluids were set to unity throughout this study.

5.1. Verification of the Mesh

Similar to previous simulations (see Figure 2), non-slip boundary conditions are prescribed on the walls, whereas periodic boundary conditions are set on the entry and exit of the computational domain.
The surrounding background fluid is initially considered at rest, with the sliding droplet dynamics acting as its only drive of motion. The final computational domain’s x-y resolution was 101 × 751 lu2, for which the width (y-direction, aligned with body force vector parallel to the grooved surface) was verified using two grid sizes: 101 × 801 lu2 and 101 × 851 lu2. When the steady sliding condition was obtained after 8000 time-steps, the relative inaccuracy of the wet length was determined. Between y-grid sizes 101 × 751 and 101 × 801 lu2, the relative error of the wet length was equal to 1.38%. A further drop to 0.71% in the relative error was attained when comparing y-grid sizes of 101 × 801 vs. 101 × 851 lu2, allowing us to adopt 801 lu as the resolution in the y-direction. Additionally, adjusting the x-grid resolution (transversal to body force vector) from 101 to 151 lu and 151 to 201 lu had little effect on the wet length, with the relative error remaining below 1%. Therefore, the final grid was set to 101 × 801 lu2.

5.2. Sliding Droplet Analysis

The Bond number, defined in Equation (11), captures the relationship between body force and surface tension:
B 0   = ρ g l 2 σ 12
where ρ is the density of the droplet [32], g is the body force per unit mass, l is the characteristic length, which is equal to the diameter of the droplet, and σ 12 is the interfacial tension between the two fluids. As specified in [25], σ 12 is equal to 68.2 × 10−4 lattice units. Hence, in our analysis, the Bond number represents the ratio between the driving body force of the fluid through the channel, which may denote either the gravity or a driving pressure gradient divided by the surface tension between the two fluids.
In the present investigation, we varied the Bond number between 0.526 and 1.29 in both Cassie–Baxter and Wenzel states for the static contact angle of 89.3° (expected on a flat surface).
The first case study corresponds to the liquid droplet with a radius of 20 lu under the g of 2.2 × 10−6 lu/ts2 and a Bond number Bo = 0.526. The contact line goes to the next ridge once the receding contact angle achieves its limited value. At the same time, the droplet deforms further until it reaches a balance between capillary and body force (see Figure 5).
Next, the body force component was varied to adjust the Bond number. At Bond number Bo = 0.704, the movement of the receding contact line from one ridge to another is observed. Once the advancing contact angle reaches the hysteresis window, the droplet motion is prompted. After the droplet motion starts, the receding and advancing contact lines jump from one ridge to the next; in the meantime, the droplet continues to slide over the grooved surface (see Figure 6).
The droplet dynamics dramatically change by increasing the Bond number to 1.29. The droplet contact region with the surface reduces until the droplet fully detaches from the bottom surface. The advancing contact line moves slower than the receding contact line. This state corresponds to a detachment regime (see Figure 7).
Changing the groove-ridge dimensions to drive the Wenzel state at the Bond number of 0.821 allowed us to determine if any substantial difference was appreciated concerning the Cassie–Baxter state. The Wenzel state increases the wall adhesion to the droplet while only deforming it slightly. The contact line remains static on the substrate (see Figure 8).
The droplet increases its deformation and starts slipping onto the surface after increasing the Bond number to 0.891. Furthermore, the advancing and receding contact lines move at nearly the same rate, representing a steady state within the fluid system (see Figure 9).
Further increase of the Bond number to 1.05 causes a droplet to stretch until it forms a neck and splits into two separate parts. The surface’s grooves-ridges still trap one part of the droplet, and the remaining part separates from the substrate and descends rapidly (see Figure 10). Additionally, an entrapment of the fluid in the space between ridges has been observed for slipping and breakup regimes. A similar observation has been noted in studies by Bhardwaj and Dalal [19].
Therefore, as the Bond number increases, the sliding droplet at Cassie–Baxter and Wenzel states experiences three regimes: stationary, slipping, and detachment. The stationary regime is reached when the contact line remains immobile on the substrate. As expected, the detachment regime occurs at a smaller critical Bond number under the Cassie–Baxter state compared to the Wenzel state.

6. Impact of the Density Ratio

Six additional simulations were performed for different wettability states to determine how the density ratio affects Wenzel/Cassie–Baxter states and the critical Bond number. In this part of the study, density ratios (Fluid 1/Fluid 2, Fluid 1: liquid; Fluid 2: background phase) of 0.5 and 4 were assessed for each wettability condition. Three different geometrical configurations associated with Wenzel/Cassie–Baxter and metastable Cassie–Baxter states were explored. Similar results as presented here for the unit density ratio were obtained and, thus, refrained from showing them in the present publication. However, one remarkable observation is that the higher density ratio requires more time-steps to stabilize the Wenzel or Cassie–Baxter states. The density ratio of 0.5 requires twice as fewer time-steps as needed for a density ratio of 1 and quadruple fewer than those required for a density ratio of 2. These results show that the density ratio is not a limitation of the current work but affects the time it takes to reach stable states.

7. Concluding Remarks

The multiphase SC-LBM model of a 2D weightless sessile droplet spreading on a horizontal grooved surface and its body force-driven slipping on a vertical configuration is presented. Results were obtained for unit density and viscosity ratios between the droplet and the surrounding fluid. Several wettability regimes were explored by varying the groove-ridge dimensions and the static contact angle. The droplet depicted the Wenzel or Cassie–Baxter states in all studied scenarios, even on the vertical wall. The sliding droplet on the vertical wall, driven by a parallel body force, portrayed dynamic contact angles affected by the groove-ridge dimensions and exhibited stationary, slipping, and breakup regimes under both states.
The Bond number was varied from 0.526 to 1.29 and determined the evolution of the sliding droplet, revealing a critical Bond number at which detachment from the wall may occur. As expected for this wetting condition, the critical transitional Bond number was smaller for the Cassie–Baxter state than for the Wenzel state in response to the lower effective contact surface between the droplet and the wall in the former.
Further analysis of the impact of the density ratio, changing it between 0.5 and 4, demonstrated that it does not affect the observed wettability regimes and their critical Bond numbers.
Future work is ongoing to investigate the sliding droplet subject to variations in the geometry and dimensions of various roughness elements. Additionally, the study will be extended to 3D configurations to explore more complex and realistic slipping conditions.

Author Contributions

Conceptualization, A.A., N.Z. and L.R.R.-S.; methodology, A.A., N.Z., A.K., B.K., E.M., L.W. and L.R.R.-S.; software, A.A., N.Z., A.K., B.K., E.M. and L.R.R.-S.; validation, N.Z., L.W. and L.R.R.-S.; formal analysis, A.A., N.Z., A.K., B.K., E.M., L.W. and L.R.R.-S.; investigation, A.A., N.Z., B.K., A.K., E.M., L.W. and L.R.R.-S.; writing—original draft preparation, A.A., N.Z., B.K., E.M., L.W. and L.R.R.-S.; writing—review and editing, A.A., N.Z., B.K., E.M., L.W. and L.R.R.-S.; visualization, A.A., N.Z., B.K., E.M., L.W. and L.R.R.-S.; supervision, L.R.R.-S.; project administration, L.R.R.-S.; funding acquisition, L.R.R.-S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Nazarbayev University PhD scholarship of Bagdagul Kabdenova and M.Sc. scholarships of Assetbek Ashirbekov and Nursultan Zhumatay, funded by the Ministry of Education and Science of the Republic of Kazakhstan. This work is funded by the NU Faculty Development Competitive Research Grant “Numerical simulation of Carbon Dioxide Sequestration in porous media using multiphase/multicomponent thermal Lattice Boltzmann Model”, grant number #021220FD0951.

Data Availability Statement

The data presented in this study are available in the displayed figures and tables within the article but can also be obtained on request from the corresponding author within a reasonable time-frame.

Acknowledgments

Open-source platform DL_MESO LBM package https://www.scd.stfc.ac.uk (accessed on 1 January 2022) is used as a base to build a new, improved model. The authors acknowledge M. Seaton for providing the original version of the code.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

SymbolMeaning
BoBond number
f Particle’s probability distribution function
f e q Equilibrium distribution function
F Force
gBody force per unit mass
h Ridge height
l Droplet diameter at initial state
r Surface roughness
s Spacing between ridges
u Particle velocity
wRidge width
θ Contact angle
ρ Density
σ Interfacial tension
χ Solid fraction
SubscriptsMeaning
a d s Adhesion
b Body
C Cohesion
i Direction
α Species

References

  1. Wijshoff, H. Drop dynamics in the inkjet printing process. Curr. Opin. Colloid Interface Sci. 2018, 36, 20–27. [Google Scholar] [CrossRef]
  2. Perazzo, A.; Tomaiuolo, G.; Preziosi, V.; Guido, S. Emulsions in porous media: From single droplet behavior to applications for oil recovery. Adv. Colloid Interface Sci. 2018, 256, 305–325. [Google Scholar] [CrossRef] [PubMed]
  3. Christodoulou, C.; Sorensen, E.; García-Muñoz, S.; Mazzei, L. Mathematical modelling of water absorption and evaporation in a pharmaceutical tablet during film coating. Chem. Eng. Sci. 2018, 175, 40–55. [Google Scholar] [CrossRef]
  4. Shi, Y.; Tang, G. Relative permeability of two-phase flow in three-dimensional porous media using the lattice Boltzmann method. Int. J. Heat Fluid Flow 2018, 73, 101–113. [Google Scholar] [CrossRef]
  5. Cassie, A.; Baxter, S. Wettability of porous surfaces. Trans. Faraday Soc. 1944, 40, 546. [Google Scholar] [CrossRef]
  6. Wenzel, R. Resistance of solid surfaces to wetting by water. Ind. Eng. Chem. 1936, 28, 988–994. [Google Scholar] [CrossRef]
  7. Quere, D.; Reyssat, M. Non-adhesive lotus and other hydrophobic materials. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 2008, 366, 1539–1556. [Google Scholar] [CrossRef]
  8. Zheng, Y.; Gao, X.; Jiang, L. Directional adhesion of superhydrophobic butterfly wings. Soft Matter 2007, 3, 178–182. [Google Scholar] [CrossRef]
  9. Tropmann, A.; Tanguy, L.; Koltay, P.; Zengerle, R.; Riegger, L. Completely superhydrophobic PDMS surfaces for microfluidics. Langmuir 2012, 28, 8292–8295. [Google Scholar] [CrossRef]
  10. Perumal, D.A.; Dass, A.K. A Review on the development of lattice Boltzmann computation of macro fluid flows and heat transfer. Alexandria Eng. J. 2015, 54, 955–971. [Google Scholar] [CrossRef]
  11. Salehi, M.S.; Afshin, H.; Firoozabadi, B. Numerical Investigation of Dynamics of Drop Motion Using Lattice Boltzmann Method. J. Mol. Liq. 2019, 296, 111926. [Google Scholar] [CrossRef]
  12. Deng, H.; Jiao, K.; Hou, Y.; Park, J.; Du, Q. A lattice Boltzmann model for multicomponent two-phase gas-liquid flow with realistic fluid properties. Int. J. Heat Mass Trans. 2019, 128, 536–549. [Google Scholar] [CrossRef]
  13. Zarei, A.; Karimipour, A.; Isfahani, A.H.M.; Tian, Z. Improve the Performance of Lattice Boltzmann Method for a Porous Nanoscale Transient Flow by Provide a New Modified Relaxation Time Equation. Physics A 2019, 535, 122453. [Google Scholar] [CrossRef]
  14. Zhang, D.; Li, S.; Jiao, S.; Shang, Y.; Dong, M. Relative Permeability of Three Immiscible Fluids in Random Porous Media Determined by the Lattice Boltzmann Method. Int. J. Heat Mass Trans. 2019, 134, 311–320. [Google Scholar] [CrossRef]
  15. Son, S.; Chen, L.; Derome, D.; Carmeliet, J. Numerical study of gravity-driven droplet displacement on a surface using the pseudopotential multiphase lattice Boltzmann model with high density ratio. Comput. Fluids 2015, 117, 42–53. [Google Scholar] [CrossRef]
  16. Suo, S.; Liu, M.; Gan, Y. An LBM-PNM framework for immiscible flow: With applications to droplet spreading on porous surfaces. Chem. Eng. Sci. 2020, 218, 115577. [Google Scholar] [CrossRef]
  17. Wang, X.; Wang, L. Numerical simulations of a droplet slipping along a filament with Lattice Boltzmann Method. Eur. J. Mech. B/Fluids 2019, 78, 230–244. [Google Scholar] [CrossRef]
  18. Chen, W.; Yang, F.; Yan, Y.; Guo, X.; Dai, R.; Cai, X. Lattice Boltzmann simulation of the spreading behavior of a droplet impacting on inclined solid wall. J. Mech. Sci. Technol. 2018, 32, 2637–2649. [Google Scholar] [CrossRef]
  19. Bhardwaj, S.; Dalal, A. Mesoscopic analysis of three-dimensional droplet displacement on wetted grooved wall of a rectangular channel. Eur. J. Mech. B/Fluids 2018, 67, 35–53. [Google Scholar] [CrossRef]
  20. Hyväluoma, J.; Koponen, A.; Raiskinmäki, P.; Timonen, J. Droplets on inclined rough surfaces. Eur. Phys. J. E 2007, 23, 289–293. [Google Scholar] [CrossRef]
  21. Zhou, W.; Yan, Y.; Liu, X.; Liu, B. Lattice Boltzmann parallel simulation of microflow dynamics over structured surfaces. Adv. Eng. Softw. 2017, 107, 51–58. [Google Scholar] [CrossRef] [Green Version]
  22. Vrancken, R.J.; Kusumaatmaja, H.; Hermans, K.; Prenen, A.M.; Pierre-Louis, O.; Bastiaansen, C.W.; Broer, D.J. Fully reversible transition from Wenzel to Cassie—Baxter states on corrugated superhydrophobic surfaces. Langmuir 2010, 26, 3335–3341. [Google Scholar] [CrossRef] [PubMed]
  23. Kang, Q.; Zhang, D.; Chen, S. Displacement of a three-dimensional immiscible droplet in a duct. J. Fluid Mech. 2005, 545, 41–66. [Google Scholar] [CrossRef]
  24. Wang, H.; Deng, Z.; Zhang, C. Droplet Detachment during Condensation Process on a Downward-facing Grooved Surface. In Proceedings of the 5th International Seminar on ORC Power Systems, Athens, Greece, 9–11 September 2019. [Google Scholar]
  25. Zhumatay, N.; Kabdenova, B.; Monaco, E.; Rojas-Solórzano, L.R. Lattice Boltzmann modeling of a gravity-driven sliding droplet under a dynamic wetting regime. Eur. J. Mech. B/Fluids 2021, 86, 198–209. [Google Scholar] [CrossRef]
  26. Dauyeshova, B.; Rojas-Solórzano, L.R.; Monaco, E. Numerical Simulation of Diffusion Process in T-Shaped Micromixer Using Shan-Chen Lattice Boltzmann Method. Comput. Fluids 2018, 167, 229–240. [Google Scholar] [CrossRef]
  27. Martys, N.S.; Chen, H. Simulation of Multicomponent Fluids in Complex Three-Dimensional Geometries by the Lattice Boltzmann Method. Phys. Rev. E 1996, 53, 743. [Google Scholar] [CrossRef]
  28. Monaco, E.; Brenner, G.; Luo, K.H. Numerical Simulation of the Collision of Two Microdroplets with a Pseudopotential Multiple-Relaxation-Time Lattice Boltzmann Model. Microfluid. Nanofluid 2014, 16, 329–346. [Google Scholar] [CrossRef]
  29. Li, Q.; Luo, K.H.; Kang, Q.J.; Chen, Q. Contact Angles in the Pseudopotential Lattice Boltzmann Modeling of Wetting. Phys. Rev. E 2014, 90, 053301. [Google Scholar] [CrossRef] [Green Version]
  30. Yagub, A.; Farhat, H.; Kondaraju, S.; Singh, T. A lattice Boltzmann model for substrates with regularly structured surface roughness. J. Comput. Phys. 2015, 301, 402–414. [Google Scholar] [CrossRef]
  31. Zheng, R.; Liu, H.; Sun, J.; Ba, Y. Droplet hysteresis investigation on non-wetting striped textured surfaces: A lattice Boltzmann study. Physics A 2014, 411, 53–62. [Google Scholar] [CrossRef]
  32. Kang, Q.; Zhang, D.; Chen, S. Displacement of a two-dimensional immiscible droplet in a channel. Phys. Fluids 2002, 14, 3203–3214. [Google Scholar] [CrossRef]
Figure 1. Main dimensions of the surface grooved-ridged roughness.
Figure 1. Main dimensions of the surface grooved-ridged roughness.
Processes 10 02356 g001
Figure 2. Initial state of the liquid droplet (red) within a background fluid (white) is depicted using computational domain and boundary conditions.
Figure 2. Initial state of the liquid droplet (red) within a background fluid (white) is depicted using computational domain and boundary conditions.
Processes 10 02356 g002
Figure 3. Grid shown near the wall and ridges (h = w = 10 lu, s = 5 lu). Solid fraction shown in blue, fluid in black.
Figure 3. Grid shown near the wall and ridges (h = w = 10 lu, s = 5 lu). Solid fraction shown in blue, fluid in black.
Processes 10 02356 g003
Figure 4. Transition from the Cassie–Baxter state to Wenzel state or so-called metastable Cassie–Baxter state (w = 10 lu, s = 6 lu, h = 5 lu). Images are taken at (a) 23,000 time-steps; (b) 33,000 time-steps; and (c) 70,000 time-steps.
Figure 4. Transition from the Cassie–Baxter state to Wenzel state or so-called metastable Cassie–Baxter state (w = 10 lu, s = 6 lu, h = 5 lu). Images are taken at (a) 23,000 time-steps; (b) 33,000 time-steps; and (c) 70,000 time-steps.
Processes 10 02356 g004aProcesses 10 02356 g004b
Figure 5. Stationary regime of the Cassie–Baxter state (Bo = 0.526; h = w = 10 lu, s = 5 lu, R = 20 lu). Images are taken at (a) initial time; (b) 2500 time-steps; (c) 5500 time-steps; (d) 7000 time-steps.
Figure 5. Stationary regime of the Cassie–Baxter state (Bo = 0.526; h = w = 10 lu, s = 5 lu, R = 20 lu). Images are taken at (a) initial time; (b) 2500 time-steps; (c) 5500 time-steps; (d) 7000 time-steps.
Processes 10 02356 g005
Figure 6. Slipping regime of the Cassie–Baxter state (Bo = 0.704; h = w = 10 lu, s = 5 lu, R = 20 lu). Images are taken at (a) initial time; (b) 5500 time-steps; (c) 6500 time-steps; (d) 12,500 time-steps; (e) 14,000 time-steps.
Figure 6. Slipping regime of the Cassie–Baxter state (Bo = 0.704; h = w = 10 lu, s = 5 lu, R = 20 lu). Images are taken at (a) initial time; (b) 5500 time-steps; (c) 6500 time-steps; (d) 12,500 time-steps; (e) 14,000 time-steps.
Processes 10 02356 g006
Figure 7. Detachment regime of the Cassie–Baxter state (Bo = 1.29; h = w = 10 lu, s = 5 lu, R = 20 lu). Images are taken at (a) initial time; (b) 2000 time-steps; (c) 6500 time-steps; (d) 9000 time-steps; (e) 10,000 time-steps.
Figure 7. Detachment regime of the Cassie–Baxter state (Bo = 1.29; h = w = 10 lu, s = 5 lu, R = 20 lu). Images are taken at (a) initial time; (b) 2000 time-steps; (c) 6500 time-steps; (d) 9000 time-steps; (e) 10,000 time-steps.
Processes 10 02356 g007
Figure 8. Stationary regime of the Wenzel state (Bo = 0.821; h = w = s = 5 lu, R = 20 lu). Images are taken at (a) initial time; (b) 2000 time-steps; (c) 5500 time-steps; (d) 8500 time-steps.
Figure 8. Stationary regime of the Wenzel state (Bo = 0.821; h = w = s = 5 lu, R = 20 lu). Images are taken at (a) initial time; (b) 2000 time-steps; (c) 5500 time-steps; (d) 8500 time-steps.
Processes 10 02356 g008
Figure 9. Slipping regime of the Wenzel state (Bo = 0.891; h = w = s = 5 lu, R = 20 lu). Images are taken at (a) initial time; (b) 2000 time-steps; (c) 5500 time-steps; (d) 12,000 time-steps; (e) 14,000 time-steps; (f) 16,000 time-steps.
Figure 9. Slipping regime of the Wenzel state (Bo = 0.891; h = w = s = 5 lu, R = 20 lu). Images are taken at (a) initial time; (b) 2000 time-steps; (c) 5500 time-steps; (d) 12,000 time-steps; (e) 14,000 time-steps; (f) 16,000 time-steps.
Processes 10 02356 g009
Figure 10. Breakup regime of the Wenzel state (Bo = 1.05; h = w = s = 5 lu, R = 20 lu). Images are taken at (a) initial time; (b) 2000 time-steps; (c) 7500 time-steps; (d) 13,000 time-steps; (e) 16,000 time-steps; (f) 18,000 time-steps.
Figure 10. Breakup regime of the Wenzel state (Bo = 1.05; h = w = s = 5 lu, R = 20 lu). Images are taken at (a) initial time; (b) 2000 time-steps; (c) 7500 time-steps; (d) 13,000 time-steps; (e) 16,000 time-steps; (f) 18,000 time-steps.
Processes 10 02356 g010
Table 1. Verification of the domain resolution based on the contact angle.
Table 1. Verification of the domain resolution based on the contact angle.
Size of the Domain (lu2)Contact Angle (°)Relative Error (%)
151 × 10169.3-
201 × 10169.40.14
251 × 10169.4~0.1
Table 2. Contact angles at Wenzel state: numerical versus analytical (w = 5 lu, s = 5 lu, h = 5 lu).
Table 2. Contact angles at Wenzel state: numerical versus analytical (w = 5 lu, s = 5 lu, h = 5 lu).
Static Contact Angle from [25]Contact Angle by SimulationContact Angle by Wenzel Equation (1)Error (%)
5329.729.930.77
6552.252.50.57
7669.369.60.43
9089.3900.78
104110.27110.380.1
115127.21127.490.22
127149.7150.10.27
Table 3. Numerical vs. analytical contact angles at Cassie–Baxter state (w = 10 lu, s = 5 lu, h = 10 lu).
Table 3. Numerical vs. analytical contact angles at Cassie–Baxter state (w = 10 lu, s = 5 lu, h = 10 lu).
Static Contact Angle from [25]Contact Angle by SimulationContact Angle by Cassie–Baxter
Equation (2)
Error (%)
53126.12126.830.56
65130.41130.120.22
76132.89133.590.52
90137.92138.590.48
104143.96144.140.12
115148.37148.830.31
127154.14154.220.05
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ashirbekov, A.; Zhumatay, N.; Kuljabekov, A.; Kabdenova, B.; Monaco, E.; Wang, L.; Rojas-Solórzano, L.R. Lattice Boltzmann Modeling of a Sessile and a Body Force-Driven Sliding Droplet over a Grooved Surface. Processes 2022, 10, 2356. https://doi.org/10.3390/pr10112356

AMA Style

Ashirbekov A, Zhumatay N, Kuljabekov A, Kabdenova B, Monaco E, Wang L, Rojas-Solórzano LR. Lattice Boltzmann Modeling of a Sessile and a Body Force-Driven Sliding Droplet over a Grooved Surface. Processes. 2022; 10(11):2356. https://doi.org/10.3390/pr10112356

Chicago/Turabian Style

Ashirbekov, Assetbek, Nursultan Zhumatay, Alibek Kuljabekov, Bagdagul Kabdenova, Ernesto Monaco, Lei Wang, and Luis R. Rojas-Solórzano. 2022. "Lattice Boltzmann Modeling of a Sessile and a Body Force-Driven Sliding Droplet over a Grooved Surface" Processes 10, no. 11: 2356. https://doi.org/10.3390/pr10112356

APA Style

Ashirbekov, A., Zhumatay, N., Kuljabekov, A., Kabdenova, B., Monaco, E., Wang, L., & Rojas-Solórzano, L. R. (2022). Lattice Boltzmann Modeling of a Sessile and a Body Force-Driven Sliding Droplet over a Grooved Surface. Processes, 10(11), 2356. https://doi.org/10.3390/pr10112356

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop