ABSTRACT
Many important microfluidic applications require the control and transport of particles immersed in a fluid. We propose a model for automatically planning good trajectories from an arbitrary point to a target in the presence of obstacles. It can be used for the manipulation of particles using actuators of mechanical or electrical type. We present the mathematical formulation of the model and a numerical method based on the optimization of travel time through the Bellman’s principle. The implementation is focused on square grids such as those built from pixelated images. Numerical simulations show that the trajectory tree produced by the algorithm successfully avoids obstacles and stagnant regions of the fluid domain.
Keywords:
Trajectory planning; manipulation of particles in microfluidics; Bellman’s principle
RESUMO
Importantes aplicações em microfluídica requerem o controle e transporte de partículas imersas em um fluido. Propomos um modelo para o planejamento automático de boas trajetórias de um ponto arbitrário até um alvo, na presença de obstáculos. Esse modelo pode ser utilizado para a manipulação de partículas fazendo uso de atuadores do tipo mecânico ou elétrico. Apresentamos a formulação matemática do modelo e um método numérico baseado na otimização do tempo de chegada através do princípio de Bellman. A implementação considera malhas quadradas, as quais são construídas a partir de imagens pixeladas. Simulações numéricas mostram que a árvore de trajetórias obtida do algoritmo evita obstáculos e regiões estagnadas do domínio fluídico com sucesso.
Palavras-chave:
Planejamento de trajetória; manipulação de partículas em microfluídica; Princípio de Bellman
1 INTRODUCTION
Particle manipulation is an important topic in microfluidics, a science that studies the behavior of small-scale fluids. A growing interest exists in the steering of particles in microfluidic systems, especially in the biological and medical fields, to analize elements of very small size such as proteins or DNA. In particular, improving the ability to manipulate and to steer particles in microfluidics can allow to better examine, control and treat biological materials. Some interesting examples are the microfluidic device presented in 1313 Y. Zhou, Y. Wang & Q. Lin. A microfluidic device for continuous-flow magnetically controlled capture and isolation of microparticles. Journal of Microelectromechanical Systems, 19(4) (2010), 743-751. or the use of optical tweezers as tools to move in a controlled way individual particles, as shown in 66 J.E. Curtis, B.A. Koss & D.G. Grier. Dynamic holographic optical tweezers. Optics communications, 207(1) (2002), 169-175..
The control of the trajectory of particles immersed in a fluid by boundary actuators is the topic presented in the articles by Chaudhary and Shapiro 44 S. Chaudhary & B. Shapiro. Arbitrary steering of multiple particles independently in an electroosmotically driven microfluidic system. IEEE Transactions on Control Systems Technology, 14(4) (2006), 669-680. and Tuval et. al. 1212 I. Tuval, I. Mezić, F. Bottausci, Y.T. Zhang, N.C. MacDonald & O. Piro. Control of particles in microelectrode devices. Physical review letters, 95(23) (2005), 1-4.. The microfluidic modeling in these works proposes the displacement of the particles from one point to another by a pre-established trajectory. The idea is to identify the location of the particle in real time so that a control algorithm calculates the demand on each actuator for the flow to carry the particle along the desired path. In Armani et. al., 11 M. Armani, S. Chaudhary, R. Probst, S. Walker & B. Shapiro. Control of microfluidic systems: two examples, results, and challenges. International Journal of Robust and Nonlinear Control, 15(16) (2005), 785-803. and 22 M. Armani, Z. Cummins, J. Gong, P. Mathai, R. Probst, C. Ropp, E. Waks, S. Walker & B. Shapiro. Feedback Control of Microflows. In “Feedback Control of MEMS to Atoms”. Springer (2012), pp. 269-319., some examples can be found of the use of control techniques. A generalization for the manipulation of particles in the three-dimensional case also exists in the study of Probst and Shapiro 1010 R. Probst & B. Shapiro. Three-dimensional electrokinetic tweezing: device design, modeling, and control algorithms. Journal of Micromechanics and Microengineering, 21(2) (2011), 027004..
In all published methods discussed above, a predefined trajectory is necessary as a datum but no mention is given about how to build such a trajectory. In fact, in their examples the trajectory is introduced manually. The present work proposes a method to automatically define a “good” trajectory γ* that can then be used as input for those particle manipulation techniques. By “good” we intend that γ* not only avoids obstacles, but also avoids parts of the fluid domain where the actuators are ineffective at producing flow (stagnant regions, or regions of low controllability). These requirements are fulfilled by time-optimal trajectories, which are the ones proposed here. In fact, we propose an algorithm that works on a rectangular grid and could be implemented directly on a pixelated image of the domain. The trajectories are discrete in the sense that they are broken lines joining neighbor grid points. This is a very practical, though not smooth, parameterization of possible paths. The resulting trajectory γ* can eventually be smoothed out before passing it to the manipulation algorithm. The combination of the proposed method with the manipulation techniques of Shapiro and coworkers (44 S. Chaudhary & B. Shapiro. Arbitrary steering of multiple particles independently in an electroosmotically driven microfluidic system. IEEE Transactions on Control Systems Technology, 14(4) (2006), 669-680.), (11 M. Armani, S. Chaudhary, R. Probst, S. Walker & B. Shapiro. Control of microfluidic systems: two examples, results, and challenges. International Journal of Robust and Nonlinear Control, 15(16) (2005), 785-803.), (22 M. Armani, Z. Cummins, J. Gong, P. Mathai, R. Probst, C. Ropp, E. Waks, S. Walker & B. Shapiro. Feedback Control of Microflows. In “Feedback Control of MEMS to Atoms”. Springer (2012), pp. 269-319. and 1010 R. Probst & B. Shapiro. Three-dimensional electrokinetic tweezing: device design, modeling, and control algorithms. Journal of Micromechanics and Microengineering, 21(2) (2011), 027004.) would allow for an integrated system of detection, planned manipulation and control of particles towards an arbitrary target.
The presentation begins by describing the general physical model and recalling the techniques used for particle steering when the desired trajectory is known. This allows us to introduce several concepts needed for trajectory planning. Then the continuous problem of optimizing trajectories (Bellman’s principle) with respect to arrival time is described, together with the proposed discretization. Finally, applications are reported in the section of numerical results. They show the good performance of the algorithm in finding suitable, sometimes counterintuitive, trajectories.
2 THE PHYSICAL MODEL
We study the manipulation of a particle in a microfluidic device in the two-dimensional case, the extension to 3D being immediate. In fact, we can see the 2D case applicable to flows between flat plates parallel to the plane xy, in the case when the plates are sufficiently close.
Let us consider therefore a square domain , for the flow of a fluid without inertia with the presence of obstacles, walls and also of actuators. As we can see in the Figure 1, the actuators act in some regions of the boundary of the domain ΓN , ΓS , ΓE and ΓW , being placed on the respective north, south, east and west borders.
At each actuator we apply a value X N , X S , X E e X W . Since the driving force of the flow is the difference between actuator values (pressure or electrical potential differences), we can fix X E to zero without losing generality and work only with the variables . Inertial effects are generally negligible in microfluidics. As a consequence, each (where 𝒜 is the set of admissible actuator values) leads to a velocity field according to the following linear relation
The matrix contains as columns the components of the velocity field at the point with X belonging to the canonical basis, i.e. and for the first column; and for the second one and and for the third one. The Figure 2 shows an example of the velocity field, for some points of the domain, generated by imposing the canonical basis.
Velocity fields generated by imposing the canonical basis. All possible velocity fields are linear combinations of these three.
The admissible set of the controls 𝒜 generates an admissible set of velocities on each point, defined by
Let P be a passive point particle whose initial position is . Each control function moves P along a trajectory that is the solution of the following problem
The goal of this work is to answer the question of how to determine the control function that moves the particle P from the position to a target set . The target set Z may be a region of the domain where the particle is analyzed or subjected to some biochemical reaction, for example. In order to completely achieve this aim we have first to solve the problem of the particle manipulation through a predefined path. As explained in the following subsection, this approach introduces an optimization problem for maximizing the possible velocity. The solution of this preliminary problem allow us to have all the tools for achieving the main goal, which is to determine a “good” trajectory that consider the arrival time optimization.
2.1 Particle manipulation along a predefined path
Let us consider the determination of the control function when the path to move the particle P from to Z is given. Let γ be the chosen path that is a smooth curve in Ω joining to . For every point , let be the unit vector tangent to γ at . In order to move P along γ it is necessary and sufficient that is parallel to for all t. Since the relation (2.1) does not depend explicitly on t, time becomes a parameter and will be omitted in what follows for conciseness. We have the following characterization:
Proposition 1. At a point the control X generates a velocity vector parallel to an arbitrary direction , if and only if
where d is the (column) array representation of. Beingof full rank, and the number of controls equal to or greater than the number of spatial dimensions, the above equation has as its solution a non-trivial subpaceof ℝn .
In particular, the maximum velocity that can be generated along direction is
and the maximizing argument is the control value that realizes the velocity .
If the admissible set of controls is a box with minimum X - and maximum X + values for each control variable, the following classical linear optimization problem results:
Linear optimization problem: Maximize c T X over (being ), with the restrictions
The maximum reached is the maximum velocity , and the maximizing argument is .
With these elements, bringing back the time variable, the problem (2.2) becomes
which will necessarily follow the curve γ at maximum speed, and has been calculated the control function arises from
The control at every instant is uniquely determined by the position of the particle and the direction in which it is desired to move it. The problem of microfluidically controlling a particle to follow a predefined trajectory was addressed by Chaudhary and Shapiro, 44 S. Chaudhary & B. Shapiro. Arbitrary steering of multiple particles independently in an electroosmotically driven microfluidic system. IEEE Transactions on Control Systems Technology, 14(4) (2006), 669-680. and the work of Mathai et al. 88 P.P. Mathai, A.J. Berglund, J.A. Liddle & B.A. Shapiro. Simultaneous positioning and orientation of a single nano-object by flow control: theory and simulations. New Journal of Physics, 13(1) (2011), 013027. , where the simplest constraints (, where ║ ∙ ║ is the Euclidean norm) are considered.
Remark 2.1. In (2.4) the vector can be multiplied by an arbitrary if the desired speed of the particle is not the maximum attainable one at each point. For example, let
If , it means that it is impossible to reach the target Z in finite time along the path γ. Otherwise , for any the control function that moves P along γ at constant speed V* is obtained by first solving
and then setting
2.2 Determining the path via arrival time optimization
When the path is not pre-defined, it can be chosen in several ways. Here we study the trajectory planning between and the Z target by minimizing the time of arrival. Other possibilities include minimizing path length, energy expenditure, and in general any cost function depending on the path and possibly on the associated control function.
For a point (controllable, i.e., such that there exists at least one function X(t) which moves P from to Z) we can consider the function defined as the minimum time required to take P from to Z. The path associated with is an optimal trajectory because it moves P to Z in minimum time and defines, at each point, the optimal direction (tangent to the optimal trajectory).
Bellman’s principle states that if a point is placed in the optimal path γ* of , then the optimal trajectory from to is the arc of γ* joining these two points, and the optimal trajectory of to the target Z is the final arc of γ*. With this principle it is possible to obtain the Hamilton-JacobiBellman equation for T , according to 33 S. Cacace, E. Cristiani & M. Falcone. Can Local Single-Pass Methods Solve Any Stationary Hamilton-Jacobi-Bellman Equation? SIAM Journal on Scientific Computing, 36(2) (2014), A570-A587.:
where is given by . We have that the associated Hamiltonian is , which gives the equation in the usual form (see 1111 J.A. Sethian & A. Vladimirsky. Ordered upwind methods for static Hamilton-Jacobi equations. Proceedings of the National Academy of Sciences, 98(20) (2001), 11069-11074.)
The fundamental property we will use in the numerical treatment of this equation is the following (a direct consequence of Bellman’s principle):
Proposition 2. Let , and C a closed curve (surface in 3D) such that is internal to C and the target Z is external to C. For each , let be the minimum time required to carry a particle from to . Then,
Moreover, if is a point of C where the minimum is attained, is the intersection of the optimal path γ* with the curve C.
If the curve C is small enough (contained in a ball of ratio ), we can consider just straight trajectories between and C, which allows to approximate the equation (2.5), according to 33 S. Cacace, E. Cristiani & M. Falcone. Can Local Single-Pass Methods Solve Any Stationary Hamilton-Jacobi-Bellman Equation? SIAM Journal on Scientific Computing, 36(2) (2014), A570-A587., by
and the optimal trajectory, locally, is given by the segment .
3 CALCULATION OF
To find the set of possible maximum velocities at a point , for each given direction , we have first to determine the admissible set of velocities , according to the equation (2.1). We consider that the system actuators can be mechanical (imposing pressure) or electrical (imposing electrode voltage).
When the velocity field is generated by pressure sources, the governing equation for the flow is the Stokes equation, that in our setting defines the following boundary value problem
where is the stress tensor and the unit normal to the boundary. We remark that the boundary conditions for this problem are essentially that pressure is imposed on the actuators and null velocity is imposed on the obstacles and walls.
When the system actuators are electrodes, the electric field is obtained from the electric potential Φ (or voltage) that solves the following problem
The electric field induces a velocity field according to the relation
where m eo is the electrosmotic constant relative to the interaction between the fluid and the wall 77 B.J. Kirby. “Micro-and nanoscale fluid mechanics: transport in microfluidic devices”. Cambridge university press (2010).. In this case, Dirichlet boundary conditions for the voltage are imposed on the actuators, and Neumann homogeneous ones on the obstacles and walls.
A detailed explanation of the numerical resolution of these problems, according to MAC (Marker and Cell) discretization, is exposed in 99 L. Meacci, F.F. Rocha, A.A. Silva, P.V. Pramiu & C.G. Buscaglia. Planejamento de trajetória para a manipulação de partículas em microfluídica. CQD-Revista Eletrônica Paulista de Matemática, 10 (2017), 19-37..
Remark 3.1. In fact, one could easily consider a combination of mechanical and electrical actuators. In this case, if a control variable X i corresponds to a mechanical actuator, the i-th column of matrix B in (2.1) would be the velocity field obtained by solving (3.1), while if it corresponds to an electrical actuator the equation to be solved would be (3.2). We have not considered problems with both kinds of actuators simultaneously to keep the exposition simpler.
The calculation of involves a resolution of the linear optimization problem of the Proposition 1 on all mesh nodes for a directions set D defined by the user.
If v1, v2 and v3 are the three vectors (column) calculated with the hydrodynamic solver on the node , and d is a direction belonging to D, the following Octave code calculates :
It is important to justify that the use of the negative sign in is because by pattern the function glpk minimizes instead of maximizing. We remark that the values X- and X+ are in the vectors Xmin and Xmax. The value of the control to generate the velocity V in the direction d is given by xopt.
The Figure 3 (a) illustrates an example of the set obtained initially with the vectors
In this Figure 3 (a), the external polygon was generated using as minimum and maximum values the vectors
When we change the vector of minimum values of the control variables to
we obtain a new set , which corresponds to the internal polygon. Notice how the polygon that represents loses its symmetry about the origin when Xmin is different from -Xmax. We remark that the number of directions used to determine the region can influence the quality of the set of possible velocities generated. The Figure 3 (b) illustrates a set given by the values
where the internal polygon was generated using 8 directions and the external polygon was generated using 320 directions.
By performing the same procedure for all nodes of the mesh, we obtain the data necessary to calculate the optimal trajectory according to the Hamilton-Jacobi-Bellman equation. The Figure 4 shows the set for some points of the domain, where in the Figure 4 (a) the actuators are of the pressure type, and in the Figure 4 (b) the actuators are of the electrode type. We can see the variability of , which is highly anisotropic and highly dependent on the geometric configuration as well as on the nature of the actuators. The regions where the polygons are smaller correspond to stagnant regions, i.e., regions where the velocity generated by the actuators hardly moves the particles at all. The regions where the polygons are bigger correspond to regions where the particles are easily transported. A “good” trajectory for manipulation avoids not just the obstacles but also these stagnant regions as much as possible.
Representation of the set of the possible maximum velocities (in each of the 32 directions considered) in some points of the domain.
4 OPTIMIZING TRAJECTORIES BY BELLMAN’S PRINCIPLE
The optimal trajectory, in terms of time, for moving a particle from a point to a target point, solves the Hamilton-Jacobi-Bellman equation. A semi-Lagrangian discretization is considered for the function T. The equation (2.6) is approximated at each node (i, j), taking the curve C (denoted as C i,j ) as the boundary of the set of the four cells that share node (i, j). The coordinates of node (i, j) are (see Figure 5). For each direction , the trajectory from to C i,j is approached by a straight line from to .
The values for the discrete approximation of T i,j are defined by
where is an interpolation of the T values at the nodes that are placed on C i,j . In the case of the direction illustrated in the Figure 5, for example, is calculated by linearly interpolating between T i+1,j e T i+1,j+1 . The direction of the optimal trajectory is that in which the minimum is attained, which is generally not parallel to ∇T because of the anisotropy of the problem.
To leave a completely discrete scheme we consider the following discretization of the set of directions D
The algorithm requires that has already been calculated at each vertex and has been stored as a table of values , corresponding to each vertex (i, j) and to each direction . T i,j is initialized to a large value (1012) on all nodes except those that belong to the target Z, which are initialized to zero. Then we iterate all the mesh nodes that do not belong to Z updating T according with the following equation
This algorithm is discussed in 55 E. Cristiani. “Numerical methods for optimal control problems: part II: local single-pass methods for stationary HJ equations”. IAC-CNR, Roma (2013). p. 21. Class notes., including error estimators for the solution in the norm of L ∞(Ω), which is important because the exact solution in general is not smooth and can not be understood in the classical sense.
Remark 4.1 (Complexity). Calculations necessary to start the previous algorithm are:
-
Computing the matrix
according to the linear relation (2.1) for the N nodes and n actuators with cost n𝒪(N a ), where 𝒪(N a ) is the cost of solving one hydrodynamic problem (tipically a ≈ from 1.3 to 1.5); -
Computing
with cost mN𝒪(n b ), obtained by the cost 𝒪(n b ) of solving one linear optimization problem (2.3) in dimension n (tipically b ≈ 2), for all N nodes and all m directions.
From (4.1) the cost of the HJB iteration is 𝒪(mN), and the number of required iterations is 𝒪(N 1/s ), s being the number of space dimensions. Summarizing, the total cost adds up to . Since , in general the HJB iterations dominate the total cost, but if the number of actuators n is large the cost of the linear optimization problems can also be significant.
5 RESULTS
In this section, we consider the set D generated by 8 directions. The accuracy of the calculation depends on the number of the considered directions, and so does the computational cost. The choice of 8 directions has the advantage of not requiring the interpolation step, since the directions considered in (4.1), with , are determined by the neighboring nodes. Therefore for each node (i, j) the optimal trajectory is a line that connects it to the neighboring mesh node at which the minimum in (4.1) is attained. As a consequence, the set of numerical trajectories produced by our algorithm has a tree structure. For this reason, to plot the numerical solution we draw all the optimal trajectories that pass through the mesh nodes, which has the geometrical appearance of a tree rooted at the target Z.
In the next simulations, we consider a domain discretized in a grid of 40 × 40 nodes, except when indicated otherwise. When the flux is induced by an electrical potential we assume an electrosmotic costant while when it is generated by a pressure source we consider viscosity . The actuators’ limit values are set to .
Figure 6 shows the trees of optimal trajectories obtained for both electrical (subfigure (a)) and mechanical (subfigure (b)) actuators.
Optimal trajectory configuration with range of actuators between −1 and 1. The destacked point is the arrival target.
One can notice that the trajectories between the two plates of the obstacle are essentially horizontal. This is because the actuators cannot generate significant vertical velocity there due to the presence of the plates. With this, it turns out that a particle that is between the plates but below or above Z must first take a path to get out of the space between the plates and then walk back in towards the target point. We find worthy of note that the optimal trajectories tend to circumvent the obstacle so as to traverse regions (more distant to the obstacle) where the admissible velocities are higher. Notice also that the trajectory trees have similarities for both types of actuators, mainly differing at the obstacles’ boundaries, where the mechanically-driven flow has zero velocity while the electrically-driven one only imposes the normal component of the velocity to be zero.
The model exhibits convergent behavior with respect to mesh refinement in all cases considered. As an example, in Figure 7 we compare the resulting trees for the electrically-driven model, with meshes 40 × 40 and 80 × 80.
Effect of mesh refinement on the resulting trajectory trees, for the electrically-driven flow.
The optimal trajectories may be counterintuitive. Let us consider the flow generated by pressure actuators setting as target Z a single point close to the obstacle as shown in Figure 8. Then we consider two possible start points, labeled with the letters A and B. In the case with A as starting point, the model designs an optimal trajectory that is a straight segment. The travel time t computed by integrating the field between A and Z along a straight line coincides with , where T is the solution of the discrete Hamilton-Jacobi-Bellman equation (4.1). On the other hand, when the starting point is B the optimal trajectory proposed by the model is not straight. It assumes a form that is not of immediate intuition. The designed path, although it is longer than the straight path, results in a shorter travel time t op = 6.7 than that obtained along the straight one t = 9.1. This advantage is due to the greater efficiency of the actuators in the region traveled by the optimal trajectory. The values to be applied to the actuators to manipulate the particle along the proposed path are easily obtained from (2.4) and are shown in Figure 9 as functions of time for the case that has B as starting point. As usual in time-optimal control, the actuator values exhibit a bang-bang-like behavior. The geometrical path is the most fundamental result since it avoids obstacles and stagnant regions. This path can then be traveled at non-optimal velocities, with smoother control functions, as discussed in Remark 2.1.
6 CONCLUSIONS
We have presented an integrated model to identify convenient trajectories for the manipulation and control of particles in microfluidics. The achievement of the above-mentioned objective involved various mathematical problems:
-
• To find an admissible set of velocities
generated by a set of admissible actuator values 𝒜 , solving -
- a Stokes problem for the velocity field induced by pressure sources;
-
- a Poisson problem for the electrical potential induced by a voltage imposed at electrodes, which in turn induces an electrosmotic velocity field;
-
• To determine the maximum admissible velocity
, through the resolution of a linear optimization problem; -
• To provide the control functions
arising from the resolution of an initial value problem for the trajectory ODE; -
• To plan the trajectory by minimizing the arrival time through the numerical implementation of Bellman’s rinciple.
The proposed method automatically indentifies a trajectory for a particle that we must manipulate from a departure point to the target in minimum time (and also provides the corresponding values to impose on the actuators). The numerical results show that the trajectories found not only avoid the obstacles, but also avoid the regions of the domain where the actuators are ineffective at producing flow. The numerical simulations show counterintuitive optimal trajectories in this kind of problem, justifying the use of mathematical models for designing automatic manipulation systems.
It should be noted that the proposed algorithm produces not just one optimal trajectory, but the whole tree of optimal paths, starting at any . This apparent overkill is very important when applying the model under realistic conditions, in particular considering external perturbations (noise). In fact, if the sensor in the control loop detects that a perturbation has taken the particle to a position outside the original optimal path, the controller can immediately switch to the path in the optimal tree that passes throught . In this way, and without any additional computation, the optimal tree can be used to correct perturbations so that the particle eventually reaches Z.
ACKNOWLEDGEMENTS
This work was carried out with the financial support of the CAPES and CNPq agencies. In particular, Luca Meacci acknowledges for the financial support the Programa de Excelência Acadêmica (Proex) - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - Brazil (CAPES) - Finance Code PROEX-9740044/D.
REFERENCES
-
1M. Armani, S. Chaudhary, R. Probst, S. Walker & B. Shapiro. Control of microfluidic systems: two examples, results, and challenges. International Journal of Robust and Nonlinear Control, 15(16) (2005), 785-803.
-
2M. Armani, Z. Cummins, J. Gong, P. Mathai, R. Probst, C. Ropp, E. Waks, S. Walker & B. Shapiro. Feedback Control of Microflows. In “Feedback Control of MEMS to Atoms”. Springer (2012), pp. 269-319.
-
3S. Cacace, E. Cristiani & M. Falcone. Can Local Single-Pass Methods Solve Any Stationary Hamilton-Jacobi-Bellman Equation? SIAM Journal on Scientific Computing, 36(2) (2014), A570-A587.
-
4S. Chaudhary & B. Shapiro. Arbitrary steering of multiple particles independently in an electroosmotically driven microfluidic system. IEEE Transactions on Control Systems Technology, 14(4) (2006), 669-680.
-
5E. Cristiani. “Numerical methods for optimal control problems: part II: local single-pass methods for stationary HJ equations”. IAC-CNR, Roma (2013). p. 21. Class notes.
-
6J.E. Curtis, B.A. Koss & D.G. Grier. Dynamic holographic optical tweezers. Optics communications, 207(1) (2002), 169-175.
-
7B.J. Kirby. “Micro-and nanoscale fluid mechanics: transport in microfluidic devices”. Cambridge university press (2010).
-
8P.P. Mathai, A.J. Berglund, J.A. Liddle & B.A. Shapiro. Simultaneous positioning and orientation of a single nano-object by flow control: theory and simulations. New Journal of Physics, 13(1) (2011), 013027.
-
9L. Meacci, F.F. Rocha, A.A. Silva, P.V. Pramiu & C.G. Buscaglia. Planejamento de trajetória para a manipulação de partículas em microfluídica. CQD-Revista Eletrônica Paulista de Matemática, 10 (2017), 19-37.
-
10R. Probst & B. Shapiro. Three-dimensional electrokinetic tweezing: device design, modeling, and control algorithms. Journal of Micromechanics and Microengineering, 21(2) (2011), 027004.
-
11J.A. Sethian & A. Vladimirsky. Ordered upwind methods for static Hamilton-Jacobi equations. Proceedings of the National Academy of Sciences, 98(20) (2001), 11069-11074.
-
12I. Tuval, I. Mezić, F. Bottausci, Y.T. Zhang, N.C. MacDonald & O. Piro. Control of particles in microelectrode devices. Physical review letters, 95(23) (2005), 1-4.
-
13Y. Zhou, Y. Wang & Q. Lin. A microfluidic device for continuous-flow magnetically controlled capture and isolation of microparticles. Journal of Microelectromechanical Systems, 19(4) (2010), 743-751.
Publication Dates
-
Publication in this collection
Sep-Dec 2018
History
-
Received
03 Oct 2017 -
Accepted
30 May 2018