Novel Approaches Concerning the Numerical Modeling of Particle and Cell Separation in Microchannels: A Review

: The demand for precise separation of particles, cells, and other biological matter has signiﬁcantly increased in recent years, leading to heightened scientiﬁc interest in this topic. More recently, due to advances in computational techniques and hardware, numerical simulations have been used to guide the design of separation devices. In this article, we establish the theoretical basis governing ﬂuid ﬂow and particle separation and then summarize the computational work performed in the ﬁeld of particle and cell separation in the last ﬁve years with an emphasis on magnetic, dielectric, and acoustic methods. Nearly 70 articles are being reviewed and categorized depending on the type of material separated, ﬂuid medium, software used, and experimental validation, with a brief description of some of the most notable results. Finally, further conclusions, future guidelines, and suggestions for potential improvement are highlighted. Coulter counter array. The electrical signal is decoded to quantify the spatial distribution of cells, and an expression histogram is produced via computational analysis of sensor signals.


Introduction
The capacity to precisely separate and manipulate suspended particles, cells, or other biological matter is a critical process in many applications of biology and medicine for disease diagnosis, genetic analysis, drug screening and delivery, tissue engineering, and therapeutics [1][2][3][4]. Recent research has been devoted to the manipulation, concentration, focusing, separation, isolation, and fractionation of particulate material such as micro-and nanoparticles, cells, liposomes, microvesicles, viruses, etc., using microfluidic devices, i.e., the so-called lab on chip devices [5][6][7]. Microfluidics-based cell/particle manipulation technologies have shown great potential due to the low sample and reagent volume consumption, low waste generation, high product purity, high sensitivity and selectivity, ease of use, and short isolation time, in addition to being able to perform particle/cell analysis at the single-cell level [8][9][10][11][12][13]. Moreover, downscaling of sorting systems enables less laboratory space, and the integration of multiple tasks and operations within the same device is possible, which not only increases the precision of the analysis but also improves the accuracy, reliability, and reproducibility of sample preparation procedures [14,15]. Thus, controlling the motion and position of these materials within fluids in microchannels possesses great importance in the previous fields, especially for the development of low-cost and portable diagnostic tools that can be used at the point of need [16][17][18].
Among the different approaches developed to separate and manipulate particles and cells, the use of active forces (i.e., magnetophoresis, dielectrophoresis, or acoustophoresis) presents several advantages over passive techniques. Passive techniques are usually restricted to separating the material based on size or compressibility, which may cause cell damage or channel obstruction; they require complex microchannels that lack reconfigurability to suit various target entities, large device footprints, long operation times, and low accuracy [14,19,20]. Due to these limitations, active sorting techniques are widely investigated as potential alternatives. Active sorting methods utilize an external field to Figure 1. Number of documents published in the last two decades (up to 2021) related to the us of numerical modeling for the microfluidic separation of particles and cells, and found in the S pus database using the keywords ''cell separation", "particle separation", "simulation", and "m crofluidics".

Principles and Numerical Modeling of the Active Microfluidic Particle Separatio
In this section, the underlying mechanisms behind active particle and cell separa are established. First, we briefly describe the different approaches used to model the p cess. Second, the main forces acting on the material during separation are introduced nally, the numerical methods most widely applied are briefly described. It should noted that the theoretical models described in this section can be implemented into dif ent commercial and custom solvers that describe the particle and cell motion as a func of multiple variables and parameters. Thus, they enable the design and optimization the separation process.

Approaches to Model the Microfluidic Separation of Suspended Particles/Cells
It should be noted that there are mainly two approaches to model the separatio particles and cells in microfluidics. First is the Lagrange method, the most widely u approach to model the active separation of particles and cells in microfluidics. This proach neglects Brownian motion and solves the Newtonian equation of motion of sin particles [77]. The second approach involves using an Eulerian approach in which pa cles are modeled collectively in terms of a time-dependent spatially varying concen tion, which is calculated by a partial differential equation that accounts for the drift fo and the Brownian diffusivity of the particles [78]. The Euler approach is generally u for nanoscale particle separation, where the effects of Brownian motion are signific Due to the size of the particles employed in most studies dealing with microfluidic se ration, which in most cases is hundreds of nanometers, the Brownian forces are neglig [61]. Thus, in this work, we only consider the Lagrangian approach, where particles modeled as discrete units, and the trajectory of each particle is estimated by applying classical Newtonian dynamics [26,52,71,79]: m p dv p dt = F d + F g + F separation where mp and vp denote the mass and velocity of the cell/particle and Fd, Fg and Fsepa represent the drag force, gravity force, and the force driving the separation, i.e., the ac force. It is to be mentioned that since the characteristic timescale of the acceleration ph of the microparticles' and cell's motion is much smaller than the time scale of the varia of the external forces, the particle always moves at a terminal velocity which can Figure 1. Number of documents published in the last two decades (up to 2021) related to the use of numerical modeling for the microfluidic separation of particles and cells, and found in the Scopus database using the keywords "cell separation", "particle separation", "simulation", and "microfluidics".

Principles and Numerical Modeling of the Active Microfluidic Particle Separation
In this section, the underlying mechanisms behind active particle and cell separation are established. First, we briefly describe the different approaches used to model the process. Second, the main forces acting on the material during separation are introduced. Finally, the numerical methods most widely applied are briefly described. It should be noted that the theoretical models described in this section can be implemented into different commercial and custom solvers that describe the particle and cell motion as a function of multiple variables and parameters. Thus, they enable the design and optimization of the separation process.

Approaches to Model the Microfluidic Separation of Suspended Particles/Cells
It should be noted that there are mainly two approaches to model the separation of particles and cells in microfluidics. First is the Lagrange method, the most widely used approach to model the active separation of particles and cells in microfluidics. This approach neglects Brownian motion and solves the Newtonian equation of motion of single particles [77]. The second approach involves using an Eulerian approach in which particles are modeled collectively in terms of a time-dependent spatially varying concentration, which is calculated by a partial differential equation that accounts for the drift force and the Brownian diffusivity of the particles [78]. The Euler approach is generally used for nanoscale particle separation, where the effects of Brownian motion are significant. Due to the size of the particles employed in most studies dealing with microfluidic separation, which in most cases is hundreds of nanometers, the Brownian forces are negligible [61]. Thus, in this work, we only consider the Lagrangian approach, where particles are modeled as discrete units, and the trajectory of each particle is estimated by applying the classical Newtonian dynamics [26,52,71,79]: where m p and v p denote the mass and velocity of the cell/particle and F d , F g and F separation represent the drag force, gravity force, and the force driving the separation, i.e., the active force. It is to be mentioned that since the characteristic timescale of the acceleration phase of the microparticles' and cell's motion is much smaller than the time scale of the variation of the external forces, the particle always moves at a terminal velocity which can be obtained from the balance of gravity, drag, and separation forces [61]. Depending on the specific application and working conditions, other forces might need to be included (such as particle-fluid or particle-particle interactions), although in most scenarios, including the drag and separation forces will suffice to estimate the separation efficiency of the device. It is widely adopted that when the volume fraction of the dispersed particles is lower than 10%, the particle-particle interaction can be neglected [27]. More specifically, particleparticle and particle-fluid interaction can be ignored for particle suspensions with a small particle volume concentration of less than 5%, which is the case in most applications [31]. Equations to develop these forces are in the following.

