Next Article in Journal
Transient Electrophoresis of a Cylindrical Colloidal Particle
Next Article in Special Issue
Modeling and Experimental Characterization of a Clutch Control Strategy Using a Magnetorheological Fluid
Previous Article in Journal
Experimental and CFD Investigation of Directional Stability of a Box-Wing Aircraft Concept
Previous Article in Special Issue
Magnetohydrodynamics Solver for a Two-Phase Free Surface Flow Developed in OpenFOAM
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Simulation of Mixing Fluid with Ferrofluid in a Magnetic Field Using the Meshless SPH Method

1
Faculty of Mechanical Engineering, University of Ljubljana, Aškerčeva cesta 6, 1000 Ljubljana, Slovenia
2
Department of Mechanical Engineering, Yasouj University, Yasouj 75918-74831, Iran
3
Institute of Metals and Technology, Lepi pot 11, 1000 Ljubljana, Slovenia
*
Author to whom correspondence should be addressed.
Fluids 2022, 7(11), 341; https://doi.org/10.3390/fluids7110341
Submission received: 16 September 2022 / Revised: 18 October 2022 / Accepted: 26 October 2022 / Published: 29 October 2022
(This article belongs to the Special Issue The Recent Advances in Magnetorheological Fluids)

Abstract

:
In this study, a numerical investigation of the effect of different magnetic fields on ferrofluid-fluid mixing processes in a two-dimensional microchannel is performed An improved version of smoothed particle hydrodynamics, SPH, by shifting particle algorithm and dummy particle boundary condition, is implemented to solve numerical continuity, ferrohydrodynamics-based momentum and mass transfer equations. SPH is formulated through the irregular arrangement of the nodes where the fields are approximated using the fifth-order Wendland kernel function. After validating the computational approach, the influence of the number (from one to three) of parallel electrical wires positioned perpendicular to the microchannel on the mixing efficiency is studied for the first time. It has originally been found that the mixing efficiency highly non-linearly depends on the Reynolds number and the number of electrical wires. For Re ≤ 20 the mixing efficiency is almost the same for two and three electrical wires and about two times higher than one electrical wire. For Re ≥ 80, the mixing efficiency of three wires is much higher than one and two electrical wires. Optimum performance of the micromixer is achieved with three electrical wires, since the mixer performs well on a broader range of Re than the other two studied cases. The outcomes of this study, obtained by a meshless method, are important for the industrial design of micromixers.

1. Introduction

The study of the mixing process of fluids in the submillimeter scales represents an interesting research subject due to the design of different microfluidic systems. Micromixing is the process of mixing at the smallest scale of fluid motion (i.e., molecular scale). Micromixing has various applications such as in microfluidic lab-on-a-chip devices, drug delivery, bioengineering, etc. [1]. There are two general types of micromixers: passive and active [2]. Passive micromixers do not need any external force for flow disturbance. In this case, the mixing process is governed primarily by diffusion and convection due to the flow’s kinetic energy. Active micromixers, on the other hand, require an external force for disturbance. Although active micromixers, compared to passive ones, are often more complex and expensive, they are essential in some of applications [3]. There are different external actuators to improve the mixing process, including acoustic field [4,5,6], thermal field [7,8], pressure-driven field [9,10,11,12,13,14,15], electric field [16,17,18,19], and magnetic field [20]. In the present work, the effect of the magnetic field on the mixing process in a channel flow will be explored. In the following, literature on the micromixing process, focused on the magnetic field, is reviewed.
Unlike the macroscale fluidic devices, where the mixing process mainly depends on the convection mechanism, mixing at the microscale mostly depends on diffusion. Moreover, due to the microscale system’s small dimensions and the fluid’s low velocity, the Reynolds number is small, and thus, the flow regime is laminar [21].
Andersson et al. [22] have numerically explored the effect of the magnetic field on a viscous ferrofluid and concluded that the presence of the magnetic field increases the skin friction. The biomagnetic fluid flow over a stretching sheet, considering a non-linear temperature-dependent magnetization, is studied in Ref. [23]. It is often assumed that the magnetic field strength is adequately strong to make the biofluid saturated, and also, the magnetization can be stated as a function of both the magnetic field intensity and temperature [24,25]. Tzirtzilakis [26] proposed a mathematical model for biomagnetic fluid dynamics to describe a laminar incompressible Newtonian blood flow under the influence of a magnetic field. He also studied the effects of the strength and the gradient of the magnetic field on the flow field properties. Tsai et al. [27] experimentally showed that by a well-positioned magnet, the mixing performance of a ferro-nanofluid and water can be improved. The mixing mechanisms of a rapid microfluidic mixer actuated by direct or alternating currents of electricity (DC or AC) have been investigated in Ref. [28]. A numerical simulation of an active mixing system has been examined in Ref. [29] by applying a hybrid gradient of a magnetic field composed of a combination of a static gradient magnetic field and an external uniform magnetic field of alternating current. The effects of the magnetic field on fluid flow and heat transfer in an aneurysm were investigated numerically by Sharifi et al. [30].
The mixing of Newtonian and non-Newtonian fluids in a three-dimensional steady-state micromixer, actuated with a magnetic field [31] and in a spiral passive micromixer [32], has been explored numerically. Lab-on-chip devices, one of the micromixers applications, were investigated by studying the effects of steady and periodic magnetic fields on their mixing performance [33]. Furthermore, magnetic drug targeting, considering the impact of a magnetic field, was introduced (e.g., see Refs. [34,35]). Many other micromixers, actuated by a magnetic field, have been presented in the literature. See reviews on micromixers with magnetic actuators [36,37,38] for more details.
The meshless methods, particularly the particle-based category, have recently received increased attention from researchers studying the mixing process [39,40,41,42]. Various types of meshless approaches, such as Moving Least Square (MLS) [43,44], particle methods [45], local radial basis function collocation method [46,47], and many others, can be found in Ref. [48]. Smoothed particle hydrodynamics (SPH), a fully Lagrangian particle-based method, was first conducted to model astronomical phenomena [49,50]. Since then, this method has been developed and applied for different fluid dynamics situations. In particular, to properly simulate problems such as free surface flows [51,52,53], multiphase flows [54,55,56], and other pertinent subjects. Due to the various applications of the mixing process and its importance, the SPH method is utilized here to study the micromixing performance under the influence of the magnetic field. Therefore, the main goal of the current work is to numerically investigate the mixing process of a laminar Newtonian viscous fluid and a ferrofluid in the presence of a differently configured magnetic field. By considering a simple relation between the magnetization M and the magnetic field strength H, the influence of the number of different conducting wires, positioned perpendicular to the microchannel on the mixing index is investigated.
The remainder of the paper is organized as follows; Section 2 describes the geometry of the micromixers. Section 3 presents the governing equations. Section 4 contains the numerical procedure. In Section 5, the results are discussed; finally, in Section 6, the conclusions are presented.

