Numerical simulation of particle motion in an ultrasound field using the lattice Boltzmann model

In this paper we investigate the motion of small particles suspended in a fluid through which an ultrasound field is propagating. The application of the lattice Boltzmann model to this problem is considered using a two dimensional model. Particles in an ultrasound field are observed to move with a mean particle motion. Further, the time-averaged force on a fixed cylinder is computed and found to be in good agreement with a theoretical expression for the radiation force. Simulations are performed with a single particle, although the approach can equally be applied for a larger number of particles.     

Dielectrophoresis of nanocolloids: A molecular dynamics study

Dielectrophoresis (DEP), the motion of polarizable particles in non-uniform electric fields, has become an important tool for the transport, separation, and characterization of microparticles in biomedical and nanoelectronics research. In this article we present, to our knowledge, the first molecular dynamics simulations of DEP of nanometer-sized colloidal particles. We introduce a simplified model for a polarizable nanoparticle, consisting of a large charged macroion and oppositely charged microions, in an explicit solvent. The model is then used to study DEP motion of the particle at different combinations of temperature and electric field strength. In accord with linear response theory, the particle drift velocities are shown to be proportional to the DEP force. Analysis of the colloid DEP mobility shows a clear time dependence, demonstrating the variation of friction under non-equilibrium. The time dependence of the mobility further results in an apparent weak variation of the DEP displacements with temperature.     

Spatial confinement of ultrasonic force fields in microfluidic channels

We demonstrate and investigate multiple localized ultrasonic manipulation functions in series in microfluidic chips. The manipulation functions are based on spatially separated and confined ultrasonic primary radiation force fields, obtained by local matching of the resonance condition of the microfluidic channel. The channel segments are remotely actuated by the use of frequency-specific external transducers with refracting wedges placed on top of the chips. The force field in each channel segment is characterized by the use of micrometer-resolution particle image velocimetry (micro-PIV). The confinement of the ultrasonic fields during single- or dual-segment actuation, as well as the cross-talk between two adjacent fields, is characterized and quantified. Our results show that the field confinement typically scales with the acoustic wavelength, and that the cross-talk is insignificant between adjacent fields. The goal is to define design strategies for implementing several spatially separated ultrasonic manipulation functions in series for use in advanced particle or cell handling and processing applications. One such proof-of-concept application is demonstrated, where flow-through-mode operation of a chip with flow splitting elements is used for two-dimensional pre-alignment and addressable merging of particle tracks.     

Microfoam formation in a capillary

The ultrasound-induced formation of bubble clusters may be of interest as a therapeutic means. If the clusters behave as one entity, i.e., one mega-bubble, its ultrasonic manipulation towards a boundary is straightforward and quick. If the clusters can be forced to accumulate to a microfoam, entire vessels might be blocked on purpose using an ultrasound contrast agent and a sound source.In this paper, we analyse how ultrasound contrast agent clusters are formed in a capillary and what happens to the clusters if sonication is continued, using continuous driving frequencies in the range 1-10 MHz. Furthermore, we show high-speed camera footage of microbubble clustering phenomena.We observed the following stages of microfoam formation within a dense population of microbubbles before ultrasound arrival. After the sonication started, contrast microbubbles collided, forming small clusters, owing to secondary radiation forces. These clusters coalesced within the space of a quarter of the ultrasonic wavelength, owing to primary radiation forces. The resulting microfoams translated in the direction of the ultrasound field, hitting the capillary wall, also owing to primary radiation forces.We have demonstrated that as soon as the bubble clusters are formed and as long as they are in the sound field, they behave as one entity. At our acoustic settings, it takes seconds to force the bubble clusters to positions approximately a quarter wavelength apart. It also just takes seconds to drive the clusters towards the capillary wall.Subjecting an ultrasound contrast agent of given concentration to a continuous low-amplitude signal makes it cluster to a microfoam of known position and known size, allowing for sonic manipulation.     

Effect of adhesion force between crack planes on subharmonic and DC responses in nonlinear ultrasound

Subharmonic and DC responses in nonlinear ultrasound have been expected as a possible means of detecting closed cracks. Recently, it has been reported that subharmonics in a closed crack markedly increases above a certain input wave amplitude. Such a phenomenon is called "threshold behavior". However, the mechanism of threshold behavior has yet to be elucidated. To clarify this, we introduced adhesion force as a short-range force into the previous analytical model, which expresses the nonlinear contact vibrations of crack planes with intense ultrasound and provides a DC displacement as an approximation of the subharmonic response. Consequently, upward convex curves of displacement against input wave amplitude above the threshold were reproduced for the first time. The validity of the derived analytical solution is discussed by comparison with experimentally observed subharmonics.     

Dynamic acoustic radiation force acting on cylindrical shells: theory and simulations

An object placed in an acoustic field is known to experience a force due to the transfer of momentum from the wave to the object itself. This force is known to be steady when the incident field is considered to be continuous with constant amplitude. One may define the dynamic (oscillatory) radiation force for a continuous wave-field whose intensity varies slowly with time. This paper extends the theory of the dynamic acoustic radiation force resulting from an amplitude-modulated progressive plane wave-field incident on solid cylinders to the case of solid cylindrical shells with particular emphasis on their thickness and contents of their hollow regions. A new factor corresponding to the dynamic radiation force is defined as Y(d) and stands for the dynamic radiation force per unit energy density and unit cross sectional surface. The results of numerical calculations are presented, indicating the ways in which the form of the dynamic radiation force function curves are affected by variations in the material mechanical parameters and by changes in the interior fluid inside the shell's hollow region. It was shown that the dynamic radiation force function Y(d) deviates from the static radiation force function for progressive waves Y(p) when the modulation frequency increases. These results indicate that the theory presented here is broader than the existing theory on cylinders.     

A numerical investigation of the resonance of gas-filled microbubbles: resonance dependence on acoustic pressure amplitude

The general Keller-Herring equation for free gas bubbles is augmented by specific terms to describe the elasticity, viscosity and thickness of the encapsulating shell in ultrasound contrast agent microbubbles. A numerical investigation that analyses the acoustic backscatter from bubbles is employed to identify resonance frequencies that can be compared, for increasing driving pressure amplitude, with linear approximations obtained via analytical considerations. Calculations for bubbles of the size employed in diagnostic ultrasound, between 2 and 6 mum diameter, that are immersed in water and blood and exposed to monochromatic insonation, causing the bubbles to undergo stable cavitation, reveal that the resonance frequency diverges from the linear approximation as the pressure amplitude is increased. The shift in resonance, to lower frequency values, is found to be more pronounced for larger bubbles with the calculated value differing by up to 40% from the linear approximation. The results of this simulation might be potentially useful in preparation of formulations of ultrasound contrast agents with the specifically desired features, such as for instance resonance frequency.     

驻波声场中单分散细颗粒的相互作用特性

利用外加声场促进悬浮在气相中的细颗粒发生相互作用,进而引起颗粒的碰撞和凝并,使得颗粒平均粒径增大、数目浓度降低,是控制细颗粒排放的重要技术途径.为探究驻波声场中单分散细颗粒的相互作用,建立包含曳力、重力、声尾流效应的颗粒相互作用模型,采用四阶经典龙格-库塔算法和二阶隐式亚当斯插值算法对模型进行求解.将数值模拟得到的颗粒声波夹带速度和相互作用过程与相应的解析解和实验结果进行对比,验证模型的准确性.进而研究颗粒初始条件和直径对相互作用特性的影响.结果表明,初始时刻颗粒中心连线越接近声波波动方向、颗粒位置越接近波腹点,颗粒间的声尾流效应就越强,颗粒发生碰撞所需要的时间就越短.研究还发现,颗粒直径对颗粒相互作用的影响取决于初始时刻颗粒中心连线偏离声波波动方向的程度.当偏离较小时,颗粒直径越大,颗粒发生碰撞所需要的时间越短;当偏离很大时,直径较小的颗粒能够发生碰撞,而直径较大的颗粒则无法发生碰撞.     

Magnetic resonance imaging of acoustic streaming: absorption coefficient and acoustic field shape estimation

In this study, magnetic resonance imaging (MRI) is used to visualize acoustic streaming in liquids. A single-shot spin echo sequence (HASTE) with a saturation band perpendicular to the acoustic beam permits the acquisition of an instantaneous image of the flow due to the application of ultrasound. An average acoustic streaming velocity can be estimated from the MR images, from which the ultrasonic absorption coefficient and the bulk viscosity of different glycerol-water mixtures can be deduced. In the same way, this MRI method could be used to assess the acoustic field and time-average power of ultrasonic transducers in water (or other liquids with known physical properties), after calibration of a geometrical parameter that is dependent on the experimental setup.     

Radiation force exerted on nanometer size non-resonant Kerr particle by a tightly focused Gaussian beam

We calculate the radiation force that is exerted by a focused continuous-wave Gaussian beam of wavelength λ on a non-absorbing nonlinear particle of radius a ? 50λ/π. The refractive index of the mechanically-rigid particle is proportional to the incident intensity according to the electro-optic Kerr effect. The force consists of two components representing the contributions of the electromagnetic field gradient and the light scattered by the Kerr particle. The focused intensity distribution is determined using expressions for the six electromagnetic components that are corrected to the fifth order in the numerical aperture (NA) of the focusing objective lens. We found that for particles with a < λ/21.28, the trapping force is dominated by the gradient force and the axial trapping force is symmetric about the geometrical focus. The two contributions are comparable with larger particles and the axial trapping force becomes asymmetric with its zero location displaced away from the focus and towards the beam propagation direction. We study the trapping force behavior versus incident beam power, NA, λ, and relative refractive index between the surrounding liquid and the particle. We also examine the confinement of a Kerr particle that exhibits Brownian motion in a focused beam. Numerical results show that the Kerr effect increases the trapping force strength and significantly improves the confinement of Brownian particles.     

Collective motion of polar active particles on a sphere

Collective motion of active particles with polar alignment is investigated on a sphere. We discussed the factors that affect particle swarm motion and define an order parameter that can show the degree of particle swarm motion. In the model, we added a polar alignment strength, along with Gaussian curvature, affecting particles swarm motion. We find that when the force exceeds a certain limit, the order parameter will decrease with the increase of the force. Combined with our definition of order parameter and observation of the model, the reason is that particles begin to move side by side under the influence of polar forces. In addition, the effects of velocity, rotational diffusion coefficient, and packing fraction on particle swarm motion are discussed. It is found that the rotational diffusion coefficient and the packing fraction have a great influence on the clustering motion of particles, while the velocity has little influence on the clustering motion of particles.     

Spin gauge fields: From Berry phase to topological spin transport and Hall effects

The paper examines the emergence of gauge fields during the evolution of a particle with a spin that is described by a matrix Hamiltonian with n different eigenvalues. It is shown that by introducing a spin gauge field a particle with a spin can be described as a spin multiplet of scalar particles situated in a non-Abelian pure gauge (forceless) field U (n). As the result, one can create a theory of particle evolution that is gauge-invariant with regards to the group Uⁿ (1). Due to this, in the adiabatic (Abelian) approximation the spin gauge field is an analogue of n electromagnetic fields U (1) on the extended phase space of the particle. These fields are force ones, and the forces of their action enter the particle motion equations that are derived in the paper in the general form. The motion equations describe the topological spin transport, pumping, and splitting. The Berry phase is represented in this theory analogously to the Dirac phase of a particle in an electromagnetic field. Due to the analogy with the electromagnetic field, the theory becomes natural in the four-dimensional form. Besides the general theory, the article considers a number of important particular examples, both known and new.     

Correlation in the velocity of a Brownian particle induced by frictional anisotropy and magnetic field

We study the motion of charged Brownian particles in an external magnetic field. It is found that a correlation appears between the components of particle velocity in the case of anisotropic friction, approaching asymptotically zero in the stationary limit. If magnetic field is smaller compared to the critical value, determined by frictional anisotropy, the relaxation of the correlation is non-oscillating in time. However, in a larger magnetic field this relaxation becomes oscillating. The phenomenon is related to the statistical dependence of the components of transformed random force caused by the simultaneous influence of magnetic field and anisotropic dissipation.     

Modeling and experimental validation of hydrodynamics in an ultrasonic batch reactor

Simulation of hydrodynamics in ultrasonic batch reactor containing immobilized enzymes as catalyst is done. A transducer with variable power and constant frequency (24 kHz) is taken as source of ultrasound (US). Simulation comprises two steps. In first step, acoustic pressure field is simulated and in second step effect of this field on particle trajectories is simulated. Simulation results are compared with experimentally determined particle trajectories using PIV Lab (particle image velocimetry). Effect of varying ultrasonic power, positioning and number of ultrasonic sources on particle trajectories is studied. It is observed that catalyst particles tend to orientate according to pattern of acoustic pressure field. An increase in ultrasonic power increases particle velocity and also brings more particles into motion. Simulation results are found to be in agreement with experimentally determined data.     

Real-Time Measurements and Modelling on Dynamic Behaviour of SonoVue Bubbles Based on Light Scattering Technology