Fluid and Gravity Forces
The forces driving the particle and cell separation inside the microchannel are balanced with the drag force on the particle in the fluid. For a spherical particle in a laminar flow regime, the Stokes drag force is [19,52,80]: where η is the fluid viscosity, r is the particle's radius, and v is the velocity vector for the particle relative to the fluid.
To compute the drag force, the fluid velocity field within the device needs to be calculated. The commercial simulation software packages used to simulate the process solve the Navier-Stokes and continuity equations using an Eulerian approach [27,79]: where v f and ρ f represent the fluid velocity and density, and P is pressure. A more sophisticated Lagrangian tracking of the particles and cells can be achieved by coupling the particles with the Eulerian solution of the Navier-Stokes equations for the carrier fluid by accounting for the particle-fluid momentum exchange between phases [27,33]. The resulting model, involving the simultaneous solutions for the Eulerian fluid flow and the Lagrangian tracking of dispersed particles, is usually referred to as the Discrete Particle Model (DPM). For Lagrangian-Eulerian models, the discrete and continuous phases are simulated using different numerical discretization and gridding schemes [27]. With that, the iterative exchange of variables between phases can only be achieved using approximated Lagrangian-Eulerian transformations and, therefore, adds for numerical errors [27].
Finally, and although it is neglected in most of the studies, the gravity force, including buoyancy, can be computed as the difference between the weight of an immersed particle and an upward buoyancy force exerted by fluid, given by [31]: where V p represents the particle volume, ρ p is the particle density, and g is the acceleration due to gravity.

Separation Forces
Different separation forces can be applied to manipulate or separate solid particles and cells from fluids. Among the most common active separation techniques are those that employ magnetic, electric, and acoustic fields. These forces are presented below.

Magnetic Forces
Magnetophoresis is one of the most promising approaches to separating particles and cells within microfluidics [81]. It can be defined as the migration of a material in the presence of a nonuniform magnetic field. The magnetic force acting on unsaturated particles or paramagnetic cells can be expressed as [27]: where ∆χ is the difference between the medium and the particle/cell magnetic susceptibility, µ 0 is the magnetic permeability of the space, and H represents the intensity of the applied magnetic field. If the magnetic susceptibility of the particle is greater than that of the carrier liquid, ∆χ > 0, and positive magnetophoresis occurs, i.e., the magnetic particle is attracted towards the high magnetic field potential. On the contrary, in negative magnetophoresis, the migration of non-magnetic (diamagnetic) particles in a magnetic fluid is carried out, and ∆χ < 0 [38]. Negative magnetophoresis can separate diamagnetic particles such as polymer beads and biological material, including cells in a paramagnetic medium such as a ferrofluid solution, without the need for magnetic labels.
It should be noted that in the majority of studies, positive magnetophoresis is employed, and thus, superparamagnetic beads are commonly the material that needs to be separated. Superparamagnetic beads are soft-magnetic materials that do not retain their magnetic moment after the applied magnetic field is removed [31,79]. The magnetization of the bead is a function of the applied magnetic field. It is a linear function of the field (χH) up to a magnetic flux of approximately 0.1-0.5 T [31,82], after which it remains nearly constant, with a saturation magnetization of M sat that is independent of the applied magnetic field [31]. In the saturation region, Equation (6) needs to be replaced by the following equation [82]: It should be noted that magnetic cell separation devices are not sensitive to factors such as ionic concentration, pH, and surface charge, allowing operation in a wide range of these parameters [31,79]. However, in most cases, separating biomaterials and cells requires the prior labeling of the material using magnetic beads. This is because biological materials (except for red blood cells (RBCs) and iron proteins or magnetostatic bacteria) usually exhibit diamagnetic properties and are not greatly affected by the applied magnetic field (unless suspended in a paramagnetic medium or ferrofluid) [31]. In these cases, the magnetic beads are specifically designed to target the cell or biomolecule that needs to be separated. Magnetic beads with diameters of 50-2500 nm, made of superparamagnetic iron oxide nanoparticles (SPIONs) embedded in a polymer matrix (e.g., polystyrene), with a net magnetic volume fraction typically less than 20% vol., are commercially available and purchased in most of the recently published research works from manufacturers [29]. Before the magnetic separation, these beads are incubated with the biological material of interest. The resulting biological material-bead complex is not spherical, making it difficult to calculate some forces, such as the drag force. However, depending on the sizes of the materials under consideration, simple approaches can be used to estimate the forces acting on the complex without introducing a significant error [31].
Finally, implementing magnetic functionality into the microfluidic devices can be done by several means. The most common approach is to place magnetic elements in the vicinity of the microchannel [79]. Such elements produce a magnetic field gradient inside the channel, which exerts an attraction force on the labeled biomaterials or the magnetic beads as they flow through the channel ( Figure 2). Microfluidic devices engineered with permanent magnets enable on-chip manipulation of magnetically labeled biomaterials which requires no power consumption [31]. Other approaches using electromagnets have been developed as well, but their disadvantages are the required circuitry and the possible temperature raise due to the Joule heating, which may pose concerns if the biological material is temperature sensitive [79].
Processes 2022, 10, x FOR PEER REVIEW 6 of 24 temperature raise due to the Joule heating, which may pose concerns if the biological material is temperature sensitive [79].

Electric Forces
Electric fields can separate particles via electrophoresis (EP) and dielectrophoresis (DEP). EP is not covered in this work, as it is not generally applied for separating particles and cells since it relies on the separation of charged materials under the influence of an electric field [45]. DEP is the most common technique used to separate particles and cells because it does not require the material to be electrically charged [45]. Thus, we will review only the simulation studies that describe particle and cell DEP.
Analogously to magnetophoresis, DEP induces a cell/particle movement in a nonuniform electric field based on the dielectric properties (electrical conductivity and permittivity) of the cell and the surrounding medium [45]. The exerted force from the nonuniform electric field on a dielectric material (with zero net charge) polarizes the particles to generate corresponding induced charges on the surface (Figure 3), which generates, in turn, dipole moments, and due to the force, the ends of the positive and negative charges are uneven [50]. This induced dipole in a nonuniform electric field causes the particle to be pulled electrostatically along the electric field gradient when it approaches the gap between the electrodes [26,45].  [79]. Reprinted with permission from Ref. [79], 2022, Elsevier.