2. Physical Models Description

A scheme of the considered micromixer is depicted in Figure 1. The microchannel depth is assumed to be large enough to ignore variations of the flow field properties perpendicular to the x-y plane. Several experimental studies proved the possibility of studying the mixing process in two dimensions [57,58]. Therefore, considering 2D models of micromixers provides a reasonably good approximation to some of their designs [59]. Dimensions of microchannel here are specified based on a unit length, U L . As Figure 1 shows, the two fluids enter the microchannel with the same velocity, U i n . The miscible fluids have the same viscosity and density, whereas their concentrations are C o = 0.0 and C 1 = 1.0 , at the entrance of the microchannel. The three electrical wires generating the magnetic field are located at fixed positions, as depicted in (Figure 1).

3. Governing Equations

In this section, the conservation laws, i.e., the continuity, momentum, and mass transport equations, as the governing equations to study the effects of different magnetic fields on the mixing process of ferrofluid, using a mixing index defined subsequently, are presented as:
1 ρ ρ t + u i x i = 0 ,
u i t + u i u j x i = 1 ρ p x i + ν 2 u j x i x i + M n H H x i ,
C t + u i C x i = α c 2 C x i x i ,
where t, x i , ρ , u i , p, ν , M n , H, C, and α c are the time, direction of the 2D Cartesian system (here are X and Y), the density of the fluid, velocity components, pressure, kinematic viscosity, dimensionless magnetic number, the magnitude of the magnetic field, concentration, and mass diffusivity, respectively. The last term on the right-hand side of Equation (2) is the magnetic force. This form of the magnetic term can be used when the assumption of M H is met [26]. The relation between the magnetization M and the magnetic field intensity H is stated as M = χ H , where M and χ are the magnetization and the magnetic susceptibility, respectively. Considering the Biot–Savart law, the components of the magnetic field H in the directions of x and y can be written as [25]:
H x ( x , y ) = I 2 π 1 ( x a ) 2 + ( y b ) 2 ( y b ) ,
H y ( x , y ) = I 2 π 1 ( x a ) 2 + ( y b ) 2 ( x a ) ,
where a and b are the x and y positions of the wires, respectively. Thus, the magnetic field intensity as a function of the electric current I can be determined as:
H ( x , y ) = H x 2 + H y 2 = I 2 π 1 ( x a ) 2 + ( y b ) 2 .
The dimensionless magnetic number M n is defined as:
M n = D 2 μ 0 χ H o 2 ν 2 ρ = D 2 B 0 M 0 ν 2 ρ ,
where μ 0 = 4 π × 10 7 N / A 2 is the magnetic permeability of the vacuum and D is the width of the channel shown in Figure 1. In addition, H 0 , B 0 = μ 0 H o and M 0 = χ H o are respectively the magnetic strength, the magnetic field induction, and the magnetization in the minimum distance between the lower wall of the microchannel and the location of the wires.
To close the listed governing equations, an equation of state is required:
P = c 0 2 ( ρ ρ 0 ) ,
in which c 0 and ρ 0 are the speed of sound within the fluid and initial density, respectively. The speed of sound is assumed to be much higher than the maximum fluid velocity. This assumption is made to decrease the fluctuations of the density. The no-slip boundary condition for velocity and the Neumann boundary condition for pressure and concentration at walls are considered.
To quantify the mixing performance of micromixers, a mixing index, based on the standard deviation of concentration, can be defined as [3]:
M I = 1 N j = 1 N ( C j C m e a n C m e a n ) 2 ,
where N, C j , and C m e a n are, respectively, the total number of particles, the concentration of particle j, and the mean concentration (here C m e a n = 0.5 ). The mixing index can be evaluated from 0, unmixed state, to 1, full mixed state. Equation (15) can be used to represent the mixing efficiency as [3]:
η m i x = ( 1 M I ) × 100 .
Higher values of the mixing efficiency η m i x indicate greater homogeneity of the mixed fluids and thus higher performance of the micromixer.

