Hostname: page-component-586b7cd67f-dlnhk Total loading time: 0 Render date: 2024-11-26T13:14:34.236Z Has data issue: false hasContentIssue false

Experimental observation of a confined bubble moving in shear-thinning fluids

Published online by Cambridge University Press:  06 December 2022

SungGyu Chun
Affiliation:
Department of Mechanical Engineering and Science, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
Bingqiang Ji
Affiliation:
Department of Mechanical Engineering and Science, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
Zhengyu Yang
Affiliation:
Department of Mechanical Engineering and Science, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
Vinit Kumar Malik
Affiliation:
Department of Mechanical Engineering and Science, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
Jie Feng*
Affiliation:
Department of Mechanical Engineering and Science, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA Materials Research Laboratory, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
*
Email address for correspondence: [email protected]

Abstract

The motion of a long gas bubble in a confined capillary tube is ubiquitous in a wide range of engineering and biological applications. While the understanding of the deposited thin viscous film near the tube wall in Newtonian fluids is well developed, the deposition dynamics in commonly encountered non-Newtonian fluids remains much less studied. Here, we investigate the dynamics of a confined bubble moving in shear-thinning fluids with systematic experiments, varying the zero-shear-rate capillary number $Ca_0$ in the range of $O(10^{-3}\unicode{x2013}10^2)$ considering the zero-shear-rate viscosity. The thickness of the deposited liquid film, the bubble speed and the bubble front/rear menisci are measured, which are further rationalized with the recent theoretical studies based on appropriate rheological models. Compared with Newtonian fluids, the film thickness decreases for both the carboxymethyl cellulose and Carbopol solutions when the shear-thinning effect dominates. We show that the film thickness follows the scaling law from Aussillous & Quéré (Phys. Fluids, vol. 12, no. 10, 2000, pp. 2367–2371) with an effective capillary number $Ca_e$, considering the characteristic shear rate in the film as proposed by Picchi et al. (J. Fluid Mech., vol. 918, no. A7, 2021, pp. 1–30). $Ca_e$ is calculated by the Carreau number and the power-law index from the Carreau–Yasuda rheological model. The shear-thinning effect also influences the bubble speed and delays the transition to the parabolic region in the bubble front and rear menisci. In particular, a high degree of undulations on the bubble surface results in an intricate rear viscosity distribution for the rear meniscus and the deviation between the experiments and theory may require a further investigation to resolve the axial velocity field. Our study may advance the fundamental understandings and engineering guidelines for coating processes involving thin-film flows and non-Newtonian fluids.

Type
JFM Papers
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
© The Author(s), 2022. Published by Cambridge University Press

1. Introduction

The transport of long gas bubbles or liquid drops in confined geometries plays an important role in many engineering and biological settings, such as enhanced oil recovery (Tran et al. Reference Tran, Ahmad, Neogi and Bai2016; Grassia Reference Grassia2019; Majeed et al. Reference Majeed, Kamal, Zhou and Solling2021), coating processes (Yu, Khodaparast & Stone Reference Yu, Khodaparast and Stone2017; Jeong et al. Reference Jeong, Kvasnickova, Boutin, Cébron and Sauret2020), drug delivery (Hernot & Klibanov Reference Hernot and Klibanov2008; Gao et al. Reference Gao, Chan, Gu, Lin, Zhang, Yeo, Alsema, Arora, Chong and Shi2016), biomechanics and biomedical devices (Clanet, Héraud & Searby Reference Clanet, Héraud and Searby2004; Chao, Jin & Fan Reference Chao, Jin and Fan2020; Ma et al. Reference Ma, Sun, Duan, van den Berg, Eijkel and Xie2020; Li et al. Reference Li, Li, Li, Chen, Li and Chen2021). When such a long bubble of length $L \gg R$ translates at a constant speed $U$ in a circular capillary of radius $R$, the bubble forms a symmetrical bullet shape, commonly called a Taylor bubble, and a thin film of liquid is generated between the bubble and capillary. Quantifying the deposition of this liquid film and its relationship with the bubble speed, fluid properties and channel geometries provides crucial information for determining the mass, momentum and heat transport in a wide range of multi-phase flow scenarios, and thus has remained as a research focus for decades.

The deposition of Newtonian fluids has been extensively studied regarding various geometries and fluid properties. Pioneering investigations on this topic were conducted by Bretherton (Reference Bretherton1961) and Taylor (Reference Taylor1961). For a long bubble translating in confined geometries with small dimensions where gravity plays a negligible role, the dynamics is characterized by the interplay between the viscosity and surface tension, as captured by the definition of the capillary number, $Ca = {\mu }U/\sigma$, where $\mu$ is the fluid viscosity and $\sigma$ is the surface tension (Aussillous & Quéré Reference Aussillous and Quéré2000; Jeong et al. Reference Jeong, Kvasnickova, Boutin, Cébron and Sauret2020). For $Ca \ll 1$, Bretherton (Reference Bretherton1961) found that the thickness of the thin liquid film $h$ scales as $h/R \sim Ca^{2/3}$ in regimes where inertia effects are negligible compared to surface tension and viscous effects. This relation was later extended to cover the range of $Ca < 2$ in the scaling analysis of Aussillous & Quéré (Reference Aussillous and Quéré2000), where the radius of curvature of the static meniscus was accounted as $R - h$ rather than $R$, yielding a semi-empirical equation as

(1.1)\begin{equation} \frac{h}{R}=\frac{1.34Ca^{2/3}}{1+2.5\times1.34Ca^{2/3}}. \end{equation}

So far, the fluid deposition by a confined bubble in non-Newtonian fluids has been much less understood, although a lot of working fluids, such as polymer solutions, colloidal suspensions and biologically relevant fluids, show non-Newtonian behaviours in many practical applications (Abishek, King & Narayanaswamy Reference Abishek, King and Narayanaswamy2015; Moreira et al. Reference Moreira, Rocha, Carneiro, Araújo, Campos and Miranda2020; Li et al. Reference Li, Li, Li, Chen, Li and Chen2021; Zhao et al. Reference Zhao, Ma, Liu, Yao, Yang and Chen2021). Resolving the hydrodynamics of a confined bubble in non-Newtonian fluids presents more complexities than in Newtonian cases because of the spatial and temporal changes of shear stress and the corresponding variations of the rheological properties. Table 1 lists a summary of the experimental, numerical and theoretical studies about the liquid film deposition dynamics in non-Newtonian fluids. In particular, many prior studies assumed a simple power-law model with a stress/shear rate relationship as $\tau ={\kappa }\dot {\gamma }^{n_p}$ (where $n_p$ is the power-law index, $\dot {\gamma }$ the shear rate and $\kappa$ the consistency factor in Pa s$^{n_p}$) to represent the shear-thinning/thickening fluids (Kamişli & Ryan Reference Kamişli and Ryan2001; de Sousa et al. Reference de Sousa, Soares, de Queiroz and Thompson2007), and suggested that the liquid film thickness scales as (Gutfinger & Tallmadge Reference Gutfinger and Tallmadge1965; Hewson, Kapur & Gaskell Reference Hewson, Kapur and Gaskell2009)

(1.2)\begin{equation} \frac{h}{R} \sim \widehat{Ca}^{2/(2n_p+1)} \quad \text{with}\ \widehat{Ca} = \kappa(U/R)^{n_p}/(\sigma/R), \end{equation}

where $\widehat {Ca}$ is a modified capillary number. A numerical prefactor on the right-hand side of (1.2) is suppressed here, but its value would depend on $n_p$ (Gutfinger & Tallmadge Reference Gutfinger and Tallmadge1965; Hewson et al. Reference Hewson, Kapur and Gaskell2009). Nevertheless, the power-law model fails to reproduce the low-shear-rate viscosity plateau and has a well-known singularity at zero-shear-rate, leading to an inaccurate velocity profile in multiple free-surface flow scenarios (Bird, Armstrong & Hassager Reference Bird, Armstrong and Hassager1987; Myers Reference Myers2005; Hewson et al. Reference Hewson, Kapur and Gaskell2009; Picchi et al. Reference Picchi, Poesio, Ullmann and Brauner2017).

