Next Article in Journal
Substituted Caffeic and Ferulic Acid Phenethyl Esters: Synthesis, Leukotrienes Biosynthesis Inhibition, and Cytotoxic Activity
Next Article in Special Issue
Solvent Supercritical Fluid Technologies to Extract Bioactive Compounds from Natural Sources: A Review
Previous Article in Journal
Protective Effect of Quercetin against Oxidative Stress-Induced Cytotoxicity in Rat Pheochromocytoma (PC-12) Cells
Previous Article in Special Issue
Pattern Evolution during Double Liquid-Vapor Phase Transitions under Weightlessness
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Measuring the Transition Rates of Coalescence Events during Double Phase Separation in Microgravity

1
Department of Physics and Astronomy, College of Charleston, Charleston, SC 29424, USA
2
CNRS, Institut de Chimie de la Matière Condensée de Bordeaux, ESEME, Université de Bordeaux, UPR 9048, F-33600 Pessac, France
3
Physique et Mécanique des Milieux Hétérogènes, UMR 7636 CNRS-ESPCI-Université Pierre et Marie Curie - Université Paris Diderot, 10 rue Vauquelin, 75005 Paris, France
*
Author to whom correspondence should be addressed.
Molecules 2017, 22(7), 1125; https://doi.org/10.3390/molecules22071125
Submission received: 15 June 2017 / Accepted: 30 June 2017 / Published: 6 July 2017
(This article belongs to the Special Issue Sub- and Supercritical Fluids and Green Chemistry)

Abstract

:
Phase transition is a ubiquitous phenomenon in nature, science and technology. In general, the phase separation from a homogeneous phase depends on the depth of the temperature quench into the two-phase region. Earth’s gravity masks the details of phase separation phenomena, which is why experiments were performed under weightlessness. Under such conditions, the pure fluid sulphur hexafluoride (SF 6 ) near its critical point also benefits from the universality of phase separation behavior and critical slowing down of dynamics. Initially, the fluid was slightly below its critical temperature with the liquid matrix separated from the vapor phase. A 0.2 mK temperature quench further cooled down the fluid and produced a double phase separation with liquid droplets inside the vapor phase and vapor bubbles inside the liquid matrix, respectively. The liquid droplets and the vapor bubbles respective distributions were well fitted by a lognormal function. The evolution of discrete bins of different radii allowed the derivation of the transition rates for coalescence processes. Based on the largest transition rates, two main coalescence mechanisms were identified: (1) asymmetric coalescences between one small droplet of about 20 μ m and a wide range of larger droplets; and (2) symmetric coalescences between droplets of large and similar radii. Both mechanisms lead to a continuous decline of the fraction of small radii droplets and an increase in the fraction of the large radii droplets. Similar coalescence mechanisms were observed for vapor bubbles. However, the mean radii of liquid droplets exhibits a t 1 / 3 evolution, whereas the mean radii of the vapor bubbles exhibit a t 1 / 2 evolution.

1. Introduction

Phase transition is a ubiquitous phenomenon in nature (e.g., the water cycle), science and technology [1]. In general, the phase separation from a homogeneous phase depends on the depth of the temperature quench into the two-phase region. The phase separation could occur through spinodal decomposition [2,3,4,5,6,7,8,9,10] or nucleation and growth [3,10,11,12,13,14,15,16]. These two processes can result in quite different morphologies. Usually, during the initial stage of the spinodal decomposition, we observe bicontinuous structures, whereas, during the nucleation and growth, we observe isolated domains in another continuous phase. The above classical theories of phase separation have some limitations. For example, the classical nucleation theory describes quite well the condensation of supercooled vapors, but gives only qualitative results when applied to supersaturated fluids. As a result, the nucleation and growth theory was further augmented—for example, with the density functional theory [17,18,19] and the diffuse interface theory [20]. Spinodal decomposition or nucleation path also depends on the volume fraction of the minority phase [9,21]. Additionally, the presence of solid walls and the wetting effects dramatically modify phase separation dynamics [22,23,24,25,26,27].
From an experimental perspective, the method most frequently used for producing phase separation is based on temperature quench in which a sample in a homogeneous state is usually heated up or cooled down. However, due to the short time scale of the phase separation and the finite quench speed, it is very difficult to study phase separation kinetics. Additionally, Earth’s gravity hampers experimental observation of fast phase separation because the vapor bubbles rise up and liquid droplets fall down, which quickly creates a flat meniscus determined by the gravity-induced coalescence of bubbles or droplets [28,29]. The microgravity environment for these experiments suppress the gravitational coalescence and allow a reasonably long time of observation of the phase-separating liquid droplets. In addition, dealing with fluids near their critical point allows critical scaling universality to be used to generalize the results to all fluids, and critical slowing down of the phase separation dynamic allows for the observation of phenomena that would otherwise be undetected due to their very fast dynamics. The DECLIC (Dispositif pour l’Étude de la Croissance et des Liquides Critiques) flight model facility is a multi-user facility for investigating critical fluids behavior and directional solidification of transparent alloys. The compact design contains three inserts, of which we refer here only to the ALI (Alice Like Insert) dedicated to the study of sulphur hexafluoride (SF 6 ) as a near-ambient temperature critical fluid. The program covers a whole characterization of SF 6 , ranging from thermodynamic quantities measurements (thermal diffusivity, heat capacity and turbidity near the critical point) to boiling effects studies [30].
This study reports on the double phase separation process that takes place after a supercritical SF 6 fluid is brought from an initial state, which is already in a two-phase state slightly below the critical temperature, to a temperature 0.2 mK below the initial value thanks to a temperature quench. Although the one-step phase separation process has been studied by several authors in weightlessness [11,12,13,24,28,31,32,33,34,35,36,37,38,39,40], double quench has only been the object of very few studies [35,41].
As discussed in detail in [41], the 0.2 mK temperature quench produced further phase separation with (1) vapor bubbles inside the majority liquid phase (the latter wetting the walls of the sample cell unit); and (2) liquid droplets inside the initially separated vapor phase. The novelty of this study is the quantitative investigation of the liquid droplets and vapor bubbles distributions, respectively, observed during the double phase separation process. We derived for the first time, to our knowledge, the transition rates that describe the evolution of different radii bins from recorded images of phase-separating systems. Two dominant coalescence mechanisms were found: (1) an asymmetric coalescence between small liquid droplets and a relatively wide range of other droplet sizes; and (2) a symmetric coalescence between large droplets of similar sizes.

2. Experimental Setup

