Cumulative distribution functions associated with bubble-nucleation processes in cavitation

Effect of Organic Ions on The Formation and Collapse of Nanometric Bubbles in Ionic Liquid/Water Solutions: A Molecular Dynamics Study

Molecular dynamics simulation is applied to investigate the effect of two ionic liquids (IL) on the nucleation and growth of (nano)cavities in water under tension and on the cavities' collapse following the release of tension. Simulations of the same phenomena in two pure water samples of different sizes are carried out for comparison. The first IL, i.e., tetra-ethylammonium mesylate ([Tea][Ms]), is relatively hydrophilic and its addition to water at 25 wt % concentration decreases its tendency to nucleate cavities. Apart from quantitative details, cavity formation and collapse are similar to those taking place in water and qualitatively follow the Rayleigh-Plesset (RP) equation. The second IL, i.e., tetrabutyl phosphonium 2,4-dimethylbenzenesulfonate ([P₄₄₄₄][DMBS]), is amphiphilic and forms nanostructured solutions with water. At 25 wt % concentrations, [P₄₄₄₄][DMBS] favors the nucleation of bubbles that tend to form at the interface between water-rich and IL-rich domains. Cavity collapse in [P₄₄₄₄][DMBS]/water solutions are greatly hindered by a shell of ions decorating the interface between the solution and the vapor phase. A similar effect is observed for the equilibration of a population of bubbles of different sizes. The drastic slowing down of the bubbles' relaxation processes suggests ways to produce long-lived nanometric cavities in the liquid phase that could be useful for nanotechnology and drug delivery.

Effects of cavitation on Kármán vortex behind circular-cylinder arrays: A molecular dynamics study

The effects of cavitation on the flow around a circular-cylinder array are studied by using a molecular dynamics simulation. Cavitation significantly affects vortex shedding characteristics. As the cavitation develops, the vibration acting on the cylinders decreases and eventually disappears. The further cavitation development generates a longer vapor region next to the cylinders, and the vortex streets are formed at further positions from the cylinders. The neighboring Kármán vortexes are synchronized in the antiphase in the absence of the cavitation. This synchronization is weakened by the cavitation, and an asymmetric wake mode can be induced. These findings help mechanical designs of fluid machinery that include cylinder arrays.

Corner-transport-upwind lattice Boltzmann model for bubble cavitation

Aiming to study the bubble cavitation problem in quiescent and sheared liquids, a third-order isothermal lattice Boltzmann model that describes a two-dimensional (2D) fluid obeying the van der Waals equation of state, is introduced. The evolution equations for the distribution functions in this off-lattice model with 16 velocities are solved using the corner-transport-upwind (CTU) numerical scheme on large square lattices (up to 6144×6144 nodes). The numerical viscosity and the regularization of the model are discussed for first- and second-order CTU schemes finding that the latter choice allows to obtain a very accurate phase diagram of a nonideal fluid. In a quiescent liquid, the present model allows us to recover the solution of the 2D Rayleigh-Plesset equation for a growing vapor bubble. In a sheared liquid, we investigated the evolution of the total bubble area, the bubble deformation, and the bubble tilt angle, for various values of the shear rate. A linear relation between the dimensionless deformation coefficient D and the capillary number Ca is found at small Ca but with a different factor than in equilibrium liquids. A nonlinear regime is observed for Ca≳0.2.

Ripening kinetics of bubbles: A molecular dynamics study

The ripening kinetics of bubbles is studied by performing molecular dynamics simulations. From the time evolution of a system, the growth rates of individual bubbles are determined. At low temperatures, the system exhibits a t^1/2 law and the growth rate is well described by classical Lifshitz-Slyozov-Wagner (LSW) theory for the reaction-limited case. This is direct evidence that the bubble coarsening at low temperatures is reaction-limited. At high temperatures, although the system exhibits a t^1/3 law, which suggests that it is diffusion-limited, the accuracy of the growth rate is insufficient to determine whether the form is consistent with the prediction of LSW theory for the diffusion-limited case. The gas volume fraction dependence of the coarsening behavior is also studied. Although the behavior of the system at low temperatures has little sensitivity to the gas volume fraction up to 10%, at high temperatures it deviates from the prediction of LSW theory for the diffusion-limited case as the gas volume fraction increases. These results show that the mean-field-like treatment is valid for a reaction-limited system even with a finite volume fraction, while it becomes inappropriate for a diffusion-limited system since classical LSW theory for the diffusion-limited case is valid at the dilute limit.