Electric Forces
Electric fields can separate particles via electrophoresis (EP) and dielectrophoresis (DEP). EP is not covered in this work, as it is not generally applied for separating particles and cells since it relies on the separation of charged materials under the influence of an electric field [45]. DEP is the most common technique used to separate particles and cells because it does not require the material to be electrically charged [45]. Thus, we will review only the simulation studies that describe particle and cell DEP.
Analogously to magnetophoresis, DEP induces a cell/particle movement in a nonuniform electric field based on the dielectric properties (electrical conductivity and permittivity) of the cell and the surrounding medium [45]. The exerted force from the nonuniform electric field on a dielectric material (with zero net charge) polarizes the particles to generate corresponding induced charges on the surface (Figure 3), which generates, in turn, dipole moments, and due to the force, the ends of the positive and negative charges are uneven [50]. This induced dipole in a nonuniform electric field causes the particle to be pulled electrostatically along the electric field gradient when it approaches the gap between the electrodes [26,45].  More specifically, at material interfaces, this polarization is revealed by the formation of macroscopic multipoles. If the excitatory field is inhomogeneous, the Coulomb forces acting on each half of the particles are unequal, which causes a net dielectrophoretic force  More specifically, at material interfaces, this polarization is revealed by the formation of macroscopic multipoles. If the excitatory field is inhomogeneous, the Coulomb forces acting on each half of the particles are unequal, which causes a net dielectrophoretic force acting on the particle, as seen in Figure 3 [61]. The magnitude and direction of the DEP force can be obtained by computing the cell's or particle's effective dipole moment [44]. The DEP force acting on the imposed electric field of an induced dipole is given by [9,44,47,80]: where where ε m is the permittivity of the medium, K(ω) is the Clausius-Mossotti factor that depends on the angular frequency of the applied field ω, E is the electric field, and ∈ * p and ∈ * m are the complex permittivities of the particle/cell and the medium, which are a function of the permittivity and conductivity of the material and the angular frequency of the applied electric field ω [15,80]. Either Direct or Alternating Current (DC or AC, respectively) sources can be employed to create a nonuniform electric field [46], but AC is chosen in most of the studies due to the low input voltage requirements, which avoids Joule heating effects and cell/particle damage [47,54].
The value of Re[K(ω)] varies in a range from −0.5 to 1 [26]. According to it, DEP is frequently classified into positive DEP (pDEP) if Re[K(ω)] > 0 and negative DEP (nDEP) if Re[K(ω)] < 0. As a result, the particles forced by pDEP response are guided toward the location of the high electric field strength (around the edges of the electrodes). In contrast, the particles forced by nDEP response are moved to the low electric gradients (away from the electrodes). In other words, if the polarizability of the particle is greater than that of the liquid, it makes the particle move towards the region with the highest electric field gradient [45]. For many materials, the real part of K(ω) is negative over a wider range of low to very high frequencies, and it is positive at very low frequencies. The frequency at which the material experiences neither p-DEP nor n-DEP is the crossover frequency [47]. This is a specific frequency at which the intrinsic properties of particles and cells can be defined [15]. For the separation or isolation of cells and biological material, it is desirable to employ nDEP so that the material is directed towards the low electric field region to minimize unwanted transfection and disruption of cell structure/membrane [11]. It should also be noted that while the electric field induces DEP forces on the particles, it also induces an electro-osmotic and electrothermal flow that might affect the separation [45].
The DEP devices consist of two main parts integrated into a substrate: the electrodes, which can be of different types and generate the nonuniform electric field inside the microchannel, and the microchannel itself, which contains the particle/cell suspension. Planar and three-dimensional (3D) electrode structures are commonly used for these applications. 3D electrodes are fabricated on the top and bottom, or sidewalls of microfluidic channels, whereas planar electrodes are commonly embedded on the bottom of microfluidic channels [15].
Usually, the separation of particles can be easily modeled once their properties are characterized. On the contrary, cells generally have a relatively complex internal structure that can affect their dielectric properties, and equivalent cells of different structures produce different results. Some can be directly equivalent to a solid sphere, and some need to establish an equivalent model for biological particles. Among them, the concentric spherical model is more common [9,15,50]. With this model, the cells are assumed to possess two concentric layers of various electric and dielectric properties [15]. In any case, the electric fields can be effectively used to analyze and separate particles and cells based on their electric properties. Indeed, DE properties of cells have been linked to changes in plasma membrane capacitance, permeability, and change in Ca + and K + ion levels in the cytoplasm, which can be related to their response to different drugs, infection, even cancer cell aggressiveness, etc.; these result in plasma membrane and cytoplasm electrical conductivity changes [44,80].

