Acoustically bound microfluidic bubble crystals

Removal of surface-attached micro- and nanobubbles by ultrasonic cavitation in microfluidics

Surface-attached micro- and nanobubbles are known for their resistance to external forces. This study experimentally and theoretically investigates their response to strong ultrasonic fields. Surface-attached micro- and nanobubbles with contact radii from 2 μm to 20 μm are generated in a microchannel and exposed to ultrasound through a vibrating glass substrate. At a driving frequency over 200 kHz up to 2 MHz tested, no significant response from the micro- and nanobubbles is observed. By contrast, at 100 kHz-200 kHz, ultrasonic cavitation bubbles appear in the microchannel and migrate toward the surface micro- and nanobubbles. Then the surface micro- and nanobubbles merge with the ultrasonic cavitation bubbles, detach from the substrate, and become free gaseous nuclei susceptible to further cavitation. Notably, the removal process leaves no observable residue. Theoretical analysis suggests that the directional migration of cavitation bubbles is driven by mutual acoustic radiation forces. This work demonstrates that ultrasonic fields can effectively remove surface micro- and nanobubbles, transforming them into free gaseous cavitation nuclei.

Acoustic manipulation of multi-body structures and dynamics

Sound can exert forces on objects of any material and shape. This has made the contactless manipulation of objects by intense ultrasound a fascinating area of research with wide-ranging applications. While much is understood for acoustic forcing of individual objects, sound-mediated interactions among multiple objects at close range gives rise to a rich set of structures and dynamics that are less explored and have been emerging as a frontier for research. We introduce the basic mechanisms giving rise to sound-mediated interactions among rigid as well as deformable particles, focusing on the regime where the particles' size and spacing are much smaller than the sound wavelength. The interplay of secondary acoustic scattering, Bjerknes forces, and micro-streaming is discussed and the role of particle shape is highlighted. Furthermore, we present recent advances in characterizing non-conservative and non-pairwise additive contributions to the particle interactions, along with instabilities and active fluctuations. These excitations emerge at sufficiently strong sound energy density and can act as an effective temperature in otherwise athermal systems.

Manipulation with sound and vibration: A review on the micromanipulation system based on sub-MHz acoustic waves

Manipulation of micro-objects have been playing an essential role in biochemical analysis or clinical diagnostics. Among the diverse technologies for micromanipulation, acoustic methods show the advantages of good biocompatibility, wide tunability, a label-free and contactless manner. Thus, acoustic micromanipulations have been widely exploited in micro-analysis systems. In this article, we reviewed the acoustic micromanipulation systems that were actuated by sub-MHz acoustic waves. In contrast to the high-frequency range, the acoustic microsystems operating at sub-MHz acoustic frequency are more accessible, whose acoustic sources are at low cost and even available from daily acoustic devices (e.g. buzzers, speakers, piezoelectric plates). The broad availability, with the addition of the advantages of acoustic micromanipulation, make sub-MHz microsystems promising for a variety of biomedical applications. Here, we review recent progresses in sub-MHz acoustic micromanipulation technologies, focusing on their applications in biomedical fields. These technologies are based on the basic acoustic phenomenon, such as cavitation, acoustic radiation force, and acoustic streaming. And categorized by their applications, we introduce these systems for mixing, pumping and droplet generation, separation and enrichment, patterning, rotation, propulsion and actuation. The diverse applications of these systems hold great promise for a wide range of enhancements in biomedicines and attract increasing interest for further investigation.

Amplification of Acoustic Forces Using Microbubble Arrays Enables Manipulation of Centimeter-Scale Objects

Manipulation of macroscale objects by sound is fundamentally limited by the wavelength and object size. Resonant subwavelength scatterers such as bubbles can decouple these requirements, but typically the forces are weak. Here we show that patterning bubbles into arrays leads to geometric amplification of the scattering forces, enabling the precise assembly and manipulation of cm-scale objects. We rotate a 1 cm object continuously or position it with 15  μm accuracy, using sound with a 50 cm wavelength. The results are described well by a theoretical model. Our results lay the foundation for using secondary Bjerknes forces in the controlled organization and manipulation of macroscale structures.

Bubble-enhanced ultrasonic microfluidic chip for rapid DNA fragmentation