Table 1. Chronological selection of previous experimental, numerical and theoretical studies on the liquid film thickness of a long bubble translating through non-Newtonian fluids in confined geometries. The non-dimensional groups appearing above are defined as follows: $Ca_0 = {\mu _0}U/\sigma$, $Wi = {\lambda }U/H$, $\widehat {Ca}_{HB} = k_{HB}(U/R)^{n_{HB}}/(\sigma /R)$ and $Ca_e = {\mu _e}U/\sigma$. Here, $\mu _0$ is the zero-shear-rate viscosity, $\lambda$ is the relaxation time, $H$ is the half of the gap in the Hele–Shaw cell, $k_{HB}$ and $n_{HB}$ are the consistency factor and the power-law index of the Herschel–Bulkely model, respectively, and $\mu _e$ is the effective viscosity defined as $\mu _e=\mu (\dot {\gamma }={U}/{h})$. Ro & Homsy (Reference Ro and Homsy1995) used $m$ and $k$ to represent the degree of shear and normal stress thinning, respectively, $S$ is the ratio of the solvent viscosity to the sum of the polymer and solvent viscosity, and $\delta$ is the ratio between $\lambda$ and the characteristic residence time in the gap. Laborie et al. (Reference Laborie, Rouyer, Angelescu and Lorenceau2017) used $a$ and $b$ are the fitting parameters and $B$ is the non-dimensional number comparing the yield stress to the capillary pressure. Zhao et al. (Reference Zhao, Ma, Liu, Yao, Yang and Chen2021) used $W$ to be the width of the rectangular microchannel.

Therefore, recent studies have focused on developing more generalized scaling laws to overcome the limitation of the power-law model. Using the Carbopol solution (1.1 wt%) which shows a strong yield-stress effect, Laborie et al. (Reference Laborie, Rouyer, Angelescu and Lorenceau2017) experimentally investigated the deposition of the yield-stress fluid in circular channels and developed a semi-empirical scaling law for the film thickness, considering the competition of the yield stress, the capillary pressure and the viscous stress (see table 1). Additionally, a recent study by Picchi et al. (Reference Picchi, Ullmann, Brauner and Poesio2021) investigated the motion of a Taylor bubble through an Ellis fluid and identified a scaling law of the film thickness with the generalized effective viscosity defined by the characteristic shear rate in the liquid film. Several recent studies also showed that the predictions using the power-law model can be erroneous at flow settings with a low-shear-rate region and have a small range of applicability compared to the more accurate Ellis and Carreau–Yasuda (C-Y) models (Moukhtari & Lecampion Reference Moukhtari and Lecampion2018; Boyko & Stone Reference Boyko and Stone2021; Picchi et al. Reference Picchi, Ullmann, Brauner and Poesio2021). However, aforementioned studies are still limited to mostly numerical or theoretical perspectives, and a systematic experimental verification of the theoretical predictions is still lacking.

In this work, we report an experimental investigation for the effect of shear-thinning rheology on the film deposition dynamics, bubble speed and bubble shape variations when a bubble is moving in a circular capillary tube filled with non-Newtonian fluids. We consider the range of the zero-shear-rate capillary number, $Ca_0 = {\mu _0}U/\sigma$ where $\mu _0$ is the zero-shear-rate viscosity for non-Newtonian fluids, over six orders of magnitude ($7\times 10^{-3} < Ca_0 < 830$). In § 2, we provide the experimental framework to measure the film thickness and bubble shape. The generalized C-Y model is used to describe the full range of rheological properties for carboxymethyl cellulose and Carbopol solutions (§ 3.1). Scaling laws for the film thickness based on different rheological models are compared (§§ 3.2 and 3.3), showing the film thickness and the bubble speed (§ 3.4) scale with an effective capillary number, $Ca_e$, containing two dimensionless numbers that describe the fluid rheology. Finally, the bubble front and rear menisci are experimentally characterized and further compared with the lubrication theory considered by Picchi et al. (Reference Picchi, Ullmann, Brauner and Poesio2021) (§ 3.5).

2. Material and experimental set-up

We use carboxymethyl cellulose (CMCell, Sigma Aldrich) and Carbopol (981, Lubrizol) solutions with different mass fractions as non-Newtonian fluids, and pure glycerin ($\mu = 0.87$ Pa s, $\sigma = 63.4\ {\rm mN}\ {\rm m}^{-1}$, Fisher Scientific) as a Newtonian fluid for a baseline comparison. CMCell consists of bentonite as a major component (Benchabane & Bekkour Reference Benchabane and Bekkour2008), while Carbopol is composed of polyacrylic acid resins (Laborie et al. Reference Laborie, Rouyer, Angelescu and Lorenceau2017). We prepare the CMCell solutions with four different mass fractions (0.5, 1.0, 1.5 and 2.0 wt%), and the Carbopol solutions with three different mass fractions (0.1, 0.2 and 0.5 wt%, neutralized using 1 M sodium hydroxide). No elastic behaviours are expected in such a low CMCell mass fraction range given the shear rates in the current experiments (Ghannam & Esmail Reference Ghannam and Esmail1997; Benchabane & Bekkour Reference Benchabane and Bekkour2008). All aqueous solutions were prepared by gradually dissolving a known weight of powders into deionized water in a cylindrical beaker with continuous stirring. The mixing was maintained for 24 h until clear and homogeneous solutions were produced. The surface tension of the working fluid is determined by the pendant drop method (Rotenberg, Boruvka & Neumann Reference Rotenberg, Boruvka and Neumann1983; Song & Springer Reference Song and Springer1996).

2.1. Experimental

Bubble motion experiments were performed in a circular glass capillary tube with a length of 500 mm and an internal radius $R = 0.47$ mm. The glass capillary was held vertical and the central length of the capillary was contained in a clear rectangular box filled with the working fluid to decrease refractive index difference compared to the glass as well as the optical distortion from the curvature of the tube wall (Yu et al. Reference Yu, Khodaparast and Stone2017; Zhao et al. Reference Zhao, Pahlavan, Cueto-Felgueroso and Juanes2018).

During the experiments, the inlet of the glass capillary tube was connected to a syringe pump (11 Pico Plus Elite, Harvard Apparatus) using a flexible connecting tube. After the glass capillary was pre-filled with a sample solution, a small volume of air was created in the connecting tube. The flow rate was then set to a very small value $\approx$10 $\mathrm {\mu }$l min$^{-1}$ to steadily transfer the long air bubble (of length $L \gg R$) from a flexible connecting tube into the inlet of the glass capillary. When the front of the bubble reached the inlet, the syringe pump was set to the targeted flow rate accordingly. The optical images of the region of interest (ROI) were recorded at the rate of 60 frames per second using a digital camera (20.9 Megapixel, D7500, Nikon) equipped with a long working distance objective lens (12$\times$ zoom lens system, Navitar); see the schematic in figure 1(a). The maximum resolution of the image in our experimental configuration is $\approx$1.3 $\mathrm {\mu }$m per pixel, leading to a maximum relative error in the film thickness measurement of less than 12 %. The average velocity of the bubble can be evaluated by tracking the gas–liquid interface at the bubble front tip using the images taken with a lower magnification. An analytical balance (ME104E, Mettler Toledo) was installed at the outlet of the glass capillary to confirm the flow rate. In the current experiment, $Bo={\rho }gR^2/\sigma < 0.04$ and $Re = {\rho }UR/\mu \ll 1$, so gravity and inertia effects are negligible (Atasi et al. Reference Atasi, Khodaparast, Scheid and Stone2017; Magnini et al. Reference Magnini, Khodaparast, Matar, Stone and Thome2019).

Figure 1. (a) Schematic of the experimental configuration. A cylindrical glass tube (with an inner radius of 0.47 mm) is filled with a sample solution (e.g. glycerin, carboxymethyl cellulose (CMCell; 0.5, 1.0, 1.5 and 2.0 wt$\%$), and Carbopol (0.1, 0.2 and 0.5 wt$\%$)). The central part of the circular glass capillary is submerged in a bath of a sample solution to match the refractive index of glass. Inset: schematic of a translating air bubble confined in a circular tube. (b) Typical experimental images of a long bubble as it translates in a circular capillary filled with glycerin (left) at $Ca_0 = 7.48 \times 10^{-2}$, CMCell (1.0 wt$\%$; centre) at $Ca_0 = 3.53 \times 10^{-2}$ and Carbopol (0.2 wt$\%$; right) solutions at $Ca_0 = 3.53 \times 10^{-2}$. Images contain the front meniscus and the middle part of the bubble, where the film thickness is uniform. The scale bar is 0.5 mm. (c) Schematic for the mass balance analysis regarding the deposition of a liquid film in a circular capillary tube. (d) Comparison of the liquid film thickness obtained by the image visualization, $h_i$ and the mass balance analysis of the liquid plug, $h_m$ (2.1) for experimental cases over the range of $8\times 10^{-3} < Ca_0 < 8\times 10^{2}$. All the data lie along the solid line with a slope of a unity.

2.2. Measurement of the liquid film thickness

The liquid film thickness around the bubble is estimated with the image visualization, in which the film profile around the bubble is obtained from the difference between the position of the bubble surface and the tube wall (see figure 1b). Since the length of the bubble is larger than the size for the field of view in the images, we use a time-strip analysis with ImageJ to ensure an accurate measurement of the uniform liquid film around the bubble.

To confirm the accuracy from the image visualization method, we also measure the film thickness using mass balance analysis. In the mass balance analysis, the film thickness is determined based on the change of the length of the liquid plug $L_p$. We ensure the initial length of the liquid plug is $L_p\approx 7$ cm, so that $L_p \gg R$. The plug advances inside the tube, deceasing $L_p$ (see figure 1c) due to the deposition of the film on the tube wall. The moving positions of the front and rear menisci of the liquid plug are analysed using ImageJ, which determines the velocity at the front and rear menisci of the liquid plug, $U_{pf}$ and $U_{pr}$, respectively, ranging from 0.2 to 50 mm s$^{-1}$. Assuming a homogeneous deposition of the liquid film near the cylindrical capillary tube wall, a mass balance on the moving plug of length $L_P$ yields the relation as (Laborie et al. Reference Laborie, Rouyer, Angelescu and Lorenceau2017)

(2.1)\begin{equation} \frac{h}{R}=1-\sqrt{1+\frac{1}{U_{pr}}\frac{\mathrm{d} L_p}{\mathrm{d} t}}=1-\sqrt{\frac{U_{pf}}{U_{pr}}}. \end{equation}

Figure 1(d) compares the liquid film thickness measured by the image visualization, $h_i$, and the mass balance analysis, $h_m$, when a bubble is translating in the glycerin, 1.0 wt$\%$ CMCell and 0.2 wt$\%$ Carbopol solutions. The measurement is conducted over a wide range of $Ca_0$ ($8\times 10^{-3} < Ca_0 < 8\times 10^{2}$). The results show that all the experimental measurement data lie on the line with a slope of unity, confirming the two methods give the same result on each experiment. We note that the deposited film thickness is observed to be invariant by rotation along the axis of the glass capillary since the shear-thinning effect dominates for the Carbopol solution with low mass fractions (figure 1b), while the experiments performed by Laborie et al. (Reference Laborie, Rouyer, Angelescu and Lorenceau2017) showed that the annular Carbopol solution film thickness was non-uniform azimuthally when the yield stress is important. In the following, we will use the image visualization to determine the film thickness, $h$.

3. Results and discussion

3.1. Rheological properties of CMCell and Carbopol solutions

Rheological measurements of the CMCell and Carbopol solutions are performed with a controlled stress rheometer (DHR-3, TA Instrument) using a parallel-plate geometry (with a diameter of 25 mm) at a controlled temperature of 25 $^{\circ }$C. Under a simple shear, the rheological properties of shear-thinning fluids are classically modelled by the power-law model $\mu =\kappa \dot {\gamma }^{n_{p}-1}$ (Kamişli & Ryan Reference Kamişli and Ryan2001; de Sousa et al. Reference de Sousa, Soares, de Queiroz and Thompson2007). However, the power-law model cannot predict the viscosity at the low-shear-rate region (Picchi et al. Reference Picchi, Poesio, Ullmann and Brauner2017, Reference Picchi, Ullmann, Brauner and Poesio2021) where the viscosity approaches to a constant value, known as the zero-shear-rate viscosity. Instead, the Ellis model (Reiner & Leaderman Reference Reiner and Leaderman1960) was proposed to capture such a viscosity plateau with a constitutive equation as

(3.1)\begin{equation} \frac{\mu}{\mu_0}=\frac{1}{1+(\tau/\tau_{1/2})^{\alpha-1}}, \end{equation}

where $\tau _{1/2}$ is the shear stress at which the viscosity is half of the Newtonian limit, while $\mu _0$ and $\alpha$ are the zero-shear-rate viscosity and the degree of shear-thinning, respectively.

Here, we consider the C-Y model, which has been used to describe emulsions, protein solutions and polymer melts (Myers Reference Myers2005; Picchi et al. Reference Picchi, Poesio, Ullmann and Brauner2017). The C-Y model is more convenient for experimental analysis since it expresses the viscosity as an explicit function of the shear rate (Carreau Reference Carreau1972; Pipe, Majmudar & McKinley Reference Pipe, Majmudar and McKinley2008; Morozov & Spagnolie Reference Morozov and Spagnolie2015). The constitutive equation of the C-Y model is $\mu =(\mu _0-\mu _\infty )(1+(\lambda \dot {\gamma })^a)^{(n_c-1)/a}+\mu _\infty$, where $\mu _0$, $\mu _\infty$, $n_c$ and $a$ are the zero-shear-rate viscosity, infinite-shear-rate viscosity, the power-law index and dimensionless parameter, respectively. Here, $\lambda$ is the inverse of a characteristic shear rate at which shear-thinning becomes apparent. Figures 2(a) and 2(b) demonstrate that the C-Y model can well capture the rheological behaviours for both the CMCell and Carbopol solutions in the range of shear rates over five orders of magnitude. We note that $\mu _\infty$ is neglected in the following discussion since $\mu _\infty /\mu _0< O(10^{-2})$ and therefore the constitutive equation for the C-Y model can be written as

(3.2)\begin{equation} \mu\approx\mu_0(1+(\lambda\dot{\gamma})^a)^{(n_c-1)/a}.\end{equation}

Figure 2. (a) Rheogram of the glycerin and carboxymethyl cellulose (CMCell) solutions with different mass fractions: viscosity $\mu$ versus shear rate $\dot {\gamma }$. (b) Rheogram of the Carbopol solutions with different mass fractions: $\mu$ versus $\dot {\gamma }$. The dashed lines represent a fitting with the Carreau–Yasuda (C-Y) model. (c) Dimensionless parameter $a$ versus the power-law index $n_c$ in the C-Y model or the degree of shear-thinning $\alpha$ in the Ellis model. (d) Ellis number $El$ versus Carreau number $Cu$ in the current experiments. Error bars are smaller than the symbols.

A very recent theoretical work by Picchi et al. (Reference Picchi, Ullmann, Brauner and Poesio2021) derived the film thickness as a function of the Ellis number $El$ and the degree of shear-thinning $\alpha$ for an Ellis fluid, where $El$ is the ratio between the characteristic shear rate of the fluid and the characteristic shear rate in the liquid film as

(3.3)\begin{equation} El=\frac{\tau_{1/2}h}{U\mu_0}. \end{equation}

To compare our experimental measurements with the lubrication theory of Picchi et al. (Reference Picchi, Ullmann, Brauner and Poesio2021), we rewrite the C-Y model with a similar form to the Ellis model as

(3.4)\begin{equation} \frac{\mu}{\mu_0}=(1+{({Cu\tilde{\dot{\gamma}}})}^a)^{(n_c-1)/a} \quad \text{with} \ Cu = \frac{\lambda{U}}{h} \ \text{and} \ \tilde{\dot{\gamma}}={\frac{\dot{\gamma}}{U/h}},\end{equation}

where the Carreau number $Cu$ is the ratio between the effective shear rate $U/h$ in the film and the cross-over strain rate $1/\lambda$ (Datt et al. Reference Datt, Zhu, Elfring and Pak2015), and $\tilde {\dot {\gamma }}$ is the dimensionless shear rate. With the rheological data in figure 2(a,b), we used the method of least squares fitting to obtain all the parameters for the Ellis and C-Y models. Notably, the value of $a$ in the C-Y model is chosen to impose $n_c=1/\alpha$ for an analogy between the Ellis and C-Y models, and we find $a\approx 0.45 n_c ^{-1}$, as shown in figure 2(c). Next, for all the experimental cases, once the bubble speed and the film thickness are determined, $El$ and $Cu$ are calculated as shown in figure 2(d). In the current work, $Cu$ is found to be inversely proportional to $El$ with an experimentally fitted relation of $Cu\approx 1.3El^{-1}$. The rheological parameters of the CMCell and Carbopol solutions are reported in table 2, and we will focus on the effects of $Cu$ and $n_c$ on the deposition dynamics of the working fluids.

Table 2. Rheological properties of the CMCell and Carbopol solutions with different mass fractions computed by the C-Y model.

We note that Carbopol solutions can exhibit both yield stress and shear-thinning behaviours. However, we use a low mass fraction of Carbopol to diminish the yield stress effect (Spiers, Subbaraman & Wilkinson Reference Spiers, Subbaraman and Wilkinson1975; Ma et al. Reference Ma, Bai, Chang, Yi, Jiao and Du2015), so that only the shear-thinning behaviour dominates. We further justify this consideration by fitting our rheological data of Carbopol with the Herschel–Bulkley model $\tau =\tau _y+k_{HB}\dot {\gamma }^{n_{HB}}$, where $\tau _y$ is the yield stress and $k_{HB}$ is the consistency factor. A dimensionless number $B=\tau _y R/\sigma$ is suggested to compare the yield stress to the capillary pressure (Deryagin & Levi Reference Deryagin and Levi1964; Laborie et al. Reference Laborie, Rouyer, Angelescu and Lorenceau2017). In our experiments, $B$ is $O (10^{-3})$ while $B$ is $O (1)$ in Laborie et al. (Reference Laborie, Rouyer, Angelescu and Lorenceau2017), thus we neglect the yield stress effect in the following discussion.

3.2. Scaling of the film thickness with the modified capillary number

The film thickness measured in the experiments with the CMCell and Carbopol solutions is shown in figures 3(a) and 3(b), respectively, as a function of the zero-shear-rate capillary number $Ca_0$, as well as the modified capillary number $\widehat {Ca}$. As shown in figure 3(a), compared to a Newtonian fluid at the same $Ca_0$, the bubble forms a thinner liquid film in both the CMCell and Carbopol solutions because of the shear-thinning effect. The thinner liquid film formed in the Carbopol solutions compared with that in the CMCell solutions at the same $Ca_0$ results from that the stronger shear-thinning effect in the Carbopol solutions than that in the CMCell solutions, as indicated by lower power-law indices of the Carbopol solutions compared with those of the CMCell solutions in table 2.

Figure 3. (a) Non-dimensional liquid film thickness h/R as a function of $Ca_0$. The black line represents prediction of (1.1). (b) $h/R$ versus $\widehat {Ca}^{2/{2n_p+1}}$, where $\widehat {Ca}=\kappa (U/R)^{n_p}/(\sigma /R)$, with $\kappa$ and $n_p$ ranging from 0.4 to 5.4 $\textrm {Pa s}^{n_p}$ and 0.16 to 0.51, respectively. The black line represents the prediction of (1.2) with a prefactor of 1. The experimental measurements are shown as open symbols, and error bars are smaller than the symbols.

With $\widehat {Ca}$ to compare the shear stress from the power-law model and the capillary pressure, (1.2) has been used to predict the film thickness (Gutfinger & Tallmadge Reference Gutfinger and Tallmadge1965; Hewson et al. Reference Hewson, Kapur and Gaskell2009). We note that the values of $\kappa$ and $n_p$ for $\widehat {Ca}$ are obtained by fitting the rheological data with a power-law model considering the range of shear rates exhibiting a shear-thinning behaviour, i.e. $\dot {\gamma } = O(10\unicode{x2013}10^3$ s$^{-1})$ for CMCell and $O(10^{-1}\unicode{x2013}10^2$ s$^{-1}$) for Carbopol. As shown in figure 3(b), we report that the data do not precisely follow the power law of $2/(2n_p+1)$, as predicted by (1.2). The results imply that all the dynamics of the coating process cannot be captured by using $\widehat {Ca}$, since the power-law model alone is not sufficient to describe the rheological behaviours of the working fluids around the bubble. In particular, the deviation of the experimental data for the CMCell solutions is larger than those for the Carbopol solutions when comparing to the prediction of (1.2), which can be attributed to the range of the effective shear rates. In the experiments, considering the effective shear rates in the film $\dot {\gamma }=U/h$, we obtain $Cu=O(10\unicode{x2013}10^2)$ for the CMCell solutions while $Cu=O(10^5)$ for the Carbopol solutions. The experiments for the CMCell solutions include low- to intermediate-shear-rate regions, while the experiments for the Carbopol solutions are performed at intermediate- to comparably high-shear-rate regions. Therefore, the viscosity plateau at low-shear-rate is required to be considered to obtain the better prediction for the case of the CMCell solutions, in addition to the power-law dependence at intermediate-shear-rates. Furthermore, in the vicinity of the uniform film thickness region, the fluid is at rest. However, the shear-thinning effect plays an important role at the bubble front meniscus. Approaching the meniscus, the film starts growing rapidly in the thickness and the shear rate also increases. However, the local shear rate will decrease again in the region of the re-circulating flow ahead of the bubble. Such a change of the local shear rate at different regions requires an accurate viscosity model for a correct representation of the flow physics, which also highlights the importance of a more realistic rheological model in free surface flow analyses.

3.3. Scaling of the film thickness with the effective capillary number

To explore the effect of the shear-thinning rheology on bubble characteristics, the following ordinary differential equation for the bubble profile has been obtained in the theoretical work by Picchi et al. (Reference Picchi, Ullmann, Brauner and Poesio2021) considering an Ellis fluid

(3.5)\begin{equation} \frac{\mathrm{d}^3\eta}{\mathrm{d}\xi^3}+\frac{3^\alpha}{(\alpha+2){El}^{\alpha-1}} \frac{\mathrm{d}^3\eta}{\mathrm{d}\xi^3} \left|\frac{\mathrm{d}^3\eta}{\mathrm{d}\xi^3}\right|^{\alpha-1}\eta^{\alpha-1}= \frac{\eta-1}{\eta^3}, \end{equation}

where $\xi =x/[h(3Ca_0)^{-1/3}]$ and $\eta =y/h$. Different from the Newtonian case, the bubble profile $\eta$ becomes a function of $\xi$, $\alpha$ and $El$. The two terms in the left-hand side of (3.5) represent the Newtonian and shear-thinning contributions, respectively. The curvature of the parabolic region needs to be matched to the curvature of the bubble spherical cap, $1/R$, and once that is done, the thickness of the uniform film region can be deduced.