Acoustic Forces
Acoustic forces originate from the scattering of imposed ultrasound waves from objects or particles suspended in a fluid. The scattering of waves occurs in a fluid when particles suspended in it are exposed to an acoustic field. Due to this scattering, second-order forces arise, known as acoustic radiation forces. The radiation forces are time-averaged forces that can effectively transport the suspended particles. The direction of such motion of a particle depends upon its properties and the surrounding fluid medium (density and compressibility) [14,16]. When such force on particles overcomes the viscous drag, gravity, and other forces, particles may exhibit acoustophoretic motion irrespective of their inherent electrical or magnetic properties [16]. Particles with varying sizes or physical parameters resulting in different acoustic properties will be exposed to different acoustic radiation pressures. They will move to specific regions of the channel (pressure nodes or antinodes) at different intervals, resulting in distinct identifiers for separation and diagnosis [2,71].
An acoustic cell sorter consists of a piezoelectric substrate with interdigital transducers (IDTs) to generate surface acoustic waves (SAWs) and a microfluidic channel bonded to the substrate, as seen in Figure 4. SAWs are mechanical waves that primarily propagate upon the surface of an elastic material. These waves consist of a longitudinal compression motion coupled with a transversal shear motion [8]. It should be noted that bulk acoustic wave devices generated by substrate-bonded bulk transducers are outside the scope of this work due to their unsuitability for microfluidic devices and high operating frequency and energy requirements [8,19,75].  When a surface standing acoustic wave (SSAW) is applied between the opposite walls of the channel (which act as reflectors) by exciting the piezoelectric elements, a harmonic standing pressure field is formed in the fluid. In most cases, a half-wavelength standing wave is usually generated in the lateral direction of the fluid channel (i.e., perpendicular to the continuous flow direction) with a pressure node in the channel center and pressure antinodes located at the channel walls [24]. Besides SSAW techniques, there exist also traveling surface acoustic waves (TSAW) sorting approaches, which change the pressure of the fluid and drive the particles to move to the other side. However, TSAW devices are not as common as SSAW devices. More specifically, TSAWs are waves that propagate in one direction and push the dispersed particles/cells in the direction of prop- When a surface standing acoustic wave (SSAW) is applied between the opposite walls of the channel (which act as reflectors) by exciting the piezoelectric elements, a harmonic standing pressure field is formed in the fluid. In most cases, a half-wavelength standing wave is usually generated in the lateral direction of the fluid channel (i.e., perpendicular to the continuous flow direction) with a pressure node in the channel center and pressure antinodes located at the channel walls [24]. Besides SSAW techniques, there exist also traveling surface acoustic waves (TSAW) sorting approaches, which change the pressure of the fluid and drive the particles to move to the other side. However, TSAW devices are not as common as SSAW devices. More specifically, TSAWs are waves that propagate in one direction and push the dispersed particles/cells in the direction of propagation. In contrast, SSAWs are formed in the volume of the microchannel as a result of the superposition of two traveling waves, which move towards each other and have some phase shift [75]. The design and electrical signal applied to IDTs completely defines the frequency, amplitude, and orientation of the generated acoustic waves; however, the microchannel dimensions and layout, together with its material and liquid phase composition, will affect the acoustic wave propagation and fluid-particle interaction [75].
The acoustic radiation force F rad exerted on compressible particles/cells suspended in an inviscid fluid arising from an acoustic field can be calculated using the Gorkov's potential approach as [5,14,16]: where where U rad is the potential in the acoustic domain, p in and v in are the first-order pressure and velocity of the incoming acoustic waves. The density and compressibility of the fluid and the particles are denoted by ρ f , ρ p , k f , and k p , respectively, and the factors f 1 and f 2 are the dimensionless scattering coefficients [5,16]. In general, for most applications, compared to smaller and less dense particles, larger and denser particles move faster towards the pressure nodes for higher acoustic radiation forces. Thus microparticles of different physical properties will end up in different lateral positions and be separated into different outlets at the end of the acoustic separation channel [24]. This formulation is valid when inviscid fluids and spherical particles (smaller than the acoustic field wavelength λ) are considered. In most cases, inviscid fluids are assumed even though viscous fluids are employed [1]. For cases where this assumption cannot be adopted, a corrected expression for f 2 must be considered (see Ref. [16]).
Regarding the particle or cell movement, the acoustophoretic contrast factor is calculated to determine the direction of the force. This is given as: If the material has a φ value higher than 0, the direction of the force is towards the nodes of the standing pressure wave (minima U rad ), whereas a negative φ implies that the force acts towards the pressure antinodes.

Numerical Methods
Different methods can be used to solve the partial differential equations, such as the force equations and particle/cell migration equation within the device. In this section, the Finite Element Method (FEM) and Finite Volume Method (FVM) are introduced since they are the most common methods and the ones integrated into commercial simulation software. The exact details of how the constituent equations are solved for each method are beyond the scope of this article. Still, a general overview is provided, along with the advantages and disadvantages of each method.

Finite Element Method
For the FEM, the computational domain is subdivided into smaller sections, i.e., elements usually in the shape of a triangle or tetrahedron. The unknown variable/parameter is then solved for each element so that the differential equation is satisfied at each boundary. Some of the core advantages of the FEM method are its ability to handle domains that include complex geometrical features, due to the use of body-fitted meshing, and its accuracy since higher-order approximations are used [83]. The FEM method is applied in various commercially available computational software packages. For example, COMSOL Multiphysics (www.comsol.com accessed on 1 June 2022) is a simulation platform that provides fully coupled multiphysics and single-physics modeling capabilities, which is based on the FEM to realize the simulation of real physical phenomena by solving partial differential equations (groups) [50]. It is applied to the research and engineering calculation of various disciplines, simulating the physical processes of science and engineering. It represents one of the most common numerical analysis software packages employed in optimizing particle and cell separation in microchannels.

