Next Article in Journal
A Multi-Part Orientation Planning Schema for Fabrication of Non-Related Components Using Additive Manufacturing
Previous Article in Journal
Zinc Carboxylate Surface Passivation for Enhanced Optical Properties of In(Zn)P Colloidal Quantum Dots
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design of Improved Flow-Focusing Microchannel with Constricted Continuous Phase Inlet and Study of Fluid Flow Characteristics

1
Key Laboratory of Metallurgical Equipment and Control Technology of Ministry of Education, Wuhan University of Science and Technology, Wuhan 430081, China
2
The State Key Laboratory of Fluid Power and Mechatronic Systems, Zhejiang University, Hangzhou 310027, China
3
College of Mechanical Engineering, Hubei University of Automotive Technology, Shiyan 442002, China
*
Author to whom correspondence should be addressed.
Micromachines 2022, 13(10), 1776; https://doi.org/10.3390/mi13101776
Submission received: 9 September 2022 / Revised: 7 October 2022 / Accepted: 17 October 2022 / Published: 19 October 2022

Abstract

:
This paper proposed an improved flow-focusing microchannel with a constricted continuous phase inlet to increase microbubble generation frequency and reduce microbubbles’ diameter. The design variables were obtained by Latin hypercube sampling, and the radial basis function (RBF) surrogate model was used to establish the relationship between the objective function (microbubble diameter and generation frequency) and the design variables. Moreover, the optimized design of the nondominated sorting genetic algorithm II (NSGA-II) algorithm was carried out. Finally, the optimization results were verified by numerical simulations and compared with those of traditional microchannels. The results showed that dripping and squeezing regimes existed in the two microchannels. The constricted continuous phase inlet enhanced the flow-focusing effect of the improved microchannel. The diameter of microbubbles obtained from the improved microchannel was reduced from 2.8141 to 1.6949 μm, and the generation frequency was increased from 64.077 to 175.438 kHz at the same capillary numbers (Ca) compared with the traditional microchannel. According to the fitted linear function, it is known that the slope of decreasing microbubble diameter with increasing Ca number and the slope of increasing generation frequency with increasing Ca number are greater in the improved microchannel compared with those in the traditional microchannel.

1. Introduction

Nowadays, microbubbles have become one of the most effective contrast agents in medical ultrasound imaging or as a carrier for targeted drug delivery [1,2,3,4,5,6,7]. Generally, the required microbubble diameter must be between 1 and 7 μm [8]. Currently, there are many ways to prepare microbubbles. Due to its high throughput and low cost, ultrasound is commonly used to generate microbubbles. Still, the ultrasound generator’s frequency, power, and pulse affect the size of the microbubbles and cannot produce uniform microbubbles of a specific size, which requires an additional filtering process to remove larger bubbles [9]. The mechanical stirring method [10] prepares microbubbles by mechanical stirrers but produces microbubbles of uneven size. The inkjet printing method [11] allows better control of the size of microbubbles but is limited to the production of liquid-filled particles. The coaxial electrofluidic atomization method [12] produces microbubbles with diameters less than 10 μm, but it also produces microbubbles of uneven size. All of the above methods produce polydisperse microbubbles in liquids, limiting the potential use of microbubbles in medicine. In contrast, microfluidics can produce highly monodisperse microbubbles [13,14].
Microfluidics has been used as a promising interdisciplinary technique for generating and manipulating monodisperse and discrete micron-sized bubbles. It is possible to reduce the size of the generated microbubbles by adjusting the flow parameters and narrowing the orifice structure to produce microbubbles with very narrow dimensions. It has the advantages of monodisperse bubble generation, easy control of bubble size, small sample volume, and low cost. Standard microfluidic techniques for generating microbubbles can be divided into three main types: the T-type, coaxial flow type, and flow-focusing type. The T-joint device has a simple structure, including liquid and gas channels that vertically intersect with the liquid channel [15]. Microbubble formation in the T-joint can be divided into three regimes: squeezing, dripping, and transiting [16]. The generation of microbubbles less than 10 μm in diameter for typical T-type structures is still problematic [17]. Coaxial flow devices usually have an exit orifice wider than the microbubble diameter. The liquid velocity is always high enough, so liquid inertia is the dominant force for microbubble formation. The coaxial flow device is characterized by a thin capillary tube supplying the gas surrounded by a coaxially arranged capillary tube supplying the liquid [18]. Flow-focusing microchannels usually consist of a central gas channel, two lateral liquid channels, and a narrow exit hole [19]. Flow-focusing microchannels use a constriction hole downstream of the dispersed phase’s inlet to modify the flow field’s hydrodynamic behavior and, thus, the formation of microbubbles. Compared with other microchannels, the microbubble generation in flow-focusing microchannels is easier to control. The bubble generation is relatively stable, efficient, and well-homogenized and has a more extensive adjustable size range. Therefore, it has attracted a lot of attention from researchers. Some experimental and numerical studies have been conducted.
Numerous visualization experiments have been conducted to investigate how microbubbles are produced in microchannels. Dietrich et al. [20] used flow-focusing microdevices to generate bubbles and measure the flow field with the micro-PIV technique. They found that the formation of bubbles depends on the geometry of the intersection and that the shear force was higher at 180° than at 60°, which indicates that the shear stress is more pronounced when the angle increases. In addition, the bubbles were more elongated in the 60° geometry than in the 180° geometry, suggesting an increase in surface tension due to deformation as the angle of connection decreases. Fu et al. [21] experimentally observed that tiny double bubbles could be formed in flow-focusing microchannels from a highly deformed gas conical jet at high liquid flow rates, whereas the gaseous jet in the downstream channel was stable. Castro-Hernández et al. [22] proposed a new mode of operation for the flow-focusing microfluidic devices, producing monodisperse bubbles with a diameter of 5 µm at a production rate of more than 105 Hz by properly adjusting the flow rates of the gas and liquid. Sheng et al. [23] first reported the generation of short bubbles at high frequency in capillary embedded step T-junction microdevices, and the mechanism of short bubble formation is mainly because of the relatively high pressure drop inside the step T-junction, which provides a much higher breakup force for the squeezing flow. The above experiments provide a visual view of microbubble generation and limited insight into the hydrodynamic mechanism of microbubble generation.
In addition to experimental methods, numerical simulations based on computational fluid dynamics methods provide detailed fluid dynamics information and become an effective tool for studying multiphase fluid flows. Jia and Zhang [24] conducted a three-dimensional numerical study of the bubble formation process in viscous fluids and found that the geo-reconstruction method performs better. Zadeh and Radespiel et al. [25] investigated the effects of two-phase flow rate, two-phase viscosity, surface tension, and geometrical channel characteristics (channel width and injection angle) on microbubble formation and its length based on the volume of fluid (VOF) method. The results showed that at a constant flow rate in the continuous phase, the bubble length significantly increased with an increasing flow rate in the dispersed phase; the bubble length increased as the injection angle increased from 45° to 90°. Chekifi et al. [26] used the coupled level set and volume of fluid (CLSVOF) method to track the interface and showed that the velocity ratio, interfacial tension, outer channel diameter, continuous phase viscosity, orifice width, and length have important effects on the bubble size and shape. At lower capillary numbers, an increase in the flow velocity ratio produces smaller bubble sizes in a shorter time, whereas an increase in the interfacial tension produces larger bubbles.
In summary, most researchers have studied the two-phase flow from its operational parameters and physical properties. Most early studies only considered the role of individual geometrical parameters without integrating the effect of geometrical parameters on microbubble generation. There was a lack of studies on the optimal design of flow-focusing microchannels. Geometry in flow-focusing devices plays an essential role in microbubbles’ size and generation frequency. In addition, these effects are necessary for the optimal design of microbubble generation devices. Therefore, it is crucial to study the geometry of flow-focusing microchannels.
This paper proposed an improved flow-focusing microchannel with constricted continuous phase entrance. A two-dimensional numerical model of microbubble generation in a flow-focusing microchannel was established using the VOF method. Four design variables of the microchannel were selected, and the radial basis function (RBF) surrogate model was used to establish the relationship between the objective function (microbubble diameter and generation frequency) and design variables. Based on this, the NSGA-II algorithm’s optimal design was carried out. Finally, the optimization results were verified by numerical simulations, and smaller microbubbles with higher generation frequency were obtained in comparison with the traditional microchannels. The fluid flow characteristics of microchannels were also further studied to investigate the effect of microchannels on microbubble generation.

