Computational fluid dynamics simulation of two

Apr 13, 2024

Scientific Reports volume 13, Article number: 9483 (2023) Cite this article

763 Accesses

1 Citations

Metrics details

In the current research work, the flow behavior of a liquid–liquid extraction (LLE) process in a serpentine microchannel was analyzed. The simulation was performed using a 3D model and the results were found to be consistent with experimental data. The impact of the flow of chloroform and water on the flow model was also examined. The data indicate that once the aqua and organic phases flow rates are low and similar, a slug flow pattern is observed. However, as the overall flow rate raises, the slug flow transforms into parallel plug flow or droplet flow. An increment in the aqua flows while maintaining a constant organic phase flow rate results in a transition from slug flow to either droplet flow or plug flow. Finally, the patterns of flow rate in the serpentine micro-channel were characterized and depicted. The results of this study will provide valuable insights into the behavior of two-phase flow patterns in serpentine microfluidic devices. This information can be used to optimize the design of microfluidic devices for various applications. Furthermore, the study will demonstrate the applicability of CFD simulation in investigating the behavior of fluids in microfluidic devices, which can be a cost-effective and efficient alternative to experimental studies.

The use of two-phase liquid–liquid (LL) systems is prevalent in chemical treatment, for instance, polymerization, nitration, chlorination, and reactive and solvent extraction1,2,3,4,5. Those procedures are mostly hampered by transport limitations, such as small mass transfer rates6,7,8. To overcome these limitations, miniaturization has been recognized as a promising method of process intensification, by reducing transport resistance and increasing transport rates9,10,11. The utilization of micro-spaces in devices can result in high heat and mass transfer rates12,13,14,15,16,17. The higher interfacial zone-to-volume fraction in micro-scale binary schemes compared to macro-scale systems results in enhanced heat and mass transfer rates and increased process efficiency, which can be higher by an order of magnitude compared to conventional systems. Additionally, the ease of scale-up developed safety, and reduced inventory requirements, specifically for systems using risky and exclusive chemicals, make microfluidic devices appropriate for a broad range of applications. The effectiveness of a specific system in LL microchannels depends greatly on the flow schemes of the two non-miscible liquids18,19,20,21.

Microfluidic flow patterns refer to the behavior of the fluid in microscale channels or devices. Three main flow, parallel, droplet and slug flow, occurs in microfluidic systems. Flow maps show graphically of these main flows versus the flow rate of two phases. Understanding microfluidic flow patterns is important for designing and optimizing microfluidic devices for specific applications. By controlling the flow pattern, researchers can manipulate the behavior of fluids in microscale channels and develop devices that can perform precise chemical reactions, separations, and detections22,23,24.

Several LL flow patterns were scrutinized in microfluidic tools based on factors such as micro-channel size and shape, physical characteristics of the liquids (for instance viscosity and surface tension), flow rate, the flow ratio of the liquids, and the wetting behavior of the micro-channel walls25,26,27. The maximum usual LL flow patterns in two-phase micro-channels include slug flow, plug flow, and droplet flow. Slug flow is favored for numerous systems due to the interior rotation inside the slugs of two phases and the diffusion among the contiguous slugs. Nonetheless, comprehensive phase split-up inside the micro-fluidic tool remains a challenge in slug flow regimes. Slug hydrodynamics, such as slug length and speed, is of significant significance as they affect the performance of microfluidic devices28,29,30,31.

The researchers have put forward various scaling laws to estimate the length of the slug in a two-phase liquid–liquid flow32,33,34,35. It has been noted in previous studies that the speed of the slug is directly related to the overall speed of the fluid36,37,38. In plug flow, separating the liquid into two phases is possible on a chip because of the creation of a stable LL interface down the center of the micro-channel, which could be achieved by using a Y-split at the microchannels exit . However, the method of transport in plug flow is limited to diffusion, leading to low system effectiveness compared to processes that utilize segmented flow39,40,41. Several research projects were performed to create flow regime maps42,43,44. These maps could be sketched according to the flow rate of binary systems, but they cannot be used as a general representation as they do not take into account all factors that impact the flow sketches. Therefore, researchers have suggested dimensionless quantities, such as Re, We, and Ca, to create general flow configuration maps45. Various mixtures of these dimensionless quantities were utilized as coordinates for generalized flow diagrams. Waelchli et al.46 employed the Buckingham Pi-theorem to find overall flow behavior for gas–liquid flows and proposed an association of Reynolds (Re) and Weber (We) numbers to generalize experimental data. Meanwhile, Cao et al.36 examined the LL flow schemes in a non-circular glass micro-channel and created comprehensive flow model transitions according to the force analysis. It has also been suggested to utilize Re and We numbers to estimate flow model transitions. Yagodinitsyna et al.47 employed an easy dimensionless analysis for LL flow to find a universal factor for flow pattern maps. They proposed the We number multiplied by the Ohnesorge number (Oh) as a novel factor to simplify the flow models of their analyzed methods. The units of physical variables used in the dimensionless analysis and The dimensionless numbers are in Tables 1 and 2.

