# Separation Criteria for Off-Axis Binary Drop Collisions

Authors: Mary D. Saroka; Nasser Ashgriz
Journal: Journal of Fluids
 (2015)

---

## Abstract

Off-axis collisions of two equal size droplets are investigated numerically.  Various governing processes in such collisions are discussed. Several commonly used theoretical models that predict the onset of separation after collision are evaluated based on the processes observed numerically. A separation criterion based on droplet deformation is found. The numerical results are used to assess the validity of some commonly used phenomenological models for drop separation after collision. Also, a critical Weber number for the droplet separation after grazing collision is reported. The effect of Reynolds number is investigated and regions of permanent coalescence and separation are plotted in a Weber-Reynolds number plane for high impact parameter collisions.


---

## Body

## 1. Introduction

Many engineering applications and natural process involve drop collisions in which the final outcome has a direct impact on the eventual success or failure of the application or process. In sprays, for example, the behavior is characterized by many small liquid drops with a particular size distribution and with each individual drop moving along a particular trajectory with a prescribed velocity. By singling out and focusing on just two drops within the spray, one can readily see that the interaction between these drops plays a part in the overall evolutional characteristics of the spray. One can further imagine that the likelihood of two drops colliding along the same axis of trajectory is less likely to occur than the situation in which the two drops collide along different trajectories. It is the latter case, typically referred to as off-axis collisions, that is the subject of this investigation.Over the years, many researchers have studied the collision dynamics of the binary drops. The early studies on drop collision were motivated by understanding of the physical processes that occur during rain fall [1–12]. Rain drops may collide with each other and break into smaller drops or coalesce and generate larger drops. For instance, Park [5] did experiments on water drops in humid environment to determine conditions for the drop separation after collision. Several phenomenological models for the drop separation after collision are provided, for instance, the rotational energy model [7] and the variation principal model for the minimum potential energy for the stretching separation [12]. The later studies were motivated by drop collision in various spray systems [13–28]. Collisions of many other liquids, such as hydrocarbons [13–18], heavy oils [19, 20], and mercury [21], in addition to water, have also been studied. The bulk of quantitative information on the drop collision is obtained through experiments [21–28], while numerical studies are used to provide detailed understanding of this complex process. The front tracking techniques [29–32] and a combination of level set and VOF methods are used to solve the dynamics of the free surfaces [33–40]. A detailed review of the droplet collision process is provided by Brenn [41]. Review of this body of research indicates that, depending upon the initial energy of the two drops and the angle at which they collide, different types of outcomes will occur. These outcomes can be categorized as bouncing, partial coalescence, permanent coalescence, separation, or shattering. There are also several subcategories for each outcome. Of these five, permanent coalescence and separation are the most frequently observed outcomes, while the other three, namely, bouncing, partial coalescence, and shattering, represent special cases in which the drops have either very small velocities, large size differences, or very high velocities.The present study is on the off-axis collision of two equal size droplets. This study is a continuation of a previous study by the authors on the head-on collision of two drops [40]. Droplet collisions in dynamically inert surroundings are considered. Therefore, bouncing collisions, which are due to air entrapment between the two approaching droplets, are not observed. In the next sections, first the problem statement and the numerical methods are described, followed by the results and analysis.


## 2. Problem Statement

A general schematic depicting the positions and velocity vectors of two drops just prior to contact is shown in Figure1. Suppose that one of the drops has a velocity defined as V1, while the other has a velocity of V2. These individual velocities are used to define a relative velocity, Vr, of the interacting drop. This is done by combining the drop velocities and their impact angle, α, through the following relation: Vr=V12+V22-2V1V2cos⁡α. Another important parameter necessary to define a particular collision involves the distance between the drops’ centers. This distance is referred to as the impact parameter, noted as x. As the impact parameter increases the amount of “interaction” is reduced, while the converse is true. Using this information, it is convenient to define various dimensionless numbers to represent the phenomena under consideration. Three dimensionless parameters, namely, a Weber number, a nondimensional impact parameter, and a drop diameter ratio, are defined based on the drop diameter, liquid density, ρ, and the coefficient of surface tension, σ: (1)We=ρDVr2σ,X=xsin⁡β-γD,where(2)sin⁡γ=V1Vrsin⁡α.In addition, Reynolds and Ohnesorge numbers are used to correlate effect of viscosity, μ, and its relationship with surface tension:(3)Re=ρDVrμ,Oh=WeRe=μρDσ.The Ohnesorge number is unique in the fact that it is only a function of fluid properties and the drop size, while both the Weber and Reynolds numbers are functions of the fluid properties, drop geometry, and impact velocity. In the present study, We, Re, and Oh numbers are calculated based on the liquid properties.Figure 1
General orientation and parameter definition for binary droplet collisions.In order to discuss different types of collision outcomes, some collision images obtained by Ashgriz and Poo [22], but not published earlier, are provided here. The experimental setup and conditions are all as described in Ashgriz and Poo [22] and will not be repeated here. Figure 2(a) shows a close to head-on collision of two water droplets. After the collision, drops spread radially on each other until the kinetic energy of the drops is transformed into the surface energy. The surface energy is large enough to push back the fluid and form a ligament. The ligament stretches until it breaks into two droplets. This is referred to as a “reflexive separation.” Ashgriz and Poo [22] have provided a theoretical model to predict the onset of this type of separation. If the collision impact parameter is increased enough, the impact process becomes complicated. Part of the drop fluid goes through the spreading and reflexive action, and another part of the drop fluid tends to stretch and pull the drops away from each other. If the collision kinetic energy is just enough to have reflexive separation, as shown in Figure 2(a), then a small increase in the impact parameter reduces this reflexive energy, and drops may not separate and stay coalesced. This process is shown in Figure 2(b), where the coalesced drop goes through significant deformation and oscillations, but drops remain coalesced. If the impact parameter is further increased, then the part of each drop which remains uninteracted with the other drop tends to stretch and pull the drops apart. This is observed in Figures 2(c), 2(d), and 2(e), and it is referred to as stretching separation. The present study only addresses stretching separation.Figure 2
Binary collision of water droplets. (a) Reflexive separation due to near head-on collision; (b) coalescence collision; (c) separation collision resulting in two drops; (d) separation collision resulting in three drops; (e) separation collision resulting in four drops; and (f) separation collision resulting in six drops.
(a)
(b)
(c)
(d)
(e)
(f)


## 3. Numerical Methodology

