Research Article  Open Access
M. I. Khodabocus, M. Sellier, V. Nock, "Slug SelfPropulsion in a Capillary Tube Mathematical Modeling and Numerical Simulation", Advances in Mathematical Physics, vol. 2016, Article ID 1234642, 16 pages, 2016. https://doi.org/10.1155/2016/1234642
Slug SelfPropulsion in a Capillary Tube Mathematical Modeling and Numerical Simulation
Abstract
A composite droplet made of two miscible fluids in a narrow tube generally moves under the action of capillarity until complete mixture is attained. This physical situation is analysed here on a combined theoretical and numerical analysis. The mathematical framework consists of the twophase flow phasefield equation set, an advectiondiffusion chemical concentration equation, and closure relationships relating the surface tensions to the chemical concentration. The numerical framework is composed of the COMSOL Laminar twophase flow phasefield method coupled with an advectiondiffusion chemical concentration equation. Through transient studies, we show that the penetrating length of the bidroplet system into the capillary tube is linear at earlytime regime and exponential at latetime regime. Through parametric studies, we show that the rate of penetration of the bidroplet system into the capillary tube is proportional to a timedependent exponential function. We also show that this speed obeys the Poiseuille law at the earlytime regime. A series of position, speedversusproperty graphs are included to support the analysis. Finally, the overall results are contrasted with available experimental data, grouped together to settle a general mathematical description of the phenomenon, and explained and concluded on this basis.
1. Introduction
When two physically distinct but miscible liquid droplets are placed in a capillary tube, a key feature that is commonly observed is the bulk motion of the composite droplet under the action of capillarity until complete mixing is attained [1]. In Figure 1 is given a freehand sketch of this scenario. This is an example of motion due to (chemical) Marangoni effect [2–4]. The motion of the bidroplet system is caused by local variations of its advancing and receding surface tensions that are in turn caused by differences in its chemical compositions. Out of this simple scenario grew two kinematical concepts, namely, the penetrating length and the penetrating rate. When both concepts are grouped together, we attain, by inference, slug selfpropulsion mechanism—a juxtaposition of two miscible liquid droplets. An important facet of this driving mechanism is that it is simple and so could be delicately adjusted to manipulate liquid droplets in smallscale fluidic devices. Earlier scientists have studied the phenomenon experimentally and have concluded that it might be used as a means of selfpropelling, transporting microscale substances at a respectable penetrating length and rate, respectively, of orders of centimeters, centimeters per second [1, 5, 6].
Bouasse [7] inquired experimentally on capillarity and interfacial phenomena. He found that a droplet predeposited in a conic tube undergoes motion towards the narrow region because of the Laplacian pressure difference induced on both its menisci. Weislogel [8] investigated the dynamics of a liquid droplet in a chemically asymmetric tube, by treating the tube to be successively hydrophilic and hydrophobic. He found the droplet to move towards the hydrophilic side. Dos Santos and Ondarçuhu [9] argued that a droplet containing a reactive substance remained in its state of motion as long as such reactive agent is operative—the dynamics of reactive wetting on solid surfaces were discussed by De Gennes [10].
Historically, it may be that the first documented study goes back to 1712, when Hauksbee [11] investigated the behaviour of a drop of oil of oranges sandwiched in between two nonparallel glass plates. In that paper, two points obviously follow in the light of this phenomenon. First, the fact that the motion of the drop reversed when the slope of those two plates was reversed led Hauksbee to believe that the direction of motion is towards the narrow region. Secondly, the fact that his experiment was reproducible in vacuum further informed Hauksbee that such motion is solely due to interfacial phenomenon.
Currently, existing (active) techniques of manipulating discrete droplets include electrowetting [12–17], dielectrophoresis [18], thermocapillary [19–22], surface acoustic wave (SAW) [23, 24], and pressuredriven basedchannel systems [25]. Thanks to these sophisticated techniques, the fields of nano and microfluidics have seen over the past decade considerable development.
Selfpropulsion phenomena as well as coalescence phenomena are found to have generated considerable interest and many features have been correctly described and explained as well. Some fundamental questions were addressed by De Gennes [10, 26, 27] and Shanahan and De Gennes [28]. Shanahan and De Gennes stressed that the motion of the system is switched on mainly due to the slight perturbation communicated to the system at the time when the second droplet is brought into close contact to the predeposited one. Furthermore, the authors pointed out that for motion to occur the system must necessarily be asymmetric and that such motion is maintained only insofar that such character is observed in the system.
Recently, there has been a growing interest in exploiting droplets selfpropulsion phenomenon, thereby avoiding the arduous fabrication procedures required [6, 29, 30]. For instance, Karpitschka and Riegler [31] gave an account on the coalescence of liquid droplets in two respective regimes, namely, the fast coalescence or delayed coalescence. In the paper of Karpitschka and Riegler, the authors pointed out that the coalescence phenomenon is predominantly governed by surface tension, regardless of whether it is fast or delayed coalescence. For instance, when the difference in the surface tensions of the working droplets is such that mNm^{−1}, delayed coalescence phenomenon is observed from the bidroplet system and, for all other cases such that mNm^{−1}, fast coalescence phenomenon is observed. In a subsequent paper, Karpitschka and Riegler [32] explore the study of delayed droplets coalescence phenomenon on a combined experimental and theoretical basis.
On the one hand, the selfpropulsion of liquid droplets has also been addressed by Sinz and Darhuber [6] and Sellier et al. [29, 30]. In the works of Sellier et al., the authors examined the selfpropulsion mechanisms in unconfined systems, that is, sessile droplets configuration. The principles of manipulating discrete droplets using active methods were considered by Darhuber, Valentino, Davis, and Troian [33, 34]. In the paper of Sinz and Darhuber, the authors addressed the selfpropelling of surfactant droplets in chemically confined systems. Therein, Sinz and Darhuber argued that this passive driving mechanism might be used to transport microscale substances.
Lunati and Or [35] undertook the problem of gravitydriven slug dynamics in capillary tubes. The authors have answered questions regarding the origin, strength, and effects of the dragforces developed at the contact lines. They also studied the kinematics of the slug, contrasting the numerical simulation results with the experiment results of Bico and Quéré [36] on falling of viscous slugs in vertical capillary tubes. Bico and Quéré deduced a set of scaling relationships characterising the steady motion of the slug. On the experimental side, the earliest laboratory experiments wherein the subject matter based on purely confined systems has been explored thoroughly were those of Bico and Quéré [1, 5]. These are the most important works in close connection with the subject at hand. In these, more work was endeavoured on the dynamics of trains of droplets in horizontal capillary tubes in both wetting and drying conditions. For instance, in the former paper of Bico and Quéré [5], which constitute their preliminary groundbased experiments, they argued mainly on the coating properties of trains of juxtaposed liquid droplets in terms of geometrical and viscous effects—another very interesting field of inquiries; they moreover drew attention on their versatility towards integrated microfluidic devices. For example, Bico and Quéré [1] have studied the response of the system in more complex geometries and displayed more features on a combined experimental and theoretical basis.
Since then, only few scientists have begun to study the subject on theoretical and numerical sides. Furthermore, studies numerically based on the slug selfpropulsion phenomenon do not appear to have received much attention. To sum up, it may be concluded that we found no scientist has come to ponder on the selfpropelling phenomenon in confined system quantitatively and to arrive at a more rational conception of it.
We have studied carefully the dynamics of a two miscible fluids in a capillary tube and, accordingly, found the law describing the position (or penetrating length) of the slug into the capillary tube is well described by the following relationship: , where , are propertydependent constants. Furthermore, it is observed at the earlytime that the slug velocity (or rate of penetration), that is, , obeys a Poiseuille law: . In this expression, designates the radius of the capillary tube; , are the surface tensions of the miscible fluids; and are the dynamic viscosity and length of the slug. We came to these laws numerically by collecting data on the slug position as a function time and the maximum travelled distance. An inspection of the form of the property (discussed hereinafter) evidences at once that small and low viscous substances are speedier than large and highly viscous ones.
In this paper, an account of the mathematical development is given. The paper is organised as follows: Section 2 is devoted to mathematical formulations and Section 3 to the numerical implementation of the mathematical framework. Results from numerical experiments and analysis are presented in Section 4: the motions of the slug with respect to time and property constants are discussed in Sections 6 and 7, respectively; the results are contrasted with the papers of Bico and Quéré [1, 5]. In Section 8, a powerlaw model describing the position of the slug with respect to time and the leading property constants involved is established. The overall results are finally argued in Section 9.
2. Mathematical Formulations
2.1. Preliminaries
First, some aspects, definitions, and notations are in order. Truly speaking, we examine here a threefluid flow problem, the properties of which are such that two of which are miscible phases but immiscible with the third one. In this context, their respective domains of definition are written as where is a point of either the set or . The scalar fields (with ) are the respective timedependent positions of the (front and rear) interfaces of the slug. Obviously, the sets and denote those timedependent interfaces. These physical domains, sets (1), are clearly illustrated in Figure 2. The choice of Figure 2 is motivated by that of Figure 1. In view of the foregoing definitions, it thus appears evident to regard merely as a twomiscibleliquid region and the (disjoint) union of its rear and front vapour phase regions. In fact, with these underlying definitions in mind, both phases can be related without loss of generality as thus where denotes the volume of its arguments. Furthermore, we characterise by the scalar fields the functionals describing the volume fraction of fluids in and , respectively. Consequently, by virtue of the above postulation, (2), it results inwhere (with ). In regard to these facts, the problem under consideration now reads an Immiscible TwoPhase Flow problem. At this point, we have made the essential geometrical features prominent and, in what follows, a twophase flow concentrationdependent surface tension problem of a viscous liquid slug in a capillary tube is undertaken.
2.2. Transport of Concentration by Diffusion and Convection
Denote by the scalar field the chemical concentration of a (highly) volatile substance filling the region initially (note that the other substance filling the region is taken to be lowly volatile). To establish the transport equation for the chemical concentration function in , we reason in the following way.
Initially, postulate the concentration within the slug to be such that and presuppose a normally imposed flux on the corresponding boundary (with ) by its external surroundings. Then, by virtue of these assertions, is simply described by the equation of chemical species transported by diffusion and convection (CDmodel, say), asfor given . In (5), the proportionality factor is the diffusivity of fluid 1 into ; the vector field denotes the outward normal unit vector corresponding to the boundary (with ). It may be pointed out in passing that this convectiondiffusion model, that is, (5), merits further attention, for it should also agree with the subsequent changes in position of the slug. This point is clarified in a later section.
2.3. Inclusion of ConcentrationDependent Surface Tension
We now discuss proper forms of interfacial forces at the boundaries , . Supposing a solution for the scalar field is granted, we can then define as the growth, (resp., decay, ) of interfacial effects in terms of at their respective boundaries , . In order to deduce such proper forms for , we shall in the first place consider the following facts: initially, the domains , are characterised to be such that in and in , and, finally, both are expected to be in equilibrium (chemically), attaining the value , say; thence, the correct conditions to be imposed upon (with , ) simply write asTaking (6) into consideration along with an additional assertion that the surface tensions evolve linearly, one arrives atThe linearity of (with ) on the scalar field is here motivated by the works of Carles and Cazabat [37, 38] on the spreading dynamics of Newtonian droplets. Therein, Carles and Cazabat [37] found that the surface tension of those mixed substances varied roughly linearly with the chemical concentration of the mixture. Equation (7) is physically reasonable, for in and in . However, out of (7) grew one important question which is somehow unclear.
Question 1. What are the magnitudes of , (with )?
The answer of this question is found in the following line of reasoning. First, it is evident thatwhere the geometrical property stands for the size of along the axis (see Figure 2). That is, the (chemical) equilibrium state is somewhat dependent upon the geometrical properties of the slug, as found in the paper of Bico and Quéré [1]. Thus, the expression which would match these conditions is , since in (with ) at time . Note that is also an admissible value but only for a particular case so to speak when , for it neglects the geometrical property of the slug, particularly, when , .
In regard to these arguments, one then concludes that the lesser the length the higher the concentration . It now remains to settle proper forms for (with ). First, we note thatHence, these suggest to consider a barycentricbased model, namely, of the form because (equilibrium surface tension) when . Clearly, there is no loss of generality in doing so, for not only does (10) agree well but it is also a physically acceptable expression. We are now well prepared to set forth the governing equations; this is the purpose of the next section.
2.4. TwoPhase Flow PhaseField Model
In Sections 2.2 and 2.3, we covered much ground on diffusive and interfacial effects. The present section takes those effects into account and models the slug selfpropulsion mechanism using the phasefield approach. In a sense, we have beforehand the motion of a twomiscibleliquid system partly in contact with an immiscible one, the vapour phase region. Therefore, the starting point is naturally based on the classical NavierStokes equations (NSmodel, say) for incompressible flow with interfacial forces coupled with the socalled CahnHilliard evolution equation (CHmodel, say). In the absence of gravitational field, those equations in vector invariant forms read as follows: where the stress tensor , constructed with the pressure and velocity fields , reads as follows:It is understood that (). In (11) and (12), the property constants characterise the mixed density and dynamic viscosity of the slug; both are defined in terms of the material properties of the two phases, and the smooth (phasefield) function is (with , ). The function defines the socalled bulk chemical potential asIn (13), is a capillary width that scales with the thickness of the boundaries , ; [m^{3}skg^{−1}], the mobility term which determines the time scale of the CahnHilliard diffusion; and , the mixing energy density. For more details see, for instance, the papers of Yue et al. [39] and Jacqmin [40] as well. We draw attention to the fact that , and, thereforeIn (14), and are operators which give the minimum and maximum of their respective arguments. Throughout the following discussion, the phasefield function is used instead of . It now remains to supplement the appropriate boundary conditions on . At the inlet and outlet boundaries a zero pressure is imposed; that is, the system containing the domains , is embedded in a still environment . Wetted wall boundary conditions are imposed along ; the contact angles are such that ( [rad]). The choice [rad] (with ) is justified from the fact that the advancing and receding boundaries of the miscible bidroplet system are convex curves, following the works of Bico and Quéré [1, 5]. More features based partly on the field are discussed in Section 3.
3. Numerical Implementation
This section furnishes the mathematical model so far discussed with numerical features to simulate the motion of the slug. Figure 3 shows a simplified computational box model with the corresponding zero level set function delineating the regions of interest.
Recall that we are concerned with the motion of a slug so in order to adapt the solution of with its motion one should remap certain property variables attached to on a phasefield basis to track the evolution of this function while following the region . In other words, this suggests to remap the scalar fields , as follows: , (with ). Moreover, since , which is equivalent to say that , the remapping is first taken into account so as to obtain , as can be easily checked by substitution. Clearly, from these remappings one observes that the convectiondiffusion equation vanishes only outside the region . Consequently, these now result in a model which spans over , while retaining its characteristics only in . Hence, invoking these mapping properties into (5) yields where the fields , are for shortness defined byFollowing (15), the problem then becomes one of finding the solution of a General Form PDE. Alternatively, one can formulate the convectiondiffusion equation describing the scalar field as follows:where the reaction term appearing in the above transport reads asThis is the problem of finding the solution of a Transport PDE. In COMSOL, the latter is called Transport of Dilute Species Interface. It may be observed in passing the fact that whenever is a point of , implies that and , thence the original convectiondiffusion equation. Contrarily, when is not a point of , implies that and , thence the restriction of the property within the domain . These latter two statements confirm very clearly the validity of the above underlying mapping procedures.
Now, some techniques which have rendered efficient aid towards implementing this model with COMSOL are in order. Henceforth, we presume the following COMSOL interfaces are activated: the phasefield and General Form PDE interfaces, respectively: the former for the NS and CHmodels and the latter for the CDmodel. Decidedly, it is observed that the conditions imposed on imply thatthence, it follows from (19), using the and operators, thatWith these features alongside, we first utilise the COMSOL and operators to construct the interfacial forces at their respective boundaries . Then, these and other fluid properties are entered into the fluid properties (with ) features. Finally, , are selected as initial interfaces. In regard to the General Form PDE interface, (15), its coefficients and flux components, namely, (with ), are entered into the appropriate default nodes, the concentration in is initialized to unity, and a zero flux boundary condition is assumed from the surroundings. This being so, several simulations were then run with an array of data extracted from the papers of Bico and Quéré [1, 5]. Extracted data from those authors are tabulated in Table 1. These simulations are discussed through several cases in the following sections.

