When two Newtonian liquid droplets are brought into contact on a solid substrate, a highly curved meniscus neck is established between the two which transforms the bihemispherically shaped fluid domain to a hemispherically shaped domain. The rate at which such topological transformation, called coalescence phenomenon, evolves results from a competition between the inertial force which resists the transformation, the interfacial force which promotes the rate, and the viscous force which arrests it. Depending on the behaviour of these forces, different scaling laws describing the neck growth can be observed, predicted theoretically, and proved numerically. The twofold objective of the present contribution is to propose a simple theoretical framework which leads to an Ordinary Differential Equation, the solution of which predicts the different scaling laws in various limits, and to validate these theoretical predictions numerically by modelling the phenomenon in the commercial Finite Element software COMSOL Multiphysics.

1. Introduction

When two Newtonian liquid droplets are brought into contact on a solid substrate, a highly curved meniscus bridge is established and thereafter driven by surface tension so as to transform the bulk flow region from bihemispherically shaped at early-time regimes to more hemispherically shaped at late-time regimes. Such topological transformation, involving the merging of two or more sessile droplets into a single one on a solid surface, is called droplet coalescence phenomenon.

Physically speaking, the coalescence process consists of three regimes, called the viscous, visco-inertia, and inertial coalescence regimes [1]. In each regime, a range of surprisingly complex behaviours can be observed, so to speak topological changes of the interface [2]. Thus, the subject of droplets coalescence is a fascinating, multifaceted subject of inquiry with many interesting phenomena to learn and discover. Indeed, the subject has intrigued many physicists [36] and mathematicians [7, 8].

In actual fact, authors of [9] appear to be the first scientists who have addressed the subject. The authors study the formation of vortex rings by droplets falling into liquids and some allied phenomena. For a miscible droplet falling from not too great a height into its bulk fluid, it is observed that the droplet descends through the bulk fluid as a ring, whereas for the case of an immiscible droplet, [9] observed the droplet to descend through the bulk fluid as a spherical droplet. Thereafter, there have been many attempts on experimental, theoretical, and numerical sides to elucidate the underlying physics governing the dynamics of coalescence.

Authors of [10] have studied the partial coalescence of a water droplet, with radius mm, with the free-surface of a bulk fluid composed of water-glycerol mixtures experimentally. For cases , mm, the authors drew attention that the coalescence process is arrested by gravitational and viscous effects and, depending on these two effects, three coalescence regimes were observed, which [10] termed gravity, inertio-capillary, and viscous regimes. Furthermore, the authors argued that the coalescence time is strongly dependent upon the Bond number and the Ohnesorge number , where , , and denote the density, dynamic viscosity, and surface tension of the droplet, respectively; and denote the properties of its vapour phase surrounding, and denotes the gravitational acceleration.

Reference [7] addressed the more general case of the coalescence of droplets surrounded by a viscous outer fluid, both analytically and numerically, using asymptotic methods and the integral equation given by [11]. Describing the initial coalescence regime by the Stokes equations, they showed that the three-dimensional solution is asymptotically equivalent to the two-dimensional one reported in the works of [1216]. When the viscosity of the outer fluid is disregarded, the radius and width of the highly curved meniscus bridge were found to obey the following laws: and , respectively, where is the time, and properties , , and are the radius, surface tension, and viscosity constant of the droplet; when the viscosity of the outer fluid is taken into account, a toroidal bubble of radius is formed by the outer fluid, and . Employing the Euler equations in the inertial regime, they showed that , where denotes the density of the coalescing droplets. (In the present study, the property will stand for ; designates the set of property constants upon which it depends.)

The paper of [7] is one classic paper which treats the coalescence of droplets both theoretically and numerically. The method used by [7] is termed conventional by [8]. The work of [8] is another classic paper on the subject. They addressed the subject using the conventional method and a different method called interface formation/disappearance method and comparing their corresponding results with experimental findings. The authors argued that their proposed model shows better agreement with experimental data than the conventional one.