The DECLIC experimental setup was described extensively in previous publications [30,41,42]. DECLIC used Alice Like Instrumentation (ALI) to carry out phase separation experiments near the critical point of SF 6 at room temperature in the weightlessness. The optical pressurized cells (or the Direct Observation Cell—DOC) contained a heater device as a transparent resistive layer appropriate for light transmission observation. The DOC was filled with SF 6 at its vapor-liquid critical point T c = 318.737 K, (45.587 C), p = 3.73 MPa, and ρ = 742.6 kg m 3 . A sketch of the cell is shown in Figure 1a (see also [41,43,44] for more details). The cylindrical DOC had an inner diameter of D = 10.6 mm and inner thickness of e = 4.115 mm, and was prepared with a mean density ρ = ρ c + 2 % . Three small thermistors were located inside the fluid volume, so that three local temperatures are measured close to the vapor–liquid interface in the microgravity environment (see Figure 1b). The light source was a 633 nm He-Ne laser.
The optical system allowed interlaced recordings of both Wide Field Of View (WFOV) images, which covers a circular area with diameter of 10.6 mm, and Narrow Field Of View (NFOV) images, which only covers a 1 mm × 1 mm area located at the center of the DOC [41,42]. Starting from an initial state, which is slightly below the critical temperature where the system is already phase-separated as shown in Figure 1b, the DOC is further cooled down by a temperature quench of 0.2 mK (see Figure 1c) that allowed the supercritical fluid to cross into a double phase separation regime in which inside the vapor bubbles we observed liquid droplets and inside the liquid matrix there were vapor bubbles (see Figure 1d).
All data analyzed here belong to ALI sequence 7 for which the temperature quench started at t = 0 s, i.e., the image reference index of i m g r e f = 1,123,710,576 with Δ = 23 frames per second. For example, the image shown in Figure 1b has an index of i m g = 1,123,708,278, which means it was taken at the actual time t = ( i m g i m g r e f ) / Δ = (1,123,708,278 − 1,123,710,576)/23 = −99.9 s from the temperature quench, i.e., approximately 100 s prior to 0.2 mK temperature quench. The image index for the snapshot shown in Figure 1d was i m g = 1,123,782,720, which corresponds to about 3136.7 s after the thermal quench.

3. Narrow Field of View (NFOV) Image Processing

Macroscopic images, or WFOV images, such as those shown in Figure 1b,d, are interlaced with microscopic images, or NFOV images, taken with a microscope objective with a magnification of ×12 (see Figure 1d for microscope location). In the NFOV images, one pixel is equivalent to 0.977 μ m.
Initially, the microscope captured images from inside a large vapor bubble that phase-separated before we applied the 0.2 mK temperature quench (as seen in Figure 1d and labeled “vapor phase”, i.e., prior to applying the 0.2 mK temperature quench). Inside this large vapor bubble, we observed coalescences between liquid droplets. As the time passes, the vapor phase recedes and the large vapor bubble seen on the right side of Figure 1d becomes smaller and smaller. As a result, at later stages, the microscope captured images from the region marked “matrix liquid” phase in Figure 1d. In this region, we also noticed coalescences as seen in Figure 2, but they take place between the vapor bubbles embedded in the liquid majority phase.
The high resolution recording and the slowing down of coalescence processes due to microgravity conditions allowed us to observe binary collisions as shown in Figure 2. Due to microgravity slowing down of all critical processes, we can follow a coalescence between two vapor bubbles for more than seven minutes (see Figure 2). Although we only focused and magnified the coalescence in the upper left corner of Figure 2a1–d1, there are other binary coalescences visible (notice the one in the upper right corner of Figure 2).

4. Results

4.1. Liquid Droplets Dynamics from NFOV Images

Inside any large vapor bubbles formed during the initial phase separation process prior to the application of the 0.2 mK thermal quench, a dynamic process of continuous condensation of liquid droplets from supercritical phase takes place. As we noticed from the WFOV images, the large vapor bubbles are in continuous and slow motion due to internal processes, such as droplets nucleation and expansion of the wetting layer at cell boundaries, Brownian motion, g-jitter, and residual steady gravity due to the fact that the sample is not at the spacecraft center of mass. As a result of the macroscopic motion, some microscopic recordings caught the slow drift of the interface between a large vapor bubble (with liquid droplets inside) and the majority liquid phase (with vapor bubbles inside).
We fitted the NFOV liquid droplets distributions both to Gauss (Figure 3a) and lognormal (Figure 3b) functions to gauge possible trends over time. We found out that there is no significant difference in the goodness of fit between the two fitting functions. However, we favor the lognormal distribution (Figure 3b) due to the presence of the asymmetric, long tail in the experimental data (see also [41] for details). Some radii bins, e.g., around 60 μ m, are significantly lower than the theoretical curves shown in Figure 3. One reason could be the fact that the 1 mm × 1 mm microscope window simply did not capture correctly the 60 μ m radii bins at that particular time. We assumed that the distribution of droplets as recorded by the 1 mm × 1 mm microscope is representative for the entire bulk distribution. However, the homogeneity assumption may not be locally valid and fluctuations are possible, especially since the system is still evolving through coalescences. Another contributing factor is the random coalescence processes that may at some particular instances deplete the droplets distributions with some specific radii, e.g., 60 μ m, while generating larger droplets.
The most noticeable result regarding droplet radii distribution is that initially the center of the distribution shifts towards larger and larger values (see Figure 3). This is consistent with the Brownian coalescence mechanism and the radii increase as t 1 / 3 (see the continuous line inside the first grayed rectangle in Figure 4c, and, for fitting details, see [41]). After the initial fast increase in the average radii, the distributions flattened, which suggests that the process reached a steady state (see the solid red circles in the middle portion of Figure 4b,c). Initially, the microscope captured images of liquid droplets inside the large vapor bubble seen in Figure 1d and schematically represented by the leftmost white-bordered square in Figure 4a. The mean radius of the distribution of liquid droplets initially shifts towards larger radii (see the grayed rectangles marked “liquid droplets” in Figure 4b,c).
At later times, we notice that the mean radius remains almost constant. While the large vapor bubble recedes, a phase separation interface between the liquid matrix and the vapor bubble comes into focus as shown schematically in the middle white-bordered square of Figure 4a and the corresponding solid red circles in Figure 4b,c. In this range of durations, the distributions are skewed due to the boundary effects induced by the moving interface line. Eventually, the vapor bubble recedes such that the interface is no longer visible in NFOV images, as shown in the rightmost white-bordered square in Figure 4a, and we only record vapor bubbles inside the liquid matrix (see the grayed rectangle marked “vapor bubbles” in Figure 4b,c). In this region, we found that the mean radius of vapor bubbles increases due to coalescences with t 1 / 2 (see the continuous line inside the last grayed rectangle in Figure 4c and also [41] for detailed discussions of possible mechanisms).
Although the mean radius of the distributions may increase over time because of coalescences, a more detailed picture of the coalescence mechanism emerged from the evolution of different radii bins. For this purpose, we binned all NFOV distributions with 10 pixels (≈10 μ m). When we focus on the initial stage of coalescence in the NFOV, i.e., for times below 1800 s and the corresponding t 1 / 3 trend in Figure 4c, it seems that the relative fraction of the droplets of radii below 40 μ m follow a descendent trend. At the same time, the relative fraction of the droplets with radii larger than 50 μ m seem to have an ascendant path (see [41] for detailed a description). This evolution could be the result of coalescence events that deplete the distribution of small droplets in favor of generating larger droplets (see the solid squares inside the first grayed rectangle of Figure 4b,c).