Darekar et al.43 investigated the flow models of standard removal schemes in Y-junction micro-channels and assessed the effectiveness of using Re, Capillary (Ca), We, and Weber multiplied by (We Oh) numbers for general flow map report. The last two quantities provided the best outcomes. It is crucial to conduct a dimensionless study to determine the appropriate grouping of dimensionless quantities to create a general flow diagram. Despite being assumed a guaranteeing alternative to traditional methods, the low output of a particular micro-channel remains a challenge, requiring the scaling-up of microfluidic devices for industrial manufacture. An answer to this issue is to develop multiple micro-channels in parallel, increasing the output while avoiding complications of traditional scaling methods48. To the best of our knowledge, there is little information on the flow model of two non-miscible mixtures in a parallel microfluidic tool, with the majority of earlier studies examining parallelized drop-based microfluidic tools for creating even emulsion drops49,50,51,52,53. Kassid et al.54 examined the flow dissemination of a LL scheme in six distinct capillaries to evaluate the mass transfer efficiency of the kerosene (+ acetic acid)—water scheme.

This study used CFD to determine the flow map in a serpentine microchannel. three regime flow in micro-channel, droplet, slug, and plug flow, were developed. This study focuses on the micro-channel to reveal their flow map at various flow rates. Furthermore, the multi-phase flow interior of the micro-channel was calculated mathematically to distinguish the flow model in the serpentine micro-channel. The outcomes of mathematical simulations were proved via experimental results. The received data from the simulations were compatible with the experimental results, indicating that the mathematical simulations accurately modeled the flow behavior in the microchannel..

The calculations for multiphase flows involving LL interfaces are done using the volume of fluid (VOF) technique55,56,57. This method is a popular technique used in CFD simulations to model the behavior of multiphase flows. The VOF method is a Eulerian method that tracks the interface between two immiscible fluids by solving the Navier–Stokes equations for each fluid phase while tracking the volume fraction of each phase in each computational cell.

The VOF method works by dividing the computational domain into a grid of small cells, and at each cell, the volume fraction of each fluid phase is tracked. The interface between the two fluids is defined as the region where the volume fraction of a phase changes from 0 to 1 or vice versa. This model is effective in monitoring the boundary between two phases that do not mix. In the VOF system, parameters such as pressure and speed are common to binary phases and reflect the average volume amounts. The correlations for conserving the mass and momentum of two incompressible, non-miscible liquids are volume-averaged and can be found in58,59,60,61,62,63,64,65:

In these correlations ρ, μ, u, ρgi, and fσ are the density of the fluid, the dynamic viscosity, the speed vector, the gravity force, and the exterior body forces, correspondingly.

The density of the mixture and viscosity is estimated via volume fraction averaging66,67,68,69,70.

where \(\mathrm{\alpha }\) is the volume percentage and digits belong to the phase. It is evident that the sum of the volume percentage of the specific phase is supposed to be hundred percent, as follows

The volume fraction in the interface ranges from 0 to 1, making it crucial to accurately track the interface within cells.

To do this, it is necessary to solve a continuity correlation for the volume percentage of unity or multiple phases. The continuity equation for the volume fraction of the i-th phase is expressed as follows:

At the gas–liquid interface, a pressure jump occurs due to the surface tension difference between both sides. This difference is considered in balance and its slope should match the added body force in the momentum balance. The pressure jump disjunction is estimated as described in66,71,72:

The wetting characteristics of the substance are analyzed by examining the water droplet contact angles on the glass side. The contact angles are measured using a standard goniometer. The physical Characteristics of the fluids are displayed in Table 3.