By introducing the effective capillary number $Ca_e$ that considers both the zero-shear-rate and the shear-thinning effects, we obtain

(3.6)\begin{equation} \frac{h}{R}=P(3Ca_0)^{2/3}=0.643(3Ca_e)^{2/3} \quad \text{with}\ Ca_e=\frac{\mu_eU}{\sigma}, \end{equation}

where $\mu _e$ is the effective viscosity defined as $\mu _e=\mu _0(P/0.643)^{3/2}$, and $P$ is the dimensionless curvature related to the capillary pressure as the second derivative of $\eta$ with respect to $\xi$, i.e. $P=\mathrm {d}^2\eta /\mathrm {d}\xi ^2$ at the limit of $\eta \gg 1$. Here, $P$ is determined by numerically solving (3.5) using the fully implicit solver $ode15i$ of Matlab and the exact initial conditions given by Picchi et al. (Reference Picchi, Ullmann, Brauner and Poesio2021), and then taking a limit when $\eta \gg 1$ toward the front meniscus. Furthermore, 0.643 is the numerical factor for the Newtonian limit (Bretherton Reference Bretherton1961). In addition, Picchi et al. (Reference Picchi, Ullmann, Brauner and Poesio2021) used the numerical results of $P$ from (3.5) to obtain a master fitting curve of $\mu _e$ as a function of $El$ and $\alpha$. Here, we revise the fitting curve with $Cu$ and $n_c$ considering the experimental rheological data ($Cu\approx 1.3{El}^{-1}$ and $n_c=1/\alpha$, § 3.1) as follows:

(3.7)\begin{equation} \frac{\mu_e}{\mu_0}=\left(\frac{P}{0.643}\right)^{3/2}= \begin{cases} 1 & \text{if} \ Cu \rightarrow 0, \\ \dfrac{10-7n_c}{4}\left(\dfrac{Cu}{1.3}\right)^{n_c-1} \equiv \varTheta & \text{if} \ Cu \gg 1. \end{cases}\end{equation}

The prefactor $(10-7{n_c})/4$ for $Cu \gg 1$ comes from the fitting curve obtained by Picchi et al. (Reference Picchi, Ullmann, Brauner and Poesio2021) at $El\rightarrow 0$. When $Cu \rightarrow 0$, (3.7) reduces to a Newtonian case. When $Cu \gg 1$, the shear-thinning effect dominates, and thus the viscosity depends on $n_c$ and $Cu$ given that $\mu _\infty$ is neglected in the current experiments. For other practical non-Newtonian fluids with non-negligible $\mu _\infty$, the limit of $Cu \rightarrow \infty$ corresponds to the high-shear-rate viscosity plateau with constant viscosity $\mu _\infty$.

For each experiment, we compute the effective shear rate in the film $U/h$, using the measured bubble speed and film thickness, and thus acquire the effective viscosity $\mu _e/\mu _0$ of the CMCell and Carbopol solutions by using the C-Y model (figure 2). We find the experimental values of $\mu _e/\mu _0$ agree well with $(P/0.634)^{3/2}$, which is numerically calculated from (3.5). Indeed, figure 4(a) provides a plot of the effective viscosity as a function of $Cu$ and $n_c$ and all the viscosity data collapse around the fitting curve of the master (3.7), as shown in figure 4(b), which suggests a universal scaling for the effective viscosity to define the effective capillary number. Therefore, the comparison of the effective capillary number $Ca_e$ obtained from (3.6) against the experimental data show good agreement over the entire range of $Ca_e$ (figure 4c). The smaller values of $\mu _e$/$\mu _0$ for the Carbopol solution demonstrate the greater extent of shear-thinning compared to the CMCell solution, which are also indicated by the higher values of $Cu$ for the Carbopol solutions, as shown in table 2. We further recast (3.6) in the following expression proposed by Aussillous & Quéré (Reference Aussillous and Quéré2000) for the range of $Ca_e$ ($10^{-3}< Ca_e<0.6$) in the current experiments as

(3.8)\begin{equation} \frac{h}{R}=\frac{1.34Ca_e^{2/3}}{1+2.5\times1.34Ca_e^{2/3}}.\end{equation}

Figure 4. (a) Non-dimensional effective viscosity $\mu _e/\mu _0$ as a function of the Carreau number $Cu$ and the power-law index $n_c$. (b) $\mu _e/\mu _0$ of the experimental data as a function of $\varTheta$. All the experimental data from the present work show good agreement with the master curve of (3.7). (c) Comparison of the effective capillary number, $Ca_e$, between the theoretical and experimental results. All the data lie along the solid line with a slope of a unity. (d) h/R as a function of $Ca_e$. The black line represents prediction of (3.8) with $Ca_e$ (Aussillous & Quéré Reference Aussillous and Quéré2000; Picchi et al. Reference Picchi, Ullmann, Brauner and Poesio2021). Error bars are smaller than the symbols.

Figure 4(d) shows a plot of the non-dimensional liquid film thickness $h/R$ as a function of $Ca_e$. The values of $h/R$ increase consistently with $Ca_e$ for the non-Newtonian fluids, as the viscous effect is increasingly important. The experimental data for both the CMCell and Carbopol solutions agree well with (3.8), including the trend in low $Ca_e$ and the saturation behaviour at relatively large $Ca_e$. Thus, we demonstrate that the liquid film thickness of the shear-thinning fluids can be estimated with the scaling law proposed by Aussillous & Quéré (Reference Aussillous and Quéré2000) using $Ca_e$. The better prediction accuracy compared with (1.2) also highlights the importance of a practical rheological model. To the best of our knowledge, our work serves as the first experimental validation of (3.8) with realistic shear-thinning fluids based on $Ca_e$, which is helpful to assess the true range of applicability for different scaling laws.

3.4. Scaling of the bubble speed with the effective capillary number

In addition, we experimentally measure the ratio of the bubble speed, $U$, to the average velocity of the fluid flowing far from the bubble, $U_\infty$. The scaling law for the ratio $U/U_\infty$ can be derived by applying the mass balance near the thin film region in a reference frame moving with the bubble (Picchi et al. Reference Picchi, Ullmann, Brauner and Poesio2021). Using $Ca_e$, we obtain

(3.9)\begin{equation} \frac{U}{U_\infty}=\frac{1}{(1-h/R)^2}=\frac{1}{\left(1- \dfrac{1.34Ca_{e}^{2/3}}{1+2.5\times1.34Ca_{e}^{2/3}}\right)^2}.\end{equation}

The experimentally obtained $U/U_\infty$ collapse well with (3.9) (figure 5), showing that $Ca_e$, as a function of $Cu$ and $n_c$, can be used to describe the evolution of the bubble speed with the rheological parameters of the shear-thinning fluids.

Figure 5. Ratio of the bubble speed to the average velocity of the fluid far from the bubble $U/U_\infty$, as a function of $Ca_e$ obtained from (3.6). Error bars are smaller than the symbols.

3.5. Characteristics of the bubble front and rear menisci

We further investigate the shape variations of the bubble translating in the CMCell solution, in which the front and rear menisci of the bubble could be identified clearly compared to those in the Carbopol solution. In this section, the bubble shape profiles near the front and rear menisci are computed by solving (3.5) as a function of $El$ and $\alpha$. Then we plot in figure 6 the numerically obtained profiles considering $Cu$ and $n_c$ to be given by the correlations $Cu \approx 1.3 El^{-1}$ and $n_c=1/\alpha$, and these correlations having already been introduced back in § 3.1. Using the experimentally obtained bubble profiles, we first identify the points at which the local film thickness increase by one pixel ($\approx 1.3\ \mathrm {\mu }$m) compared with the uniform film thickness. Such a thickness increase corresponds to ${\approx }1.7$$8.5\,\%$ of the uniform thickness of the deposited film. Then, these points are overlapped with those corresponding to the same thickness increase in the numerically obtained bubble profiles, as shown in figure 6. We note that the bubble profiles at the front and the rear menisci are solved separately by integrating (3.5) with a different set of boundary conditions. At the bubble front, we assume that the thin film region extends to $\xi \rightarrow -\infty$ (i.e. $\eta (-\infty )=1$), and the front meniscus is obtained by integrating (3.5) towards positive $\xi$. However, at the rear meniscus, the thin film region is at $\xi \rightarrow +\infty$, and the profile at the rear meniscus is obtained by integrating (3.5) towards negative $\xi$ starting from the boundary condition $\eta (+\infty )=1$ (Bretherton Reference Bretherton1961; Picchi et al. Reference Picchi, Ullmann, Brauner and Poesio2021).

Figure 6. Bubble front meniscus as a function of (a) the Carreau number $Cu$ with the power-law index $n_c = 0.48$ and (b) $n_c$ with $Cu = 13.8$. Bubble rear meniscus as a function of (c) the Carreau number $Cu$ with the power-law index $n_c = 0.48$ and (d) $n_c$ with $Cu = 13.8$.

For the bubble front meniscus, figure 6(a) shows good agreement between the experimental results and the numerical predictions (Picchi et al. Reference Picchi, Ullmann, Brauner and Poesio2021) of the shape changes when $Cu$ increases at fixed $n_c=0.48$. Although the bubble maintains the rounded shape similar to that in the Newtonian fluids, the delayed transition from the uniform film to the parabolic region characterized by a constant dimensionless curvature becomes significant due to the higher effective shear rate in the film. Considering the viscosity field obtained from the numerical simulations by Moreira et al. (Reference Moreira, Rocha, Carneiro, Araújo, Campos and Miranda2020), the appearance of high viscosity in the film is due to the almost stagnant liquid, while in the axis of the channel, it is due to a low velocity gradient. In between, as $Cu$ can be interpreted as the ratio between the representative shear rate in the film to the onset of the shear-thinning effects, larger values of $Cu$ (i.e. lower $El$) indicate stronger shear-thinning effects and thus the weight of the second term on the left-hand side of (3.5) increases, resulting in the delayed transition to the parabolic region for the bubble shape. A similar trend is observed when decreasing $n_c$ at fixed $Cu=13.8$, as shown in figure 6(b). The decrease of $n_c$ indicates stronger shear-thinning effects and thus the transition to the parabolic region is also expected to be delayed. We note that the numerical solution of (3.5) starts to deviate from the experimental data as $Ca_0 > 0.3$ (figure 6b) since the lubrication approximation will no longer strictly hold for the relatively large $Ca_0$. We note that such a delayed transition from the thin film to the parabolic profile is also theoretically observed in the case where a charged oil droplet moves through a charged capillary. When the electrostatic interaction between the capillary wall and the droplet surface is attractive, the visco-electro-osmotic balance might not only reduce the film thickness, but also delay the transition because of the cooperation of the electro-osmotic and capillary pressure (Grassia Reference Grassia2020, Reference Grassia2022).

For the bubble rear meniscus, (3.5) is solved following the initial and boundary conditions given by Picchi et al. (Reference Picchi, Ullmann, Brauner and Poesio2021). The general observation of the rear meniscus is similar to that of the Newtonian case, where the bubble profile exhibits one main crest and one main valley (Magnini et al. Reference Magnini, Beisel, Ferrari and Thome2017) with a high degree of undulations, as shown in figure 6(c,d). Although the experimental results do not agree with the numerical solutions very well, the general trends are consistent with the results of Picchi et al. (Reference Picchi, Ullmann, Brauner and Poesio2021), considering the effect of the shear-thinning rheology. As the shear-thinning effect becomes more important with increasing $Cu$ or decreasing $n_c$, the bubble profile stretches along $\xi$. Unlike the front meniscus, where the viscosity profile is regular, the axial velocity gradient near the rear meniscus should be considered when computing the shear rate due to the undulations (Picchi et al. Reference Picchi, Ullmann, Brauner and Poesio2021). Therefore, an accurate description of the viscosity field at the rear meniscus of the bubble requires a further correction for the axial derivative of the velocity in future work.

4. Conclusion

In this work, we provide an experimental framework to study the motion of a long bubble translating in a circular capillary tube filled with non-Newtonian shear-thinning fluids. The Carreau–Yasuda rheological model is used to describe the rheological properties of the CMCell and Carbopol solutions, with full consideration for the viscosity plateaus at the very low- or high-shear-rates and the shear-thinning behaviour at the intermediate-shear-rates. We show that the deposited film thickness and the bubble speed cannot be scaled by the modified capillary number based on the simple power-law rheological model. Instead, the extended Bretherton's law holds well if the effective capillary number is considered, as a function of the Carreau number and power-law index in the Carreau–Yasuda rheological model. In addition, we investigate the shear-thinning effect on the variation of the bubble profile near the front and rear menisci. Based on a recent theoretical work by Picchi et al. (Reference Picchi, Ullmann, Brauner and Poesio2021), we systematically compare the experimental measurements to the numerical prediction for the bubble profile. Stronger shear-thinning effect, indicated by large Carreau number and smaller power-law index, delays the transition from the uniform film to the parabolic region at the bubble front while stretches the undulations at the rear meniscus. The numerical prediction of the bubble profile works well for the bubble front but with less accuracy for the bubble rear given the complexity of the velocity field. We believe our results serve as an experimental validation of the recent modelling approach, which provides confidence in applying these models to a variety of problems involving lubrication and coating flows with shear-thinning fluids. The influences of other rheological properties, such as viscoelasticity and the resulting coupling with the channel geometry, on the deposition dynamics will be the focus of our future investigation.

Acknowledgements

We acknowledge Professor R.H. Ewoldt and Y. Wang in Mechanical Science and Engineering at the University of Illinois at Urbana-Champaign for fruitful discussion about the rheology analysis. Rheological experiments were carried out in part in the Materials Research Laboratory Central Research Facilities, University of Illinois.

Funding

We acknowledge the support from American Chemical Society Petroleum Research Fund Grant No. 61574-DNI9 (to J.F.).

Declaration of interests

The authors report no conflict of interest.

References

REFERENCES