4.2. Vapor Bubbles Dynamics from NFOV Images

Inside the liquid matrix, a dynamic process of phase separation of vapor bubbles from supercritical phase takes place. As in the previous subsection, we also fitted the NFOV vapor bubbles distributions both to Gauss and lognormal functions to gauge possible trends over time (not shown, but see Figure 3 for similar trends seen in liquid droplets case and also [41]). Both the Gaussian and lognormal fit of the vapor bubbles distributions show similar trends. While the large vapor bubble (with liquid droplets inside) receded from the field of view, we initially recorded vapor bubbles distributions in the presence of a moving phase separation interface for times below approximately 5000 s. As we discussed in the previous subsection, such distributions (marked with solid circles in Figure 4b,c) are skewed due to the finite size effects near the phase separation interface. For this reason, we did not consider them when analyzing the trend of the mean radius versus time (see the continuous lines in Figure 4c that show the t 1 / 3 power law for liquid droplets evolution and t 1 / 2 power law for vapor bubbles evolution, respectively. For fitting details see also [41]).
After the phase separation interface is out of the field of view, we recorded images from the bulk of the liquid matrix that contained vapor bubbles undergoing coalescences. The most noticeable result regarding vapor bubbles radii distribution is that the center of the distribution shifts towards larger and larger values (see the second grayed rectangle in Figure 4b,c). This is consistent with a faster than Brownian coalescence mechanism, possibly a directional motion of vapor bubbles due to a composition Marangoni force [45], possible g-jitter and remaining steady gravity due to the fact that the sample is not at the spacecraft center of mass [41], and the radii increase as t 1 / 2 (see the continuous line inside the second grayed rectangle in Figure 4c).
As in the preceding subsection, a more detailed picture of the coalescence mechanism emerges when we investigate the evolution of different radii bins. During the initial stage of coalescences in the NFOV (for times below 5000 s), it seems that the dynamic is similar to the one observed across the phase separation interface inside the large vapor bubble. During the ascending trend of the average radius (see the grayed rectangles marked “vapor bubbles” in Figure 4b,c) the fraction of the small droplets below 30 μ m seem to decay over time and slightly increase for intermediate 40–50 μ m bubbles. Larger bubbles with radii from 60 to 90 μ m seem to also decrease their relative fraction, whereas large droplets over 100 μ m increase their relative contribution to the distributions. This evolution could be the result of coalescence events and has different dynamics than previous measurements in liquid droplets (see the previous subsection).

5. Theoretical Modeling

From a mathematical perspective, if we assume that a new droplet can only occur due to a binary coalescence between two droplets of radii r i ( t ) and r j ( t ) then the resultant droplet has a radius r k ( t + 1 ) = ( r i ( t ) 3 + r j ( t ) 3 ) 1 / 3 as in Figure 5a. Since we are limited by the image resolution, it is natural to consider only radii with integer values in pixels. In the following, we consider that measured radii are binned with a 1-pixel (≈0.977 μ m) increment. Therefore, by a coalescence event between a droplet of radius r i ( t ) = 1 pixel and r j ( t ) = 1 pixel it results in a droplet of radius r k ( t + 1 ) = 1.26 pixels, which will be counted in the 1-pixel bin at the time step t + 1 (see Figure 5a). Such an event occurs with a transition rate f 11 . Similarly, the droplets in the bin of size 2 pixels at time step t + 1 could only be created by collisions between r i ( t ) = 1 pixel and r j ( t ) = 2 pixels (see Figure 5a). The transition rate for such a collision is f 12 . The bin r 3 ( t + 1 ) could be generated by binary collisions r i ( t ) = 1 pixel and r j ( t ) = 3 pixels, or r j ( t ) = 2 pixels and r j ( t ) = 3 pixels, or r i ( t ) = 2 pixels and r j ( t ) = 2 pixels (see Figure 5a). Each possible combination has it own collision transition rate, i.e., f 13 , f 23 , and f 22 , respectively. Similarly, the 4-pixel bin could be generated by binary collisions r i ( t ) = 1 pixel and r j ( t ) = 4 pixels, or r i ( t ) = 2 pixels and r j ( t ) = 4 pixels, r i ( t ) = 3 pixels and r j ( t ) = 4 pixels, or r i ( t ) = 3 pixels and r j ( t ) = 3 pixels, and so forth.
As a reminder, 1 pixel W F O V 12.04 μ m in WFOV images and 1 pixel N F O V 0.977 μ m in NFOV images. The two-dimensional state space of all possible binary events that generate a droplet of radius r k at iteration time t + 1 is determined by all possible pairs as shown in Figure 5. The iteration time step is equal to the image acquisition sampling rate. The contour lines in Figure 5 separate points in integer value (pixel) radii. For example, the first contour line in Figure 5 goes between 0.5 and 1.5 pixels and only includes one possible data point, i.e., r i ( t ) = 1 pixel and r j ( t ) = 1 pixel, which could generate a droplet of radius ( 1 3 + 1 3 ) 1 / 3 = 1.26 pixels that rounds to a radius of 1 pixel at time t + 1 .
Similarly, the 3-pixel bin (with radii between 2.5 and 3.5 pixels) marked with a “3” inside the shaded area of Figure 5b, can only result from the following integer radii combinations: r i ( t ) = 2 pixels and r j ( t ) = 2 pixels, which gives r k ( t + 1 ) = ( 2 3 + 2 3 ) 1 / 3 2.52 pixels, or r i ( t ) = 1 pixel and r j ( t ) = 3 pixels, which gives r k ( t + 1 ) = ( 1 3 + 3 3 ) 1 / 3 3.04 pixels, or r i ( t ) = 2 pixels and r j ( t ) = 3 pixels, which gives r k ( t + 1 ) = ( 2 3 + 3 3 ) 1 / 3 3.27 pixels. All of the above three possible integer radii coalescence events would be counted in the 3-pixel bin at time step t + 1 .
To summarize, the following recursive relationships describe the evolution of the number n i ( t ) of droplet with the radius r i ( t ) by accounting for all possible binary collisions that could generate such a droplet:
Δ n 1 ( t ) = f 11 n 1 ( t ) 2 , Δ n 2 ( t ) = f 12 n 1 ( t ) n 2 ( t ) , Δ n 3 ( t ) = f 13 n 1 ( t ) n 3 ( t ) + f 22 n 2 ( t ) 2 + f 23 n 2 ( t ) n 3 ( t ) , Δ n 4 ( t ) = f 14 n 1 ( t ) n 4 ( t ) + f 24 n 2 ( t ) n 4 ( t ) + f 34 n 3 ( t ) n 4 ( t ) + f 33 n 3 ( t ) 2 , Δ n 5 ( t ) = f 15 n 1 ( t ) n 5 ( t ) + f 25 n 2 ( t ) n 5 ( t ) + f 35 n 3 ( t ) n 5 ( t ) + f 45 n 4 ( t ) n 5 ( t ) + f 44 n 4 ( t ) 2 ,
where Δ n i ( t ) = n i ( t + 1 ) n i ( t ) is the finite difference increase in the relative frequency distribution of droplets of radius r i = i pixels during one time step, f i j represents the transition rates of a stable collision process between droplets of radii r i = i pixels and r j = j pixels with i j . To avoid cluttering the recursive formulas above, we only mention that n 6 ( t + 1 ) is determined by the set of transition rates f 16 , f 26 , f 36 , f 46 , f 45 , f 55 ; n 7 ( t + 1 ) is determined by the set f 17 , f 27 , f 37 , f 47 , f 56 ; n 8 ( t + 1 ) is determined by the set f 18 , f 28 , f 38 , f 48 , f 56 , f 57 , f 66 ; n 9 ( t + 1 ) is determined by the set f 19 , f 29 , f 39 , f 49 , f 59 , f 58 , f 68 , f 78 , f 77 ; etc.
Although for each possible cluster of size r k there are multiple possible binary collisions, e.g., a droplet of radius four pixels could be generated by four distinct binary collision events according to Equation (1), we could always select the largest (most probable) value of all f i j values that could generate a droplet of radius r k . By connecting the states with the most likely (largest) binary collisions transition rates for each cluster shown in Figure 5, we obtained the phase space path of the system, which depends on many parameters, such as temperature, depth of the thermal quench that led to phase separation, volume fraction, etc.
We notice from Equation (1) that the transition rates f 11 ( t ) = ( n 1 ( t + 1 ) n 1 ( t ) ) / n 1 2 ( t ) and f 12 ( t ) = ( n 2 ( t + 1 ) n 2 ( t ) ) / ( n 1 ( t ) n 2 ( t ) ) are first order recursions as they only involve current state time t and the immediate next state at time t + 1 . However, for larger radii r k 2 (in pixels), there is a multiplicity of possible ways of generating a droplet of radius r k ( t + 1 ) that requires higher order recursions. For example, n 3 ( t + 1 ) n 3 ( t ) = f 13 n 1 ( t ) n 3 ( t ) + f 22 n 2 ( t ) 2 + f 23 n 2 ( t ) n 3 ( t ) is determined by three transition rates f 13 ( t ) , f 22 ( t ) and f 23 ( t ) and, therefore, required three successive distributions as initial conditions to solve the recursion:
n 3 ( t + 1 ) n 3 ( t ) = f 13 ( t ) n 1 ( t ) n 3 ( t ) + f 22 ( t ) n 2 ( t ) 2 + f 23 ( t ) n 2 ( t ) n 3 ( t ) , n 3 ( t + 2 ) n 3 ( t + 1 ) = f 13 ( t + 1 ) n 1 ( t + 1 ) n 3 ( t + 1 ) + f 22 ( t + 1 ) n 2 ( t + 1 ) 2 + f 23 ( t + 1 ) n 2 ( t + 1 ) n 3 ( t + 1 ) , n 3 ( t + 3 ) n 3 ( t + 2 ) = f 13 ( t + 2 ) n 1 ( t + 2 ) n 3 ( t + 2 ) + f 22 ( t + 2 ) n 2 ( t + 2 ) 2 + f 23 ( t + 2 ) n 2 ( t + 2 ) n 3 ( t + 2 ) .
Although in principle the transitions rates f i j ( t ) in Equation (2) could change over time to reflect the changes in thermophysical conditions (temperature quenches, external pressure, etc.), we assumed that, for the purpose of estimating them from the radii distributions at successive images, they remain quasiconstant. A possible justification for such an assumption is the fact that all our images are acquired on the thermal plateau (see Figure 1c). Additionally, even if the transition rates f i j change over time, we assumed that they do not change dramatically over, for example, three successive images. As a result, the transition rates f i j for a 3-pixel bin are solutions of the linear system:
n 3 ( t + 1 ) n 3 ( t ) n 3 ( t + 2 ) n 3 ( t + 1 ) n 3 ( t + 3 ) n 3 ( t + 2 ) = n 1 ( t ) n 3 ( t ) n 2 ( t ) 2 n 2 ( t ) n 3 ( t ) n 1 ( t + 1 ) n 3 ( t + 1 ) n 2 ( t + 1 ) 2 n 2 ( t + 1 ) n 3 ( t + 1 ) n 1 ( t + 2 ) n 3 ( t + 2 ) n 2 ( t + 2 ) 2 n 2 ( t + 2 ) n 3 ( t + 2 ) f 13 f 22 f 23 .
As we notice from Equation (1), for larger droplets, the multiplicity associated with the possible binary collisions increases quickly, which requires longer recursions.
In the WFOV images, we identified droplets at the lower end of an image resolution limit of a 1-pixel radius only in about 10% of analyzed images, which gives quite a space matrix when solving Equation (1). As a result, if we were to use the measured n 1 ( t ) for predicting n j ( t + 1 ) , then the solutions of Equation (1) could not always be determined. Therefore, we dropped any contribution of n 1 ( t ) from all recursions shown in Equation (1).
In general, both in the WFOV and NFOV analysis, we identified two dominant tendencies based on the clustering of the transition rates. One is an asymmetric coalescences that favors coalescences between small droplets, usually, r i 20–30 μ m and a broad range of other radii r j (30,120) μ m (see Figure 6 and Figure 7). The other cluster favors almost symmetrical coalescences between droplets of similar radii such as r i 80 μ m and a broad range of other similar radii, usually, r j (60,90) μ m.
The landscape of the positive transition rates f i j for WFOV contains different combinations of radii r i and r j (see Figure 6) with the above-mentioned two dominant clusters: one that generates collisions between small radius droplets r i 20–30 μ m and a wide range of other radii r j (20,120) μ m; and another subset with a radius around r i 70 μ m that generated binary collisions with similarly large droplets with r j (80,100) μ m.
For the NFOV images, we investigated separately the collision mechanism for the liquid droplets inside the large vapor bubble (see also Figure 7a) from the collision mechanism for the vapor bubbles inside the liquid matrix (see Figure 7b). As radii distributions change over time, the collision landscape of the positive transition rates f i j for different combinations of radii r i and r j would also change over time. Therefore, we again averaged over specific frames and plotted the contours of equal transition rates on a logarithmic scale (see Figure 7).
For the NFOV images of liquid droplet collisions inside the large vapor bubble, we found that the maximum positive transition rates (marked with solid circles in Figure 7a) have a distinctive pattern that favors large transition rates f i j for binary collisions between small droplets r i 20 μ m and a wide range of large droplets with r j (30,90) μ m. However, as we notice from Figure 7a, there are two dominant collision processes: (1) between droplets of almost equal size of 20 μ m and 30 μ m; and (2) between droplets of 30 μ m and 80 μ m. Both collision mechanisms suggest that the droplets of small radii will rapidly disappear in favor of larger droplets. The second group of relatively highly probable collisions is around r i 60–70 μ m that produces coalescences with similarly large droplets with r j 70–80 μ m (see Figure 7a). The second highly probable mechanism leads to a fast growth of the fraction of large droplets (as seen in Figure 7b).
For the NFOV images of vapor bubbles inside the liquid matrix, we found that the largest transition rates (marked with filled circles in Figure 7b) show a significantly different behavior, i.e., the binary collisions are primarily driven by asymmetric size bubbles such as 30 μ m and 80 μ m (see Figure 7b). At the same time, similar size collisions, however, are more evenly distributed across all sizes compared to liquid droplets case. For example, binary collisions between vapor bubbles of 30 μ m and 40 μ m have almost the same transition probability as 60 μ m and 70 μ m collisions.