4. Numerical Procedure

As mentioned before, SPH is one of the frequently used Lagrangian meshless methods for modeling complex phenomena in physics. The method is based on particle theory, in which continuous partial differential equations are approximated by algebraic equations. A comprehensive description of the theory and applications of SPH can be found in Ref. [60].
The function of interpolation for a scalar field A, the pressure and concentration, or a vector field A , the velocity and magnetic field, of particle a, which is located at position r a , in the SPH framework are defined respectively as [61]:
A a b V b A b w a b ,
A a b V b A b w a b ,
in which V b is the volume of particle b, w a b = w h ( r a r b ) is the function of interpolation and h is the smoothing length. The Wendland kernel function of 5th order is considered in the present work [62]:
w h ( q ) = α w h 2 ( 1 q 2 ) 4 ( 2 q + 1 ) 0 q 2 0 2 q ,
where α w = 7 4 π is a constant, q = r a r b h , ra and rb are the position vectors of particles a and b, respectively.
The operators such as gradient (grad), divergence (div), and laplacian (div(grad)) for a scalar (A) and a tensor (A) are expressed as follows [61,63]:
{ ( grad A ) a } b V b A b a B a . w a b ´ ,
{ ( div A ) a } b V b A b a ( B a . w a b ´ ) ,
div { grad { A a } } b 2 V b A b a r a b e ab . ( B a . w a b ´ ) ,
in which e ab , w a b ´ , and B a are respectively the unit vector between particles a and b, the derivative of the kernel function, and the factor of renormalization, which is determined as [61]:
B a = [ b V b ( r a b ) w a b ´ ] 1 .
Implementation of boundary conditions represents a challenge in the SPH method. Various techniques have been developed to treat the boundary conditions properly [62,63,64]. In this work, the Dummy particle method [65] is used for wall boundaries, and the open channel flow conditions are considered for inflow and outflow boundaries [64]. For further information on these techniques, see the work of Lee et al. [65]. Different remedies were also proposed for solving the non-uniformity distribution in SPH [66,67]. We have used the particle shifting algorithm for this purpose [68]. Furthermore, the prediction-correction scheme is used for the time stepping algorithm, in which a variable time-step size is utilized as [69]:
Δ t = ε t min ( δ min c + V max , δ min 2 ν , δ min g , δ min 2 α c ) ,
where ϵ t is a constant ( 0 < ε t < 1 ), δ min is the minimum distance between two particles, g is gravity acceleration and V max is the maximum particle velocity.

5. Results and Discussions

5.1. Validation

In order to verify and check the accuracy of the computational approach, the problem of the mixing process of two fluids in a rectangular microchannel is solved and the results are compared to those of Kim et al. [70], which were obtained by the Lattice Boltzmann method (LBM). LBM, which originated in lattice gas automata, is another class of computational fluid mechanics that uses streaming and collision (relaxation) processes to simulate a fluid density on a lattice rather than directly implementing the Navier–Stokes equations. In contrast, SPH is a meshless method that directly solves the Navier–Stokes equations to simulate fluid flow in a Lagrangian framework. The Reynolds number R e = 240 and the Schmidt number S c = 10 , defined as the ratio of the kinematic viscosity to the mass diffusivity, are the characteristic dimensionless numbers of this problem. A total of 14,000 computational particles are used here to simulate the problem. Figure 2 shows the variations of the concentration versus the non-dimensional characteristic length of the channel. A good agreement between the LBM and SPH results is observed.

5.2. Mixing Process under the Effect of Magnetic Fields