In the literature, [4] has given a combined theoretical and experimental account of the coalescence of water droplets on a solid substrate surrounded by an atmosphere of nitrogen saturated with water. For small contact angle and large contact angle , the bidroplet undergoes very fast coalescence. Moreover, the relaxation time characterising the coalescence flow is found to be (or, according to [4], ), and the bulk capillary relaxation time, , where , , and are the radius, surface tension, and viscosity constant of the bidroplet system.

From a hydrodynamical viewpoint, [3] found that the time evolution of the characteristic scale of that small liquid bridge—where the meniscus of the bidroplet system is highly curved—obeys a time-dependent scaling/power law: , where . For instance, in the inertial coalescence regime it was found that , whereas in the viscous coalescence regime . The confirmation of these laws as regards the growth of the property in those coalescence regimes was also reported by earlier scientists. Reference [2] addressed the early-time evolution of the highly curved meniscus bridge experimentally. They found that , where is the inverse of the curvature of the meniscus bridge; by comparison, [7].

On increasing the viscosity of the fluid, [17] has shown the linear dependence of on time. Reference [18] has studied the influence of the geometrical properties of low viscosity liquid droplets experimentally. For contact angles , they deduced , while for their analysis yielded to . Most interestingly, they demonstrated the exponent scaling law in the limit that , thereby unifying the coalescence of liquid droplets and freely suspended droplets in the inertial regime. The exponent scaling law was also argued by [19] who attributed it to the self-similar description of the shape of inviscid pinchoff of a liquid droplet in a still air.

On the other hand, [20] has considered that the influence of surface tension gradient on droplet coalescence was addressed on a combined experimental and numerical basis. The situation considered by these authors was the coalescence of a droplet deposited on the surface of a miscible liquid reservoir. Their results reveal three distinct coalescence regimes based on the reservoir-to-droplet surface tension ratio . For instance, the situation results in the ejection of a small daughter droplet from the top of the top of the coalescing droplet. Contrarily, it was only when that total coalescence could be observed.

Reference [21] has addressed the early-time coalescence of viscous droplets on a flat, wettable substrate. The authors have shown in terms of their respective radii , initial heights , surface tensions , and dynamics viscosities that , which led them to conclude that the evolution of is highly sensible to the geometrical properties of the liquid droplets. Reference [22] considered the coalescence of air bubbles and salt water droplets systems in silicone oils experimentally. The proportionality constant and the -exponent of the scaling law are analysed; the authors argued that the evolution of is independent of the viscosity of its outer fluid, changing only the proportionality constant by a factor roughly equal to and leaving the -exponent unaffected.

The above references reveal extensive studies of the characters of the property with respect to some properly/carefully chosen pair of scalar fields , paying little attention to the characters of the coalescence flow fields in three-dimensional space and on cutting planes and leaving the simultaneous establishment of these power law regimes from one and only one generic equation untouched. Moreover, should an attempt be made to study a class of physical situations ranging from droplets spreading to self-propulsion and coalescence, the several modelling frameworks proposed by those scientists are inconvenient to employ without reconstructing either the mathematical framework or the numerical framework, though the outcomes of their works are no doubt striking. Thus, on the one hand, we must admit that a -evolution equation which permits the deduction of those power laws has never been derived and, on the other hand, the models proposed by those authors make it difficult, if not impossible, to display the global characters of the coalescence flow fields as well as to adapt from one physical situation to another.

In this part, we develop a sound framework which upon variations of an interdependent leading parameter permits to study either the coalescence of two liquid droplets or the spreading of a single droplet on a three-dimensional solid substrate. Taking advantage of the proposed geometrical model, we then establish through theoretical analysis that -evolution equation and derive from it the power law relationships (, , ).

To demonstrate that the proposed model might be adopted with advantage, the , , -exponent scaling laws are proved numerically, and the characters of the coalescence flow fields are illustrated in three-dimensional space and on cutting planes.

2. Theory

2.1. Mathematical Formulations

To formulate the problem mathematically, we reason in the following manner. Suppose two liquid droplets , coalesce on a solid substrate , at some time they will occupy a master domain satisfying the geometrical conditionFollowing [8], equality follows from (1) only when the size of the interfacial circular sector centered at the point (see Figure 1) is zero, and strict inequality whenever it is nonzero, as is easily understood. In the present work, the domain , which can vary in time, will be termed coalescence domain, though the term bidroplet system is also attributed to; the set (with ) denotes the period of time over which the study is carried out. Mathematically, the coalescence domain is bounded by two disjoint sets called the interface and the footprint of the coalescence domain. In Figure 1, we illustrate the coalescence of two viscous liquid droplets , conditioned by zero-flux of momentum, mass and energy across the (free) surface and the (contact) surface ; the null flux vectors (with , , ) stand for those boundary conditions. Figure 2 exposes to view the projection of the coalescence domain in the -plane (top) and in the -plane (bottom). Points , designate the rear and centroid of the bidroplet system; and, points , the centroids of its first and second droplets, both taken at the instant of time . Quantities , (with ) are the initial radii of the coalescing droplets, and the characteristic size at their corresponding points of contact. In particular, it is assumed in Figure 2 that the coalescence domain is initially bihemispherical, units in length in the direction of , with respective radii , (with , by hypothesis), and satisfying the condition . (We agree that necessarily implies that and, implies that the coalescing droplets have not yet been brought into contact; see Figure 2.)

The dynamic state of the coalescence domain is characterised completely by the following macroscopic vector and scalar fields quantities: the velocity and the pressure . We draw attention that, the vector field describes the coalescence flow within the coalescence domain .

The governing equations describing these property variables in are the following generic momentum and mass transport equations of a divergence-free flow field:In the above system of balance laws, the scalar field denotes the pressures field and the vector field denotes the density of volume force due to the surface tension. The tensor denotes the incompressible Newtonian viscous stress tensor, defined asIn (2), the property constants , are, respectively, the density and dynamic viscosity of the mixed fluid occupying the domain . The operator designates the tensor product operator; for given vector field variables , , their tensor product yields to a matrix with the following entries: (with , ).

To model the vector field , we employ the interfacial scalar field into the following alternative expressionThe detailed derivation of (4) is found in [23]. The vector is the unit outer vector normal to and (, since, for , ), respectively. The operator is the surface gradient operator on , and the matrix stands for the identity matrix.

2.2. Boundary and Initial Conditions

At the instant of time the distributions of the velocity field , the pressure field and the chemical concentration should be supplemented in , and thereafter conditioned on boundaries and for all . These are discussed in the following two sections.

The dynamic boundary conditions imposed on the interface for the velocity field are the Neumann-type conditionsIn the above conditions, the property designates the external pressure applied on the interface ; is the force per unit area due to local curvature of that interface. The term on the right-hand side of the second equation designates the tangential stress associated with gradients in the surface tension coefficient. The proper boundary conditions to be considered on the contact surface are now discussed. At some time thereafter, the footprint of the coalescence domain is expected to deform from a bidisc shaped domain to a monodisc shaped domain. Therefore, a wetted wall boundary condition is considered, the solid-fluid boundary conditions. Thus, if the scalar quantity designates the Navier-slip coefficient [24], considering the Navier-slip-with-friction boundary condition yieldsIn (5), (6), the vector designates the tangent vectors to boundaries and . The following initial condition is imposed on the vector field at :The momentum and mass transport equations along with the foregoing set of boundary and initial conditions are solved using the COMSOL Moving Mesh Interface, discussed in what follows.

2.3. Moving Coalescence Domain Explained

At this stage, we have beforehand the mathematical framework and its respective domain of definition. Bearing in mind that the bulk motion of the coalescence domain results from topological changes of its interface from one instant of time to another, to adapt the mathematical descriptions with the movement of the coalescence domain we employ the moving mesh method, which we shall now throw some light upon. For more details, see the works of [23, 25, 26].

Let the interface of the coalescence domain designates a material point-set consisting of material points at the instant of time . Under this condition, let this coalescence domain at some time be described as a spatial point-set consisting of points . On the other hand, the coalescence domain being subjected to smooth deformation from one instant to another, there follows thatwhere . Having these material and spatial vector fields at hand, to relate the topological changes of the interface to its bounded bulk fluid, it suffices to equate the associated velocities on the said interface. Thus, if the vector field describes the movement of the interface, there results inNote that, we write instead of to arrest attention that the interface is undeformed initially.

3. Numerical Implementation

In this section, the construction of the physical domain , described in the -space by Figure 1 and in the -plane by Figure 2, and the numerical implementation of those equations and conditions, namely, (1) to (9), describing its dynamic state is outlined.

The bidroplet system—that is, the coalescence domain—is built as a three-dimensional bihemispherical computational domain sitting on a -plane solid substrate. The solid substrate is defined as a rectangular -plane sized units in length in the direction of (with , ). To realise these geometrical constructions, the radii , , and position vectors (with , , ) must be supplied, and are tabulated in Table 1. They are based on the following line of reasoning.

Initially, the bidroplet system is assumed symmetrical with respect to the -plane. On the solid substrate, the position vector denotes its centroid of the coalescence domain ; the position vectors , denote the centroids of its first and second droplets, respectively. These descriptions are clearly represented in Figure 2. To stress the dependence of the property on those geometrical quantities, we consider the following expressions:where , . Furthermore, to stress that the property is a relatively small quantity initially, we impose the condition ; in regard to the freehand sketches of Figure 2, this means that . We note that for all ; thus the coalescence phenomenon is a relative process, since it is studied through the evolution of the relative property . Consequently, the expansion rate at that region where the curvature is very large is very dependent upon the magnitude of . Clearly, the smaller be the size of the property , the higher be the curvature at the common point of contact of two droplets, and, consequently, the faster be the coalescence rate at early-time.

Topologically, the case for which prescribes a one-point interface, while for the other cases for which , all prescribe infinitely many-points interface. Since those interfaces disappear for all , the first case might well be termed the one-point model and the other three cases infinitely many-points model.

Granted the above descriptions, the following set of equations is implemented: the COMSOL Laminar Two-Phase Flow Moving Mesh equations set, having appropriate interfaces of their initial and boundary conditions. This model, as its name implies, is used to describe the spatiotemporal descriptions of the variable field function . To adapt the COMSOL -equation model with the physical situation under consideration, the COMSOL default model in question is expanded so as to suit the -set of initial and boundary conditions given in earlier (theoretical) sections. The values attributed to those conditions are reported in Table 1.

A series of CFD experiments are ran on this ground and their results are presented in what follows.

4. Results and Discussion

The results to which the proposed model led on a combined theoretical and numerical basis are discussed here.

4.1. Preliminaries

The highly curved meniscus is a fundamental entity. On the grounds of earlier works [3, 18, 19, 21, 27], it can be characterised as follows. By hypothesis, let the -dimensional sequence denotes the property-set of positive property constants, upon which rests the character of the coalescence dynamics. Further, let be a -dependent property, where the subscript stands for the associated -exponent law. Then, following earlier works one can postulate the following expression:Given , the fundamental units of length and time, the dimensional unit of the property is deduced as follows: We observe that , therefore, , where . Moreover, if by properties , , are meant the characteristic scaling of the scalar fields , , , the dimensionless form of (11) reads: , where for some , meaning that the law is independent of the working Newtonian fluids. Consequently, sufficient is that our CFD framework is capable of proving coalescence laws. These statements motivate us to study the coalescence of the bihemispherically shaped droplet system.

5. The Laws of Predicted Theoretically