Tensile Strength of Liquids: Equivalence of Temporal and Spatial Scales in Cavitation

It is well known that strain rate and size effects are both important in material failure, but the relationships between them are poorly understood. To establish this connection, we carry out molecular dynamics (MD) simulations of cavitation in Lennard-Jones and Cu liquids over a very broad range of size and strain rate. These studies confirm that temporal and spatial scales play equivalent roles in the tensile strengths of these two liquids. Predictions based on smallest-scale MD simulations of Cu for larger temporal and spatial scales are consistent with independent simulations, and comparable to experiments on liquid metals. We analyze these results in terms of classical nucleation theory and show that the equivalence arises from the role of both size and strain rate in the nucleation of a daughter phase. Such equivalence is expected to hold for a wide range of materials and processes and to be useful as a predictive bridging tool in multiscale studies.

Spontaneous cavitation in a Lennard-Jones liquid: Molecular dynamics simulation and the van der Waals-Cahn-Hilliard gradient theory

The process of bubble nucleation in a Lennard-Jones (LJ) liquid is studied by molecular dynamics (MD) simulation. The bubble nucleation rate J is determined by the mean life-time method at temperatures above that of the triple point in the region of negative pressures. The results of simulation are compared with classical nucleation theory (CNT) and modified classical nucleation theory (MCNT), in which the work of formation of a critical bubble is determined in the framework of the van der Waals-Cahn-Hilliard gradient theory (GT). It has been found that the values of J obtained in MD simulation systematically exceed the data of CNT, and this excess in the nucleation rate reaches 8-10 orders of magnitude close to the triple point temperature. The results of MCNT are in satisfactory agreement with the data of MD simulation. To describe the properties of vapor-phase nuclei in the framework of GT, an equation of state has been built up which describes stable, metastable and labile regions of LJ fluids. The surface tension of critical bubbles γ has been found from CNT and data of MD simulation as a function of the radius of curvature of the surface of tension R*. The dependence γ(R*) has also been calculated from GT. The Tolman length has been determined, which is negative and in modulus equal to ≈(0.1 - 0.2) σ. The paper discusses the applicability of the Tolman formula to the description of the properties of critical nuclei in nucleation.

Molecular dynamics simulations of bubble nucleation in dark matter detectors

Bubble chambers and droplet detectors used in dosimetry and dark matter particle search experiments use a superheated metastable liquid in which nuclear recoils trigger bubble nucleation. This process is described by the classical heat spike model of F. Seitz [Phys. Fluids (1958-1988) 1, 2 (1958)PFLDAS0031-917110.1063/1.1724333], which uses classical nucleation theory to estimate the amount and the localization of the deposited energy required for bubble formation. Here we report on direct molecular dynamics simulations of heat-spike-induced bubble formation. They allow us to test the nanoscale process described in the classical heat spike model. 40 simulations were performed, each containing about 20 million atoms, which interact by a truncated force-shifted Lennard-Jones potential. We find that the energy per length unit needed for bubble nucleation agrees quite well with theoretical predictions, but the allowed spike length and the required total energy are about twice as large as predicted. This could be explained by the rapid energy diffusion measured in the simulation: contrary to the assumption in the classical model, we observe significantly faster heat diffusion than the bubble formation time scale. Finally we examine α-particle tracks, which are much longer than those of neutrons and potential dark matter particles. Empirically, α events were recently found to result in louder acoustic signals than neutron events. This distinction is crucial for the background rejection in dark matter searches. We show that a large number of individual bubbles can form along an α track, which explains the observed larger acoustic amplitudes.

Simple improvements to classical bubble nucleation models

We revisit classical nucleation theory (CNT) for the homogeneous bubble nucleation rate and improve the classical formula using a correct prefactor in the nucleation rate. Most of the previous theoretical studies have used the constant prefactor determined by the bubble growth due to the evaporation process from the bubble surface. However, the growth of bubbles is also regulated by the thermal conduction, the viscosity, and the inertia of liquid motion. These effects can decrease the prefactor significantly, especially when the liquid pressure is much smaller than the equilibrium one. The deviation in the nucleation rate between the improved formula and the CNT can be as large as several orders of magnitude. Our improved, accurate prefactor and recent advances in molecular dynamics simulations and laboratory experiments for argon bubble nucleation enable us to precisely constrain the free energy barrier for bubble nucleation. Assuming the correction to the CNT free energy is of the functional form suggested by Tolman, the precise evaluations of the free energy barriers suggest the Tolman length is ≃0.3σ independently of the temperature for argon bubble nucleation, where σ is the unit length of the Lennard-Jones potential. With this Tolman correction and our prefactor one gets accurate bubble nucleation rate predictions in the parameter range probed by current experiments and molecular dynamics simulations.

Cavitation in a metallic liquid: homogeneous nucleation and growth of nanovoids

Large-scale molecular dynamics (MD) simulations are performed to investigate homogeneous nucleation and growth of nanovoids during cavitation in liquid Cu. We characterize in detail the atomistic cavitation processes by following the temporal evolution of cavities or voids, analyze the nucleation behavior with the mean first-passage time (MFPT) and survival probability (SP) methods, and discuss the results against classical nucleation theory (CNT), the Tolman equation for surface energy, independent calculation of surface tension via integrating the stress profiles, the Johnson-Mehl-Avrami (JMA) growth law, and the power law for nucleus size distributions. Cavitation in this representative metallic liquid is a high energy barrier Poisson processes, and the steady-state nucleation rates obtained from statistical runs with the MFPT and SP methods are in agreement. The MFPT method also yields the critical nucleus size and the Zeldovich factor. Fitting with the Tolman's equation to the MD simulations yields the surface energy of a planar interface (~0.9 J m⁻²) and the Tolman length (0.4-0.5 Å), and those values are in accord with those from integrating the stress profiles of a planar interface. Independent CNT predictions of the nucleation rate (10(33 - 34) s(-1) m(-3)) and critical size (3-4 Å in radius) are in agreement with the MFPT and SP results. The JMA law can reasonably describe the nucleation and growth process. The size distribution of subcritical nuclei appears to follow a power law with an exponent decreasing with increasing tension owing to coupled nucleation and growth, and that of the supercritical nuclei becomes flattened during further stress relaxation due to void coalescence.

Bubble evolution and properties in homogeneous nucleation simulations

We analyze the properties of naturally formed nanobubbles in Lennard-Jones molecular dynamics simulations of liquid-to-vapor nucleation in the boiling and the cavitation regimes. The large computational volumes provide a realistic environment at unchanging average temperature and liquid pressure, which allows us to accurately measure properties of bubbles from their inception as stable, critically sized bubbles, to their continued growth into the constant speed regime. Bubble gas densities are up to 50% lower than the equilibrium vapor densities at the liquid temperature, yet quite close to the gas equilibrium density at the lower gas temperatures measured in the simulations: The latent heat of transformation results in bubble gas temperatures up to 25% below those of the surrounding bulk liquid. In the case of rapid bubble growth-typical for the cavitation regime-compression of the liquid outside the bubble leads to local temperature increases of up to 5%, likely significant enough to alter the surface tension as well as the local viscosity. The liquid-vapor bubble interface is thinner than expected from planar coexistence simulations by up to 50%. Bubbles near the critical size are extremely nonspherical, yet they quickly become spherical as they grow. The Rayleigh-Plesset description of bubble-growth gives good agreement in the cavitation regime.

Direct simulations of homogeneous bubble nucleation: Agreement with classical nucleation theory and no local hot spots

We present results from direct, large-scale molecular dynamics simulations of homogeneous bubble (liquid-to-vapor) nucleation. The simulations contain half a billion Lennard-Jones atoms and cover up to 56 million time steps. The unprecedented size of the simulated volumes allows us to resolve the nucleation and growth of many bubbles per run in simple direct micro-canonical simulations while the ambient pressure and temperature remain almost perfectly constant. We find bubble nucleation rates which are lower than in most of the previous, smaller simulations. It is widely believed that classical nucleation theory (CNT) generally underestimates bubble nucleation rates by very large factors. However, our measured rates are within two orders of magnitude of CNT predictions; only at very low temperatures does CNT underestimate the nucleation rate significantly. Introducing a small, positive Tolman length leads to very good agreement at all temperatures, as found in our recent vapor-to-liquid nucleation simulations. The critical bubbles sizes derived with the nucleation theorem agree well with the CNT predictions at all temperatures. Local hot spots reported in the literature are not seen: Regions where a bubble nucleation event will occur are not above the average temperature, and no correlation of temperature fluctuations with subsequent bubble formation is seen.

Amyloid fibril disruption by ultrasonic cavitation: nonequilibrium molecular dynamics simulations

We describe the disruption of amyloid fibrils of Alzheimer's amyloid-β peptides by ultrasonic cavitation. For this purpose, we performed nonequilibrium all-atom molecular dynamics simulations with sinusoidal pressure and visualized the process with movies. When the pressure is negative, a bubble is formed, usually at hydrophobic residues in the transmembrane region. Most β-strands maintain their secondary structures in the bubble. When the pressure becomes positive, the bubble collapses, and water molecules crash against the hydrophilic residues in the nontransmembrane region to disrupt the amyloid. Shorter amyloids require longer sonication times for disruption because they do not have enough hydrophobic residues to serve as a nucleus to form a bubble. These results agree with experiments in which monodispersed amyloid fibrils were obtained by ultrasonication.

Single-bubble dynamics in pool boiling of one-component fluids

We numerically investigate the pool boiling of one-component fluids with a focus on the effects of surface wettability on the single-bubble dynamics. We employed the dynamic van der Waals theory [Phys. Rev. E 75, 036304 (2007)], a diffuse-interface model for liquid-vapor flows involving liquid-vapor transition in nonuniform temperature fields. We first perform simulations for bubbles on homogeneous surfaces. We find that an increase in either the contact angle or the surface superheating can enhance the bubble spreading over the heating surface and increase the bubble departure diameter as well and therefore facilitate the transition into film boiling. We then examine the dynamics of bubbles on patterned surfaces, which incorporate the advantages of both hydrophobic and hydrophilic surfaces. The central hydrophobic region increases the thermodynamic probability of bubble nucleation while the surrounding hydrophilic region hinders the continuous bubble spreading by pinning the contact line at the hydrophobic-hydrophilic intersection. This leads to a small bubble departure diameter and therefore prevents the transition from nucleate boiling into film boiling. With the bubble nucleation probability increased and the bubble departure facilitated, the efficiency of heat transfer on such patterned surfaces is highly enhanced, as observed experimentally [Int. J. Heat Mass Transfer 57, 733 (2013)]. In addition, the stick-slip motion of contact line on patterned surfaces is demonstrated in one-component fluids, with the effect weakened by surface superheating.

Toward a molecular theory of homogeneous bubble nucleation: II. Calculation of the number density of critical nuclei and the rate of nucleation

In the present paper, we develop a method to calculate the rate of homogeneous bubble nucleation within a superheated L-J liquid based on the (n,v) equilibrium embryo free energy surface introduced in the first paper (DOI: 10.1021/jp404149n). We express the nucleation rate as the product of the concentration of critical nuclei within the metastable liquid phase and the relevant forward rate coefficient. We calculate the forward rate coefficient of the critical nuclei from their average lifetime as determined from MD simulations of a large number of embryo trajectories initiated from the transitional region of the metastable liquid configuration space. Therefore, the proposed rate coefficient does not rely on any predefined reaction coordinate. In our model, the critical nuclei belong to the region of the configuration space where the committor probability is about one-half, guaranteeing the dynamical relevance of the proposed embryos. One novel characteristic of our approach is that we define a limit for the configuration space of the equilibrium metastable phase and do not include the configurations that have zero committor probability in the nucleation free energy surface. Furthermore, in order to take into account the transitional degrees of freedom of the critical nuclei, we develop a simulation-based approach for rigorously mapping the free energy of the (n,v) equilibrium embryos to the concentration of the critical nuclei within the bulk metastable liquid phase.

Phase diagram and universality of the Lennard-Jones gas-liquid system

The gas-liquid phase transition of the three-dimensional Lennard-Jones particles system is studied by molecular dynamics simulations. The gas and liquid densities in the coexisting state are determined with high accuracy. The critical point is determined by the block density analysis of the Binder parameter with the aid of the law of rectilinear diameter. From the critical behavior of the gas-liquid coexisting density, the critical exponent of the order parameter is estimated to be β = 0.3285(7). Surface tension is estimated from interface broadening behavior due to capillary waves. From the critical behavior of the surface tension, the critical exponent of the correlation length is estimated to be ν = 0.63(4). The obtained values of β and ν are consistent with those of the Ising universality class.