In this study, three-dimensional micro-channel geometry was chosen for the flow map analysis. The geometry and mesh are shown in Fig. 1. The microchannel has a smooth profile with a surface roughness average of 0.22 μm. The area of the cross-section is approximately 0.13 mm2 with a hydraulic diameter of 0.32 mm. The length of each microchannel (mixing part) is 102 mm. Three-dimensional structured elements were used for mesh generation due to their ability to handle complex geometries and reduce errors. The computational domain was initially established using components with standard side lengths of 5, 3, 2, and 1 µm (mesh numbers of 495,000, 950,000, 1,920,000, and 3,850,000). After explaining the flow field and volume percentage correlations, the slug length of the simulation was compared to the experimental data and plotted. The outcomes are presented in Table 4, with comparative errors in the estimation of the slug length for the 1 and 2 µm elements being less than 1%. Thus, networks with average lengths of 2 µm were chosen as the threshold for the numerical simulation.

(a) Geometry of micro-channel (b) creating the grid.

In order to assess the quality of the mesh used in the simulation, the aspect ratio and skewness of the mesh elements were analyzed. The aspect ratio measures the elongation of each mesh element, and a high aspect ratio can result in inaccurate simulations. The skewness measures the deviation from a regular shape, and a high skewness can lead to numerical instabilities. The results of the analysis showed that the aspect ratio of most elements was below 3, indicating that the mesh was not excessively elongated. Additionally, the skewness was within the acceptable range of less than 0.5. Based on these results, it was concluded that the mesh quality was acceptable and suitable for use in the simulation.

For the two-phase flow simulations, the uniform inlet velocity boundary conditions for both liquid phases are implemented. At the outlet, it is assumed that the Pressure outlet boundary condition for liquid and gas. At the walls, a no-slip boundary condition for liquid phases is imposed. The finite volume method was used for the numerical solutions. The SIMPLE algorithm was selected to compute the pressure–velocity coupling. The second-order upwind discretization was used for momentum. The convergence criterion of the residual errors was set to 10–4. For this purpose, the finite volume method of the ANSYS FLUENT software package was used to obtain numerical solutions.

The diagrams that depict the flow models of different LL schemes, known as flow maps, demonstrate how the flow rate affects the flow regime. Figures 2, 3, 4 and 5 present the different flow regimes—slug flow, droplet flow and plug flow—in a serpentine microfluidic channel as validated by the experimental findings of Asadi et al.45. The simulation outcomes match well with the experimental work.

The slug flow regime (a) simulation (b) experimental work45.

The droplet flow regime (a) simulation (b) experimental work45.

The plug flow regime (a) simulation (b) experimental work45.

3D chloroform volume fraction (a) slug flow (b) droplet flow (c) parallel flow.

Figure 2 illustrates the occurrence of slug flow once the flow rates of the aqua and organic phases are somewhat small and similar. In this figure, the flow rates are 100 and 100 µl per minute for the aqua and organic phases, respectively. As depicted, the organic phase enters the main channel first and takes up a large portion of its cross-section, causing the continuous phase to be blocked to a significant extent. These results in increased drag force being applied to the interface, leading to the gradual and complete entry of the organic phase into the major channel over time. The pressure gradient created in the forming clump and the drag force acting on the interface counteract the surface tension force, causing the dispersed phase to separate from the Y-shaped inlet of the micro-channel. With these two forces dominating the surface tension force, the dispersed phase becomes separated from the Y-junction and forms a clump. As the aqueous phase flows back into its designated entrance, the clump is completely detached and moves down the main channel. This process repeats in an alternating manner. The size of the clumps produced by the two-phase flow and the physical characteristics of the fluids used can be altered.

As the overall flow rate rises, slug flow transforms into either plug flow or droplet flow. The flow scheme depends on the form and organic phase rates. If the flow rate of the aqua phase is reduced while the organic phase rate is increased, the flow model will switch from slug flow to droplet flow, as shown in Fig. 3a with an aqua flow rate of 600 and organic flow rate of 30 µl per minute. If the aqua flow rate remains constant while the organic phase rate increases, the resulting flow model will always be plug flow, as illustrated in Fig. 4a with aqua and organic flow rates both at 500 µl per minute. Figure 5a-c show the chloroform volume fraction slug flow, droplet flow and parallel flow, respectively.

The flow maps in Fig. 6 demonstrate the impact of flow rate on flow regimes in a liquid–liquid system. At a moderate state, comparable flow rates for the aqua and organic phase, slug flow are observed. As the overall flow rate rises, the flow regime shifts to plug flow or droplet flow. If the aqua flow rate is maintained constant and the organic phase flow rate increases, the flow model transforms from slug to droplet or plug flow, depending on the type and flow rate of the organic phase. On the other hand, if the organic flow rate is increased whereas keeping the aqua flow rate uniform, the consequential flow model is permanent plug flow, as demonstrated in Fig. 6.

The flow models of chloroform-water based on computational work.

Figures 7, 8 and 9 show flow pattern maps of the two-phase systems with Re, Ca and We numbers as the coordinates based on computational work, respectively. As can be seen from Fig. 7, at higher Re numbers of the organic phase, the flow moves towards the parallel flow. At high Re numbers of the aqueous phase and low Re numbers of the organic phase, the flow approaches are droplet flow. At equal Re numbers as well as low Re numbers, we will have slug flow.

Flow pattern maps of the two-phase systems with Re as the coordinates based on computational work.

Flow pattern maps of the two-phase systems with Ca as the coordinates based on computational work.

Flow pattern maps of the two-phase systems with We as the coordinates based on computational work.

Figures 8 and 9 of the diagram are similar to what happened for the Re number, with the difference that in Fig. 9, the range of the Weber number covers a larger range than the two dimensionless numbers Re and Ca.

In this paper, a new model was presented, based on which the flow pattern in the serpentine microchannel can be predicted. In addition, this research examines the flow behavior of liquid–liquid extraction using chloroform and water in a serpentine microchannel. A 3D model was employed to calculate the flow behavior in the serpentine microchannel, and the results were compatible with experimental data. The impact of the flow rate of chloroform and water on the flow pattern was analyzed, revealing that slug flow occurs at low and comparable flow rates, but transitions to parallel or droplet flow as the total flow rate improves. Increasing the aqua flow rate while maintaining the organic phase flow rate stable lead to a shift from slug flow to either droplet flow or plug flow. The flow patterns in the serpentine microchannel were also depicted. In addition, flow pattern maps of the two-phase systems with Re, Ca and We numbers as the coordinates based on computational work are present. Resuts show that at higher Re, Ca and We numbers of the organic phase, the flow moves towards the parallel flow. At high Re, Ca and We numbers of the aqueous phase and low Re, Ca and We numbers of the organic phase, the flow approaches are droplet flow.

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Ghaini, A., Mescher, A. & Agar, D. W. Hydrodynamic studies of liquid–liquid slug flows in circular microchannels. Chem. Eng. Sci. 66, 1168–1178 (2011).

Article CAS Google Scholar

Yang, C., Zhang, J. & Huang, Z. Numerical study on cavitation–vortex–noise correlation mechanism and dynamic mode decomposition of a hydrofoil. Phys. Fluids 34, 125105 (2022).

Article ADS CAS Google Scholar

Sattari-Najafabadi, M. & Esfahany, M. N. Hexavalent chromium extraction from aqueous solutions in a liquid–liquid slug flow microreactor. Chem. Eng. Process. Process Intensif. 157, 108156 (2020).

Article CAS Google Scholar

Cao, Z., Wu, Z., Sattari Najafabadi, M. & Sunden, B. In Heat Transfer Summer Conference. (2017), V002T010A013 (American Society of Mechanical Engineers).

Zhang, Y. et al. Analysis of grinding mechanics and improved predictive force model based on material-removal and plastic-stacking mechanisms. Int. J. Mach. Tools Manuf 122, 81–97 (2017).

Article Google Scholar

Santana, H., Silva, J., Aghel, B. & Ortega-Casanova, J. Review on microfluidic device applications for fluids separation and water treatment processes. SN Appl. Sci. 2, 1–19 (2020).

Article Google Scholar

Sattari-Najafabadi, M., Sundén, B., Wu, Z. & Nasr Esfahany, M. In International Conference on Nanochannels, Microchannels, and Minichannels. (2016), V001T003A003 (American Society of Mechanical Engineers).

Yu, Y. et al. Targeting loop3 of sclerostin preserves its cardiovascular protective action and promotes bone formation. Nat. Commun. 13, 4241 (2022).

Article ADS CAS PubMed PubMed Central Google Scholar

Wang, L. et al. Therapeutic aptamer targeting sclerostin loop3 for promoting bone formation without increasing cardiovascular risk in osteogenesis imperfecta mice. Theranostics 12, 5645 (2022).