4. Results and Discussion
In this section we analyse the results of our simulations based on the data given in Table 1. The data appearing in Table 1 are those extracted from the papers of Bico and Quéré [1, 5]. This made the need of a condensed notation imperative. Define in compact form those property constants, namely, , upon which the numerical experiments were carried out as , and henceforward view the surface tension difference between both menisci as a function of the set of parameters and of time , respectively; that is,Let the position and speed/velocity, that is, and , respectively, of the slug be functions of the formsso that the corresponding differentials with respect to time write as From now on, we can then characterise as the instantaneous rate of change but , the net change in the property with respect to interfacial variation and , the net rate of change in interfacial forces with respect to time. In particular, by is meant the time derivative of following the dynamics of , where , . We will come back to those terms in later sections. In other words, these bring forward the key factors that do have significant effects on the evolution of the scalar fields , . Suppose the following baseline properties are granted: , and important questions related to those properties might then be addressed. For instance, it would be fruitful to extend theoretical understanding based on the following inequalities: Henceforth, we primarily aim at answering two sets of questions. In the first place, we aim to answer the following questions.
Question 2 (given ). What are the kinematics of the slug in steady and unsteady regimes?
Question 3 (given ). How does the slug evolve when is perturbed to ?
And, in the second place, we aim to answer the following question.
Question 4 (given and , resp.). Are the positions of the slug predictable for all ?
In the above last question, the scalar field denotes the position of the miscible fluids system after selfpropulsion and the corresponding time. Thus, by the pairs of scalars , are meant the ending points of its path in the plane. In the following sections we shed lights on the above three questions by advancing a series of graphical results and mechanical explanations.
5. The PhaseField Function
The proposed methodological approach rests on the phasefield technique. Knowing the subsequent positions of the phasefield function when time evolves, one can take advantage if it calculates the subsequent position of the slug. Thus, the discussion commences by computing the scalar field . To do so, we use the set of property constants tabulated in Table 1.
The time evolution of the phase function is illustrated at the sequence of times [s] in Figure 4. We draw attention that is graphed by the remapping in order to have instead of .
6. Case 1: Motion of the Slug versus Time
The following two cases concern the rate of penetration of the slug into the capillary tube in steady and unsteady regimes. The steady case/regime refers to the case for which the positiontime graph of the slug is linear and the unsteady case/regime to other cases.
6.1. Case 1.1: Steady Regime
Generally, on a Poiseuille flow hypothesis together with the assertion that the slug is sufficiently long, say, its path results in an integral expression of the formTo settle (25), we have taken into account the following statements [1]: the Poiseuille flow induced within the region is of the form , where the Laplacian pressure exerted at relative to that exerted at reads as . The velocity field is averaged in , and the resulting expression equated to , yielding , since , thence the antiderivative of (25).
Consequently, it follows that there exists a functional such that for given , . Numerical simulation was carried out for the steady case in the hope of establishing the character of . The earlytime results confirm that , and the position of the slug obeys the following linear law at earlytime:where [s], as can be easily seen in Figure 5. To graph Figure 5, the COMSOL Cut Line data set is used to create a line parallel to the axis of the capillary tube, having its ending points at , , where denotes the length of the capillary tube along the axis. This allows to track the subsequent positions occupied by the advancing and receding boundaries, that is, , , since it will cut at every instant of time these boundaries at precisely two and only two (boundary) points (with ). Setting , we then construct a timeposition sequence . Representing these data points by circular markers (), the corresponding approximation by solid line, and that of Bico and Quéré [5] by dotted line, one obtains Figure 5. Henceforth, whenever the study is carried out with respect to time, it is to be understood that the plotting referred to in the plane is based on that approach, unless the study is carried out parametrically, in which case we shall then draw attention upon.
In the sequel, we will show through numerical experiments that . Thus, by (26) it would appear that after some time , the slug will have travelled at a constant velocity a distance . Clearly, (26) appears to show that whenever , a situation which is reasonable mathematically but unreasonable physically; this is simply because the mathematical description using the linear approximation will at some time (unsteady case) fail to give physical agreement. In that case a nonlinear law (as will be shown in the sequel to exist) should be sought.
The above statements suggest that there exists a scalar , characterising the distance at which the slug stops, such that , since we know from the nature of the problem that the difference in interfacial forces, which was initially very high, that is, (), owing to the initial conditions in (with ), will have decreased to a vanishingly small value when the equilibrium of has been established. Thus, this demonstrates the existence of a temporal (dependant) property such that mechanical equilibrium is attained via interfacial effects; namely,Clearly, there is no loss of generality in assigning, for instance, this particular case to an initially prescribed set, say , where . Based on the said facts, we ran the numerical simulation with a large final time [s] to demonstrate that (27) is satisfied.
Surprisingly, we observed with a slug bearing the characteristics [mm], [mPas] (with ), and [m^{2}s]—the diffusion constant for alcohol fluid 1 in water fluid 2 [41, 42]—thatThe graphical evolution of the surface tension function with respect to time is shown in Figure 6. When the surface tensions () are plotted on the same plane, one obtains Figure 7, where curved arrows appearing in that plane denote the direction of convergence of these functions with respect to time.
A very interesting feature is observed from Figure 7. We have [Nm^{−1}] at time [s]. Substituting the values tabulated in Table 1 into the barycentricmodel equation (10) for the pairs of scalar fields , , one arrives at [Nm^{−1}] (with ), thereby validating the form of (10).
Out of the graphical results displayed in Figures 6 and 7, it would appear that (26) does not hold in both the early and latetime regimes, more precisely, in the steady (or earlytime) regime as well as in the unsteady (or latetime) regime. The next section discusses the dynamics of the slug in the unsteady regime.
6.2. Case 1.2: Unsteady Regime
When the period over which the numerical experiments were extended to [s], the subsequent positions of the slug graphed exponentially, as can be easily shown in Figure 8. Hence, the path of the slug could fairly well be described by an exponentialtype equation of the formIn (29), the fields , are dependent scalars; that is, and . Furthermore, since the slug moves steadily only for a short period of time (roughly, 10 [s]), consequently, there does necessarily exist , such that, for given ,where . Here, it is to be understood that is the place at which the slug stops whereas simply characterises its steady regime. Also, stand, respectively, for the initial time [s], the time [s], and the time [s] at which the slug stops. But designates the period over which the numerical experiments are carried out. With these characteristics at hand, the path of the slug then writes asFrom this, one concludes that the position of the slug is merely a superposition arising from two contributions: the steady and unsteady regimes; the former is characterised by and the latter by , respectively. Following Figure 5, it is obvious that the equality holds only for [s]. Thus, for early times, the position of the slug can be approximated aswhich is precisely the same expression in the steady case, being put for . It now remains to demonstrate that () coincides with the result of Bico and Quéré [5].
For two typical slugs of lengths , [mm], respectively, we combined (29) with numerical data in the leastsquares sense. We found for [mm] and for [mm]. Thence, with the former, we deduce [mms] and with the latter [mms^{−1}]. All these quantities are arranged carefully in Table 2. Figure 8 is intended to exhibit these features graphically using (31).