Granted the proposed three-dimensional coalescence domain (see Figure 1), and the accompanying two-dimensional geometrical model (see Figure 2) wherein are reflected its geometrical properties in the -plane and in the -plane, one can predict the time evolution of the property in the power law regimes. In the present section, the theoretical prediction of the property in those regimes is undertaken. Let it be given , , and , and consider the elementary domain defined bywith boundary , at the instant of time (see Figure 3 for geometrical insight). Note that, the elementary domain may be conceived of as a degree sector of -order thick aligned with the -plane, and having its footprint centered at the point . The region is characterised geometrically by the following length and width at time : , . Following [7], the coalescence flow is driven by the time evolution of . One can therefore postulate thatThe placement of the velocity vector fields , is illustrated in Figure 3. The coalescence dynamics being dominated radially, consequently, the radial component of the equation of motion describing the evolution of the region —from the standpoint of Eulerian description—might at some time be written in the formThe establishment of (14) is based on these chains of reasoning. In earlier sections, the standard basis was used. Here, use of the cylindrical basis will be made to prove (14). Thus, every point on the basis can be represented by on the basis . For , it is postulated that the rate of change of linear momentum equates the differences of capillary and viscous forces across boundary , i.e.,With this view in mind, one then reasons as follows: The curve traced by a -plane cutting the boundary is described by an equation of the formConsequently, the curvature of that curve yieldsOn the other hand, an -volume of and an -surface of gives:respectively. Note in passing that the plane and are clearly illustrated in Figure 4. Bearing in mind that the property is a -dependent scalar field, the rate of change of the linear momentum of along , say, then results inThe capillary force difference sustained across the interface is given byThat due to viscous force is given byIn (21), the tensor , with -dependent entries , where is the row index and the column index, stands for the viscous stress tensor. Substituting (19) through (21) into (15) yieldsFor points , it is asserted that on the basis of postulate (13), yieldingApproximating the last term of (22) by (23) and then setting the triple , (14) follows. It is to be noted here that is a -dependent sequence because of the dependence of on , and on , as is easily understood. One is now well prepared to discuss the so-called -exponent laws, reals.

Granted , after its substitution into (14), one arrives at the following algebraic expression:For that matter of analytical inquiries, setrespectively. Mathematically, it is proposed to predict the power law behaviour of the property based upon the magnitudes of . Since, (24) is a three-term equation, the predictions are argued on three possible pairwise terms.

5.1. Case 1

Suppose that , and . The latter may then be neglected. Thus, for every fixed set of property constants, the following conditions must be fulfilled, namely, . The law describing this situation therefore writes

5.2. Case 2

Let it be assumed that , and , respectively. Then, analogous to the above statements, one obtains . Thus, is in this case of the form

5.3. Case 3

Conceive the orders of magnitudes of with respect to themselves to be as follows: , and . In this case, one deduces . Consequently, the scaling law describing the time evolution of the property writesIt is to me pointed out that by is meant the order of magnitude of the property (with , , ) with respect to time.

For liquid droplets whose viscosities are not high enough, and having contact angles , practical experiments have shown that [18, 19]. Contrarily, for liquids whose viscosities are high enough and having contact angles , many authors have established the law [18, 21, 27]. On the other hand, in purely viscous coalescence regime, the law has also been deduced experimentally [3].

In the next section are proved the power law regimes numerically, thereby giving agreement to experimental data.

6. The Laws of Proved Numerically

Having established the power laws , , theoretically, the chief aim of the present section is to validate them numerically. Before doing so, the groundwork necessary for such endeavours is first laid.

Given the triple -dependent proportionality constants , where . Write to stress that the prefactor referred to is the one associated with the -exponent law and, let quantities , stand for its theoretical and numerical predictions, respectively. Thus, to express an -order discrepancy between the theoretical predictions and the simulation results, one can write and to negate the latter one can write . The question arises: Which of those -exponent laws, , should be computed first? The answer to this question lies in the following paragraph.

Suppose two droplets , having the same size coalesce and time evolves, the bidroplet system then will experience a viscous-inertial coalescence regime which is preceded by a viscous coalescence regime and followed by an inertial coalescence regime, as is easily understood. This shows the existence of two overlapping regimes. Moreover, if the set designates the interval of time over which the law holds for some , , then it is clear that the time intervals corresponding to those three coalescence regimes interrelate themselves as follows:These open sets are termed the coalescence times of the corresponding -regimes, . For some , , it may be said that those two overlapping regimes occur over the time intervalsIn view of the above discussion, it follows, in chronological order, that the establishment of -exponent law must precede that of the -exponent law and follow that of the -exponent law. Indeed, if one is to express for every the law by making use of the open set , the results would resolve into three algebraic expressions of the forms:Carrying out the discussion a step further, one can generalise the form of the scalar field with respect to . Supposed granted , where , and a carefully chosen set of -dependent proportionality constants, with the further supposition that , i.e., . Then, it is no error to state that the time evolution of in the most general case is a series function of the form:Leaving (32) out of consideration, thought it appears to demonstrate an important law, hence, might be a fairer, if not the fairest, approximation in practice, the following sections are devoted to the establishment of the laws , where .