Article CAS PubMed PubMed Central Google Scholar

Liu, Y. et al. A microfluidic biosensor for rapid detection of Salmonella typhimurium based on magnetic separation, enzymatic catalysis and electrochemical impedance analysis. Chin. Chem. Lett. 33, 3156–3160 (2022).

Article CAS Google Scholar

Huang, L. et al. A microfluidic device for accurate detection of hs-cTnI. Chin. Chem. Lett. 32, 1555–1558 (2021).

Article CAS Google Scholar

Sattari-Najafabadi, M., Esfahany, M. N., Wu, Z. & Sundén, B. Hydrodynamics and mass transfer in liquid-liquid non-circular microchannels: Comparison of two aspect ratios and three junction structures. Chem. Eng. J. 322, 328–338 (2017).

Article CAS Google Scholar

Sadeghi, A., Amini, Y., Saidi, M. H. & Yavari, H. Shear-rate-dependent rheology effects on mass transport and surface reactions in biomicrofluidic devices. AIChE J. 61, 1912–1924 (2015).

Article CAS Google Scholar

Sadeghi, A., Amini, Y., Saidi, M. H. & Chakraborty, S. Numerical modeling of surface reaction kinetics in electrokinetically actuated microfluidic devices. Anal. Chim. Acta 838, 64–75 (2014).

Article CAS PubMed Google Scholar

Yang, M. et al. Predictive model for minimum chip thickness and size effect in single diamond grain grinding of zirconia ceramics under different lubricating conditions. Ceram. Int. 45, 14908–14920 (2019).

Article CAS Google Scholar

Liu, S., Niu, B., Zong, G., Zhao, X. & Xu, N. Adaptive neural dynamic-memory event-triggered control of high-order random nonlinear systems with deferred output constraints. IEEE Trans. Autom. Sci. Eng. (2023).

Ban, Y. et al. 2D/3D multimode medical image alignment based on spatial histograms. Appl. Sci. 12, 8261 (2022).

Article CAS Google Scholar

Amini, Y. et al. Optimization of liquid-liquid extraction of calcium with a serpentine microfluidic device. Int. Commun. Heat Mass Transf. 140, 106551 (2023).

Article CAS Google Scholar

Abdollahi, P. & Moosavian, M. A. Microfluidic solvent extraction of calcium: Modeling and optimization of the process variables. Sep. Purif. Technol. 231, 115875 (2020).

Article CAS Google Scholar

Zhao, Y., Niu, B., Zong, G., Xu, N. & Ahmad, A. M. Event-triggered optimal decentralized control for stochastic interconnected nonlinear systems via adaptive dynamic programming. Neurocomputing 539, 126163 (2023).

Article Google Scholar

Yang, B. et al. Motion prediction for beating heart surgery with GRU. Biomed. Signal Process. Control 83, 104641 (2023).

Article Google Scholar

Asadi-Saghandi, H., Karimi-Sabet, J., Ghorbanian, S. & Moosavian, S. M. A. Liquid–liquid extraction of calcium in a scaled-out microfluidic device: Process intensification using a crown ether-ionic liquid system. Chem. Eng. Process. Process Intensif. 183, 109261 (2022).

Article Google Scholar

Cheng, F., Liang, H., Niu, B., Zhao, N. & Zhao, X. Adaptive neural self-triggered bipartite secure control for nonlinear MASs subject to DoS attacks. Inf. Sci. 631, 256–270 (2023).

Article Google Scholar

Gao, Z., Hong, S. & Dang, C. An experimental investigation of subcooled pool boiling on downward-facing surfaces with microchannels. Appl. Therm. Eng. 226, 120283 (2023).

Article CAS Google Scholar

Dessimoz, A.-L., Cavin, L., Renken, A. & Kiwi-Minsker, L. Liquid–liquid two-phase flow patterns and mass transfer characteristics in rectangular glass microreactors. Chem. Eng. Sci. 63, 4035–4044 (2008).

Article CAS Google Scholar

Zhao, Y. et al. A microfluidic surface-enhanced Raman scattering (SERS) sensor for microRNA in extracellular vesicles with nucleic acid-tyramine cascade amplification. Chin. Chem. Lett. 33, 2101–2104 (2022).

Article CAS Google Scholar