DNA fragmentation is an essential process in developing genetic sequencing strategies, genetic research, as well as for the diagnosis of diseases with a genetic signature like cancer. Efficient on-chip DNA fragmentation protocols would be beneficial to process integration and open new opportunities for microfluidics in genetic applications. Here we present an acoustic microfluidic chip comprising an array of ultrasound-actuated microbubbles located at dedicated positions adjacent to a channel containing the DNA sample solution. The efficiency of the on-chip DNA fragmentation process arises mainly from tensile forces generated by acoustic streaming near the oscillating bubble interfaces, as well as a synergistic effect of streaming stress and ultrasonic cavitation. Acoustic microstreaming and the pressure distribution in the DNA channel were assessed by finite element simulation. We characterized the bubble-enhanced effect by measuring gene fragment size distributions with respect to different ultrasound parameters. For optimized on-chip conditions, purified lambda (λ) DNA (48.5 kbp) could be disrupted to fragments with an average size of 2 kbp after 30 s and down to 300 bp after 90 s. Mouse genomic DNA (1.4 kbp) fragmentation size decreased to 500 bp in 30 s and reduced further to 250 bp in 90 s. Bubble-induced fragmentation was more than 3 times faster than without bubbles. On-chip performance and process yield were found to be comparable to a sophisticated high-end commercial system. In this view, our new bubble-enhanced microfluidic approach is a promising tool for current and next generation sequencing platforms with high efficiency and good capacity. Moreover, the availability of an efficient on-chip DNA fragmentation process opens perspectives for implementing full molecular protocols on a single microfluidic platform.

Spectrally wide acoustic frequency combs generated using oscillations of polydisperse gas bubble clusters in liquids

Acoustic frequency combs leverage unique properties of the optical frequency comb technology in high-precision measurements and innovative sensing in optically inaccessible environments such as under water, under ground, or inside living organisms. Because acoustic combs with wide spectra would be required for many of these applications but techniques of their generation have not yet been developed, here we propose an approach to the creation of spectrally wide acoustic combs using oscillations of polydisperse gas bubble clusters in liquids. By means of numerical simulations, we demonstrate that clusters consisting of bubbles with precisely controlled sizes can produce wide acoustic spectra composed of equally spaced coherent peaks. We show that under typical experimental conditions, bubble clusters remain stable over time, which is required for a reliable recording of comb signals. We also demonstrate that the spectral composition of combs can be tuned by adjusting the number and size of bubbles in a cluster.

Feedback-controlled microbubble generator producing one million monodisperse bubbles per second

Monodisperse lipid-coated microbubbles are a promising route to unlock the full potential of ultrasound contrast agents for medical diagnosis and therapy. Here, we present a stand-alone lab-on-a-chip instrument that allows microbubbles to be formed with high monodispersity at high production rates. Key to maintaining a long-term stable, controlled, and safe operation of the microfluidic device with full control over the output size distribution is an optical transmission-based measurement technique that provides real-time information on the production rate and bubble size. We feed the data into a feedback loop and demonstrate that this system can control the on-chip bubble radius (2.5 μm-20 μm) and the production rate up to 10⁶ bubbles/s. The freshly formed phospholipid-coated bubbles stabilize after their formation to a size approximately two times smaller than their initial on-chip bubble size without loss of monodispersity. The feedback control technique allows for full control over the size distribution of the agent and can aid the development of microfluidic platforms operated by non-specialist end users.

Modeling Acoustic Cavitation Using a Pressure-Based Algorithm for Polytropic Fluids

A fully coupled pressure-based algorithm and finite-volume framework for the simulation of the acoustic cavitation of bubbles in polytropic gas-liquid systems is proposed. The algorithm is based on a conservative finite-volume discretization with collocated variable arrangement, in which the discretized governing equations are solved in a single linear system of equations for pressure and velocity. Density is described by the polytropic Noble-Abel stiffened-gas model and the interface between the interacting bulk phases is captured by a state-of-the-art algebraic Volume-of-Fluid (VOF) method. The new numerical algorithm is validated using representative test-cases of the interaction of acoustic waves with the gas-liquid interface as well as pressure-driven bubble dynamics in infinite and confined domains, showing excellent agreement of the results obtained with the proposed algorithm compared to linear acoustic theory, the Gilmore model and high-fidelity experiments.

Acoustic levitation with optimized reflective metamaterials