6. Conclusions

There are two major novelties of this study: (1) we recorded during the same thermal plateau both the evolution of liquid droplets statistics immersed in a vapor bubble and the evolution of vapor bubbles immersed in the liquid matrix; and (2) based on the respective statistics we computed the transition rates that govern the evolution of individual radii. Both the WFOV and NFOV image analysis showed that the droplet distributions could be fitted with a lognormal function and the peak of the distributions shifts towards larger droplet radii over time. By following the temporal evolution of narrow bins of droplet radii (1-pixel increment in WFOV and 10 pixels increment in NFOV), we noticed that small radii distributions reduce their contribution over time, whereas the large radii distributions contribution increase. Since we only observed binary coalescences in all images, we derived a mathematical model that describes that rate of change of droplet radii distributions. By recursively solving for the transition rates, we identified two characteristic coalescence mechanisms: (1) dominantly asymmetric coalescences between one droplet of small radius (2 pixels in WFOV) and a broad range of other radii (between 2–10 pixels in WFOV); and (2) dominantly symmetric coalescences between one droplet of large radius (eight pixels in WFOV) and similar radii (between 6–8 pixels in WFOV).
Our novel approach to extracting the transition rates from images of liquid droplets and vapor bubbles, respectively, would allow us in the future to derive a realistic phenomenological model of droplet distributions, their evolution, and their dependence on different thermophysical parameters. It is known that, for large droplet diameters and binary coalescences, the rate of change of the number of domains was estimated as d N / d t = N 2 Σ p ( r ) V × n d Σ , where N is the number of domains per unit volume, p ( r ) is the pair distribution function of the tubes/droplets, V is the relative velocity of the tubes/droplets, and n is the outward normal to the collision surface Σ [31,46]. In these references, it is found that the long tail of the droplet size distribution, such as those found in both our WFOV and NFOV images can be reasonably modeled mathematically for radii over a natural cutoff with a power law, i.e., d N ( r ) / d r = A r θ [31,47]. Molecular dynamic simulations showed that, for a large volume fractions of the minority phase, the distribution of droplets versus their corresponding diameters becomes wider over time (see Figure 7b in [47]). Our study opens the possibility of modeling droplet dynamics not only in the limit of large radii, but for all radii, and could offer a more complete quantitative view of coalescence processes during phase separation.

Acknowledgments

A.O. acknowledges a research and development grant from the College of Charleston. Y.G., C.L. and D.B. acknowledge a research grant from Centre National d’Études Spatiales (CNES), and wish to thank the CNES and NASA teams involved in the DECLIC project. We are grateful to the anonymous reviewers for their constructive and pertinent comments.

Author Contributions

A.O. conducted image analysis and wrote the manuscript. Y.G., C.L. and D.B. are the Principal Investigators for the joint NASA/CNES project. They designed the experiment and collected the experimental data. They also reviewed the manuscript.

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses or interpretation of data; in the writing of the manuscript; nor in the decision to publish the results.