Mao, S. et al. A chemo-mechanical switchable valve on microfluidic chip based on a thermally responsive block copolymer. Chin. Chem. Lett. 33, 3083–3086 (2022).

Article CAS Google Scholar

Yao, C., Liu, Y., Xu, C., Zhao, S. & Chen, G. Formation of liquid–liquid slug flow in a microfluidic T-junction: Effects of fluid properties and leakage flow. AIChE J. 64, 346–357 (2018).

Article CAS Google Scholar

Zhang, Q., Liu, H., Zhao, S., Yao, C. & Chen, G. Hydrodynamics and mass transfer characteristics of liquid–liquid slug flow in microchannels: The effects of temperature, fluid properties and channel size. Chem. Eng. J. 358, 794–805 (2019).

Article CAS Google Scholar

Li, P., Yang, M. & Wu, Q. Confidence interval based distributionally robust real-time economic dispatch approach considering wind power accommodation risk. IEEE Trans. Sustain. Energy 12, 58–69 (2020).

Article ADS Google Scholar

Xiang, J. et al. Design and thermal performance of thermal diode based on the asymmetric flow resistance in vapor channel. Int. J. Therm. Sci. 191, 108345 (2023).

Article Google Scholar

Xu, J. H., Li, S., Tan, J. & Luo, G. Correlations of droplet formation in T-junction microfluidic devices: From squeezing to dripping. Microfluid. Nanofluid. 5, 711–717 (2008).

Article CAS Google Scholar

Xu, J., Luo, G., Li, S. & Chen, G. Shear force induced monodisperse droplet formation in a microfluidic device by controlling wetting properties. Lab Chip 6, 131–136 (2006).

Article CAS PubMed Google Scholar

Tsaoulidis, D. & Angeli, P. Effect of channel size on liquid-liquid plug flow in small channels. AIChE J. 62, 315–324 (2016).

Article CAS Google Scholar

Si, Z., Yang, M., Yu, Y. & Ding, T. Photovoltaic power forecast based on satellite images considering effects of solar position. Appl. Energy 302, 117514 (2021).

Article Google Scholar

Cao, Z., Wu, Z. & Sundén, B. Dimensionless analysis on liquid-liquid flow patterns and scaling law on slug hydrodynamics in cross-junction microchannels. Chem. Eng. J. 344, 604–615 (2018).

Article CAS Google Scholar

Moradi, R., Mosavat, M., Gerdroodbary, M. B., Abdollahi, A. & Amini, Y. The influence of coolant jet direction on heat reduction on the nose cone with Aerodome at supersonic flow. Acta Astronaut. 151, 487–493 (2018).

Article ADS CAS Google Scholar

Gerdroodbary, M. B., Amini, Y., Ganji, D. & Takam, M. R. The flow feature of transverse hydrogen jet in presence of micro air jets in supersonic flow. Adv. Space Res. 59, 1330–1340 (2017).

Article ADS Google Scholar

Kashid, M. N. & Agar, D. W. Hydrodynamics of liquid–liquid slug flow capillary microreactor: Flow regimes, slug size and pressure drop. Chem. Eng. J. 131, 1–13 (2007).

Article CAS Google Scholar

Tang, F., Wang, H., Chang, X.-H., Zhang, L. & Alharbi, K. H. Dynamic event-triggered control for discrete-time nonlinear Markov jump systems using policy iteration-based adaptive dynamic programming. Nonlinear Anal. Hybrid Syst 49, 101338 (2023).

Article MathSciNet Google Scholar

Sabernaeemi, A., Barzegar Gerdroodbary, M., Salavatidezfouli, S. & Valipour, P. Influence of stent-induced vessel deformation on hemodynamic feature of bloodstream inside ICA aneurysms. Biomech. Model. Mechanobiol. 22, 1–15 (2023).

Zhao, Y., Chen, G. & Yuan, Q. Liquid-liquid two-phase flow patterns in a rectangular microchannel. AIChE J. 52, 4052–4060 (2006).

Article CAS Google Scholar

Darekar, M., Singh, K. K., Mukhopadhyay, S. & Shenoy, K. T. Liquid–liquid two-phase flow patterns in Y-junction microchannels. Ind. Eng. Chem. Res. 56, 12215–12226 (2017).

Article CAS Google Scholar

Cheng, Y., Niu, B., Zhao, X., Zong, G. & Ahmad, A. M. Event-triggered adaptive decentralised control of interconnected nonlinear systems with Bouc–Wen hysteresis input. Int. J. Syst. Sci. 54(6), 1–14 (2023).

