Molecular dynamics simulation of shock-induced microscopic bubble collapse

Characterization of the Dynamic Behavior of Multinanobubble System under Shock Wave Influence

Shockwave-induced changes in nanobubbles cause cavitation erosion and membrane damage but can also be applied to biocarrier transport. Currently, research focuses on single nanobubbles; however, in reality, nanobubbles usually appear as a multibubble system. Therefore, this study proposes a method based on cutting and replicating to construct a multibubble model. This method can be widely applied to molecular dynamics (MD) models and enhance the customization capabilities of MD models. The dynamic behavior of a multinanobubble system with different numbers and arrangements of nanobubbles is investigated with the MD method under the influence of shock waves in a liquid argon system. The study also explores the range of influence between nanobubbles. The results show that in the case of two nanobubbles, when the distance between the bubbles is constant, the smaller the angle between the direction of the shock wave and the line connecting the bubbles, the greater is the influence between nanobubbles, and the moment of collapse of the nanobubbles farther away from the shock wave is slower. When three nanobubbles are arranged with a right offset, after the first bubble collapses, the effect on the other two bubbles is similar to the changes in bubbles when the angle of arrangement is 30° or 60°. Under a different arrangement, the change of shock wave velocity on the nanobubble size only affects its collapse time and contraction collapse rate. When the shock wave with a radian of about 2.87 or greater than 2.87 touches the bubbles, the collapse of the second nanobubble will not be affected.

The impact of low-velocity shock waves on the dynamic behaviour characteristics of nanobubbles

Low-velocity shock wave-induced contraction and expansion of nanobubbles can be applied to biocarriers and microfluidic systems. Although experiments have been conducted to study the application effects, the dynamic behavior characteristics of nanobubbles remain unexplored. In this work, we utilize molecular dynamics (MD) simulations to investigate the dynamic behavior characteristics of nanobubbles influenced by low-velocity shock waves in a liquid argon system. The DBSCAN (Density-Based Spatial Clustering of Applications with Noise) machine learning method is used to calculate the equivalent radius of nanobubbles. Two statistical methods are then utilized to predict the time series changes in the equivalent radius of nanobubbles without rebound shock waves. The piston velocity is analyzed using the bisection method to obtain the critical impact states of the nanobubble. The results show that at the low velocity shock wave (piston velocity of 0.1 km s-1), the shock wave pressure is small, the non-vacuum nanobubbles contract and expand in a circular shape, and the gas particles inside the bubble are not dispersed. In contrast, the vacuum nanobubbles collapse directly. As the shock wave rebounds upon impact, it triggers periodic contraction and expansion of the nanobubbles. The predictions indicate that the equivalent radius will vary within a small range according to the pre-predicted values in the absence of the rebound shock wave. Nanobubbles are present in four critical impact states: dispersed gaps, multiple smaller bubbles, two split bubbles, and a concave bubble.

Effect of Shock-Wave-Mediated Collapse on Nanobubble Characteristics

To enhance the comprehension of the cavitation mechanism and explore its practical use in industrial production, this study developed models involving oxygen, varying bubble radii, and bubble quantities. This study uses molecular dynamics simulations coupled with the momentum mirror method to examine the collapse characteristics of bubbles during ultrasonic cavitation. The investigation uncovers patterns in the fluctuation of the maximum local density of water molecules, released pressure, and temperature. The findings demonstrate that, when oxygen-containing bubbles collapse at identical radii, the local density is notably higher and diminishes more rapidly. Moreover, the changes in the shape exhibit greater regularity. During the bubble collapse, a depression forms on the bubble's surface, coinciding with a notable surge in local density around the depression. As bubble radii and quantities increase, so does the local density along with a concurrent rise in the maximum pressure. Intriguingly, the model demonstrates the lowest pressure at Z = 35 Å accompanied by the emergence of a small crescent-shaped region with a reduced density. Throughout the pressure ascension phase, the rate of the maximum pressure change escalates with an increase in the number of bubbles. Conversely, during the pressure descent phase, the rate of the maximum pressure change diminishes with a growing number of bubbles. However, it is important to note that the maximum pressure does not exhibit a direct correlation with the number of bubbles. Ultimately, this study provides valuable technical guidance and a theoretical foundation for the integration of ultrasonic cavitation in industrial production processes.

Near-Wall Cavitation Effect: A Molecular Dynamics Study

Cavitation has been the subject of abundant studies, but the internal mechanism of cavitation is less well known. In this article, a microlevel near-wall model was established by using LAMMPS to present the process of cavitation effect. The results of molecular dynamics simulation revealed the fluctuation process of the liquid near the wall with the change in pressure. Molecular dynamics was also used to evaluate the void volume fraction and density distribution of the system. The results exhibited that the cavitation process can be divided into two stages: the initial cavitation stage and the rapid growth stage. Based on these results, the effects of wettability and initial system temperature on the near-wall cavitation effect were demonstrated. The results indicated that the hydrophobic near-wall forms a gas layer to weaken the density fluctuation, while the hydrophilic wall is opposite. Increasing the temperature could positively affect molecular motion and cavitation. This work provides a theoretical basis for further exploration of the cavitation effect.