This invention was made with support of the United States Government. The United States Government has certain rights in the invention.
1. Field of the Invention
The invention described herein is directed to the manipulation of multiple objects suspended in fluids within microfluidic devices through the application of complex force fields. More specifically, the present invention is a method for determining and subsequently applying a set of signals to one or more actuators on a microfluidic device so as to respectively apply a corresponding force on each of one or more objects contained therein to thereby manipulate the position, velocity, shape, orientation, and/or distribution thereof.
2. Description of the Prior Art
Microfabrication techniques have been used for over a decade to produce a variety of submillimeter mechanical structures. For example, the new fabrication techniques have led the way to the production of Micro-Electro-Mechanical Systems (MEMS) in which microscopic machinery, sensors, actuators, and electronic circuitry are assembled on, and in many cases etched from, a common substrate, such as silicon. Many of these micromachined devices have enjoyed a wide range of applicability in such fields as chemical and biological research.
In certain technological fields, such as the aforementioned chemical and biological research, the physical scale of certain domains of interest have motivated the development of sophisticated equipment capable of manipulating microscopic objects, both individually and in selected groups. One device widely used in this area is the optical tweezer (also known as laser tweezers) which uses laser light to manipulate objects of molecular size scales. Optical tweezers create an optical trap on an object through light scattering forces and light intensity gradient forces of a focused laser beam. The forces combine to hold the object in the center of the focused laser spot. The trapped object may be repositioned by moving the focused laser spot as desired.
Optical tweezers are effective in manipulating certain types of objects, but suffer several shortcomings which prevent their implementation in a wider range of applications. Optical tweezers can be used to manipulate particles, provided there is a difference in the index of refraction of the particle and that of the surrounding medium, but have yet to be used effectively to redistribute fluids. Manipulation of fluids on a microscopic scale has become a useful method in delivering chemical agents to individual cells to observe how the cells react thereto.
Another shortcoming of optical tweezers is that they are large and generally expensive pieces of equipment. In typical applications, a laser tweezer will consist of one or more lasers, a microscope, and high-quality focusing optics to produce each optical trap.
Manipulation of submillimeter objects by electrophoresis and other methods using an applied electric field have been used for many years. Electrophoresis has been widely used for separating and sorting particles into bands in accordance with particle size and inherent electric charge. Gel electrophoresis, for example, whereby an electric field is applied to molecules suspended in a porous gel, is used in the field of genetics for DNA profiling. However, while electrophoresis is useful for applying a force on certain particles, such as molecules, the process is not operative on objects immovable by an electric field, such as objects made of a dielectric material. Additionally, electrophoresis may be used to sort particles according to a charge/size ratio into sorting bins located along a straight path, but does not provide a means for steering the objects toward alternative locations not on the path.
Sorting of minute particles is a prevalent requirement in many research and biochemical fields, and many means for performing this task are widely available. For example, one common application passes a stream of particles suspended in an electrolyte through a small aperture over which an electric field is applied. A particle in the aperture displaces an amount of electrolyte equal to its own volume. In accordance with the Coulter principle, the volume displaced changes the impedance of the aperture and is measured as a voltage pulse, the height of which is proportional to the volume of electrolyte displaced, i.e., the volume of the particle. The particles may then be sorted by size by deflecting different sized particles into a corresponding sorting location or bin, by some mechanism such as an optical tweezer. The Coulter particle sorter illustrates the benefits of device implementation on a micromachined platform, i.e., to measure a change in impedance in a channel or aperture caused by the presence of a microscopic particle, the channel or aperture is required to be formed on a size scale comparable with the size of the particle.
An illustrative example of another cell sorting device constructed by micromachined techniques is provided by the journal article, “An Integrated Microfabricated Cell Sorter,” by Anne Fu, et al. (Analytical Chemistry, Vol. 74, No. 11, Jun. 1, 2002). The referenced cell sorter implements a network of microvalves and micropumps for controlling the movement of cells suspended in a fluid after the classification thereof by controlling the surrounding fluid flow. Cells within the device are classified by means of fluorescence of the cell resulting from excitation by an argon laser. Various valves and pumps are activated in accordance with one of a number of predetermined patterns so as to direct a particle to a destination sorting bin by directing the flow of the suspending fluid along one of a number of predetermined paths. However, the geometry of the fluid channel and the pump and valve configuration allow only a limited control over the motion of any particular cell. Moreover, the configuration does not afford simultaneous parallel control of multiple particles within the fluid. For example, the device does not contemplate directing different objects toward each other.
Particle placement and sorting are not the only tasks for which microscopic object manipulation means are desired. Many applications require the manipulation of fluids on a microscopic scale for purposes of, for example, mixing, dosing, and delivering small quantities of drug to individual cells. In other applications, objects such as strands require shape orientation or conformation. For example, in certain applications, DNA strands may need to be “unwrapped” to expose certain structural features for study. To perform these functions, complex motion of multiple particles, strand segments and surfaces is necessary. However, simultaneous arbitrary control of the trajectories of multiple objects presents a challenging controller design problem.
One device for controlling the motion of multiple objects is the Universal Planar Manipulator (UPM) developed by Dan Reznik, formerly of the University of California at Berkeley. Objects to be manipulated are placed on a rigid, horizontally oriented plate. The plate is coupled to one or more actuators which vibrates the plate in the horizontal plane. The objects are moved by means of frictional forces selectively overcome or engaged by the acceleration of the vibrating plate.
The control of motion of the objects on the UPM is achieved through a closed loop configuration consisting of a camera, for photographing the horizontal plate and the objects thereon, a set of motors for vibrating the plate and a computer for a) determining the positions of the objects at each sampling interval, b) computing the forces to be applied to each object so that the object follows its predesignated trajectory, c) computing the motion of the plate which will bring about all of the required forces, and d) applying a signal to each actuator so as to move the plate in the required manner. The process is repeated periodically according to a predetermined sampling schedule.
The UPM control method determines, at each sample period, a center of rotation (COR) about which the plate is to be rotated and the magnitude (i.e., duration) of the rotation. By strategically placing the COR at each sampling period, the required forces are generated, in a time-averaged sense, so as to move the objects in their respectively assigned trajectories.
Whereas the UPM illustrates that parallel control of multiple objects on a common medium is possible, its method of control cannot be applied to systems where gravity has much less influence on the objects than do other forces. For example, in fluidic realms, the effects of turbulence and fluid viscosity are as significant as those due to gravity. Fluid flow, in general, is a complex process presenting exceptional control challenges. Some of the unwanted effects of turbulence may be mitigated by controlling the fluid on a small size scale where the momentum of the fluid reaches negligibility. However, the control of fluid flow by acceleration (i.e., by relying on gravitational forces) on such size scales becomes highly impractical.
As shown by consideration of the shortcomings of the prior art, there is an apparent need for parallel control of multiple objects suspended or immersed in a fluid such that each object follows an arbitrary trajectory.
One aspect of the present invention provides a method for transporting objects suspended in a fluid respectively along corresponding trajectories. A microfluidic receptacle is provided to contain the fluid and the objects suspended therein. The microfluidic receptacle includes a plurality of actuators installed thereon for creating a force field within the microfluidic device. The force field respectively imparts a corresponding force on each of the objects. The microfluidic receptacle further includes a sensor for determining at least the location of each object therein. A plurality of force fields defining forces on the objects responsive to a set of actuation signals is determined for the microfluidic receptacle. The method of the present invention then determines, at each sampling interval, a destination point on a trajectory corresponding to each object. The method selects a set of force fields from the plurality of force fields for producing the forces on each object to transport it along its corresponding trajectory. A plurality of actuation signals corresponding to the set of fields is selected and respectively applied to each actuator so as to produce the total force field. Once the signals have been applied to the actuators, the method is repeated until all of the objects have traversed their corresponding trajectories.
Another aspect of the present invention provides a method for transporting a plurality of particles suspended in a fluid respectively along a corresponding one of a plurality of trajectories. A microfluidic receptacle is provided to receive the fluid in which the particles are suspended. The microfluidic receptacle includes a plurality of fluid actuators installed thereon for respectively applying a corresponding force on the fluid. A plurality of fluid flow fields defining the fluid flow responsive to a set of actuation signals is determined for the microfluidic receptacle. The method of the present invention then determines, at each sampling interval, a destination point on a trajectory corresponding to each particle. The method selects a set of fluid flow fields from the plurality of fluid flow fields for producing the fluid flow to transport each particle along its trajectory. A plurality of actuation signals corresponding to the set of fluid flow fields is selected and respectively applied to each fluid actuator so as to produce the flow. Once the signals have been applied to the fluid actuators, the method is repeated until all of the particles have traversed their corresponding trajectories.
In a further aspect of the present invention, a method is provided for sorting a plurality of particles suspended in a fluid in accordance with an attribute possessed by the particles. In this aspect of the present invention, the microfluidic receptacle includes a plurality of fluid actuators, a sensor to detect the position and attribute of each particle within the fluid, and a plurality of sorting bin locations for receiving particles having a corresponding attribute. A plurality of fluid flow fields defining the fluid flow responsive to a set of actuation signals is determined for the microfluidic receptacle. The method of the present invention then determines the attribute of each particle via the sensor at each sampling interval. The method then establishes a plurality of trajectories, one for each particle, which directs the particle to the sorting bin location receiving particles of the associated attribute. The method selects a set of fluid flow fields from the plurality of fluid flow fields that produces a fluid flow such that each particle is moved toward its corresponding sorting bin location. The method then selects the actuation signals to be applied to the actuators so that the fluid flow is produced and, once the signals have been applied, repeats the method until each of the particles has arrived at its corresponding sorting bin location.
In a further aspect of the present invention, a method is provided for conforming a strand from a first conformation to a second conformation, where the strand is suspended in a fluid. The microfluidic receptacle is provided with a plurality of fluid actuators and a sensor to detect a position of a plurality of strand segments forming the strand. A plurality of fluid flow fields defining the fluid flow responsive to a set of actuation signals is determined for the microfluidic receptacle. Once the respective position of each of the strand segments has been determined by sensor, a plurality of trajectories, one for each of the strand segments, is established and is directed to a corresponding segment position of that segment in the second conformation of the strand. The method selects a set of fluid flow fields from the plurality of fluid flow fields for that produces a flow that moves the strand segments toward the second conformation segment position. A plurality of actuation signals corresponding to the selected fluid flow fields is applied to the fluid actuators, and the method is repeated until each of the strand segments has arrived to its strand segment location in the second strand conformation.
In yet another aspect of the present invention, a method is provided for redistributing a first volume of fluid from a first distribution to a second distribution, where the first volume of fluid is immersed within a second volume of fluid and is separated therefrom by at least one fluid interface. The microfluidic receptacle is provided with the plurality of fluid actuators and a sensor to detect a position of a first plurality of segments defining the at least one fluid interface. A plurality of fluid flow fields defining the fluid flow responsive to a set of actuation signals is determined for the microfluidic receptacle. Each of the first plurality of interface segments is located via the sensor and a plurality of trajectories is established, each directed to a corresponding interface segment location defining the second distribution of fluid. A set of fluid flow fields is selected from the plurality of fluid flow fields that produces a fluid flow such that the plurality of segments is directed along its corresponding trajectory. A plurality of actuation signals corresponding to the selected fluid flow fields is applied to the plurality of fluid actuators and the method is repeated until the first plurality of interface segments has arrived at the corresponding one of the plurality of interface segments of the second distribution of the first fluid.
Prior to discussing the exemplary embodiments of the present invention, it is believed as beneficial to briefly define certain terminology as used throughout this Application for Patent. For purposes of the following, a “force field” is intended to mean a set of forces in a region of space imparting a total force vector on an object as a function of the object's position in the region of space. In this sense, force fields include, but are not limited to, electric fields in a region, magnetic fields in a region and fluid flow fields in a region.
The term “fluid” is meant to refer to non-solid media including, but not limited to, gases, liquids and gels.
The term “microfluidic” refers to properties of, and processes on, fluids constrained to regions of a physical size scale wherein inertial effects of the fluid are much less than the viscous effects thereof, i.e., a flow having a vanishingly small Reynolds number. Additionally, gravitational forces on the fluid at the physical size scale are negligible. Typically, the size scale of a microfluidic process or system is less than 1 mm.
An “object” is to be understood as an inanimate (i.e., not self-propelling) particle of arbitrary shape, an inanimate chain or strand or a region of space occupied by a fluid.
Object “conformation” is to be understood as any aspect of an object's internal or external shape or orientation. Additionally, conformation refers to an object's internal or external physical state (e.g. temperature, shear stress, etc.).
An “actuator” refers to a means for generating a force field within a region of space. This includes, but is not limited to, electrodes forming electric and magnetic fields, means for applying forces on fluids by, for example, pressure, electrical, magnetic or electromagnetic fields, thermal processes, surface tension, or electro-osmosis.
A “force field eigenmode” (or alternatively, simply “eigenmode”) is a characteristic force field within a microfluidic channel or receptacle as constrained by the channel geometry, actuation type, and actuator placement.
Referring now to
In the example of
At periodic time intervals, i.e., once every sampling period, the positions of each of the particles within the fluid channel 102 are sensed by a particle position sensor, such as by exemplary methods discussed in paragraphs that follow. From each of these particle locations, a local particle trajectory is computed to set each particle on a path towards its corresponding sorting bin location 104a-104c. By means of the present invention, it is then determined what local force field vector would cause a corresponding force on the particle so that the particle is conducted along its corresponding trajectory. The force may be, for example, an electric or magnetic field or may be a flow vector of the fluid in contact with the object. For purposes of description, the latter of the fields is assumed.
As is shown in the Figure, microfluidic device 100 is adapted to receive a plurality of actuators 106a-106h. The actuators 106a-106h are independently operated to produce a local fluid current flow about its corresponding location. It is an object of the present invention to respectively apply an appropriate signal to each actuator 106a-106h to produce an underlying current flow of the fluid suspending particles 120, 130, 140 so that the force applied to each particle is approximately equal to the force required to transport that particular particle along its corresponding trajectory.
It should be clear from the simplified diagram of
By way of example, the system of
The present invention provides a method by which multiple objects within a microfluidic receptacle may be respectively transported along corresponding trajectories. A microfluidic receptacle of a given geometry and having a given actuation mechanism (e.g., a complex electric field generated by a plurality of electrodes or a complex fluid flow generated by pumps or other means described below) is capable of supporting a set of force fields therein. As will be shown in paragraphs that follow, selecting a subset of supported force fields provides stable, parallel control of the motion of the objects within the microfluidic receptacle.
To implement the present invention, a model of the forces on the objects contained within a microfluidic system provides insight as to hoe the force fields supported by the device may be used in the control method of the present invention. Once the supported force fields have been ascertained, they may be subsequently selected and combined to control the motion of the objects. An exemplary embodiment of the present invention is now presented in a microfluidic system implementing an electro-osmotic actuation mechanism.
The fluid flow in any microfluidic system having a minimum device length well above the mean free path of the molecules composing the fluid can be accurately modeled by the Navier-Stokes equations:
where v* is the fluid velocity field, ρ* is the density, P* is the pressure, η* is the viscosity. Here the asterisk denotes dimensional quantities.
For the microfluidic devices of interest, inertial effects are negligible compared to the effects of viscosity and Equation (2) may then be reduced to:
It is thus apparent from Equation (3) that a change in the fluid flow field may be brought about by either a gradient in the pressure of the fluid or by viscous coupling to a moving object (since inertia is minimal at these scales). It is the choice of actuation and a corresponding choice in microfluidic channel geometry that defines the controllability of the fluid flow. Equation (3) may thus be used to derive a set of control equations once the geometry and actuation mechanism for a particular microfluidic application have been chosen.
An exemplary embodiment of the electro-osmotic microfluidic device is illustrated in
When an aqueous (polar) fluid, such as water, is introduced to the microfluidic chamber 220, an electrical double layer of ions, commonly referred to as the Debye layer, is formed at the plate/liquid interface as is shown in
When a small object 210 is placed in the microfluidic chamber 220 of microfluidic device 200, the trajectory thereof follows the local fluid flow of the suspended fluid. Thus, as will be detailed in paragraphs that follow, the trajectory of object 210 may be controlled by selective application of an electric potential on one or more of the plurality of electrodes 203. The electric field produces a corresponding fluid velocity field serving as the exemplary force field.
The Debye layer thickness, hD, (10 nm in most cases) is very small compared to the chamber dimensions (1 cm×1 cm×0.05 mm in the exemplary device), the boundary conditions at the walls are accurately captured by the velocity slip conditions,
where ε* is the permittivity of the fluid and ξ* is the zeta potential at the wall. The electric potential Φ* satisfies the Laplace equation
∇(ε*∇Φ*)=0 (5)
with the boundary conditions of the applied voltage at the ith control electrode given by
Φ*(δDi)=ui* (6)
where ui* is the electric potential of the ith electrode, and δDi denotes the electrode surface.
If the viscosity and surface properties are uniform and the fluid velocity at the inlets is given by (7):
then the quasi-steady state solution to partial differential equations (1) and (2) is simply given by
showing that the velocity profile is uniform in the vertical z direction (perpendicular to plates 205 and 207). Though the condition (7) is not satisfied at the surface of the electrodes (at the electrode surface the velocity should vanish to zero due to no-slip boundary conditions), the flow field relaxes to satisfy this condition within a few lengths of the channel height from the electrode surfaces and so (8) can be used to predict the fluid velocity in the domain except very near to the electrode surfaces.
The condition on how slowly the electrode voltages can be varied for the fluid flow to maintain a quasi-steady state should be determined to ensure controllability. The time required for the fluid flow in the exemplary device filled with water (ρ*=1000 kg/m3, v*=1 mm/s, h*=0.05 mm, η*=0.001 Ns/m2) to reach the steady state of (8) in response to a step voltage input has been determined to be on the order of 1 μs. Thus, if the electric field is varied such that the time period of the highest frequency input voltage is much greater than 1 μs, the velocity of the fluid at all times may be given by (8).
The motion of an object within the fluid follows the motion of the fluid local to the object. Thus, at any time t, the velocity of the jth object at position pj*=(xj*, yj*) is given by
For convenience, the dimensionalized equations above may be non-dimensionalized as follows. Since the numerical values of x*, y*, |pj*| are of the order of L*; and Φ*, u* are of the order of |umax*| and |vmax*| is of the order of |vmax*|, the relations
can be used to write the equations of the system in dimensionless variables as
∇2Φ=0 (11)
v=∇Φ (12)
{dot over (p)}j=v(pj) (13)
with boundary conditions
Φ(δDi)=ui (14)
The non-dimensionalized system of equations (11), (12), and (14) are linear and contain no time derivatives, so at any time t the velocity field v can be expressed as a superposition of velocity fields ∇Φi(1, 2, . . . , n) as follows
where Φi solves
∇2Φi=0 (16)
with boundary conditions
Φ(δDi)=1, Φ(δDj)=0, j≠i (17)
where n is the number of electrodes and u=[u1 u2 . . . un] represents the vector of electrode voltages. An exemplary velocity field ∇Φ is shown in
Now, since u and u−c (c being some arbitrary constant) generate the same velocity field v, a vector u may always be chosen such that a particular electrode voltage is always 0. In other words any achievable velocity field can be expressed as a linear superposition of any n−1 fields out of ∇Φi (i=1, 2, . . . , n) which constitute a linearly independent set. Hence we can rewrite (15) as
If at time t, the objects are at positions p1=(x1, y1), p2=(x2, y2), . . . , pm=(xm, ym), then the velocity of the jth object is given by
Let v1D, v2D, . . . , vmD be the desired object velocities at time t. It is desired to proportionally combine the velocity fields ∇Φi so that the fluid velocities at p1, p2, . . . , pm are as close to v1D, v2D, . . . , vmD as possible. Such a voltage vector uopt can be obtained by solving the least squares problem
the analytical solution of which is given by
uopt=[AT(p)A(p)]−1A(p)vD. (22)
When 2m is less than n−1 (i.e., when the twice the number of objects is less than one less than the number of actuators), Equation (22) constitutes an underdetermined system of linear equations which has multiple solutions. The least square optimal solution is that with the smallest norm which allows the required velocities to be reached with minimal applied voltage. Ideally, the optimal solution generates the required velocity field by combining the lowest order eigenmodes of the available force field eigenmodes. The higher eigenmodes require large voltage components in the solution uopt while having an insignificant impact on vD. As a result, small changes in vD introduces large changes in uopt. This instability imposes significant obstacles in the control of the fluid and thereby the motion of the objects suspended therein.
Increasing the number of objects imposes other restrictions on the velocity field, i.e., higher eigenmodes of the velocity field must be evoked. Thus, attempting to overcome the control deficiencies by adding actuators to the system fails to do so. To stabilize the system, only a subset of the force field eigenmodes are implemented in the control method, by way of the present invention, as will now be discussed.
Methods for suppressing the higher eigenmodes and to stabilize the least squares solution are widely known. For example, Tikhonov regularization works by solving the modified problem,
Another optimization method is TSVD (truncated singular value decomposition), which as the name suggests, works by expressing the solution space as a superposition of only a finite number the lower eigenmodes and discarding the higher eigenmodes of the system. In certain embodiments of the present invention, TSVD is utilized to select the set of implemented force field eigenmodes. From the solution of Equation (15),
on a set of discrete q×q rectangular grid points ri in the domain, where q is chosen such that the grid may resolve the fluid velocity field sufficiently. Thus,
The ith force field (fluid flow) eigenmode is then given by
Ei=Σj=1nNji∇Φj(x,y)
The TSVD forces the column vectors M and N to satisfy:
∥Mi∥2=1 and ∥Ni∥2=1 (28)
Thus, the application of unit voltage vector Ni produces a unit velocity field vector Mi amplified by σii. The ratio σ11/σii indicates the strength of the first mode of the velocity field in comparison to the ith mode when a voltage vector of identical strength is applied to the corresponding electrodes of the exemplary device in both cases. This ratio is illustrated in
As previously discussed, the control of the motion of objects within the microfluidic device is stabilized by selecting a subset of possible eigenmodes for implementation. If the subset of is chosen as k eigenmodes, k≦m, a new matrix, Ñ can be constructed from the first k columns of N. Then, a new matrix, Ã, may be formed from the original matrix A such that,
Ã(p)=A(p)Ñ[ÑTÑ]−1ÑT, (29)
and, then,
{dot over (p)}=Ã(p)·u. (30)
The actuator signals for the stabilized control method is then given by,
uD=[ÃT(p)Ã(p)]−1Ã(p)vD (31)
where uD is the stabilized actuator input corresponding to the desired force field, vD.
A further design consideration lies in the number of actuators to implement. As previously stated, only the lower order fluid flow eigenmodes contribute significantly to the fluid flow. As such, a smaller number of actuators (electrodes in the exemplary device) need be installed. This is shown in the graph of
Having now described an exemplary microfluidic device, a system for implementing the method of the present invention is described with reference to
Microfluidic device 200 is coupled to sensor 730 by which the positions of objects within the microfluidic device are located and the attributes of the objects may be ascertained. The present invention may be implemented by any object sensor means capable of locating the objects to within a desired accuracy. Such sensor means includes, but is not limited to, microscopic cameras, fiber optics, electromagnetic sensors and thermal sensors. The sensor should, of course, be appropriate to determine the property of the object upon which control thereof is conditioned.
Coupling of the sensor 730 to microfluidic device 200 need not be a physical connection (as indicated by the dashed line). For example, sensor 730 may be a digital camera fitted with optics appropriate to view the objects in microfluidic device 200. In certain embodiments of the present invention, the optics may be part of a microscope, the camera being coupled to the ocular port thereof. In such embodiments, it should be apparent that one or more of the plates 205, 207 of microfluidic device 200 be transparent to provide optical access to the objects in the microfluidic chamber 220. Images captured by the camera may be processed by known image processing techniques to determine an object's position, conformation, or other attribute, an exemplary embodiment of which will be described in paragraphs below.
Actuators 203a-203d are respectively coupled to an amplifier or transducer of amplifier/transducer stage 720. Amplifier/transducer section 720 converts electrical signals from computer 710, into signals appropriate to the actuator type. For example, if actuators 203a-203d are driven by the pressure of a gas, amplifier/transducer section 720 controls the gas pressure for each actuator responsive to an electrical signal supplied by computer 710. In other embodiments, such as the electro-osmotic device described hereinabove, amplifier/transducer section 720 conditions the electrical signal from computer 710 to an electrical signal having an electric potential appropriate to the corresponding actuator.
Computer 710 is configured as the controller of the system of
As is illustrated in
In certain embodiments of the present invention, system controller 715 is implemented as a Time Varying Linear Quadratic Regulator (TVLQR) such as is well-known in the optimal control art. The TVLQR may be adapted to compensate for errors in position of the objects being controlled by system controller 715.
As was developed in paragraphs hereinabove, the desired position of an object pD(t+Δt) is set by an input to the actuators by the relationship
pD(t+Δt)=A[pD(t+Δt)]u(t+Δt) (32)
As is the case with any physical system, there will be errors in the object's position after the actuation cycle has been completed. This positional error, denoted hereinafter as pE(t), results from shortcomings of the model to accurately portray the physical process of the fluid flow within the microfluidic device, deviations from nominal in the actual actuator signals, movement of the object through Brownian motion, as well as other system noise sources.
If the position of the object is represented as a sum of the desired position and the positional error, i.e.,
pT(t)=pD(t)+pE(t), (33)
then the position of the object at some future time t+Δt is given by
{dot over (p)}T(t)=pD(t)+pE(t)+A[pD(t+Δt)+pE(t+Δt)][uD(t+Δt)+uE(t+Δt)] (34)
where,
{dot over (p)}T(t)=pT(t+Δt). (35)
The vector uE(t) represents the actuation signal that would bring about the positional error pE(t).
In a manner consistent with optimal control theory, the positional error of the objects under control of system controller 715 may be modeled by
{dot over (p)}E(t)=F(t)pE(t)+G(t)uE(t) (36)
where F(t) and G(t) contain the coefficients to the differential equation characterizing the error behavior of the open loop system. Using this model for the positional error, {dot over (p)}E(t), it is desired to correct the positional deviation from an initial non-zero state pE(t0) to a final state pE(tf)=0, where tf is the time within which the positional error is to be extinguished. The corrective actuation signal, uE(t) can be found by minimizing the cost function
where Q(t) and R(t) are positive definite matrices and are chosen to insure realizability of uE(t) and pE(t). The actuation signal that solves this optimization problem is known to be,
uE(t+Δt)=−K(t)pE(t) (38)
where,
K(t)=Q−1(t)G(t)S(t) (39)
The matrix S(t) is the solution to the Ricatti differential equation,
{dot over (S)}(t)=−S(t)F(t)−FT(t)S(t)+S(t)G(t)R−1(t)GT(t)S(t)−Q(t) (40)
S(tf)=0, (41)
where
{dot over (S)}(t)=S(t+Δt) (42)
A block diagram of an exemplary set of processing functions of computer 710 is depicted in
As shown at block 830 of
The desired trajectories are provided to actuation processor 840, which determines the force field required to transport the N objects along their respective trajectories. Actuation processor 840 provides at its output a vector of actuation signals that, when applied to the actuation system of the microfluidic device, produces the computed force field.
At prescribed sampling intervals, an image of the microfluidic chamber 220 is acquired via sensor 730 as shown at block 810. In certain embodiments of the present invention, the image data are in the form of pixels, each having a value corresponding to physical properties of the fluid and objects suspended therein at the corresponding pixel location. For example, in the case where the image is constructed of optical data, the value of each pixel may correspond to optical transmissivity of the fluid (and objects suspended therein) at the location in the microfluidic chamber 220 corresponding to the pixel location.
The captured image is transferred to image processor 820 by which the N objects are located and characterized. Any known image processing technique capable of isolating multiple objects of varying conformation, density, distribution, etc. may be utilized by the present invention and an exemplary embodiment will be described with reference to
The actual object locations as determined by the image processor 820 are compared with the desired trajectories computed at block 830 at summing node 860, the output of which is a positional error for each of the N objects within the microfluidic chamber 220. The positional error for each object is used to compute a corresponding corrective trajectory at block 870. The corrective trajectory is one which is directed from the actual location of the corresponding object to its desired position. The corrective trajectories are provided to eigenmode processor 880 which determines the underlying fluid flow to move each object along its corrective trajectory and provides at its output the corresponding actuation signals. The actuation signals of eigenmode processor 840 and those of eigenmode processor 880 are combined at summing node 850 to produce a set of total actuation signals which are subsequently applied to the actuators of the microfluidic device.
Referring now to
Initially, Δt=0 and the actuation signal uD(t) is applied to the actuators as described by block 930. At the next sampling period, i.e., t+Δt (set by delay unit 935), the positions of all the objects in the microfluidic chamber 220 are determined by the sensor 730 to produce a position vector pT(t). The output of summing node 945 is the error position vector pE(t)=pT(t)−pD(t). The error position vector is applied to the dynamic gain block 950 to produce the error actuation signal for the next sample period, uE(t+Δt). The error actuation signal is applied to the desired actuation signal in subsequent operations to produce the total actuation signal uT(t+Δt), which is applied to the actuation block 930.
Referring now to
At block 1040, the Ricatti differential equation is solved for S(t). In certain embodiments of the present invention, such as when the object trajectories have been precalculated prior to the instantiation of the method, the solution for S(t) may be determined off-line and stored for each of a plurality of predetermined regions of microfluidic chamber 220. The regions for which S(t) are solved are determined in accordance with its adjacency to each object trajectory.
Once a solution for the Ricatti differential equation has been computed in block 1040, flow is transferred to block 1045 in which the feedback gain K(t) is computed. Once again, values for K(t) may be computed off-line when the object trajectories were predetermined.
In block 1050, a corrective actuation signal is computed from the positional error of each object and, in block 1055, an open loop actuation signal is determined from the desired object position. It should be clear to one of ordinary skill in the art, that the order in which the corrective voltage and the open loop voltage are computed does not affect the outcome of the method of the present invention.
As shown at block 1060, the total actuation signal is computed as the sum of the open loop actuation signal and the corrective actuation signal and is applied to the actuators of the microfluidic device. Flow is transferred to block 1065 at which it is determined if the next sample period has arrived. If not, the method is suspended until the next sample period. Once the next cycle time has arrived, the method repeats at block 1015 until each object is in its final position as determined at block 1020.
Referring now to
In certain embodiments of the present invention, an image captured by camera 1440 is passed through a threshold filter to convert the image to a binary form. As is shown in
As shown in
An exemplary image processing method is illustrated by way of the flow chart of
Once the objects have been detected, they may be identified and distinguished from one another by attribute. Each object may then be labeled, as shown at block 1880, by attribute or as an object previously identified, such as for purposes of object tracking. The positions of the center of the objects are identified in block 1890 and are transmitted to the object control algorithm. The object position may be determined by known techniques in the image processing art.
By way of the method of the present invention, operations on several different types of objects may be accomplished, as shown in
Although the invention has been described herein in conjunction with specific embodiments thereof, many alternatives, modifications, and variations will be apparent to those skilled in the art. The present invention is intended to embrace all such alternatives, modifications, and variations that fall within the spirit and broad scope of the appended claims.
Number | Date | Country | |
---|---|---|---|
60500777 | Sep 2003 | US |