No doubt, it is no exaggeration to say that our model, verified numerically in those two cases, namely, for , [mm], as depicted in Figure 8, may be extended to analyse the movement of small volumes of viscous liquids in capillary tubes of square, rectangular cross sections and so forth. Hence, given the set , that pair of scalar fields designating the point at which the slug stops for an initially given speed, one can predict a priori the position as well as the rate of penetration of the slug into the capillary tube. Further, if the focus is on rapid motion, one should evidently consider a small slug, instead of a big one. Contrarily, if the primary goal is to achieve large penetrating length, one might endeavour to take a large slug; the trueness of these two statements is reflected in Figure 8. The following sections enable further physical relationships to be grasped.
7. Case 2: Motion of the Slug versus Geometrical and Physical Properties
The position of the slug is such that . So far, we studied this position function with respect to time, keeping the elements of the parameter vector fixed at those values which are tabulated in Table 1. Here it is that this function is studied parametrically. We agree that parametric studies are best undertaken when the property in question is timeinvariant. We have, on the one hand, and ; see (30). And, on the other hand, the slug moves according to Poiseuille flow only in the earlytime regime. As a consequence, it seems, therefore, legitimate to carry out the parametric studies in the earlytime regime, which are what we proposed to consider hereinafter.
Let designate the dynamic viscosity of the slug (mixture of the two miscible liquid droplets). On a barycentricbased basis we may then consider , in which case the velocity of the slug, following (29), writes asWe understand that (33) is the derivative of (25) with respect to ; that is, . If (33) holds true, so are the following proportionalities: , , and . These are what we are going to show through a series of numerical experiments in the following sections.
7.1. Case 2.1: Motion of the Slug versus the Radius of the Capillary Tube
In this section the motion of the slug is studied in a series of capillary tubes having the following radii: [mm] (with ). The properties of the slug are unaltered. We aim at establishing an algebraic expression relating the velocity of the slug to the radius of the capillary tube. For every , a numerical experiment is carried out and the velocity of the slug is calculated. When, for each , the ratio is analysed, it was found thatfrom which we conclude that . The corresponding pairs of values are plotted on a plane; see Figure 9. Consequently, the dependence of the functional finally restricted to evolves linearly with respect to the radius of the capillary tube according to the following linear relationship:It is understood that the scalar fields , are partially dependant constants (see the character of (33) to understand this point). In other words, an increase in by results in a linear increase in , conversely, which manifests the nature of the size of the capillary tube. The next section examines the influence of the geometrical property on the kinematics of the slug.
7.2. Case 2.2: Motion of the Slug versus Its Length
We herein explore the effect of the geometrical property (size of the slug) on the slug motion. We ran for the sequence [mm] (with , ) several numerical experiments, keeping . The accompanying table (Table 3) compares the typical parameters for the slug selfpropulsion and its laboratory analogs extracted from the paper of Bico and Quéré [5].