Abishek, S., King, A.J.C. & Narayanaswamy, R. 2015 Dynamics of a Taylor bubble in steady and pulsatile co-current flow of Newtonian and shear-thinning liquids in a vertical tube. Intl J. Multiphase Flow 74, 148164.CrossRefGoogle Scholar
Atasi, O., Khodaparast, S., Scheid, B. & Stone, H.A. 2017 Effect of buoyancy on the motion of long bubbles in horizontal tubes. Phys. Rev. Fluids 2 (9), 094304.CrossRefGoogle Scholar
Aussillous, P. & Quéré, D. 2000 Quick deposition of a fluid on the wall of a tube. Phys. Fluids 12 (10), 23672371.CrossRefGoogle Scholar
Benchabane, A. & Bekkour, K. 2008 Rheological properties of carboxymethyl cellulose (CMC) solutions. Colloid Polym. Sci. 286 (10), 11731180.CrossRefGoogle Scholar
Bird, R.B., Armstrong, R.C. & Hassager, O. 1987 Dynamics of Polymeric Liquids. Vol. 1: Fluid Mechanics. John Wiley and Sons Inc.Google Scholar
Boehm, M.W., Sarker, S. & Koelling, K. 2011 An experimental investigation of two-phase coating flow within microchannels: the effect of coating fluid rheology. Microfluid Nanofluid 10 (6), 11751183.CrossRefGoogle Scholar
Boyko, E. & Stone, H.A. 2021 Flow rate–pressure drop relation for shear-thinning fluids in narrow channels: approximate solutions and comparison with experiments. J. Fluid Mech. 923, R5.CrossRefGoogle Scholar
Bretherton, F.P. 1961 The motion of long bubbles in tubes. J. Fluid Mech. 10 (2), 166188.CrossRefGoogle Scholar
Carreau, P.J. 1972 Rheological equations from molecular network theories. Trans. Soc. Rheol. 16 (1), 99127.CrossRefGoogle Scholar
Chao, C., Jin, X. & Fan, X. 2020 Evolution of thin-liquid films surrounding bubbles in microfluidics and their impact on the pressure drop and fluid movement. Langmuir 36 (49), 1510215111.CrossRefGoogle ScholarPubMed
Clanet, C., Héraud, P. & Searby, G. 2004 On the motion of bubbles in vertical tubes of arbitrary cross-sections: some complements to the dumitrescu–Taylor problem. J. Fluid Mech. 519, 359376.CrossRefGoogle Scholar
Datt, C., Zhu, L., Elfring, G.J. & Pak, O.S. 2015 Squirming through shear-thinning fluids. J. Fluid Mech. 784, R1.CrossRefGoogle Scholar
Deryagin, B.V. & Levi, S.M. 1964 Film Coating Theory: The Physical Chemistry of Coating Thin Layers on a Moving Support. Focal Press.Google Scholar
Gao, Y., Chan, C., Gu, Q., Lin, X., Zhang, W., Yeo, D., Alsema, A., Arora, M., Chong, M. & Shi, P. 2016 Controlled nanoparticle release from stable magnetic microbubble oscillations. NPG Asia Mater. 8 (4), e260.CrossRefGoogle Scholar
Gauri, V. & Koelling, K.W. 1999 The motion of long bubbles through viscoelastic fluids in capillary tubes. Rheol. Acta 38 (5), 458470.CrossRefGoogle Scholar
Ghannam, M.T. & Esmail, M.N. 1997 Rheological properties of carboxymethyl cellulose. J. Appl. Polym. Sci. 64 (2), 289301.3.0.CO;2-N>CrossRefGoogle Scholar
Grassia, P. 2019 Motion of an oil droplet through a capillary with charged surfaces. J. Fluid Mech. 866, 721758.CrossRefGoogle Scholar
Grassia, P. 2020 Viscous and electro-osmotic effects upon motion of an oil droplet through a capillary. J. Fluid Mech. 899, A31.CrossRefGoogle Scholar
Grassia, P. 2022 Electro-osmotic and viscous effects upon pressure to drive a droplet through a capillary. Proc. R. Soc. A 478 (2258), 20210801.CrossRefGoogle ScholarPubMed
Gutfinger, C. & Tallmadge, J.A. 1965 Films of non-Newtonian fluids adhering to flat plates. AIChE J. 11 (3), 403413.CrossRefGoogle Scholar
Hernot, S. & Klibanov, A. 2008 Microbubbles in ultrasound-triggered drug and gene delivery. Adv. Drug Deliv. Rev. 60 (10), 11531166.CrossRefGoogle ScholarPubMed
Hewson, R.W., Kapur, N. & Gaskell, P.H. 2009 A model for film-forming with Newtonian and shear-thinning fluids. J. Non-Newtonian Fluid Mech. 162 (1–3), 2128.CrossRefGoogle Scholar
Jeong, D.H., Kvasnickova, A., Boutin, J.-B., Cébron, D. & Sauret, A. 2020 Deposition of a particle-laden film on the inner wall of a tube. Phys. Rev. Fluids 5 (11), 114004.CrossRefGoogle Scholar
Kamişli, F. & Ryan, M.E. 2001 Gas-assisted non-Newtonian fluid displacement in circular tubes and noncircular channels. Chem. Engng Sci. 56 (16), 49134928.CrossRefGoogle Scholar
Laborie, B., Rouyer, F., Angelescu, D.E. & Lorenceau, E. 2017 Yield-stress fluid deposition in circular channels. J. Fluid Mech. 818, 838851.CrossRefGoogle Scholar
Li, Z., Li, G., Li, Y., Chen, Y., Li, J. & Chen, H. 2021 Flow field around bubbles on formation of air embolism in small vessels. Proc. Natl Acad. Sci. 118 (26), e2025406118.CrossRefGoogle ScholarPubMed
Ma, Y.-C., Bai, F.-Q., Chang, Q., Yi, J.-M., Jiao, K. & Du, Q. 2015 An experimental study on the atomization characteristics of impinging jets of power law fluid. J. Non-Newtonian Fluid Mech. 217, 4957.CrossRefGoogle Scholar
Ma, Y., Sun, M., Duan, X., van den Berg, A., Eijkel, J. & Xie, Y. 2020 Dimension-reconfigurable bubble film nanochannel for wetting based sensing. Nat. Commun. 11 (1), 814.CrossRefGoogle ScholarPubMed
Magnini, M., Beisel, A.M., Ferrari, A. & Thome, J.R. 2017 Pore-scale analysis of the minimum liquid film thickness around elongated bubbles in confined gas-liquid flows. Adv. Water Resour. 109, 8493.CrossRefGoogle Scholar
Magnini, M., Khodaparast, S., Matar, O.K., Stone, H.A. & Thome, J.R. 2019 Dynamics of long gas bubbles rising in a vertical tube in a cocurrent liquid flow. Phys. Rev. Fluids 4 (2), 023601.CrossRefGoogle Scholar
Majeed, T., Kamal, M.S., Zhou, X. & Solling, T. 2021 A review on foam stabilizers for enhanced oil recovery. Energy Fuels 35 (7), 55945612.CrossRefGoogle Scholar
Moreira, A.I., Rocha, L.A.M., Carneiro, J., Araújo, J.D.P., Campos, J.B.L.M. & Miranda, J.M. 2020 Isolated taylor bubbles in co-current with shear thinning cmc solutions in microchannels–a numerical study. Processes 8 (2), 242.CrossRefGoogle Scholar
Morozov, A. & Spagnolie, S.E. 2015 Introduction to complex fluids. In Complex Fluids in Biological Systems, pp. 3–52. Springer.CrossRefGoogle Scholar
Moukhtari, F.E. & Lecampion, B. 2018 A semi-infinite hydraulic fracture driven by a shear-thinning fluid. J. Fluid Mech. 838, 573605.CrossRefGoogle Scholar
Myers, T.G. 2005 Application of non-Newtonian models to thin film flow. Phys. Rev. E 72 (6), 066302.CrossRefGoogle ScholarPubMed
Picchi, D., Poesio, P., Ullmann, A. & Brauner, N. 2017 Characteristics of stratified flows of Newtonian/non-Newtonian shear-thinning fluids. Intl J. Multiphase Flow 97, 109133.CrossRefGoogle Scholar
Picchi, D., Ullmann, A., Brauner, N. & Poesio, P. 2021 Motion of a confined bubble in a shear-thinning liquid. J. Fluid Mech. 918, A7.CrossRefGoogle Scholar
Pipe, C.J., Majmudar, T.S. & McKinley, G.H. 2008 High shear rate viscometry. Rheol. Acta 47 (5), 621642.CrossRefGoogle Scholar
Reiner, M. & Leaderman, H. 1960 Deformation, strain, and flow. Phys. Today 13 (9), 47.CrossRefGoogle Scholar
Ro, J.S. & Homsy, G.M. 1995 Viscoelastic free surface flows: thin film hydrodynamics of hele-shaw and dip coating flows. J. Non-Newtonian Fluid Mech. 57 (2–3), 203225.CrossRefGoogle Scholar
Rotenberg, Y., Boruvka, L. & Neumann, A. 1983 Determination of surface tension and contact angle from the shapes of axisymmetric fluid interfaces. J. Colloid Interface Sci. 93 (1), 169183.CrossRefGoogle Scholar
Song, B. & Springer, J. 1996 Determination of interfacial tension from the profile of a pendant drop using computer-aided image processing: 2. Experimental. J. Colloid Interface Sci. 184 (1), 7791.Google Scholar
Sontti, S. & Atta, A. 2017 CFD analysis of Taylor bubble in a co-flow microchannel with Newtonian and non-Newtonian liquid. Ind. Engng Chem. Res. 56 (25), 74017412.Google Scholar
de Sousa, D.A., Soares, E.J., de Queiroz, R.S. & Thompson, R.L. 2007 Numerical investigation on gas-displacement of a shear-thinning liquid and a visco-plastic material in capillary tubes. J. Non-Newtonian Fluid Mech. 144 (2–3), 149159.CrossRefGoogle Scholar
Spiers, R.P., Subbaraman, C.V. & Wilkinson, W.L. 1975 Free coating of non-Newtonian liquids onto a vertical surface. Chem. Engng Sci. 30 (4), 379395.CrossRefGoogle Scholar
Taylor, G.I. 1961 Deposition of a viscous fluid on the wall of a tube. J. Fluid Mech. 10 (2), 161165.Google Scholar
Tran, T.Q., Ahmad, M.A., Neogi, P. & Bai, B. 2016 A single pore model for displacement of heavy crude oil with carbon dioxide. SPE J. 21 (03), 08640872.CrossRefGoogle Scholar
Yamamoto, T., Suga, T., Nakamura, K. & Mori, N. 2004 The gas penetration through viscoelastic fluids with shear-thinning viscosity in a tube. Trans. ASME J. Fluids Engng 126 (2), 148152.CrossRefGoogle Scholar
Yu, Y.E., Khodaparast, S. & Stone, H.A. 2017 Armoring confined bubbles in the flow of colloidal suspensions. Soft Matt. 13 (15), 28572865.CrossRefGoogle ScholarPubMed
Zhao, Q., Ma, H., Liu, Y., Yao, C., Yang, L. & Chen, G. 2021 Hydrodynamics and mass transfer of taylor bubbles flowing in non-Newtonian fluids in a microchannel. Chem. Engng Sci. 231, 116299.CrossRefGoogle Scholar
Zhao, B., Pahlavan, A.A., Cueto-Felgueroso, L. & Juanes, R. 2018 Forced wetting transition and bubble pinch-off in a capillary tube. Phys. Rev. Lett. 120 (8), 084501.CrossRefGoogle Scholar
Figure 0