6.1. The Viscous Regime Law

In this section, the numerical proof of the viscous coalescence power law is proved. How the numerical study has been undertaken may be understood through the following lines.

Numerical experiments is first run. The geometrical properties used to construct the computational domain are those given in Table 1. The physical properties of the working fluid occupying the domain as well as the initial and boundary conditions imposed on that domain are those tabulated in Table 1. The said power law is then established in the following manner:so that the pair of scalar fields generates a sequence of points in the -plane. Obviously, the dimensional units are such that and (and , ). For every , [ms], the numerical value of is extracted. This gives the necessary data required to describe the time evolution of the geometrical property . By means of (33), the pairs of scalar fields are collected in the viscous regime. When plotted in the -plane, one obtains the graphical result exposed in Figure 5.

In Figure 5, circular markers are used to represent the numerical data, and the solid line the corresponding fitted power law. This power law function is shaped approximately as follows:From Figure 5, one sees that only for . The critical instant of time at which is observed at [s]; thus, approximately [s].

6.2. The Viscous-Inertial Regime Law

The viscous-inertial coalescence law is demonstrated numerically in the present section. The procedure is similar to that discussed in the above section, differing only in the upper bound of the time . That is, to demonstrate numerically the law , the property values exposed in Table 1 are used, and set [s] (i.e., ). Furthermore, since knowing the pair of scalar fields yields the desired law in the -plane numerically, to settle the -exponent law one then proceed in likely the same way as in the earlier section (see Section 6.1).

In Figure 6 are graphed the graphical results. Circular markers are used to represent the numerical data, and the solid line the equation corresponding fitted power law. That fitted power law curve is shaped asIt is understood that only for . The critical instant of time at which occurs at [s]; thus, approximately [s].

6.3. The Inertial Regime Law

The power law is here validated numerically. To do so numerical experiments were first run. The geometrical properties, initial and boundary conditions imposed, and the physical properties of the working fluid are those tabulated in Table 1. The said power law is then established in exactly the same way as in the foregoing two sections (see Sections 6.1 and 6.2).

In Figure 7 are exposed to view their graphical results. Circular markers are used to represent the numerical data, and the solid line the corresponding fitted power law. Through numerical experiments, it is believed that the following scalar function—representative of that solid line curve—gives good agreement:In regard to (36), it follows that . The power-the law is thus proved numerically, demonstrated experimentally by [18, 21, 27]. It may be remarked that only for . The critical instant of time at which occurs at [s]; thus, approximately [s]. The values of the scalar fields (, , ) are observed to be roughly equal to these values:Moreover, from (36), (34), (35), one also observes thatIn what follows the -dependent sequence are analysed.

6.4. Of the Characters of

When the literature was surveyed, one found no more than two elements of the sequence have gained attention, yet their establishment in terms of the property constants involved in the problems they ultimately lie have been correctly settled. With this view in mind, it is proposed to discuss in the present section the characters of through analysis.

The analysis commences by relating the leading property associated with the interfacial coalescence hydrodynamics to that associated with the bulk coalescence hydrodynamics, that is, it is proposed to relate to on . It is clear that , where , . Therefore, taking its derivative, one obtainsOn the other hand, the leading property constants being , , , , , consequently, . Given the characteristic scales , one rescales the properties appearing in (39) as follows: , , , and , yieldingIt now remains to express the right-hand side of (40) in terms of the property vector . To do this, one reasons according to the -regimes, where . Analysing the set of properties in terms of , the following orders of magnitude follow:These being so, substituting (41) into (40), one findsMoreover, inserting the property constants of Table 1 into (42) and, thereafter, dividing each element of by the corresponding element of , the orders of magnitude of the ratio of to yieldIn (43), it is understood that quantities , stand for the theoretical and numerical estimation, respectively, of the scalar field . Expressing (34), (35) and (36) in terms of (42) givesEquation (44) at once justifies the reasonableness of the foregoing chains of reasoning and, this completes the discussion of the present section concerning the -dependent evolution of the field . In the next section, it is proposed to enter upon the discussion of the characters of the curve along a carefully chosen one-dimensional boundary.

