Next Article in Journal
MEMS-Based Wavelength-Selective Bolometers
Previous Article in Journal
Evaluating Nanoparticles in Preclinical Research Using Microfluidic Systems
Previous Article in Special Issue
Hydrophobicity Influence on Swimming Performance of Magnetically Driven Miniature Helical Swimmers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Driving Forces of the Bubble-Driven Tubular Micromotor Based on the Full Life-Cycle of the Bubble

1
Hubei Key Laboratory of Theory and Application of Advanced Materials Mechanics, Department of Mechanics and Engineering Structure, Wuhan University of Technology, Wuhan 430070, China
2
“111” Project Laboratory of Biomechanics and Tissue Repair, Bioengineering College, Chongqing University, Chongqing 400044, China
*
Author to whom correspondence should be addressed.
Micromachines 2019, 10(6), 415; https://doi.org/10.3390/mi10060415
Submission received: 24 April 2019 / Revised: 10 June 2019 / Accepted: 14 June 2019 / Published: 21 June 2019
(This article belongs to the Special Issue Microswimmer)

Abstract

:
Micromotors show many advantages in practical applications, including small size, large push-to-weight ratio, and low power consumption. Micromotors have been widely used in a variety of applications, including cell manipulation, payload delivery, and removal of toxic components. Among them, bubble-driven micromotors have received great attention due to their large driving force and high speed. The driving force of the bubble-driven micromotor movement comes from the four stages of the life cycle of the bubble: nucleation, growth, slip, and ejection. At present, investigators are still unclear about the driving mechanism of the bubble-driven micromotors, the source of the driving force being still especially controversial. In response to this problem, this paper combines the mass transfer model, hydrodynamic theory, and numerical simulation to explain the driving force generated by the various stages of the life-cycle of the bubble. A mass transfer model was used to calculate the driving force of the motor contributed by the bubble nucleation and slip stage. Based on equilibrium of force and conservation of energy, a theoretical model of the driving force of the tubular micromotor in the growth and ejection stage of the bubble was established. The results show that the driving force contributed by the bubble in the nucleation and the slip stage is rather small. However, the stage of bubble growth and ejection provide most of the driving force. On further evaluating the effect of the bubble driving force on the motor speed, it was found that the growth stage plays a major role in the motion of the bubble-driven micromotor. The micromotor velocity based on the driving forces of the full life-cycle of bubbles agrees well with the experimental results.

1. Introduction

Micromotors are micro-scale structures that exhibit many advantages in practical applications, including small size, large thrust-to-weight ratio, and low power consumption [1]. Micromotors are widely used in many fields, including cell manipulation [2,3], payload transport [4,5,6], and removal of pollutants [7,8]. Bubble-driven micromotors, which can be divided into Janus particles and tubular structures, have received great attention due to their large driving force and high speed [9]. Tubular micromotors perform typically better than Janus micromotors in locomotion. Moreover, the bubble-driven micromotor does not need to be driven by an external force, has a simple structure, and the propulsion of the micromotor is hardly affected by the ion concentration. The bubble-driven micromotor can be driven to move at a fast speed of more than 100 μm/s [10,11]. For example, a micromotor with zinc as the inner wall and polyaniline as the outer wall can reach a speed of 1000 μm/s. A similar speed can be achieved with a micromotor with Pt/Co/Ti as the inner wall [12].
The bubbles are generated by the reaction between the inner wall of the micromotor and the solution. Among them, hydrogen peroxide is the most commonly used fuel. The driving force of the bubble-driven micromotors comes from the entire life cycle of the bubbles, from the generation to the detachment [9,13]. The life cycle of bubble dynamics includes four stages, namely the nucleation, growth, slip, and ejection within the motor [14,15]. First under the action of a catalytic metal, hydrogen peroxide decomposes to generate gas and form bubbles. The bubbles nucleate at the inner surface of the micromotor wall to form spherical bubbles. Second, after the nucleation of the bubble, the outer boundary of the bubble moves outward under the action of the internal pressure of the bubble. Then the radius of the bubble becomes larger, and the bubble enters the growth stage [16]. Third, the bubble grows to a certain stage, moves to the opening under the action of the fluid, and enters the slip stage. Finally, the bubble moves to the opening of the motor and is ejected away from the tubular motor under different forces.
Investigators are still deep in controversy about the source of the driving force of bubble-driven tubular micromotors. Some authors believe that ejection contributes the most to micromotor motion, while other researchers believe that bubble growth plays a major role. These studies have generally focused on the effects of bubbles on the micromotor motion in a single stage. Solovev et al. [17] suggested that the driving force mainly from bubble ejection accelerates the tubular micromotor to 2 mm/s. Li et al. [12] considered that the driving force of the bubble-driven motor mainly comes from the force at the bubble growth stage. Manjare et al. [18] also reported that bubble growth mainly generated the driving force. Subsequently, this group gave the opposite conclusion. They found that the driving force of bubble growth is much smaller than that of bubble ejection [19]. Fomin’s group [20] proposed that both bubble growth and ejection produced driving forces and established the corresponding theoretical models for prediction. The experimental results showed that the theoretical prediction velocity is two orders of magnitude smaller than the measured velocity. That is, the model underestimates the driving force generated by the bubble. Gibbs et al. [21] suggested that the growth of bubbles provides a driving force for the Janus micromotor.
As for the debate about the source of the driving force of bubble-driven micromotors, it is very important to quantify the driving force generated at each stage of the life cycle of the bubbles. Our previous work established a model of micromotor motion [22] based on geometric asymmetry and fluid viscosity, and predicted the resistance of the cone micromotor immersed in the flow field based on the hydrodynamic model [23]. We also discussed strategies to make fast, efficient bubble-driven micromotors [24], and micromotors powered by biocompatible fuel [25]. The paper further proposed a model of driving force based on the tubular micromotor of the full life cycle of the bubble. In the bubble nucleation stage, the mass transfer model was used to propose the driving force model of the tubular micromotor in the nucleation stage. For the bubble growth stage, the force on the bubble changes as the bubble radius changes. Based on this, a driving force model of the bubble growth stage was proposed. In the slip stage of the bubble, the mass transfer model was used to establish the driving force of the motor. At the ejection stage, a driving force model of the tubular micromotor in the stage was established based on the force. In this paper, the driving force models of bubble nucleation, growth, slip, and ejection are established respectively, and the influence of bubble driving force on the motor motion is further evaluated.