The results of the numerical simulation, using the SPH method, for the mixing process of ferrofluid under the influence of different magnetic fields follow. This arrangement contributes to the forced mixing due to the external magnetic field force. For the problem under consideration, the magnetic number M n = 1000 , saturation magnetization M s = 60   A / m , magnetization at the location of wires M o = 60 , Schmidt number S c = 1 and different Reynolds numbers, 1 R e 100 , are considered. Furthermore, the effects of different numbers of electrical wires with various Reynolds numbers are explored; to this end, the contours of velocity, concentration, and particle distribution for three electrical wires at R e = 15 are drawn in Figure 3 and Figure 4. It is worth noting that the results of R e = 15 are presented here as a sample of the results of the other Reynolds numbers. Figure 3 suggests that by increasing the number of wires, the production of the vortices also increases. The vortex shedding phenomenon leads to an increase in the interface between two fluids and ferrofluid, and subsequently, the mixing efficiency is enhanced. The main reason for this enhancement is due to the increased contact time of the fluids due to the vortex shedding. Figure 4 indicates that an increase in the number of wires leads to a better homogeneity for the fluid in a specific desired time. It is worth mentioning that the results are presented at t * = 0.014 , where t * = t U i n / U L is the dimensionless time. Moreover, considering both the contours of concentration and particle distribution, the particle dispersion at the outlet of the channel is improved, which favorably influences the mixing index. The mixing efficiencies at Re = 15 for one, two, and three wires are 61.85%, 72.22%, and 77.16%, respectively, confirming this fact.
Figure 5 shows a comprehensive view of the effects of the number of electrical wires on the mixing efficiency in a wide range of Reynolds numbers. This figure can be approximately divided into some small intervals; for the interval of 0 R e 20 , the mixing efficiencies for the case of two and three wires are better than in the case of one wire; moreover, the three-wires case has higher efficiencies, compared to the case of two wires. This trend can be interpreted as the vortex shedding intensification due to an increase in the number of wires so that the two fluid elements are more in contact and have more time to become homogenized. Interestingly, there is an intersection point around R e = 20 for all three cases. In the interval of 20 R e 52 , however, the effect of the number of electrical wires is opposed to that of the beginning interval; furthermore, the behavior of the one wire case is almost parabolic, and surprisingly this case has better efficiencies, compared to the other cases. In addition, it is indicated that the performance of the three-wires case is the worst, in this interval. As the Reynolds number exceeds 55, the mixing efficiency of the one wire case significantly decreases. In contrast, this behaviour for the case of two wires is observed at the Reynolds number of about 70. The three-wires case performance slightly increases in the interval of 60 R e 100 . For Reynolds numbers higher than 80, the mixing performance of the three-wires case increases compared to the two other cases. The overall result indicates that the number of electrical wires positively affects the performance of the micromixing process. However, it is fair to say that the mixing efficiency for the case of three wires has a more stable behavior (fewer ups and downs) for different Reynolds numbers. If one tends to consider a micromixing device working in a wide range of Reynolds numbers, it is more suitable to use the three electrical wire case.
Figure 6 shows the variation of the mixing index with the dimensionless length of the microchannel for different Reynolds numbers; from which it can be concluded that for the cases of R e = 15 , 30 , 45 , the effect of the number of electrical wires at the end of the microchannel is less significant, compared to the other Reynolds numbers. For the case of R e = 1 , the influence of increasing the number of electrical wires on the mixing performance is more effective, in comparison with the other Reynolds numbers, at the end of the microchannel (Figure 6). The mixing index values at R e = 100 for the cases of one and two wires are almost identical, except at the end of the channel, whereas the three-wires case clearly indicates a significant improvement in the mixing process (Figure 6).
Generally speaking, for a wide range of Reynolds numbers, the case of a micromixer with three wires is optimal, among all the examined cases. The three-wires case maintains the behavior of mixing efficiency at a high level for different Reynolds numbers.

6. Conclusions

In the present work, the mixing process of ferrofluid in a microchannel was studied numerically using the meshless SPH method. The effects of increasing the number of magnetic fields on the mixing index for a wide range of applicable Reynolds numbers have been explored. After validating the computational approach, the results for the problem under consideration were presented. The most important findings of the current investigation can be summarized as:
  • By increasing the number of electrical wires, for some Reynolds numbers, the vortex shedding intensifies and causes an increase in the interface between the two fluids, which leads to a better mixing efficiency;
  • The homogeneity of the mixing is improved by increasing the number of wires;
  • The micromixer with three wires has optimum performance for a range of Reynolds numbers compared to the cases of one and two wires;
  • For R e = 15 , 30 , 45 , the variation of efficiency for different numbers of wires is not as significant as that of the other Reynolds numbers at the outlet of the channel;
  • If one aims to design a micromixer working in a wide range of Reynolds numbers, the results show that the case with three magnetic wires is recommended.
The findings of the represented study, performed by a meshless method, are important for the industrial design of micromixers. The future extension of this study will focus on mixing non-Newtonian fluids in domains with geometrically moving boundaries and channels with passive mixers, for which SPH is the most suitable approach. This research also represents the power of meshless methods in industrial applications.

Author Contributions

Conceptualization, M.A. (Mohsen Abdolahzadeh) and A.T.; validation, M.A.(Mohsen Abdolahzadeh); formal analysis, M.A. (Mohsen Abdolahzadeh) and A.T.; resources, M.A.(Mohsen Abdolahzadeh); data curation, M.A.(Mohsen Abdolahzadeh); writing—original draft preparation, M.A.(Mohsen Abdolahzadeh), A.T., and M.A.(Mehrdad Ahmadinejad); writing—review and editing, A.T. and B.Š.; supervision, B.Š.; project administration, B.Š.; funding acquisition, B.Š. All authors have read and agreed to the published version of the manuscript.