Finite Volume Method
The FVM is a subsection of the finite differences method (FDM) commonly used in many commercial CFD software packages. In FVM, the domain is split into rectangular volumes, i.e., cells with their edges aligned to the coordinate system used. The solved quantities are stored in the center of each volume for non-vectors (e.g., temperature) and at the edges of each volume for vectors (e.g., velocity). Some of the advantages of the FVM are its computational simplicity which makes it computationally possible to solve for large domains and ease of implementation [83]. There are a plethora of commercially available software packages that are based on the FVM, such as PHOENICS (www.cham.co.uk accessed on 1 June 2022), ANSYS Fluent (www.ansys.com accessed on 1 June 2022), FLOW-3D (www.flow3d.com accessed on 1 June 2022), and STAR-CCM+ (https://www.plm. automation.siemens.com/global/en/products/simcenter/STAR-CCM.html accessed on 1 June 2022). Most notably, OpenFOAM, an open-source software, uses the FVM to solve partial differential equations [58]. The software packages mentioned above have also been used to model particle separation in microfluidics.

Mesh Dependency Analysis
Independent of the software employed, it is required to assess the accuracy of the numerical results by performing a mesh dependency study. That, in principle, consists of running a given problem with a progressively refined mesh, i.e., gradually increasing the number of elements or cells until the difference between the results of the mesh independent domain and the immediately finer domain is less than a given value. Usually, errors below 5% are considered, but some applications might require more accurate simulations, and finer meshes are employed to ensure that the error stays below 1% [70]. The mesh dependency study must be performed as the last step of problem setup before proceeding to the parametric analysis.

Recent Simulation Studies
This section reviews the most recent simulation studies on the active separation of particles and cells inside microchannels. We start with magnetophoresis studies, followed by dielectrophoresis, and end with acoustophoresis. Table 1 presents the main simulation studies performed in the last 5 years on the microfluidic separation of particles and cells using magnetic means. Most of these works used a Lagrangian-Eulerian (particle-based) model, but other approaches have been explored, too [27], and COMSOL Multiphysics is the preferred platform to simulate the process. As seen in the table, one of the main challenges of magnetophoretic microfluidic devices is their typically small volume capacity, which implies the use of relatively low flow rates. Nevertheless, numerical methods can be employed to estimate the device performance and optimize the working conditions in advance of fabrication so that high magnetic field gradients and flow rates are achieved. For example, Castillo-Torres et al. [29], using a combination of modeling and experimental techniques, designed a microfluidic device able to separate bacteria magnetically conjugated to magnetic disks at flow rates higher than 100 µL/s using a magnet array. They achieved around 94% efficiency in the bacterial detection/separation, analyzing 100 mL of water samples in less than 15 min. Moreover, Golozar et al. modeled and optimized a magnetophoretic chip to isolate cells conjugated to magnetic beads [31]. Their design consisted of a continuous-flow microfluidic platform that contained locally engineered magnetic field gradients able to separate different cells previously conjugated to magnetic beads at flow rates up to 200 mL/h. The simulations were used to predict the trajectory and capturing region of the cells and to analyze the effect of the magnetic bead size, cell size, number of beads per cell, and flow rate on the separation performance. On the other hand, Civelekoglu et al. [36], with the aid of numerical simulations, introduced a cytometry technique that computed the surface expression of immunomagnet-ically labeled cells by electrically tracking their trajectory under a magnetic field gradient on a microfluidic chip with a throughput of >500 cells/min. More specifically, immunomagnetically labeled cells were deflected under a magnetic field gradient according to their magnetic load (or surface expression) and detected at the different outlets of the device ( Figure 5). Their device allowed direct expression profiling of target subpopulations from non-purified samples. It was demonstrated via application to measure epithelial cell adhesion molecule (EpCAM) expression on human breast cancer cells.

Magnetic Separation
Magnetic particles (40- On the other hand, Civelekoglu et al. [36], with the aid of numerical simulations, introduced a cytometry technique that computed the surface expression of immunomagnetically labeled cells by electrically tracking their trajectory under a magnetic field gradient on a microfluidic chip with a throughput of >500 cells/min. More specifically, immunomagnetically labeled cells were deflected under a magnetic field gradient according to their magnetic load (or surface expression) and detected at the different outlets of the device ( Figure 5). Their device allowed direct expression profiling of target subpopulations from non-purified samples. It was demonstrated via application to measure epithelial cell adhesion molecule (EpCAM) expression on human breast cancer cells. The separation of magnetic particles from fluids within microchannels has also been addressed for its implementation within nanoparticle synthesis procedures. The cleaning process of synthesized nanoparticles is usually performed by multi-step procedures involving filtrations, centrifugations, or manually washing with the help of a powerful magnet; these procedures are not only time-consuming but also inefficient as a large number of nanoparticles are often lost, and traces of solvents, surfactants, or byproducts remain after the process [37]. To overcome these limitations, Cardoso et al. [37] used a microfluidic device for the automated and continuous cleaning and separation of the material. Using numerical simulation and experimental testing, their microfluidic device cleaned and recovered more than 99% of iron oxide nanoparticles, obtaining a higher recovery than manual procedures (94%). Other studies have also focused on the sorting of the nanoparticles after synthesis to obtain fractions uniform in size [43]. Moreover, numerical simulation has also been employed to optimize devices where negative magnetophoresis occurs. Indeed, Munaz et al. [38] designed a device where different diamagnetic particles suspended in a ferrofluid are separated via negative magnetophoresis, as seen in Table 1. The authors used a numerical model to optimize the magnet array configuration and the effect of flow rate ratios, and the concentration distribution of ferrofluid in the device.
It should be noted that numerical models have also been applied to optimize the separation of cells/biomolecules based on their intrinsic magnetization (without using magnetic beads as labels). The separation of intrinsically magnetic (paramagnetic) cells or biomaterial is more challenging because of the orders of magnitude weaker magnetizations observed in these materials (much weaker than what can typically be achieved in immunomagnetically labeled cells) [34,84]. For example, Sun et al. [34] employed a separation system called magnetic deposition microscopy (MDM) to not only separate cells, but to deposit cells in specific locations on slides for further microscopic analysis based on their intrinsic magnetizations. These authors characterized the MDM system using finite The separation of magnetic particles from fluids within microchannels has also been addressed for its implementation within nanoparticle synthesis procedures. The cleaning process of synthesized nanoparticles is usually performed by multi-step procedures involving filtrations, centrifugations, or manually washing with the help of a powerful magnet; these procedures are not only time-consuming but also inefficient as a large number of nanoparticles are often lost, and traces of solvents, surfactants, or byproducts remain after the process [37]. To overcome these limitations, Cardoso et al. [37] used a microfluidic device for the automated and continuous cleaning and separation of the material. Using numerical simulation and experimental testing, their microfluidic device cleaned and recovered more than 99% of iron oxide nanoparticles, obtaining a higher recovery than manual procedures (94%). Other studies have also focused on the sorting of the nanoparticles after synthesis to obtain fractions uniform in size [43]. Moreover, numerical simulation has also been employed to optimize devices where negative magnetophoresis occurs. Indeed, Munaz et al. [38] designed a device where different diamagnetic particles suspended in a ferrofluid are separated via negative magnetophoresis, as seen in Table 1. The authors used a numerical model to optimize the magnet array configuration and the effect of flow rate ratios, and the concentration distribution of ferrofluid in the device.
It should be noted that numerical models have also been applied to optimize the separation of cells/biomolecules based on their intrinsic magnetization (without using magnetic beads as labels). The separation of intrinsically magnetic (paramagnetic) cells or biomaterial is more challenging because of the orders of magnitude weaker magnetizations observed in these materials (much weaker than what can typically be achieved in immunomagnetically labeled cells) [34,84]. For example, Sun et al. [34] employed a separation system called magnetic deposition microscopy (MDM) to not only separate cells, but to deposit cells in specific locations on slides for further microscopic analysis based on their intrinsic magnetizations. These authors characterized the MDM system using finite element and computational fluid mechanics software and predicted the separation performance of magnetic spores (Bacillus spores). Analogously, Kasetsirikul et al. [41] designed a device to separate malaria-infected RBCs, which show paramagnetism and suffer positive magnetophoresis when exposed to magnetic fields. Their numerical simulations suggest that the infected RBCs could be separated from healthy RBCs in a continuous flow using magnet arrays.
Numerical models have also been used to analyze the mixing of co-flowing fluids while magnetophoresis takes place. Gomez-Pastora et al. [32] employed a Lagrangian computational model for designing continuous-flow magnetophoretic microsystems for blood detoxification, assuring complete particle/biomaterial capture without compromising the quality of the blood (avoiding diffusion or loss of blood). The model was used to study the particle trajectory, the flow pattern and mass transfer between streams, and the analysis of multiple variables and parameters (flow rates, beads and magnets dimensions, fluid viscosities, etc.) was numerically and experimentally performed. In a later study [33], these authors provided a detailed analysis of the bead motion and its effect on fluid flow while magnetophoresis takes place ( Figure 6). Different scenarios were modeled and compared: (i) one-way coupling wherein momentum is transferred from the fluid to the beads, which are treated as point particles, (ii) two-way coupling wherein the beads are treated as point particles, and momentum is transferred from the bead to the fluid and vice versa, and (iii) two-way coupling taking into account the effects of bead volume in the fluid displacement. It was demonstrated that although there is little difference in the bead trajectories for the three scenarios, there is significant variation in the flow fields, especially when high magnetic forces are applied to the beads. Therefore, employing an accurate full flow-focused model that considers the effects of the bead motion and volume on the flow field when high magnetic forces are employed is recommended; however, when the beads are subjected to medium or low magnetic forces, computationally inexpensive models can be safely employed. Even though the experimental validation was not performed, similar studies proved that the fluid flow is perturbated under high magnetic fields and gradients when performing the continuous-flow magnetophoresis of ferrofluid droplets [85]. element and computational fluid mechanics software and predicted the separation performance of magnetic spores (Bacillus spores). Analogously, Kasetsirikul et al. [41] designed a device to separate malaria-infected RBCs, which show paramagnetism and suffer positive magnetophoresis when exposed to magnetic fields. Their numerical simulations suggest that the infected RBCs could be separated from healthy RBCs in a continuous flow using magnet arrays.
Numerical models have also been used to analyze the mixing of co-flowing fluids while magnetophoresis takes place. Gomez-Pastora et al. [32] employed a Lagrangian computational model for designing continuous-flow magnetophoretic microsystems for blood detoxification, assuring complete particle/biomaterial capture without compromising the quality of the blood (avoiding diffusion or loss of blood). The model was used to study the particle trajectory, the flow pattern and mass transfer between streams, and the analysis of multiple variables and parameters (flow rates, beads and magnets dimensions, fluid viscosities, etc.) was numerically and experimentally performed. In a later study [33], these authors provided a detailed analysis of the bead motion and its effect on fluid flow while magnetophoresis takes place ( Figure 6). Different scenarios were modeled and compared: (i) one-way coupling wherein momentum is transferred from the fluid to the beads, which are treated as point particles, (ii) two-way coupling wherein the beads are treated as point particles, and momentum is transferred from the bead to the fluid and vice versa, and (iii) two-way coupling taking into account the effects of bead volume in the fluid displacement. It was demonstrated that although there is little difference in the bead trajectories for the three scenarios, there is significant variation in the flow fields, especially when high magnetic forces are applied to the beads. Therefore, employing an accurate full flow-focused model that considers the effects of the bead motion and volume on the flow field when high magnetic forces are employed is recommended; however, when the beads are subjected to medium or low magnetic forces, computationally inexpensive models can be safely employed. Even though the experimental validation was not performed, similar studies proved that the fluid flow is perturbated under high magnetic fields and gradients when performing the continuous-flow magnetophoresis of ferrofluid droplets [85].

Dielectrophoretic Separation
The use of dielectrophoretic forces is very popular and one of the most widely used in the microfluidic separation of micron-sized materials. Table 2 presents the numerical

Dielectrophoretic Separation
The use of dielectrophoretic forces is very popular and one of the most widely used in the microfluidic separation of micron-sized materials. Table 2 presents the numerical studies published in the last 5 years dealing with the modeling and optimization of DEP in microdevices. As can be seen from the table, COMSOL Multiphysics is the preferred software to perform these studies. This software combines different modules (the laminar flow module, AC/DC module, and particle trajectory tracking module) valid to simulate the separation process [50,52]. More specifically, the laminar flow module is used to calculate the distribution of the flow field in the main channel; the AC/DC module is used to calculate the distribution of the electric field in the main channel; the particle trajectory tracking module coupling the influence of the flow field and electric field on the particle trajectory is used to calculate the trajectory of the particles in the channel [50]. These models have been experimentally validated in most of the studies, and good agreement between experiments and simulations has been achieved [49]. Table 2. Simulation studies dealing with particle and cell dielectrophoresis in microchannels.  No [69] The reported simulation studies have optimized the electrode configuration, applied electric frequency, voltage amplitude, fluid flow rate, conductivity, etc. [9,26,45,46]. It can be seen in Table 2 that most of the studies involving the separation of cells deal with blood separation and are performed using nDEP. For example, Nguyen et al. [26] carried out theoretical calculations and simulations to define the favorable parameters in the electric field operation of the microchip before its application to separate circulating tumor cells (CTCs) from RBCs, which is one of the most crucial procedures for subsequent detection and treatment of cancer.

Material to Be Separated
Similarly, Zhang and Chen performed a numerical study to find the optimum conditions (medium parameters, voltage, and frequency) for the DEP separation of blood cells [9]. In a different study [54], CTCs were separated from RBCs via DEP using microchannels with flow obstacles. A finite element simulation was conducted to predict the cell trajectories and revealed that the device with rectangular obstacles showed the best performance, achieving 100% separation efficiency and purity for cancerous cells in the specific outlets [54]. On the other hand, the optimization of electrode location and operating voltages were numerically performed by Kumar et al. [46]. These authors simulated several conditions and voltages for the separation of bacteria from blood cells and achieved effective separation of cells smaller than 1 µm using voltages below 20 V. It should be noted that high voltages are required to separate small particles via DEP. However, operating at higher voltages could lead to localized Joule heating, which may induce changes in conductivity and permittivity of the particles near the electrodes, and may create turbulence due to electrothermal fluid flow [46]. Furthermore, PDMS and polymethyl methacrylate (PMMA), such as deformable substrates used for microchannel manufacturing, may be affected by the induced heating at high voltages [46]. Thus, to overcome such electro thermal problems, optimization of electrode location and operating voltages need to be performed, and numerical models are powerful techniques to carry out such studies in advance of fabrication.
DEP has also been studied for the separation of particles [3,49,52]. Tirapu-Azpiroz et al. [49], using a FEM solver, simulated the dielectrophoretic separation of polystyrene (PE) beads using a capillary-driven microdevice. More specifically, the model was used to optimize the electrode layout and voltage configuration. These authors demonstrated the separation of 10 µm and 5 µm beads experimentally, with~98% efficiency. Chen et al. [52] employed three-dimensional (3D) electrodes composed of Ag powder (particle size of 10 nm) and PDMS for the separation of PE particles via DEP. Through theoretical calculations, numerical simulations, and experimental verification, the role of their new composite 3D microelectrodes was tested for the separation of 5, 10, and 20 µm PE beads.
Moreover, electric fields have been used to analyze different particles and cells according to their dielectric signatures. Differences in DEP responses of individual cells are based on differences in their structural and dielectric properties that might be linked to different intracellular processes. Fikar et al. [44] simulated and experimentally employed a novel method called distributed dielectrophoretic cytometry (2DEP cytometry) to measure dielectrophoresis forces acting on live cells in a microfluidic device. This method was employed to quantify the DE signatures of immortalized myelogenous leukemia cell lines. Analogously, Ettehad et al. [15] developed a microfluidic device that enabled the manipulation and characterization of live and dead yeast cells using dielectrophoretic forces. This is possible because of the different size and dielectric properties between live and dead cells. These authors, using FEM modeling, optimized the electrode configuration of the devices and selected the optimum operating conditions. The final device and their experimental approach successfully demonstrated the feasibility of the technology to trap live yeast cells inside the device and purify them from dead cells (Figure 7). nm) and PDMS for the separation of PE particles via DEP. Through theoretical calculations, numerical simulations, and experimental verification, the role of their new composite 3D microelectrodes was tested for the separation of 5, 10, and 20 µ m PE beads.
Moreover, electric fields have been used to analyze different particles and cells according to their dielectric signatures. Differences in DEP responses of individual cells are based on differences in their structural and dielectric properties that might be linked to different intracellular processes. Fikar et al. [44] simulated and experimentally employed a novel method called distributed dielectrophoretic cytometry (2DEP cytometry) to measure dielectrophoresis forces acting on live cells in a microfluidic device. This method was employed to quantify the DE signatures of immortalized myelogenous leukemia cell lines. Analogously, Ettehad et al. [15] developed a microfluidic device that enabled the manipulation and characterization of live and dead yeast cells using dielectrophoretic forces. This is possible because of the different size and dielectric properties between live and dead cells. These authors, using FEM modeling, optimized the electrode configuration of the devices and selected the optimum operating conditions. The final device and their experimental approach successfully demonstrated the feasibility of the technology to trap live yeast cells inside the device and purify them from dead cells (Figure 7). Finally, it should be noted that the nonuniform electric field can be generated using an array of electrodes or an array of insulating posts. Traditionally, DEP studies employed different microelectrode configurations, causing a nonuniform electrical field distribution inside the microchip. Advantages include a low voltage application for operating, but disadvantages exist, such as the need for complex microfabrication techniques and loss of functionality due to fouling [61]. The use of insulator-based DEP can overcome some of these disadvantages and has been the subject of numerical studies. For example, Aghaamoo et al. [61] used a numerical model to investigate the performance of different design aspects of insulating posts (pots geometry, electrode setups) and optimized the working conditions for the separation of CTCs from white blood cells (WBCs).

Acoustophoretic Separation
Different studies have addressed the simulation of particle and cell acoustophoretic separation with high accuracy. Table 3 presents the numerical studies that have been Finally, it should be noted that the nonuniform electric field can be generated using an array of electrodes or an array of insulating posts. Traditionally, DEP studies employed different microelectrode configurations, causing a nonuniform electrical field distribution inside the microchip. Advantages include a low voltage application for operating, but disadvantages exist, such as the need for complex microfabrication techniques and loss of functionality due to fouling [61]. The use of insulator-based DEP can overcome some of these disadvantages and has been the subject of numerical studies. For example, Aghaamoo et al. [61] used a numerical model to investigate the performance of different design aspects of insulating posts (pots geometry, electrode setups) and optimized the working conditions for the separation of CTCs from white blood cells (WBCs).

Acoustophoretic Separation
Different studies have addressed the simulation of particle and cell acoustophoretic separation with high accuracy. Table 3 presents the numerical studies that have been performed in the last 5 years to optimize the particle and cell acoustophoretic separation. As occurred with DEP, COMSOL Multiphysics is the software preferred since it integrates several modules that allow the accurate simulation of the process (pressure acoustics, laminar flow, particle tracing, etc.) [1,8,76]. More specifically, FEM codes have been used to optimize several variables and parameters that affect the process, such as the applied voltage or the effect of the IDTs location (distance from the channel and angle), number of IDT fingers, microchannel's height, the main flow velocity, and flow rate, etc. [8,19,71,74]. As seen in Table 3, the separation of cells via acoustophoresis, mostly blood cells, has been successfully modeled and optimized using numerical models. Moreover, acoustophoresis has been used not only to manipulate cells within microchannels but also to capture them in specific regions covered by antibodies. For example, Gupta et al. [16] placed antibody-coated polystyrene layers in the nodal regions of an acoustic field to capture infected RBCs with malaria parasites from plasma (the fluid condition was considered stationary to avoid cell lysis or cell-antibody bonds rupture due to shear stress). These layers have acoustic properties quite similar to the fluids involved (i.e., they are "acoustically transparent" as they do not impede wave propagation) and are biocompatible. Their device allowed bidirectional migration of the suspended cells to the biofunctionalized surfaces. Thus, simultaneously capturing infected erythrocytes on both layers could be feasible. The authors made use of FEM to model the pressure field as well as the motion of erythrocytes under the influence of acoustic radiation, drag, and gravitational forces. Based on the simulations and the trajectories of the erythrocytes, the location of the functionalized surface could be calculated before the fabrication of the device.
Furthermore, numerical models can also be implemented to account for the effect of channel elements (materials and orientations) in the SAW path. The total internal reflection (TIR) at the channel-fluid interface affects the entire fluid domain, where diffractive interference patterns arise from the imposition of a channel-bounded traveling SAW. TIR occurs when a wavefront propagates between domains with different sound speeds. It can cause undesired particle aggregation at channel walls and might require a complex design of channel or addition of sheath flows if not taken into account [72]. Since the walls are an essential component of these devices, it is important to couple particle actuation to the channel geometry/material rather than just the underlying SAW, allowing for highly localized patterning and focusing activities that can be incorporated by shaping the channel features [1,5]. For example, Richard et al. [1] used a numerical FEM model to study the effect of the wall material on the process (in this case, SU-8, which has a higher acoustic impedance than that of water, and this results in a more complex acoustic behavior than regular PDMS channels due to increased reflections at the fluid-solid interface). After ob-taining the numerical solution for the pressure field, the authors optimized the separation of RBCs from platelets in such a device (Figure 8).
design of channel or addition of sheath flows if not taken into account [72]. Since the walls are an essential component of these devices, it is important to couple particle actuation to the channel geometry/material rather than just the underlying SAW, allowing for highly localized patterning and focusing activities that can be incorporated by shaping the channel features [1,5]. For example, Richard et al. [1] used a numerical FEM model to study the effect of the wall material on the process (in this case, SU-8, which has a higher acoustic impedance than that of water, and this results in a more complex acoustic behavior than regular PDMS channels due to increased reflections at the fluid-solid interface). After obtaining the numerical solution for the pressure field, the authors optimized the separation of RBCs from platelets in such a device (Figure 8). Baek et al. performed a series of simulations accounting not only for the particle properties but also its surrounding component parameters [72]. To this end, they investigated the role played by the sound speed mismatch between the fluid and the channel in an SSAW device. More specifically, they determined the positions of patterned particles by analyzing the effects of the acoustic properties of the channel and the fluid on the pressure field, acoustic streaming field, and particle destination through both computational and experimental techniques [72]. On the other hand, Collins et al. [5] demonstrated (using an analytical model though) how curved and straight PDMS channel interfaces affect particle separation and tried to develop a comprehensive theory of channel wall interactions to examine the full range of channel wall orientations.
Finally, it should be noted that most theoretical/numerical studies on the acoustophoretic motion of a particle caused by acoustophoresis assume forces applied on small rigid particles (particle size <<< acoustic wavelength) suspended within an ideal fluid [70]. A few works have addressed the compressibility of the particle and used formulations beyond the ideal fluid assumption [70]. Studies dealing with the simulation of bigger particles (particle size ≈ wavelength) are scarce. For example, Samandari et al. [70] used a Baek et al. performed a series of simulations accounting not only for the particle properties but also its surrounding component parameters [72]. To this end, they investigated the role played by the sound speed mismatch between the fluid and the channel in an SSAW device. More specifically, they determined the positions of patterned particles by analyzing the effects of the acoustic properties of the channel and the fluid on the pressure field, acoustic streaming field, and particle destination through both computational and experimental techniques [72]. On the other hand, Collins et al. [5] demonstrated (using an analytical model though) how curved and straight PDMS channel interfaces affect particle separation and tried to develop a comprehensive theory of channel wall interactions to examine the full range of channel wall orientations.
Finally, it should be noted that most theoretical/numerical studies on the acoustophoretic motion of a particle caused by acoustophoresis assume forces applied on small rigid particles (particle size <<< acoustic wavelength) suspended within an ideal fluid [70]. A few works have addressed the compressibility of the particle and used formulations beyond the ideal fluid assumption [70]. Studies dealing with the simulation of bigger particles (particle size ≈ wavelength) are scarce. For example, Samandari et al. [70] used a FEM simulation approach to simulate with high resolution the acoustophoresis of particles at high frequencies in a microchannel. More specifically, they used a FEM model to amend analytical calculations of the acoustic radiation force arising from an imposed standing ultrasound field. To calculate the force and the resultant deformation-induced on the particles, they implemented a solid acoustic interaction (ASI) approach. A particle-tracking scheme was used to obtain the separation efficiency and the patterning of the particles under relevant conditions. They demonstrated a significant mismatch between the force obtained from the ASI approach and previous analytical predictions that considered small particle assumptions when working at high manipulation frequencies due to both different ARF values and the variation in force maps in multidirectional wave propagation.

Conclusions and Further Directions
In this work, we have summarized the numerical modeling of the active separation of particles and cells within microchannels. Specifically, separation using magnetophoretic, dielectrophoretic, and acoustophoretic forces has been analyzed. First, we have established the different approaches used to calculate the particle/cell migration inside the device, emphasizing the Lagrangian approach, where the particle trajectory is calculated as a function of the external forces acting on the material. The formulation of the drag force, gravitational force, and the active forces responsible for the separation are also formulated. For active forces, equations to calculate the particle movement in magnetophoresis, dielectrophoresis, and acoustophoresis are developed. Different methods to solve the equations are also introduced, such as the FEM and FVM, and the importance of performing a mesh dependency study is highlighted. Finally, the different studies dealing with the numerical modeling of the process are presented. We demonstrated the effectiveness of numerical modeling in optimizing the different variables and parameters that affect the process in advance of fabrication. Also, we reviewed the models developed in the different studies, the software employed, and the conditions simulated to offer the readers different guidelines about the capabilities of the CFD software. Nevertheless, there are different directions where improvement can be done.
On the one hand, and as seen in the sections above, COMSOL Multiphysics ® (versions 5 and above, founded by Dr. Svante Littmarck and Farhad Saeidi, 100 District Avenue, Burlington, MA 01903, USA)is the preferred software to simulate the process. This multipurpose software package is user-friendly and contains the modules required to simulate the separation. It is specifically designed to solve multi-applied physics and engineering phenomena as a simulation tool for electrical, mechanical, fluid flow, chemical, and other applications within one package [39,40]. Other software platforms have also been used, but they require the implementation of the fields and forces by the user using customized codes, which limits their use to more experienced CFD engineers. For example, DEP has been modeled using Ansys Fluent, but a User Defined Function (UDF) needs to be written in C++ and integrated into the Ansys Workbench to solve for the electric field and electric potential distribution within the computational domain and also to calculate the DEP forces at any location depending on the fields and particle/cell parameters [63]. Analogously, for modeling magnetophoresis, some authors have used Ansys Fluent or Flow-3D [32,33], which require the implementation of the magnetic field, gradient, and force by using external subroutines linked to the main solver. Future efforts should be directed towards implementing modules that allow the calculation of electric/magnetic fields and gradients and the calculation of forces acting on the particulate material.
On the other hand, most of the studies have neglected particle-particle and particlefluid interactions. It should be noted that a precise calculation of the interparticle effects as well as the effects that the particles have on the fluid field might be required for some applications. For example, magnetic dipole-dipole interactions or electrostatic interactions between particulate material could be beneficial for the separation process as they may speed up the aggregation of the material and/or magnify the force acting on it [82]. On the other hand, modifying the flow patterns as the material is separated might have detrimental consequences for the process if mixing is undesired and could underestimate the value of the force required for separation [33]. These effects are not easily implemented in commercial CFD software programs. Therefore, future work should focus on improving the accuracy of these models to apply them to the processes where particle-particle and particle-fluid effects cannot be neglected.
It should also be noted that although the experimental validation is carried out in most studies, it is not properly quantified. The difference between the experimental and theoretical results (for example, the difference in particle position values or the particle recovery) is only presented in a few works. Thus, future efforts should also be directed toward quantifying the accuracy of the theoretical models by performing proper experimentation.
In conclusion, the potential of CFD software to accurately model microscopic and macroscopic systems while integrating a plethora of different physical models into one simulation, and along with advances in computer hardware and algorithms, only favor numerical modeling. We believe that the future of this technology is bright, and should prove paramount in the development and optimization of the design of novel cell and particle separation medical devices.