The dynamic behaviour of SonoVue microbubbles, a new generation ultrasound contrast agent, is investigated in real time with light scattering method. Highly diluted SonoVue microbubbles are injected into a diluted gel made of xanthan gum and water. The responses of individual SonoVue bubbles to driven ultrasound pulses are measured. Both linear and nonlinear bubble oscillations are observed and the results suggest that SonoVue microbubbles can generate strong nonlinear responses. By fitting the experimental data of individual bubble responses with Sarkar＇s model, the shell coating parameter of the bubbles and dilatational viscosity is estimated to be 7.0 nm-s-Pa.     

Motion of a group of microparticles in a viscoelastic medium under the action of acoustic radiation force

We theoretically and experimentally substantiate the method of detecting microcalcifications in mammary gland tissue. Calcium salts accumulate in soft tissues, primarily forming clusters of individual microparticles. We study the motion of solid microparticles distributed in a viscoelastic medium. Displacement of particles is caused by the radiation force occurring as a consequence of energy scattering and absorption of an ultrasound beam focused in the particle region. The radiation force acts over the course of 200 μs, after which the medium with distributed particles relaxes to the initial state. Motion of the medium is tracked by the cross-correlation method with short probing pulses following at a frequency of 5 kHz. The presence of solid microparticles leads to a change in the character of motion of the medium after pulsed ultrasound action. The amplitude and duration of displacements increases in comparison to the homogeneous medium, and the motion character itself becomes significantly complicated.     

Frequency dependence of the acoustic radiation force acting on absorbing cylindrical shells

The frequency dependence of the radiation force function Y(p) for absorbing cylindrical shells suspended in an inviscid fluid in a plane incident sound field is analysed, in relation to the thickness and the content of their interior hollow region. The theory is modified to include the effect of hysteresis type absorption of compressional and shear waves in the material. The results of numerical calculations are presented for two viscoelastic (lucite and phenolic polymer) materials, with the hollow region filled with water or air indicating how damping and change of the interior fluid inside the shell's hollow region affect the acoustic radiation force. The acoustic radiation force acting on cylindrical lucite shells immersed in a high density fluid (in this case mercury) and filled with water in their hollow region, is also studied.     

Mode-switching: A new technique for electronically varying the agglomeration position in an acoustic particle manipulator

Acoustic radiation forces offer a means of manipulating particles within a fluid. Much interest in recent years has focussed on the use of radiation forces in microfluidic (or “lab on a chip”) devices. Such devices are well matched to the use of ultrasonic standing waves in which the resonant dimensions of the chamber are smaller than the ultrasonic wavelength in use. However, such devices have typically been limited to moving particles to one or two predetermined planes, whose positions are determined by acoustic pressure nodes/anti-nodes set up in the ultrasonic standing wave. In most cases devices have been designed to move particles to either the centre or (more recently) the side of a flow channel using ultrasonic frequencies that produce a half or quarter wavelength over the channel, respectively.It is demonstrated here that by rapidly switching back and forth between half and quarter wavelength frequencies - mode-switching - a new agglomeration position is established that permits beads to be brought to any arbitrary point between the half and quarter-wave nodes. This new agglomeration position is effectively a position of stable equilibrium. This has many potential applications, particularly in cell sorting and manipulation. It should also enable precise control of agglomeration position to be maintained regardless of manufacturing tolerances, temperature variations, fluid medium characteristics and particle concentration.     

A Modeling Study on Particle Dispersion in Wall-Bounded Turbulent Flows

Three physical mechanisms which may affect dispersion of particle's motion in wall-bounded turbulent flows, including the effects of turbulence, wall roughness in particle-wall collisions, and inter-particle collisions, are numerically investigated in this study. Parametric studies with different wall roughness extents and with different mass loading ratios of particles are performed in fully developed channel flows with the Eulerian-Lagrangian approach. A low-Reynolds-number $k-\epsilon$ turbulence model is applied for the solution of the carrier-flow field, while the deterministic Lagrangian method together with binary-collision hard-sphere model is applied for the solution of particle motion. It is shown that the mechanism of inter-particle collisions should be taken into account in the modeling except for the flows laden with sufficiently low mass loading ratios of particles. Influences of wall roughness on particle dispersion due to particle-wall collisions are found to be considerable in the bounded particle-laden flow. Since the investigated particles are associated with large Stokes numbers, i.e., larger than $\mathcal{O}(1)$, in the test problem, the effects of turbulence on particle dispersion are much less considerable, as expected, in comparison with another two physical mechanisms investigated in the study.