The invention concerns a novel acoustofluidic device to separate acoustically active particles from fluids comprising a novel device arrangement for improved acoustic pressure and particle velocity; and a method of separating particles from a fluid comprising use of same.
Acoustophoresis is the separation of particles using acoustic wave. It has been known that acoustic wave can exert forces on particles in the acoustic field which can be determined by the particles' volume, density and compressibility. The pressure profile in a standing acoustic wave contains areas of net zero pressure at the nodes and of maximum at the anti-nodes. Depending on the acoustic contract factor of the particles, they will be trapped at the pressure nodes or pressure anti-nodes of the standing acoustic wave.
A wide range of acoustofluidic devices integrating a channel and acoustic transducer have been developed for applications in biochemistry and biomedicine. Acoustofluidic technology involves the use of acoustic radiation force to generate acoustic pressure across a flow channel. Recently, acoustofluidic manipulation of microparticles (i.e. particles with dimensions between about 0.1 and about 1000 μm) such as bacteria, blood cells, circulating tumour cells (CTCs) and extracellular vesicles, has received attention in biochemical, biophysical, and biomedical areas due to its biocompatible, versatile, contactless and label-free advantages.
Devices based on surface acoustic waves (SAWs) producing acoustophoretic motion independent of the acoustic impedance ratio at the boundary which allows the water-filled channel to be made from either acoustically hard (e.g. silicone, pyrex) or soft materials (e.g. polydimethylsiloxane (PDMS)). SAWs versatility also enables droplet actuation with free boundary conditions. SAW based acoustophoretic devices integrate an acoustic source and channel, where the acoustic energy is coupled to fluid in the channel via the soft and/or hard walls defining said channel.
The principle of the separation by SAW devices is driven by the primary acoustic radiation force Frad, and acoustic streaming drag force Fdrag induced by the acoustic waves. Due to the attenuation of the acoustically soft channel material, acoustic waves propagating inside the channel are absorbed into the channel material resulting in acoustic energy loss. Techniques have been explored to manipulate the particle in the channel and controlling the movement along the vertical direction by adjusting the input power from the transducers, however, smaller particles require higher input power. Further, the height of channel further limits the size of samples that can be processed. Increasing the input power of the SAW device may be able to compensate the loss but the induced Joule heat on the interdigital transducers (IDTs) can damage the piezoelectric substrate, such as lithium niobate (LiNbO3) which has high electro-mechanical coupling coefficient but poor thermal conductivity. Thus, the maximum power received by the SAW device is typically constrained by thermal stress produced by the IDTs on the substrate.
Therefore, conventional acoustophoresis devices have had limited efficacy due to several factors including inefficient heat dissipation and weak mechanical sustainability. Improved acoustophoresis devices using innovative acoustic structure are therefore desirable.
In this work, we have developed an alternative acoustofluidic device wherein we have employed a channel having a roof comprised of an active acoustic source. By doing so, significantly reduced acoustic energy loss and increased acoustic pressure inside the channel is achievable leading to vastly improved channel flow velocities compared to conventional devices. Further, the provision of an active acoustic source can be carefully controlled as an adjustable actuator, which can further increase the acoustic energy density in the channel to enhance acoustic manipulation of particles.
According to a first aspect of the invention there is provided an acoustofluidic device comprising:
Reference herein to an IDT refers to a transducer comprising two interlocking comb-shaped arrays of metallic electrodes (in the fashion of a zipper), also known as interdigital electrodes (IDEs). These metallic electrodes are deposited on the surface of a piezoelectric substrate to form a periodic structure and, upon application of radio frequency (RF) voltage, convert electric signals to surface acoustic waves (SAW) by generating periodically distributed mechanical forces via a piezoelectric effect causing the substrate to expand and contract. In the case of the present arrangement, when an RF voltage is applied to the at least one IDT this generates a single travelling SAW within the piezoelectric substrate that propagates toward the channel positioned adjacent thereto. Without wishing to be bound by theory, internal reflection of the SAW from the channel wall reflects the SAW towards the IDT in a counter direction such that the outgoing SAW and reflected SAW travel towards one another wherein they interfere to generate a standing surface acoustic wave (SSAW) field of periodic pressure nodes and pressure antinodes in channel.
More preferably, the device comprises at least a pair of interdigitated transducers (IDTs) deposited on the surface of a piezoelectric substrate to form at least one SSAW transducer wherein the at least one channel is positioned between said at least one pair of IDTs. As is known to those skilled in the art, SSAW transducers comprise a piezoelectric substrate patterned with at least two interdigital transducers IDTs
According to this preferred embodiment of the invention, the at least one channel is positioned between said pair of IDTs at the point wherein the SSAW field is generated. While the SSAWs transmit along the fluid/solid interface of the channel(s) as a transverse wave, upon entering the fluid this becomes a longitudinal wave and causes a pressure field inside the fluid. As a result, particles suspended in the fluid are subjected to lateral acoustic forces due to the acoustic radiation and pressure fluctuations, which allow manipulation of the particles by changing the parameter of the IDTs such as input power, frequency, length, number of electrodes, and spacing between two IDTs.
Reference herein to a piezoelectric substrate refers to any material that exhibits a piezoelectric effect, that is the internal generation of electrical charge resulting from an applied mechanical force (and so also exhibit the reverse piezoelectric effect, that is the internal generation of a mechanical strain resulting from an applied electrical field). Examples include, but are not limited to, polyvinylidene difluoride (PVDF), Gallium Nitride (GaN), Aluminium nitride (AlN), Silicon carbide (SiC), Aluminum Gallium Nitride (AlGaN), Langasite (La3Ga5SiO14), Gallium orthophosphate (GaPO4), a Lithium niobate (LiNbO3), Lithium tantalate (LiTaO3), Barium titanate (BaTiO3), Lead zirconate titanate (Pb[ZrxTi1-x]O3 with 0≤x≤1, or more commonly known as PZT), Potassium niobate (KNbO3), Sodium tungstate (Na2WO3), and Zinc oxide (ZnO). In a preferred embodiment of the invention, said piezoelectric substrate is Lithium niobate (LiNbO3)
In a preferred embodiment, the longitudinal axis of said channel is substantially orthogonal with respect to said IDT(s). As will be appreciated, in this arrangement as fluid flows through said channel, acoustic wave forms generated by the IDT(s) are substantially transverse to the fluid flow path, thus exposing any particles present in the fluid to lateral pressure force permitting generation of fluid flow paths that allow separation of particles according to shape and/or size. Alternatively, the longitudinal axis of said channel is provided at an angle with respect to said IDT(s) (and therefore SAW generated by same), which results in the generation of a tapered SSAW within the channel to facilitate particle separation as the particles flow through the channel. Preferably, said angle is between 0 and 90 degrees or any 1 degree increment therebetween.
In a preferred embodiment, said channel floor is configured to functionally couple with the at least one IDT or SSAW transducer such that the travelling wave or SAWs, respectively, is propagated across same. As will be appreciated, in this preferred embodiment said floor of the channel is provided by the piezoelectric substrate wherein the walls of the channel are bonded to the surface of the substrate to provide a channel.
Alternatively, said channel floor and/or walls can be made from any suitable material that can house a fluid mixture and permits coupling of the acoustic wave energy from the piezoelectric substrate to the fluid in the channel. As is known by those skilled in the art, silicon, glass, or metal materials are commonly used for acoustophoresis because the rigid channel walls provide a near ideal acoustic boundary against the sample fluid, enhancing the required standing wave resonance. Such suitable materials for channel floor and/or wall include, but are not limited to, medical grade plastics, such as polycarbonates or polymethyl methacrylates, polyphenylsulfone (PPS), glass, silicone, ceramic, elastomers, thermoset polyester (TPE), poly-methyl methacrylate (PMMA), polystyrene (PS), polycarbonate (PC), poly-ethylene glycol diacrylate (PEGDA), teflons, polyurethane (PU), paper, hydrogels, pyrex and polydimethyl siloxane (PDMS). Most ideally, said channel walls and/or floor are PDMS and silicone. In yet a further preferred embodiment, the material may be configured to be transparent to permit visualisation of the flow of fluid during operation of the device.
Reference herein to an acoustic wave source refers to any means for generating acoustic wave energy wherein said acoustic wave energy is transmitted into the fluid of the channel. Without wishing to be bound by theory, it is believed that by providing a roof defined by an acoustic wave source one can further generate and enhance a vertical pressure field in the channel, in addition to the longitudinal pressure field induced by the substantially orthogonal IDT(s), producing another standing wave thereby increasing the channel pressure field. Further, the vertical channel pressure field also permits manipulation of particles in the vertical direction of the channel in addition to laterally.
In a preferred embodiment, said acoustic wave source is provided as a further at least one interdigitated transducer (IDT) deposited on the surface of a piezoelectric substrate or standing surface acoustic wave (SSAW) transducer as defined herein. Preferably, said acoustic wave source is provided as a further standing surface acoustic wave (SSAW) transducer as defined herein. As will be appreciated by those skilled in the art, in this arrangement the channel roof and floor is provided by two substantially opposing SSAW transducers which, with the channel walls, define the channel therebetween (SAW-SAW). In this arrangement one can achieve stronger pressure gradients in the channel, thus achieving increased flow rates compared to conventional arrangements known in the art. Further, by tuning the phase difference between the waves generated by the roof and floor SSAW transducers, the acoustic pressure distribution can be controlled permitting greater particle manipulation.
In a preferred embodiment of this arrangement, the IDT(s) deposited on the surface of the piezoelectric substrate and the acoustic wave source are configured such that, in use, a phase difference of between about Δφ=π/2 and Δφ=3π/2 exists between the acoustic wave(s) originating in the piezoelectric substrate (e.g. the floor SSAW transducer) and the acoustic wave(s) originating in the roof of the channel (e.g. the roof SSAW transducer). Most ideally, the phase difference is about Δφ=π, which has been found to produce four symmetrical pressure anti-nodes to form good particle trajectories.
In an alternative embodiment, and more ideally, said acoustic wave source is provided as a bulk acoustic wave (BAW) piezoelectric transducer producing bulk acoustic waves (BAWs). In this arrangement (called BAW-SAW) the channel floor is provided by the at least one IDT or SSAW transducer of the substrate and the roof is defined by a bulk acoustic (BAW) piezoelectric transducer. It has been found that in this particular arrangement one can also achieve stronger pressure gradients in the channel of several orders of magnitude greater than conventional arrangements known in the art. RF signals drive both the BAW and SSAW transducers to produce a combined acoustic energy in the channel. Changing the input voltage of the BAW and SSAW transducers can vary the integrated acoustic field. In a preferred embodiment of this arrangement, the vibration amplitude of the BAW transducer is at least two times the vibration amplitude of the SSAW transducer, more preferably at least five times, and most ideally at least ten times which has been found to produce four symmetrical pressure anti-nodes to form good particle trajectories.
In all embodiments, advantageously it has been found that the pressure gradients achievable can be further increased by cooling the at least one SSAW transducer such that increased input voltages can be applied and thus counter any excess heating of same.
In a preferred embodiment of this arrangement, said BAW piezoelectric transducer is a piezoelectric ceramic such as, but not limited to, PZT, LiNbO3, or the like.
In a preferred embodiment, said channel has a width to height ratio of between about 10:1 and 1:1. More preferably, said channel has a width to height ratio of between about 6:1 and 3:1. In this arrangement, it has been found that where the channel dimensions are proportionally greater in width than height, maximum pressure field can be achieved across the entire cross-section of the channel allowing more careful particle manipulation.
In yet a further preferred embodiment, said channel has a width between about 10-1000 μm including every 1 μm therebetween. More preferably, said channel has a width between about 100-750 μm, and more preferably still between about 300-700 μm, and most preferably between about 400-650 μm.
In a yet further preferred embodiment, said channel has a height between about 1-250 μm including every 1 μm therebetween. More preferably, said channel has a height between about 25-200 μm, and most preferably between about 100-150 μm.
Channels of the following dimensions: (i) 600 μm (W)×125 μm (H); or (ii) 450 μm (W)×120 (H) are particularly suitable for use in the device of the present invention.
In yet a further preferred embodiment of the invention, the channel comprises at least one inlet configured to introduce a fluid into a proximal end portion of the channel. Additionally, or alternatively, the channel comprises at least one outlet which is located at a downstream portion of the channel positioned substantially along the longitudinal axis of the channel. Ideally, the channel comprises at least two outlets. As will be appreciated by those skilled in the art, in this embodiment, a fluid can be introduced through the first inlet and flowed through the at least one channel. Through exposure to the acoustic pressure generated by the SSAW transducer and/or acoustic wave source, particles present in the fluid can be separated into different specific outlets of the channel according to particle shape and/or size. Advantageously, the standing wave from the SSAW transducers and acoustic wave source can control movement of particles both laterally and vertically in channel. Preferably, said inlet(s) and/or outlet(s) are branched to permit separation of particles into different flow streams.
In yet a further preferred embodiment, said inlet(s) and/or outlet(s) comprise tubing to permit flow of a fluid into the inlet(s) and/or out of the outlet(s).
Preferably, said device comprises a pump to control flow rate of fluid through the inlet(s)/channel(s)/outlet(s).
In yet a further preferred embodiment still, said device comprises a plurality of channels in fluid communication with one another. Preferably, the channels are connected in series, so that each channel shares a connection with at least another channel. More preferably, each channel is connected via tubing. More preferably still each channel is functionally coupled with at least one IDT deposited on the surface of a piezoelectric substrate or a SSAW transducer such that each channel can separate different particles with respect to one another according to the standing wave generated for each respective channel. In this arrangement multi-stage particle separation can be achieved.
In yet a further preferred embodiment, the SSAW transducers and/or acoustic wave source can be operated in phase with each other, or operated out of phase with each other. Each SSAW transducer and/or acoustic wave source of the present disclosure may have individual electrical attachments (e.g. electrodes), so that each SSAW transducer and/or acoustic wave source can be individually controlled for frequency and power. Configuration allows for not only the generation of a multi-dimensional acoustic standing wave, but also improved control of the acoustic standing wave. In this way, it is possible to drive individual transducers with arbitrary phasing and/or different or variable frequencies and/or in various out-of-phase modes.
In yet a further preferred embodiment each SSAW transducer and/or acoustic wave source can generate a resonance frequency, or a mean resonance frequency, of between about 100 kHz to 1000 MHz and more preferably between about 1 MHz to 60 MHz.
In yet a further preferred embodiment, it may be required to modulate the frequency or voltage amplitude of the standing wave. This may be done by amplitude modulation and/or by frequency modulation.
As will be appreciated by those skilled in the art, the device disclosed therein herein can be used to separate acoustically active particles from fluids. For example, particles and cells (e.g., target particles) can be removed from a fluid based on the target particles' acoustic properties with respect to the fluid in which they are contained. The fluids can be biological based (e.g., a bodily fluid such as blood) or non-biological based (e.g., waste water). For example acoustic focusing of cells and particles is a technique that can be used in cytometric applications. Acoustic focusing can be implemented in devices for purifying and enriching samples prior to analysis of use of the samples for various applications such as prior to therapeutic injection or diagnosis. In some embodiments, a purified or enriched sample can be integrated into a conventional flow cytometer for further analysis. The acoustic manipulation of particles described herein can be used in clinical applications, requiring the separation of micro- and/or nano-particles. In some embodiments, the invention can be used to manipulate, separate and/or enrich viruses, cells, cell clusters, organisms, tissues, bacteria, exosomes, platelets, parasites, worms, nanotubes, fibres, beads, zebrafish, apoptotic bodies, microvesicles, lipoproteins, liposomes, aerosols, droplets and other nanoparticle and/or microparticle components in biological fluids. In some embodiments, the invention can be used to separate two different sizes of cells. In some embodiments, the invention can be used to separate two different sized cancer cells or to separate cancer cells or disease infected cells (e.g. pathogen infected cells) from healthy cells.
According to a second aspect of the invention, there is provided a method for separating a mixture of particles comprising use of the device as defined herein. As would be appreciated by those skilled in the art, such a method typically comprises suspending a mixture of acoustically active particles in a liquid flow stream and flowing said flow stream through the channel of the device of the first aspect invention, thereby exposing the flow stream in the channel to a standing acoustic wave field to affect acoustic fluid relocation of said acoustically active particles.
Throughout the description and claims of this specification, the words “comprise” and “contain” and variations of the words, for example “comprising” and “comprises”, mean “including but not limited to” and do not exclude other moieties, additives, components, integers or steps. Throughout the description and claims of this specification, the singular encompasses the plural unless the context otherwise requires. In particular, where the indefinite article is used, the specification is to be understood as contemplating plurality as well as singularity, unless the context requires otherwise.
All references, including any patent or patent application, cited in this specification are hereby incorporated by reference. No admission is made that any reference constitutes prior art. Further, no admission is made that any of the prior art constitutes part of the common general knowledge in the art.
Preferred features of each aspect of the invention may be as described in connection with any of the other aspects.
Other features of the present invention will become apparent from the following examples. Generally speaking, the invention extends to any novel one, or any novel combination, of the features disclosed in this specification (including the accompanying claims and drawings). Thus, features, integers, characteristics, compounds or chemical moieties described in conjunction with a particular aspect, embodiment or example of the invention are to be understood to be applicable to any other aspect, embodiment or example described herein, unless incompatible therewith.
Moreover, unless stated otherwise, any feature disclosed herein may be replaced by an alternative feature serving the same or a similar purpose.
The Invention will now be described by way of example only with reference to the Examples below and to the following Figures wherein:
(A) An acoustofluidic configuration comprising PDMS channel sidewalls, a piezoelectric substrate top plate and a piezoelectric substrate bottom plate. A single IDT is deposited on the surface of each of the top and bottom plates, wherein the IDTs of the top and bottom plates are positioned adjacent to and on the same side of the channel sidewalls.
(B) An acoustofluidic configuration comprising PDMS channel sidewalls, a piezoelectric substrate top plate and a piezoelectric substrate bottom plate. A single IDT is deposited on the surface of each of the top and bottom plates, wherein the IDTs of the top and bottom plates are positioned adjacent to and on the opposite side of the channel sidewalls.
(C) An acoustofluidic configuration comprising PDMS channel sidewalls, a piezoelectric substrate top plate and a piezoelectric substrate bottom plate. A single IDT is deposited on the surface of the top plate, and a pair of IDTs forming a SSAW transducer are deposited on the surface of the bottom plate.
(D) An acoustofluidic configuration comprising PDMS channel sidewalls, a top BAW transducer and a piezoelectric bottom plate. A single IDT is deposited on the surface of the bottom plate, adjacent to the channel sidewalls.
(E) An acoustofluidic configuration comprising PZT (or any other piezoelectric substrate) sidewalls, a top BAW transducer and a piezoelectric bottom plate. A single IDT is deposited on the surface of the bottom plate, adjacent to the channel sidewalls. Alternatively, a pair of IDTs forming a SSAW transducer may be deposited on the surface of the bottom plate.
Table 1. Parameters used in numerical analysis at T=25° C.
Referring to the figures and, firstly, to
As is shown in
As will be appreciated, to separate particles in a fluid it is required to provide fluid flow path through said channel [2] through which a fluid containing particles to be separated can flow. In
According to the invention, as shown, a single channel [2] is provided to allow a single staged particle separation stage. However, in alternative arrangements it is envisaged that multiple stages of separation can be achieved. For example, a single channel [2] positioned between multiple staggered pairs of IDTs [3b] may be provided wherein each pair of IDTs [3b] produces a separate SSAW thereby providing differing fields of separation along the flow path. Alternatively, the device [1] comprises a plurality of channels [2] (not shown) in fluid communication with one another. Preferably, the channels are connected in series, so that each channel shares a connection with at least another channel, wherein the outlet of a first channel forms the inlet for the second, and so on. More preferably, each channel is connected via tubing. More preferably still each channel comprises a SSAW transducer such that each channel can separate different particles with respect to one another according to the standing wave generated for each respective channel. In this arrangement multi-stage particle separation can be achieved.
Turning to the channel, referring to
With reference to the channel roof [2d], an acoustic wave source [7] refers to any means for generating acoustic wave energy wherein said acoustic wave energy is transmitted into the fluid of the channel [2]. In this way it is believed that by providing a roof [2d] defined by an acoustic wave source [7] one can further generate a vertical pressure field in the channel, in addition to the longitudinal pressure field induced by the orthogonal IDTs, producing another standing wave thereby increasing the channel pressure field. Further, the vertical channel pressure field also permits manipulation of particles in the vertical direction of the channel in addition to laterally.
In a first embodiment, as shown in
Alternatively, in a second embodiment as shown in
Therefore, as will be appreciated, in all embodiments the SSAW transducers [3] and/or acoustic wave source [7a, 7b] can be operated in phase with each other, or operated out of phase with each other depending on the configuration. Each SSAW transducer [3] and/or acoustic wave source [7a, 7b] of the present disclosure may have individual electrical attachments (e.g. electrodes), so that each SSAW transducer and/or acoustic wave source can be individually controlled for frequency and power. Configuration allows for not only the generation of a multi-dimensional acoustic standing wave, but also improved control of the acoustic standing wave. In this way, it is possible to drive individual transducers with arbitrary phasing and/or different or variable frequencies and/or in various out-of-phase modes.
As will be appreciated, the channel [2] takes the form predominantly of a longitudinal channel whose dimensions (height and width) can vary according to the nature of the fluid to be flowed therethrough, the number of particles to be separated, or respective number of inlet and outlet channels, for example. It has been found that where the channel dimensions are proportionally greater in width than height, maximum pressure field can be achieved across the entire cross-section of the channel allowing more careful particle manipulation. Ideally, said channel has a width between about 400-650 μm and a height between about 100-150 μm, although variations outside these ranges are possible and within the spirit of the invention.
The flow dynamics and particle separation of the devices according to the invention are described in the following examples.
Methods
When acoustic wave is applied to a suspension of particles, the scattering of the wave on the particles will exert an acoustic radiation force that can be utilised to manipulate the particles. To understand how the design of the acoustofluidic devices affects distribution patterns of the particles in a channel, we performed numerical simulations for different design scenarios. In all cases, the channel length is significantly longer than the height and the width, and the acoustic waves are perpendicular to the longitudinal direction. Thus, the fluid flow and particle movement in the channel are investigated as two-dimensional problems.
A. Governing Equations for Fluid Flow
For very dilute suspensions, influences of the particles on the bulk fluid flow can be neglected as long as the particle size is significantly smaller than the dimension of the channel and the acoustic wavelength. Thus, the governing equations for the bulk fluid flow are,
where {circumflex over (ρ)} is the fluid density, the bold letter {circumflex over (v)} is the vector of fluid velocity, {circumflex over (p)} is the fluid pressure, η and ηb are the shear viscosity and bulk viscosity, respectively.
In our devices, there is no fluid flow before application of acoustic waves. As a result, the fluid density and pressure, ρ0 and p0, are uniform and time-independent. When acoustic waves propagate through the fluid, they cause small perturbations in the density, pressure, and velocity fields, which can be expressed as,
{circumflex over (p)}=p
0
+{circumflex over (p)}
1
+{circumflex over (p)}
2+ . . . (4)
{circumflex over (ρ)}=ρ0+{circumflex over (ρ)}1+{circumflex over (p)}2+ . . . (5)
{circumflex over (v)}={circumflex over (v)}
1
+{circumflex over (v)}
2+ . . . (6)
where the subscripts 1 and 2 indicate the first and the second order terms, respectively. Higher order terms are neglected in the simulations. Additionally, we assume that {circumflex over (p)}1 is proportional to {circumflex over (ρ)}1,
{circumflex over (p)}
1
=c
0
2{circumflex over (ρ)}1 (7)
where c0 is a constant and approximately equal to the speed of sound in the fluid. Substituting Eq. (4) through (7) into Eq. (1) and Eq. (2) yields the continuity and momentum equations for the first- and second-order terms:
For periodic perturbations, the time average of the Eq. (10) and (11) become
where X denotes the temporal average of X over an oscillation period. To solve the first order equations, Eq. (8) and Eq. (9) are first combined to obtain the governing equation for {circumflex over (p)}1.
In the study, we assume the first-order fields of the density, pressure, and velocity to be harmonic time dependence, i.e.
{circumflex over (ρ)}1(r,t)=ρ1(r)eiωt (15)
{circumflex over (p)}
1(r,t)=eiωt (16)
{circumflex over (v)}
1(r,t)=v1(r)eiωt (17)
where ω=2πf, which is the angular frequency, and f is the wave frequency. Substituting Eq. (15) through (17) into Eq. (9) and (14) yields,
Eq. (19) can be solved with specific boundary conditions (see descriptions below) to obtain p1, which can be substituted into Eq. (18) to determine the first order velocity, v1(r). In the current study, however, we determined v1(r) using an approximate method. Based on the Helmholtz decomposition theorem, a vector field can be separated into two terms: irrotational and solenoidal. Previous studies have shown that the second term for v1 is negligible in the bulk fluid. It is important only within the boundary layer around a solid surface. Thus, we assumed v1 to be irrotational in the bulk fluid, which can be expressed as the gradient of a velocity potential, v1=∇ϕ1. Substituting this relationship into Eq. (18) yields,
This is the equation used to calculate v1 after solving the governing equation for p1 (i.e., Eq. (19)).
B. Governing Equation for Acoustophoretic Trajectories of Particles
Once the first order acoustic pressure p1 and velocity v1 are obtained, we can determine the time-averaged acoustic radiation force Frad on a spherical particle, which leads to the net movement of the particles besides local oscillation. Frad is the sum of the second-order pressure and the first-order momentum flux integrated over the particle surface,
F
rad=∂Ωn·[{circumflex over (p)}2I+ρ0({circumflex over (v)}1{circumflex over (v)}1)]dA (21)
where ∂Ω is a fixed surface in the bulk fluid around the particle. If the particle radius is much smaller than the wave length, an analytical expression of the force has been derived by Settnes and Bruus,
where the asterisk denotes the complex conjugate of the quantity, a is the particle radius, and κ0 is the isentropic compressibility of the fluid defined as
where V and s are the volume and entropy of the fluid. After neglecting second and higher order terms, κ0=1/(ρ0c02). The scattering coefficients f1 and f2 are calculated by
where ρp and κp is the mass density and compressibility of the particle, respectively. δ is called the viscous penetration depth, which characterises the boundary layer thickness.
Apart from Frad, particles also experience drag force from the viscous fluid due to the relative movement of the particle with respective to the fluid. Since the time-averaged streaming velocity is {circumflex over (v)}2, the time-averaged drag force is,
F
drag=6πηa({circumflex over (v)}2−vp) (28)
where a is the particle radius and vp is the particle velocity vector. Applying the Newton's second law of motion to the particle yields,
where mp is the mass of the particle. In most experimental setup, the particle acceleration time is much shorter than the time scale of experimental observation. As a result, we can neglect the acceleration term in Eq. (29) to obtain an expression for vp:
Once {circumflex over (v)}2 is determined by numerically solving Eq. (12) and Eq. (13) with the boundary conditions described below, Eq. (30) can be used to calculate the particle velocity.
C. Model Configurations and Boundary Conditions
C1. Model Configurations
The four different model configurations considered in this study are shown in
To further increase the acoustic energy pressure in the channel for enhanced manipulation of particles, we developed two new models of acoustofluidic structures as shown in
C2. Boundary Conditions
To solve the first order pressure field, we employ the impendence or lossy-wall boundary condition at the water-PDMS interface, due to partial absorption of the acoustic energy by PDMS:
where ρm and cm are the mass density of PDMS and the speed of sound in PDMS, respectively. n is the normal vector of the solid boundary surface. The same lossy-wall boundary condition also applies to the water-glass interface shown in
where ρg and cg are the mass density of glass and the speed of sound in glass, respectively. To derive the boundary condition at the water-LiNbO3 interface, we considered the LiNbO3 substrate to be actuated by the SSAW, and ignored the wave decay along the propagation path in the substrate because of the short path length. Thus, the displacement and the velocity of the substrate in the z direction at the interface are,
û=u
0
{e
i(ωt-k
y)
+e
i[ωt-k
(w
-y)]} (33)
{circumflex over (v)}=iωu
0
{e
i(ωt-k
y)
+e
i[ωt-k
(w
-y)]} (34)
where u0, t, ks, y and w0 denote displacement amplitude, time, wave number, location on y-axis and channel width, respectively. The continuity of the displacement in the z direction requires the z component of the velocity to be continuous. Using Eq. (20), the boundary condition for p1 at the water-LiNbO3 interface is,
where ez is the unit vector in the z direction. In the design of the SAW-SAW configuration, the same attenuation boundary condition as that shown in Eq. (35) applies to both the top and the bottom boundaries.
For the water-PZT interface, we simulated the design in which PZT vibrated only in the z direction. Thus, the displacement and the velocity of the substrate at this interface are,
û=u
T
e
iωt (36)
{circumflex over (v)}=iωu
T
e
iωt (37)
where uT denotes the maximum displacement amplitude of the PZT surface which is controlled by the applying RF voltage. Again, the continuity of the displacement in the normal direction at the interface requires the normal velocity to be continuous. Using Eq. (20), the boundary condition for p1 at the water-PZT interface is,
In the design with the BAW-SAW configuration, Eq. (38) and Eq. (34) apply to the top and the bottom boundaries respectively.
D. Numerical Simulations
Computational mesh with maximum element size length db at the domain boundary and 10 db in the bulk of the domain is reasonable to capture the physics of the model. We use an illustrative mesh with db=20δ, where δ is the viscous penetration depth defined in Eq. (27). For verifying the correctness of this solution, an investigation of the mesh-convergence is required. We compare a series of meshes with decreasing mesh element size length and define a relative convergence function C(g) for a solution g with respect to a reference solution gref taken to be the solution for the smallest value of db.
where we use a reference solution gref with db=0.4δ, which resulted in 2.2×105 elements.
To simulate the flow and particle distribution patterns, we first solve Eq. (19) to determine the first order pressure field p1 and use it to calculate the velocity field v1 with Eq. (20). The results are substituted into Eq. (12) and Eq. (13) to solve for the time averaged, second order velocity field {circumflex over (v)}2. Finally, the particle velocity and trajcetories are calculated with Eq. (30).
Results
Acoustofluidic Field and Particle Trajectories in the SAW-PDMS and SAW-Glass
The time-averaged second-order velocities v2 in the SAW-PDMS and SAW-Glass are given in
Acoustofluidic Field and Particle Trajectories in the SAW-SAW
In the SAW-Glass channel, the reflected wave from the water-glass interface interacts with the leaky wave in the water produced by the bottom SSAW transducer to produce a pseudo-standing wave (PSW) on the z direction. The PSW can be further improved and controlled by using another SSAW or BAW transducer such as PZT to replace the glass positioned on the top of the channel (
Varying the phase difference Δφ between the top and bottom SSAW transducers can redistribute the pressure gradients and alter the pressure amplitude in the channel, due to the phase shift results in an interchange in the position of the nodes and the anti-nodes. By sweeping the Δφ, the maximum acoustic pressure and pressure gradients are shown in
Acoustofluidic Field and Particle Trajectories in the BAW-SAW
By replacing the top SSAW transducer by the BAW transducer, such as PZT, as shown in
The particle trajectories in the BAW-SAW configuration are simulated in
A comprehensive comparison amongst traditional SAW-PDMS, hybrid SAW-Glass, novel SAW-SAW and BAW-SAW structures have been presented in the study. The model of the SAW-SAW transducers has notably increased the 10-μm particle velocity to 573 μm/s, comparing to the velocity of 10.4 μm/s in the state-of-the-art hybrid SAW-Glass configuration. The active acoustic generation by the top SAW transducer instead of the passive glass reflection employs the same actuation boundary condition as the bottom IDT used in most of acoustofluidic devices. The BAW-SAW transducer whose piezoelectric material (e.g. PZT) can produce much greater vibration to incorporate with the bottom SSAW transducer allowing significantly stronger acoustic resonance generated in the channel. The 10-μm particles migrate at the maximum velocity of 3,310 μm/s when the BAW actuation amplitude is 10 times larger than the SAW actuation. The future work is to manufacture the SAW-SAW and BAW-SAW acoustofluidic chips to verify the model system and numerical analysis.
Number | Date | Country | Kind |
---|---|---|---|
1904601.0 | Apr 2019 | GB | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/EP2020/058531 | 3/26/2020 | WO | 00 |