References

  1. Akgerman, A. Supercritical Fluids in Environmental Remediation and Pollution Prevention. In Supercritical Fluids; Abraham, M.A., Sunol, A.K., Eds.; American Chemical Society: Washington, DC, USA, 1997; Chapter 15; pp. 208–231. [Google Scholar]
  2. Bates, F.S.; Wiltzius, P. Spinodal decomposition of a symmetric critical mixture of deuterated and protonated polymer. J. Chem. Phys. 1989, 91, 3258–3274. [Google Scholar] [CrossRef]
  3. Beysens, D. Critical Phenomena. In Materials Sciences in Space, a Contribution to the Scientific Basis of Space Processing; Feuerbacher, B., Hamacher, H., Naumann, R., Eds.; Springer: Berlin, Germany, 1986; p. 191. [Google Scholar]
  4. Binder, K. Spinodal Decomposition. In Materials Science and Technology. A Comprehensive Treatment; Cahn, R.W., Haasen, P., Kramer, E., Eds.; Phase Transformations in Materials; Wiley-VCH: Weinheim, Germany, 1991; Chapter 7; Volume 5, pp. 405–471. [Google Scholar]
  5. Frisch, H.; Lebowitz, J.L. The Equilibrium Theory of Classical Fluids; Benjamin: New York, NY, USA, 1964. [Google Scholar]
  6. Zappoli, B.; Daniel, B.; Garrabos, Y. Heat Transfers and Related Effects in Supercritical Fluids; Springer: Dordrecht, The Netherlands, 2015; Volume 108. [Google Scholar]
  7. Abyzov, A.S.; Schmelzer, J.W.P. Nucleation versus spinodal decomposition in confined binary solutions. J. Chem. Phys. 2007, 127, 114504. [Google Scholar] [CrossRef] [PubMed]
  8. Cumming, A.; Wiltzius, P.; Bates, F.S.; Rosedale, J.H. Light-scattering experiments on phase-separation dynamics in binary fluid mixtures. Phys. Rev. A 1992, 45, 885. [Google Scholar] [CrossRef] [PubMed]
  9. Perrot, F.; Guenoun, P.; Baumberger, T.; Beysens, D.; Garrabos, Y.; Le Neindre, B. Nucleation and growth of tightly packed droplets in fluids. Phys. Rev. Lett. 1994, 73, 688–691. [Google Scholar] [CrossRef] [PubMed]
  10. Schmelzer, J.W.P.; Abyzov, A.S.; Moller, J. Nucleation versus spinodal decomposition in phase formation processes in multicomponent solutions. J. Chem. Phys. 2004, 121, 6900–6917. [Google Scholar] [CrossRef] [PubMed]
  11. Beysens, D.; Straub, J.; Turner, D. Phase Transitions and Near-Critical Phenomena. In Space, a European Perspective; Walter, H., Ed.; Springer: Berlin, Germany, 1987; p. 221. [Google Scholar]
  12. Domb, C.; Lebowitz, J.L.; Lebowitz, J. Phase Transitions and Critical Phenomena; Academic Press: London, UK, 2001; p. 320. [Google Scholar]
  13. Onuki, A. Phase Transition Dynamics; Cambridge University Press: Cambridge, UK, 2002. [Google Scholar]
  14. Schmelzer, J.W.P. Comments on the Nucleation Theorem. J. Colloid Interface Sci. 2001, 242, 354–372. [Google Scholar] [CrossRef]
  15. Schmelzer, J.W.P.; Schmelzer, J.J.; Gutzow, I.S. Reconciling Gibbs and van der Waals: A new approach to nucleation theory. J. Chem. Phys. 2000, 112, 3820. [Google Scholar] [CrossRef]
  16. Schmelzer, J.W.P.; Schmelzer, J.J. Kinetics of condensation of gases: A new approach. J. Chem. Phys. 2001, 114, 5180. [Google Scholar] [CrossRef]
  17. Cahn, W.J.; Hilliard, E.J. Free Energy of a Nonuniform System. I. Interfacial Free Energy. J. Chem. Phys. 1958, 29, 258–267. [Google Scholar] [CrossRef]
  18. Cahn, J.; Hilliard, J. Free energy of a non-uniform system III: Nucleation in a two-component incompressible fluid. J. Chem. Phys. 1959, 31, 688–699. [Google Scholar] [CrossRef]
  19. Cahn, J. On spinodal decomposition. Acta Metall. 1961, 9, 795–801. [Google Scholar] [CrossRef]
  20. Granasy, L. Nucleation and spinodal decomposition. Solid State Phenom. 1997, 56, 67. [Google Scholar] [CrossRef]
  21. Oprisan, A.; Oprisan, S.A.; Hegseth, J.J.; Garrabos, Y.; Lecoutre-Chabot, C.; Beysens, D. Universality in early-stage growth of phase-separating domains near the critical point. Phys. Rev. E 2008, 77, 051118. [Google Scholar] [CrossRef] [PubMed]
  22. Bonn, D.; Bertrand, E.; Meunier, J.; Blossey, R. Dynamics of Wetting Layer Formation. Phys. Rev. Lett. 2000, 84, 4661–4664. [Google Scholar] [CrossRef] [PubMed]
  23. Fenistein, D.; Bonn, D.; Rafai, S.; Wegdam, G.; Meunier, J.; Parry, A.; TelodaGama, M. What controls the thickness of wetting layers near bulk criticality? Phys. Rev. Lett. 2002, 89, 096101. [Google Scholar] [CrossRef] [PubMed]
  24. Guenoun, P.; Beysens, D.; Robert, M. Dynamics of wetting and phase separation. Phys. Rev. Lett. 1990, 65, 2406–2409. [Google Scholar] [CrossRef] [PubMed]
  25. Hegseth, J.; Garrabos, Y.; Nikolayev, V.S.; Lecoutre-Chabot, C.; Wunenburger, R.; Beysens, D. Gas wets a solid wall in orbit. Int. J. Thermophys. 2002, 23, 89–101. [Google Scholar] [CrossRef]
  26. Lipowsky, R.; Huse, D.A. Diffusion-Limited Growth of Wetting Layers. Phys. Rev. Lett. 1986, 57, 353. [Google Scholar] [CrossRef] [PubMed]
  27. Steiner, U.; Klein, J. Growth of Wetting Layers from Liquid Mixtures. Phys. Rev. Lett. 1996, 77, 2526–2529. [Google Scholar] [CrossRef] [PubMed]
  28. Beysens, D.; Garrabos, Y. The phase transition of gas and liquids. Physica A 2000, 281, 361–380. [Google Scholar] [CrossRef]
  29. Beysens, D.; Forgacs, G.; Glazier, J.A. Cell sorting is analogous to phase ordering in fluids. Proc. Natl. Acad. Sci. USA 2000, 97, 9467–9471. [Google Scholar] [CrossRef] [PubMed]
  30. Pont, G.; Barde, S.; Blonde, D.; Zappoli, B.; Garrabos, Y.; Lecoutre, C.; Beysens, D.; Hicks, M.; Hegde, U.; Hahn, I.; et al. DECLIC, soon two years of successful operations. In Proceedings of the 62nd International Astronautical Congress, Cape Town, South Africa, 3–7 October 2011; International Astronautical Federation: Paris, France, 2011; pp. 1–12. [Google Scholar]
  31. Beysens, D.; Guenoun, P.; Sibille, P.; Kumar, A. Dimple and nose coalescences in phase-separation processes. Phys. Rev. E 1994, 50, 1299–1302. [Google Scholar] [CrossRef]
  32. Beysens, D. Kinetics and morphology of phase separation in fluids: The role of droplet coalescence. Physica A 1997, 239, 329–339. [Google Scholar] [CrossRef]
  33. Binder, K. Theory for the dynamics of “clusters.” II. Critical diffusion in binary systems and the kinetics of phase separation. Phys. Rev. B 1977, 15, 4425–4447. [Google Scholar] [CrossRef]
  34. Hegseth, J.; Nikolayev, V.; Beysens, D.; Garrabos, Y.; Chabot, C. Growth and Morphology of Phase Separating Supercritical Fluids (GMSF), Boiling in Subcritical Fluids, and Critical Fluctuations. In Proceedings of the Fourth Microgravity Fluid Physics and Transport Phenomena Conference, Cleveland, OH, USA, 12–14 August 1998; National Center for Microgravity Research on Fluids and Combustion: Cleveland, OH, USA, 1998; Volume NASA/CP-1999-208526/SUPPL1, pp. 184–189. [Google Scholar]
  35. Klein, H.; Schmitz, G.; Woermann, D. Temperature propagation in near-critical fluids prior to and during phase separation. Phys. Rev. A 1991, 43, 4562–4563. [Google Scholar] [CrossRef] [PubMed]
  36. Oprisan, A. Fluctuations, Phase Separation and Wetting Films near Liquid-Gas Critical Point. Ph.D. Thesis, University of New Orleans, New Orleans, LA, USA, 2006. [Google Scholar]
  37. Perrot, F.; Beysens, D.; Garrabos, Y.; Froumlhlich, T.; Guenoun, P.; Bonetti, M.; Bravais, P. Morphology transition observed in a phase separating fluid. Phys. Rev. E 1999, 59, 3079. [Google Scholar] [CrossRef]
  38. Roy, S.; Das, S.K. Effects of domain morphology on kinetics of fluid phase separation. J. Chem. Phys. 2013, 139, 044911. [Google Scholar] [CrossRef] [PubMed]
  39. Sagui, C.; O’Gorman, D.S.; Grant, M. Nucleation, growth and coarsening in phase-separating systems. Scanning Microsc. 1998, 12, 3–8. [Google Scholar]
  40. Tanaka, H. New mechanisms of droplet coarsening in phase-separating fluid mixtures. J. Chem. Phys. 1997, 107, 3734. [Google Scholar] [CrossRef]
  41. Oprisan, A.; Garrabos, Y.; Lecoutre, C.; Beysens, D. Pattern Evolution during Double Liquid-Vapor Phase Transitions under Weightlessness. Molecules 2017, 22, 947. [Google Scholar] [CrossRef] [PubMed]
  42. Marcout, R.; Raymond, G.; Martin, B.C.G.; Zappoli, B.; Duclos, F.; Barde, S.; Beysens, D.; Garrabos, Y.; Lecoutre, C.; Billia, B.; et al. DECLIC: A facility to investigate fluids and transparent materials in microgravity conditions in ISS. In Proceedings of the 57th International Astronautical Congress, Valencia, Spain, 2–6 October 2006; American Institute of Aeronautics and Astronautics: Reston, VA, USA, 2006. [Google Scholar]
  43. Beysens, D.; Pichavant, G.; Chatain, D.; Nikolayev, V.; Lecoutre, C.; Garrabos, Y. Non Marangoni motion of bubble under temperature gradient. In Proceedings of the 62nd International Astronautical Congress, Cape Town, South Africa, 3–7 October 2011; International Astronautical Federation: Paris, France, 2011. [Google Scholar]
  44. Beysens, D.; Chatain, D.; Nikolayev, V.S.; Ouazzani, J.; Garrabos, Y. Possibility of long-distance heat transport in weightlessness using supercritical fluids. Phys. Rev. E 2010, 82, 061126. [Google Scholar] [CrossRef] [PubMed]
  45. Shimizu, R.; Tanaka, H. A novel coarsening mechanism of droplets in immiscible fluid mixtures. Nat. Commun. 2015, 6, 7407. [Google Scholar] [CrossRef] [PubMed]
  46. Oprisan, A.; Hegseth, J.; Smith, G.; Lecoutre, C.; Garrabos, Y.; Beysens, D.A. Dynamics of wetting layer and Marangoni convection in microgravity. Phys. Rev. E 2011, 84, 021202. [Google Scholar] [CrossRef] [PubMed]
  47. Martula, D.; Hasegawa, T.; Lloyd, D.; Bonnecaze, R. Coalescence-induced coalescence of inviscid droplets in a viscous fluid. J. Colloid Interface Sci. 2000, 232, 241–253. [Google Scholar] [CrossRef] [PubMed]
