Next Article in Journal
Numerical Investigation of Special Heat Transfer Phenomenon in Wire-Wrapped Fuel Rod of SFR
Next Article in Special Issue
Mixing Performance of the Modified Tesla Micromixer with Tip Clearance
Previous Article in Journal
High Responsivity Thermopile Sensors Featuring a Mosaic Structure
Previous Article in Special Issue
Microfluidic Applications of Artificial Cilia: Recent Progress, Demonstration, and Future Perspectives
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Hydrodynamic and Thermal Behaviour of Non-Newtonian-Nanofluid Mixing in a Chaotic Micromixer

1
Gas Turbine Joint Research Team, University of Djelfa, Djelfa 17000, Algeria
2
Department of Industrial and Production Engineering, Jashore University of Science and Technology, Jashore 7408, Bangladesh
3
Institute of Nuclear Power Engineering, Bangladesh University of Engineering and Technology, Dhaka 1000, Bangladesh
4
Mechanical Engineering Department, Ziane Achour University of Djelfa, Djelfa 17000, Algeria
5
Department of Mechanical Engineering, Inha University, 100 Inha-ro, Michuhol-gu, Incheon 22212, Korea
*
Authors to whom correspondence should be addressed.
Micromachines 2022, 13(6), 933; https://doi.org/10.3390/mi13060933
Submission received: 27 April 2022 / Revised: 30 May 2022 / Accepted: 2 June 2022 / Published: 11 June 2022
(This article belongs to the Special Issue Analysis, Design and Fabrication of Micromixers, Volume II)

Abstract

:
Three-dimensional numerical investigations of a novel passive micromixer were carried out to analyze the hydrodynamic and thermal behaviors of Nano-Non-Newtonian fluids. Mass and heat transfer characteristics of two heated fluids have been investigated to understand the quantitative and qualitative fluid faction distributions with temperature homogenization. The effect of fluid behavior and different Al2O3 nanoparticles concentrations on the pressure drop and thermal mixing performances were studied for different Reynolds number (from 0.1 to 25). The performance improvement simulation was conducted in intervals of various Nanoparticles concentrations (φ = 0 to 5%) with Power-law index (n) using CFD. The proposed micromixer displayed a mixing energy cost of 50–60 comparable to that achieved for a recent micromixer (2021y) in terms of fluid homogenization. The analysis exhibited that for high nanofluid concentrations, having a strong chaotic flow enhances significantly the hydrodynamic and thermal performances for all Reynolds numbers. The visualization of vortex core region of mass fraction and path lines presents that the proposed design exhibits a rapid thermal mixing rate that tends to 0.99%, and a mass fraction mixing rate of more than 0.93% with very low pressure losses, thus the proposed micromixer can be utilized to enhance homogenization in different Nano-Non-Newtonian mechanism with minimum energy.

1. Introduction

Heat and mass transfer produced by different fluid concentrations is generally applied in high mixing processes that affect the dynamic transport of chemical species within physical chaotic advection. Mixing of non-Newtonian fluids process of micromixers occurs in multiple applications of various devices and industrial applications that have significant utility in bioengineering fields [1], chemical engineering [2], extraction, polymerization techniques and cleaning systems [3,4,5,6,7], etc. To enhance the kinematic behaviour [8] and fluid homogenization performance, many researchers contribute various solutions using chaotic advection inside microchannels [9,10] while reducing pressure drops [11,12,13,14].
Xia et al. [15] investigated experimentally and numerically different passive micromixers with multi-layer crossing channels. Their proposed micromixer can give a fast mixing (more than 95%) at minimum Reynolds number (not exceeding 40), for example, a high mixing index of 0.96 was found for very low flow regimes. The same geometry was selected by Hossain et al. [16] to carry out a parametric study. To develop the mixing efficiency, the effect of the geometrical parameter was analyzed at (not exceeding 40) low Reynolds numbers.
The mixing enhancement of split and recombination micromixers using Newtonian fluid flow with low regimes was investigated numerically by Hossain et al. [17]. They evaluated the homogenization performance and the pressure-drop characteristics compared to the TLCCM micromixer which was selected by Xia et al. [15], in a range of Reynolds numbers between 0.2 and 120. Raza et al. [18,19] developed the same problem to create an unbalanced split and recombination micromixer (1.5 mm). They measured numerically the mixing efficiency within various cross-sections via low flow regime. They found that the optimal configuration has 88% mixing efficiency. A parametric study using five geometric parameters of the micromixer was performed by Naas et al. [8]. They presented the effect of a low Reynolds number on mixing efficiency. Their proposed short micromixer gives mixing rates greater than 99% throughout the Reynolds number range in a short length of three-dimensional micromixer, which is rarely achieved in early micromixer configurations.
The effects of thermal condition of non-Newtonian fluids rheology were investigated by Naas et al. for C-shape chaotic geometry [20,21] and Multi-layer micromixers [22]. They selected a highly efficient heat mixing system and better thermodynamic performance, for rheology index ranging from 1 to 0.49. The results showed that the decrease in heat transfer flow of shear-thinning flows was higher than that of Newtonian fluids.
In order to enhance the heat transfer efficiency, several researchers carried out numerically [23,24,25,26,27,28,29,30] and experimentally [31,32] the thermo-physical enhancement of various applied thermal engineering systems. They used effective nanoparticles such as: Al2O3 [23,24], H2O/SWCNT [29], CuO [28,29,30]. The authors proved experimentally that the nanofluids have significantly more powerful thermal conductivities than the same flow with zero nano concentration.
Inside tubes, Xuan and Li [33] investigated experimentally the thermal performance of nanofluid swirling under wall heat transfer. Compared with pure water, the heat transfer coefficient was augmented more than 39% with the Cu nanoparticles concentration of 2.0%.
Within rectangular microchannels, convection laminar regimes with nano-shear thinning flow have been carried out by Esmaeilnejad et al. [34]. Their obtained results presented that, with Peclet number of 700 and 4% particle concentration, the thermal coefficient will decease around 27.2% and the pressure drop will increase roughly 50.7%. Evangelos et al. [35] studied the mixing performances of heated water with the effect of an electromagnetic field. For different mixers, Pouya [36] investigated numerically, the heat transfer enhancement and mixing quality of a hybrid nanofluid. They found that for high Reynolds numbers and a constant frequency, the mixing rate increased over time.
In other works [37], the authors created a chaotic geometry called C-shape to enhance the heat transfer characteristics of non-Newtonian nanofluid. They found that the development of chaotic advection causes more enhancement to thermal rate in comparison with pumping power.
In recent years, Antar and Kacem [24] analyzed the enhancement of forced convection flow within a straight pipe under the effect of nano-particles of Al2O3subjected to a fixed wall of heat. They examined the effects of both fluid behavior index and nanofluid concentrations on the gradient temperature and heat transfer rates. Their obtained results indicated that, for various Peclet and Brinkman numbers, the heat transfer rate increases when nanoparticles concentration.
Due to the diminutive measurement of micromixers, it’s difficult to enhance the molecular diffusion operation. While the improvement in dynamic behaviors is limited, the injection of nanofluid concentrations is a suitable manner in which to enhance mass transfer and thermal mixing efficiency within a passive micromixer. Therefore, our contribution aims to improve high mixing performances for Nano-non-Newtonian fluids in terms of mass and heat transfer inside the optimum multi-layer micromixer as it was recently used by Naas et al. [8] beneath a very low Reynolds number. Various nanoparticles concentrations with different values of fluid behavior index were proposed to investigate the chaotic flow formation and thermal mixing performances within the suggested micromixer. In order to get important energy efficiency, the homogenization of the fluid indexes and mixing energy cost will be appraised.