The simplest and most commonly used acoustic levitator is comprised of a transmitter and an opposing reflecting surface. This type of device, however, is only able to levitate objects along one direction, at distances multiple of half of a wavelength. In this work, we show how a customised reflective acoustic metamaterial enables the levitation of multiple particles, not necessarily on a line and with arbitrary mutual distances, starting with a generic input wave. We establish a heuristic optimisation technique for the design of the metamaterial, where the local height of the surface is used to introduce delay patterns to the reflected signals. Our method stands for any type and number of sources, spatial resolution of the metamaterial and system's variables (i.e. source position, phase and amplitude, metamaterial's geometry, relative position of the levitation points, etc.). Finally, we explore how the strength of multiple levitation points changes with their relative distance, demonstrating sub-wavelength field control over levitating polystyrene beads into various configurations.

Asymmetricity and sign reversal of secondary Bjerknes force from strong nonlinear coupling in cavitation bubble pairs

Most of the current applications of acoustic cavitation use bubble clusters that exhibit multibubble dynamics. This necessitates a complete understanding of the mutual nonlinear coupling between individual bubbles. In this study, strong nonlinear coupling is investigated in bubble pairs which is the simplest case of a bubble-cluster. This leads to the derivation of a more comprehensive set of coupled Keller-Miksis equations (KMEs) that contain nonlinear coupling terms of higher order. The governing KMEs take into account the convective contribution that stems from the Navier-Stokes equation. The system of KMEs is numerically solved for acoustically excited bubble pairs. It is shown that the higher-order corrections are important in the estimation of secondary Bjerknes force for closely spaced bubbles. Further, asymmetricity is witnessed in both magnitude and sign reversal of the secondary Bjerknes force in weak, regular, and strong acoustic fields. The obtained results are examined in the light of published scientific literature. It is expected that the findings reported in this paper may have implications in industries where there is a requirement to have a control on cavitation and its effects.

Acoustic manipulation of particles in a cylindrical cavity: Theoretical and experimental study on the effects of boundary conditions

Precise manipulation of microparticles in microchannels is a primary technique for numerous lab-on-a-chip bioengineering research and applications, as it determines the chip's functions and analytical results. Acoustic manipulation, using the acoustic radiation force, is a compact, versatile and contactless manipulation technique, which can be easily integrated with other microfluidic components. It is our main purpose to report the effect of boundary condition of a cylindrical microfluidic cavity on the acoustic particles' manipulation. A device consisting of a cylindrical cavity in a silicon wafer with three kinds of top boundary conditions (rigid, soft, and imperfect rigid boundary) has been built. The corresponding distributions of acoustic radiation force are analyzed analytically and numerically. Experiments are performed with 2.5 μm radius polystyrene microspheres in the cavity covered by three reflective layers (340 μm-thick glass, 400 μm-thick PDMS, and 660 μm-thick glass film), respectively, which specify the three different boundary conditions at the top of the cavity. It is demonstrated that the boundary condition of a cavity influences the acoustic radiation force and the stable positions of particles, and this is in agreement with the theoretical predictions. Thus, the effects of boundary conditions need to be considered for precise acoustic manipulation.

Effects of Nonlinear Propagation of Focused Ultrasound on the Stable Cavitation of a Single Bubble

Many biomedical applications such as ultrasonic targeted drug delivery, gene therapy, and molecular imaging entail the problems of manipulating microbubbles by means of a high-intensity focused ultrasound (HIFU) pressure field; namely stable cavitation. In high-intensity acoustic field, bubbles demonstrate translational instability, the well-known erratic dancing motion, which is caused by shape oscillations of the bubbles that are excited by their volume oscillations. The literature of bubble dynamics in the HIFU field is mainly centered on experiments, lacking a systematic study to determine the threshold for shape oscillations and translational motion. In this work, we extend the existing multiphysics mathematical modeling platform on bubble dynamics for taking account of (1) the liquid compressibility which allows us to apply a high-intensity acoustic field; (2) the mutual interactions of volume pulsation, shape modes, and translational motion; as well as (3) the effects of nonlinearity, diffraction, and absorption of HIFU to incorporate the acoustic nonlinearity due to wave kinematics or medium—all in one model. The effects of acoustic nonlinearity on the radial pulsations, axisymmetric modes of shape oscillations, and translational motion of a bubble, subjected to resonance and off-resonance excitation and various acoustic pressure, are examined. The results reveal the importance of considering all the involved harmonics and wave distortion in the bubble dynamics, to accurately predict the oscillations, translational trajectories, and the threshold for inertial (unstable) cavitation. This result is of interest for understanding the bubble dynamical behaviors observed experimentally in the HIFU field.