6.5. Of the Characters of

If the velocity norm of the interacting droplets systems is plotted on the curved boundary , the curve traced by the intersection of the interface with the plane , one obtains Figure 8. The placement of arrows on the curves of Figure 8 serve merely to indicate the evolution of time evolution of these curves. In Figure 8, the graphs of the velocity norm correspond to the sequence of times [ms]. The nature of these curves are now explained.

Inspecting Figure 8, one sees that the highly curved meniscus induces an -order singularity into the coalescence flow in the early-time regime, peaking to magnitudes of the order of mm·s inside that highly curved meniscus zone (indicated by dashed, vertical lines). It is perhaps because of this striking effect that some authors have characterised such motion as singular motion [7, 22].

It is interesting to see the spatiotemporal distributions of the field variable . This is the purpose of the next section.

7. The Coalescence Flow Fields Illustrated

Having discussed the growth of the time-dependent function in the inertial, inertial-viscous and viscous regimes, the evolutions of the vector field function with respect to are now reported. In few words, the present endeavours is an attempt to lay additional stress on the qualitative aspects of the proposed model.

When the coalescence domain undergoes coalescence phenomenon, as soon as its pair of droplets are brought into contact, a mass movement is started, thereby causing the corresponding droplets to fuse. It is not necessary to go into details; sufficient is to demonstrate that the proposed model is capable of describing the underlying motion phenomena and the latter are exposed in the space and on the planes (with , ).

In Figures 9 and 10, the time evolutions of the velocity norm of the coalescence domain at times , , , [ms] are exposed in the space. The corresponding two-dimensional plots of the scalar field are reported, respectively, on the planes (with , ) given in Figures 1114, where solid lines stand for the iso- contours.

8. Closing Remarks

In this part, a three-dimensional modelling methodology has been proposed to study droplets coalescence phenomenon. On the whole, the statements arrived at conclusively are:

When, upon variations of geometrical parameters, the footprint and the free-surface of the bidroplet system reduce to the footprint and the free-surface , respectively, of a monodroplet system , the corresponding boundary conditions are by default adapted into the COMSOL interfaces. As a result, with one numerical model, its property constants can be so arranged so as to compute many physical situations, a task which one would have undertaken had the proposed model been unversatile.

Taking advantage of the geometrical features of the proposed model, it is possible to deduce through theoretical analysis the scaling laws (with , , ), laws already reported in the (experimental) literature.

Carried out numerical proofs of these laws, is observed in the visco-dominated coalescence flow, in the inertia-visco-dominated coalescence flow, and in the inertia-dominated coalescence flow.

Through scaling analysis, it is found that the dependence of the sequence of prefactors on can be approximated as follows: , and .

The flow field variables illustrated by , plots clearly demonstrate that the expected features associated with droplets coalescence phenomenon are well captured by the model.

On the whole, yet the approach in the present analysis has been to present the leading features of the proposed model by considering the coalescence that results under an idealized situation, it is no way contradictory to declare that it can be used to model more complex physical situations. For example, numerical experiments on droplets coalescence phenomenon can be carried out under appropriate CFD conditions to display useful information concerning mixing process driven by Marangoni flow.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.


The current address of M. I. Khodabocus and M. Sellier is Mechanical Engineering Department, University of Canterbury, Private Bag 4800, Christchurch 8140, New Zealand. The current address of V. Nock is Electrical and Computer Engineering Department, University of Canterbury, Private Bag 4800, Christchurch 8140, New Zealand.

Conflicts of Interest

The authors declare that they have no conflicts of interest.