The ability to separate a particle/fluid mixture into its separate components is desirable in many applications. Physical size exclusion filters can be used for this purpose, where the particles are trapped on the filter and the fluid flows through the filter. Examples of physical filters include those that operate by tangential flow filtration, depth flow filtration, hollow fiber filtration, and centrifugation. However, physical filters can be complicated to work with. As the filter fills up, filtration capacity is reduced. Also, using such filters requires periodic stopping to remove the filter and obtain or clear the particles trapped thereon.
Acoustophoresis is the separation of particles and secondary fluids from a primary or host fluid using high intensity acoustic standing waves, and without the use of membranes or physical size exclusion filters. It has been known that high intensity standing waves of sound can exert forces on particles in a fluid when there is a differential in both density and/or compressibility, otherwise known as the acoustic contrast factor. The pressure profile in a standing wave contains areas of local minimum pressure amplitudes at its nodes and local maxima at its anti-nodes. Depending on the density and compressibility of the particles, they will be trapped at the nodes or anti-nodes of the standing wave. The higher the frequency of the standing wave, the smaller the particles that can be trapped due the pressure of the standing wave.
The use of a single planar acoustic standing wave has been used to accomplish this separation process. However, this single, planar wave tends to trap the particles or secondary fluid in a manner such that they can only be separated from the primary fluid by turning off the planar standing wave. This does not allow for continuous operation. Also, the amount of power that is needed to generate the acoustic planar standing wave tends to heat the primary fluid through waste energy.
Conventional acoustophoresis devices have thus had limited efficacy due to several factors including heat generation, limits on fluid flow, and the inability to capture different types of materials. Improved acoustophoresis devices using three dimensional acoustic standing waves combined with improved fluid dynamics which allow the acoustophoresis particle separation process to be a continuous process driven by gravity are presented herein.
The present disclosure relates to acoustophoretic systems with improved fluid dynamics that can be used to improve separation of particles from a particle/fluid mixture. A new mixture with an increased concentration of particles can be obtained, or the particles themselves can be obtained. In more specific embodiments, the particles are biological cells, such as Chinese hamster ovary (CHO) cells, NS0 hybridoma cells, baby hamster kidney (BHK) cells, and human cells.
Briefly, the systems described herein all include a substantially vertical flow path of the mixture through the acoustic chamber in order to improve separation of particles/secondary fluid from a primary fluid using fluid dynamics. The vertical flow path reduces velocity non-uniformities in the acoustic chamber resulting from gravity forces. In additional systems, a dump diffuser is used to make the incoming flow more uniform, so that the efficiency of the acoustophoretic system is maximized. In other systems, baffles (i.e. laminar plates) are used as collection surfaces for the particles/secondary fluid as the primary fluid flows over them.
Disclosed in various embodiments are acoustophoresis devices comprising: a housing having a sidewall that defines an acoustic chamber; at least one upper outlet in the acoustic chamber; a shallow wall leading to at least one lower concentrate outlet below the acoustic chamber; at least one ultrasonic transducer located on the sidewall of the housing; at least one reflector located on the sidewall of the housing opposite the at least one ultrasonic transducer; and at least one diffuser inlet for the inflow of fluid mixture into the acoustic chamber.
The at least one diffuser inlet may be part of a dump diffuser. The at least one diffuser inlet may be located at a height between 0% and 100% of the height of the acoustic chamber, as measured from the bottom. The at least one diffuser inlet can be in the shape of holes or slots that provide a flow direction normal to the axial direction of an acoustic standing wave generated by the at least one ultrasonic transducer. The shallow wall concentrate outlet can serve to collect particle clusters, and may have an angle of 60° or less relative to a horizontal plane. The device may have a circular or rectangular cross-section in a flow direction.
In some particular embodiments, the at least one diffuser inlet is a plurality of diffuser inlets located about the housing so that the device has a vertical plane or line of symmetry and the inflow of fluid is symmetrical. The plurality of diffuser inlets can be fed by a plurality of dump diffusers.
In other particular embodiments, the at least one diffuser inlet is located at an upper end of the acoustic chamber along with the at least one upper outlet. A wall separates the at least one diffuser inlet from the at least one upper outlet such that fluid flows vertically downwards from the at least one diffuser inlet, then horizontally, then vertically upwards to the at least one upper outlet.
Also disclosed are methods of separating a host fluid from a second fluid or particulate carried in suspension. The methods comprise flowing a mixture of the host fluid and the second fluid or particulate in suspension through an acoustophoresis device in a uniform fashion. The device comprises: a housing having a sidewall that defines an acoustic chamber; at least one upper outlet in the acoustic chamber; a shallow wall leading to at least one lower concentrate outlet below the acoustic chamber; at least one ultrasonic transducer located on the sidewall of the housing, the transducer including a piezoelectric material driven by a voltage signal to create a multi-dimensional acoustic standing wave in the acoustic chamber; at least one reflector located on the sidewall of the housing opposite the at least one ultrasonic transducer; and at least one diffuser inlet in the sidewall of the acoustic chamber above the shallow wall concentrate outlet. The methods further comprise capturing smaller particles of the second fluid or particulate in the multi-dimensional acoustic standing wave to cluster and continuously gravity separate the second fluid or particulate from the host fluid; wherein the gravity separated second fluid or particulate subsequently falls into the at least one lower concentrate outlet. The mixture of the host fluid and the second fluid or particulate can be a slurry.
Disclosed in other embodiments herein are acoustophoresis devices, comprising: a housing defining a primary flow channel between an inlet end and an outlet end of the housing; at least one ultrasonic transducer located on a side of the housing; at least one reflector located on the side of the housing opposite the at least one ultrasonic transducer; and at least one baffle extending into the primary flow channel.
The acoustophoresis device can further include a port on a side of the housing between the at least one transducer and the at least one reflector connecting to a secondary flow channel, the at least one baffle extending into both the primary flow channel and the secondary flow channel. The secondary flow channel may be angled relative to the primary flow channel. In particular embodiments, the angle between the primary flow channel and the secondary flow channel is from about 15° to 90°. At least one baffle is generally angled relative to the primary flow channel, but is not angled relative to the secondary flow channel. A collection well may be located in the secondary flow channel downstream of the at least one baffle.
The acoustophoresis device may have a plurality of baffles, arranged in a variety of configurations. The plurality of baffles can be arranged in parallel to each other. The plurality of baffles can extend into the primary flow channel in at least two different distances. The baffles in the plurality of baffles can be arranged in at least two different angles relative to the primary flow channel. The baffles in the plurality of baffles may all have the same length, or may vary in length.
In some embodiments, the primary flow channel is substantially U-shaped. In others, the inlet end and the outlet end are located on opposite ends of the housing, i.e. the primary flow channel is straight.
Also disclosed are methods of separating a second fluid or a particulate from a host fluid, comprising flowing a mixture of the host fluid and the second fluid or particulate through an acoustophoresis device. The device comprises: an acoustic chamber having at least one inlet and at least one outlet that define a primary flow channel; an ultrasonic transducer located on a wall of the acoustic chamber, the transducer including a piezoelectric material driven by a voltage signal to create a multi-dimensional acoustic standing wave in the acoustic chamber; and at least one reflector located on a wall of the acoustic chamber opposite the at least one ultrasonic transducer; and at least one baffle extending into the primary flow channel. The methods further comprise capturing smaller particles of the second fluid or particulate in the multi-dimensional acoustic standing wave to clump, aggregate, and coalesce and continuously gravity separate the second fluid or particulate from the host fluid; wherein the gravity separated second fluid or particulate subsequently contact the at least one baffle and are directed to a collection outlet.
The particulate may be Chinese hamster ovary (CHO) cells, NS0 hybridoma cells, baby hamster kidney (BHK) cells, insect cells or human cells such as stem cells and T-cells. The mixture may be continuously flowed through the acoustic chamber. The standing wave may have an axial force and a lateral force, the lateral force being at least the same order of magnitude as the axial force.
These and other non-limiting characteristics are more particularly described below.
The following is a brief description of the drawings, which are presented for the purposes of illustrating the exemplary embodiments disclosed herein and not for the purposes of limiting the same.
The present disclosure may be understood more readily by reference to the following detailed description of desired embodiments and the examples included therein. In the following specification and the claims which follow, reference will be made to a number of terms which shall be defined to have the following meanings.
Although specific terms are used in the following description for the sake of clarity, these terms are intended to refer only to the particular structure of the embodiments selected for illustration in the drawings, and are not intended to define or limit the scope of the disclosure. In the drawings and the following description below, it is to be understood that like numeric designations refer to components of like function. Furthermore, it should be understood that the drawings are not to scale.
In the figures, interior surfaces are designated by dashed lines in cross-sectional views, unless otherwise noted.
The singular forms “a,” “an,” and “the” include plural referents unless the context clearly dictates otherwise.
As used in the specification and in the claims, the open-ended transitional phrases “comprise(s),” “include(s),” “having,” “contain(s),” and variants thereof require the presence of the named ingredients/steps and permit the presence of other ingredients/steps. These phrases should also be construed as disclosing the closed-ended phrases “consist of” or “consist essentially of” that permit only the named ingredients/steps and unavoidable impurities, and exclude other ingredients/steps.
All numerical values used herein include values that are the same when reduced to the same number of significant figures and values that differ by less than the experimental error of conventional techniques for measuring that value.
All ranges disclosed herein are inclusive of the recited endpoint and independently combinable (for example, the range of “from 2 grams to 10 grams” is inclusive of the endpoints, 2 grams and 10 grams, and all the intermediate values).
The terms “substantially” and “about” can be used to include any numerical value that can vary without changing the basic function of that value. When used with a range, they also disclose the range defined by the absolute values of the two endpoints, e.g. “about 2 to about 4” also discloses the range “from 2 to 4.” The term “about” may refer to plus or minus 10% of the indicated number.
It should be noted that many of the terms used herein are relative terms. For example, the terms “upper” and “lower” are relative to each other in location, i.e. an upper component is located at a higher elevation than a lower component in a given orientation, but these terms can change if the device is flipped. The terms “inlet” and “outlet” are relative to a fluid flowing through them with respect to a given structure, e.g. a fluid flows through the inlet into the structure and flows through the outlet out of the structure. The terms “upstream” and “downstream” are relative to the direction in which a fluid flows through various components, i.e. the flow fluids through an upstream component prior to flowing through the downstream component. It should be noted that in a loop, a first component can be described as being both upstream of and downstream of a second component.
The terms “horizontal” and “vertical” are used to indicate direction relative to an absolute reference, i.e. ground level. However, these terms do not require structures to be absolutely parallel or absolutely perpendicular to each other, or to the earth. The terms “upwards” and “downwards” are also relative to an absolute reference; an upwards flow is always against the gravity of the earth.
The term “parallel” should be construed in its lay sense of two surfaces that maintain a generally constant distance between them, and not in the strict mathematical sense that such surfaces will never intersect when extended to infinity.
The present application refers to “the same order of magnitude.” Two numbers are of the same order of magnitude if the quotient of the larger number divided by the smaller number is a value less than 10.
The acoustophoretic separation technology of the present disclosure employs ultrasonic acoustic standing waves to trap, i.e., hold stationary, particles or a secondary fluid in a host fluid stream. The particles or secondary fluid collect at the nodes or anti-nodes of the multi-dimensional acoustic standing wave, depending on the particles' or secondary fluid's acoustic contrast factor relative to the host fluid, forming clusters that eventually fall out of the multi-dimensional acoustic standing wave when the clusters have grown to a size large enough to overcome the holding force of the multi-dimensional acoustic standing wave (e.g. by coalescence or agglomeration). This is an important distinction from previous approaches where particle trajectories were merely altered by the effect of the acoustic radiation force. The scattering of the acoustic field off the particles results in a three dimensional acoustic radiation force, which acts as a three-dimensional trapping field. The acoustic radiation force is proportional to the particle volume (e.g. the cube of the radius) when the particle is small relative to the wavelength. It is proportional to frequency and the acoustic contrast factor. It also scales with acoustic energy (e.g. the square of the acoustic pressure amplitude). For harmonic excitation, the sinusoidal spatial variation of the force is what drives the particles to the stable axial positions within the standing waves. When the acoustic radiation force exerted on the particles is stronger than the combined effect of fluid drag force and buoyancy and gravitational force, the particle is trapped within the acoustic standing wave field. This results in concentration, agglomeration and/or coalescence of the trapped particles. The strong lateral forces create rapid clustering of particles. Relatively large solids of one material can thus be separated from smaller particles of a different material, the same material, and/or the host fluid through enhanced gravitational separation.
One specific application for the acoustophoresis device is in the processing of bioreactor materials. It is important to be able to filter all of the cells and cell debris from the expressed materials that are in the fluid stream. The expressed materials are composed of biomolecules such as recombinant proteins or monoclonal antibodies, and are the desired product to be recovered. Through the use of acoustophoresis, the separation of the cells and cell debris is very efficient and leads to very little loss of the expressed materials. This is an improvement over current filtration processes (depth filtration, tangential flow filtration, and the like), which show limited efficiencies at high cell densities, so that the loss of the expressed materials in the filter beds themselves can be up to 5% of the materials produced by the bioreactor. The use of mammalian cell cultures including Chinese hamster ovary (CHO), NS0 hybridoma cells, baby hamster kidney (BHK) cells, and human cells has proven to be a very efficacious way of producing/expressing the recombinant proteins and monoclonal antibodies required of today's pharmaceuticals. The filtration of the mammalian cells and the mammalian cell debris through acoustophoresis aids in greatly increasing the yield of the bioreactor. The acoustophoresis process, through the use of multidimensional acoustic waves, may also be coupled with a standard filtration process upstream or downstream, such as depth filtration using diatomaceous earth, tangential flow filtration (TFF), or other physical filtration processes.
In this regard, the contrast factor is the difference between the compressibility and density of the particles and the fluid itself. These properties are characteristic of the particles and the fluid themselves. Most cell types present a higher density and lower compressibility than the medium in which they are suspended, so that the acoustic contrast factor between the cells and the medium has a positive value. As a result, the axial acoustic radiation force (ARF) drives the cells, with a positive contrast factor, to the pressure nodal planes, whereas cells or other particles with a negative contrast factor are driven to the pressure anti-nodal planes. The radial or lateral component of the ARF is larger than the combined effect of fluid drag force and gravitational force. The radial or lateral component drives the cells/particles to planes where they can cluster into larger groups, which will then gravity separate from the fluid.
As the cells agglomerate at the nodes of the standing wave, there is also a physical scrubbing of the cell culture media that occurs whereby more cells are trapped as they come in contact with the cells that are already held within the standing wave. This generally separates the cells from the cell culture media. The expressed biomolecules remain in the nutrient fluid stream (i.e. cell culture medium).
Desirably, the ultrasonic transducer(s) generate a three-dimensional or multi-dimensional acoustic standing wave in the fluid that exerts a lateral force on the suspended particles to accompany the axial force so as to increase the particle trapping and clumping capabilities of the standing wave. Typical results published in literature state that the lateral force is two orders of magnitude smaller than the axial force. In contrast, the technology disclosed in this application provides for a lateral force to be higher, up to the same order of magnitude as the axial force.
For three-dimensional acoustic fields, Gor'kov's formulation can be used to calculate the acoustic radiation force Fac applicable to any sound field. The primary acoustic radiation force Fac is defined as a function of a field potential U,
FA=−∇(U),
where the field potential U is defined as
and f1 and f2 are the monopole and dipole contributions defined by
where p is the acoustic pressure, μ is the fluid particle velocity, is the ratio of cell density ρp to fluid density ρf, σ is the ratio of cell sound speed cp to fluid sound speed cf, o is the volume of the cell, and < > indicates time averaging over the period of the wave. Gor'kov's formulation applies to particles smaller than the wavelength. For larger particle sizes, Ilinskii provides equations for calculating the 3D acoustic radiation forces for any particle size. See Ilinskii, Acoustic Radiation Force on a Sphere in Tissue, The Journal of the Acoustical Society of America, 132, 3, 1954 (2012), which is incorporated herein by reference.
Perturbation of the piezoelectric crystal in an ultrasonic transducer in a multimode fashion allows for generation of a multidimensional acoustic standing wave. A piezoelectric crystal can be specifically designed to deform in a multimode fashion at designed frequencies, allowing for generation of a multi-dimensional acoustic standing wave. The multi-dimensional acoustic standing wave may be generated by distinct modes of the piezoelectric crystal such as the 3×3 mode that would generate multidimensional acoustic standing waves. A multitude of multidimensional acoustic standing waves may also be generated by allowing the piezoelectric crystal to vibrate through many different mode shapes. Thus, the crystal would excite multiple modes such as a 0×0 mode (i.e. a piston mode) to a 1×1, 2×2, 1×3, 3×1, 3×3, and other higher order modes and then cycle back through the lower modes of the crystal (not necessarily in straight order). This switching or dithering of the crystal between modes allows for various multidimensional wave shapes, along with a single piston mode shape to be generated over a designated time.
It is also possible to drive multiple ultrasonic transducers with arbitrary phasing. In other words, the multiple transducers may work to separate materials in a fluid stream while being out of phase with each other. Alternatively, a single ultrasonic transducer that has been divided into an ordered array may also be operated such that some components of the array will be out of phase with other components of the array.
It may be necessary, at times, due to acoustic streaming, to modulate the frequency or voltage amplitude of the standing wave. This may be done by amplitude modulation and/or by frequency modulation. The duty cycle of the propagation of the standing wave may also be utilized to achieve certain results for trapping of materials. In other words, the acoustic beam may be turned on and shut off at different frequencies to achieve desired results.
The lateral force of the total acoustic radiation force (ARF) generated by the ultrasonic transducers of the present disclosure is significant and is sufficient to overcome the fluid drag force at high linear velocities up to 1 cm/s and beyond. For example, linear velocities through the devices of the present disclosure can be a minimum of 4 cm/min for separation of cells/particles, and can be as high as 1 cm/sec for separation of oil/water phases. Flow rates can be a minimum of 25 mL/min, and can range as high as 40 mL/min to 270 mL/min, or even higher. This is true for batch reactors, fed-batch bioreactors and perfusion bioreactors.
The present disclosure relates to acoustophoretic devices that have improved fluid dynamics for separation of particles/secondary fluid from a primary or host fluid. A U-turn may be present in the flow path of such devices. Baffles (i.e. laminar plates) can also be utilized to increase the concentration of the particles/secondary fluid after being separated from the primary fluid.
The fluid dynamic parameters that are affecting the acoustophoresis separation process may be defined in terms of Reynolds numbers. For instance, the particle Reynolds numbers for particle sizes from one micron to 250 μm with a flow velocity of 4 cm/min to 25 cm/min are from about 0.0005 to 1. The chamber Reynolds numbers for a 1 inch length chamber at a flow rate of 4 cm/min to a 3 inch length chamber at 12 cm/min are about from 10 to 150. The Reynolds numbers based on the inlet tubing diameters of 0.12 inches to 0.31 inches at 30 mL per minute to 540 mL per minute are from about 200 to 1500.
If desired, multiple standing waves from multiple ultrasonic transducers can also be used, which allows for multiple separation stages. For example, in a flow path that runs past two ultrasonic transducers, the first transducer (and its standing wave) will collect a certain amount of the particles, and the second transducer (and its standing wave) will collect additional particles that the first transducer was not able to hold. This construction can be useful where the particle/fluid ratio is high (i.e. large volume of particles), and the separation capacity of the first transducer is reached. This construction can also be useful for particles that have a bimodal or greater size distribution, where each transducer can be optimized to capture particles within a certain size range.
Briefly, in
Referring first to the inlet/outlet module 400, a wall 440 is located between the inlet port 432 and the outlet port 434. Due to the presence of the wall, fluid flows from the inlet port down through the rest of the device and then to the outlet port. The wall essentially divides the primary flow channel into two separate sub-channels, one sub-channel 480 being identified by the inlet port 432 and the wall, and the other sub-channel 482 being identified by the outlet port 434 and the wall. The cross-sectional area of the flow channel for the inlet port can be smaller than, equal to, or greater than the cross-sectional area of the flow channel for the outlet port.
Also visible is a first retainer wall 442 adjacent the inlet port and a second retainer wall 444 adjacent the outlet port. Incoming fluid must flow over the first retainer wall 442 before flowing downwards into the ultrasonic transducer module 200. Similarly, fluid coming back upwards must flow over the second retainer wall 444 before exiting through the outlet port 434. This construction provides a means by which the turbulence of incoming fluid can be reduced, so that the particles trapped in the acoustic standing wave in the ultrasonic transducer module are not disrupted or washed out of the standing wave before aggregating to a sufficient size.
As depicted here, in some embodiments, the wall 440 is spaced apart from the upper end of the housing. This gap 446 forms and acts as a pressure relief passage between the inlet port 432 and the outlet port 434, for example in case the flow path is inadvertently blocked.
An ultrasonic transducer 240 is located on one side of the housing of the ultrasonic transducer module 200, and the reflector 242 is located on the side of the housing opposite the ultrasonic transducer. It should be noted that the ultrasonic transducer is directly adjacent to the primary flow channel, and would be directly exposed to any fluid passing through the primary flow channel.
The particles within the ultrasonic standing wave collect, clump, or agglomerate, and eventually grow to a size where gravity overcomes the acoustic force of the standing wave, and the particle aggregates then fall/sink into the collection well module 300. The collection well module includes a well 330 that tapers downwards in cross-sectional size to a vertex 334. A drain line 340 connects the vertex 334 to a port 342 where the concentrated particles can be drawn out of the well.
The three modules of
This ultrasonic transducer module 900 has a housing 902 with a first end 904, second end 906, and four side walls 910, 912, 914, 916. A primary flow channel 930 is defined between an inlet end and an outlet end, which in this case correspond to the first end and the second end of the housing, as represented by circular openings 932, 934. Put another way, this housing is a flow chamber with an inlet and an outlet. The first end 904 and the second end 906 of the housing can be considered as defining a z-axis. The sides of the housing on which the ultrasonic transducer (not shown) and the reflector (not shown) would be located are represented by square openings 936, 938, and can be considered as defining a y-axis.
In the ultrasonic transducer of
As seen in
In one mode of operation illustrated in
In
In
Finally, in
The acoustophoretic devices discussed above can be made from any suitable material. Such suitable materials include medical grade plastics, such as polycarbonates or polymethyl methacrylates, or other acrylates. It is generally desirable for the material to be somewhat transparent, so that a clear window can be produced and the internal flow channels and flow paths can be seen during operation of the acoustophoresis device/system.
Various coatings may be used on the internal flow channels of the modules. Such coatings include epoxies, for example epichlorohydrin bisphenol crosslinked with an amine or a polyamide; or polyurethane coatings, for example a polyester polyol crosslinked with aliphatic isocyanates. Such coatings are useful for producing a smooth surface and/or reducing surface tension, permitting cells to slide better under the influence of gravity along the flow channel surface and into desired locations (such as collection well modules).
The flow rate of the acoustophoretic device must be controlled so that gravity can act on particle aggregates. In this regard, it is contemplated that the particle/fluid mixture passing in/out of the flow path in the acoustophoretic device through the inlet/outlet modules or the port module can flow at rates of up to about 100 milliliters per minute (ml/min). By way of comparison, the flow rate out of the collection well modules through the ports is much less, from about 3 ml/min up to about 10 ml/min.
Some explanation of the ultrasonic transducers used in the devices of the present disclosure may be helpful as well. In this regard, the transducers use a piezoelectric crystal, usually made of PZT-8 (lead zirconate titanate). Such crystals may have a 1 inch diameter and a nominal 2 MHz resonance frequency. Each ultrasonic transducer module can have only one crystal, or can have multiple crystals that each act as a separate ultrasonic transducer and are either controlled by one or multiple amplifiers.
Screws 88 attach an aluminum top plate 82a of the housing to the body 82b of the housing via threads. The top plate includes a connector 84 for powering the transducer. The top surface of the PZT crystal 86 is connected to a positive electrode 90 and a negative electrode 92, which are separated by an insulating material 94. The electrodes can be made from any conductive material, such as silver or nickel. Electrical power is provided to the PZT crystal 86 through the electrodes on the crystal. Note that the crystal 86 has no backing layer or epoxy layer. Put another way, there is an air gap 87 in the transducer between aluminum top plate 82a and the crystal 86 (i.e. the air gap is completely empty). A minimal backing 58 and/or wear plate 50 may be provided in some embodiments, as seen in
The transducer design can affect performance of the system. A typical transducer is a layered structure with the ceramic crystal bonded to a backing layer and a wear plate. Because the transducer is loaded with the high mechanical impedance presented by the standing wave, the traditional design guidelines for wear plates, e.g., half wavelength thickness for standing wave applications or quarter wavelength thickness for radiation applications, and manufacturing methods may not be appropriate. Rather, in one embodiment of the present disclosure the transducers, there is no wear plate or backing, allowing the crystal to vibrate in one of its eigenmodes (i.e. near eigenfrequency) with a high Q-factor. The vibrating ceramic crystal/disk is directly exposed to the fluid flowing through the flow chamber.
Removing the backing (e.g. making the crystal air backed) also permits the ceramic crystal to vibrate at higher order modes of vibration with little damping (e.g. higher order modal displacement). In a transducer having a crystal with a backing, the crystal vibrates with a more uniform displacement, like a piston. Removing the backing allows the crystal to vibrate in a non-uniform displacement mode. The higher order the mode shape of the crystal, the more nodal lines the crystal has. The higher order modal displacement of the crystal creates more trapping lines, although the correlation of trapping line to node is not necessarily one to one, and driving the crystal at a higher frequency will not necessarily produce more trapping lines.
In some embodiments, the crystal may have a backing that minimally affects the Q-factor of the crystal (e.g. less than 5%). The backing may be made of a substantially acoustically transparent material such as balsa wood, foam, or cork which allows the crystal to vibrate in a higher order mode shape and maintains a high Q-factor while still providing some mechanical support for the crystal. The backing layer may be a solid, or may be a lattice having holes through the layer, such that the lattice follows the nodes of the vibrating crystal in a particular higher order vibration mode, providing support at node locations while allowing the rest of the crystal to vibrate freely. The goal of the lattice work or acoustically transparent material is to provide support without lowering the Q-factor of the crystal or interfering with the excitation of a particular mode shape.
Placing the crystal in direct contact with the fluid also contributes to the high Q-factor by avoiding the dampening and energy absorption effects of the epoxy layer and the wear plate. Other embodiments may have wear plates or a wear surface to prevent the PZT, which contains lead, contacting the host fluid. This may be desirable in, for example, biological applications such as separating blood. Such applications might use a wear layer such as chrome, electrolytic nickel, or electroless nickel. Chemical vapor deposition could also be used to apply a layer of poly(p-xylylene) (e.g. Parylene) or other polymers or polymer films. Organic and biocompatible coatings such as silicone or polyurethane are also usable as a wear surface.
Initially, when a suspension is flowing through the system with primarily small micron sized particles, it is necessary for the acoustic radiation force to balance the combined effect of fluid drag force and buoyancy force for a particle to be trapped in the standing wave. In
Particle size growth continues until the buoyancy force becomes dominant, which is indicated by a second critical particle size, Rc2, at which size the particles will rise or sink, depending on their relative density with respect to the host fluid. At this size, acoustic forces are secondary, gravity/buoyancy forces become dominant, and the particles naturally drop out of the host fluid. As the particles rise or sink, they no longer reflect the acoustic radiation force, so that the acoustic radiation force then increases. Not all particles will drop out, and those remaining particles will continue to grow in size as well. This phenomenon explains the quick drops and rises in the acoustic radiation force beyond size Rc2. Thus,
The size, shape, and thickness of the transducer determine the transducer displacement at different frequencies of excitation, which in turn affects particle separation efficiency. Typically, the transducer is operated at frequencies near the thickness resonance frequency (half wavelength). Gradients in transducer displacement typically result in more places for particles to be trapped. Higher order modal displacements generate three-dimensional acoustic standing waves with strong gradients in the acoustic field in all directions, thereby creating equally strong acoustic radiation forces in all directions, leading to multiple trapping lines, where the number of trapping lines correlate with the particular mode shape of the transducer.
To investigate the effect of the transducer displacement profile on acoustic trapping force and particle separation efficiencies, an experiment was repeated ten times, with all conditions identical except for the excitation frequency. Ten consecutive acoustic resonance frequencies, indicated by circled numbers 1-9 and letter A on
As the emulsion passed by the transducer, the trapping lines of oil droplets were observed and characterized. The characterization involved the observation and pattern of the number of trapping lines across the fluid channel, as shown in
The effect of excitation frequency clearly determines the number of trapping lines, which vary from a single trapping line at the excitation frequency of acoustic resonance 5 and 9, to nine trapping lines for acoustic resonance frequency 4. At other excitation frequencies four or five trapping lines are observed. Different displacement profiles of the transducer can produce different (more) trapping lines in the standing waves, with more gradients in displacement profile generally creating higher trapping forces and more trapping lines.
In the present systems, the system is operated at a voltage such that the particles and particle clusters are trapped in the ultrasonic standing wave. The particles and clusters are collected in well-defined trapping lines. Each trapping line is aligned with the main direction of the acoustic standing wave. Particles and clusters in the trapping lines are separated by half a wavelength. Within each pressure nodal plane of the standing wave, the particles are trapped at very specific points, typically the minima of the acoustic radiation potential. The axial component of the acoustic radiation force drives the particles, with a positive contrast factor, to the pressure nodal planes, whereas particles with a negative contrast factor are driven to the pressure anti-nodal planes. The radial or lateral component of the acoustic radiation force is the force that traps the particles in the standing wave, clumps or clusters them into tightly packed clusters, which then gravity separate when the clusters reach a critical size. In systems using typical transducers, the radial or lateral component of the acoustic radiation force is typically several orders of magnitude smaller than the axial component of the acoustic radiation force. It therefore has two limitations. It has very weak trapping capabilities of particles and moreover, it cannot generate tightly enough packed clusters that will separate out due to gravity. The lateral force in the present devices can be significant, on the same order of magnitude as the axial force component, and is sufficient to overcome the fluid drag force at linear velocities of up to 1 cm/s. The strong clustering capability leads to rapid formation of clusters which continuously separate out from the host fluid through gravity/buoyancy separation.
The three-dimensional acoustic standing waves are the result of superposition of the vibration modes of the crystal. Three dimensional force gradients are generated in every nodal plane of the standing wave. Multiple particle clusters are formed along a line in the axial direction of the standing wave, as illustrated in
The overall system thus operates as shown in
Referring now to
This device 500 includes a symmetrical, dual dump diffuser, plenum inlet configuration. Here, two dump diffusers 530 are placed on opposite sides of the device. Each dump diffuser has a plenum/chamber with an upper end 532 and a lower end 534. The plenum volume provides flow diffusion and dramatically reduces incoming flow non-uniformities. An inlet flow port 542 is located at the upper end 532, and at least one flow outlet 544 is located at the lower end of the plenum. These plenum flow outlets 544 also pass through the acoustic chamber sidewall 508, and can be considered flow inlets into the acoustic chamber. These acoustic chamber flow inlets can be in the form of slots or a line of holes, and they are placed above the bottom of the acoustic chamber. In embodiments, the acoustic chamber inlets are located above the chamber bottom 501 at a height 515 that is between 0% and 100% of the height 517 of the acoustic chamber, and more particularly between 5% and 25% of the height of the acoustic chamber. The diffuser outlets 544 provide a flow direction normal to the axial direction of the acoustic standing waves generated by the ultrasonic transducer. The acoustic chamber inlets are also arranged so that they are in opposing locations, so that the horizontal velocity will decrease to zero in the center of the acoustic chamber.
Each dump diffuser includes an entrance port 542 into which the mixture of host fluid/second fluid or particulate flows (the X refers to the flow direction into the paper). This eliminates downward flow in the acoustic chamber. The mixture fills up the plenum in the dump diffuser and then flows horizontally into the acoustic chamber inlets 544 and enters the acoustic chamber, where the mixture flows vertically upwards and out of the upper outlet 512. The dump diffuser reduces/eliminates flow pulsations and flow non-uniformities that result from pumps, hosing and horizontal inlet flow where gravity effects dominate. The acoustic chamber inlets 544 then bring the heavier mixture into the acoustic chamber above the bottom of the chamber (line 501) and below the ultrasonic transducer and the nodal clusters that form in the ultrasonic standing waves. This minimizes any disturbances of the clusters set up by inflowing material.
The vertical plane or line of symmetry 505 is aligned with gravity forces. Also shown are flow streamlines which are desirably symmetrical, since this minimizes non-uniformities, eddy disturbances, circulation, and disturbance of clusters falling through outlet 514 to be collected. Symmetry also maximizes gravity forces in the inlet flow distribution and particle collection process. Because it is heavier than the permeate exiting at the top of the device, the (relatively) heavy incoming mixture comes in near the bottom of the acoustic chamber, spreads out across the bottom of the chamber due to gravity forces, and provides near uniform velocity profiles from bottom to top. The horizontal velocity of the mixture will decrease to zero as it approaches the center of the acoustic chamber due to the dual opposing inlet flows. This assures minimum interference between the chamber flow and dropping particle clusters. A uniform velocity provides the best separation and collection results because the lateral acoustic forces have to overcome particle drag for the clusters to grow and continuously drop out of the acoustic chamber. This also eliminates the need for an inlet flow distributor.
As the particle clusters drop out, the axial acoustic forces associated with the standing wave will keep the clusters intact. This assures rapid dropping of the clusters with high terminal velocities, on the order of 1 cm/sec. This is extremely fast compared to the chamber flow velocities, which are on the order of 0.1 cm/sec to 0.3 cm/sec. The shallow wall angle means the cylindrical particle clusters have to drop only a very short distance before they exit the acoustic chamber, so that little dispersion of the clusters occurs. Ideally, the system operates with 3 to 12 crystal vibration nodes per square inch of transducer. The symmetry, minimum flow disturbance in the central collection region, and shallow collector walls provide good collection without the need for baffles/laminar plates.
The flow distributors of
In embodiments, the pulsed voltage signal driving the transducer can have a sinusoidal, square, sawtooth, or triangle waveform; and have a frequency of 500 kHz to 10 MHz. The pulsed voltage signal can be driven with pulse width modulation, which produces any desired waveform. The pulsed voltage signal can also have amplitude or frequency modulation start/stop capability to eliminate streaming.
The transducer(s) is/are used to create a pressure field that generates forces of the same order of magnitude both orthogonal to the standing wave direction and in the standing wave direction. When the forces are roughly the same order of magnitude, particles of size 0.1 microns to 300 microns will be moved more effectively towards regions of agglomeration (“trapping lines”). Because of the equally large gradients in the orthogonal acoustophoretic force component, there are “hot spots” or particle collection regions that are not located in the regular locations in the standing wave direction between the transducer 240 and the reflector 242. Hot spots are located at the minima of acoustic radiation potential. Such hot spots represent particle collection locations.
The comparisons between numerical and experimental results are excellent for the multidimensional acoustic standing waves of the present disclosure, and demonstrate the accuracy of the models.
As discussed above, the lateral force can be increased by driving the transducer in higher order mode shapes, as opposed to a form of vibration where the crystal effectively moves as a piston having a uniform displacement. The types of waves thus generated can be characterized as composite waves, with displacement profiles that are similar to leaky symmetric (also referred to as compressional or extensional) Lamb waves. The waves are leaky because they radiate into the water layer, which result in the generation of the acoustic standing waves in the water layer. Symmetric Lamb waves have displacement profiles that are symmetric with respect to the neutral axis of the piezoelectric element, which causes multiple standing waves to be generated in a 3-D space. These higher order modes of vibration can include modes (1,1), (1,2), (2,1), (2,2), (2, 3), or (m, n), where m and n are 1 or greater. The acoustic pressure is proportional to the driving voltage of the transducer. The electrical power is proportional to the square of the voltage.
The transducer is typically a thin piezoelectric plate, which is operated in the (3,3) mode, with electric field in the z-axis and primary displacement in the z-axis, as shown in
In biological applications, it is contemplated that all of the parts of the system (e.g. the acoustic chamber, tubing leading to and from the acoustic chamber, the transducer, dump diffusers, baffles, etc.) can be separated from each other and be disposable. The frequency of the transducers may also be varied to obtain optimal effectiveness for a given power.
The following examples are provided to illustrate the devices and processes of the present disclosure. The examples are merely illustrative and are not intended to limit the disclosure.
Three different devices were made and tested for their cell separation efficiency. The baseline device was similar to
The devices were tested using a 1.5 volume % yeast solution. Each device used two transducers (arranged vertically one on top of the other) with 50 volts to the top transducer and 60 volts to the bottom transducer, operating at a transducer frequency 2.23 MHz. The flowthrough rate was about 270 mL/min, with about 8.1 mL/min concentrate drawoff. The cell separation efficiency was measured as the percentage reduction in turbidity between the feed and the permeate. Table 1 below provides the results:
Next, the 2DD45 device was tested at different flowrates, and the clarification was again measured. Table 2 shows the results:
The results in Table 2 showed that good cell separation could be obtained at up to three times the flowrate of the baseline device when using the dual dump diffuser design.
The present disclosure has been described with reference to exemplary embodiments. Obviously, modifications and alterations will occur to others upon reading and understanding the preceding detailed description. It is intended that the present disclosure be construed as including all such modifications and alterations insofar as they come within the scope of the appended claims or the equivalents thereof.
This application claims priority to U.S. Provisional Patent Application Ser. No. 62/020,088, filed on Jul. 2, 2014, and to U.S. Provisional Patent Application Ser. No. 62/154,672, filed on Apr. 29, 2015, the disclosures of which are hereby fully incorporated by reference in its entirety.
Number | Name | Date | Kind |
---|---|---|---|
2667944 | Crites | Feb 1954 | A |
3555311 | Weber | Jan 1971 | A |
4055491 | Porath-Furedi | Oct 1977 | A |
4158629 | Sawyer | Jun 1979 | A |
4165273 | Azarov et al. | Aug 1979 | A |
4173725 | Asai et al. | Nov 1979 | A |
4204096 | Barcus et al. | May 1980 | A |
4398325 | Piaget et al. | Aug 1983 | A |
4666595 | Graham | May 1987 | A |
4699588 | Zinn et al. | Oct 1987 | A |
4743361 | Schram | May 1988 | A |
4759775 | Peterson et al. | Jul 1988 | A |
4983189 | Peterson et al. | Jan 1991 | A |
5225089 | Benes et al. | Jul 1993 | A |
5371729 | Manna | Dec 1994 | A |
5395592 | Bolleman et al. | Mar 1995 | A |
5431817 | Braatz et al. | Jul 1995 | A |
5443985 | Lu et al. | Aug 1995 | A |
5452267 | Spevak | Sep 1995 | A |
5484537 | Whitworth | Jan 1996 | A |
5527460 | Trampler et al. | Jun 1996 | A |
5560362 | Sliwa, Jr. et al. | Oct 1996 | A |
5594165 | Madanshetty | Jan 1997 | A |
5604301 | Mountford et al. | Feb 1997 | A |
5626767 | Trampler et al. | May 1997 | A |
5688405 | Dickinson et al. | Nov 1997 | A |
5711888 | Trampler et al. | Jan 1998 | A |
5831166 | Kozuka et al. | Nov 1998 | A |
5902489 | Yasuda et al. | May 1999 | A |
5912182 | Coakley et al. | Jun 1999 | A |
5951456 | Scott | Sep 1999 | A |
6090295 | Raghavarao et al. | Jun 2000 | A |
6166231 | Hoeksema | Dec 2000 | A |
6216538 | Yasuda et al. | Apr 2001 | B1 |
6205848 | Faber et al. | Jun 2001 | B1 |
6332541 | Coakley et al. | Dec 2001 | B1 |
6391653 | Letcher et al. | May 2002 | B1 |
6487095 | Malik et al. | Nov 2002 | B1 |
6592821 | Wada et al. | Jul 2003 | B1 |
6649069 | DeAngelis | Nov 2003 | B2 |
6763722 | Fjield et al. | Jul 2004 | B2 |
6881314 | Wang et al. | Apr 2005 | B1 |
6929750 | Laurell et al. | Aug 2005 | B2 |
6936151 | Lock et al. | Aug 2005 | B1 |
7010979 | Scott | Mar 2006 | B2 |
7061163 | Nagahara et al. | Jun 2006 | B2 |
7081192 | Wang et al. | Jul 2006 | B1 |
7093482 | Berndt | Aug 2006 | B2 |
7108137 | Lal et al. | Sep 2006 | B2 |
7150779 | Meegan, Jr. | Dec 2006 | B2 |
7186502 | Vesey | Mar 2007 | B2 |
7191787 | Redeker et al. | Mar 2007 | B1 |
7331233 | Scott | Feb 2008 | B2 |
7340957 | Kaduchak et al. | Mar 2008 | B2 |
7373805 | Hawkes et al. | May 2008 | B2 |
7541166 | Belgrader et al. | Jun 2009 | B2 |
7601267 | Haake et al. | Oct 2009 | B2 |
7673516 | Janssen et al. | Mar 2010 | B2 |
7837040 | Ward et al. | Nov 2010 | B2 |
7846382 | Strand et al. | Dec 2010 | B2 |
7968049 | Takahashi et al. | Jun 2011 | B2 |
8080202 | Takahashi et al. | Dec 2011 | B2 |
8256076 | Feller | Sep 2012 | B1 |
8266950 | Kaduchak et al. | Sep 2012 | B2 |
8273253 | Curran | Sep 2012 | B2 |
8273302 | Takahashi et al. | Sep 2012 | B2 |
8309408 | Ward et al. | Nov 2012 | B2 |
8319398 | Vivek et al. | Nov 2012 | B2 |
8334133 | Fedorov et al. | Dec 2012 | B2 |
8387803 | Thorslund et al. | Mar 2013 | B2 |
8592204 | Lipkens et al. | Nov 2013 | B2 |
8679338 | Rietman et al. | Mar 2014 | B2 |
20020134734 | Campbell et al. | Sep 2002 | A1 |
20030195496 | Maguire | Oct 2003 | A1 |
20030209500 | Kock et al. | Nov 2003 | A1 |
20030230535 | Affeld et al. | Dec 2003 | A1 |
20040016699 | Bayevsky | Jan 2004 | A1 |
20050196725 | Fu | Sep 2005 | A1 |
20060037915 | Strand et al. | Feb 2006 | A1 |
20070272618 | Gou et al. | Nov 2007 | A1 |
20070284299 | Xu et al. | Dec 2007 | A1 |
20080011693 | Li | Jan 2008 | A1 |
20080217259 | Siversson | Sep 2008 | A1 |
20090029870 | Ward et al. | Jan 2009 | A1 |
20090053686 | Ward et al. | Feb 2009 | A1 |
20090098027 | Tabata et al. | Apr 2009 | A1 |
20090178716 | Kaduchak et al. | Jul 2009 | A1 |
20090194420 | Mariella, Jr. et al. | Aug 2009 | A1 |
20090045107 | Ward et al. | Dec 2009 | A1 |
20090295505 | Mohammadi et al. | Dec 2009 | A1 |
20100000945 | Gavalas | Jan 2010 | A1 |
20100078384 | Yang | Apr 2010 | A1 |
20100124142 | Laugharn et al. | May 2010 | A1 |
20100192693 | Mudge et al. | Aug 2010 | A1 |
20100193407 | Steinberg et al. | Aug 2010 | A1 |
20100206818 | Leong et al. | Aug 2010 | A1 |
20100255573 | Bond et al. | Oct 2010 | A1 |
20100317088 | Radaelli et al. | Dec 2010 | A1 |
20100323342 | Gonzalez Gomez et al. | Dec 2010 | A1 |
20100330633 | Walther et al. | Dec 2010 | A1 |
20110024335 | Ward et al. | Feb 2011 | A1 |
20110092726 | Clarke | Apr 2011 | A1 |
20110123392 | Dionne et al. | May 2011 | A1 |
20110154890 | Holm et al. | Jun 2011 | A1 |
20110166551 | Schafer | Jul 2011 | A1 |
20110262990 | Wang et al. | Oct 2011 | A1 |
20110281319 | Swayze et al. | Nov 2011 | A1 |
20110309020 | Rietman et al. | Dec 2011 | A1 |
20120088295 | Yasuda et al. | Apr 2012 | A1 |
20120325727 | Dionne et al. | Dec 2012 | A1 |
20120328477 | Dionne et al. | Dec 2012 | A1 |
20120329122 | Lipkens et al. | Dec 2012 | A1 |
20130175226 | Coussios et al. | Jul 2013 | A1 |
20130277316 | Dutra et al. | Oct 2013 | A1 |
20130284271 | Lipkens et al. | Oct 2013 | A1 |
20140011240 | Lipkens et al. | Jan 2014 | A1 |
Number | Date | Country |
---|---|---|
30 27 433 | Feb 1982 | DE |
0 292 470 | Nov 1988 | EP |
1 254 669 | Nov 2002 | EP |
2 420 510 | May 2006 | GB |
WO 198707178 | Dec 1987 | WO |
WO 02072234 | Sep 2002 | WO |
WO 2009111276 | Sep 2009 | WO |
WO 2009144709 | Dec 2009 | WO |
WO 2010024753 | Apr 2010 | WO |
WO 2011023949 | Mar 2011 | WO |
WO 2011027146 | Mar 2011 | WO |
WO 2011161463 | Dec 2011 | WO |
Entry |
---|
International Search Report for PCT Application No. PCT/US2015/039125. |
Alvarez et al.; Shock Waves, vol. 17, No. 6, pp. 441-447, 2008. |
Benes et al.; Ultrasonic Separation of Suspended Particles, 2001 IEEE Ultrasonics Symposium; Oct. 7-10, 2001; pp. 649-659; Atlanta, Georgia. |
Castro; Tunable gap and quantum quench dynamics in bilayer graphene; Jul. 13, 2010; Mathematica Summer School. |
Cravotto et al.; Ultrasonics Sonochemistry, vol. 15, No. 5, pp. 898-902, 2008. |
Garcia-Lopez, et al; Enhanced Acoustic Separation of Oil-Water Emulsion in Resonant Cavities. The Open Acoustics Journal. 2008, vol. 1, pp. 66-71. |
Hill et al.; Ultrasonic Particle Manipulation; Microfluidic Technologies for Miniaturized Analysis Systems, Jan. 2007, pp. 359-378. |
Kuznetsova et al.; Microparticle concentration in short path length ultrasonic resonators: Roles of radiation pressure and acoustic streaming; Journal of the Acoustical Society of America, American Institute of Physics for the Acoustical Society of America, vol. 116, No. 4, Oct. 1, 2004, pp. 1956-1966, DOI: 1.1121/1.1785831. |
Latt et al.; Ultrasound-membrane hybrid processes for enhancement of filtration properties; Ultrasonics sonochemistry 13.4 (2006): 321-328. |
Lipkens et al.; Frequency sweeping and fluid flow effects on particle trajectories in ultrasonic standing waves; Acoustics 08, Paris, Jun. 29-Jul. 4, 2008. |
Lipkens et al.; Prediction and measurement of particle velocities in ultrasonic standing waves; J. Acoust. Soc. Am., 124 No. 4, pp. 2492 (A) 2008. |
Lipkens et al.; Separation of micron-sized particles in macro-scale cavities by ultrasonic standing waves; Presented at the International Congress on Ultrasonics, Santiago; Jan. 11-17, 2009. |
Lipkens et al.; The effect of frequency sweeping and fluid flow on particle trajectories in ultrasonic standing waves; IEEE Sensors Journal, vol. 8, No. 6, pp. 667-677, 2008. |
Lipkens et al., Macro-scale acoustophoretic separation of lipid particles from red blood cells, The Journal of the Acoustical Society of America, vol. 133, Jun. 2, 2013, p. 045017, XP055162509, New York, NY. |
Meribout et a.; An Industrial-Prototype Acoustic Array for Real-Time Emulsion Layer Detection in Oil Storage Tanks; IEEE Sensors Journal, vol. 9, No. 12, Dec. 2009. |
Nilsson et al.; Review of cell and particle trapping in microfluidic systems; Department of Measurement Technology and Industrial Electrical Engineering, Div. of Nanobiotechnology, Lund University, P.O. Box 118. Lund, Sweden, Analytica Chimica Acta 649, Jul. 14, 2009, pp. 141-157. |
Pangu et al.; Droplet transport and coalescence kinetics in emulsions subjected to acoustic fields; Ultrasonics 46, pp. 289-302 (2007). |
Ponomarenko et al.; Density of states and zero Landau level probed through capacitance of graphene; Nature Nanotechnology Letters, Jul. 5, 2009; DOI: 10.1038/NNANO.2009.177. |
Seymour et al, J. Chem. Edu., 1990, 67(9), p. 763, published Sep. 1990. |
Annex to Form PCT/ISA/206—Communication Relating to the Results of the Partial International Search Report, dated Jul. 18, 2013. |
European Search Report of European Application No. 11769474.5 dated Oct. 10, 2012. |
International Search Report and Written Opinion dated Dec. 20, 2011, for corresponding PCT application No. PCT/US2011/032181. |
International Search Report and Written Opinion dated Feb. 27, 2012, for PCT application No. PCT/US2011/040787. |
International Search Report and Written Opinion of International Application No. PCT/US2013/037404 dated Jun. 21, 2013. |
International Search Report and Written Opinion of International Application No. PCT/US2013/050729 dated Sep. 25, 2013. |
International Search Report for corresponding PCT Application Serial No. PCT/US2014/015382 dated May 6, 2014. |
Phys. Org. “Engineers develop revolutionary nanotech water desalination membrane.” Nov. 6, 2006. http://phys.org/news82047372.html. |
“Proceedings of the Acoustics 2012 Nantes Conference,” Apr. 23-27, 2012, Nantes, France, pp. 278-282. |
Sony New Release: <http://www.sony.net/SonyInfo/News/Press/201010/10-137E/index.html>. |
Number | Date | Country | |
---|---|---|---|
20160002070 A1 | Jan 2016 | US |
Number | Date | Country | |
---|---|---|---|
62020088 | Jul 2014 | US | |
62154672 | Apr 2015 | US |