Shape Oscillation of a Single Microbubble in an Ultrasound Field

The shape oscillation of a single two-dimensional nitrogen microbubble in an ultrasound field is numerically investigated. The Navier–Stokes equations are solved by using the finite-volume method combined with the volume-of-fluid model. The numerical results are in good accordance with experimental and theoretical results reported. According to the analyses of the shape oscillation process, the bubble deformation period is twice the driving acoustic pressure period and the shape oscillation is mainly caused by the change of interface velocity. The vortexes produced due to velocity variations lead to the deformation of the bubble interface.

Acoustofluidic Measurements on Polymer-Coated Microbubbles: Primary and Secondary Bjerknes Forces

The acoustically-driven dynamics of isolated particle-like objects in microfluidic environments is a well-characterised phenomenon, which has been the subject of many studies. Conversely, very few acoustofluidic researchers looked at coated microbubbles, despite their widespread use in diagnostic imaging and the need for a precise characterisation of their acoustically-driven behaviour, underpinning therapeutic applications. The main reason is that microbubbles behave differently, due to their larger compressibility, exhibiting much stronger interactions with the unperturbed acoustic field (primary Bjerknes forces) or with other bubbles (secondary Bjerknes forces). In this paper, we study the translational dynamics of commercially-available polymer-coated microbubbles in a standing-wave acoustofluidic device. At increasing acoustic driving pressures, we measure acoustic forces on isolated bubbles, quantify bubble-bubble interaction forces during doublet formation and study the occurrence of sub-wavelength structures during aggregation. We present a dynamic characterisation of microbubble compressibility with acoustic pressure, highlighting a threshold pressure below which bubbles can be treated as uncoated. Thanks to benchmarking measurements under a scanning electron microscope, we interpret this threshold as the onset of buckling, providing a quantitative measurement of this parameter at the single-bubble level. For acoustofluidic applications, our results highlight the limitations of treating microbubbles as a special case of solid particles. Our findings will impact applications where knowing the buckling pressure of coated microbubbles has a key role, like diagnostics and drug delivery.

Sound-mediated stable configurations for polystyrene particles

Here we report an experimental observation of the self-organization effect of polystyrene particles formed by acoustically induced interaction forces. Two types of stable configurations are observed experimentally: one is mechanically equilibrium and featured by nonzero interparticle separations, and the other corresponds to a closely packed assembly, which is created by strong attractions among the aggregated particles. For the former case involving two or three particles, the most probable interparticle separations (counted for numerous independent initial arrangements) agree well with the theoretical predictions. For the latter case, the number of the final stable configurations grows with the particle number, and the occurrence probability of each configuration is interpreted by a simple geometric model.

Acoustic force measurements on polymer-coated microbubbles in a microfluidic device

This work presents an acoustofluidic device for manipulating coated microbubbles, designed for the simultaneous use of optical and acoustical tweezers. A comprehensive characterization of the acoustic pressure in the device is presented, obtained by the synergic use of different techniques in the range of acoustic frequencies where visual observations showed aggregation of polymer-coated microbubbles. In absence of bubbles, the combined use of laser vibrometry and finite element modelling supported a non-invasive measurement of the acoustic pressure and an enhanced understanding of the system resonances. Calibrated holographic optical tweezers were used for direct measurements of the acoustic forces acting on an isolated microbubble, at low driving pressures, and to confirm the spatial distribution of the acoustic field. This allowed quantitative acoustic pressure measurements by particle tracking, using polystyrene beads, and an evaluation of the related uncertainties. This process facilitated the extension of tracking to microbubbles, which have a negative acoustophoretic contrast factor, allowing acoustic force measurements on bubbles at higher pressures than optical tweezers, highlighting four peaks in the acoustic response of the device. Results and methodologies are relevant to acoustofluidic applications requiring a precise characterization of the acoustic field and, in general, to biomedical applications with microbubbles or deformable particles.

Acoustic streaming produced by a cylindrical bubble undergoing volume and translational oscillations in a microfluidic channel