Funding

This research was made possible by support from Yasouj University and University of Ljubljana under core funding P2-0162 and project J2-4477.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The first author is grateful for the support from the Slovenian Grant Agency ARRS in the framework of the Young Researcher Program.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jeong, G.S.; Chung, S.; Kim, C.B.; Lee, S.H. Applications of micromixing technology. Analyst 2010, 135, 460–473. [Google Scholar] [CrossRef] [PubMed]
  2. 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]
  3. Nguyen, N.T. Micromixers: Fundamentals, Design and Fabrication; William Andrew: Norwich, NY, USA, 2011. [Google Scholar]
  4. Lin, Y.; Gao, C.; Gao, Y.; Wu, M.; Yazdi, A.A.; Xu, J. Acoustofluidic micromixer on lab-on-a-foil devices. Sens. Actuators B Chem. 2019, 287, 312–319. [Google Scholar] [CrossRef]
  5. Van Phan, H.; Coşkun, M.B.; Şeşen, M.; Pandraud, G.; Neild, A.; Alan, T. Vibrating membrane with discontinuities for rapid and efficient microfluidic mixing. Lab A Chip 2015, 15, 4206–4216. [Google Scholar] [CrossRef]
  6. Ang, K.M.; Yeo, L.Y.; Hung, Y.M.; Tan, M.K. Amplitude modulation schemes for enhancing acoustically-driven microcentrifugation and micromixing. Biomicrofluidics 2016, 10, 054106. [Google Scholar] [CrossRef] [Green Version]
  7. Tan, H. Numerical study of a bubble driven micromixer based on thermal inkjet technology. Phys. Fluids 2019, 31, 062006. [Google Scholar] [CrossRef]
  8. Huang, K.R.; Chang, J.S.; Chao, S.D.; Wung, T.S.; Wu, K.C. Study of active micromixer driven by electrothermal force. Jpn. J. Appl. Phys. 2012, 51, 047002. [Google Scholar] [CrossRef]
  9. Niu, X.; Lee, Y.K. Efficient spatial-temporal chaotic mixing in microchannels. J. Micromech. Microeng. 2003, 13, 454. [Google Scholar] [CrossRef] [Green Version]
  10. Du, M.; Ma, Z.; Ye, X.; Zhou, Z. On-chip fast mixing by a rotary peristaltic micropump with a single structural layer. Sci. China Technol. Sci. 2013, 56, 1047–1054. [Google Scholar] [CrossRef]
  11. Wang, X.; Ma, X.; An, L.; Kong, X.; Xu, Z.; Wang, J. A pneumatic micromixer facilitating fluid mixing at a wide range flow rate for the preparation of quantum dots. Sci. China Chem. 2013, 56, 799–805. [Google Scholar] [CrossRef]
  12. Abbas, Y.; Miwa, J.; Zengerle, R.; Von Stetten, F. Active continuous-flow micromixer using an external braille pin actuator array. Micromachines 2013, 4, 80–89. [Google Scholar] [CrossRef] [Green Version]
  13. Xia, Q.; Zhong, S. Liquid mixing enhanced by pulse width modulation in a Y-shaped jet configuration. Fluid Dyn. Res. 2013, 45, 025504. [Google Scholar] [CrossRef]
  14. Abdolahzadeh, M.; Tayebi, A.; Omidvar, P. Mixing process of two-phase non-newtonian fluids in 2D using smoothed particle hydrodynamics. Comput. Math. Appl. 2019, 78, 110–122. [Google Scholar] [CrossRef]
  15. Yao, L.B.; Wu, W.; Wu, X.S.; Chu, G.W.; Luo, Y.; Sun, B.C. Intensification of micromixing efficiency in a spinning disk reactor: Experimental investigation. Chem. Eng. Process. Process Intensif. 2021, 166, 108500. [Google Scholar] [CrossRef]
  16. Usefian, A.; Bayareh, M.; Shateri, A.; Taheri, N. Numerical study of electro-osmotic micro-mixing of Newtonian and non-Newtonian fluids. J. Braz. Soc. Mech. Sci. Eng. 2019, 41, 238. [Google Scholar] [CrossRef]
  17. Usefian, A.; Bayareh, M. Numerical and experimental study on mixing performance of a novel electro-osmotic micro-mixer. Meccanica 2019, 54, 1149–1162. [Google Scholar] [CrossRef]
  18. Kumar, V.; Paraschivoiu, M.; Nigam, K. Single-phase fluid flow and mixing in microchannels. Chem. Eng. Sci. 2011, 66, 1329–1373. [Google Scholar] [CrossRef]
  19. Zhou, T.; Wang, H.; Shi, L.; Liu, Z.; Joo, S.W. An enhanced electroosmotic micromixer with an efficient asymmetric lateral structure. Micromachines 2016, 7, 218. [Google Scholar] [CrossRef] [Green Version]
  20. Neuringer, J.L.; Rosensweig, R.E. Ferrohydrodynamics. Phys. Fluids 1964, 7, 1927–1937. [Google Scholar] [CrossRef]
  21. Cai, G.; Xue, L.; Zhang, H.; Lin, J. A review on micromixers. Micromachines 2017, 8, 274. [Google Scholar] [CrossRef]
  22. Andersson, H.; Valnes, O. Flow of a heated ferrofluid over a stretching sheet in the presence of a magnetic dipole. Acta Mech. 1998, 128, 39–47. [Google Scholar] [CrossRef]
  23. Kafoussias, N.; Tzirtzilakis, E. Biomagnetic fluid flow over a stretching sheet with non linear temperature dependent magnetization. Z. Für Angew. Math. Und Phys. ZAMP 2003, 54, 551–565. [Google Scholar] [CrossRef]
  24. Tzirtzilakis1, E.; Sakalis, V.; Kafoussias, N.; Hatzikonstantinou, P. Biomagnetic fluid flow in a 3D rectangular duct. Int. J. Numer. Methods Fluids 2004, 44, 1279–1298. [Google Scholar] [CrossRef]
  25. Tzirtzilakis, E.; Kafoussias, N. Three-dimensional magnetic fluid boundary layer flow over a linearly stretching sheet. J. Heat Transf. 2010, 132, 011702. [Google Scholar] [CrossRef]
  26. Tzirtzilakis, E. A mathematical model for blood flow in magnetic field. Phys. Fluids 2005, 17, 077103. [Google Scholar] [CrossRef]
  27. Tsai, T.H.; Liou, D.S.; Kuo, L.S.; Chen, P.H. Rapid mixing between ferro-nanofluid and water in a semi-active Y-type micromixer. Sens. Actuators A Phys. 2009, 153, 267–273. [Google Scholar] [CrossRef]
  28. Wen, C.Y.; Liang, K.P.; Chen, H.; Fu, L.M. Numerical analysis of a rapid magnetic microfluidic mixer. Electrophoresis 2011, 32, 3268–3276. [Google Scholar] [CrossRef]
  29. Cao, Q.; Han, X.; Li, L. An active microfluidic mixer utilizing a hybrid gradient magnetic field. Int. J. Appl. Electromagn. Mech. 2015, 47, 583–592. [Google Scholar] [CrossRef]
  30. Sharifi, A.; Yekani Motlagh, S.; Badfar, H. Investigation of the effects of two parallel wires’ non-uniform magnetic field on heat and biomagnetic fluid flow in an aneurysm. Int. J. Comput. Fluid Dyn. 2018, 32, 248–259. [Google Scholar] [CrossRef]
  31. Bayareh, M.; Usefian, A.; Ahmadi Nadooshan, A. Rapid mixing of Newtonian and non-Newtonian fluids in a three-dimensional micro-mixer using non-uniform magnetic field. J. Heat Mass Transf. Res. 2019, 6, 55–61. [Google Scholar]
  32. Tokas, S.; Zunaid, M.; Ansari, M.A. Non-Newtonian fluid mixing in a Three-Dimensional spiral passive micromixer. Mater. Today Proc. 2021, 47, 3947–3952. [Google Scholar] [CrossRef]
  33. Saadat, M.; Shafii, M.B.; Ghassemi, M. Numerical investigation on mixing intensification of ferrofluid and deionized water inside a microchannel using magnetic actuation generated by embedded microcoils for lab-on-chip systems. Chem. Eng. Process. Process Intensif. 2020, 147, 107727. [Google Scholar] [CrossRef]
  34. Badfar, H.; Motlagh, S.Y.; Sharifi, A. Numerical Simulation of Magnetic Drug Targeting to the Stenosis Vessel Using Fe3O4 Magnetic Nanoparticles Under the Effect of Magnetic Field of Wire. Cardiovasc. Eng. Technol. 2020, 11, 162–175. [Google Scholar] [CrossRef] [PubMed]
  35. Sharifi, A.; Motlagh, S.Y.; Badfar, H. Numerical investigation of magnetic drug targeting using magnetic nanoparticles to the Aneurysmal Vessel. J. Magn. Magn. Mater. 2019, 474, 236–245. [Google Scholar] [CrossRef]
  36. Chen, X.; Zhang, L. A review on micromixers actuated with magnetic nanomaterials. Microchim. Acta 2017, 184, 3639–3649. [Google Scholar] [CrossRef]
  37. Cao, Q.; Fan, Q.; Chen, Q.; Liu, C.; Han, X.; Li, L. Recent advances in manipulation of micro-and nano-objects with magnetic fields at small scales. Mater. Horiz. 2020, 7, 638–666. [Google Scholar] [CrossRef]
  38. Meng, X.; Bachmann, M.; Artinov, A.; Rethmeier, M. The influence of magnetic field orientation on metal mixing in electromagnetic stirring enhanced wire feed laser beam welding. J. Mater. Process. Technol. 2021, 294, 117135. [Google Scholar] [CrossRef]
  39. Orthmann, J.; Kolb, A. Temporal Blending for Adaptive SPH; Computer Graphics Forum; Wiley Online Library: Malden, NJ, USA, 2012; Volume 31, pp. 2436–2449. [Google Scholar]
  40. Robinson, M.; Cleary, P.; Monaghan, J. Analysis of mixing in a twin cam mixer using Smoothed Particle Hydrodynamics. AIChE J. 2008, 54, 1987–1998. [Google Scholar] [CrossRef]
  41. Lenaerts, T.; Dutré, P. Mixing Fluids and Granular Materials; Computer Graphics Forum; Wiley Online Library: Malden, NJ, USA, 2009; Volume 28, pp. 213–218. [Google Scholar]
  42. Abdolahzadeh, M.; Tayebi, A.; Omidvar, P. Thermal effects on two-phase flow in 2D mixers using SPH. Int. Commun. Heat Mass Transf. 2021, 120, 105055. [Google Scholar] [CrossRef]
  43. Tayebi, A.; Shekari, Y.; Heydari, M.H. A meshless method for solving two-dimensional variable-order time fractional advection–diffusion equation. J. Comput. Phys. 2017, 340, 655–669. [Google Scholar] [CrossRef]
  44. Liu, F.; Zhang, K.; Liu, Z. Three-dimensional MLS-based numerical manifold method for static and dynamic analysis. Eng. Anal. Bound. Elem. 2019, 109, 43–56. [Google Scholar] [CrossRef]
  45. Tayebi, A.; Jin, Y.C. Development of moving particle explicit (MPE) method for incompressible flows. Comput. Fluids 2015, 117, 1–10. [Google Scholar] [CrossRef]
  46. Hon, Y.C.; Šarler, B.; Yun, D.F. Local radial basis function collocation method for solving thermo-driven fluid-flow problems with free surface. Eng. Anal. Bound. Elem. 2015, 57, 2–8. [Google Scholar] [CrossRef]
  47. Mramor, K.; Vertnik, R.; Šarler, B. Application of the local RBF collocation method to natural convection in a 3D cavity influenced by a magnetic field. Eng. Anal. Bound. Elem. 2020, 116, 1–13. [Google Scholar] [CrossRef]
  48. Liu, G.R.; Gu, Y.T. An Introduction to Meshfree Methods and Their Programming; Springer Science & Business Media: Dordrecht, The Netherland, 2005. [Google Scholar]
  49. Gingold, R.A.; Monaghan, J.J. Smoothed particle hydrodynamics: Theory and application to non-spherical stars. Mon. Not. R. Astron. Soc. 1977, 181, 375–389. [Google Scholar] [CrossRef]
  50. Lucy, L.B. A numerical approach to the testing of the fission hypothesis. Astron. J. 1977, 82, 1013–1024. [Google Scholar] [CrossRef]
  51. Jiang, T.; Ouyang, J.; Li, Q.; Ren, J.; Yang, B. A corrected smoothed particle hydrodynamics method for solving transient viscoelastic fluid flows. Appl. Math. Model. 2011, 35, 3833–3853. [Google Scholar] [CrossRef]
  52. Kiara, A.; Hendrickson, K.; Yue, D.K. SPH for incompressible free-surface flows. Part II: Performance of a modified SPH method. Comput. Fluids 2013, 86, 510–536. [Google Scholar] [CrossRef]
  53. Monaghan, J.J. Simulating free surface flows with SPH. J. Comput. Phys. 1994, 110, 399–406. [Google Scholar] [CrossRef]
  54. Cleary, P.W. Modelling confined multi-material heat and mass flows using SPH. Appl. Math. Model. 1998, 22, 981–993. [Google Scholar] [CrossRef] [Green Version]
  55. Monaghan, J.; Kocharyan, A. SPH simulation of multi-phase flow. Comput. Phys. Commun. 1995, 87, 225–235. [Google Scholar] [CrossRef]
  56. Aristodemo, F.; Federico, I.; Veltri, P.; Panizzo, A. Two-phase SPH modelling of advective diffusion processes. Environ. Fluid Mech. 2010, 10, 451–470. [Google Scholar] [CrossRef]
  57. Gobby, D.; Angeli, P.; Gavriilidis, A. Mixing characteristics of T-type microfluidic mixers. J. Micromech. Microeng. 2001, 11, 126. [Google Scholar] [CrossRef]
  58. Erickson, D.; Li, D. Microchannel flow with patchwise and periodic surface heterogeneity. Langmuir 2002, 18, 8949–8959. [Google Scholar] [CrossRef]
  59. Bedram, A.; Moosavi, A.; Hannani, S.K. Analytical relations for long-droplet breakup in asymmetric T junctions. Phys. Rev. E 2015, 91, 053012. [Google Scholar] [CrossRef] [Green Version]
  60. Violeau, D. Fluid Mechanics and the SPH Method: Theory and Applications; Oxford University Press: Oxford, UK, 2012. [Google Scholar]
  61. Bonet, J.; Lok, T.S. Variational and momentum preservation aspects of smooth particle hydrodynamic formulations. Comput. Methods Appl. Mech. Eng. 1999, 180, 97–115. [Google Scholar] [CrossRef]
  62. Monaghan, J.J. Smoothed particle hydrodynamics. Rep. Prog. Phys. 2005, 68, 1703. [Google Scholar] [CrossRef]
  63. Monaghan, J.J.; Kos, A. Solitary waves on a Cretan beach. J. Waterw. Port Coastal Ocean Eng. 1999, 125, 145–155. [Google Scholar] [CrossRef]
  64. Abdolahzadeh, M.; Tayebi, A.; Mansouri Mehryan, M. Numerical Simulation of Mixing in Active Micromixers Using SPH. Transp. Porous Media 2022, 1–18. [Google Scholar] [CrossRef]
  65. Lee, E.S.; Moulinec, C.; Xu, R.; Violeau, D.; Laurence, D.; Stansby, P. Comparisons of weakly compressible and truly incompressible algorithms for the SPH mesh free particle method. J. Comput. Phys. 2008, 227, 8417–8436. [Google Scholar] [CrossRef]
  66. Shadloo, M.S.; Zainali, A.; Sadek, S.H.; Yildiz, M. Improved incompressible smoothed particle hydrodynamics method for simulating flow around bluff bodies. Comput. Methods Appl. Mech. Eng. 2011, 200, 1008–1020. [Google Scholar] [CrossRef]
  67. Xu, R.; Stansby, P.; Laurence, D. Accuracy and stability in incompressible SPH (ISPH) based on the projection method and a new approach. J. Comput. Phys. 2009, 228, 6703–6725. [Google Scholar] [CrossRef]
  68. Tayebi, A.; Dehkordi, B.G. Development of a piso-sph method for computing incompressible flows. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2014, 228, 481–490. [Google Scholar] [CrossRef]
  69. Morris, J.P. Analysis of Smoothed Particle Hydrodynamics with Applications; Monash University Australia: Melbourne, Australia, 1996. [Google Scholar]
  70. Kim, Y.D.; An, S.J.; Maeng, J.S. Numerical analysis of the fluid mixing behaviors in a microchannel with a circular cylinder and an oscillating stirrer. J. Korean Phys. Soc. 2007, 50, 505–513. [Google Scholar]
