Multiscale interactions of liquid, bubbles and solid phases in ultrasonic fields revealed by multiphysics modelling and ultrafast X-ray imaging

The volume of fluid (VOF) and continuous surface force (CSF) methods were used to develop a bubble dynamics model for the simulation of bubble oscillation and implosion dynamics under ultrasound. The model was calibrated and validated by the X-ray image data acquired by ultrafast synchrotron X-ray. Coupled bubble interactions with bulk graphite and freely moving particles were also simulated based on the validated model. Simulation and experiments quantified the surface instability developed along the bubble surface under the influence of ultrasound pressure fields. Once the surface instability exceeds a certain amplitude, bubble implosion occurs, creating shock waves and highly deformed, irregular gas-liquid boundaries and smaller bubble fragments. Bubble implosion can produce cyclic impulsive stresses sufficient enough to cause µs fatigue exfoliation of graphite layers. Bubble-particle interaction simulations reveal the underlying mechanisms for efficient particle dispersion or particle wrapping which are all strongly related to the oscillation dynamics of the bubbles and the particle surface properties.

Modeling and Analysis of the Two-Dimensional Axisymmetric Acoustofluidic Fields in the Probe-Type and Substrate-Type Ultrasonic Micro/Nano Manipulation Systems

The probe-type and substrate-type ultrasonic micro/nano manipulation systems have proven to be two kinds of powerful tools for manipulating micro/nanoscale materials. Numerical simulations of acoustofluidic fields in these two kinds of systems can not only be used to explain and analyze the physical mechanisms of experimental phenomena, but also provide guidelines for optimization of device parameters and working conditions. However, in-depth quantitative study and analysis of acoustofluidic fields in the two ultrasonic micro/nano manipulation systems have scarcely been reported. In this paper, based on the finite element method (FEM), we numerically investigated the two-dimensional (2D) axisymmetric acoustofluidic fields in the probe-type and substrate-type ultrasonic micro/nano manipulation systems by the perturbation method (PM) and Reynolds stress method (RSM), respectively. Through comparing the simulation results computed by the two methods and the experimental verifications, the feasibility and reasonability of the two methods in simulating the acoustofluidic fields in these two ultrasonic micro/nano manipulation systems have been validated. Moreover, the effects of device parameters and working conditions on the acoustofluidic fields are clarified by the simulation results and qualitatively verified by the experiments.

Numerical simulation of single bubble dynamics under acoustic travelling waves

The objective of this paper is to apply CLSVOF method to investigate the single bubble dynamics in acoustic travelling waves. The Naiver-Stokes equation considering the acoustic radiation force is proposed and validated to capture the bubble behaviors. And the CLSVOF method, which can capture the continuous geometric properties and satisfies mass conservation, is applied in present work. Firstly, the regime map, depending on the dimensionless acoustic pressure amplitude and acoustic wave number, is constructed to present different bubble behaviors. Then, the time evolution of the bubble oscillation is investigated and analyzed. Finally, the effect of the direction and the damping coefficient of acoustic wave propagation on the bubble behavior are also considered. The numerical results show that the bubble presents distinct oscillation types in acoustic travelling waves, namely, volume oscillation, shape oscillation, and splitting oscillation. For the splitting oscillation, the formation of jet, splitting of bubble, and the rebound of sub-bubbles may lead to substantial increase in pressure fluctuations on the boundary. For the shape oscillation, the nodes and antinodes of the acoustic pressure wave contribute to the formation of the "cross shape" of the bubble. It should be noted that the direction of the bubble translation and bubble jet are always towards the direction of wave propagation. In addition, the damping coefficient causes bubble in shape oscillation to be of asymmetry in shape and inequality in size, and delays the splitting process.

A numerical and experimental study of acoustic micromixing in 3D microchannels for lab-on-a-chip devices

This paper reports a numerical and experimental study of acoustic streaming and micromixing in polydimethylsiloxane microchannels. The mixing between two fluids flowing in microchannels was evaluated through the following conditions: (1) using a 28 μm thick β-poly(vinylidene fluoride) (β-PVDF) as a piezoelectric transducer actuated with a 24 Vpp and 40 MHz sinusoidal voltage; (2) using different flow rates. The results suggest that the mixing length increases as the flow rate increases and that the acoustic streaming phenomenon leads to a reduction on the mixing length. The good qualitative agreement between numerical and experimental results is a valuable indicator to predict the mixing performance of microfluidic devices, for improving biological fluid analysis in diagnosis lab-on-a-chip devices.

Improving acoustic streaming effects in fluidic systems by matching SU-8 and polydimethylsiloxane layers

This paper reports the use of acoustic waves for promoting and improving streaming in tridimensional polymethylmethacrylate (PMMA) cuvettes of 15mm width×14mm height×2.5mm thickness. The acoustic waves are generated by a 28μm thick poly(vinylidene fluoride) - PVDF - piezoelectric transducer in its β phase, actuated at its resonance frequency: 40MHz. The acoustic transmission properties of two materials - SU-8 and polydimethylsiloxane (PDMS) - were numerically compared. It was concluded that PDMS inhibits, while SU-8 allows, the transmission of the acoustic waves to the propagation medium. Therefore, by simulating the acoustic transmission properties of different materials, it is possible to preview the acoustic behavior in the fluidic system, which allows the optimization of the best layout design, saving costs and time. This work also presents a comparison between numerical and experimental results of acoustic streaming obtained with that β-PVDF transducer in the movement and in the formation of fluid recirculation in tridimensional closed domains. Differences between the numerical and experimental results are credited to the high sensitivity of acoustic streaming to the experimental conditions and to limitations of the numerical method. The reported study contributes for the improvement of simulation models that can be extremely useful for predicting the acoustic effects of new materials in fluidic devices, as well as for optimizing the transducers and matching layers positioning in a fluidic structure.

Evaluation of the successive approximations method for acoustic streaming numerical simulations

This work evaluates the successive approximations method commonly used to predict acoustic streaming by comparing it with a direct method. The successive approximations method solves both the acoustic wave propagation and acoustic streaming by solving the first and second order Navier-Stokes equations, ignoring the first order convective effects. This method was applied to acoustic streaming in a 2D domain and the results were compared with results from the direct simulation of the Navier-Stokes equations. The velocity results showed qualitative agreement between both methods, which indicates that the successive approximations method can describe the formation of flows with recirculation. However, a large quantitative deviation was observed between the two methods. Further analysis showed that the successive approximation method solution is sensitive to the initial flow field. The direct method showed that the instantaneous flow field changes significantly due to reflections and wave interference. It was also found that convective effects contribute significantly to the wave propagation pattern. These effects must be taken into account when solving the acoustic streaming problems, since it affects the global flow. By adequately calculating the initial condition for first order step, the acoustic streaming prediction by the successive approximations method can be improved significantly.