2. Materials and Methods

2.1. Bubble Nucleation

According to the classical theory of bubble nucleation [26], the free energy required for the nucleation of bubbles on the inner wall of the micromotor is
W = σ L G A L G + ( σ S G σ S L ) A S G + Δ G V V G + x ( μ G μ L )
The subscripts S, L, and G represent the solid state, liquid state, and gaseous three-stage state, respectively, where σLG is the surface tension coefficient σ of the hydrogen peroxide solution required in the calculation. The equilibrium state is shown in Figure 1. According to Young’s equation of the contact angle
σ S L = σ S G + σ L G cos ( π θ )
and
m = Z R = cos ( π θ ) = cos θ = σ S L σ S G σ L G
The liquid and gas contact area ALG and the solid and gas contact area ASG are respectively
A L G = 2 π R ( R Z ) = 2 π R 2 ( 1 m )
A S G = π r 2 = π R 2 ( 1 m 2 )
Gas volume is
V G = 1 3 π ( 3 R H ) H 2 = 1 3 π R 3 ( 2 3 m + m 3 ) = 4 3 π R 3 F
The change of free energy per unit volume can be obtained according to the Laplace equation,
Δ G V = P L P G = 2 σ R
The difference in chemical potential is
μ G μ L = K T ln P G P V = K T ln P L + 2 σ / R P L + 2 σ / R C
The ideal gas equation PV = xKT can be applied to obtain the free energy required for bubble nucleation
W = 4 3 π R 2 σ F + 4 3 π R 3 σ F ( P L + 2 σ R ) ln P L + 2 σ / R P L + 2 σ / R C
The critical radius Rc can be obtained by the following formula
R c = 2 σ Δ G V = 2 σ P P G
P is the standard atmospheric pressure, and PG is the internal pressure of the nucleation bubble.
As shown in Figure 2, the volume of the bubble is partially Vc = 4π R c 3 /3. Assuming that the volume of fluid discharged along the left orifice is fVc, the volume of fluid discharged along the right orifice is (1 − f)Vc. As can be seen from the figure, due to the incompressibility of the fluid, the volume Vc of the nucleation bubble is not the volume of the fluid in the tube. The fluid discharged from the tube flows out from the two openings of the tube, respectively, wherein the ratio of the volume of the left outflow fluid to the total flow is f. The parameter f is obtained from FLUENT 18.0 (ANSYS).
The fluid flow rates at the left and right openings are
v 1 = d ( f V c ) d t ,   v 2 = d [ ( 1 f ) V c ] d t
In the nucleation stage, the driving force of the tubular micromotor is
F jet 1 = d ( m v ) d t = d d t ( m 1 v 1 m 2 v 2 ) = d d t [ ρ l f V c f V c A 1 ρ l ( 1 f ) V c ( 1 f ) V c A 2 ]
Substituting into the tubular micromotor two-port area A1, A2, and the function of bubble and oxygen production frequency q
A 1 = π R min 2 = π ( R max L tan δ ) 2 ,   A 2 = π R max 2 ,   R b = ( 3 q t 4 π ) 1 3 [ 0 , R c ]
Therefore, the driving force in the stage is
F jet 1 = ρ l π C H 2 O 2 2 n 2 L 2 cos 2 δ ( 2 R max L tan δ ) 2 [ f 2 π ( R max L tan δ ) 2 ( 1 f ) 2 π R max 2 ]
where CH2O2 is the concentration of hydrogen peroxide and f is a function of time t.

2.2. Bubble Growth