2. Micromixer Structure and Problem Statement

A novel micromixer was suggested in this study, namely a Two-Layer Crossing Modified geometry (TLCM) which was used recently and firstly by Naas et al. [8] to achieve High-Mixing Processing under the effect of nanofluid concentrations for power-law non-Newtonian fluid [24], see Figure 1. The micromixer is constructed of a couple of twisted channels; the higher and lower channels are arranged with a periodic chamber. The following mixing components occur in reconstructed forms of several grooves. Dimensions detail is presented in Table 1, where D is the chamber diameter, d is the grooves diameter, dhyd is the hydraulic diameter, l is the distance between the inlets and L* is the length of the geometry. The same nano-non-Newtonian fluid of Antar and Kacem [23,24] and Santra et al. [36] was proposed as working fluid, represented in Table 2 and Table 3 respectively.
Heated and cold Nano-non-Newtonian fluids flow from the inlets with a uniform velocity profile. The outsides are deemed adiabatic and the other boundaries include no-slip. Pressure outlet condition is awarded at the outflow section. A numerical CFD code Fluent© numerically solved the governing equations, which are given by the following expressions [22,38]:
d i v V = 0
where   V is the velocity vector.
V ¯ ¯ V = 1 ρ n f P + d i v τ
where σ (Pa) is shear stress and P is the pressure.
ρ n f c n f V T =   λ n f Δ T
where ρ n f , λ n f and T are the density, the conductivity and the temperature of the nanofluid, respectively.
The constitutive connection among the shear rate γ ˙ (s−1) and shear stress τ (Pa) can be characterized by a simple power-law equation:
τ = m γ ˙ n
where, m (Pas−1) is fluid consistency index and n is the fluid behavior index.
The apparent viscosity is:
μ n f = k γ ˙ n 1
The advised boundary conditions are:
-
Uniform velocity profile imposed to the inlets flow, and the temperatures equal to Tmin = 300 k for one inlet, and the other one Tmax = 330 k
-
No-slip conditions within the solid walls.
-
Pressure outlet condition is considered at the outlet section flow.

2.1. Mass Transfer Characteristics of the Chaotic Flows

Metzner and Reed [39,40] defined the generalized Reynolds number (Reg) for power-law fluids as follows:
Re g = ρ nf u 2 n d hyd n [ 8 n 1 ( b + a n ) n m ]
  • a Geometric parameter ( a = 0.2121)
  • b Geometric parameter ( b = 0.6771)
  • ρ nf Density [ kg m 3 ]
  • u Average speed [ m s 1 ]
  • d hyd Hydraulic diameter [m]
  • m Consistency index (N sn m−2)
  • n Rheological behavior index of the fluid
Mixing efficiency of two Nano-non-Newtonian fluids inside the micromixer is characterized by:
MM i = 1 σ σ 0  
  • σ   is the standard deviation measured by post-CFD.
  • σ 0 is the maximum standard deviation.
σ 2 = 1 N i = 1 N ( C i C ¯ ) 2
For fully mixed fluids, the standard deviation rate was taken from the minimum values and maximum for unmixed fluids. N presents the total number of nodes. C ¯ is the average mass fraction. The maximum standard deviation σ 2 is measured by:
σ 0 2 = ( 1 C ¯ )
The mixing energy cost (MEC) is used to decide the effectiveness of the micromixer and is established by merging the pressure losses and the mixing degree, as follows [41]:
MEC = Δ P × Q MMi
where ΔP and Q are the pressure drop and the flow rate (m3/s) along the geometry, respectively.