Sample Availability: Samples of the compounds are not available from the authors.
Figure 1. (a) a schematic cross section (not to scale) of the Direct Observation Cell (DOC) dedicated to studying sulfur hexafluoride (SF 6 ) phase separation in Dispositif pour l’Étude de la Croissance et des Liquides Critiques (DECLIC) with Alice Like Instrumentation (ALI). experiments; (b) a wide field of view (WFOV) image of the DOC at equilibrium in the two-phase range recorded 100 s prior to the application of 0.2 mK temperature quench; (c) the DOC was slightly below the critical temperature when the cooling down temperature quench of 0.2 mK was applied at t = 0 s; (d) as a result of the 0.2 mK thermal quench, a double phase separation takes place with liquid droplets nucleating inside the vapor phase and vapor bubbles forming inside the liquid matrix. A microscope magnifies a 1 mm × 1 mm square area at the center of the DOC.
Figure 1. (a) a schematic cross section (not to scale) of the Direct Observation Cell (DOC) dedicated to studying sulfur hexafluoride (SF 6 ) phase separation in Dispositif pour l’Étude de la Croissance et des Liquides Critiques (DECLIC) with Alice Like Instrumentation (ALI). experiments; (b) a wide field of view (WFOV) image of the DOC at equilibrium in the two-phase range recorded 100 s prior to the application of 0.2 mK temperature quench; (c) the DOC was slightly below the critical temperature when the cooling down temperature quench of 0.2 mK was applied at t = 0 s; (d) as a result of the 0.2 mK thermal quench, a double phase separation takes place with liquid droplets nucleating inside the vapor phase and vapor bubbles forming inside the liquid matrix. A microscope magnifies a 1 mm × 1 mm square area at the center of the DOC.
Molecules 22 01125 g001
Figure 2. Binary collisions between vapor bubbles in narrow filed of view (NFOV) images at successive times (a) 6607.3 s; (b) 6879.7 s; (c) 6970.4 s; and (d) 7061.2 s. The highlighted white-bordered rectangle in panes (a1d1) show two vapor bubbles of slightly different radii that approached each other (a1,a2) and then form a continuous bubble (b1,b2), which slowly changes into an oblongated shape (c1,c2) and finally becomes spherical again (d1,d2). We also magnified the region of interest for better visualization of coalescence (see panels a2d2).
Figure 2. Binary collisions between vapor bubbles in narrow filed of view (NFOV) images at successive times (a) 6607.3 s; (b) 6879.7 s; (c) 6970.4 s; and (d) 7061.2 s. The highlighted white-bordered rectangle in panes (a1d1) show two vapor bubbles of slightly different radii that approached each other (a1,a2) and then form a continuous bubble (b1,b2), which slowly changes into an oblongated shape (c1,c2) and finally becomes spherical again (d1,d2). We also magnified the region of interest for better visualization of coalescence (see panels a2d2).
Molecules 22 01125 g002
Figure 3. Histograms of microscopic views at successive times 620.9 s (solid squares), 1437.9 s (solid circles), and 2704.6 s (solid triangles). The continuous lines show the corresponding Gauss (a) and lognormal (b) fits. For image at 620.9 s (solid squares) the center of Gauss distribution is at x c G = (28.7 ± 0.4) μ m and of lognormal is at x c L N = (29.5 ± 0.8) μ m. The corresponding standard deviations are w G = (8.4 ± 0.5) μ m and w L N = (0.28 ± 0.03) μ m. For the snapshot at 1437.9 s (solid circles), x c G = (36.1 ± 1.1) μ m and x c L N = (37.0 ± 1.8) μ m with the standard deviations and w G = (12.3 ± 2.4) μ m and w L N = (0.33 ± 0.06) μ m. For the image at 2704.6 s (solid triangles), x c G = (38.6 ± 1.0) μ m and x c L N = (39.9 ± 1.0) μ m with w G = (9.8 ± 1.5) μ m and w L N = (0.29 ± 0.04) μ m.
Figure 3. Histograms of microscopic views at successive times 620.9 s (solid squares), 1437.9 s (solid circles), and 2704.6 s (solid triangles). The continuous lines show the corresponding Gauss (a) and lognormal (b) fits. For image at 620.9 s (solid squares) the center of Gauss distribution is at x c G = (28.7 ± 0.4) μ m and of lognormal is at x c L N = (29.5 ± 0.8) μ m. The corresponding standard deviations are w G = (8.4 ± 0.5) μ m and w L N = (0.28 ± 0.03) μ m. For the snapshot at 1437.9 s (solid circles), x c G = (36.1 ± 1.1) μ m and x c L N = (37.0 ± 1.8) μ m with the standard deviations and w G = (12.3 ± 2.4) μ m and w L N = (0.33 ± 0.06) μ m. For the image at 2704.6 s (solid triangles), x c G = (38.6 ± 1.0) μ m and x c L N = (39.9 ± 1.0) μ m with w G = (9.8 ± 1.5) μ m and w L N = (0.29 ± 0.04) μ m.
Molecules 22 01125 g003
Figure 4. (a) a schematic representation of the sample cell unit (not at scale) containing one large vapor bubble inside the liquid matrix. The microscope view shifts over time and at early stages captures liquid droplets distributions from inside the large vapor bubble (see the leftmost white-bordered square). Later on, the microscope captures the vapor–liquid interface and a mixture of distributions (see the middle white-bordered square). At late stages, the microscope view only captures vapor bubbles embedded in the liquid matrix (see the rightmost white-bordered square). The location of the center of the Gaussian (b), respectively, lognormal (c) distributions shift over time. At first, the distribution of liquid droplets shifts towards larger radii (solid squares inside the grayed “liquid droplets” area). At some point, the interface line between the large vapor bubble and the liquid matrix drifts into focus between 2981.2 s and 5790.3 s (solid red circles) and significantly changes the trend of the distributions due to boundary effects. While the best fit is still a lognormal distribution, the center of the distribution seems to remain constant and even moves towards smaller radii. At a later time, after the vapor bubble completely receded from focus, the distribution of liquid droplets also shows an increase in its mean radius (see the grayed rectangle marked “liquid droplets”). The continuous lines in (c) indicate a t 1 / 3 evolution for liquid droplet radii and a t 1 / 2 evolution for vapor droplets (see [41] for data fitting details and a discussion regarding possible coalescence mechanisms).
Figure 4. (a) a schematic representation of the sample cell unit (not at scale) containing one large vapor bubble inside the liquid matrix. The microscope view shifts over time and at early stages captures liquid droplets distributions from inside the large vapor bubble (see the leftmost white-bordered square). Later on, the microscope captures the vapor–liquid interface and a mixture of distributions (see the middle white-bordered square). At late stages, the microscope view only captures vapor bubbles embedded in the liquid matrix (see the rightmost white-bordered square). The location of the center of the Gaussian (b), respectively, lognormal (c) distributions shift over time. At first, the distribution of liquid droplets shifts towards larger radii (solid squares inside the grayed “liquid droplets” area). At some point, the interface line between the large vapor bubble and the liquid matrix drifts into focus between 2981.2 s and 5790.3 s (solid red circles) and significantly changes the trend of the distributions due to boundary effects. While the best fit is still a lognormal distribution, the center of the distribution seems to remain constant and even moves towards smaller radii. At a later time, after the vapor bubble completely receded from focus, the distribution of liquid droplets also shows an increase in its mean radius (see the grayed rectangle marked “liquid droplets”). The continuous lines in (c) indicate a t 1 / 3 evolution for liquid droplet radii and a t 1 / 2 evolution for vapor droplets (see [41] for data fitting details and a discussion regarding possible coalescence mechanisms).
Molecules 22 01125 g004
Figure 5. (a) binary collisions between droplets of radii r i ( t ) and r j ( t ) lead to a larger droplet of radius r k ( t + 1 ) with a collision transition rate f i j . As the resultant droplet radius (binned in increments of one pixel) increased, the multiplicity of possible combinations increases. The number inside the circles symbolizing droplets represents their radii (in pixels) before coalescence; (b) the two-dimensional state space of the system at time step t shows all possible pairs of liquid droplet radii r i ( t ) and r j ( t ) < r i ( t ) (below continuous diagonal line) that could lead to a new droplet of radius r k at iteration t + 1 . Each shaded contour plot covers a 1-pixel bin. For example, the shaded contour label with three pixels shows that it is the result of three possible types of coalescence events between droplets of following sizes: two pixels and two pixels, or one pixel and three pixels, or two pixels and three pixels.
Figure 5. (a) binary collisions between droplets of radii r i ( t ) and r j ( t ) lead to a larger droplet of radius r k ( t + 1 ) with a collision transition rate f i j . As the resultant droplet radius (binned in increments of one pixel) increased, the multiplicity of possible combinations increases. The number inside the circles symbolizing droplets represents their radii (in pixels) before coalescence; (b) the two-dimensional state space of the system at time step t shows all possible pairs of liquid droplet radii r i ( t ) and r j ( t ) < r i ( t ) (below continuous diagonal line) that could lead to a new droplet of radius r k at iteration t + 1 . Each shaded contour plot covers a 1-pixel bin. For example, the shaded contour label with three pixels shows that it is the result of three possible types of coalescence events between droplets of following sizes: two pixels and two pixels, or one pixel and three pixels, or two pixels and three pixels.
Molecules 22 01125 g005
Figure 6. A logarithmic contour plot of the positive transition rates for WHOV images show solid circles proportional to the relative power law exponent (see the color legend). The logarithmic contour plot shows two dominant mechanisms for generating droplets: (1) small droplets that collide with droplets of any radius; and (2) preferential coalescences of large droplets with similar size.
Figure 6. A logarithmic contour plot of the positive transition rates for WHOV images show solid circles proportional to the relative power law exponent (see the color legend). The logarithmic contour plot shows two dominant mechanisms for generating droplets: (1) small droplets that collide with droplets of any radius; and (2) preferential coalescences of large droplets with similar size.
Molecules 22 01125 g006
Figure 7. (a) the NFOV images of liquid droplets collisions inside the large vapor bubble. The landscape of the positive transition rates is dominated by two mechanisms: (1) low and similar radii (20 μ m with 30 μ m); and (2) asymmetric sizes at 30 μ m and 80 μ m; (b) NFOV images of vapor bubbles inside the liquid matrix. The transition rates for binary collisions are dominated by 30 μ m and 80 μ m sizes. At the same time, similar size collisions are much more likely than in the case of liquid droplets. All contour plots are shown on a logarithmic scale. The size of the solid circle is proportional to the transition probability exponent.
Figure 7. (a) the NFOV images of liquid droplets collisions inside the large vapor bubble. The landscape of the positive transition rates is dominated by two mechanisms: (1) low and similar radii (20 μ m with 30 μ m); and (2) asymmetric sizes at 30 μ m and 80 μ m; (b) NFOV images of vapor bubbles inside the liquid matrix. The transition rates for binary collisions are dominated by 30 μ m and 80 μ m sizes. At the same time, similar size collisions are much more likely than in the case of liquid droplets. All contour plots are shown on a logarithmic scale. The size of the solid circle is proportional to the transition probability exponent.
Molecules 22 01125 g007

Share and Cite

MDPI and ACS Style

Oprisan, A.; Garrabos, Y.; Lecoutre, C.; Beysens, D. Measuring the Transition Rates of Coalescence Events during Double Phase Separation in Microgravity. Molecules 2017, 22, 1125. https://doi.org/10.3390/molecules22071125

AMA Style

Oprisan A, Garrabos Y, Lecoutre C, Beysens D. Measuring the Transition Rates of Coalescence Events during Double Phase Separation in Microgravity. Molecules. 2017; 22(7):1125. https://doi.org/10.3390/molecules22071125

Chicago/Turabian Style

Oprisan, Ana, Yves Garrabos, Carole Lecoutre, and Daniel Beysens. 2017. "Measuring the Transition Rates of Coalescence Events during Double Phase Separation in Microgravity" Molecules 22, no. 7: 1125. https://doi.org/10.3390/molecules22071125

APA Style

Oprisan, A., Garrabos, Y., Lecoutre, C., & Beysens, D. (2017). Measuring the Transition Rates of Coalescence Events during Double Phase Separation in Microgravity. Molecules, 22(7), 1125. https://doi.org/10.3390/molecules22071125

Article Metrics

Back to TopTop