Figure 1. Schematic of the physical models of the micromixers.
Figure 1. Schematic of the physical models of the micromixers.
Fluids 07 00341 g001
Figure 2. Comparison of the variation of concentration versus the dimensionless length of the channel in the work of Kim et al. [70].
Figure 2. Comparison of the variation of concentration versus the dimensionless length of the channel in the work of Kim et al. [70].
Fluids 07 00341 g002
Figure 3. Velocity contours for different number of wires at R e = 15 .
Figure 3. Velocity contours for different number of wires at R e = 15 .
Fluids 07 00341 g003
Figure 4. Concentration (left) and particle contours (right) for a different number of wires at R e = 15 .
Figure 4. Concentration (left) and particle contours (right) for a different number of wires at R e = 15 .
Fluids 07 00341 g004
Figure 5. Comparison of the mixing efficiency versus the Reynolds numbers for a different number of electrical wires.
Figure 5. Comparison of the mixing efficiency versus the Reynolds numbers for a different number of electrical wires.
Fluids 07 00341 g005
Figure 6. Variations of the mixing index versus the dimensionless length of the microchannel for various Reynolds numbers.
Figure 6. Variations of the mixing index versus the dimensionless length of the microchannel for various Reynolds numbers.
Fluids 07 00341 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abdolahzadeh, M.; Tayebi, A.; Ahmadinejad, M.; Šarler, B. Numerical Simulation of Mixing Fluid with Ferrofluid in a Magnetic Field Using the Meshless SPH Method. Fluids 2022, 7, 341. https://doi.org/10.3390/fluids7110341

AMA Style

Abdolahzadeh M, Tayebi A, Ahmadinejad M, Šarler B. Numerical Simulation of Mixing Fluid with Ferrofluid in a Magnetic Field Using the Meshless SPH Method. Fluids. 2022; 7(11):341. https://doi.org/10.3390/fluids7110341

Chicago/Turabian Style

Abdolahzadeh, Mohsen, Ali Tayebi, Mehrdad Ahmadinejad, and Božidar Šarler. 2022. "Numerical Simulation of Mixing Fluid with Ferrofluid in a Magnetic Field Using the Meshless SPH Method" Fluids 7, no. 11: 341. https://doi.org/10.3390/fluids7110341

APA Style

Abdolahzadeh, M., Tayebi, A., Ahmadinejad, M., & Šarler, B. (2022). Numerical Simulation of Mixing Fluid with Ferrofluid in a Magnetic Field Using the Meshless SPH Method. Fluids, 7(11), 341. https://doi.org/10.3390/fluids7110341

Article Metrics

Back to TopTop