2.2. Thermal Characteristics of the Chaotic Flows

Heat transfer coefficient, h , different inlet temperature is assumed as:
h = q ( T b T w )
where, Tb (k) is the mean bulk temperature fluid, Tw (k) is the perimeter average wall temperature and q″ (w/m2) is the wall heat flux.
These two temperatures are defined as:
  T w ( s ) = 1 P P T w dp
T b ( s ) = 1 AU i A V · n T · dA
The mean heat transfer coefficient, h mean , is defined as:
h mean = 1 L 0 L h ( s ) ds
The Thermal Mixing index ( TM i ) utilized by Naas et al. [22] of two fluids (hot and cold) is given by the following equation:
TM i = 1 s | T T avg | ds ( T max T min ) s
where, Tavg = (Tmin + Tmax)/2. T min = 300 K and T max = 330 .
K The values of TM i range from 0 for the unmixed flow case, to 1 for fully mixed flow.
In order to confirm the efficiency of the scalar temperature T which was measured between two values; (Ta, Tb) at a given section, a probability density function PDF % (T) was calculated which was equal to the number of nodes within (Ta, Tb) divided by the total number of nodes [22].
All governing equations in the present work were solved in a laminar flow by using ANSYS Fluent 16© CFD software [42], which works by the (FVM) finite volumes method. The SIMPLEC scheme was chosen for pressure and velocity coupling. To determine the mass and momentum equations a second-order upwind scheme was selected. The computations were ensured and simulated to be converged at 10−7 of root mean square (RMS) residual values. Non-Newtonian power-law fluids were done as a working fluid for various Al2O3 nanoparticles concentrations.

3. Mesh Sensitivity Test

To check the sensitivity of the CFD results, a quantitative grid study was carried outby varying the total number of cells. Using an unstructured mesh with uniform tetrahedral cells, four mesh grids were analyzed ranging from 100,000 to 800,000.
The evolutions of the velocity along the X-axis of the outlet section with a Reynolds number equal to 10 are shown in Figure 2. It can be observed that the outlet velocity rates were sensitive to the mesh, except for the mesh densities with 700,000 and 800,000 nodes where no important difference is seen. As a consequence, the 700,000 nodes grid is selected as the favorable mesh for the analyses (Figure 2 and Figure 3).

4. Results and Discussion

Kinematic and thermal behaviors of mixing Nano-non-Newtonian flows are investigated in detail within novel micromixers, which are compared with potential micromixers used recently in the literature. Heat, mass transfer and fluid mixing process are investigated for several low generalized Reynolds number ranging between 0.1 to 25.

4.1. CDF Validation Case

The numerical solution procedure of mass transfer performances has been reported and validated thoroughly by comparing the present results with the results of Xia et al. [15] and Hossain et al. [16], for Reynolds number equal to 0.2, quantitative and quantitative comparison for mixing efficiency index and mass transfer contours for different successive cross-sectional planes, as shown in Figure 4.
An excellent agreement is seen to exist between the present numerical values and the literature values of mass mixing index. Based on these comparisons, it is perhaps reasonable to conclude that the present results are reliable to within ±0.4%. Deviations of this order are not at all uncommon in numerical studies and arise due to the differences in the flow schematics, problem formulations, grid and/or domain sizes, discretization schemes and numerical methods.
Moreover, a quantitative numerical validation was carried out with those obtained by Ning et al. (2016) [43], and the results present a heat transfer rate as a function of various Reynolds numbers for non-Newtonian cases. The comparison is satisfactory and revealed good agreements among results which are shown in Figure 5.
In addition, we added another validation with those obtained by Jibo et al. [44]. The problem was analyzed for multiphase fluids of mixing enhancement inside micromixers, see Figure 6, where the relative error with the numerical results is less than 1%.

4.2. Mass Transfer and Fluid Mixing Processing

Flow visualization of mass fraction between two nanofluids is shown in Figure 7. Different cases of fluid concentrations are subjected within the fluid pattern to understand the development of visual mixing inside the new configuration. We remark that the fluid mixing augments with the Reg for all fluid behavior index n and there is no exactly black region in the second unites. Therefore, the new micromixers have a quick mass transfer performance.
Figure 8 shows the evaluation of the mass mixing efficiency of two fluids as a function of Reg for various cases of nanofluid concentrations (φ = 0.5 to 5%). As the flow homogenization was only performed by the nanoparticles, the prevalence of Reynolds number is very low (Re < 5) and the flow behaviors were not effective for improving mixing (molecular diffusion dominates). When the generalized Reynolds fluid number increased, the homogenization was more effective, and the mixing intensity developed rapidly. Note that in large Reynolds numbers, the nano-particles are more effective, and the mixing intensity develops quickly so the most select mixing state is reached when the concentration increases. Moreover, it is observed that the proposed micromixer shows a 2.22% enhancement of mixing intensity when the fluid behavior index decreases to 0.46, as compared to the case of n = 0.88.
The homogeneity achieved over midlines at the exit of the several cases of nano-fluid concentrations inside the micromixer at a given Reg is demonstrated by the standard deviation like mass transfer profiles as shown in Figure 5. Significant changes in the values of the mass transfer for all cases of Nano-non-Newtonian fluids can be observed. For the case of φ = 5%, the results indicates that the best degree of mixing done at the exit of the micromixer.
To explain the structure of the fluid flow inside the proposed micromixer, a visualization of vectors and streamlines of the mass fraction and vortex core regions are shown in Figure 9. The kinematic behavior has an important role in enhancing homogenization for nanofluids. As can be seen in Figure 10, the micromixer has a single strong vortex region inside each corner which develops the mixing rate inwardly of the geometry, whereas it has a low-pressure drop near the outlet part. Furthermore, we can see that the flow is more chaotic and dynamic due to the structure and the curve of the configuration. Moreover, it is remarkable that the pathline inside the selected new micromixer produces a reversed flow pattern and strong secondary flows are created which can enhance more mass transfer efficiency and ensure excellent quality of homogenization.