Table 1. Chronological selection of previous experimental, numerical and theoretical studies on the liquid film thickness of a long bubble translating through non-Newtonian fluids in confined geometries. The non-dimensional groups appearing above are defined as follows: $Ca_0 = {\mu _0}U/\sigma$, $Wi = {\lambda }U/H$, $\widehat {Ca}_{HB} = k_{HB}(U/R)^{n_{HB}}/(\sigma /R)$ and $Ca_e = {\mu _e}U/\sigma$. Here, $\mu _0$ is the zero-shear-rate viscosity, $\lambda$ is the relaxation time, $H$ is the half of the gap in the Hele–Shaw cell, $k_{HB}$ and $n_{HB}$ are the consistency factor and the power-law index of the Herschel–Bulkely model, respectively, and $\mu _e$ is the effective viscosity defined as $\mu _e=\mu (\dot {\gamma }={U}/{h})$. Ro & Homsy (1995) used $m$ and $k$ to represent the degree of shear and normal stress thinning, respectively, $S$ is the ratio of the solvent viscosity to the sum of the polymer and solvent viscosity, and $\delta$ is the ratio between $\lambda$ and the characteristic residence time in the gap. Laborie et al. (2017) used $a$ and $b$ are the fitting parameters and $B$ is the non-dimensional number comparing the yield stress to the capillary pressure. Zhao et al. (2021) used $W$ to be the width of the rectangular microchannel.

Figure 1

Figure 1. (a) Schematic of the experimental configuration. A cylindrical glass tube (with an inner radius of 0.47 mm) is filled with a sample solution (e.g. glycerin, carboxymethyl cellulose (CMCell; 0.5, 1.0, 1.5 and 2.0 wt$\%$), and Carbopol (0.1, 0.2 and 0.5 wt$\%$)). The central part of the circular glass capillary is submerged in a bath of a sample solution to match the refractive index of glass. Inset: schematic of a translating air bubble confined in a circular tube. (b) Typical experimental images of a long bubble as it translates in a circular capillary filled with glycerin (left) at $Ca_0 = 7.48 \times 10^{-2}$, CMCell (1.0 wt$\%$; centre) at $Ca_0 = 3.53 \times 10^{-2}$ and Carbopol (0.2 wt$\%$; right) solutions at $Ca_0 = 3.53 \times 10^{-2}$. Images contain the front meniscus and the middle part of the bubble, where the film thickness is uniform. The scale bar is 0.5 mm. (c) Schematic for the mass balance analysis regarding the deposition of a liquid film in a circular capillary tube. (d) Comparison of the liquid film thickness obtained by the image visualization, $h_i$ and the mass balance analysis of the liquid plug, $h_m$ (2.1) for experimental cases over the range of $8\times 10^{-3} < Ca_0 < 8\times 10^{2}$. All the data lie along the solid line with a slope of a unity.

Figure 2

Figure 2. (a) Rheogram of the glycerin and carboxymethyl cellulose (CMCell) solutions with different mass fractions: viscosity $\mu$ versus shear rate $\dot {\gamma }$. (b) Rheogram of the Carbopol solutions with different mass fractions: $\mu$ versus $\dot {\gamma }$. The dashed lines represent a fitting with the Carreau–Yasuda (C-Y) model. (c) Dimensionless parameter $a$ versus the power-law index $n_c$ in the C-Y model or the degree of shear-thinning $\alpha$ in the Ellis model. (d) Ellis number $El$ versus Carreau number $Cu$ in the current experiments. Error bars are smaller than the symbols.

Figure 3

Table 2. Rheological properties of the CMCell and Carbopol solutions with different mass fractions computed by the C-Y model.

Figure 4

Figure 3. (a) Non-dimensional liquid film thickness h/R as a function of $Ca_0$. The black line represents prediction of (1.1). (b) $h/R$ versus $\widehat {Ca}^{2/{2n_p+1}}$, where $\widehat {Ca}=\kappa (U/R)^{n_p}/(\sigma /R)$, with $\kappa$ and $n_p$ ranging from 0.4 to 5.4 $\textrm {Pa s}^{n_p}$ and 0.16 to 0.51, respectively. The black line represents the prediction of (1.2) with a prefactor of 1. The experimental measurements are shown as open symbols, and error bars are smaller than the symbols.

Figure 5

Figure 4. (a) Non-dimensional effective viscosity $\mu _e/\mu _0$ as a function of the Carreau number $Cu$ and the power-law index $n_c$. (b) $\mu _e/\mu _0$ of the experimental data as a function of $\varTheta$. All the experimental data from the present work show good agreement with the master curve of (3.7). (c) Comparison of the effective capillary number, $Ca_e$, between the theoretical and experimental results. All the data lie along the solid line with a slope of a unity. (d) h/R as a function of $Ca_e$. The black line represents prediction of (3.8) with $Ca_e$ (Aussillous & Quéré 2000; Picchi et al.2021). Error bars are smaller than the symbols.

Figure 6

Figure 5. Ratio of the bubble speed to the average velocity of the fluid far from the bubble $U/U_\infty$, as a function of $Ca_e$ obtained from (3.6). Error bars are smaller than the symbols.

Figure 7

Figure 6. Bubble front meniscus as a function of (a) the Carreau number $Cu$ with the power-law index $n_c = 0.48$ and (b) $n_c$ with $Cu = 13.8$. Bubble rear meniscus as a function of (c) the Carreau number $Cu$ with the power-law index $n_c = 0.48$ and (d) $n_c$ with $Cu = 13.8$.