Article MathSciNet CAS Google Scholar

Asadi-Saghandi, H., Karimi-Sabet, J., Ghorbanian, S. & Moosavian, S. M. A. Dimensionless analysis on liquid–liquid two-phase flow patterns in a numbered-up microfluidic device. Chem. Eng. J. 429, 132428 (2022).

Article CAS Google Scholar

Waelchli, S. & von Rohr, P. R. Two-phase flow characteristics in gas–liquid microreactors. Int. J. Multiph. Flow 32, 791–806 (2006).

Article CAS MATH Google Scholar

Yagodnitsyna, A. A., Kovalev, A. V. & Bilsky, A. V. Flow patterns of immiscible liquid-liquid flow in a rectangular microchannel with T-junction. Chem. Eng. J. 303, 547–554 (2016).

Article CAS Google Scholar

Jeong, H.-H., Issadore, D. & Lee, D. Recent developments in scale-up of microfluidic emulsion generation via parallelization. Korean J. Chem. Eng. 33, 1757–1766 (2016).

Article CAS Google Scholar

Schroen, K., Bliznyuk, O., Muijlwijk, K., Sahin, S. & Berton-Carabin, C. C. Microfluidic emulsification devices: From micrometer insights to large-scale food emulsion production. Curr. Opin. Food Sci. 3, 33–40 (2015).

Article Google Scholar

Nisisako, T. & Torii, T. Microfluidic large-scale integration on a chip for mass production of monodisperse droplets and particles. Lab Chip 8, 287–293 (2008).

Article CAS PubMed Google Scholar

Holtze, C. Large-scale droplet production in microfluidic devices—an industrial perspective. J. Phys. D Appl. Phys. 46, 114008 (2013).

Article ADS Google Scholar

Qu, M. et al. Laboratory study and field application of amphiphilic molybdenum disulfide nanosheets for enhanced oil recovery. J. Petrol. Sci. Eng. 208, 109695 (2022).

Article CAS Google Scholar

Hariri, S., Poueinak, M. M., Hassanvand, A., Gerdroodbary, M. B. & Faraji, M. Effects of blood hematocrit on performance of endovascular coiling for treatment of middle cerebral artery (MCA) aneurysms: Computational study. Interdiscip. Neurosurg. 32, 101729 (2023).

Article Google Scholar

Kashid, M. N., Gupta, A., Renken, A. & Kiwi-Minsker, L. Numbering-up and mass transfer studies of liquid–liquid two-phase microstructured reactors. Chem. Eng. J. 158, 233–240 (2010).

Article CAS Google Scholar

Jahromi, P. F., Karimi-Sabet, J., Amini, Y. & Fadaei, H. Pressure-driven liquid-liquid separation in Y-shaped microfluidic junctions. Chem. Eng. J. 328, 1075–1086 (2017).

Article Google Scholar

Amini, Y., Shadman, M. M., Ghazanfari, V. & Hassanvand, A. Enhancement of immiscible fluid mixing using passive micromixers to increase the performance of liquid–liquid extraction. Int. J. Mod. Phys. C 2350149, 1–20 (2023).

Google Scholar

Amini, Y., Mokhtari, M., Haghshenasfard, M. & Gerdroodbary, M. B. Heat transfer of swirling impinging jets ejected from Nozzles with twisted tapes utilizing CFD technique. Case Stud. Therm. Eng. 6, 104–115 (2015).

Article Google Scholar

Ghazanfari, V., Imani, M., Shadman, M. M. & Zahakifar, F. Numerical study on the thermal performance of the shell and tube heat exchanger using twisted tubes and Al2O3 nanoparticles. Prog. Nucl. Energy 155, 104526 (2023).

Article CAS Google Scholar

Jianhua, X. et al. Heat Transfer performance and structural optimization of a novel micro-channel heat sink. Chin. J. Mech. Eng. 35, 1–12 (2022).

Google Scholar

Wu, Z. et al. Computationally driven design of low dielectric-loss epoxy resin for medium-frequency transformers. J. Phys. D Appl. Phys. 56, 184001 (2023).

Article ADS Google Scholar

Li, Z. et al. Analysis of surface pressure pulsation characteristics of centrifugal pump magnetic liquid sealing film. Intern. Flow Mech. Mod. Hydraul. Mach 16648714, 124 (2023).