4.3. Heat Transfer and Thermal Mixing Processing

Thermal mixing fluids between hot and cold Nano-non Newtonian fluids are measured for various cases of nanofluid concentrations (φ = 0.5 to 5%). For this, the cold fluid is injected in one inlet at 300K and the heated fluid is injected in the other inlet at 330 K, see Figure 11. In this part, we investigate thermal mixing performances and mining energy coast in the present micromixer, and its performances will be compared to other strong micromixers. Figure 11 shows a top view of the temperature contours in three cases of nanofluid concentration (φ = 0, 2.5 and 5%) for Reg ranging from 0.1 to 25. For a given value of nanofluid concentration or fluid behavior index, the thermal homogenization quality is more vigorous when the Reg is more important. The development of the generalized number highly improves the dynamic of the movement, the kinematic of the fluid nanoparticles varies considerably and the mixing intensity will be improved.
The influence of the nanofluid concentration loading on the improved thermal efficiency is analyzed, see Figure 12, Figure 13 and Figure 14. Figure 12 shows the evolution of heat transfer coefficient for several generalized Reynolds numbers with different nano-fluid concentrations (φ = 0.5 to 5%). As the nanofluid concentration ranges from 0.5% to 5%, the heat transfer augments from 122 W/m2 K to 222 W/m2 K for Reg = 25. It can be resolved that other devices besides thermal conductivity increase can be efficient for thermal flow enhancement. Thermal nanofluid performance is under the effect of various parameters. It appears that nanofluid behavior may develop the temperature gradient within the flow and thus increase the heat of the nanofluid, see Figure 13. Figure 14 shows the evaluation of the thermal mixing efficiency of two fluids for different Reg for various cases of nanofluid concentrations ( φ = 0.5   to   5 % ) .
Due to the heat transfer being done only by conduction, it is indicated that the mixing quality is not affected by the fluid concentrations for low Reg numbers. Therefore, the rheology fluid index n has no effect on the mode of energy transfer in the flow. When Reg augments, the flow is chaotic which leads to greater mixing. This behavior enhances greatly the thermal efficiency as a function of Reg number for all cases of nanofluid concentrations. For Reg numbers equal to 25, the mixing degree is close to 1 and the quality of mixing is perfect. The degree of thermal mixing has an insignificant difference as a function of the fluid concentrations considered, so the fluid behavior index n has the same effect on the TMI as described above.
Table 4 presents the probability density function of temperature distribution for different cross-sections. As it is known, the greatest PDF is 100% for T = 315 K. Despite the fluid passing through the micromixer, the fluids are well mixed and tend to be homogenized (T= 315 K) under the effect of the chaotic behavior of the flow. For all cases of various nanofluid concentrations, the temperature distribution from the third plane (P3) is examined in a small variety of the order of two Kelvin, where the best of this arrangement is compared to the coveted mixing fluid temperature, 315 K.
To choose the better Nano-non-Newtonian fluid concentration with high mixing efficiency and lower cost of mixing energy, Figure 15 presents Mixing Energy Cost under the effect of generalized Reynolds number ranging (Reg = 0.1 to 25) for various nanofluid concentrations (φ = 0.5 to 5%). The flow of low Reg has qualitatively the same behavior in terms of cost of the mixing energy, due to the fluid being more viscous and the secondary flows not yet active. When the Reg overrun the value of 5, the difference enhances obviously for both the kinematic and thermal process. As can be seen, when the Reg increases, the MEC become larger inside the micromixers because the fluid flow turns out more sheared and agitated. Moreover, the fluid flow concentration of 0.5% displays low rates of energy cost compared to the other cases.
The development of Mixing Energy Costis illustrated in Table 5 for various Reynolds numbers. For cases of Reynolds numbers, it can be clearly remarked that the present micromixer has the most economical mixing energy cost compared to the micromixer studied recently by Embarek et al. in [39]. In addition, to compare quantitatively the mixing index, pressure-losses, and the mixing energy cost, Table 6 shows different recent works (2017, 2019, 2021y) for these parameters at fixed Reynolds number (Re = 70). It is shown that the proposed micromixer has the lower pressure drops and the higher mixing index value of 99.99%, and the best mixing energy cost.

5. Conclusions

A numerical CFD work was performed to analyze the flow characteristics of Nano-non-Newtonian shear-thinning fluid inside a novel micromixer. This study outlines the hydrodynamic and thermal mixing performances for different generalized Reynolds numbers with nanofluid concentrations (φ = 0.5 to 5%). According to this study, the following points can be given:
-
Effects of generalized Reynolds numbers on the hydrodynamic behavior of non-Newtonian flow were improved within the proposed micromixers.
-
Strong secondary flows are created inside the micromixer to enhance the mixing quality for all cases of nanofluid concentration.
-
The non-Newtonian fluid of φ = 0.5% exhibits low mass transfer compared to the preferable nanofluid concentration case (φ = 5%).
-
As Reynolds number increases, the flow visualization of both heat and mass transfer revealed that the vortex created in the micromixer had more vigorous intensity.
-
Higher rates of generalized Reynolds number have more effects to increase both mass and thermal homogenization rates.
-
The heat transfer coefficient increases from 122 W/m2K to 222 W/m2K when φ rises from 0.5 to 5%.
-
The cases of low fluid behavior index n have a more effective improvement in the mixing efficiency than the other cases
-
For all of Reg, high intensity thermal and fluid mixing is obtained for high nanofluid concentration (φ = 5%).

Author Contributions

Conceptualization, Numerical Simulation, N.T.T.; Writing Original Draft, S.H.; Conceptualization, Methodology, Writing Original Draft. A.H.K.; Writing—review and editing, K.-Y.K.; Funding, Supervision, T.M.; Numerical Simulation, Writing Original Draft. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The work was partially supported by the UGC research fund 2020–2021. Authors are gratefully acknowledging this support. Authors also acknowledge the support of Jashore University of Science and Technology.

Conflicts of Interest

The authors declare there is no conflict of interest.

References

  1. Lee, C.Y.; Fu, L.M. Recent advances and applications of micromixers. Sens. Actuators B Chem. 2018, 259, 677–702. [Google Scholar] [CrossRef]
  2. Park, T.; Lee, S.; Seong, G.H.; Choo, J.; Lee, E.K.; Kim, Y.S.; Ji, W.H.; Hwang, S.Y.; Gweon, D.G.; Lee, S. Highly sensitive signal detection of duplex dye-labelled DNA oligonucleotides in a PDMS microfluidic chip: Confocal surface-enhanced Raman spectroscopic study. Lab Chip 2005, 5, 437–442. [Google Scholar] [CrossRef] [PubMed]
  3. Rapp, B.E.; Gruhl, F.J.; Länge, K. Biosensors with label-free detection designed for diagnostic applications. Anal. Bioanal. Chem. 2010, 398, 2403–2412. [Google Scholar] [CrossRef] [PubMed]
  4. Rahman, M.; Rebrov, E. Microreactors for gold nanoparticles synthesis: From faraday to flow. Processes 2014, 2, 466–493. [Google Scholar] [CrossRef]
  5. Stone, H.A.; Stroock, A.D.; Ajdari, A. Engineering flows in small devices. Annu. Rev. Fluid Mech. 2004, 36, 381–411. [Google Scholar] [CrossRef]
  6. Squires, T.M.; Quake, S.R. Microfluidics: Fluid physics at the nanoliter scale. Rev. Mod. Phys. 2005, 77, 977–1026. [Google Scholar] [CrossRef]
  7. Hessel, V.; Löwe, H.; Schönfeld, F. Micromixers—A review on passive and active mixing principles. Chem. Eng. Sci. 2005, 60, 2479–2501. [Google Scholar] [CrossRef]
  8. Naas, T.T.; Hossain, S.; Aslam, M.; Rahman, A.; Hoque, A.S.M.; Kim, K.-Y.; Islam, S.M.R. Kinematic measurements of novel chaotic micromixers to enhance mixing performances at low Reynolds numbers: Comparative study. Micromachines 2021, 12, 364. [Google Scholar] [CrossRef]
  9. Raza, W.; Hossain, S.; Kim, K.Y. A Review of Passive Micromixers with a Comparative Analysis. Micromachines 2020, 11, 455. [Google Scholar] [CrossRef]
  10. Ingham, C.J.; van Hylckama Vlieg, J.E.T. MEMS and the microbe. Lab Chip 2008, 8, 1604–1616. [Google Scholar] [CrossRef]
  11. Schulte, T.H.; Bardell, R.L.; Weigl, B.H. Microfluidic technologies in clinical diagnostics. Clin. Chim. Acta 2002, 321, 1–10. [Google Scholar] [CrossRef]
  12. Razzacki, S.Z.; Thwar, P.K.; Yang, M.; Ugaz, V.M.; Burns, M.A. Integrated microsystems for controlled drug delivery. Adv. Drug Deliv. Rev. 2004, 56, 185–198. [Google Scholar] [CrossRef] [PubMed]
  13. Hossain, S.; Lee, I.; Kim, S.M.; Kim, K.Y. A micromixer with two-layer serpentine crossing channels having excellent mixing performance at low Reynolds numbers. Chem. Eng. J. 2017, 327, 268–277. [Google Scholar] [CrossRef]
  14. Afzal, A.; Kim, K.Y. Passive split and recombination micromixer with convergent-divergent walls. Chem. Eng. J. 2012, 203, 182–192. [Google Scholar] [CrossRef]
  15. Xia, H.M.; Wan, S.Y.M.; Shu, C.; Chew, Y.T. Chaotic micromixers using two-layer crossing channels to exhibit fast mixing at low Reynolds numbers. Lab Chip 2005, 5, 748–755. [Google Scholar] [CrossRef] [PubMed]
  16. Hossain, S.; Kim, K.Y. Parametric investigation on mixing in a micromixer with two-layer crossing channels. Springerplus 2016, 5, 794. [Google Scholar] [CrossRef]
  17. Hossain, S.; Kim, K.Y. Mixing analysis in a three-dimensional serpentine split-and-recombine micromixer. Chem. Eng. Res. Des. 2015, 100, 95–103. [Google Scholar] [CrossRef]
  18. Raza, W.; Hossain, S.; Kim, K.Y. Effective mixing in a short serpentine split-and-recombination micromixer. Sens. Actuators B Chem. 2018, 258, 381–392. [Google Scholar] [CrossRef]
  19. Raza, W.; Kim, K.Y. Unbalanced Split and Recombine Micromixer with Three-Dimensional Steps. Ind. Eng. Chem. Res. 2020, 59, 3744–3756. [Google Scholar] [CrossRef]
  20. Naas, T.T.; Lasbet, Y.; Benzaoui, A.; Loubar, K. Characterization of Pressure Drops and Heat Transfer of Non-Newtonian Power-Law Fluid Flow Flowing in Chaotic Geometry. Int. J. Heat Technol. 2016, 34, 251–260. [Google Scholar] [CrossRef]
  21. Naas, T.T.; Lasbet, Y.; Aidaoui, L.; Ahmed, L.B.; Khaled, L. High performance in terms of thermal mixing of non-Newtonian fluids using open chaotic flow: Numerical investigations. Therm. Sci. Eng. Prog. 2020, 16, 100454. [Google Scholar] [CrossRef]
  22. Naas, T.T.; Kouadri, A.; Khelladi, S.; Laib, L. Thermal mixing performances of shear thinning non-Newtonian fluids inside two-layer crossing channels micromixer using entropy generation method: Comparative study. Chem. Eng. Process. 2020, 156, 108096. [Google Scholar]
  23. Antar, T.; Kacem, M. Theoretical investigation of laminar flow convective heat transfer in a circular duct for a non-Newtonian nanofluid. Appl. Therm. Eng. 2017, 112, 1027–1039. [Google Scholar]
  24. Antar, T.; Kacem, M. Analytical solution by Laplace-ritz variational method for non-Newtonian nanofluid inside a circular tube. Int. J. Mech. Sci. 2018, 135, 596–608. [Google Scholar]
  25. Mehryan, S.A.M.; Mohammad, G.; Mohammad, V.; Seyed, M.H.Z.; Nima, S.; Obai, Y.; Ali, J.C.; Hani, A. Non-Newtonian phase change study of nano-enhanced n-octadecane comprising mesoporous silica in a porous medium. Appl. Math. Model. 2021, 97, 463–482. [Google Scholar] [CrossRef]
  26. Nadeem, S.; Rashid, M.; Noreen, S.A. Non-orthogonal stagnation point flow of a nano non-Newtonian fluid towards a stretching surface with heat transfer. Int. J. Heat Mass Transf. 2013, 57, 679–689. [Google Scholar] [CrossRef]
  27. Liang, C.; Nawaz, M.; Hajra, K.; Alaoui, M.K.; Abdellatif, S.; Chuanxi, L.; Hamid, A. Flow and heat transfer analysis of elastoviscoplastic generalized non-Newtonian fluid with hybrid nano structures and dust. Int. Commun. Heat Mass Transf. 2021, 126, 105275. [Google Scholar]
  28. Kourosh, J.; Habib, K.; Abbas, K. Numerical study of heat transfer enhancement of non-Newtonian nanofluid in porous blocks in a channel partially. Powder Technol. 2021, 383, 270–279. [Google Scholar]
  29. Kumar, M.S.; Raju, C.S.K.; El-Sayed, M.S.; Ebrahem, A.A.; Bilal, S.; Harri, J. A Comprehensive Physical Insight about Enhancement in Thermo Physical Features of Newtonian Fluid Flow by Suspending of Metallic oxides of Single Wall Carbon Nano Tube Structures. Surf. Interfaces 2021, 23, 100838. [Google Scholar] [CrossRef]
  30. Amin, S.; Seyed, S.A.; Askari, I.B.; Muhammad, H.A. Numerical investigation of the effect of corrugation profile on the hydrothermal characteristics and entropy generation behavior of laminar forced convection of non-Newtonian water/CMC-CuO nanofluid flow inside a wavy channel. Int. Commun. Heat Mass Transf. 2021, 121, 105–117. [Google Scholar]
  31. Yang, J.C.; Li, F.C.; Zhou, W.W.; He, Y.R.; Jiang, B.C. Experimental investigation on the thermal conductivity and shear viscosity of viscoelastic-fluid-based nanofluids. Int. J. Heat Mass Transf. 2012, 55, 3160–3166. [Google Scholar] [CrossRef]
  32. Li, F.C.; Yang, J.C.; Zhou, W.W.; He, Y.R.; Jiang, B.C. Experimental study on the characteristics of thermal conductivity and shear viscosity of viscoelastic-fluid-based nanofluids containing multi-walled carbon nanotubes. Acta 2013, 556, 47–53. [Google Scholar]
  33. Xuan, Y.; Li, Q. Investigation on convective heat transfer and flow features of nanofluids. J. Heat Transf. 2003, 125, 151–155. [Google Scholar] [CrossRef]
  34. Esmaeilnejad, A.; Aminfar, H.; Neistanak, M.S. Numerical investigation of forced convection heat transfer throughmicrochannels with non-Newtonian nanofluids. Int. J. Therm. Sci. 2014, 75, 76–86. [Google Scholar] [CrossRef]
  35. Evangelos, K.; Christos, L.; Theodoros, K.; Ioannis, S. Mixing of Particles in Micromixers under Different Angles and Velocities of the Incoming Water. Proceeding 2018, 2, 577. [Google Scholar]
  36. Pouya, B. Numerical assessment of heat transfer and mixing quality of a hybrid nanofluid in a microchannel equipped with a dual mixer. Int. J. 2021, 12, 100111. [Google Scholar]
  37. Mehdi, B.; Nima, M.; Masoud, A. Development of chaotic advection in laminar flow of a non-Newtonian nanofluid: A novel application for efficient use of energy. Appl. Therm. Eng. 2017, 124, 1213–1223. [Google Scholar]
  38. Santra, A.K.; Sen, S.; Chakraborty, N. Study of heat transfer due to laminar flow of copper—Water nanofluid through two isothermally heated parallel plates. Int. J. Therm. Sci. 2009, 48, 391–400. [Google Scholar] [CrossRef]
  39. Delplace, F.; Leuliet, J.C. Generalized Reynolds number for the flow of power-law fluids in cylindrical ducts of arbitrary cross-section. Chem. Eng. J. 1995, 56, 33–37. [Google Scholar] [CrossRef]
  40. Metzne, A.B.; Reed, J.C. Flow of non-Newtonian fluids—Correlation of the laminar transition, and turbulent-flow regions. AIChE J. 1955, 1, 434–440. [Google Scholar] [CrossRef]
  41. Embarek, D.; Samir, L.; Amar, K.; Naas, T.T.; Sofiane, K.; Abdelylah, B. High hydrodynamic and thermal mixing performances of efficient chaotic micromixers: A comparative study. Chem. Eng. Process.-Process Intensif. 2021, 164, 108394. [Google Scholar]
  42. ANSYS, Inc. ANSYS Fluent User’s Guide; Release 15, 2013; ANSYS, Inc.: Canonsburg, PA, USA, 2013; pp. 1–814. [Google Scholar]
  43. Li, S.-N.; Zhang, H.-N.; Li, X.-B.; Li, Q.; Li, F.-C.; Qian, S.; Joo, S.W. Numerical study on the heat transfer performance of non-Newtonian fluid flow in a manifold microchannel heat sink. Appl. Therm. 2017, 115, 1213–1225. [Google Scholar] [CrossRef]
  44. Jibo, W.; Guojun, L.; Xinbo, L.; Fang, H.; Xiang, M. A Micromixer with Two-Layer Crossing Microchannels Based on PMMA Bonding Process. Int. J. Chem. React. Eng. 2019, 17, 20180265. [Google Scholar] [CrossRef]
  45. Gidde, R.R. Concave wall-based mixing chambers and convex wall-based constriction channel micromixers. Int. J. Environ. Anal. Chem. 2019, 101, 561–583. [Google Scholar] [CrossRef]