After the bubbles nucleate on the tube wall, the generated gas continues to fill into the bubbles and the bubbles gradually grow, as shown in Figure 3.
The force of a bubble along the pipe wall (a axis), perpendicular to the pipe wall (p axis) and parallel to the flow direction (x axis), perpendicular to the flow direction (y axis) is as follows
{ Σ F a = F s a + F q s cos δ + F d u a ( F b + F s l ) sin δ Σ F p = F s p + F q s sin δ + F d u p + ( F b + F s l ) cos δ + F h + F c p
{ F x = Σ F a cos δ + Σ F p sin δ F y = Σ F p cos δ Σ F a sin δ
Fs is the surface tension, Fqs is the fluid resistance of the bubble, Fdu is the unstable force of the asymmetric growth of the bubble, Fb is the buoyancy, Fsl is the shearing force, Fh is the force caused by the dynamic pressure, and Fcp is the contact pressure. The subscript a, p denotes a semi-cone angle of a conical micromotor along the wall of the tube and perpendicular to the wall of the tube.
Define contact angle function γ(φ)
r ( ϕ ) = β + ( α β ) ϕ π , 0 ϕ π
At the same time, the contact angle function needs to meet the conditions γ(0) = β, γ(π) = α, γ′(0) = γ′(π) = 0. Substituting the contact angle function, the projection of the surface tension on the coordinate axis
{ F s a = 0 π d w σ cos γ cos ϕ d ϕ = d w σ π ( α β ) π 2 ( α β ) 2 ( sin α + sin β ) F s p = 0 π d w σ sin γ d ϕ = d w σ π α β ( cos β cos α )
At the beginning of the bubble growth, the bubble is spherical. Then dw = 2Rsinθ. The volume of the bubble is expressed as V = 1 / 3 π R 3 ( 2 + 3 cos θ cos 3 θ ) .
While the growth of the bubble is in accordance with V = qt [27,28], and the bubble radius is equal to
R = γ t n = [ 3 q π ( 2 + 3 cos θ cos 3 θ ) ] 1 / 3 t 1 / 3
The contact angle θ = 30° and the semi-cone angle φ = 4.5° are applied [29].
The fluid resistance of the bubble is [30]
F q s = ( 2 3 + 1 12 R e ) 6 π μ V y R ,   R e = ρ l V y 2 R μ
where Re is the fluid Reynolds number, ρl is the fluid density, μ is the fluid viscosity, R is the bubble radius, and Vy is the fluid velocity at the bubble center point.
Substituting the average diameter to obtain the axial velocity distribution inside the conical tube
V y = 8 Q π h 2 32 Q π h 4 y 2 = 8 Q π ( 2 R max L tan δ ) 2 32 Q π ( 2 R max L tan δ ) 4 y 2
Substituting Vy into the projection of surface tension on the coordinate axis, and considering the fluid resistance of small bubbles is 2/3 of the same size solid particles [31], the fluid resistance of the bubble is expressed as
F q s = 2 3 6 π μ R [ 8 Q π ( 2 R max L tan δ ) 2 32 Q π ( 2 R max L tan δ ) 4 ( R cos θ h 2 ) 2 ]
During the growth of the bubble, the deformation of the bubble results in an asymmetrical growth force which can be obtained through a series of simplifications
F d u = Γ ρ l ( R R ¨ + 3 2 R ˙ 2 ) d A = π R 2 ρ l ( R R ¨ + 3 2 R ˙ 2 ) = 1 18 π ρ l ( 3 q π ( 2 + 3 cos θ cos 3 θ ) ) 4 / 3 t 2 / 3
Bubble buoyancy
F b = V ( ρ l ρ b ) g = 1 3 π R 3 ( 2 + 3 cos θ cos 3 θ ) ( ρ l ρ b ) g
where ρb is the density of the gas within the bubble, namely the oxygen density.
The fluid dynamic pressure Fh mainly acts on the circular surface where the bubble contacts the tube wall [32], namely
F h = 9 8 ρ l V y 2 π d w 2 4 = 9 8 ρ l π d w 2 4 [ 8 Q π ( 2 R max L tan δ ) 2 32 Q π ( 2 R max L tan δ ) 4 ( R cos θ h 2 ) 2 ] 2
The contact pressure is
F c p = π d w 2 4 2 σ r r = π d w 2 4 2 σ 5 R
The shear lift force is [33]
F s l = 2 3 81.2 μ V y R 2 κ 1 2 / ν 1 2
where μ is the dynamic viscosity, kinematic viscosity is ν = μ/ρl, velocity gradient κ = | d V y d y | = | 64 Q π ( 2 R max L tan δ ) 4 y | .
A conversion factor λ needs to be set as the force of the bubble acting on the tubular micromotor. Wang et al. [34] used experimental and molecular dynamics simulations to find that the chemical energy conversion efficiency of bubble-driven motors is about 10−10, but the conversion efficiency of bubbles was not studied. Fomin et al. [20] then studied the bubble-driven tubular micromotor and found that the conversion force of the bubble is converted to a tubular motor driving force with a conversion coefficient of about 1/30. The driving force of the tubular micromotor in the stage of bubble growth is
F jet 2 = λ F x

2.3. Bubble Slip

When the bubble grows up to the moment of contact with the tube wall, the relationship between the bubble radius Ri and the inner diameter Hi of the tube wall at the nucleation point is H i = 2 R i cos δ , or H i = 2 X i tan δ + 2 R min .
If the bubble is kept in the micromotor, there exists the expression XiC/2 ≥0, Xi + C/2 ≤L.
As shown in Figure 4a, the volume of the fluid in the tube is divided into two parts. Then the left and right spheres and the volume of the trough are
V 1 = π 3 h 1 2 ( 3 R 1 h 1 ) = π 3 [ ( X i C / 2 ) tan δ + R min cos ( θ + δ ) ] 3 [ 2 3 sin ( θ + δ ) + sin 3 ( θ + δ ) ] V 2 = π 3 h 2 2 ( 3 R 2 h 2 ) = π 3 [ ( X i + C / 2 ) tan δ + R min cos ( θ δ ) ] 3 [ 2 3 sin ( θ δ ) + sin 3 ( θ δ ) ] V 3 = π 3 C ( H 1 2 / 4 + H 1 H 2 / 4 + H 2 2 / 4 ) = π 3 C [ R 1 2 cos 2 ( θ + δ ) + R 1 R 2 cos ( θ + δ ) cos ( θ δ ) + R 2 2 cos 2 ( θ δ ) ] = π 3 C [ R 1 2 cos 2 ( θ + δ ) + R 1 R 2 ( cos 2 θ sin 2 δ ) + R 2 2 cos 2 ( θ δ ) ]
The total volume of the bubble is V = q t = C h d n S t , the expression is
V = V 1 + V 2 + V 3 = π 3 R 1 3 [ 2 3 sin ( θ + δ ) + sin 3 ( θ + δ ) ] π 3 R 2 3 [ 2 3 sin ( θ δ ) + sin 3 ( θ δ ) ] + π 3 C [ R 1 2 cos 2 ( θ + δ ) + R 1 R 2 ( cos 2 θ sin 2 δ ) + R 2 2 cos 2 ( θ δ ) ]
The left and right fluid volumes are
{ V l = π 3 X i ( R min 2 4 + R min H 1 4 + H 1 2 4 ) V 1 V r = π 3 ( L X i ) ( H 2 2 4 + R max H 2 4 + R max 2 4 ) V 2
The fluid flow rates at the left and right openings are
v l = 1 A 1 d V l d t ,   v r = 1 A 2 d V r d t
Then the driving force of the tubular micromotor is
F jet 3 = ρ l 1 A 1 ( d V l d t ) 2 ρ l 1 A 2 ( d V r d t ) 2

2.4. Bubble Ejection

The asymmetric deformation of the bubble due to buoyancy is not considered at the beginning, as shown in Figure 5.
The radius of the two spheres is R1 and R2, respectively, and the volume of the two spheres is V1 and V2. R 1 = R max / cos δ , H 1 = R 1 R max tan δ = R max / cos δ R max tan δ
V 1 = π 3 H 1 2 ( 3 R 1 H 1 ) = π R max 3 3 cos 3 δ ( sin 3 δ 3 sin δ + 2 )
and R 2 = R max / sin θ , H 2 = R 2 + R max / tan θ = R max / sin θ + R max / tan θ
V 2 = π 3 H 2 2 ( 3 R 2 H 2 ) = π R max 3 3 sin 3 θ ( cos 3 θ + 3 cos θ + 2 )
Total bubble volume is
V = V 1 + V 2 = π R max 3 3 [ sin 3 δ 3 sin δ + 2 cos 3 δ + cos 3 θ + 3 cos θ + 2 sin 3 θ ]
The bubble moves outward to the opening, where θ ( π / 2 , θ t ) , θt is the angle at which the bubble is ejected. The main change in bubble volume is reflected in the V2 part,
d V 2 = π R max 3 ( 1 cos θ ) 2 sin 4 θ d θ
The forces applied to the stage of bubble ejection are shown in Figure 6.
The surface tension
F s x = F s sin θ + F s cos δ = 2 π σ R max ( cos δ sin θ )
If the bubble growth obeys R = γin, the equivalent bubble radius R′ can be found.
F q s = ( 2 3 + 1 12 R e ) 6 π μ u s R ,   R e = ρ l u s 2 R μ
The unsteady force of asymmetric growth of bubbles is
F d u = π R 2 ρ l ( R R ¨ + 3 2 R ˙ 2 ) = 1 18 π ρ l ( 3 q 4 π ) 4 / 3 t 2 / 3
If only considering the buoyancy of the portion of the opening that is exposed to the volume of the V2 portion,
F b = ( ρ l ρ b ) g V 2
In this stage, the bubble is only in contact with the opening, and the stage radius dw is substantially close to zero. Consequently, the flow pressure, contact pressure, and shear lift force are zero
F h = 0 ,   F c p = 0 ,   F s l = 0
The resultant force of the bubble in the fluid flow direction and the vertical direction is
{ F x = F s + F q s + F d u x F y = F b + F d u y
According to the flow rate and speed of the openings of the tubular micromotor, the driving force can be obtained
F jet 4 = ρ l 1 A 1 ( d V d t ) 2 ρ l 1 A 2 ( d V d t ) 2
During the stage of bubble ejection, the velocity of the fluid flowing through the opening is approximately the same as the velocity of the ejected bubble, which is,
d V d t = π R max 2 v b
where vb is the speed at which the bubble is ejected from the opening.
Assuming that the bubble and the micromotor wall are smooth and frictionless, the work done due to the difference in surface tension on the bubble is all converted into the kinetic energy of the bubble, namely
1 2 m b v b 2 = W S 1 W S 2
The bubble is ejected from the opening at a speed of vb
{ W S 1 = 0 L X i + C / 2 F S 1 d x = π σ cos ( θ + δ ) [ ( L X i ) 2 tan δ + 2 R min ( L X i + C / 2 ) ] W S 2 = 0 L X i C / 2 F S 2 d x = π σ cos ( θ δ ) [ ( L X i ) 2 tan δ + 2 R min ( L X i C / 2 ) ]
The speed when the bubble is ejected is
v b = 2 π σ ρ g V g { cos θ cos δ R min C 2 sin θ sin δ [ ( L X i ) 2 tan δ + 2 R min ( L X i ) ] } 1 2
Then driving force during the stage could be expressed as
F jet 4 = 4 π 3 σ 2 ρ l R max 2 ( R max 2 R min 2 ) ρ g 2 V g 2 R min 2 { cos θ cos δ R min C 2 sin θ sin δ [ ( L X i ) 2 tan δ + 2 R min ( L X i ) ] }

2.5. Micromotor Velocity

Here we show a simplified version of the micromotor motion theory based on driving force and drag force. A detailed theoretical study of drag force can be found in Wang et al. [23]. According to Newton’s second law, the movement of a tubular micromotor immersed in a fluid can be expressed as
m v ˙ = Σ F = F jet F d
where Fd is the drag force experienced by the micromotor. The relationship between Fd and Reynolds number, micromotor geometry, and the drag coefficient is obtained by the computational fluid dynamics software FLUENT 18.0 (ANSYS).
The mass of the micromotor m is
m = ρ V j = ρ [ π 3 L ( R max 2 + R max R min + R min 2 )    π 3 L ( R max 2 + R max R min + R min 2 ) ]
The drag force is [23]
F d = π 4 μ v R max [ ( a e b λ + c ξ + d ξ 2 ) tan δ + f ξ ]
Solving the differential equation of Equation (51), the expression of the motion velocity of the tubular micromotor could be obtained on the basis of driving force and drag force.

3. Results and Discussion

3.1. Bubble Nucleation and Growth

The relationship between the nucleation free energy change W and the nucleus radius R is shown in Figure 7.
When the bubble nucleates on the inner wall of the micromotor, the volume fraction of fluid discharged from the left end is f = 0.358. The computation is fulfilled by FLUENT. The relationship between the moving speed of the tubular micromotor and the solution surface tension coefficient, the viscosity, concentration, motor tube length, opening radius, and the semi-cone angle were calculated.
The data in Figure 8 is fitted to obtain the relationship between the driving force and the surface tension coefficient, Fjet2 = f (σ2). As can be seen from Figure 8b, when the viscosity changes from 0.1 mPa·s to 4.0 mPa·s, the driving force remains substantially unchanged. Thus, the fluid viscosity almost has no influence on the driving force in the stage. When the concentration of the hydrogen peroxide solution in Figure 9a changes, the speed does not change much, and the difference between the maximum and minimum values is only 0.62%. Basically, it can be considered that the change in the concentration of the peroxide solution almost has no effect on the driving force. Fitting the curve in Figure 9b to obtain the relationship between the driving force and the opening radius F jet 2 = f ( R max 3 ) .The relationship with the micromotor length and semi-cone angle is linear, as demonstrated by Figure 10a,b. θ is the contact angle of the bubble with the micromotor wall, and φ is the inclination angle of the bubble under the action of the flow rate. The front and back contact angle α and β of the surface tension in Equation (18) are determined by both contact angels. Furthermore, Equation (19) used for calculation of the bubble radius R is related to the contact angle θ. Considering Equations (20)–(27) are all related to the bubble radius R, the formula of the driving force could be expressed as a function of θ and φ, namely f(θ) and g(φ).
Combing the relationship between Fjet2 and σ as well as Rmax, the final expression of the driving force in the stage is
F jet 2 = χ ρ l σ 2 R max 3 tan δ f ( θ ) g ( φ )
where χ is a conversion coefficient, which means that the resultant force of the fluid and solid motor acting on the bubble along the direction of the fluid flow cannot be completely converted to the driving force of the micromotor. The conversion coefficient χ can be obtained from experimental results in the literature [20]. It should be noted that Equation (53) is only used to characterize the relationship between the driving force and other parameters. In fact, the driving force is obtained from Equations (16) and (28).

3.2. Driving Force of the Full Life-Cycle of the Bubbles

Taking a peroxide solution with a concentration of 5% as an example, the density is ρl = 103 kg/m3, the viscosity μ = 0.9 mPa⋅s, the surface tension coefficient σ = 30 mN/m. The geometry of the tubular micromotor has: length L = 100 μm, larger opening radius Rmax = 10 μm and semi-cone angle δ = 3.2°. The driving forces at each stage are shown in Table 1. A phase diagram of the bubble in the full life cycle of bubbles can be found in Figure 11.
As can be seen from Table 1, the bubbles contribute to the movement of the tubular micromotor at various stages of micromotor motion. The driving forces contributed by bubble nucleation and slip are insignificant. It is worth mentioning that, the driving force contributed by the nucleation stage is almost negligible. The driving forces contributed by bubble growth and ejection are prominent, as can be seen from Table 1. Bubbles provide most of the driving force during the stage of growth and ejection. Our results are a quantitative comparison of the driving forces generated by bubbles in the four stages of the full life cycle of bubbles, which is more comprehensive than the previous reports on the individual stages of the research [17,18,19]. The four stages of the bubble life cycle are strictly distinguished and compared quantificationally separately. A similar conclusion provided by ref. [12] is from experimental observation and is only a qualitative conclusion while our results further provide a quantitative description to arrive at the conclusion that the driving force is mainly from the growth stage. For example, the driving force of the growth stage accounts for the ratio of the total driving force. The results also confirmed results reported by Fomin et al. that both bubble growth and ejection generate the driving force [20]. Furthermore, our results used detailed computations to reveal that bubble growth and ejection provide the primary driving force for micromotors.

3.3. The Influence of Driving Force on Speed in Each Stage

Taking the Reynolds number Re = 0.0125 in Table 2 as an example, FLUENT. calculates the pressure distribution cloud diagram and streamline diagram on the tubular motor as shown in Figure 12.
The results show that the pressure near the upstream ring is larger and the pressure at the tail of the motor tube becomes smaller. The drag forces of all the walls of the tubular micromotor were monitored separately, including the upstream ring, the backflow ring, the outer wall, and the inner wall. Drag forces of various areas are shown in Table 2.
The results demonstrated that the drag forces of the outer wall were the smallest. The drag forces on the region of the upstream and backflow ring were small, and the drag forces on the inner wall were the largest. Therefore, the drag forces of the tubular motor mainly come from the inner wall of the micromotor body. Considering the drag forces of the inner wall of the tube, we obtained a more accurate drag force for conical tubular micromotors
F d = π 4 μ v R max [ ( a e b λ + c ξ + d ξ 2 ) tan δ + f ξ ]
The effect of the driving force generated by the bubble at each stage on the motor speed was evaluated. The dimensions and fluid parameters for a given tubular micromotor are as follows: 5% hydrogen peroxide solution, density ρl = 1130 kg/m3, viscosity μ = 0.9 mPa⋅s, surface tension coefficient σ = 30 mN/m; micromotor length L = 100 μm, larger opening radius Rmax = 10 μm, semi-cone angle δ = 3.2°; solution density ρg = 1.33, contact angle α = 45°, β = 36°, angle of inclination φ = 4.5°. The viscosity μ = 0.9 mPa⋅s that we used falls into the general viscosity of the fluid ranges from 0.1–4 mPa·s. The generation rate of bubbles in the chemical reaction is n = 9.8 × 10−4 m/s. For example, when the random nucleation point is in the middle of the tube, the stage of bubble ejection has a shorter duration. The velocity of the tubular micromotor in the various stages of the bubble is as shown in Figure 13.
As can be seen in Figure 13, the nucleation stage takes the longest time and the ejection stage takes the shortest time during the entire bubble cycle. The stage of bubble nucleation generates driving forces but has a limited contribution to the motion. The growth stage plays a major role in the micromotor motion, as seen in Figure 13. The speed of the motor rapidly increases in the stage. However, the driving force in the slip stage drops rapidly and the speed of the micromotor is accelerated to a higher level than the previous stage. The drag forces in the slip stage are also larger, even exceeding the driving force temporally. Therefore, the micromotor exhibits a decrease in speed in the stage. On the contrary, a decrease in the speed reduces fluid resistance, causing the speed to decrease more and more slowly. The bubble finally approaches an equilibrium state until it leaves the opening and bursts [18].
It can be seen from Figure 13, the velocity of the bubble in the ejection stage is increased, but the increase is smaller than the growth stage. Although the driving force in the ejection stage is large, the duration of the stage is too short to make a comparable contribution to the micromotor motion. Among the four stages of the bubble life cycle, the most prominent important contribution to micromotor motion is the growth stage, as seen in Figure 13 and the most insignificant contribution to micromotor motion is the nucleation stage. The contribution of the whole bubble cycle to the velocity of the tubular micromotor is averaged by integrating the velocity of the whole life-cycle of the bubble using the formula. The calculations show that the average velocity is about 6.49 μm/s when the bubbles nucleate in the middle of the micromotor.
Bubbles are produced in large quantities at a certain frequency to propel the micromotor ahead [35]. A further consideration is given to the effect of the generation of numerous bubbles on the motor motion. The frequency of generated bubbles used in the analysis was the same as that of Li et al. [5]. When the semi-cone angle of the tubular micromotor increased from 0° to 5°, the frequency of bubbles decreased from 97 Hz to 54 Hz which can be obtained by linear interpolation. When the semi-cone angle is 3.2°, the bubble generation frequency should be 69.48 Hz. Micromotor geometry used in the analysis has length L = 100 μm, the larger opening radius Rmax = 10 μm, semi-cone angle δ = 3.2°. When moving in 5% peroxide solution, the average speed of motion is 450.93 μm/s which is in accordance with 431 μm/s reported by Li et al [5].
We further validated our theory of driving forces with experimental measurements. The sample was placed in an aqueous solution of hydrogen peroxide, and the speed of the micromotor was measured to compare with the theoretical results. Figure 14 shows the images of optical microscopy and scanning electron microscopy of the tubular micromotor prepared by ourselves. The parameters for the micromotor are as follows: 5% hydrogen peroxide solution, density ρl = 1130 kg/m3, viscosity μ = 0.9 mPa⋅s, surface tension coefficient σ = 30 mN/m, micromotor length L = 100 μm, larger opening radius Rmax = 10 μm, semi-cone angle δ = 3.2°, solution density ρg = 1.33, contact angle α = 45°, β = 36°, angle of inclination φ = 4.5°.
Because the micromotor used in the experiment has a small size, the bubbles quickly contact the wall after nucleation and then go to the slip stage. In this case, the life-cycle of the bubbles appears as three stages in the absence of the growth stage. The velocity of the micromotor in the life-cycle can be found in Figure 15. The average velocity of the micromotor in a single bubble cycle computed by our theory is 7.39 μm/s. The experimentally measured frequency of bubbles is 75 Hz. The average speed obtained from theory is 554.25 μm/s. Moreover, the average speed obtained from the measured trajectory is 540 μm/s. Consequently, the experimental results agree very well with the theoretical prediction.

4. Conclusions

The dynamical behavior of the bubble plays a crucial role throughout the movement of the bubble-driven tubular micromotor. In this paper, based on the unclear problem of the driving mechanism of the bubble-driven micromotor, the theoretical model of the driving force of the motor in different stages was established by using mass transfer model, hydrodynamic theory, and numerical simulation. Based on the bubble nucleation theory, a theoretical calculation model of the driving force of the tubular micromotor in the nucleation stage of the bubble was determined. Based on the principle of force balance, the theoretical calculation model of the driving force of the tubular micromotor in the growth stage of the bubble tube was established. Based on the mass transfer theory, a theoretical model of the driving force of the motor in the slip stage of the bubble tube was derived. Based on the principle of force and energy conservation, a theoretical model of the driving force of the tubular micromotor in the ejection stage of the bubble was formed. The results show that the driving forces produced by bubble nucleation and the slip are insignificant. The driving force contributed by the bubble nucleation is the smallest and is almost negligible. The stage of bubble growth and ejection provide most of the driving force. Further, comparing the contribution of the driving force at each stage to the motor speed, it was found that the growth stage plays a major role in the locomotion of the bubble-driven micromotor. The micromotor speed based on the driving force of the bubble life-cycle is in consistent with the experimental results.

Author Contributions

Y.L. implemented the analysis and wrote the manuscript; X.G. and Q.C. contributed to the writing and discussion of the paper; K.X. and C.W. made some analysis of the results; Z.W. proposed the idea of the manuscript.

Funding

This study was supported by grants from the National Natural Science Foundation of China (11602181), the Fundamental Research Funds for the Central Universities (WUT: 2018IB005), the Visiting Scholar Foundation of Key Laboratory of Biorheological Science and Technology(Chongqing University), Ministry of Education (Grant Number: CQKLBST-2018-006, CQKLBST-2018-009), the open foundation of Hubei Key Laboratory of Theory and Application of Advanced Materials Mechanics (Grant Number: TAM201813), the Independent Innovation Foundation of Wuhan University of Technology (WUT: 2019IⅡ075GX).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gao, W.; Uygun, A.; Wang, J. Hydrogen-bubble-propelled zinc-based microrockets in strongly acidic media. J. Am. Chem. Soc. 2012, 134, 897–900. [Google Scholar] [CrossRef] [PubMed]
  2. Balasubramanian, S.; Kagan, D.; Jack Hu, C.-M.; Campuzano, S.; Lobo-Castañon, M.J.; Lim, N.; Kang, D.Y.; Zimmerman, M.; Zhang, L.; Wang, J. Micromachine-enabled capture and isolation of cancer cells in complex media. Angew. Chem. 2011, 123, 4247–4250. [Google Scholar] [CrossRef]
  3. Yoshizumi, Y.; Okubo, K.; Yokokawa, M.; Suzuki, H. Programmed transport and release of cells by self-propelled micromotors. Langmuir 2016, 32, 9381–9388. [Google Scholar] [CrossRef] [PubMed]
  4. Wu, Z.; Esteban-Fernández de Ávila, B.; Martín, A.; Christianson, C.; Gao, W.; Thamphiwatana, S.K.; Escarpa, A.; He, Q.; Zhang, L.; Wang, J. Rbc micromotors carrying multiple cargos towards potential theranostic applications. Nanoscale 2015, 7, 13680–13686. [Google Scholar] [CrossRef] [PubMed]
  5. Li, L.; Wang, J.; Li, T.; Song, W.; Zhang, G. Hydrodynamics and propulsion mechanism of self-propelled catalytic micromotors: Model and experiment. Soft Matter 2014, 10, 7511–7518. [Google Scholar] [CrossRef] [PubMed]
  6. Li, J.; Angsantikul, P.; Liu, W.; Esteban-Fernández de Ávila, B.; Thamphiwatana, S.; Xu, M.; Sandraz, E.; Wang, X.; Delezuk, J.; Gao, W.; et al. Micromotors spontaneously neutralize gastric acid for ph-responsive payload release. Angew. Chem. Int. Ed. 2017, 56, 2156–2161. [Google Scholar] [CrossRef] [PubMed]
  7. Zhang, Q.; Dong, R.; Wu, Y.; Gao, W.; He, Z.; Ren, B. Light-driven au-wo3@c janus micromotors for rapid photodegradation of dye pollutants. ACS Appl. Mater. Interfaces 2017, 9, 4674–4683. [Google Scholar] [CrossRef] [PubMed]
  8. Eskandarloo, H.; Kierulf, A.; Abbaspourrad, A. Nano- and micromotors for cleaning polluted waters: Focused review on pollutant removal mechanisms. Nanoscale 2017, 9, 13850–13863. [Google Scholar] [CrossRef]
  9. Hayakawa, M.; Onoe, H.; Nagai, K.H.; Takinoue, M. Influence of asymmetry and driving forces on the propulsion of bubble-propelled catalytic micromotors. Micromachines 2016, 7, 229. [Google Scholar] [CrossRef]
  10. Liang, Y.; Wang, H.; Yao, D.; Chen, Y.; Deng, Y.; Wang, C. Transportation and release of janus micromotors by two-stage rocket hydrogel. J. Mater. Chem. A 2017, 5, 18442–18447. [Google Scholar] [CrossRef]
  11. Yang, F.; Manjare, M.; Zhao, Y.; Qiao, R. On the peculiar bubble formation, growth, and collapse behaviors in catalytic micro-motor systems. Microfluid. Nanofluid. 2017, 21, 6. [Google Scholar] [CrossRef]
  12. Li, J.; Huang, G.; Ye, M.; Li, M.; Liu, R.; Mei, Y. Dynamics of catalytic tubular microjet engines: Dependence on geometry and chemical environment. Nanoscale 2011, 3, 5083–5089. [Google Scholar] [CrossRef] [PubMed]
  13. Gallino, G.; Zhu, L.; Gallaire, F. The hydrodynamics of a micro-rocket propelled by a deformable bubble. Fluids 2019, 4, 48. [Google Scholar] [CrossRef]
  14. Kulkarni, A.A.; Joshi, J.B. Bubble formation and bubble rise velocity in gas−liquid systems:  A review. Ind. Eng. Chem. Res. 2005, 44, 5873–5931. [Google Scholar] [CrossRef]
  15. Bari, S.A.; Hallett, J. Nucleation and growth of bubbles at an ice–water interface. J. Glaciol. 1974, 13, 489–520. [Google Scholar] [CrossRef]
  16. Kundumani-Sridharan, V.; Dyukova, E.; Hansen, D.E.; Rao, G.N. 12/15-lipoxygenase mediates high-fat diet-induced endothelial tight junction disruption and monocyte transmigration: A new role for 15(s)-hydroxyeicosatetraenoic acid in endothelial cell dysfunction. J. Biol. Chem. 2013, 288, 15830–15842. [Google Scholar] [CrossRef]
  17. Solovev, A.A.; Yongfeng, M.; Esteban, B.U.A.; Gaoshan, H.; Schmidt, O.G. Catalytic microtubular jet engines self-propelled by accumulated gas bubbles. Small 2010, 5, 1688–1692. [Google Scholar] [CrossRef]
  18. Manjare, M.; Yang, B.; Zhao, Y.P. Bubble driven quasioscillatory translational motion of catalytic micromotors. Phys. Rev. Lett. 2012, 109, 128305. [Google Scholar] [CrossRef]
  19. Manjare, M.; Yang, B.; Zhao, Y.P. Bubble-propelled microjets: Model and experiment. J. Phys. Chem. C 2013, 117, 4657–4665. [Google Scholar] [CrossRef]
  20. Fomin, V.M.; Hippler, M.; Magdanz, V.; Soler, L.; Sanchez, S.; Schmidt, O.G. Propulsion mechanism of catalytic microjet engines. IEEE Trans. Robot. 2014, 30, 40–48. [Google Scholar] [CrossRef]
  21. Gibbs, J.G.; Zhao, Y.P. Autonomously motile catalytic nanomotors by bubble propulsion. Appl. Phys. Lett. 2009, 94, 759. [Google Scholar] [CrossRef]
  22. Wang, Z.; Chi, Q.; Liu, L.; Liu, Q.; Bai, T.; Wang, Q. A viscosity-based model for bubble-propelled catalytic micromotors. Micromachines 2017, 8, 198. [Google Scholar] [CrossRef] [PubMed]
  23. Wang, Z.; Chi, Q.; Bai, T.; Wang, Q.; Liu, L. A dynamic model of drag force for catalytic micromotors based on navier–stokes equations. Micromachines 2018, 9, 459. [Google Scholar] [CrossRef] [PubMed]
  24. Liu, L.; Bai, T.; Chi, Q.; Wang, Z.; Xu, S.; Liu, Q.; Wang, Q. How to make a fast, efficient bubble-driven micromotor: A mechanical view. Micromachines 2017, 8, 267. [Google Scholar] [CrossRef] [PubMed]
  25. Chi, Q.; Wang, Z.; Tian, F.; You, J.A.; Xu, S. A review of fast bubble-driven micromotors powered by biocompatible fuel: Low-concentration fuel, bioactive fluid and enzyme. Micromachines 2018, 9, 537. [Google Scholar] [CrossRef]
  26. Blander, M.; Katz, J.L. Bubble nucleation in liquids. Adv. Colloid Interface Sci. 1975, 10, 1–32. [Google Scholar] [CrossRef]
  27. Thorncroft, G.E.; Klausner, J.F.; Mei, R. An experimental investigation of bubble growth and detachment in vertical upflow and downflow boiling. Int. J. Heat Mass Transf. 1998, 41, 3857–3871. [Google Scholar] [CrossRef]
  28. Staniszewski, B.E. Nucleate Boiling Bubble Growth and Departure; Massachusetts Institute of Technology: Cambridge, MA, USA, 2011. [Google Scholar]
  29. Ishida, N.; Inoue, T.; Minoru Miyahara, A.; Higashitani, K. Nano bubbles on a hydrophobic surface in water observed by tapping-mode atomic force microscopy. Langmuir 2000, 16, 6377–6380. [Google Scholar] [CrossRef]
  30. Mei, R.; Klausner, J.F. Unsteady force on a spherical bubble at finite reynolds number with small fluctuations in the free: Tream velocity. Phys. Fluids A Fluid Dyn. 1991, 4, 613–631. [Google Scholar]
  31. Happel, J.; Brenner, H. Low Reynolds Number Hydrodynamics, with Special Applications to Particulate Media; Springer Science & Business Media: Berlin, Germany, 1973. [Google Scholar]
  32. Klausner, J.F.; Mei, R.; Bernhard, D.M.; Zeng, L.Z. Vapor bubble departure in forced convection boiling. Int. J. Heat Mass Transf. 1993, 36, 651–662. [Google Scholar] [CrossRef]
  33. Saffman, P.G. The lift on a small sphere in a slow shear flow. J. Fluid Mech. 2006, 22, 385–400. [Google Scholar] [CrossRef]
  34. Wang, W.; Chiang, T.Y.; Velegol, D.; Mallouk, T.E. Understanding the efficiency of autonomous nano- and microscale motors. J. Am. Chem. Soc. 2013, 135, 10557–10565. [Google Scholar] [CrossRef] [PubMed]
  35. Orozco, J.; Jurado-Sanchez, B.; Wagner, G.; Gao, W.; Vazquez-Duhalt, R.; Sattayasamitsathit, S.; Galarnyk, M.; Cortes, A.; Saintillan, D.; Wang, J. Bubble-propelled micromotors for enhanced transport of passive tracers. Langmuir 2014, 30, 5082–5087. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Gas-liquid–solid interface and contact angle of the bubble.
Figure 1. Gas-liquid–solid interface and contact angle of the bubble.
Micromachines 10 00415 g001
Figure 2. A schematic diagram of the bubble in the nucleation stage.
Figure 2. A schematic diagram of the bubble in the nucleation stage.
Micromachines 10 00415 g002
Figure 3. Various forces during the growth stage of the bubble.
Figure 3. Various forces during the growth stage of the bubble.
Micromachines 10 00415 g003
Figure 4. The contact of the bubble with the inner wall of the micromotor. (a) The volume of several parts as the bubble is in full contact with the inner wall of the micromotor. (b) The geometry of the bubble and the forces it experiences.
Figure 4. The contact of the bubble with the inner wall of the micromotor. (a) The volume of several parts as the bubble is in full contact with the inner wall of the micromotor. (b) The geometry of the bubble and the forces it experiences.
Micromachines 10 00415 g004
Figure 5. A schematic diagram of the ejection stage of the bubble at the opening of the micromotor.
Figure 5. A schematic diagram of the ejection stage of the bubble at the opening of the micromotor.
Micromachines 10 00415 g005
Figure 6. Various forces in the stage of bubble ejection of a conical micromotor.
Figure 6. Various forces in the stage of bubble ejection of a conical micromotor.
Micromachines 10 00415 g006
Figure 7. Relationship between nucleation free energy change and the nucleus radius of the bubble.
Figure 7. Relationship between nucleation free energy change and the nucleus radius of the bubble.
Micromachines 10 00415 g007
Figure 8. Relationship between driving force and (a) surface tension coefficient (b) fluid viscosity.
Figure 8. Relationship between driving force and (a) surface tension coefficient (b) fluid viscosity.
Micromachines 10 00415 g008
Figure 9. Relationship between driving force and (a) concentration of hydrogen peroxide CH2O2 and (b) opening radius.
Figure 9. Relationship between driving force and (a) concentration of hydrogen peroxide CH2O2 and (b) opening radius.
Micromachines 10 00415 g009
Figure 10. Relationship between driving force and (a) micromotor length (b) tangent semi-cone angle.
Figure 10. Relationship between driving force and (a) micromotor length (b) tangent semi-cone angle.
Micromachines 10 00415 g010
Figure 11. A phase diagram of bubbles in the stage of (a) bubble nucleation (b) bubble growth (c) bubble slip (d) bubble ejection, where red represents the gas and blue represents the fluid.
Figure 11. A phase diagram of bubbles in the stage of (a) bubble nucleation (b) bubble growth (c) bubble slip (d) bubble ejection, where red represents the gas and blue represents the fluid.
Micromachines 10 00415 g011
Figure 12. Tubular micromotor (a) pressure distribution and (b) streamline.
Figure 12. Tubular micromotor (a) pressure distribution and (b) streamline.
Micromachines 10 00415 g012
Figure 13. The velocity of the micromotor in the full life-cycles of bubbles, namely bubble nucleation, growth, slip, and ejection. Inset: ejection stage.
Figure 13. The velocity of the micromotor in the full life-cycles of bubbles, namely bubble nucleation, growth, slip, and ejection. Inset: ejection stage.
Micromachines 10 00415 g013
Figure 14. Optical microscopy and scanning electron microscopy of tubular micromotor.
Figure 14. Optical microscopy and scanning electron microscopy of tubular micromotor.
Micromachines 10 00415 g014
Figure 15. The velocity of the tubular micromotor in the stage of bubble nucleation, slip, and ejection.
Figure 15. The velocity of the tubular micromotor in the stage of bubble nucleation, slip, and ejection.
Micromachines 10 00415 g015
Table 1. The driving forces, average speed, and duration of the bubble-driven micromotor at different stages.
Table 1. The driving forces, average speed, and duration of the bubble-driven micromotor at different stages.
StagesEquationDriving ForceAverage Speed (μm/s)Duration (ms)
Nucleation(14)0.039586.808 × 10−40–4.2
Growth(56)299510114.2–5.25
Slip(33)18.83191.45.25–8.97
Ejection(49)236015.658.97–8.98
Table 2. Drag forces of the various surfaces of the tubular motor.
Table 2. Drag forces of the various surfaces of the tubular motor.
Surface AreaUpstream RingBackflow RingOuter WallInner Wall
Drag force (nN)0.0280.0430.0020.191
Drag force (nN)0.264

Share and Cite

MDPI and ACS Style

Lin, Y.; Geng, X.; Chi, Q.; Wang, C.; Wang, Z. Driving Forces of the Bubble-Driven Tubular Micromotor Based on the Full Life-Cycle of the Bubble. Micromachines 2019, 10, 415. https://doi.org/10.3390/mi10060415

AMA Style

Lin Y, Geng X, Chi Q, Wang C, Wang Z. Driving Forces of the Bubble-Driven Tubular Micromotor Based on the Full Life-Cycle of the Bubble. Micromachines. 2019; 10(6):415. https://doi.org/10.3390/mi10060415

Chicago/Turabian Style

Lin, Yongshui, Xinge Geng, Qingjia Chi, Chunli Wang, and Zhen Wang. 2019. "Driving Forces of the Bubble-Driven Tubular Micromotor Based on the Full Life-Cycle of the Bubble" Micromachines 10, no. 6: 415. https://doi.org/10.3390/mi10060415

APA Style

Lin, Y., Geng, X., Chi, Q., Wang, C., & Wang, Z. (2019). Driving Forces of the Bubble-Driven Tubular Micromotor Based on the Full Life-Cycle of the Bubble. Micromachines, 10(6), 415. https://doi.org/10.3390/mi10060415

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