Figure 10 gives some graphical illustrations on this basis. On a careful perusal of Table 3 and Figure 10, respectively, one concludes that the agreement between simulations and experiments is good. As a result, the present study puts forward an important point: an increase in (by [mm], say) weakens the surface tension difference, that is, . Hence, one concludes that a slug of length moves faster than the one whose length slightly differs, for instance, by . Furthermore, perusing Figure 10, we conclude that the dependence of the functional restricted to on evolves on a loglog scale according toAs pointed out so far, it is understood that the scalar fields , are partially dependant constants. From Figures 10 and 11, it may be remarked in passing thatConverting every loglog sizespeed pair into the usual pair and, further, letting , be the corresponding numerical, experimental values, the tabulation of all these property constants yields Table 4.

7.3. Case 2.3: Motion of the Slug versus Diffusivity
This section discusses the effects of diffusivity constant . A way of interpreting the latter may be understood on the following basis. At leading order, we deduced from the constitutive equation of the following scaling (diffusion length): , which gave insight of how far the concentration has propagated inside the slug in time for a given diffusivity constant . In Figure 12 is graphed on a loglog scale the speed of the slug versus the diffusivity constant . When every loglog diffusivityspeed pair is converted into the usual pair and, further, to stress that the dependence of the speed of the slug on the diffusivity has not been studied numerically, the tabulation of these property constants yields Table 5.