Figure 1. (a): 3-D view design of the proposed micromixers, (b): geometric parameters.
Figure 1. (a): 3-D view design of the proposed micromixers, (b): geometric parameters.
Micromachines 13 00933 g001
Figure 2. Evolutions of outlet velocities at the exit mid-line of X coordinates.
Figure 2. Evolutions of outlet velocities at the exit mid-line of X coordinates.
Micromachines 13 00933 g002
Figure 3. Structure of the generated mesh grid.
Figure 3. Structure of the generated mesh grid.
Micromachines 13 00933 g003
Figure 4. Quantitative validation for mixing efficiency index (a), and quantitative validation of local mass fraction contours for various planes (b) With the results of Xia et al., reproduced with permission from [15] and Hossain et al. reproduced with permission from [16].
Figure 4. Quantitative validation for mixing efficiency index (a), and quantitative validation of local mass fraction contours for various planes (b) With the results of Xia et al., reproduced with permission from [15] and Hossain et al. reproduced with permission from [16].
Micromachines 13 00933 g004
Figure 5. Heat transfer coefficient Vs Reynolds number for non-Newtonian case with Li et al. reproduced with permission from [43].
Figure 5. Heat transfer coefficient Vs Reynolds number for non-Newtonian case with Li et al. reproduced with permission from [43].
Micromachines 13 00933 g005
Figure 6. Evolution of mixing rates Vs Reynolds number with Jibo et al. [44].
Figure 6. Evolution of mixing rates Vs Reynolds number with Jibo et al. [44].
Micromachines 13 00933 g006
Figure 7. Mass fraction contours at various Reynolds numbers with different fluid concentrations (φ = 0.5 to 5%).
Figure 7. Mass fraction contours at various Reynolds numbers with different fluid concentrations (φ = 0.5 to 5%).
Micromachines 13 00933 g007
Figure 8. Development of mass mixing performance for different Reynolds numbers with variation cases of nanofluid concentration (φ = 0.5 to 5%).
Figure 8. Development of mass mixing performance for different Reynolds numbers with variation cases of nanofluid concentration (φ = 0.5 to 5%).
Micromachines 13 00933 g008
Figure 9. Evolution of the mass standard deviation along the (a): exit X-line and (b) exit Y-line of different cases of nano-fluid concentration (φ = 0.5 to 5%) at fixed Reg = 25.
Figure 9. Evolution of the mass standard deviation along the (a): exit X-line and (b) exit Y-line of different cases of nano-fluid concentration (φ = 0.5 to 5%) at fixed Reg = 25.
Micromachines 13 00933 g009
Figure 10. Streamlines and vectors of the mass fraction with vortex core region.
Figure 10. Streamlines and vectors of the mass fraction with vortex core region.
Micromachines 13 00933 g010
Figure 11. Qualitative representation of Temperature contours for different generalized Reynolds number at the horizontal middle section of each nano-fluid concentrations (φ = 0.5 to 5%).
Figure 11. Qualitative representation of Temperature contours for different generalized Reynolds number at the horizontal middle section of each nano-fluid concentrations (φ = 0.5 to 5%).
Micromachines 13 00933 g011
Figure 12. Evolutions of heat transfer coefficient for different generalized Reynolds numbers with various nanofluid concentrations (φ = 0.5 to 5%).
Figure 12. Evolutions of heat transfer coefficient for different generalized Reynolds numbers with various nanofluid concentrations (φ = 0.5 to 5%).
Micromachines 13 00933 g012
Figure 13. Evolutions of local heat transfer coefficient as a function along the geometry with φ = 1 and 5%.
Figure 13. Evolutions of local heat transfer coefficient as a function along the geometry with φ = 1 and 5%.
Micromachines 13 00933 g013
Figure 14. Improvement of thermal mixing performance for different generalized Reynolds numbers for variation cases of nanofluid concentrations (φ = 0.5 to 5%).
Figure 14. Improvement of thermal mixing performance for different generalized Reynolds numbers for variation cases of nanofluid concentrations (φ = 0.5 to 5%).
Micromachines 13 00933 g014
Figure 15. Development of Mixing Energy Cost for several generalized Reynolds numbers (Reg = 0.1 to 25) for various nanofluid concentrations (φ = 0.5 to 5%).
Figure 15. Development of Mixing Energy Cost for several generalized Reynolds numbers (Reg = 0.1 to 25) for various nanofluid concentrations (φ = 0.5 to 5%).
Micromachines 13 00933 g015
Table 1. Dimensions information of the proposed micromixer.
Table 1. Dimensions information of the proposed micromixer.
W0.2 mm
l*0.8 mm
D*0.2 mm
d*0.1 mm
dhyd0.22 mm
D0.6 mm
L4.5 mm
Table 2. Rheological parameters of Al2O3 Nano-non-Newtonian, reproduced with permission from [23,24].
Table 2. Rheological parameters of Al2O3 Nano-non-Newtonian, reproduced with permission from [23,24].
φ % m (N sn m−2)N
0.50.001870.88
1.00.002300.83
1.50.002830.78
2.00.003470.730
2.50.004260.680
3.00.005350.625
3.50.006410.580
4.00.007500.540
4.50.008760.500
5.00.010200.460
Table 3. Thermal parameters of Al2O3 Nano-non-Newtonian, reproduced with permission from [23,24].
Table 3. Thermal parameters of Al2O3 Nano-non-Newtonian, reproduced with permission from [23,24].
( φ   % ) ρ (Kg/m3) C p   ( J / Kg   k ) k (w/m k)
0.51013.14165.90.6248
1.01027.94148.80.6367
1.51042.84131.70.6488
2.01057.64114.60.6610
2.51072.54097.50.6734
3.01087.44080.50.6859
3.51102.24063.40.6987
4.01117.14046.30.7116
4.51131.94029.20.7246
5.01146.84012.10.7379
Table 4. Local Temperature PDF ranging from 314 K to 316 K for different cross-sections.
Table 4. Local Temperature PDF ranging from 314 K to 316 K for different cross-sections.
φ (%) P1 = 0.93 mmP2 = 1.82 mmP3 = 2.62 mmP4 = 3.47 mmOutlet
0.55%31.45%99.99%99.99%99.99%
27.83%26.23%99.99%99.99%99.99%
47.85%24.43%99.99%99.99%99.99%
512.24%30.49%99.99%99.99%99.99%
Table 5. Quantitative comparison of the Energy Mixing Cost for varied generalized Reynolds number (Re = 1, 5, 15 and 15) with that obtained by Embarek et al. reprinted with permission from [41] in 2021y.
Table 5. Quantitative comparison of the Energy Mixing Cost for varied generalized Reynolds number (Re = 1, 5, 15 and 15) with that obtained by Embarek et al. reprinted with permission from [41] in 2021y.
ReMEC (μW) of Embarek [39]MEC (μW) of Present WorkEvolution Parentage
10.009200.006430%
50.236000.19418%
152.573021.54840%
3012.997.4345%
Table 6. Comparison of the pressure losses, mixing degree and energy mixing cost with different recent works at fixed Reynolds number (Re = 70).
Table 6. Comparison of the pressure losses, mixing degree and energy mixing cost with different recent works at fixed Reynolds number (Re = 70).
Micromixer, YearUnitLength (mm)ΔP (Pa)MIMEC (μW)
Split and recombination micromixer, 2012 [14]610.66165000.83210.25
Convergent and divergent micromixer, 2017 [41]611140000.94209.55
rectangular obstacles micromixer 2019 [45]58.0125420.98179.75
triangular obstacles micromixer 2019 [45]58.0114850.93175.38
teardrop obstacles micromixer 2019 [45]58.0109080.94163.22
Elongation micromixer with 2021 [41]43.580020.99111.48
Present Two-layer Modified micromixer 54.5834550.9955.67
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tayeb, N.T.; Hossain, S.; Khan, A.H.; Mostefa, T.; Kim, K.-Y. Evaluation of Hydrodynamic and Thermal Behaviour of Non-Newtonian-Nanofluid Mixing in a Chaotic Micromixer. Micromachines 2022, 13, 933. https://doi.org/10.3390/mi13060933

AMA Style

Tayeb NT, Hossain S, Khan AH, Mostefa T, Kim K-Y. Evaluation of Hydrodynamic and Thermal Behaviour of Non-Newtonian-Nanofluid Mixing in a Chaotic Micromixer. Micromachines. 2022; 13(6):933. https://doi.org/10.3390/mi13060933

Chicago/Turabian Style

Tayeb, Naas Toufik, Shakhawat Hossain, Abid Hossain Khan, Telha Mostefa, and Kwang-Yong Kim. 2022. "Evaluation of Hydrodynamic and Thermal Behaviour of Non-Newtonian-Nanofluid Mixing in a Chaotic Micromixer" Micromachines 13, no. 6: 933. https://doi.org/10.3390/mi13060933

APA Style

Tayeb, N. T., Hossain, S., Khan, A. H., Mostefa, T., & Kim, K. -Y. (2022). Evaluation of Hydrodynamic and Thermal Behaviour of Non-Newtonian-Nanofluid Mixing in a Chaotic Micromixer. Micromachines, 13(6), 933. https://doi.org/10.3390/mi13060933

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