Abstract
The spreading kinetics of Herschel-Bulkley fluids on horizontal solid substrates were theoretically studied. The equations of film thickness were derived in both gravitational and capillary regimes. The dynamic contact angle for the capillary regime was also derived. Finally, a limiting result for the case of τ0 = 0 was obtained, which was compared with the known solution for validation. The results show that the yield behavior of the fluids had a significant impact on the spreading kinetics in both cases. Only when stress was larger than the yield stress, would substantial flow occur. The spreading zone was divided into two parts by the yield surface: sheared zone and yield zone, which was completely different from common Newtonian fluids or power-law fluids. The thickness of the yield zone mainly depended on yield stress and pressure gradient along the z-direction. According to the final evolution, both the film thickness and dynamic contact angle were affected not only by the power-law index but also by the yield behavior.
Introduction
Many studies have addressed the spreading of Newtonian and power-law fluids over solid surfaces by experiments, simulations, and theories [1–7]. For Newtonian fluids, Tanner law is considered to be the classic theory [3]. For power-law fluids, Starov built a hydrodynamic model for the spreading of a complete wetting case [4]. Wang experimentally investigated both complete wetting and partial wetting cases of spreading, which agreed with the theory proposed by Starov [5-6]. Carré and Eustach proposed a simplified theoretical analysis and experimentally observed the drop shape [7].
Besides power-law fluids, Bingham and Herschel-Bulkley fluids also belong to the common non-Newtonian type, which present yielding behavior. For Bingham fluids, the spreading dynamics both in impact and spin-coating cases have been analyzed in detail. In the impact case, the results indicated that the film shape mainly depends on initial yield stress, power-law index, and the impact force [8–13]. In the spin-coating case, the analysis showed that film thickness is not homogeneous, the yield behavior, power-law index, and the centrifugal force are the key factors on spreading dynamics [14–17].
Because of the complex rheological model of Herschel-Bulkley fluids, a minority of theories of complete wetting have been proposed. In this paper, we will theoretically analyze the spreading of Heschel–Bulkley fluids in a complete wetting case based on some proper assumptions. We begins with the simple spreading schematic of Newtonian fluids and the Navier–Stokes equation. According to the analytical approach of the spreading of power-law fluids, the film thickness equations are derived in both gravitational and capillary regimes and the dynamic contact angle is also discussed.
Assumptions
To simplify the problem, the following assumptions are made [18, 19]. 1) The spreading process belongs to the complete wetting case. 2) The fluids are incompressible, which will ensure the Navier-Stokes equation is workable. 3) The gas viscosity is ignored. 4) The film is much thinner than the horizontal length, so the flow can be converted into 2D. 5) The Reynolds number is small enough to ignore inertial influence. And 6) complete wetting is applied here to ensure a small contact angle.
Theory
The Constitutive Model of Herschel-Bulkley Fluids
The constitutive model of Herschel-Bulkley fluids is similar to Bingham fluids or power-law fluids. The specific constitutive equation is given as follows:where τ is shear stress, k is viscosity coefficient,is the shear rate, n is the power-law index, and τ0 is the initial yield stress.
Thin Film Equation
The spreading schematic of the advancing system of Newtonian or power-law fluids is shown in Figure 1. The derivation of the Herschel-Bulkley case was also started in Figure 1. In Figure 1, u is the horizontal velocity of the fluid, which is affected by time t and position x, z. U is the framing moving velocity of the contact line, θ(U) is the dynamic contact angle, and h is the height of film thickness of any position x.
FIGURE 1
With the proposed assumption (3), the pressure along the z-direction will meet the following equation.
Applying the Young-Laplace equation to the free surface (z = h), there iswhere pG represents the atmospheric pressure and σ is surface tension.
Integrating Eq. (2) with respect to the boundary condition listed in Eq. (3), the pressure can be further calculated as
When z = 0, the stress on the horizontal substrate can be obtained from Eq. (4) as
Consequently, when the initial yield stress τ0 is larger than τs, substantial flow will not occur. Thus the approximate condition for substantial flow is
Substantial flow will be discussed below. Owing to the initial yield stress, there must be a yield surface at z = h0. Thus the spreading zone will be divided into two zones, the significant shear exists only below the yield surface. The shear rate is zero at or above the yield surface. The viewpoint is consistent with the spreading of Bingham fluids proposed by Liu [20]. Therefore, the actual flow schematic of the advancing systems of Herschel-Bulkley fluids is shown in Figure 2. The dashed line separates the zone into sheared and yield zones respectively. The schematic figure is completely different from the Newtonian or power-law case.
FIGURE 2
For the Navier–Stokes equation, it can be written as where p is the pressure and x is the distance of one point from the z-axis. According to Eq. (1), when the stress is larger than the yield initial stress, the viscosity can be expressed as
Substituting Eq. (8) into Eq. (7), there is
Based on the boundary condition of no shear at the yield surface (h(x) = h0, ∂u/∂z = 0), integrating Eq. (9) with respect to z and the following equation can be acquired.With the boundary condition of no-slip at the solid surface (z = 0, u = 0), the expression of u can be obtained by integrating Eq. (10).
Equation (11) shows that the maximum velocity up occurs at the yield surface z = h0 that can be obtained as
When z > h0, the velocity is still equal to the maximum velocity up. According to mass conservation, there iswhere q is the flow of fluids.
For Herschel-Bulkley fluids, q is made up of two parts as
The first term of the right side represents the shear zone below the yield surface and its detailed expression is
Then the flow q will be acquired as follows:
In Fig. 2, the velocity along the x-direction is parabolic, which can also be explained by Eqs. (11)–(12). When it refers to the velocity, the case here is similar to classical Poiseuille flow. Both of them are driven by the pressure gradient, and they have the same initial and boundary conditions, thus the relationship between h and h0 can be obtained based on our previous work [21–22] as
Substituting Eq. (17) into Eq. (16), a new formula of flow q is generated as
Finally, substituting Eq. (18) into Eq. (13), the following equation is obtained as
According to Eq. (16), the pressure mainly depends on the gravity and capillary forces. In the following part, two limiting regimes will be considered further. One is only the gravitational regime taken into consideration by ignoring the capillary force, ρgR2>>σ. The other is only the capillary force regime considered, ignoring the gravity, σ>>ρgR2.
Gravitational Spreading Regime
The capillary force is ignored here, so Eq. (19) can be simplified into the following equation.
When the contact line moves at velocity U, the above equation can be converted into the following form by introducing a new variable .
Integrating the above equation with respect to ξ, there isWith the boundary condition: ξ = 0, h = 0, Eq. (22) can be solved as
Equation (23) addresses the fact that film thickness in capillary spreading is determined by initial yield stress τ0 and power index n.
Capillary Spreading Regime
In the condition, the gravitational action is ignored, then Eq. (19) can be simplified as follows
Introducing the above variable ξ and h0 into Eq. (24),
Integrating the above equation with respect to ξ,
Then it is supposed that , so Eq. (26) can be translated intoWith the boundary condition: ξ = 0, h = 0, Eq. (27) can be solved as
Further, the film thickness equation can be obtained by replacing h0,
Based on Eq. (17), the variable λ can be calculated as
In this way, the equation of film thickness h can be obtained, which depends on the power index n and initial yield stress τ0.
Dynamic Contact Angle
The inclination angle at x = ξ can be calculated by differentiating film thickness h as shown in Eq. (29) for the capillary spreading regime.
In previous studies, many researchers found that the local microscopic contact angle can not be measured directly for boundary conditions. Thus dynamic contact angle θ is taken as the replacement, which is equal to the inclination angle at x = xm. Furthermore, the dynamic contact angle can be obtained as
To obtain Eq. (32), it is considered that when the angle θ is small enough, it is approximately equal to sin θ or tan θ.
Example
In a previous study [21, 22], the typical Herschel-Bulkley model was obtained as follows
The rheological figure is shown in Figure 3. The model was taken as an example to explain the spreading results.
FIGURE 3
Figure 4 gives the film thickness in the gravitational regime. The area below the dashed line is the sheared zone. In Figure 5, the figure of the contact angle vs. the moving velocity at different film thicknesses is presented. It can be found that a thinner film will lead to a smaller contact angle.
FIGURE 4
FIGURE 5
Discussion and Conclusion
The spreading of Herschel-Bulkley fluids over solid substrates for a complete wetting case was theoretically studied here. The film thicknesses in both the gravitational and capillary regimes were derived and the dynamic contact angles were also explored. There was an obvious difference between other kinds of fluids (Newtonian fluids and power-law fluids). The flow was dramatically affected by the yield behavior of Herschel-Bulkley fluids. According to the analysis, only when the stress on the substrate was larger than the initial yield stress τ0, would the flow be appreciable. Even if appreciable flow occurred, the spreading zone was also divided into two zones: one was the shear zone and the other was the yield zone. The dividing surface was named yield surface at z = h0, which does not exist in common fluids such as Newtonian fluids or power-law fluids. In “Example” section, the detailed numerical examples are given to explain the results further.
It is noteworthy that when initial yield stress τ0 = 0, the fluids become power-law fluids. The film thicknesses of the gravitational regime and the capillary regime will be obtained as follows
Equations (34) and (35) represent the gravitational case and the capillary case, respectively. Also, the dynamic contact angle θ of the capillary regime will be
Substituting Eq. (35) into Eq. (36), then the specific equation of the dynamic contact angle can be obtained as
It can be found that all of them are in agreement with the results obtained by Wang [5]. Another noticeable feature of the analysis is that some of the equations have a collective limitation , which means the theoretical results are not applicable for Bingham fluids, although the fluid corresponds to it when n = 1.
Funding
This work was financially supported by the Young and Middle-Aged Scientific Research Backbone Training Project of Nantong Institute of Technology (Grant No. ZQNGG204), the Key University Science Research Project of Jiangsu Province (Grant No. 18KJA460006), the Priority Discipline Construction Program of Jiangsu Province (Grant No. 2016-9), the Top-Notch Academic Programs Project of Jiangsu Higher Education Institutions (Grant No. 2020-9), and the Key R&D plan of Jiangsu Province (Grant No. BE2018010-4).
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
JZ: idea, HG: method, JS: calculation, BL: calculation, JZ: calculation, WW: validation and writing.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2020.609926/full#supplementary-material.
Glossary
- g
Gravitational acceleration
- h
The height of film thickness at position x
- h0
The film thickness of yield surface
- k
Viscosity coefficient
- n
The power-law index
- p
Pressure
- pG
Atmospheric pressure
- q
Flow of fluids
- t
Time
- u
Velocity along the x-direction
- up
Maximum velocity of u
- U
Frame moving velocity
- x/z
Values at two directions in the Cartesian coordinate system
Shear rate
- θ
Dynamic contact angle
- λ
Variable
- ξ
Variable, ξ = x-Ut
- ρ
Density
- σ
Surface tension
- τ
Stress
- τ0
Initial yield stress
- τs
Stress on the substrate
Nomenclature Symbols
References
1.
LiangZPWangXDLeeDJ. Spreading dynamics of power-law fluid droplets. Phys-Condens Mat (2009) 21:3543–9. 10.1088/0953-8984/21/46/464117
2.
NeogiPYbarraRM. The absence of a rheological effect on the spreading of small drops. J Chem Phys (2001) 115:7811–3. 10.1063/1.1415455
3.
TannerLH. The spreading of silicone oil drops on horizontal surfaces. J Phys D Appl Phys (1979) 12:1473. 10.1088/0022-3727/12/9/009
4.
StarovVMTyatyushkinANVelardeMG. Spreading of non-Newtonian liquids over solid substrates. J Colloid Interface Sci (2003) 257:284–90. 10.1016/s0021-9797(02)00034-6
5.
ZhdanovXDLeeDJPengXF. Spreading dynamics and dynamic contact angle of non-Newtonian fluids. Langmuir (2007) 23:8042–7. 10.1021/la0701125
6.
LaiXDZhangYLeeDJ. Spreading of completely wetting or partially wetting power-law fluid on solid surface. Langmuir (2007) 23:9258–62. 10.1021/la700232y
7.
PengAEustacheF. Spreading kinetics of shear-thinning fluids in wetting and dewetting modes. Langmuir (2000) 16:2936–41. 10.1021/la991021d
8.
GermanGBertolaV. The spreading behaviour of capillary driven yield-stress drops. Colloid Surf Physicochem Eng Aspect (2010) 366:18–26. 10.1016/j.colsurfa.2010.05.019
9.
GermanGBertolaV. Impact of shear-thinning and yield-stress drops on solid substrates. J Phys Condens Matter (2009) 21:375111. 10.1088/0953-8984/21/37/375111
10.
NigenS. Experimental investigation of the impact of an (apparent) yield-stress material. Atomization Sprays (2005) 15:103–18. 10.1615/atomizspr.v15.i1.60
11.
SaïdiAMartinCMagninA. Influence of yield stress on the fluid droplet impact control. J Non-Newtonian Fluid Mech (2010) 165:596–606. 10.1016/j.jnnfm.2010.02.020
12.
LuuL-HForterreY. Drop impact of yield-stress fluids. J Fluid Mech (2009) 632:301–27. 10.1017/s0022112009007198
13.
KimEBaekJ. Numerical study of the parameters governing the impact dynamics of yield-stress fluid droplets on a solid surface. J Non-Newtonian Fluid Mech (2012) 173–174:62–71. 10.1016/j.jnnfm.2012.02.005
14.
BurgessSLWilsonSDR. Spin‐coating of a viscoplastic material. Phys Fluids (1996) 8(9):2291–7. 10.1063/1.869016
15.
JenekheSASchuldtSB. Flow and film thickness of Bingham plastic liquids on a rotating disk. Chem Eng Commun (1985) 33:135–47.
16.
TabuteauHBaudezJCChateauX. Flow of a yield stress fluid over a rotating surface. Rheol Acta (2007) 46:341–55. 10.1080/00986448508911165
17.
TsamopoulosJAChenMEBorkarAV. On the spin coating of viscoplastic fluids. Rheol Acta (1996) 35:597–615. 10.1007/bf00396510
18.
FermigierMJenfferP. An experimental investigation of the dynamic contact angle in liquid–liquid systems. J Colloid Interface Sci (1991) 146:226–41. 10.1016/0021-9797(91)90020-9
19.
ChenJ-DWadaN. Edge profiles and dynamic contact angles of a spreading drop. J Colloid Interface Sci (1992) 148:207–22. 10.1016/0021-9797(92)90129-a
20.
LiuKFMeiCC. Approximate equations for the slow spreading of a thin sheet of Bingham plastic fluid. Phys Fluids (1998) 2(1):30–6.
21.
WuWHuangXFangC. An improved MRT-LBM for Herschel-Bulkley fluids with high Reynolds number, Numer Heat Tr B-Fund (2017) 72(6):1–12. 10.1080/10407790.2017.1409521
22.
WuWHuangXLiYFangauCJiangX. A modified LBM for non-Newtonian effect of cement paste flow in 3D printing. Rpj (2019) 25(1):22–9. 10.1108/rpj-06-2017-0124
Summary
Keywords
spreading, Herschel-Bulkley, film thickness, dynamic contact angle, non-newtonian
Citation
Zhang J, Gu H, Sun J, Li B, Jiang J and Wu W (2020) Spreading Kinetics of Herschel-Bulkley Fluids Over Solid Substrates. Front. Phys. 8:609926. doi: 10.3389/fphy.2020.609926
Received
24 September 2020
Accepted
17 November 2020
Published
11 December 2020
Volume
8 - 2020
Edited by
Nuno A. M. Araújo, University of Lisbon, Portugal
Reviewed by
Shahin Rouhani, Sharif University of Technology, Iran
Martin Kröger, ETH Zürich, Switzerland
Updates
Copyright
© 2020 Zhang, Gu, Sun, Li, Jiang and Wu.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jie Zhang, zhangjienit@163.com
This article was submitted to Mathematical and Statistical Physics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.