2. Model Description

2.1. Structural Design of Flow-Focusing Microchannels

The flow-focusing microchannel is a typical axisymmetric microchannel, consisting of a dispersed phase inlet and two continuous phase inlets to form focusing holes. The microbubble generation is achieved by passing the dispersed and continuous phases through the focusing holes. The geometry of the flow-focusing microchannel is shown in Figure 1. The geometry is fully defined by six parameters, namely the orifice width ( W d ), continuous phase inlet constriction width ( W c 1 ), continuous phase inlet width ( W c 2 ), outlet width ( W o u t ), constriction angle ( α ), and expansion angle ( β ), where the orifice width is 2 μm and the continuous phase inlet width is 2 W d = 4 μm. To facilitate the uniform adjustment of the structural parameters, the dispersed phase inlet width is kept the same as the orifice width, and the continuous phase inlet constriction width and outlet width are normalized to the orifice width as shown in Equations (1) and (2).
W ¯ c 1 = W c 1 W d
W ¯ o u t = W o u t W d
where W ¯ c 1 is the normalized continuous phase inlet constriction width, and W ¯ o u t is the normalized outlet width. The gas phase inlet is L 1 = 2 W d upstream from the orifice, and the liquid phase inlet is L 2 = 2 W c 2 from the constriction so that the flow is fully developed before entering the orifice. The outlet channel is located at L 3 = 15 W d from the end of the exit expansion angle to ensure that the two-phase flow is fully developed in the downstream microchannel. The continuous phase inlet constriction width, outlet width, constriction angle, and expansion angle are selected as design variables, and the design variables and design space are shown in Table 1.
In this paper, water was selected as the continuous phase (its density and viscosity coefficient is ρ c = 998.2 kg/m3,   μ c = 1.003 mPa·s), nitrogen is the dispersed phase (its density and viscosity coefficient is ρ d = 1.138 kg/m3,   μ d = 0.01663 mPa·s), and the surface tension coefficient between the two phases is σ = 0.072 N/m. The Ca number ( Ca = μ υ / σ ) is defined for the aqueous phase (continuous phase), where all physical properties of both phases remain constant, so the Ca number varies only with the flow rate. The straight passageway transports the gas phase (N2), and the liquid phase flow (water) penetrates into the main passageway in the connection area. At the junction, the gas phase transforms into bubbles. At the end of the fragmentation process, the generated bubbles flow in the straight channel to the outlet.

2.2. Numerical Model

The Reynolds number Re is less than 10 for all cases in the flow-focusing system [27]. Therefore, the flow can be considered a laminar flow. In the simulations in this paper, the VOF method was used to track the interface between the continuous and dispersed phases. Since the density and viscosity of the two-phase fluid are constant, it can be regarded as an incompressible isotropic Newtonian fluid, so the Navier–Stokes equations describing the conservation of its momentum control equations are shown below.
· ν = 0
ν t + ( ν · ) ν = 1 ρ P + μ ρ ( ν + ν T ) + 1 ρ F
where ν is the velocity vector of the fluid; t is the time; P is the pressure; ρ and μ are the density and viscosity of the fluid; ( ν · ) ν is the inertial force per unit volume of fluid; P is the pressure gradient per unit volume of fluid; μ ( ν + ν T ) denotes the viscous force per unit volume of fluid; and F is the momentum source term associated with the surface tension.
The capture of the two-phase interface motion can be characterized by calculating the distribution of the continuous and dispersed phase volume fractions φ c and φ d in a grid cell, where φ c = 1 and φ d = 0 for the continuous phase and φ d = 1 and φ c = 0 for the dispersed phase. Therefore, the two intersecting interfaces in a calculation cell depend on the magnitude of φ c and φ d , taking values between 0 and 1. In the two-phase mixing unit, the calculation of the two-phase mixing density and viscosity in Equations (3) and (4) can be obtained from Equations (5) and (6).
ρ = φ d ρ d + ( 1 φ d ) ρ c
μ = φ d μ d + ( 1 φ d ) μ c
In addition, the volume fraction φ d can be obtained by solving the continuous equation for the volume fraction:
φ d t + ν φ d = 0
In the calculation grid adjacent to the wall surface, the two phases of liquid and gas are in contact with the wall surface at a fixed angle, and a contact angle model can be used to characterize the effect of wall wettability. The contact angle θ can be defined as the angle between the liquid–gas and solid–liquid phase interfaces at the intersection of the solid wall surface, liquid phase, and gas phase. The normal unit vector in the calculation cell close to the wall surface in this model can be obtained by processing Equation (8):
n ^ = n ^ w cos θ + t ^ w sin θ
where n ^ w and t ^ w represent the wall unit normal vector and tangential vector, respectively.
In addition to the control equations, the relevant boundary conditions of the flow-focusing microchannel are as follows:
  • Inlet: Both phases are uniform inlet, and the flow velocity of the continuous and dispersed phase inlets are set to 1 m/s and 0.2 m/s, respectively, to ensure that microbubbles can be generated in the dripping regime.
  • Outlet: The ambient pressure is used as the outlet fluid’s reference pressure, equal to 0 Pa.
  • Wall: The contact angle is set to 50°, and the channel wall is set to no-slip boundary condition.
In this study, a double-precision pressure solver was used for numerical simulations. The pressure implicit algorithm with operator splitting (PISO) is used to implement the pressure level coupling, which is based on a higher order of the approximate relationship between pressure and velocity corrections. Pressure staggering option (PRESTO) is chosen as the pressure interpolation algorithm; the momentum term is in the second-order windward form, and the volume fraction is in the geo-reconstruct format. Due to the high accuracy of curvature calculation, the geometric reconstruction format is used for interface interpolation. The convergence criterion is set to 1 × 10−6. In the simulation, the variable step size is automatically adjusted according to the criterion that the global Courant number is less than 0.25, thus maintaining the stability and convergence of the solution.