Google Scholar

Liu, M. et al. Three-dimensional modeling of heart soft tissue motion. Appl. Sci. 13, 2493 (2023).

Article CAS Google Scholar

Sheidani, A. et al. Influence of the coiling porosity on the risk reduction of the cerebral aneurysm rupture: Computational study. Sci. Rep. 12, 19082 (2022).

Article ADS CAS PubMed PubMed Central Google Scholar

Saei Moghaddam, M. & Barzegar Gerdroodbary, M. Analytical study of heat and mass transfer in axisymmetric unsteady flow by ADM. J. Comput. Appl. Res. Mech. Eng. (JCARME) 11, 151–163 (2021).

Google Scholar

Edalatpour, A., Hassanvand, A., Gerdroodbary, M. B., Moradi, R. & Amini, Y. Injection of multi hydrogen jets within cavity flameholder at supersonic flow. Int. J. Hydrog. Energy 44, 13923–13931 (2019).

Article CAS Google Scholar

Amini, Y., Karimi-Sabet, J. & Esfahany, M. N. Experimental and numerical study of multiphase flow in new wire gauze with high capacity structured packing. Chem. Eng. Process. 108, 35–43 (2016).

Article CAS Google Scholar

Amiri, H. A., Shafaghat, R., Alamian, R., Taheri, S. M. & Safdari Shadloo, M. Study of horizontal axis tidal turbine performance and investigation on the optimum fixed pitch angle using CFD: A case study of Iran. Int. J. Numer. Methods Heat Fluid Flow 30, 206–227 (2020).

Article Google Scholar

Hashemipour, N. et al. Experimental and simulation investigation on separation of binary hydrocarbon mixture by thermogravitational column. J. Mol. Liq. 268, 791–806 (2018).

Article CAS Google Scholar

Rahimi, M., Aghel, B. & Alsairafi, A. Experimental and CFD studies on using coil wire insert in a proton exchange membrane fuel cell. Chem. Eng. Process. 49, 689–696 (2010).

Article CAS Google Scholar

Poueinak, M. M. et al. Computational study of blood hemodynamic in ICA aneurysm with coiling embolism. Int. J. Mod. Phys. C 34, 2350138 (2023).

Article Google Scholar

Sadeh, A., Kazemi, A., Bahramkhoo, M. & Barzegar Gerdroodbary, M. Computational study of blood flow inside MCA aneurysm with/without endovascular coiling. Sci. Rep. 13, 4560 (2023).

Article ADS CAS PubMed PubMed Central Google Scholar

Amini, Y. & Nasr-Esfahany, M. CFD simulation of the structured packings: A review. Sep. Sci. Technol. 54, 2536–2554 (2019).

Article CAS Google Scholar

Download references

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Nuclear Fuel Cycle Research School, Nuclear Science and Technology Research Institute, Tehran, Iran

Younes Amini, Valiyollah Ghazanfari, Mehran Heydari, Mohammad Mahdi Shadman, A. Gh. Khamseh & Mohammad Hassan Khani

Department of Polymer Engineering, Faculty of Engineering, Lorestan University, Khorramabad, Iran

Amin Hassanvand

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

Y.A.: wrote the main manuscript text and prepared figures and tabels. Simulation, All authors reviewed the manuscript.V.G.: wrote the main manuscript text and prepared figures and tabels. All authors reviewed the manuscript.M.H.: wrote the main manuscript text and prepared figures and tabels. All authors reviewed the manuscriptM.M.S.: wrote the main manuscript text and prepared figures and tabels. All authors reviewed the manuscriptA.G.K.: wrote the main manuscript text and prepared figures and tabels. All authors reviewed the manuscriptM.H.K.: wrote the main manuscript text and prepared figures and tabels. All authors reviewed the manuscriptA.H.: wrote the main manuscript text and prepared figures and tabels. All authors reviewed the manuscript.

Correspondence to Younes Amini.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Amini, Y., Ghazanfari, V., Heydari, M. et al. Computational fluid dynamics simulation of two-phase flow patterns in a serpentine microfluidic device. Sci Rep 13, 9483 (2023). https://doi.org/10.1038/s41598-023-36672-6

Download citation

Received: 21 April 2023

Accepted: 07 June 2023

Published: 10 June 2023

DOI: https://doi.org/10.1038/s41598-023-36672-6

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.