1. Introduction
Studies of droplets squeezing through small or narrow constrictions have been extensively performed for various applications, such as microfluidics or biomedical engineering [
1]. Furthermore, extensive research has delved into the movement of biological cells navigating narrow passages within the human body, encompassing phenomena such as our previous modeling work on the filtration of red blood cells within the spleen [
2], the process of white blood cell diapedesis [
3], and the intricate journey of cancer cells during metastasis [
4]. In our recent work, we conducted extensive microfluidic experiments and combined them with multiscale simulations to understand the critical conditions for red blood cells’ passage through splenic inter-endothelial slits (IESs) and the physical mechanisms that control their transit dynamics [
2]. Although our computationally predicted transit time matches well with our experimental measurements, the simulations are expensive and less insightful than analytical theory. Additionally, in another recent work, we developed an analytical theory for a droplet passing through a circular pore [
5], but the idealized circular pore is significantly different from the IES geometry we studied experimentally and computationally [
2]. In this study, we extended our analytical theory of a circular pore to a slit for the transit time of a droplet.
In general cases under both in vivo and in vitro conditions, the speed at which droplets or cells pass through constrictions is a critical factor, yet its quantitative characterization remains unclear. For instance, red blood cells’ transient slowdown as they transit through splenic IESs can lead to capture and subsequent destruction by macrophages of the immune system. An in vitro example is that of droplets flowing through microfluidic devices, where droplet velocity is also important for effective manipulation and separation processes. Another instance where flow velocity is important is in the flow of ferrofluids into ducts of circular or rectangular cross-sections, where additional parameters like pressure drop, constriction cross-section size, and curvature can significantly impact the efficiency of heat exchangers and mixers [
6].
Numerous experimental studies have been conducted to understand the dynamics of droplets, vesicles, and cells as they transit through constrictions. In 1967, Gregersen et al. investigated the passage of red blood cells through small pores created in paper [
7]. The development of versatile microfluidic technology 20 years or so has vastly expanded the opportunities to explore the behavior of soft objects squeezing through constrictions of various size and geometry under controlled flows. Examples of microfluidic-based experiment are Ma et al., who studied the flow patterns within droplets navigating through rectangular microchannels [
8], and Wang et al., who discussed the dynamic behavior of viscoelastic droplets in Y-shaped capillary channels [
9]. Gambhire et al., followed by Moreau et al., examined red blood cells as they pass through submicron-wide slits designed to mimic the splenic IESs, highlighting the crucial role of red blood cells’ mechanical properties for efficient transit [
2,
10]. Ren et al. created a microfluidic device to study the biophysical characteristics and transit times of cells navigating cyclically through constrictions [
11]. Finally, to relate transit time and object viscosity, Khan et al. studied cancer cells and viscous droplets flowing in confining microchannels [
12]. Still, not all experiments are performed at the micrometer scale. For example, Chen et al. sought to discover the effects of droplet size (in the millimeter range), constriction dimensions, and applied pressure drop on the trajectory of droplets in irrigation and agricultural applications [
13]. Nevertheless, despite all this extensive work, how the transit dynamics of these soft objects precisely relate to the setup parameters (applied pressure drop or constriction geometry) and object intrinsic properties remains unclear.
Apart from the above experimental approaches, the transit of soft objects through small apertures has been studied both from a theoretical and a computational point of view. First, theoretically, Zhang et al. developed analytical models for investigating the pressure and minimum impulse associated with a droplet’s passage through a circular pore under a constant flow rate [
14,
15,
16]. However, in many applications, pressure control takes precedence over flow rate regulation. In addition, Jensen et al. employed an energy-based approach to examine the behavior of a bubble as it squeezes through a short constriction [
17], while Marmur conducted an analysis of droplet penetration through a capillary under the influence of gravity [
18]. In the case of vesicles, Gompper and Kroll developed an analytical solution to model the mobility of vesicles as they squeeze through a cylindrical pore past a threshold driving field strength [
19]. When considering red blood cells, Waugh et al. formulated an analytical model to estimate transit times during egress in bone marrow [
20]. When studying Stoke’s flow, Dassios and Vafeas developed a 3D concentric sphere model for particles in creeping flows and derived analytical expressions for the velocity, the total pressure, the angular velocity, and the stress tensor fields [
21]. In addition, we recently developed an analytical study of droplet transit through a circular pore [
5].
In addition to analytical models, a wide array of numerical methods has been harnessed to study the passage of droplets, vesicles, and cells through microfluidic pores, as reviewed by Zhang et al. in 2014 [
22]. For instance, Barthes-Biesel made pioneering contributions by employing boundary integral simulations to investigate the passage of vesicles and capsules through constrictions in an axisymmetric configuration [
23]. Similarly, Zinchenko and Davis applied a similar approach to explore the three-dimensional scenario of a droplet passing through constrictions formed between spheres [
24].
Despite extensive prior research, understanding how the transit time of an object through a constriction changes with the constriction dimensions or the applied pressure drop and with the object intrinsic properties such as viscosity or surface tension remains unclear. Existing studies, such as those by Zhang et al. [
14,
15,
16], often assume constant flow rates, while real-world applications frequently involve constant pressure drops. Additionally, microscale applications primarily exhibit creeping flows instead of finite Reynolds number flows [
14,
15,
16]. Moreover, real-world constrictions often have rectangular cross sections rather than circular ones. In this study, we developed analytical models for droplets passing through slit-like constrictions under constant pressure conditions in order to derive the exact solution of transit time. Our work has two novelties: first, it is an analytical study, rather than a computational study, where it is rather trivial to change the constriction geometry; second, from the application perspective, the slit configuration is much more closely related to real in vivo situation like that of IESs, and our analytical approach gives a much more precise relationship between transit time and slit dimensions than a circular pore theory. Altogether, our results give much more insight into explaining the dynamics of cells passing through slits, which we previously observed in our experiments and simulations [
2].
2. Theory of a Droplet Squeezing through an Infinitely Wide Slit
Our model seeks to mathematically describe the forces acting on the droplet in each stage and at different locations.
Figure 1a describes the position of the droplet at the boundaries of different stages and
Figure 1b labels pressures
P1 to
P7 at different locations of the system during stage II. We take into account the lubrication layer between the droplet and the slit walls, and we consider the droplet and lubrication layer flows inside the slit to be planar Poiseuille flow and Couette flow, respectively. Forces due to surface tension are calculated at the spherical and/or elliptical ends of the droplet. The flow entering and exiting the slit is estimated with Roscoe’s solution with modifications to address the presence of multiple phases. The parameters that are used in the equations include the dimensions of an infinitely wide slit of thickness
T and length
L, the viscosity of the droplet
ηd and external fluid
η0, the surface tension
and volume
of the droplet, and the pressure drop
that drives flow through the slit.
2.1. Contributions to Pressure Drop
We aimed to calculate the total transit time when the prescribed total pressure drop
is constant during the transit process. The total pressure drop
between the two sides along the dashed red path in
Figure 1b can be grouped into four terms as:
where the first term
is the pressure term corresponding to pressure drops across the droplet membrane on both sides of the droplet. The second and third terms are the Poiseuille-related pressure drops, and the last term is the Roscoe flow term of a flow passing through an infinitely thin slit [
25].
We assume that the width of the droplet inside the slit is constant, the same as the initial diameter of the droplet, given by
. The flow rate
Q due to pressure differential ∆
Pηi and viscosity
η inside the slit is represented by:
which is the flow created between two infinitely long parallel plates, or planar Poiseuille flow. The equation can be rearranged to express the pressure drop due to viscous friction inside the slit:
The pressure drop due to the viscous friction in the spherical part outside the slit can be found using Roscoe’s solution for the flow through a slit (Roscoe’s extension of Sampson flow [
25]).
Thus, the total pressure drop of a droplet in an infinitely wide slit is represented as:
2.2. Relationships between Pressure Drop and Lubrication Layer, Velocity Profile, and Flow Rate
As the droplet progresses through the slit, it does not touch the slit walls. Instead, the external fluid forms a lubrication layer between the droplet and each wall as shown in
Figure 2, where
h is the thickness of each layer. For a 2D droplet in a Hele–Shaw flow in a channel, it was shown that
h is related to the channel height
T and to the capillary number
Ca as:
Let
be the velocity of the droplet at the droplet–lubrication layer interface and
be the maximum velocity of the droplet. The shear stress is continuous across the interface, so that:
Note that this is exactly the same as the circular pore case (see Equation (7) in [
5]).
To prepare to derive an analytical solution later, we define:
where
for ideal Poiseuille flow and
for ideal plug flow.
Let
Q1 represent the flow rate of the droplet and
Q2 represent the flow rate of a single lubrication layer. We calculate the flow rates as:
Because in a planar Poiseuille flow, the average velocity is 2/3 of the maximum velocity, we have:
This leads to our final equation, giving the total pressure drop
as function of the droplet velocities and flow rates:
where
and
are the mean curvatures of the right and left ends of the droplet, respectively. Like the Sampson term shown in our previous work [
5], the Roscoe term can be split into two terms of two half problems (the last two terms in Equation (14)). Equations (7) and (14) can be solved for
and
.
2.3. Transit Times
Here, we determine the durations of the droplet transit stages and of its total transit through the slit, based on the pressure drops calculated above. To do so, we adapted the procedures used in the circular pore case [
5]. For stages I and V, the time needed for the droplet to form or retract a curved half-ellipsoidal shaped droplet head is very small, and so they are considered negligible for this study.
During stage II, the constant total pressure
tot is:
where
when
is approximated in a rectangular slit with dimensions
and
[
26] and
. After solving Equation (15) with Equation (7) together, we obtain:
This nonlinear ODE can be solved by integrating from the initial condition to . The unknown time then can be solved numerically.
The flow rate during stage III can be found using the pressure drop:
where
and
, and
is the volume of the droplet exiting the slit.
The ODE of can be solved with initial condition . By numerically integrating from t = 0 to an unknown time where , we can then numerically solve for .
Since stage II and stage IV are similar,
is found the same as in stage II except each appearance of
is replaced with
.
where
and
. If we denote
, we have:
and
with
This is a nonlinear ODE
which can be solved by integrating from the initial condition of
to
. This unknown time
can then be solved numerically.
The transit time
t1 and
t5 of stages I and V are negligible as these are very short processes. Thus, the total transit time is:
2.4. Relationship between Critical Pressure and Surface Tension
For given values of surface tension, droplet viscosity, and slit dimensions, the critical pressure
c is the minimum total pressure drop for which the droplet is able to transit completely through the slit, which is determined by the slit thickness
T, the droplet volume
, width
W, and its surface tension
σ as:
For given values of pressure drop, droplet viscosity, and slit dimensions, the critical surface tension σc is the maximum surface tension for which the droplet is able to transit completely through the slit, which can be solved in Equation (23) by setting as the given pressure drop value.
2.5. Procedure to Obtain Numerical Solutions for Finite Surface Tension Cases
While an analytical solution can be found when the surface tension term is absent, the differential equations in the above sections cannot be integrated manually when the surface tensions of both the right and left parts of the droplet are considered. In order to find the transit time when the surface tension is non-zero, we used the ODE solver ode45. In cases where ode45 failed to integrate an equation, the ode23 solver was used to perform the calculations needed to generate the figures in our study. We also applied the event function to calculate the transit times when certain conditions are met, such as in Equation (16).
4. Conclusions and Discussion
This paper presents analytical approaches to understand the passage of a droplet through a narrow slit. These approaches can have many applications in fields ranging from understanding the flow of biofuel and ethanol droplets [
27] to improving the properties of agricultural pesticides in emulsion form [
28]. Additionally, this model of droplet flow can be more broadly applied to liquid flows in general, with applications in improving insulation systems [
29] and enhancing the understanding of co-flowing gas-liquid flows [
30]. Compared to our previous work [
5], there are several novelties. First, we found the available analytical solutions to different terms in our current theory of the slit. Notably, we used the analytical theory by Roscoe in 1949 [
25] on flow passing a slit, which was the counterpart of the Sampson theory of a circular pore [
31]. This is very significant in terms of fundamental transport theory in fluid mechanics. Second, from the application perspective, the current work is much more closely related to our recent work on the study of red blood cells passing through IESs in the spleen [
2]. The new result in our current study, such as Equation (28), gives more precise relationships between transit times (total and for each stage) and slit dimensions than a circular pore theory. This will give more insights into explaining the dynamics of cells passing through slits observed in our experiments and simulations, because it explicitly gives relationships which cannot be obtained in numerical simulations. Third, for the surface tension terms, we incorporated the recent result from Darvishzadeh et al. [
26] for applying the Young–Laplace equation to a rectangular opening, which is quite different from applying the Young–Laplace equation to a circular pore. Fourth, for the lubrication thickness, we utilized the result for a droplet inside a Hele–Shaw flow, which is different from a droplet in a circular tube. Finally, besides the parametric studies of transit times, we also showed detailed temporal evolution of the droplet of volumes in
Figure 5, which was not carried out in the previous work.
Although we considered more realistic geometry in this study, there remain some limitations. While our expression evaluates
as a sphere,
in the beginning of stage III is actually a hemispherical cap on the right side. As a result, the expression for
appears larger and our analytical model does not perform as well when the pressure drop approaches
. This effect can be seen as otherwise inexplicable sudden flow rate changes in
Figure 5 not associated with the droplet initially or fully crossing a boundary. The sudden change in curvature calculation between stages appears to be a sharp change in flow rate between stages II and III for
VII and between stages III and IV for
VIII. In addition, the calculation for critical pressure, although it would theoretically be found with the initial conditions of stage II, it is ultimately determined by the initial conditions in stage III due to assuming a spherically shaped droplet end. It will be more accurate to consider a spherical cup to evaluate the left and right spheres in this study; however, this will significantly increase the complexity of the analytical expressions. Because this study focused on total transit time, and the fact that the time where the right and left caps are not approximately spherical is relatively short, our model was sufficiently accurate.