2.3. Grid Independence and Model Validation

Usually, the smaller the grid is, the more accurate the calculation results are; however, the smaller the grid is, the more computational resources are consumed. In order to find the most economical mesh size, the unstructured mesh is used to mesh the 2D geometric model of the flow-focusing microchannel. Since different design variables affect the number of meshes, seven scale sizes ( Δ h = W d / 10 , W d / 15 , W d / 20 , W d / 25 , W d / 30 , W d / 35 , W d / 40 ) are selected for the meshing from cases 1 to 7. As shown in Figure 2, the variations between the minimum and maximum values of the microbubble diameter and microbubble detachment time are only 0.24% and 0.55%, respectively, from cases 4 to 7, indicating that the grid cells meet sufficient accuracy requirements. Therefore, to balance the accuracy of the calculation results and computational cost, Δ h = W d / 25 grid cells are chosen for meshing the geometric channels in the subsequent simulations.
In addition, to verify the accuracy of the current numerical simulation calculation method, the numerical simulation results were compared with the experimental results in [28], as shown in Figure 3. The figure shows that the microbubble formation results obtained from the numerical simulation agree well with the experimental results. The average error between the simulated and experimental values is 6.27%, which is within the acceptable range. Therefore, it can be concluded that the numerical simulation method used in this paper is accurate and feasible for characterizing the two-phase flow patterns within the flow-focusing microchannels.

3. Research Methods

3.1. Design of Sampling Point

The Latin hypercube sampling method was proposed in 1979 by M.D. Mc Kay, R.J. Beckman, and others as a constrained and homogeneous sampling method [29]. Its purpose is to ensure that the sampling points cover all sampling areas. Suppose that the sampling size is n ; the number of variables is m ; x i is the variable, and each variable is divided into n intervals of equal probability,   x i , l = x i , 0 < x i , 1 < < x i , n = x i , u ; and x i , l and x i , u are the upper and lower bounds of the variable x i , respectively. Extract a random sample from each interval of the variable x i , combining the n points extracted from x 1 with the n points in x 2 in a random and nonrepetitive manner until finally combining them with the n sample points in the variable x m . The sample points are designed according to the minimum distance maximum criterion, and the minimum distance d should satisfy
min ( d ) = m a x ( m i n ( d i j ) )
where d i j is the distance between sampling points x i and x j :
d i j = μ = 1 m | x i μ x j μ | 2
The study used a total of four design variable parameters. Latin hypercube sampling was used to choose the 40 sampling sites that would provide a large enough sample of data. After the initial simulation test, a sampling point matrix of 4 × 40 was obtained, with each row representing a set of test parameters and each column representing a variable. The sample of the collected data obtained are shown in Table A1 in Appendix A.

3.2. Surrogate Model

A radial basis function is a neural network that includes an input layer of all influences, a hidden layer, and an output layer consisting of all responses [30,31,32]. The transfer function between the first two layers is a radial basis function, whereas the transfer function between the hidden and output layers is linear. RBF neural networks are characterized by fast training speed and compact networks.
A surrogate model with design variable parameters ( W ¯ c 1 , W ¯ o u t , α , β ) and the objective function after the response of each sample point was developed, where the optimization metrics included the diameter ( D ) and frequency ( F ). In this work, the four parameters of the design variables were used as variables in the input layer, and the response of the output layer was obtained by transforming the hidden layer. The RBF topology of this work is shown in Figure 4.

3.3. NSGA-II