A theoretical model is developed for acoustic streaming generated by a cylindrical bubble confined in a fluid channel between two planar elastic walls. The bubble is assumed to undergo volume and translational oscillations. The volume oscillation is caused by an imposed acoustic pressure field and generates the bulk scattered wave in the fluid gap and Lamb-type surface waves propagating along the fluid-wall interfaces. The translational oscillation is induced by the velocity field of an external sound source such as another bubble or an oscillatory fluid flow. The acoustic streaming is assumed to result from the interaction of the volume and the translational modes of the bubble oscillations. The general solutions for the linear equations of fluid motion and the equations of acoustic streaming are calculated with no restrictions on the ratio between the viscous penetration depth and the bubble size. Approximate solutions for the limit of low viscosity are provided as well. Simulations of streamline patterns show that the geometry of the streaming resembles flows generated by a source dipole, while the vortex orientation is governed by the driving frequency, bubble size, and the distance of the bubble from the source of translational excitation. Experimental verification of the developed theory is performed using data for streaming generated by bubble pairs.

Effect of surface waves on the secondary Bjerknes force experienced by bubbles in a microfluidic channel

An analytical expression is derived for the secondary Bjerknes force experienced by two cylindrical bubbles in a microfluidic channel with planar elastic walls. The derived expression takes into account that the bubbles generate two types of scattered acoustic waves: bulk waves that propagate in the fluid gap with the speed of sound and Lamb-type surface waves that propagate at the fluid-wall interfaces with a much lower speed than that of the bulk waves. It is shown that the surface waves cause the bubbles to form a bound pair in which the equilibrium interbubble distance is determined by the wavelength of the surface waves, which is much smaller than the acoustic wavelength. Comparison of theoretical and experimental results demonstrates good agreement.

Governing Principles of Alginate Microparticle Synthesis with Centrifugal Forces

A controlled synthesis of polymeric particles is becoming increasingly important because of emerging applications ranging from medical diagnostics to self-assembly. Centrifugal synthesis of hydrogel microparticles is a promising method, combining rapid particle synthesis and the ease of manufacturing with readily available laboratory equipment. This method utilizes centrifugal forces to extrude an aqueous polymer solution, sodium alginate (NaALG) through a nozzle. The extruded solution forms droplets that quickly cross-link upon contact with aqueous calcium chloride (CaCl2) solution to form hydrogel particles. The size distribution of hydrogel particles is dictated by the pinch-off behavior of the extruded solution through a balance of inertial, viscous, and surface tension stresses. We identify the parameters dictating the particle size and provide a numerical correlation predicting the average particle size. Furthermore, we create a phase map identifying different pinch-off regimes (dripping without satellites, dripping with satellites, and jetting), explaining the corresponding particle size distributions, and present scaling arguments predicting the transition between regimes. By shedding light on the underlying physics, this study enables the rational design and operation of particle synthesis by centrifugal forces.

Lamb-type waves generated by a cylindrical bubble oscillating between two planar elastic walls

The volume oscillation of a cylindrical bubble in a microfluidic channel with planar elastic walls is studied. Analytical solutions are found for the bulk scattered wave propagating in the fluid gap and the surface waves of Lamb-type propagating at the fluid-solid interfaces. This type of surface wave has not yet been described theoretically. A dispersion equation for the Lamb-type waves is derived, which allows one to evaluate the wave speed for different values of the channel height h. It is shown that for h<λ_t, where λ_t is the wavelength of the transverse wave in the walls, the speed of the Lamb-type waves decreases with decreasing h, while for h on the order of or greater than λ_t, their speed tends to the Scholte wave speed. The solutions for the wave fields in the elastic walls and in the fluid are derived using the Hankel transforms. Numerical simulations are carried out to study the effect of the surface waves on the dynamics of a bubble confined between two elastic walls. It is shown that its resonance frequency can be up to 50% higher than the resonance frequency of a similar bubble confined between two rigid walls.

Exploring bubble oscillation and mass transfer enhancement in acoustic-assisted liquid-liquid extraction with a microfluidic device

We investigated bubble oscillation and its induced enhancement of mass transfer in a liquid-liquid extraction process with an acoustically-driven, bubble-based microfluidic device. The oscillation of individually trapped bubbles, of known sizes, in microchannels was studied at both a fixed frequency, and over a range of frequencies. Resonant frequencies were analytically identified and were found to be in agreement with the experimental observations. The acoustic streaming induced by the bubble oscillation was identified as the cause of this enhanced extraction. Experiments extracting Rhodanmine B from an aqueous phase (DI water) to an organic phase (1-octanol) were performed to determine the relationship between extraction efficiency and applied acoustic power. The enhanced efficiency in mass transport via these acoustic-energy-assisted processes was confirmed by comparisons against a pure diffusion-based process.

A high-power ultrasonic microreactor and its application in gas-liquid mass transfer intensification

The combination of ultrasound and microreactor is an emerging and promising area, but the report of designing high-power ultrasonic microreactor (USMR) is still limited. This work presents a robust, high-power and highly efficient USMR by directly coupling a microreactor plate with a Langevin-type transducer. The USMR is designed as a longitudinal half wavelength resonator, for which the antinode plane of the highest sound intensity is located at the microreactor. According to one dimension design theory, numerical simulation and impedance analysis, a USMR with a maximum power of 100 W and a resonance frequency of 20 kHz was built. The strong and uniform sound field in the USMR was then applied to intensify gas-liquid mass transfer of slug flow in a microfluidic channel. Non-inertial cavitation with multiple surface wave oscillation was excited on the slug bubbles, enhancing the overall mass transfer coefficient by 3.3-5.7 times.

Meshless bubble filter using ultrasound for extracorporeal circulation and its effect on blood

A bubble filter with no mesh structure for extracorporeal circulation using ultrasound was developed. Hemolysis was evaluated by measuring free hemoglobin (FHb). FHb in 120 mL of bovine blood was measured in acoustic standing-wave fields. With a sound pressure amplitude of 60 kPa at driving frequencies of 1 MHz, 500 kHz and 27 kHz for 15 min. FHb values were 641.6, 2575 and 8903 mg/dL, respectively. Thus, hemolysis was inhibited with higher driving frequencies when the same sound pressure amplitude was applied. An ultrasound bubble filter with a resonance frequency of 1 MHz was designed. The filtering characteristics of the flowing microbubbles were investigated with a circulation system using bovine blood with a flow rate of 5.0 L/min. Approximately 99.1% of microbubbles were filtered with 250 kPa and a flow of 5.0 L/min. Hemolysis decreased as the sound pressure decreased; FHb values were 225.8 and 490.7 mg/dL when using 150 and 200 kPa, respectively.

Stable tridimensional bubble clusters in multi-bubble sonoluminescence (MBSL)

In the present work, stable clusters made of multiple sonoluminescent bubbles are experimentally and theoretically studied. Argon bubbles were acoustically generated and trapped using bi-frequency driving within a cylindrical chamber filled with a sulfuric acid aqueous solution (SA85w/w). The intensity of the acoustic pressure field was strong enough to sustain, during several minutes, a large number of positionally and spatially fixed (without pseudo-orbits) sonoluminescent bubbles over an ellipsoidally-shaped tridimensional array. The dimensions of the ellipsoids were studied as a function of the amplitude of the applied low-frequency acoustic pressure (PAc(LF)) and the static pressure in the fluid (P0). In order to explain the size and shape of the bubble clusters, we performed a series of numerical simulations of the hydrodynamic forces acting over the bubbles. In both cases the observed experimental behavior was in excellent agreement with the numerical results. The simulations revealed that the positionally stable region, mainly determined by the null primary Bjerknes force (F→Bj), is defined as the outer perimeter of an axisymmetric ellipsoidal cluster centered in the acoustic field antinode. The role of the high-frequency component of the pressure field and the influence of the secondary Bjerknes force are discussed. We also investigate the effect of a change in the concentration of dissolved gas on the positional and spatial instabilities through the cluster dimensions. The experimental and numerical results presented in this paper are potentially useful for further understanding and modeling numerous current research topics regarding multi-bubble phenomena, e.g. forces acting on the bubbles in multi-frequency acoustic fields, transient acoustic cavitation, bubble interactions, structure formation processes, atomic and molecular emissions of equal bubbles and nonlinear or unsteady acoustic pressure fields in bubbly media.

Self-organizing microfluidic crystals

We consider how to design a microfluidic system in which suspended particles spontaneously order into flowing crystals when driven by external pressure. Via theory and numerics, we find that particle-particle hydrodynamic interactions drive self-organization under suitable conditions of particle morphology and geometric confinement. Small clusters of asymmetric "tadpole" particles, strongly confined in one direction and weakly confined in another, spontaneously order in a direction perpendicular to the external flow, forming one dimensional lattices. Large suspensions of tadpoles exhibit strong density heterogeneities and form aggregates. By rationally tailoring particle shape, we tame this aggregation and achieve formation of large two-dimensional crystals.