We consider two drops composed of the same fluid with constant properties. The flow inside the drops is solved using the standard conservation of mass and momentum equations:(4)∇·V=0,ρ∂V∂t+ρV·∇V=-∇p+μ∇2V+Fb,where V is the velocity vector, p is the pressure, and Fb is any volume force acting on the drops. At the free surface, the boundary condition requires that both mass and momentum are conserved. For drops surrounded by a less dense/viscous fluid or gas, the viscous stresses at the surface are zero. Additionally, assume a constant surface tension coefficient allows the stress boundary condition to reduce to Laplace’s equation ps=σκ, where ps is the pressure jump across the surface and κ is the surface curvature. Details of the numerical simulation have been discussed in earlier publications [40, 42, 43] and are not repeated here.Two drops with a prescribed distance separating their respective centers are considered. In all cases, this distance was defined to occur in one direction only. That is to say that for drops translating in one direction (i.e.,z-direction) their centers were placed in the same plane (called the defining plane) composed of this direction and one of the other two remaining directions (i.e., y-z or x-z plane). The third direction (i.e., x or y, resp.) was held constant for both drops. One might imagine that this third direction could also represent an off-set, thus indicating a second impact parameter. However, this case was not considered during this investigation. Figure 3 provides an example of the defining plane for an off-axis collision. This plane represents the y-z plane and is useful in observing the internal flow patterns as the two drops collide and evolve. Although similar planes such as X-Z or X-Y planes could be defined, they would only represent a specific region of the collision and are therefore not as comprehensive.Figure 3
Definition of they-z plane. Original two drops translating in the ± z-direction.The distance between the drops’ centers was used to compute the impact parameter. For this investigation, impact parameters ranging from 0 to 1 were considered, with the bulk of the simulations focusing on collisions forX>0.5 (i.e., high impact parameters). In order to cause the drops to collide, each drop was given a nonzero initial velocity. These individual drop velocities were subsequently used to compute the relative velocity, Vr. All simulations were conducted using water drops for Weber numbers ranging from 10 to 60 and Reynolds numbers ranging from 135 to 1639. Since various drop sizes ranging from 10 μm to 2 mm were used in the course of this investigation, a wide range of relative velocities ranging from 0.5 m/s to 21 m/s were considered. For very small drops of 10 μm, relative velocities as high as 21 m/s were used, while for the larger drops of 2 mm, velocities less than 1 m/s were predominately used. The fluid properties, density, ρ, viscosity, μ, and the coefficient of the surface tension, σ, were kept constant at 1000 kg/m3, 10−6 m2/s, and 0.073 N/m, respectively.In order to resolve the flow within the drops, 10 cells per radius were used. A computational grid of 70 × 70 × 180 cells in thex-, y-, and z-directions was used to allow the collision to evolve without hitting any walls. The use of 180 cells in the z-direction was necessary, since this was the primary direction in which the drops approached each other. Higher Weber numbers than those considered were not simulated, since these types of collisions required a larger computational grid to fully capture the evolutional dynamics of the collision.


## 4. Results and Discussion

Figures4–6 show typical surface evolutions of water drops colliding at impact parameters greater than zero. Figure 4 shows a coalescence collision of two 300 μm drops colliding at We = 25, Re = 740, and X = 0.7. Figure 5 shows results for a similar case as in Figure 4 but with a slightly larger relative velocity, where two 300 μm drops collide at We = 30, Re = 811, and X = 0.7. It is shown that the boundary between the coalescence and separation at the impact parameter of X = 0.7 occurs at a Weber number between 25 < We < 30. A small increase in the impact We number causes the drops to separate after collision. Comparing Figures 4 and 5 with each other, it is clear that, in Figure 4, the surface energy pulls the ligament formed at stage 5 back together, and drops coalesce. In Figure 5, the ligament formed at stage 10 is longer than that in Figure 4. In this case, the pinching effect, which is governed by the Rayleigh type effect, is faster than the forces that tend to coalesce the drops. Therefore, the ligament in Figure 5 breaks before it comes back forming a coalesced droplet.Figure 4
Coalescence collision of two 300μm water droplets with We = 25, Re = 740, and X = 0.7. Numbers in the figure indicate time = (0.1 ms)i, where i = image number.Figure 5
Separation collision of two 300μm water droplets with We = 30, Re = 810, and X = 0.7. The numbers in figure indicate time = (0.05 ms)i, where i = image number.Figure 6
Stretching separation with one satellite droplet for the collision of two 700μm water droplets with We = 60, Re = 1751, and X = 0.7. The numbers in the figure indicate time = (0.2 ms)i, where i = image number.A general description of an off-axis collision process is as follows. As the drops come into contact with one another, a small region of interaction develops, while the remaining portion of the drops continues on their original course. This stretches the collided region. As the combined mass continues to stretch, the interaction region allows a portion of the fluid at the ends to drain back towards the middle of the system. This in turn causes the entire mass to rotate in a clockwise direction when the top drop is to the right of the bottom drop (opposite orientation results in counterclockwise rotation). Fluid within the mass continues to adjust itself in this manner until a liquid bridge between the two ends is formed. Depending upon the amount of initial energy present in the system and the magnitude of the impact parameter, the rotating mass will either contract resulting in permanent coalescence or separate without or with satellite drops. Figure6 shows collision of two 700 μm drops colliding at We = 60, Re = 1751, and X = 0.7. In this case, a long ligament is formed, which later results in a satellite droplet in addition to the two main droplets. It can be seen in images 16 to 24 that the rate of contraction of the inner ligament is faster than the end pinching of the drops. Therefore, only one satellite is formed. If the inner ligament was slightly longer, end pinching could result in three satellite droplets.


### 4.1. Internal Flow Patterns

There has been a lot of debate about the existence of internal circulations in droplets of a spray. Since it is very difficult to see the internal flow of a moving droplet, the main way to determine this effect has been through numerical modeling. The present study puts forward a concept for the formation of the internal circulations in the droplets. Figure7 shows the internal flow patterns for the drop collision of Figure 6 (We = 60, Re = 1751, and X = 0.7) and for the plane described in Figure 3. The first image in Figure 7 shows the two drops just after contact. Immediately after the collision a small region of interaction is produced. In this region the fluid is redirected away from the collision point. As this occurs, the outer (or noninteracting) portion continues to translate unimpeded in the original direction (as shown in images 2–4). In the fifth image, in the center of the combined mass some of the flow begins to be redirected “radially” while the bulk of the motion still occurs in original translating direction. However, surface tension compensates and the radial flow is redirected, thus causing the overall flow to be directed upward/downward and away from the center. This action causes the stretching phenomena observed in off-axis collisions. The bulk of the fluid eventually ends up in the bulbous ends as indicated by the streamlines in images 7–9. Image 10 indicates that separation has occurred with the production of two parent drops and a ligament. It is interesting to note several observations at this point. First, review of early images does not provide any evidence of the pinching mechanism in this defining plane. Therefore, the pinching phenomena are the result of radial type flow. Second, the parent drops exhibit flow patterns representing translating flow. Unlike the scenario for the head-on collisions [40] these parent drops reflect similar flow patterns of the original two drops which causes them to continue to translate after separation occurs. Lastly, the ligament collapses into a satellite drop producing the symmetrical flow pattern observed. Therefore, main droplets after collision have unidirectional internal flows, whereas satellite droplets may have multidirectional and internal flow circulations.Figure 7
Internal flow patterns for they-z plane. Results are for a binary water drop collision with a diameter of 700 μm and We = 60, Re = 1751.


### 4.2. The Upper Boundary between Coalescence and Separation

The upper boundary between a coalescence and separation collision occurs for grazing collisions, that is,X = 1. As soon as the surfaces of two drops touch, the liquid flows between them. Depending on the collision conditions, the two drops may touch and temporarily coalesce, and then separate, or remain coalesced. The We number corresponding to this condition is referred to as the upper critical We number, Wecr. Table 1 shows the numerically determined upper critical Weber number for equally sized water drops as a function of drop diameter. It is clear that the upper critical Weber number increases with decreasing diameter. The pressure inside smaller drops is larger than those in larger drops; therefore, the flow is pushed into the liquid bridge between the two drops faster in the smaller drops than for the larger drops. The larger the mass of the bridge, the more difficult it is to break the drops. Therefore, the separation occurs at higher Weber numbers for the smaller drop collisions. The curve fit to the data of Table 1 is Wecr=-2.7lnD+36. As the drop size increases, the effect of change in drop size becomes less important. The upper critical We number approaches We = 18 for large drop sizes.Table 1
Upper critical Weber numbers for two equal size water droplets (X=1).
Diameterμm
We
cr
Re
10
32.5
154
75
22
347
150
21
480
300
20.3
665
700
18.3
966
2000
18
1622


### 4.3. Stretching Separation Boundary

The boundary between permanent coalescence and separation for off-axis collision of two equal size droplets for different collision impact parameters is obtained. Figure8 presents the results of this investigation for equal size drop collisions with a diameter of 300 μm (Figure 8(a)) and 700 μm (Figure 8(b)). Outcomes of permanent coalescence, separation without satellite production, and separation with one satellite produced are plotted for a range of Weber numbers and impact parameters. Also included in the figures are the boundaries defined by the analytical models proposed by Ashgriz and Poo [22], Arkhipov et al. [12], Brazier-Smith et al. [7], and Park [5]. The model by Ashgriz and Poo [22] is referred to as AP/EXP model for the remainder of the discussion. Similarly, the other models are referred to as A, B-S, and P, respectively. The results show that, for a constant impact parameter, as the Weber number is increased, the outcome changes from coalescence to separation without satellites to separation with a satellite. Similarly, as the Weber number is held constant, the same progression is observed provided there is enough energy initially present in the drops (i.e., at low Weber numbers, regardless of impact parameter, coalescence will always occur).Figure 8
Coalescence-fragmentation curves: (a) 300μm drop collisions and (b) 700 μm drop collisions. Symbols represent numerical data while lines represent boundary models presented by previous researchers.
(a)
(b)The numerical division between the coalescence and separation occurs at higher Weber and impact parameter combinations than those observed by the AP/EXP and B-S models. A possible explanation for this behavior is attributed to the existence of trapped air in the gap between the colliding drops. In the experiments, the impact parameter is determined by measuring the lines of trajectory and the corresponding separating angle [22]. Depending upon when this is done (i.e., a few diameters from the drop generators or just prior to impact) the value of the impact parameter can be quite different. As the two drops travel through air, their trajectories remain relatively constant. However, as they come closer and closer, the surrounding gas becomes trapped in between them, which results in the slippage of the drops on each other. This may alter the original impact parameter based on the drop trajectories prior to collision. The effect of such a slippage is shown in Figure 9 by a shift in data. In this figure, the solid line represents the boundary using the AP/EXP model while the symbols represent the numerical results for the stretching separation without satellite production. The dashed line represents a linear translation of the AP/EXP model and is meant to convey the concept of impact parameter slippage. For example, at We = 31, the impact parameter on the AP/EXP boundary is X = 0.44, while that numerically is Xs = 0.64. Thus, due to the presence of the trapped air, the impact parameter at the time of collision could actually be higher than that previously thought.Figure 9
Impact parameter deviation for the same Weber numbers.X is a point on the experimental boundary, while Xs represents the impact parameter observed numerically. The solid line represents the boundary model by Ashgriz and Poo [22], while the dashed line represents a linear translation added to their model.The presence of the trapped air has three different effects on the drops, which ultimately cause an increase (or slippage) in the effective impact parameter. These effects are (i) altered trajectories, (ii) altered internal flow fields, and (iii) shape deformation. The trapped air conveys an external force on the drops, which may alter their trajectories, as shown in Figure10. The second effect is shown in Figure 11, where drops moving in air may develop internal flows. The internal flow changes the drop dynamics after collision and may change the collision outcome. Finally, the surrounding air may cause drop deformation, as shown in Figure 12. As a result of this deformation the region of contact may increase as opposed to the case in which the drops stay perfectly spherical. This increased contact region may result in a larger interaction region and thus more mixing. It is, therefore, concluded that, in determining the collision impact parameter, one has to be worried about the slippage effect that may occur due to the air entrapment.Figure 10
Effect of trapped air between colliding drops. Experimental photo from Ashgriz (unpublished material from Ashgriz and Poo [22]). The trapped air alters the collision impact parameter based on droplet trajectories.Figure 11
Internal flow of two drops approaching each other in (a) vacuum and (b) in air.
(a)
(b)Figure 12
Two cases show the effect of trapped air in deforming the drop prior to collision and merging [unpublished photos from Ashgriz and Poo experiments [22]].
(a)
(b)Figure8 also compares the numerical results with models by Arkhipov et al. (A) [12] and Park (P) [5]. Arkhipov assumed that after the collision a spherical drop is formed that rotates with a constant angular velocity, while Park assumed that the combined mass rotates as two connected spheres. The numerical results, as evident in Figure 5, show that, for the drops to separate after collision, they rotate as two drops that are connected, and not as a combined spherical mass. Therefore, Park’s model is closer to numerical observation of the collision dynamics. However, the drop shapes and masses are altered and are not exactly the same as the original droplets. This becomes particularly true as the Weber number is increased. However, at low Weber numbers, We < 25, the numerical results do not correlate with any of the models considered. This is most likely explained by an incorrect assessment of the amount of angular velocity produced during the evolution. One could imagine that the amount of angular velocity would in fact not be constant (as assumed by Arkhipov et al.) and be a function of the impact parameter instead.


### 4.4.  Deformation

The deformation variations in each direction for an off-axis collision of two water drops with a Weber number of 30 and a Reynolds number of 810 are provided in Figures13–16. For clarity, the deformations are grouped according to either low impact parameters, X<0.5, or high impact parameters, X>0.5. It is convenient to nondimensionalize the amount of deformation by dividing by the drop’s initial diameter. In doing so, the following definitions are used to represent the dimensionless deformation in each direction: (5)η=LxD,γ=LyD,ψ=LzD,where, Lx is the amount of deformation in the x-direction, Ly is the amount of deformation in the y-direction, and Lz is the amount of deformation in the z-direction. Since off-axis collisions occur on different trajectories, the evolution process occurs in all three dimensions, resulting in separate behavior for each direction. In the next three subsections, the behavior in each direction is analyzed.Figure 13
x-direction deformation for We = 30 and Re = 810. (a) Low impact parameters, X < 0.5, and (b) high impact parameters, X ≥ 0.5.
(a)
(b)Figure 14
y-direction deformation for We = 30 and Re = 810. (a) Low impact parameters, X < 0.5, and (b) high impact parameters, X ≥ 0.5.
(a)
(b)Figure 15
z-direction deformation for We = 30 and Re = 810. (a) Low impact parameters, X < 0.5, and (b) high impact parameters, X ≥ 0.5.
(a)
(b)Figure 16
Maximumz-deformation for various Weber numbers and impact parameters. (a) 300 μm drops and (b) 700 μm drops. Solid symbols, coalescence, and open symbols, separation. Boundary line, ψmax = 0.006 We + 3.
(a)
(b)As observed with head-on collisions [40], three major periods of the collision process are spreading and production of a disk, period I, contraction of the disk, period II, and stretching along the collision axis (z-direction), period III.


### 4.5.x-Direction (Minor Axis 1)

Figure13 contains the dimensionless deformation in the x-direction. Review of this figure during period I provides the following information. First, it is observed that, at t∗=0, η increases as the impact parameter, X, increases. This behavior is a direct consequence of the impact parameter definition, since the drop’s separating distance is defined in this direction. Second, as the impact parameter is increased the amount of spread decreases. Therefore, as the region of interaction becomes smaller due to the increased impact parameter, less fluid is allowed to flow in this direction. Additionally, at X~0.6, no spread is observed and the deformation actually becomes drainage or contraction for higher impact parameters. This would seem to imply that, at this (X) and higher impact parameters, fluid is actually pulled out of this direction. The amount of time this period takes to occur is relatively constant regardless of impact parameter.Period II was defined to represent the action of fluid contraction of the disk formed during period I. When the variation of impact parameter is considered, the following observations can be made. First, the amount of contraction is relatively constant until a value ofη~1.3 is obtained at the end of this period for X<0.3. Essentially, this means that, for low impact parameters, the redirected fluid momentum in this direction dominates the flow which in turn causes the disk to contract inward. Second, as the impact parameter is increased until X~0.6, more contraction occurs until a minimum value of η~0.75 is reached. In this range of 0.3<X<0.6, a competition between the amount of fluid momentum and the surface tension effect becomes evident, with the fluid momentum effect decreasing as the impact parameter is increased (i.e., due to the reduction in the interaction region). For impact parameters greater than 0.7, this minimum deformation value is maintained. Therefore, at these high impact parameters, the amount of deformation does not change since the surface tension is now the controlling factor.


### 4.6.y-Direction (Minor Axis 2)

The dimensionless deformation in they-direction is provided in Figure 14. As with the x-direction variation, period I behavior still exists regardless of impact parameter variation. However, the amount of spread decreases as the impact parameter increases. Again, this is seen as a consequence of the reduction of the interaction region. Unlike in the x-direction, spread is observed up until X~0.8 (instead of X~0.6 as seen in the x-direction). At higher impact parameters, the amount of deformation shifts to contraction instead. This indicates that more mixing occurs in the y-direction than in the x-direction, since the fluid is still able to spread outward for higher impact parameters.During period II the overall amount of contraction remains constant (since similar slopes are observed) for low impact parameters ofX<0.3 (i.e., γ at the end of this period decreases from 0.85 to 0.7). Therefore, even though more deformation is observed at the end of this period, the control parameter is still fluid momentum as it is for the x-direction case. The smaller γ values are the consequence of smaller γmax values observed during period I. For impact parameters in the range of 0.4–0.6, the amount of deformation at the end of this period is constant and has a value of γ~0.6 indicating a shift in dominance from fluid momentum to surface tension. As the impact parameter is further increased, X>0.7, less overall deformation occurs since the amount of initial spread seen during period I has decreased and surface tension effect is in control.


### 4.7.z-Direction (Major Axis)

In general, as the impact parameter increases, the region of interaction between the two drops decreases causing the intensity of deformation to shift from thex- and y-directions to the z-direction. This is clearly observed in the z-deformation curves of Figure 15. For an impact parameter of 0, the maximum deformation obtained is 2.7 times the drop’s initial diameter and occurs at tc∗=6.7. As the impact parameter is increased the maximum height obtained decreases and occurs earlier in time until an impact parameter of 0.3 is reached. At this impact parameter the maximum deformation is 1.8 times the diameter and occurs at tc∗=4.1. As the impact parameter is increased further the maximum deformation increases linearly until ψmax=3 times the diameter is obtained. As with head-on collisions, this value defines a limiting value in which either permanent coalescence or separation occurs. This is observed at impact parameters of 0.8 and 0.9, in which the combined mass separates by the pinching mechanism explained earlier. In the case of the impact parameter of 0.7, the contraction rate is greater than the pinching mechanism and stretching-coalescence is observed.The behavior of this maximum deformation parameter,ψmax, as a function of impact parameter and Weber number is plotted in Figure 16. Figure 16(a) represents the results for collisions with drop diameters of 300 μm and Figure 16(b) is for collision with diameters of 700 μm. Open symbols indicate the collision outcome is coalescence while solid symbols indicate separation collision. Review of this figure indicates several aspects. First, the behavior of ψmax is linear over the range of Weber numbers simulated and for illustrative purposes a straight line is fitted through each set of impact data. The slope of each impact parameter data set, m, is computed and presented in Table 2. The subscripts 300 and 700 refer to the drop size. Also included in the table is the change in slope, δm, which indicates the difference in the slope as the drop size is decreased from 700 μm to 300 μm.Table 2
Slopes of lines in Figure16.
X
m
D
=
300
μ
m
m
D
=
700
μ
m
δ
m
0.4
0.018
n/a
—
0.5
0.038
0.046
0.008
0.6
0.043
0.051
0.008
0.7
0.058
0.063
0.005
0.8
0.073
0.078
0.005
0.9
0.088
0.087
0.001
0.95
n/a
0.096
—Regardless of drop size, the slope of the lines increases with the increase in the impact parameter. This is attributed to the reduction in the interaction region which in turn causes the bulk of the fluid motion to continue traveling in the originalz-direction. It is also interesting to note that the drop size has a small impact on the final stretch length. This is evident by the increased slope values for the 700 μm drop collisions as opposed to 300 μm drop collision occurring at the same impact parameter. Additionally, the slopes for X=0.9 are very similar and as a result this indicates that the size effect is no longer a factor in the overall stretching process.A last comment concerning Figure16 concerns the boundary line (represented by a thick dashed line) indicating the regions of coalescence and separation. According to Lord Rayleigh’s capillary instability model [44], this boundary has a value of ψmax=π. However, for off-axis collisions, this boundary was empirically determined to be(6)ψmax=0.006We+3,where the slope reflects the impact parameter effect.


### 4.8. Dissipation Based Models

At the beginning of the collision, both drops have a known amount of surface,SEo, and kinetic energies, KEo. These energies are computed as follows:(7)SEo=2πDi2σ,KEo=π24ρDi3Vr2,where the ratio of the kinetic energy to the surface energy is defined as We∗ = We/48 [22]. Each of these energies is computed at each time step using the following definitions:(8)SEn=σ∑p=1NcellsSAp,KEn=12∑ijk=1Ncellsfijkρvcelluijk2+vijk2+wijk2,where SAp is the polygon area in the computational domain and n is the time index, Ncells is number of computational cells, fijk is volume fraction of the liquid in the cell, and uijk,vijk, and wijk are the velocity components in each direction. Therefore, using these relationships each energy variation is computed at each time step during the collisional process. Additionally, by using these energies and conservation of energy an estimation of the change in total energy or dissipation energy, ϕ, is obtained as follows:(9)ϕ=KEo+SEo-KEn-SEn.The kinetic and surface energies are scaled with the initial surface energy as α=KEn/SEo and β=SEn/SEo. The initial surface energy is chosen as the scaling energy since for a same-size drop collision, it is independent of the impact energy. Since the amount of dissipated energy is simply the combined transient kinetic and surface energy subtracted from the initial total energy, a nondimensional dissipation energy is defined as(10)Φ=1-KEn+SEnKEo+SEo.


#### 4.8.1. Kinetic Energy

The dimensionless kinetic energy variation for the binary collision of two water drops for We = 30 and Re = 810 is provided in Figure17. As the impact parameter is increased, the maximum kinetic energy, αmax, decreases from 0.25 for X=0 to 0.22 for 0<X<0.6. It should be noted that, for this range of impact parameters, the final outcome was coalescence, while higher impact parameters (X>0.7) resulted in separation. This maximum kinetic energy decrease is the result of more energy conversion to surface energy (due to the increased surface area accompanying an increase in impact parameter). Additionally, whereas head-on collisions exhibit defined peak and valley behavior corresponding to times of maximum spread and contraction, an increase in the impact parameter causes the change in kinetic energy to undergo a more chaotic variation. It is not until the impact parameter is increased to a value greater than 0.6 that no oscillatory behavior (of any kind) is observed. Therefore, if the kinetic energy does not peak early in the evolutional sequence, separation will likely occur.Figure 17
Dimensionless kinetic energy,α = KEn/SEo for We = 30 and Re = 810. (a) Low impact parameters, X < 0.5, and (b) high impact parameters, X ≥ 0.5.
(a)
(b)


#### 4.8.2. Surface Energy

Figure18 provides the dimensionless surface energy variations for the binary collision case under consideration. The maximum surface energy, βmax, decreases from a value of 1.325 for X=0 to 1.025 for the range, 0<X<0.6. This behavior parallels the behavior seen for the kinetic energy. At impact parameters greater than 0.6, the absence of this peak behavior at this time in the evolution is conspicuous. In fact at this time, βmax~1 and as noted earlier in the previous discussion, all the remaining collisions at higher impact parameters resulted in separation. Additionally, the variation in surface energy becomes more oscillatory or chaotic for impact parameters up to 0.6. At these low and intermediate impact parameters, various surface deformations occur as the combined mass tries to stabilize itself to its minimum energy state.Figure 18
Dimensionless surface energy,β = SEn/SEo for We = 30 and Re = 810. (a) Low impact parameters, X < 0.5, and (b) high impact parameters, X ≥ 0.5.
(a)
(b)


### 4.9. Reynolds Number Effect

The Reynolds number effect is presented based on the collision of two drops with impact parameters ofX=1 and X=0.9. For the grazing collision (X=1), 25 different water drop collisions with drop sizes of 10 μm, 75 μm, 150 μm, 300 μm, 700 μm, and 2 mm are simulated and plotted in Figure 19(a). Note that the Weber and Reynolds numbers are computed using fluid properties, relative velocity, and drop size. By holding the fluid properties and relative velocity constant, the only way to produce different We/Re combinations is to model different drop diameters. Therefore, for this range of drop diameters, the Ohnesorge number varies from 0.0370 to 0.00443. Similarly, the results of 26 different water drop collisions with an impact parameter of 0.9 are provided in Figure 19(b).Figure 19
Weber versus Reynolds number. Water drop collisions. (a)X=1, boundary between coalescence and separation is depicted using a power fitted line. (b) X=0.9, in this case the boundary between coalescence and separation is shown using a spline fitted line.
(a)
(b)A clear boundary between coalescence and separation is observed for both impact parameters. For the grazing collisions (X=1), the boundary was produced by fitting a power function between the data points. This boundary indicates that as the Re number approached zero, the We number approached We→60. This implies that all collisions with this impact parameter and a Weber number greater than 60 would result in separation. By lowering the impact parameter to 0.9 (and thus creating a region of interaction), a minimum Reynolds number between 126 and 100 is observed as the We number approaches zero. Therefore, this interacting region (which allows portions of each drop to mix) increases the viscous forces while reducing the surface tension forces, such that permanent coalescence will occur. Also, the We number approaches 22 for X=0.9 and 18 for X=1 as the Re number is increased. For head-on collision, it was shown that as the Re number increases, the We number approached We→35 [40]. Therefore, the Reynolds number effect is only important for binary drop collisions with a Weber number less than 35 for X=0 and decreases as the impact parameter is increased.


## 4.1. Internal Flow Patterns

There has been a lot of debate about the existence of internal circulations in droplets of a spray. Since it is very difficult to see the internal flow of a moving droplet, the main way to determine this effect has been through numerical modeling. The present study puts forward a concept for the formation of the internal circulations in the droplets. Figure7 shows the internal flow patterns for the drop collision of Figure 6 (We = 60, Re = 1751, and X = 0.7) and for the plane described in Figure 3. The first image in Figure 7 shows the two drops just after contact. Immediately after the collision a small region of interaction is produced. In this region the fluid is redirected away from the collision point. As this occurs, the outer (or noninteracting) portion continues to translate unimpeded in the original direction (as shown in images 2–4). In the fifth image, in the center of the combined mass some of the flow begins to be redirected “radially” while the bulk of the motion still occurs in original translating direction. However, surface tension compensates and the radial flow is redirected, thus causing the overall flow to be directed upward/downward and away from the center. This action causes the stretching phenomena observed in off-axis collisions. The bulk of the fluid eventually ends up in the bulbous ends as indicated by the streamlines in images 7–9. Image 10 indicates that separation has occurred with the production of two parent drops and a ligament. It is interesting to note several observations at this point. First, review of early images does not provide any evidence of the pinching mechanism in this defining plane. Therefore, the pinching phenomena are the result of radial type flow. Second, the parent drops exhibit flow patterns representing translating flow. Unlike the scenario for the head-on collisions [40] these parent drops reflect similar flow patterns of the original two drops which causes them to continue to translate after separation occurs. Lastly, the ligament collapses into a satellite drop producing the symmetrical flow pattern observed. Therefore, main droplets after collision have unidirectional internal flows, whereas satellite droplets may have multidirectional and internal flow circulations.Figure 7
Internal flow patterns for they-z plane. Results are for a binary water drop collision with a diameter of 700 μm and We = 60, Re = 1751.


## 4.2. The Upper Boundary between Coalescence and Separation

The upper boundary between a coalescence and separation collision occurs for grazing collisions, that is,X = 1. As soon as the surfaces of two drops touch, the liquid flows between them. Depending on the collision conditions, the two drops may touch and temporarily coalesce, and then separate, or remain coalesced. The We number corresponding to this condition is referred to as the upper critical We number, Wecr. Table 1 shows the numerically determined upper critical Weber number for equally sized water drops as a function of drop diameter. It is clear that the upper critical Weber number increases with decreasing diameter. The pressure inside smaller drops is larger than those in larger drops; therefore, the flow is pushed into the liquid bridge between the two drops faster in the smaller drops than for the larger drops. The larger the mass of the bridge, the more difficult it is to break the drops. Therefore, the separation occurs at higher Weber numbers for the smaller drop collisions. The curve fit to the data of Table 1 is Wecr=-2.7lnD+36. As the drop size increases, the effect of change in drop size becomes less important. The upper critical We number approaches We = 18 for large drop sizes.Table 1
Upper critical Weber numbers for two equal size water droplets (X=1).
Diameterμm
We
cr
Re
10
32.5
154
75
22
347
150
21
480
300
20.3
665
700
18.3
966
2000
18
1622


## 4.3. Stretching Separation Boundary

The boundary between permanent coalescence and separation for off-axis collision of two equal size droplets for different collision impact parameters is obtained. Figure8 presents the results of this investigation for equal size drop collisions with a diameter of 300 μm (Figure 8(a)) and 700 μm (Figure 8(b)). Outcomes of permanent coalescence, separation without satellite production, and separation with one satellite produced are plotted for a range of Weber numbers and impact parameters. Also included in the figures are the boundaries defined by the analytical models proposed by Ashgriz and Poo [22], Arkhipov et al. [12], Brazier-Smith et al. [7], and Park [5]. The model by Ashgriz and Poo [22] is referred to as AP/EXP model for the remainder of the discussion. Similarly, the other models are referred to as A, B-S, and P, respectively. The results show that, for a constant impact parameter, as the Weber number is increased, the outcome changes from coalescence to separation without satellites to separation with a satellite. Similarly, as the Weber number is held constant, the same progression is observed provided there is enough energy initially present in the drops (i.e., at low Weber numbers, regardless of impact parameter, coalescence will always occur).Figure 8
Coalescence-fragmentation curves: (a) 300μm drop collisions and (b) 700 μm drop collisions. Symbols represent numerical data while lines represent boundary models presented by previous researchers.
(a)
(b)The numerical division between the coalescence and separation occurs at higher Weber and impact parameter combinations than those observed by the AP/EXP and B-S models. A possible explanation for this behavior is attributed to the existence of trapped air in the gap between the colliding drops. In the experiments, the impact parameter is determined by measuring the lines of trajectory and the corresponding separating angle [22]. Depending upon when this is done (i.e., a few diameters from the drop generators or just prior to impact) the value of the impact parameter can be quite different. As the two drops travel through air, their trajectories remain relatively constant. However, as they come closer and closer, the surrounding gas becomes trapped in between them, which results in the slippage of the drops on each other. This may alter the original impact parameter based on the drop trajectories prior to collision. The effect of such a slippage is shown in Figure 9 by a shift in data. In this figure, the solid line represents the boundary using the AP/EXP model while the symbols represent the numerical results for the stretching separation without satellite production. The dashed line represents a linear translation of the AP/EXP model and is meant to convey the concept of impact parameter slippage. For example, at We = 31, the impact parameter on the AP/EXP boundary is X = 0.44, while that numerically is Xs = 0.64. Thus, due to the presence of the trapped air, the impact parameter at the time of collision could actually be higher than that previously thought.Figure 9
Impact parameter deviation for the same Weber numbers.X is a point on the experimental boundary, while Xs represents the impact parameter observed numerically. The solid line represents the boundary model by Ashgriz and Poo [22], while the dashed line represents a linear translation added to their model.The presence of the trapped air has three different effects on the drops, which ultimately cause an increase (or slippage) in the effective impact parameter. These effects are (i) altered trajectories, (ii) altered internal flow fields, and (iii) shape deformation. The trapped air conveys an external force on the drops, which may alter their trajectories, as shown in Figure10. The second effect is shown in Figure 11, where drops moving in air may develop internal flows. The internal flow changes the drop dynamics after collision and may change the collision outcome. Finally, the surrounding air may cause drop deformation, as shown in Figure 12. As a result of this deformation the region of contact may increase as opposed to the case in which the drops stay perfectly spherical. This increased contact region may result in a larger interaction region and thus more mixing. It is, therefore, concluded that, in determining the collision impact parameter, one has to be worried about the slippage effect that may occur due to the air entrapment.Figure 10
Effect of trapped air between colliding drops. Experimental photo from Ashgriz (unpublished material from Ashgriz and Poo [22]). The trapped air alters the collision impact parameter based on droplet trajectories.Figure 11
Internal flow of two drops approaching each other in (a) vacuum and (b) in air.
(a)
(b)Figure 12
Two cases show the effect of trapped air in deforming the drop prior to collision and merging [unpublished photos from Ashgriz and Poo experiments [22]].
(a)
(b)Figure8 also compares the numerical results with models by Arkhipov et al. (A) [12] and Park (P) [5]. Arkhipov assumed that after the collision a spherical drop is formed that rotates with a constant angular velocity, while Park assumed that the combined mass rotates as two connected spheres. The numerical results, as evident in Figure 5, show that, for the drops to separate after collision, they rotate as two drops that are connected, and not as a combined spherical mass. Therefore, Park’s model is closer to numerical observation of the collision dynamics. However, the drop shapes and masses are altered and are not exactly the same as the original droplets. This becomes particularly true as the Weber number is increased. However, at low Weber numbers, We < 25, the numerical results do not correlate with any of the models considered. This is most likely explained by an incorrect assessment of the amount of angular velocity produced during the evolution. One could imagine that the amount of angular velocity would in fact not be constant (as assumed by Arkhipov et al.) and be a function of the impact parameter instead.


## 4.4.  Deformation

The deformation variations in each direction for an off-axis collision of two water drops with a Weber number of 30 and a Reynolds number of 810 are provided in Figures13–16. For clarity, the deformations are grouped according to either low impact parameters, X<0.5, or high impact parameters, X>0.5. It is convenient to nondimensionalize the amount of deformation by dividing by the drop’s initial diameter. In doing so, the following definitions are used to represent the dimensionless deformation in each direction: (5)η=LxD,γ=LyD,ψ=LzD,where, Lx is the amount of deformation in the x-direction, Ly is the amount of deformation in the y-direction, and Lz is the amount of deformation in the z-direction. Since off-axis collisions occur on different trajectories, the evolution process occurs in all three dimensions, resulting in separate behavior for each direction. In the next three subsections, the behavior in each direction is analyzed.Figure 13
x-direction deformation for We = 30 and Re = 810. (a) Low impact parameters, X < 0.5, and (b) high impact parameters, X ≥ 0.5.
(a)
(b)Figure 14
y-direction deformation for We = 30 and Re = 810. (a) Low impact parameters, X < 0.5, and (b) high impact parameters, X ≥ 0.5.
(a)
(b)Figure 15
z-direction deformation for We = 30 and Re = 810. (a) Low impact parameters, X < 0.5, and (b) high impact parameters, X ≥ 0.5.
(a)
(b)Figure 16
Maximumz-deformation for various Weber numbers and impact parameters. (a) 300 μm drops and (b) 700 μm drops. Solid symbols, coalescence, and open symbols, separation. Boundary line, ψmax = 0.006 We + 3.
(a)
(b)As observed with head-on collisions [40], three major periods of the collision process are spreading and production of a disk, period I, contraction of the disk, period II, and stretching along the collision axis (z-direction), period III.


## 4.5.x-Direction (Minor Axis 1)

Figure13 contains the dimensionless deformation in the x-direction. Review of this figure during period I provides the following information. First, it is observed that, at t∗=0, η increases as the impact parameter, X, increases. This behavior is a direct consequence of the impact parameter definition, since the drop’s separating distance is defined in this direction. Second, as the impact parameter is increased the amount of spread decreases. Therefore, as the region of interaction becomes smaller due to the increased impact parameter, less fluid is allowed to flow in this direction. Additionally, at X~0.6, no spread is observed and the deformation actually becomes drainage or contraction for higher impact parameters. This would seem to imply that, at this (X) and higher impact parameters, fluid is actually pulled out of this direction. The amount of time this period takes to occur is relatively constant regardless of impact parameter.Period II was defined to represent the action of fluid contraction of the disk formed during period I. When the variation of impact parameter is considered, the following observations can be made. First, the amount of contraction is relatively constant until a value ofη~1.3 is obtained at the end of this period for X<0.3. Essentially, this means that, for low impact parameters, the redirected fluid momentum in this direction dominates the flow which in turn causes the disk to contract inward. Second, as the impact parameter is increased until X~0.6, more contraction occurs until a minimum value of η~0.75 is reached. In this range of 0.3<X<0.6, a competition between the amount of fluid momentum and the surface tension effect becomes evident, with the fluid momentum effect decreasing as the impact parameter is increased (i.e., due to the reduction in the interaction region). For impact parameters greater than 0.7, this minimum deformation value is maintained. Therefore, at these high impact parameters, the amount of deformation does not change since the surface tension is now the controlling factor.


## 4.6.y-Direction (Minor Axis 2)

The dimensionless deformation in they-direction is provided in Figure 14. As with the x-direction variation, period I behavior still exists regardless of impact parameter variation. However, the amount of spread decreases as the impact parameter increases. Again, this is seen as a consequence of the reduction of the interaction region. Unlike in the x-direction, spread is observed up until X~0.8 (instead of X~0.6 as seen in the x-direction). At higher impact parameters, the amount of deformation shifts to contraction instead. This indicates that more mixing occurs in the y-direction than in the x-direction, since the fluid is still able to spread outward for higher impact parameters.During period II the overall amount of contraction remains constant (since similar slopes are observed) for low impact parameters ofX<0.3 (i.e., γ at the end of this period decreases from 0.85 to 0.7). Therefore, even though more deformation is observed at the end of this period, the control parameter is still fluid momentum as it is for the x-direction case. The smaller γ values are the consequence of smaller γmax values observed during period I. For impact parameters in the range of 0.4–0.6, the amount of deformation at the end of this period is constant and has a value of γ~0.6 indicating a shift in dominance from fluid momentum to surface tension. As the impact parameter is further increased, X>0.7, less overall deformation occurs since the amount of initial spread seen during period I has decreased and surface tension effect is in control.


## 4.7.z-Direction (Major Axis)

In general, as the impact parameter increases, the region of interaction between the two drops decreases causing the intensity of deformation to shift from thex- and y-directions to the z-direction. This is clearly observed in the z-deformation curves of Figure 15. For an impact parameter of 0, the maximum deformation obtained is 2.7 times the drop’s initial diameter and occurs at tc∗=6.7. As the impact parameter is increased the maximum height obtained decreases and occurs earlier in time until an impact parameter of 0.3 is reached. At this impact parameter the maximum deformation is 1.8 times the diameter and occurs at tc∗=4.1. As the impact parameter is increased further the maximum deformation increases linearly until ψmax=3 times the diameter is obtained. As with head-on collisions, this value defines a limiting value in which either permanent coalescence or separation occurs. This is observed at impact parameters of 0.8 and 0.9, in which the combined mass separates by the pinching mechanism explained earlier. In the case of the impact parameter of 0.7, the contraction rate is greater than the pinching mechanism and stretching-coalescence is observed.The behavior of this maximum deformation parameter,ψmax, as a function of impact parameter and Weber number is plotted in Figure 16. Figure 16(a) represents the results for collisions with drop diameters of 300 μm and Figure 16(b) is for collision with diameters of 700 μm. Open symbols indicate the collision outcome is coalescence while solid symbols indicate separation collision. Review of this figure indicates several aspects. First, the behavior of ψmax is linear over the range of Weber numbers simulated and for illustrative purposes a straight line is fitted through each set of impact data. The slope of each impact parameter data set, m, is computed and presented in Table 2. The subscripts 300 and 700 refer to the drop size. Also included in the table is the change in slope, δm, which indicates the difference in the slope as the drop size is decreased from 700 μm to 300 μm.Table 2
Slopes of lines in Figure16.
X
m
D
=
300
μ
m
m
D
=
700
μ
m
δ
m
0.4
0.018
n/a
—
0.5
0.038
0.046
0.008
0.6
0.043
0.051
0.008
0.7
0.058
0.063
0.005
0.8
0.073
0.078
0.005
0.9
0.088
0.087
0.001
0.95
n/a
0.096
—Regardless of drop size, the slope of the lines increases with the increase in the impact parameter. This is attributed to the reduction in the interaction region which in turn causes the bulk of the fluid motion to continue traveling in the originalz-direction. It is also interesting to note that the drop size has a small impact on the final stretch length. This is evident by the increased slope values for the 700 μm drop collisions as opposed to 300 μm drop collision occurring at the same impact parameter. Additionally, the slopes for X=0.9 are very similar and as a result this indicates that the size effect is no longer a factor in the overall stretching process.A last comment concerning Figure16 concerns the boundary line (represented by a thick dashed line) indicating the regions of coalescence and separation. According to Lord Rayleigh’s capillary instability model [44], this boundary has a value of ψmax=π. However, for off-axis collisions, this boundary was empirically determined to be(6)ψmax=0.006We+3,where the slope reflects the impact parameter effect.


## 4.8. Dissipation Based Models

At the beginning of the collision, both drops have a known amount of surface,SEo, and kinetic energies, KEo. These energies are computed as follows:(7)SEo=2πDi2σ,KEo=π24ρDi3Vr2,where the ratio of the kinetic energy to the surface energy is defined as We∗ = We/48 [22]. Each of these energies is computed at each time step using the following definitions:(8)SEn=σ∑p=1NcellsSAp,KEn=12∑ijk=1Ncellsfijkρvcelluijk2+vijk2+wijk2,where SAp is the polygon area in the computational domain and n is the time index, Ncells is number of computational cells, fijk is volume fraction of the liquid in the cell, and uijk,vijk, and wijk are the velocity components in each direction. Therefore, using these relationships each energy variation is computed at each time step during the collisional process. Additionally, by using these energies and conservation of energy an estimation of the change in total energy or dissipation energy, ϕ, is obtained as follows:(9)ϕ=KEo+SEo-KEn-SEn.The kinetic and surface energies are scaled with the initial surface energy as α=KEn/SEo and β=SEn/SEo. The initial surface energy is chosen as the scaling energy since for a same-size drop collision, it is independent of the impact energy. Since the amount of dissipated energy is simply the combined transient kinetic and surface energy subtracted from the initial total energy, a nondimensional dissipation energy is defined as(10)Φ=1-KEn+SEnKEo+SEo.


### 4.8.1. Kinetic Energy

The dimensionless kinetic energy variation for the binary collision of two water drops for We = 30 and Re = 810 is provided in Figure17. As the impact parameter is increased, the maximum kinetic energy, αmax, decreases from 0.25 for X=0 to 0.22 for 0<X<0.6. It should be noted that, for this range of impact parameters, the final outcome was coalescence, while higher impact parameters (X>0.7) resulted in separation. This maximum kinetic energy decrease is the result of more energy conversion to surface energy (due to the increased surface area accompanying an increase in impact parameter). Additionally, whereas head-on collisions exhibit defined peak and valley behavior corresponding to times of maximum spread and contraction, an increase in the impact parameter causes the change in kinetic energy to undergo a more chaotic variation. It is not until the impact parameter is increased to a value greater than 0.6 that no oscillatory behavior (of any kind) is observed. Therefore, if the kinetic energy does not peak early in the evolutional sequence, separation will likely occur.Figure 17
Dimensionless kinetic energy,α = KEn/SEo for We = 30 and Re = 810. (a) Low impact parameters, X < 0.5, and (b) high impact parameters, X ≥ 0.5.
(a)
(b)


### 4.8.2. Surface Energy

Figure18 provides the dimensionless surface energy variations for the binary collision case under consideration. The maximum surface energy, βmax, decreases from a value of 1.325 for X=0 to 1.025 for the range, 0<X<0.6. This behavior parallels the behavior seen for the kinetic energy. At impact parameters greater than 0.6, the absence of this peak behavior at this time in the evolution is conspicuous. In fact at this time, βmax~1 and as noted earlier in the previous discussion, all the remaining collisions at higher impact parameters resulted in separation. Additionally, the variation in surface energy becomes more oscillatory or chaotic for impact parameters up to 0.6. At these low and intermediate impact parameters, various surface deformations occur as the combined mass tries to stabilize itself to its minimum energy state.Figure 18
Dimensionless surface energy,β = SEn/SEo for We = 30 and Re = 810. (a) Low impact parameters, X < 0.5, and (b) high impact parameters, X ≥ 0.5.
(a)
(b)


## 4.8.1. Kinetic Energy

The dimensionless kinetic energy variation for the binary collision of two water drops for We = 30 and Re = 810 is provided in Figure17. As the impact parameter is increased, the maximum kinetic energy, αmax, decreases from 0.25 for X=0 to 0.22 for 0<X<0.6. It should be noted that, for this range of impact parameters, the final outcome was coalescence, while higher impact parameters (X>0.7) resulted in separation. This maximum kinetic energy decrease is the result of more energy conversion to surface energy (due to the increased surface area accompanying an increase in impact parameter). Additionally, whereas head-on collisions exhibit defined peak and valley behavior corresponding to times of maximum spread and contraction, an increase in the impact parameter causes the change in kinetic energy to undergo a more chaotic variation. It is not until the impact parameter is increased to a value greater than 0.6 that no oscillatory behavior (of any kind) is observed. Therefore, if the kinetic energy does not peak early in the evolutional sequence, separation will likely occur.Figure 17
Dimensionless kinetic energy,α = KEn/SEo for We = 30 and Re = 810. (a) Low impact parameters, X < 0.5, and (b) high impact parameters, X ≥ 0.5.
(a)
(b)


## 4.8.2. Surface Energy

Figure18 provides the dimensionless surface energy variations for the binary collision case under consideration. The maximum surface energy, βmax, decreases from a value of 1.325 for X=0 to 1.025 for the range, 0<X<0.6. This behavior parallels the behavior seen for the kinetic energy. At impact parameters greater than 0.6, the absence of this peak behavior at this time in the evolution is conspicuous. In fact at this time, βmax~1 and as noted earlier in the previous discussion, all the remaining collisions at higher impact parameters resulted in separation. Additionally, the variation in surface energy becomes more oscillatory or chaotic for impact parameters up to 0.6. At these low and intermediate impact parameters, various surface deformations occur as the combined mass tries to stabilize itself to its minimum energy state.Figure 18
Dimensionless surface energy,β = SEn/SEo for We = 30 and Re = 810. (a) Low impact parameters, X < 0.5, and (b) high impact parameters, X ≥ 0.5.
(a)
(b)


## 4.9. Reynolds Number Effect

The Reynolds number effect is presented based on the collision of two drops with impact parameters ofX=1 and X=0.9. For the grazing collision (X=1), 25 different water drop collisions with drop sizes of 10 μm, 75 μm, 150 μm, 300 μm, 700 μm, and 2 mm are simulated and plotted in Figure 19(a). Note that the Weber and Reynolds numbers are computed using fluid properties, relative velocity, and drop size. By holding the fluid properties and relative velocity constant, the only way to produce different We/Re combinations is to model different drop diameters. Therefore, for this range of drop diameters, the Ohnesorge number varies from 0.0370 to 0.00443. Similarly, the results of 26 different water drop collisions with an impact parameter of 0.9 are provided in Figure 19(b).Figure 19
Weber versus Reynolds number. Water drop collisions. (a)X=1, boundary between coalescence and separation is depicted using a power fitted line. (b) X=0.9, in this case the boundary between coalescence and separation is shown using a spline fitted line.
(a)
(b)A clear boundary between coalescence and separation is observed for both impact parameters. For the grazing collisions (X=1), the boundary was produced by fitting a power function between the data points. This boundary indicates that as the Re number approached zero, the We number approached We→60. This implies that all collisions with this impact parameter and a Weber number greater than 60 would result in separation. By lowering the impact parameter to 0.9 (and thus creating a region of interaction), a minimum Reynolds number between 126 and 100 is observed as the We number approaches zero. Therefore, this interacting region (which allows portions of each drop to mix) increases the viscous forces while reducing the surface tension forces, such that permanent coalescence will occur. Also, the We number approaches 22 for X=0.9 and 18 for X=1 as the Re number is increased. For head-on collision, it was shown that as the Re number increases, the We number approached We→35 [40]. Therefore, the Reynolds number effect is only important for binary drop collisions with a Weber number less than 35 for X=0 and decreases as the impact parameter is increased.


## 5. Conclusion

An investigation of off-axis collisions of water drops across various impact parameters is presented. The results show that as the impact parameter is increased to values greater than zero (off-axis collision), the region of interaction between the drops is reduced. Two different types of break-up mechanisms are observed. The ligament holding the two drops may go through pinching process governed by the Rayleigh instability, or the ligament may be pulled until it is broken. The pulling (tearing) breakup occurs at high impact parameter and high Weber number combinations. In these situations, the region of interaction is insufficient to allow adequate mixing between the two colliding drops. Therefore, the unmixed portions of the drops continue their travel in their respective original directions, thus tearing the ligament holding the two drops.A study of the internal flow field of droplets shows that the droplets generated after collision may have different types of internal flows. Collision with high impact parameters may result in main drops and satellite drops. The internal flow field of the main drops seems to be similar to those of the parent drops, generally unidirectional. However, the internal flow field of the satellite drops may have internal circulations and may be multidirectional. Therefore, the internal flow field of droplets in dilute sprays, which do not have significant droplet collisions, may be unidirectional. However, the internal flow fields of dense sprays, which may have significant number of high impact droplet collisions and formation of satellite droplets, may be multidirectional and have internal circulations.The upper critical Weber number, describing the boundary between the coalescence and separation for grazing collisions (X=1), is determined for water drop collisions with drop diameters ranging from 10 μm to 2 mm. These critical Weber numbers are found to be a function of the drop diameter. The smaller drops have higher We number values, that is, Wecr=32.5 for D=10 μm, whereas Wecr=18 for D=2 mm. This is due to the higher pressures that exist in smaller drops, which quickly form a liquid bridge between the drops as soon as they touch. Smaller drops form relatively larger bridges, which requires higher impact We numbers to cause separation. However, as the drop size increases, the effect of change in drop size becomes less important. The upper critical We number approaches We = 18 as the drop size is increased. It should be noted that the above conditions are relevant for cases in which the surrounding gas is dynamically neutral, or collisions occur in vacuum. Otherwise, in determining the collision impact parameter, one has to be worried about the slippage effect that may occur due to the air entrapment in between the two drops.A deformation criterion based on the drop deformation after collision is obtained. This criterion indicates that the drops separate if the length of the ligament formed after collision is approximately 3 times the diameter of the original drop. The maximum length of the ligament formed after the collision is found to vary linearly with the impact parameter.Plots of the collisional outcomes in the Re-We plane are provided for the grazing collisions (X=1) and collisions with an impact parameter of X=0.9. As the Weber number increases, the Re for the separation approaches zero for an impact parameter of 1 and approached Re → 100 for collisions with an impact parameter of 0.9. It was shown [40] that Re → 1000 for an impact parameter of X=0. Therefore, it can be concluded that collisions below this critical Reynolds number will always result in coalescence regardless of the Weber number. Also, the Reynolds number does have an effect on the collisional outcome as the impact parameter is reduced. However, this latter influence is small for drops with diameters larger than 75 μm and for impact parameters of 1 and 0.9.



---

*Source: 405696-2015-05-25.xml*