When multiple objectives exist, especially competing for objective functions, single-objective optimization methods cannot meet the requirements. Therefore, there is a need for a multi-objective optimization method that can simultaneously optimize multiple objective functions. Genetic algorithms have good applications in solving multi-objective problems.
Kalyanmoy Deb et al. proposed an efficient and fast global optimization algorithm in 2001 inspired by the Darwinian evolutionary theory [33]. The nondominated sorting genetic algorithm II (NSGA-II) randomly generates the initial population based on the objective function and then performs nondominated sorting and crowding distance calculation. The new generation is obtained through crossover, selection, and mutation by genetic algorithms. These progeny populations are combined with the initial parent population for nondominated sorting and calculation of the crowding distance for each individual. The best individuals are selected based on sorting and crowding distance to form the initial population of the next generation. This is a cyclic process until a certain number of evolutionary generations are reached, or the convergence criterion is satisfied, and the whole optimization process ends.
NSGA-II has good robustness and convergence. It ensures a uniform distribution of non-inferior optimal solutions, which are less likely to fall into local optima during the search process. It is especially suitable for dealing with nonlinear problems that are difficult to be solved by traditional search methods. In this paper, flow-focusing microchannel optimization was a complex mixed-variable nonlinear problem. Therefore, this paper used the NSGA-II to solve the multi-objective optimization problem. The flow chart of multi-objective optimization based on RBF and NSGA-II is shown in Figure 5.
Suppose f 1 ( X ) = D ,   f 2 ( X ) = F , and the objective function satisfies the following constraints:
{ M i n i m i z e f 1 ( x 1 , x 2 , x 3 , x 4 ) M a x i m i z e f 2 ( x 1 , x 2 , x 3 , x 4 )
subject   to { 0.5 x 1 1.5 1.5 x 2 3.5 10 x 3 30 5 x 4 15

3.4. TOPSIS

In multi-objective optimization, selecting the final optimal solution from the Pareto frontier is important. The ideal point of the Pareto frontier is the optimal value for each objective, and the non-ideal point is the worst value for each objective. The technique for order preference by similarity to an ideal solution (TOPSIS) method can select the final optimal solution from the Pareto front [34].
The dimensions and scales of different objective spaces in multi-objective optimization problems should be unified, since the TOPSIS method used utilizes the Euclidean dimensionless objective X i j n , which is defined as
X i j n = X i j i = 1 m ( X i j ) 2 2 ( i = 1 , 2 , , m ; j = 1 , 2 , , n )
where X i j is the point on the Pareto front, m is the solution, and n is the target number.
In the TOPSIS method, the spatial distance between each solution on the Pareto boundary and the ideal point represented by d i + is determined as follows:
d i + = j = 1 n ( X i j X j i d e a l ) 2 2
where F j i d e a l is the ideal value of the jth objective obtained in the single-objective optimization.
The solution distance from the non-ideal point is denoted by d i . Therefore,
d i = j = 1 n ( X i j X j N o n i d e a l ) 2 2
where F j N o n i d e a l is the non-ideal value of the jth objective obtained in the single-objective optimization.
The calculation of relative intimacy is shown as follows:
r i = d i d i + + d i
In the TOPSIS method, all solutions are ranked, and the final solution required is the solution with the maximum value; therefore, the compromise of the final selected solution is
A b e s t = m a x ( r i )

4. Results and Discussion

4.1. Surrogate Model Accuracy Analysis

The RBF neural network training model uses 30 of the 40 sets of data from the sample space that the Latin hypercube sampled, while the post-training prediction model uses the remaining 10 sets. This post-training prediction model then necessitates an analysis of the microchannel design variables and the fitting accuracy of the objective function. The coefficient R 2 in Equation (18) was used in the paper to assess the fitting accuracy of the surrogate model. The maximum value of R 2 is 1, and the closer its value to 1, the better the fitting accuracy. The approximate model R 2 ≥ 0.9 was generally considered for engineering applications [35].
R 2 = 1 i = 1 M ( y i y ^ i ) 2 i = 1 M ( y i y ¯ i ) 2
where y i is the ith simulated value of the model, y ^ i is the ith predicted value of the model, y ¯ i is the mean of the simulated values, and M is the number of samples.
It can be seen in Figure 6 the relationship between the simulated and predicted values of the flow-focusing microchannels by linearly fitting them; the evaluation coefficients R 2 are 0.940 and 0.931 for the two targets, respectively, based on the calculation of the above equations. The RBF predictions of the two objective functions are close enough to the computational fluid dynamics (CFD) test values, indicating that the developed RBF model is reliable and sufficient to design flow-focusing microchannel structures.
In sensitivity analysis, the input variables of Pearson and Spearman linear correlation analysis are correlated with a certain objective function [36,37]. The results of the sensitivity analysis of the model are shown in Figure 7. As can be seen from Figure 7, the constriction angle is positively correlated with the microbubble diameter. It has the highest effect on the microbubble diameter (79%), followed by the expansion angle (69%) and the continuous phase inlet constriction width (46%). The outlet width has the least effect on the diameter (−19%), negatively correlated with the diameter.
For the generation frequency of microbubbles, the constriction angle was the most important influencing factor (−87%), negatively correlated with the generation frequency. This is followed by the expansion angle (−67%) and the continuous phase inlet constriction width (−38%). The outlet width has the least effect on the frequency (31%), positively correlated with the generation frequency.

4.2. Optimization Design Results

The optimization aims to obtain the best target values (highest generation frequency and minimum diameter). The diameter is measured by the microbubbles produced. The formation frequency is calculated from the time difference between the microbubbles’ separation from the orifice and the moment of the previous bubble separation. It can be seen from Figure 7 that all objective functions vary with different design parameters, and there exist design parameters corresponding to the optimal objective function. The Pareto optimal solution obtained by NSGA-II is necessary to represent when the problem to be addressed is targeted [38]. To detect the Pareto optimal frontier, the NSGA-II algorithm is used to calculate the optimal value. Table 2 shows the initial settings of the NSGA-II algorithm.
Figure 8 shows the Pareto optimal frontier diagram between the objective functions (red dots are optimal solutions). Each frontier solution represents the trade-off between the two objective functions. It is worth noting that the points on the Pareto front are optimal and are not dominated by each other. Each point corresponds to a specific structural parameter and optimization objective. To select the appropriate design point from many optimized design points to meet the requirements of practical engineering applications, TOPSIS is used to obtain the structural parameters corresponding to the final optimization results from the Pareto front, as shown in Figure 8. The objective values of the compromise solution are 1.6877 and 177.507, and the corresponding combinations of input parameters are W ¯ c 1 = 0.9, W ¯ o u t = 2.6, α = 10, and β = 5. The results of the error analysis of the predicted results with the numerical simulation show that the relative errors of the diameter and frequency are 0.43% and 1.17%, respectively, as shown in Table 3. This indicates that the predicted values are consistent with the CFD calculation results. The results show that the multi-objective optimization method based on RBF, NSGA-II, and TOPSIS proposed in this paper is effective.

4.3. Performance Evaluation of Flow-Focusing Microchannels

Most research efforts have focused on fluid properties and operating conditions for traditional flow-focusing microchannels. In this paper, improved performance was achieved by enhancing the flow-focusing microchannel, and the performance of the improved microchannel is compared with that of the traditional microchannel [22]. To be consistent with the structural parameters of this paper, the orifice width of the traditional microchannel is 2 μm, and the continuous phase inlet width is set to 4 μm.

4.3.1. Flow Regimes in Flow-Focusing Microchannels

Figure 9 and Figure 10 show the time evolution of the fluid dynamics information during microbubble generation in two flow-focusing microchannels. The figures include two flow states: dripping and squeezing regime, showing different fluid flow properties and microbubble generation performance. Microbubble generation mechanisms are essential for understanding both microchannels; therefore, these effects are discussed in both microchannels using flow field profiles.
Figure 9 shows the time evolution of the hydrodynamic information during microbubble generation in the improved microchannel, where the flow regime is in the dripping regime. The process of microbubble generation in the dripping regime can be divided into three stages: growing, necking, and pinch-off. The initial phase of the first microbubble generation can be seen in the dripping regime (see t = 5.96–9.26 μs in Figure 9), where the front part of the dispersed phase threads mainly stretches toward the main downstream channel without significant expansion in the direction of the side channels. As the thread continues to extend into the main downstream channel, the transverse interface of the thread becomes concave. The necking phase (see t = 9.26–10.06 μs in Figure 9) begins with the formation of a neck, which contributes to necking and thread extension due to the enhanced squeezing of the continuous phase under the action of the constriction inlet, in combination with the viscous shearing effect of the continuous phase. After the necking stage, a pinch-off stage occurs (see t = 10.06–10.36 μs in Figure 9), where high pressure appears on both sides of the neck, and the pressure difference is generated with the neck further promoting the necking. It can be seen that the fragmentation and generation of microbubbles occur at the junction of two phases. Due to the improved continuous phase inlet channel in this paper, it will be found that the microbubble diameter is smaller than the channel width. In this flow state, microbubbles’ generation is consistent, stable, and uniform.
Figure 10 shows the time evolution of the hydrodynamic information during microbubble generation in a traditional microchannel where the flow regimes are in the squeezing regime. The process of microbubble generation in the squeezing regime can be divided into four stages: expanding, squeezing, necking, and pinch-off. During the expanding stage (see t = 10.09–21.07 μs in Figure 10), the front part of the dispersed phase continuously swells with the continuous injection of the dispersed phase and remains spherical at the intersection. As the front area increases, the gap between the interface and corners of the main and side channels decreases, thus limiting the continuous phase flow into the main downstream channel. A squeeze will form in the dispersed phase, and the squeezing phase begins. During the squeezing stage (see t = 21.07–22.77 μs in Figure 10), the thread tip starts to thin and gradually moves downward under continuous phase squeezing, blocking the downstream exit of the intersection. The convex surface of the front gradually flattens out, followed by the necking stage. At this stage (see t = 22.77–23.37 μs in Figure 10), the thread plane starts to change to a concave surface, and a neck is formed at the intersection. Under the squeezing pressure of the continuous phase, the neck continuously shrinks while the thread stretches downstream and then enters the pinch-off phase (see t = 23.37–23.57 μs in Figure 10). As the continuous phase constantly squeezes the neck, the neck rapidly contracts, and pinch-off occurs. The fragmentation and generation of microbubbles occur near the intersection of the aqueous and gas phases. The generated microbubbles become clogged in the outlet channel and encounter higher resistance to move downstream.
The time-dependent pressure and velocity contours near the intersection allow us to understand the generation of microbubbles under different microchannels. It is clear from the pressure contours (Figure 9 and Figure 10) that the pressure in the main channel is decreasing in sequence. According to the Laplace pressure theory, the surface tension between the two phases causes the pressure to jump at the curved interface, and the internal pressure of the microbubble is higher than the internal pressure of the surrounding stable phase. It can be seen that the gas continuously enters the main channel; as the water phase constantly squeezes the gas, there is a pressure drop; this pressure drop plays an essential role in the fragmentation of microbubbles. In the microbubble pinch-off phase, a high-pressure region can be observed near the separation point (see t = 23.57 μs in Figure 10b).
For the velocity contours, the maximum velocity in the improved and traditional microchannels during fragmentation occurs during the pinch-off stage (see t = 10.36 μs in Figure 9c and t = 23.37 μs in Figure 10c). Comparing the two microchannels, it can be seen that the maximum velocity reached in the traditional microchannel is 13.9 m/s. The velocity at the orifice is higher in the improved microchannel, and the outlet channel is higher in the traditional microchannel.
The evolution of the thread profile helps understand the fluid dynamics behind the generation of microbubbles in different microchannels. Figure 11 shows the variation of tip length over time during the formation of microbubbles in both microchannels; L t i p is the distance from the threaded tip to the end of the main inlet channel in the connection region (see Figure 1).
For the improved microchannel, the L t i p tip increases as the dispersed phase threads stretch downstream during the growth phase in the dripping regime. At t = 8.26 μs, it can be seen that the enhancement of the continuous phase squeeze leads to an acceleration of the downstream motion of the thread tip. After the first microbubble is generated, a significant retraction of the remaining thread tip can be seen. After 10.46 μs, the tip of the dispersed phase gradually increases and enters the next cycle of microbubble generation. For the traditional microchannel, the change in the L t i p tip is similar to that of the improved microchannel, and the slope of the change in the L t i p tip is the same in both microchannels before the inflection point acceleration occurs. After 21.67 μs, the tip of the thread rapidly stretches downstream, and after the neck pinch-off, it can be seen that the tip of the remaining thread and the tail of the newly formed microbubble quickly retract due to the interfacial tension (t = 23.47–23.77 μs), after which the next cycle of microbubble generation begins. Compared with the traditional microchannel, the maximum L t i p and minimum retraction L t i p of the improved microchannel are both smaller.

4.3.2. Pressure Changes during Microbubble Formation

The contribution of pressure in different microchannels to the fragmentation mechanism during microbubble generation can be further investigated by observing the pressure gradient with time near the orifice [39]. The pressure variation near the orifice helps to understand the fundamental physical process of microbubble generation. Figure 12 shows the variation of the two-phase (dispersed phase and continuous phase) pressure with time for the two microchannels. P c is the pressure at the continuous phase channel connection point c, and P d is measured at point d before the main channel connection point.
In both microchannels, pressure changes exhibit similar behavior on the time scale. The improved microchannel (Figure 12a) has P d > P c at the very beginning, and as the gas enters, it will appear that P c > P d . Initially, the gas continuously enters the main channel, and P c remains approximately stable. As the gas continues to enter, the pressure accumulates, P c rises, and the pressure rise slows down at 9.26 μs when the microbubble neck begins to contract. At t = 10.36 μs, the gas ruptured to form microbubbles, followed by a rapid decrease in P c . For P d , first, the gas keeps entering the main channel, and P d keeps increasing. After the gas enters the crossover, it is subjected to resistance from the liquid phase, and the pressure basically remains constant. Before 9.26 μs, P d starts to be smaller than P c , and there is a necking trend. At 9.26 μs, the neck starts to constrict, and P d starts to decrease. P d rapidly increased after microbubble fragmentation, and then P c and P d started to fluctuate. In the traditional microchannel (Figure 12b), it can be seen that, initially, P c > P d , and as the gas continues to enter, P d first increases and then slowly decreases. P c remains stable at first and then starts to increase until it is in line with P d . In the squeezing stage, P d = P c . After t = 22.77 μs, P c starts to be larger than P d after the start of necking. At t = 23.57 μs, the gas ruptures to form microbubbles, and P d rapidly rises and P d > P c . After that, P d and P c fluctuate. It is clear that the average pressure value in the improved microchannel is lower than in the traditional one.
The pressure fluctuations at the inlet of the main channel shown at the p-point during microbubble generation in traditional and improved microchannels are given in Figure 13. The curve indicates that the pressure remains constant at the beginning and then periodically fluctuates. Usually, during the formation of each microbubble, the pressure increases and then decreases, so it can be seen from the graph that there are many peaks in the pressure fluctuations, and each peak represents the production of one microbubble. The pressure at the p-point in the improved microchannel is lower than that in the traditional microchannel, and the pressure more rapidly changes in the improved microchannel.

4.3.3. Microbubble Diameter and Generation Frequency

Figure 14a shows the variation of the microbubble diameter for different capillary numbers Ca , where D denotes the area-equivalent circle diameter. The figure shows that the microbubble diameter decreases with increasing Ca number. This reduction is due to the injection of the continuous phase at a higher rate, resulting in smaller microbubbles. It is clear from the figure that the microbubbles under the traditional flow-focusing microchannel are larger than the channel width, whereas the microbubbles under the improved flow-focusing microchannel are smaller than the channel width. The smallest microbubbles were from the improved flow-focusing microchannel at Ca = 0.0167, whereas the largest microbubbles were observed in the blocked state under the traditional flow-focusing microchannel at Ca = 0.0056. By fitting the microbubble diameter with the number of capillaries, the slope of the linear fitting function of the traditional microchannel is smaller than that of the improved microchannel. It can be seen that the diameter of the microbubble generated by the improved microchannel will be smaller with the increase in the number of capillaries.
Figure 14b shows the frequency of microbubble generation for different capillary numbers. The generation frequency is inversely proportional to the microbubble fragmentation time, which is the moment when the gas phase is separated from the two lateral channels by the water phase. When the water flow is injected into the connection at a higher velocity, the higher the Ca number, the shorter the microbubble fragmentation time, and vice versa, the higher the microbubble generation frequency. The microbubble generation frequency was fitted to the capillary number. Moreover, the slope of the increase in generation frequency with increasing Ca number in the improved microchannel was much larger than that in the traditional microchannel, as seen from the fitted linear function. At the same Ca number, the generation frequency of the improved microchannel is greater than that of the traditional microchannel. At Ca = 0.0167, the maximum frequency is 200 kHz in the improved microchannel and 68.617 kHz in the traditional microchannel. At Ca = 0.0056, the minimum frequency in the improved microchannel is 90.09 kHz, and the minimum frequency in the traditional microchannel is 42.765 kHz. As shown in Table 3, the constriction of the continuous phase inlet enhanced the flow-focusing effect; the microbubble diameter of the improved microchannel decreased from 2.8141 μm to 1.6949 μm, and the microbubble generation frequency increased from 64.077 kHz to 175.438 kHz compared with the traditional microchannel when Ca = 0.0138.

5. Conclusions

In this paper, numerical simulations and optimization algorithms investigated the fluid flow characteristics of microbubble generation in an improved flow-focusing microchannel with a constricted continuous phase inlet. A two-dimensional numerical model of microbubble generation in a flow-focusing microchannel was developed using the VOF method. Four design variables were obtained by the Latin hypercube sampling method, and the RBF surrogate model was used to establish the relationship between the objective function and design variables. The R 2 of the microbubble diameter and generation frequency were 0.940 and 0.931, respectively, which indicated that the RBF surrogate model could respond well to the relationship between the design variables and corresponding objective function. The results of the sensitivity analysis showed that the contraction angle had the most significant effect on the diameter with a positive correlation, and the outlet width had the least impact on the diameter with a negative correlation. The contraction angle has the most significant effect on the generation frequency with a negative correlation, followed by the expansion angle, continuous phase inlet contraction width, and outlet width, where the outlet width is positively correlated with the generation frequency. Then the optimal design of the NSGA-II algorithm was carried out, and TOPSIS decision analysis was performed on the obtained Pareto solution set to obtain the final optimal result of the microbubble diameter of 1.6949 μm and generation frequency of 175.438 kHz. Finally, the fluid flow characteristics during microbubble generation in the improved and traditional flow-focusing microchannels were investigated. It was found that the flow regime of microbubbles generated by the improved microchannel is dripping, whereas the flow regime of microbubbles caused by the traditional microchannel is squeezing. The pressure at the orifice in the improved flow-focusing microchannel reaches its maximum earlier in the whole microbubble formation cycle, which contributes to the separation of microbubbles more than the traditional microchannel. Constriction of the continuous phase inlet enhances the flow-focusing effect. When the Ca numbers are the same, the generated microbubbles in the improved flow-focusing microchannels have smaller diameters and higher frequencies; the microbubble diameters in both microchannels decrease with increasing Ca numbers, and the generation frequencies increase with increasing Ca numbers. According to the fitted linear function, it is known that the slope of the decrease in the microbubble diameter with increasing Ca number and the slope of the increase in microbubble generation frequency with increasing Ca number are greater in the improved microchannel compared with the traditional one. We hope that the results of the study of flow-focusing microchannels will provide guidance for the design of microbubble generators for high frequencies. In this paper, the reliability of the microchannel simulations was confirmed by correct grid independence checks and model validation, and the analytical procedure and numerical results used help to provide guidance for the design of microbubble generators at high frequencies. The necessary experiments are important for the study of the microbubble generation characteristics, and we will make experiments a focus of future research.

Author Contributions

Conceptualization, Z.W.; methodology, W.D.; software, W.D., Y.F. and J.C.; validation, Y.F., J.W. and H.W.; investigation, Y.F., J.W. and H.W.; resources, Z.W.; writing—original draft preparation, W.D.; writing—review and editing, W.D.; supervision, Z.W. and J.C.; project administration, Z.W.; funding acquisition, Z.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (51875419), Open Foundation of the State Key Laboratory of Fluid Power and Mechatronic Systems (GZKF-202122), Natural Science Foundation of Hubei Province (2021CFB495), and Hubei Higher Education Research Project (2015233).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Latin hypercube sampling method for generating 40 data samples within CFD.
Table A1. Latin hypercube sampling method for generating 40 data samples within CFD.
No.Design VariablesObjective Functions
W ¯ c 1 W ¯ o u t α β D F
11.22.21191.8939138.888
21.41.526122.1013112.994
30.92.623121.7746156.250
41.12.51381.8194149.253
50.92.71061.7054173.913
60.83.430141.9273134.228
71.42.824102.0886112.994
80.93.518131.8610143.884
91.03.020141.8219149.253
100.82.317131.9165136.986
111.22.016101.9321136.986
121.32.124122.0200124.223
131.43.32382.0623114.942
140.83.32591.8178149.253
151.02.514141.7962151.515
161.41.71752.0995116.279
170.62.22851.8103152.671
180.51.62171.8769140.845
191.22.429151.9197133.333
201.53.11462.1888103.626
211.12.72291.8060151.515
221.02.025111.7564158.730
231.21.815121.9391134.228
240.81.828101.8421147.058
251.53.227152.2303101.010
260.63.02271.8623144.927
270.72.61981.8363147.058
280.72.41981.8535144.927
291.32.921112.0171124.223
301.31.61271.9769128.205
310.73.215101.9033138.888
320.51.91271.8604140.845
331.32.92061.9706128.205
340.62.81361.8337147.058
350.93.418131.8764143.884
360.71.929131.9444131.578
371.12.126111.8311149.253
381.01.716141.8022150.375
390.63.127111.9577130.718
401.12.31191.8286149.253

References

  1. Keller, M.W.; Feinstein, S.B.; Watson, D.D. Successful left ventricular opacification following peripheral venous injection of sonicated contrast agent: An experimental evaluation. Am. Heart J. 1987, 114, 570–575. [Google Scholar] [CrossRef]
  2. Wei, K.; Jayaweera, A.R.; Firoozan, S.; Linka, A.; Skyba, D.M.; Kaul, S. Quantification of myocardial blood flow with ultrasound-induced destruction of microbubbles administered as a constant venous infusion. Circulation 1998, 97, 473–483. [Google Scholar] [CrossRef]
  3. van Wamel, A.; Kooiman, K.; Harteveld, M.; Emmer, M.; ten Cate, F.J.; Versluis, M.; de Jong, N. Vibrating microbubbles poking individual cells: Drug transfer into cells via sonoporation. J. Control. Release 2006, 112, 149–155. [Google Scholar] [CrossRef]
  4. Dixon, A.J.; Dhanaliwala, A.H.; Chen, J.L.; Hossack, J.A. Enhanced intracellular delivery of a model drug using microbubbles produced by a microfluidic device. Ultrasound Med. Biol. 2013, 39, 1267–1276. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Duncanson, W.J.; Abbaspourrad, A.; Shum, H.C.; Kim, S.H.; Adams, L.L.; Weitz, D.A. Monodisperse gas-filled microparticles from reactions in double emulsions. Langmuir 2012, 28, 6742–6745. [Google Scholar] [CrossRef] [Green Version]
  6. Abbaspourrad, A.; Duncanson, W.J.; Lebedeva, N.; Kim, S.H.; Zhushma, A.P.; Datta, S.S.; Dayton, P.A.; Sheiko, S.S.; Rubinstein, M.; Weitz, D.A. Microfluidic fabrication of stable gas-filled microcapsules for acoustic contrast enhancement. Langmuir 2013, 29, 12352–12357. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Zhang, M.-J.; Wang, W.; Xie, R.; Ju, X.-J.; Liu, L.; Gu, Y.-Y.; Chu, L.-Y. Microfluidic fabrication of monodisperse microcapsules for glucose-response at physiological temperature. Soft Matter 2013, 9, 4150–4159. [Google Scholar] [CrossRef]
  8. Cosgrove, D. Ultrasound contrast agents: An overview. Eur. J. Radiol. 2006, 60, 324–330. [Google Scholar] [CrossRef]
  9. Zhao, Y.Z.; Liang, H.D.; Mei, X.G.; Halliwell, M. Preparation, characterization and in vivo observation of phospholipid-based gas-filled microbubbles containing hirudin. Ultrasound Med. Biol. 2005, 31, 1237–1243. [Google Scholar] [CrossRef] [PubMed]
  10. Soetanto, K.; Chan, M. Fundamental studies on contrast images from different-sized microbubbles: Analytical and experimental studies. Ultrasound Med. Biol. 2000, 26, 81–91. [Google Scholar] [CrossRef]
  11. Böhmer, M.R.; Schroeders, R.; Steenbakkers, J.A.M.; de Winter, S.H.P.M.; Duineveld, P.A.; Lub, J.; Nijssen, W.P.M.; Pikkemaat, J.A.; Stapert, H.R. Preparation of monodisperse polymer particles and capsules by ink-jet printing. Colloids Surf. A Physicochem. Eng. Asp. 2006, 289, 96–104. [Google Scholar] [CrossRef]
  12. Farook, U.; Stride, E.; Edirisinghe, M.J.; Moaleji, R. Microbubbling by co-axial electrohydrodynamic atomization. Med. Biol. Eng. Comput. 2007, 45, 781–789. [Google Scholar] [CrossRef]
  13. Stride, E.; Edirisinghe, M. Novel preparation techniques for controlling microbubble uniformity: A comparison. Med. Biol. Eng. Comput. 2009, 47, 883–892. [Google Scholar] [CrossRef] [PubMed]
  14. Wang, S.; Dhanaliwala, A.H.; Chen, J.L.; Hossack, J.A. Production rate and diameter analysis of spherical monodisperse microbubbles from two-dimensional, expanding-nozzle flow-focusing microfluidic devices. Biomicrofluidics 2013, 7, 14103. [Google Scholar] [CrossRef] [Green Version]
  15. Peng, H.; Xu, Z.; Chen, S.; Zhang, Z.; Li, B.; Ge, L. An easily assembled double T-shape microfluidic devices for the preparation of submillimeter-sized polyacronitrile (PAN) microbubbles and polystyrene (PS) double emulsions. Colloids Surf. A Physicochem. Eng. Asp. 2015, 468, 271–279. [Google Scholar] [CrossRef]
  16. Fu, T.; Ma, Y.; Funfschilling, D.; Zhu, C.; Li, H.Z. Squeezing-to-dripping transition for bubble formation in a microfluidic T-junction. Chem. Eng. Sci. 2010, 65, 3739–3748. [Google Scholar] [CrossRef]
  17. Xu, J.H.; Li, S.W.; Wang, Y.J.; Luo, G.S. Controllable gas-liquid phase flow patterns and monodisperse microbubbles in a microfluidic T-junction device. Appl. Phys. Lett. 2006, 88, 133506. [Google Scholar] [CrossRef]
  18. Lorenceau, E.; Sang, Y.Y.C.; Höhler, R.; Cohen-Addad, S. A high rate flow-focusing foam generator. Phys. Fluids 2006, 18, 097103. [Google Scholar] [CrossRef] [Green Version]
  19. Dollet, B.; van Hoeve, W.; Raven, J.P.; Marmottant, P.; Versluis, M. Role of the channel geometry on the bubble pinch-off in flow-focusing devices. Phys. Rev. Lett. 2008, 100, 034504. [Google Scholar] [CrossRef] [Green Version]
  20. Dietrich, N.; Poncin, S.; Midoux, N.; Li, H.Z. Bubble formation dynamics in various flow-focusing microdevices. Langmuir 2008, 24, 13904–13911. [Google Scholar] [CrossRef]
  21. Fu, T.; Ma, Y.; Funfschilling, D.; Li, H.Z. Bubble formation and breakup mechanism in a microfluidic flow-focusing device. Chem Eng Sci 2009, 64, 2392–2400. [Google Scholar] [CrossRef]
  22. Castro-Hernandez, E.; van Hoeve, W.; Lohse, D.; Gordillo, J.M. Microbubble generation in a co-flow device operated in a new regime. Lab Chip 2011, 11, 2023–2029. [Google Scholar] [CrossRef] [PubMed]
  23. Sheng, L.; Chen, Y.; Deng, J.; Luo, G. High-frequency formation of bubble with short length in a capillary embedded step T-junction microdevice. AIChE J. 2021, 67, e17376. [Google Scholar] [CrossRef]
  24. Jia, H.; Zhang, P. Bubble formation in viscous fluids by a microfluidic flow-focusing junction: A computational study. Microfluid. Nanofluidics 2020, 24, 85. [Google Scholar] [CrossRef]
  25. Zadeh, S.A.; Radespiel, R. Numerical Study on Bubble Formation in a Microchannel Flow-Focusing Device Using the VOF Method. In Proceedings of the ASME 2011 9th International Conference on Nanochannels, Microchannels, and Minichannels, Edmonton, AB, Canada, 19–22 June 2011; Volume 2, pp. 419–430. [Google Scholar] [CrossRef]
  26. Chekifi, T.; Boukraa, M.; Aissani, M. DNS using CLSVOF method of single micro-bubble breakup and dynamics in flow focusing. J. Visual.-Jpn. 2021, 24, 519–530. [Google Scholar] [CrossRef]
  27. Yu, W.; Liu, X.; Zhao, Y.; Chen, Y. Droplet generation hydrodynamics in the microfluidic cross-junction with different junction angles. Chem. Eng. Sci. 2019, 203, 259–284. [Google Scholar] [CrossRef]
  28. Fu, T.; Funfschilling, D.; Ma, Y.; Li, H.Z. Scaling the formation of slug bubbles in microfluidic flow-focusing devices. Microfluid. Nanofluidics 2010, 8, 467–475. [Google Scholar] [CrossRef]
  29. McKay, M.D.; Beckman, R.J.; Conover, W.J. A Comparison of Three Methods for Selecting Values of Input Variables in the Analysis of Output from a Computer Code. Technometrics 1979, 21, 239–245. [Google Scholar] [CrossRef]
  30. Tian, W.; Jiang, C.; Ni, B.; Wu, Z.; Wang, Q.; Yang, L. Global sensitivity analysis and multi-objective optimization design of temperature field of sinter cooler based on energy value. Appl. Therm. Eng. 2018, 143, 759–766. [Google Scholar] [CrossRef]
  31. Peng, Y.; Parsian, A.; Khodadadi, H.; Akbari, M.; Ghani, K.; Goodarzi, M.; Bach, Q.-V. Develop optimal network topology of artificial neural network (AONN) to predict the hybrid nanofluids thermal conductivity according to the empirical data of Al2O3—Cu nanoparticles dispersed in ethylene glycol. Phys. A Stat. Mech. Its Appl. 2020, 549, 124015. [Google Scholar] [CrossRef]
  32. Giwa, S.O.; Sharifpur, M.; Goodarzi, M.; Alsulami, H.; Meyer, J.P. Influence of base fluid, temperature, and concentration on the thermophysical properties of hybrid nanofluids of alumina–ferrofluid: Experimental data, modeling through enhanced ANN, ANFIS, and curve fitting. J. Therm. Anal. Calorim. 2021, 143, 4149–4167. [Google Scholar] [CrossRef]
  33. Deb, K.; Pratap, A.; Agarwal, S.; Meyarivan, T. A fast and elitist multiobjective genetic algorithm: NSGA-II. IEEE Trans. Evol. Comput. 2002, 6, 182–197. [Google Scholar] [CrossRef] [Green Version]
  34. Sayyaadi, H.; Mehrabipour, R. Efficiency enhancement of a gas turbine cycle using an optimized tubular recuperative heat exchanger. Energy 2012, 38, 362–375. [Google Scholar] [CrossRef]
  35. Bagherzadeh, S.A.; D’Orazio, A.; Karimipour, A.; Goodarzi, M.; Bach, Q.-V. A novel sensitivity analysis model of EANN for F-MWCNTs–Fe3O4/EG nanofluid thermal conductivity: Outputs predicted analytically instead of numerically to more accuracy and less costs. Phys. A Stat. Mech. Its Appl. 2019, 521, 406–415. [Google Scholar] [CrossRef]
  36. Hauke, J.; Kossowski, T. Comparison of Values of Pearson’s and Spearman’s Correlation Coefficients on the Same Sets of Data. Quageo 2011, 30, 87–93. [Google Scholar] [CrossRef] [Green Version]
  37. Wu, H.; Bagherzadeh, S.A.; D’Orazio, A.; Habibollahi, N.; Karimipour, A.; Goodarzi, M.; Bach, Q.-V. Present a new multi objective optimization statistical Pareto frontier method composed of artificial neural network and multi objective genetic algorithm to improve the pipe flow hydrodynamic and thermal properties such as pressure drop and heat transfer coefficient for non-Newtonian binary fluids. Phys. A Stat. Mech. Its Appl. 2019, 535, 122409. [Google Scholar] [CrossRef]
  38. Tucci, M.; Crisostomi, E.; Giunta, G.; Raugi, M. A Multi-Objective Method for Short-Term Load Forecasting in European Countries. IEEE Trans. Power Syst. 2015, 31, 3537–3547. [Google Scholar] [CrossRef] [Green Version]
  39. Bashir, S.; Rees, J.M.; Zimmerman, W.B. Investigation of pressure profile evolution during confined microdroplet formation using a two-phase level set method. Int. J. Multiph. Flow 2014, 60, 40–49. [Google Scholar] [CrossRef]
Figure 1. Schematics for improved flow-focusing microchannel.
Figure 1. Schematics for improved flow-focusing microchannel.
Micromachines 13 01776 g001
Figure 2. Grid independence test.
Figure 2. Grid independence test.
Micromachines 13 01776 g002
Figure 3. Comparison of numerical results of dimensionless microbubble size with experimental results.
Figure 3. Comparison of numerical results of dimensionless microbubble size with experimental results.
Micromachines 13 01776 g003
Figure 4. RBF neural network topology.
Figure 4. RBF neural network topology.
Micromachines 13 01776 g004
Figure 5. Multi-objective optimization flow chart.
Figure 5. Multi-objective optimization flow chart.
Micromachines 13 01776 g005
Figure 6. Comparison between prediction and reality: (a) D and (b) F.
Figure 6. Comparison between prediction and reality: (a) D and (b) F.
Micromachines 13 01776 g006
Figure 7. Sensitivity analysis.
Figure 7. Sensitivity analysis.
Micromachines 13 01776 g007
Figure 8. Pareto optimal front.
Figure 8. Pareto optimal front.
Micromachines 13 01776 g008
Figure 9. Time evolution of hydrodynamic information during microbubble generation in improved microchannel (Ca = 0.0139, Ud = 0.2 m/s): (a) evolution of the interface, (b) evolution of pressure field, and (c) evolution of velocity field.
Figure 9. Time evolution of hydrodynamic information during microbubble generation in improved microchannel (Ca = 0.0139, Ud = 0.2 m/s): (a) evolution of the interface, (b) evolution of pressure field, and (c) evolution of velocity field.
Micromachines 13 01776 g009
Figure 10. Time evolution of hydrodynamic information during microbubble generation in traditional microchannel (Ca = 0.0139, Ud = 0.2 m/s): (a) evolution of the interface, (b) evolution of pressure field, and (c) evolution of velocity field.
Figure 10. Time evolution of hydrodynamic information during microbubble generation in traditional microchannel (Ca = 0.0139, Ud = 0.2 m/s): (a) evolution of the interface, (b) evolution of pressure field, and (c) evolution of velocity field.
Micromachines 13 01776 g010
Figure 11. Variation of tip length with time during formation of microbubbles in two microchannels.
Figure 11. Variation of tip length with time during formation of microbubbles in two microchannels.
Micromachines 13 01776 g011
Figure 12. Pressure curve: (a) improved microchannel; (b) traditional microchannel (Ca = 0.0139, Ud = 0.2 m/s).
Figure 12. Pressure curve: (a) improved microchannel; (b) traditional microchannel (Ca = 0.0139, Ud = 0.2 m/s).
Micromachines 13 01776 g012
Figure 13. Variation of pressure at p-point with time in traditional and improved microchannels.
Figure 13. Variation of pressure at p-point with time in traditional and improved microchannels.
Micromachines 13 01776 g013
Figure 14. (a) Microbubble diameter (D) and (b) generation frequency (F) as a function of capillary number (Ca) (Ud = 0.2 m/s).
Figure 14. (a) Microbubble diameter (D) and (b) generation frequency (F) as a function of capillary number (Ca) (Ud = 0.2 m/s).
Micromachines 13 01776 g014
Table 1. Design variables and design space.
Table 1. Design variables and design space.
Limit W ¯ c 1 W ¯ o u t α β
Upper1.53.53015
Lower0.51.5105
Table 2. Parameters setting of NSGA-II.
Table 2. Parameters setting of NSGA-II.
ParametersValues
Population size12
Number of generations100
Crossover probability0.9
Crossover distribution index10.0
Mutation distribution index20.0
Table 3. Comparison between numerical and RBF predicted values.
Table 3. Comparison between numerical and RBF predicted values.
DesignDiameter (μm)Frequency (kHz)
Improved microchannel1.6949175.438
Traditional microchannel2.814164.077
Relative error (%)
Im (simulation and predictive)0.431.17
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Z.; Ding, W.; Fan, Y.; Wang, J.; Chen, J.; Wang, H. Design of Improved Flow-Focusing Microchannel with Constricted Continuous Phase Inlet and Study of Fluid Flow Characteristics. Micromachines 2022, 13, 1776. https://doi.org/10.3390/mi13101776

AMA Style

Wang Z, Ding W, Fan Y, Wang J, Chen J, Wang H. Design of Improved Flow-Focusing Microchannel with Constricted Continuous Phase Inlet and Study of Fluid Flow Characteristics. Micromachines. 2022; 13(10):1776. https://doi.org/10.3390/mi13101776

Chicago/Turabian Style

Wang, Zhaohui, Weibing Ding, Yiwei Fan, Jian Wang, Jie Chen, and Hongxia Wang. 2022. "Design of Improved Flow-Focusing Microchannel with Constricted Continuous Phase Inlet and Study of Fluid Flow Characteristics" Micromachines 13, no. 10: 1776. https://doi.org/10.3390/mi13101776

APA Style

Wang, Z., Ding, W., Fan, Y., Wang, J., Chen, J., & Wang, H. (2022). Design of Improved Flow-Focusing Microchannel with Constricted Continuous Phase Inlet and Study of Fluid Flow Characteristics. Micromachines, 13(10), 1776. https://doi.org/10.3390/mi13101776

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop