1. Introduction
Liquid droplets are commonly found in nature. At the gas–liquid interface of a droplet, the cohesion of molecules generates the interfacial tension, which has a great influence on the behavior of the droplet. In a gravitational field, an independently existing droplet is subjected to the combined effect of gravity and interfacial tension, and can be suspended at the lower end of a water pipe, showing the shape of an ellipsoid. Under microgravity, the morphology of the droplet is different from that in the gravity field. Because the droplet tends to be subjected to zero gravity, it appears to be suspended in the air in the form of a sphere under the action of interfacial tension. The interfacial tension plays an important supporting role to the droplet model, maintaining it in equilibrium, migration, and motion. Microgravity science focuses on the study of the laws of motion of matter under microgravity (low-gravity) environments and the effects of gravity changing on the laws of motion. Microgravity fluid physics has been focusing on the interfacial tension for a long time to carry out research on the process and mechanism of heat and mass transport, especially on the study of secondary processes that are difficult to realize on the ground, as well as on the fluid physics problems associated with manufacturing of materials in space and biotechnology. One of the typical scientific problems in microgravity fluid physics is the thermocapillary migration of droplets. The equilibrium, deformation, vibration and interaction of immiscible droplets in the mother liquor are affected by the liquid–liquid interfacial tension, and the construction of a droplet migration model is scientifically important for the study of interfacial tension mechanism. With the development of microgravity science, droplet dynamics has gradually become a hot spot of research and has a more profound meaning compared to classical fluid mechanics. The study of the droplet thermocapillary migration mechanism has important application value in industries such as space alloy preparation, glass preparation, multiphase separation, crystal growth, droplet combustion and fire safety.
Young, Goldstein and Block [
1] proposed the YGB theory in 1959, which initiated the theoretical exploration of droplet thermocapillary migration. In 1990, Barton [
2] theoretically analyzed the velocity field and temperature field of droplet thermocapillary migration. For the study of thermocapillary migration of double droplets, Keh [
3] theoretically analyzed the thermocapillary migration process of axisymmetric double droplets to obtain the velocity field and droplet migration velocity in the mother liquor.
For the problem of the thermocapillary migration of droplets, research work mostly focuses on numerical simulation. Brady [
4] analyzed the thermocapillary migration process of three-dimensional deformable droplets and found that the shape of the cell wall ultimately affects the droplet migration process. Loewenberg [
5] calculated the thermocapillary migration velocities of two axisymmetric droplets with different radius ratios. Berejnov [
6] gave the evolution of the center distance between two drops over time and discussed the thermocapillary migration process and the interaction of two drops with different radius ratios. Yin [
7,
8] investigated the thermocapillary migration process of a single droplet with no deformation and discussed the influences of dimensionless parameters and initial conditions on the migration, and also investigated the thermocapillary migration of a two-droplet pair with different radius ratios and the interaction between the two droplets, and analyzed the migration velocities of the large succeeding droplet and the small leading droplet.
Compared with theoretical studies and numerical simulations, there are fewer experiments related to droplet migration. Barton [
9] found that the migration rate is proportional to the droplet radius in the experiments of single-droplet thermocapillary migration, which is the same trend predicted by the YGB theory. Hähnel [
10] observed the process of large droplets overtaking small droplets in the experiments of double-droplet thermocapillary migration. Our project team, Kang et al. [
11] experimentally investigated the thermocapillary migration of two droplets and the interactions between them, and found that the small leading droplets appeared to migrate in an inclined “8” trajectory. In recent years, Kang et al. [
12,
13] observed the migration processes of single and double droplets in real time using digital holographic interferometry. They discussed the thermal wake problem, analyzed the relationship between the droplet thermocapillary migration velocity and the temperature field around the droplet, and revealed the physical law of the same-phase two-droplet thermocapillary migration and interaction; they carried out the experiment of coaxial three-phase two-droplet thermocapillary migration, observed that the small droplet surpassed the large droplet, and explored the influence of the viscosity coefficient between two droplets of different phases on the migration. Vincent and Delville [
14] experimentally studied the thermocapillary migration of the growing droplet in the microchannel by limiting the front advance of the growing droplet by local laser, and found that the laser can control the flow of the droplet in the microchannel, acting as a “soft door”. Vikas Chaurasiya et al. [
15] explored a mathematical model which accounts for the effect of convection in the frozen and vapor regions, addressed gaps in the modeling of sublimation of a humid porous body. For a specific form of the velocity profile, an exact solution of the current problem is obtained via similarity technique.
Based on the above work, in this paper, we construct a physical model of droplet migration in various multiphysics fields, build a synchronized measurement platform of digital holographic interference and particle image velocimetry, and observe the thermocapillary migration process of a droplet in a temperature gradient field in real time, so as to study the mechanisms that involved multiphysics coupling from a completely new viewpoint. Simultaneous measurements of temperature and velocity fields of single-droplet migration, same-phase two-droplet migration and different-phase two-droplet migration have been conducted. Not only has the droplet thermal wake been observed, but also a region with a zero-velocity field at the rear of the droplet, four vortex cell structures, and the off-axis behavior in two-droplet migration have been discovered for the first time. Numerical simulations of droplet migration have been carried out, and the results of the simulations are in line with those of the experiments. The driving mechanism of droplet migration by the thermocapillary force and mother liquor flow is explored, and the coupling correlation between the temperature field and velocity field during the droplet migration process as well as the influence factors on the droplet migration state are analyzed.
3. Results and Analysis of Droplet Migration Experiments
3.1. Principle of Droplet Thermocapillary Migration
The droplet migration problem is a complex process driven by the coupling of multiple physical fields. The droplet is submerged in another immiscible mother liquor, and a temperature gradient is added to the mother liquor, which causes the temperature inhomogeneity at the droplet interface and leads to the inhomogeneity of the interfacial tension, thus driving the flow of the fluids on both sides of the interface, and the droplet will move due to the reaction force of the shear stress of the mother liquor, which is called the thermocapillary migration of the droplet. Previous studies have concluded that in the process of droplet thermocapillary migration, convective cell elements are formed inside the droplet, and the around-flow is formed outside the droplet, as shown in
Figure 8. We have obtained results that are different from this conclusion by simultaneous measurements of the velocity field and the temperature field during the droplet migration process, and have gained a brand new understanding of the process of droplet migration.
The Maragoni number (hereafter referred to as Ma number) represents the ratio of thermal convective transport to thermal diffusion and is defined as follows:
3.2. Single-Droplet Migration
In the experiment of single-droplet migration in a uniform temperature gradient field, the mother liquor is 30 cSt silicone oil and the droplet is a mixture of water and ethanol to ensure density matching. The temperature gradient of the mother liquor Γ = 0.57 K/mm, d = 6 mm, Ar = 0.15, Ma = . The droplet migration temperature and velocity fields are observed in real time using a synchronized measurement platform developed in this study.
The hologram at t = 0 and the hologram at a certain moment in the droplet migration process are collected, and the original holograms are processed to obtain the phase distribution map that reflects the temperature field distribution and the temperature field distribution (unwrapped phase map) at the corresponding moments, as shown in
Figure 9. Throughout the whole process of droplet migration, from the digital holographic interferometric measurements, it is found that there exists a thermal boundary layer near the interface of the droplet, and the droplet tail appears as a thermal wake.
The droplet is surrounded by a low-temperature zone, where the temperature around the droplet is lower than in the mother liquor of the same height, and a thermal boundary layer exists. The main reason is because the effect of convective transport is strong. In the experiment, the Ma number is large, and the thermal convection effect is relatively strong, while the thermal conduction effect is relatively weak. The droplet will migrate to the high-temperature end due to the thermocapillary effect, and the droplet will have a lower temperature than the surrounding mother liquor. However, the heat used to increase the internal temperature of the droplet is mainly supplied by the mother liquor through thermal conduction, rather than through thermal radiation and thermal convection. This heat transfer mechanism results in a relatively weak heat exchange between the interior and exterior of the droplet. This ultimately results in a temperature difference between the interior and exterior of the droplet. A thermal boundary layer exists near the droplet interface to maintain the continuity of heat flux.
The presence of thermal wake can also be observed clearly at the rear of the droplet. The actual temperature isotherms at the rear of the droplet are bent, neither straight nor symmetrical with respect to the horizontal axis passing through the droplet center, which is mainly due to the thermal convective transport. The droplet affects the surrounding temperature field, and the droplet migration bends the isotherms of the originally flat background temperature field. In the vicinity of the front stationary point of the droplet, the isotherms are bent to the high-temperature end, and the isotherms are gathered near the front stationary point of the droplet, which makes the temperature gradient increase in the mother liquor near and above the front stationary point; in the vicinity of the rear stationary point of the droplet, there exists a long thermal wake at the rear of the droplet, and the isotherms are more sparse in the region of the thermal wake below the rear stationary point of the droplet, which makes the temperature gradient smaller in the mother liquor near and below the rear stationary point.
The coupling between the velocity field and temperature field during the droplet migration process is investigated, and the results of the simultaneous PIV measurements are shown in
Figure 10. It is found that the phenomena of boundary layer and wake also exist in the velocity field, and the PIV results are consistent with the digital holographic interferometry results. The complex temperature field in the region of the droplet thermal wake in turn leads to a region with a zero-velocity field at the rear of the droplet, a result that is new to us and has not been reported in previous studies. The temperature field is strongly coupled with the velocity field during the droplet migration process. The interfacial tension drives the flow of mother liquor to generate the velocity field, the velocity field causes the temperature field to change, and the temperature field acts on the droplet migration, resulting in a complex structure of temperature field and flow field. This complex structure has a direct impact on the droplet migration state, affecting the equilibrium of the droplets and leading to migration instability.
In addition to the phenomena of thermal boundary layer, thermal wake, velocity boundary layer and velocity wake, a new understanding of the flow field structure is obtained. Four vortex cells are found to be distributed in the mother liquor during the thermocapillary migration process of a single droplet in the experiments of simultaneous measurements of temperature field and velocity field, as shown in
Figure 11. Two larger vortex cells exist in the mother liquor above the droplet, and two smaller vortex cells exist in the mother liquor in the tail region of the droplet, which is consistent with the results of numerical simulation. This is different from the previous results (see
Figure 8), and the coupling between the temperature field and velocity field in the droplet migration process should be investigated to re-recognize and explore the mechanism of surface tension effect of droplet migration. Previous theoretical studies have mainly investigated the thermocapillary behavior of droplets in an infinite flow field, but in a finite flow field, the boundary conditions of the flow field affect the flow of mother liquor. When the droplet moves vertically upward driven by the thermocapillary force, the mother liquor around the droplet also moves upward with the droplet, and due to the finite flow field and the conservation of mass of the mother liquor, the mother liquor farther away from the droplet moves downward, so a region in the mother liquor with zero velocity appears inside the droplet wake region. And the position of this region of zero velocity continues to move upward in the process of droplet migration. When the region is farther away from the wall of the liquid pool, the horizontal velocity of the mother liquor is no longer zero, and there is vertical motion in the mother liquor in the opposite direction at this time; the vortex cells in the mother liquor begin to form in the region of the droplet tail.
3.3. Same-Phase Two-Droplet Migration
In the experiment of same-phase two-droplet migration in a uniform temperature gradient field, the temperature gradient is 0.57 K/mm. The leading droplets of 7 mm diameter are injected into the mother liquor first, and then the succeeding droplets of 7.8 mm diameter are injected afterwards. Ma
lead =
, and Ma
suc =
. Since the diameter of the succeeding droplet is relatively large, according to the YGB theory, the migration velocity of the droplet increases with the increase in the diameter, so the migration velocity of the succeeding droplet is faster than that of the leading droplet, and the process of the succeeding droplet catching up with the leading droplet can be observed in the experiment. Through the analysis of the velocity field and the corresponding phase field, it is found that the large succeeding droplet will cause the deviation of the small leading droplet from the axis of the succeeding droplet during the process of chasing; when the deviation between the two droplet axes is large, the phase field will be changed from the axisymmetric state to the asymmetric state; the asymmetric temperature field causes the two droplet axes to continue to shift, as shown in
Figure 12. The reason for this phenomenon is that the two-droplet migration has a more complicated temperature field structure and flow field structure, which leads to the oscillation of droplet migration, and the double-droplet migration deviates from the coaxial motion.
From the initial phase field, it can be seen that the initial temperature field has a vertical uniform temperature gradient. When the succeeding droplets are farther away from the leading droplets, the temperature gradient in the mother liquor is still relatively uniform; when the succeeding droplets are about to catch up with the leading droplets, the direction of the temperature gradient in the mother liquor is obviously tilted. The tilted temperature gradient leads to the tilt of the surface tension gradient; as a result, the direction of the surface tension on the droplets will change. Since the diameter of the succeeding droplet is larger than that of the leading droplet, the migration velocity of the succeeding droplet is greater than that of the leading droplet, and it can be seen from
Figure 12a,b when the temperature field appears to be tilted, the velocity field of the mother liquor around the surfaces of the two droplets is also tilted. Regions of zero velocity are found inside the droplet wake region in both the two-droplet migration and single-droplet migration.
3.4. Heterophase Two-Droplet Migration
In the experiment of different-phase two-droplet migration in a uniform temperature gradient field, the temperature gradient is 0.57 K/mm. A leading droplet with a diameter of 8 mm and a glycerol–ethanol mixture (
ν = 8.7) as the working medium is injected into the mother liquor first, and then a succeeding droplet with a diameter of 3.8 mm and a water-ethanol mixture (
ν = 2.8) as the working medium is injected into the mother liquor afterwards. Ma
lead =
, and Ma
suc =
. According to the YGB theory, the droplet motion velocity is proportional to the droplet diameter and inversely proportional to the liquid viscosity. In this working condition, the velocity of the small succeeding droplet is greater than that of the large leading droplet, and the small succeeding droplet will chase and overtake the large leading droplet. It is found in this process that when the small succeeding droplet overtakes the large leading droplet, the moving directions of both the small succeeding droplet and the large leading droplet are tilted, as shown in
Figure 13.
5. Conclusions
In this paper, we construct a physical model of droplet migration in various multiphysics fields, and set up a synchronized optical measurement platform by using digital holographic interferometry and particle image velocimetry based on the development of the fluid physics experimental cabinet on the Chinese space station, in order to observe the thermocapillary migration of droplets in temperature gradient field in real time, and to study the coupling mechanism of multiphysics fields from a new viewpoint. The temperature and velocity fields of single-droplet migration, same-phase two-droplet migration, and different-phase two-droplet migration are measured experimentally, and the mechanism of droplet thermal wake is discussed and analyzed. By measuring the velocity field of droplet migration, (1) a region with a zero-velocity field at the rear of the droplet is firstly found; (2) four vortex structures are firstly observed and verified by numerical simulation; (3) the off-axis behavior of double-droplet migration are found for the first time. The physical phenomena are analyzed and discussed from the mechanical point of view, and the coupling correlation of multiphysics fields and the influencing factors are analyzed.
The study of droplet thermocapillary migration due to the uniform temperature field is of great significance in revealing the role of thermocapillary force in the interfacial behavior. Liquid rocket combustion systems, space environment cooling systems, and the space welding are among the many applications of space technology that involve droplet thermocapillary migration. In future long-duration space flights, sewage treatment and food processing in space are issues that astronauts must consider, which involve problems of two-phase or multiphase separation. However, the droplet behavior cannot be manipulated more precisely in this paper, but it still lays theoretical and experimental foundations for the use of thermocapillary force to lead the interfacial behavior.