Bubble Evolution Research Articles

Modeling the distribution characteristics of vapor bubbles in cavitating flows

Dispersed vapor bubbles are the dominant rheology in cloud cavitation, and their size distribution is directly associated with cavitation noise and erosion. However, the numerical resolution of large numbers of dispersed bubbles remains a challenge. In this work, we establish a new cavitation model based on the population balance equation (PBE) that can predict the size distribution and spatiotemporal evolution of bubbles within cloud cavitation under different cavitation numbers. An expression for the phase transition source term without empirical parameters is derived based on the bubble size distribution (BSD) function, enabling the coupling of mass transfer in the governing equations with the PBE cavitation model. The cavitation model is solved alongside the Eulerian homogeneous mixture flow. The mass transfer between water and vapor, and the bubble coalescence and breakup under turbulent flows, are modeled to determine the BSD. The numerical model is carefully validated through comparisons with experimental results for cavitation flows on a wedge-shaped flat plate, and good agreement is achieved with respect to the pressure distribution, void fraction, and BSD. This confirms that our proposed cavitation model can accurately predict the void fraction and BSD within the cloud cavitation region.

He bubbles inhibition and abnormal hardening in GH3535 alloy by niobium element addition

GH3535 and 2 wt.% Nb-modified GH3535 (Nb-GH3535) alloys were irradiated with 1.2 MeV He ion at fluences of 1 × 1016, 3 × 1016 and 1 × 1017 ions/cm2 at 750 °C. Microstructural changes, including the formation of He bubbles and superlattice phases, as well as the irradiation-induced hardening, were investigated by TEM and nanoindentation. Notably, Nb-GH3535 exhibited a smaller size and lower number density of He bubbles, especially at higher fluences, indicating the inhibitory effect of Nb on bubble evolution. First-principles calculations revealed that the presence of Nb inhibits the combination of He atoms and vacancies, reducing their diffusion rate and consequently impeding the nucleation and growth of He bubbles. Furthermore, the addition of Nb facilitates the irradiation-induced precipitation of the superlattice phase, leading to a higher degree of total hardening in Nb-GH3535 alloy.

Numerical simulation of the equatorial plasma bubble: the effect of seeding by the vertical winds and random background noise perturbations

A wide variety of small-amplitude waves widely exist in the ionosphere and have significant effects on the evolution of equatorial plasma bubbles. In this paper, we simulated equatorial plasma bubbles (EPB) seeded by vertical neutral wind perturbations with wavelengths of 125 km and 250 km, and compared the morphology characteristics of plasma bubble structures with those under random noise perturbations in the background density. The numerical results showed that both vertical winds and random background noise perturbations can contribute to the growth of plasma bubbles, and the perturbations under additional random background noise can promote the growth of the plasma bubble structures faster. Additionally, several processes of the nonlinear behavior of bifurcated EPB structures, including bifurcation, pinching, and small-scale turbulent structures, were successfully obtained. Our simulation captured supersonic flows within the low-density plasma structures characterized by vertical velocities of about 1.5 km/s, which is consistent with experimental studies found in the literature.

Effect of water head on the permeability of foam-conditioned sands: Experimental and analytical investigation

The water permeability of conditioned soils is one of the most essential properties for Earth Pressure Balance (EPB) tunnelling in coarse-grained soils. Permeability tests are conducted to study the influence of water heads on the permeability of foam-conditioned sands. The initial permeability coefficient of foam-conditioned sands increases with the water head, while the stable permeability coefficient and the initial stable period duration decrease. Meanwhile, a novel analytical model is proposed to estimate the initial permeability coefficient. In this model, the effect of the water head on the initial permeability coefficient is incorporated by calculating void ratios of the foam and effective diameters of foam bubbles under different water pressures. Experimental results are in close agreement with analytical solutions, indicating the excellent performance of the proposed calculation method. In addition, the physical mechanisms of how the water head affects the permeability of foam-conditioned sands are discussed from the contraction and evolution of foam bubbles.

Shrinkage and coalescence behavior of helium bubbles in Fe9Cr1.5W0.4Si F/M steel under the coupling effect of compression stress and temperature

As important irradiation defects, helium bubbles will threaten the stability and reliability of reactor structural materials. Extensive research has been conducted on the dynamic response of helium bubbles under thermal annealing, irradiation, or tension, but the behavior and mechanisms of helium bubble evolution under compression remain to be understood. Here, by adopting in-situ cyclic nanocompression testing on Fe9Cr1.5W0.4Si F/M steel containing helium bubbles at 673 K, we present direct experimental evidence for the shrinkage and coalescence behavior of helium bubbles under the coupling effect of compression stress and temperature for the first time. Unlike the stress-free condition where the bubble size shows an increasing trend, compressive stress causes the bubble to shrink. Moreover, the large-sized bubbles show a greater change in size than small-sized bubbles. At 673 K, there is no thermal motion of the bubbles in F/M steel, and it is the deformation and mutual squeezing of bubbles under compression that causes the coalescence of tightly arranged bubbles. Furthermore, helium bubbles can act as dislocation sources to emit dislocations under compression, and the emitted dislocations can in turn envelope around the helium bubbles, thus hindering bubble coalescence and leading to the formation of deformed bubbles.

Formation of Hierarchical Zn/N‐doped Carbon Hollow Nanofibers towards Dendrite‐Free Zn Metal Anodes

AbstractRechargeable aqueous Zn metal batteries are considered as promising next‐generation energy storage device because of their intrinsic safety and low cost. Nonetheless, the poor lifespan stemming from uncontrollable dendrite growth and inevitable parasitic reactions of metallic Zn impede their large‐scale application. Herein, the formation of the hierarchical lotus root‐like Zn/N‐doped carbon hollow nanofibers‐based paper is reported decorated with interconnected Zn/N‐doped carbon nanocages (denoted as LRZCF@ZCC) as a freestanding host with multifunctionality for Zn metal anode. Of note, the interconnected porous structure mitigates the internal stress to ensure robust structure. Besides, the concave surface of inner LRZCF@ZCC channel facilitates the nucleation of Zn2+ to geometrically guide the Zn deposition. Meanwhile, the well‐distributed metallic Zn and N‐doped carbon species can serve as the zincophilic sites to reduce the nucleation overpotential and avoid the evolution of hydrogen bubbles. As expected, the LRZCF@ZCC exhibits low polarization for the controllable Zn deposition with small voltage hysteresis. Both the symmetric cell and the full cell using this unique host demonstrate long‐term cycling duration and enhanced rate capability.

Bubble Management for Electrolytic Water Splitting by Surface Engineering: A Review.

During electrocatalytic water splitting, the management of bubbles possesses great importance to reduce the overpotential and improve the stability of the electrode. Bubble evolution is accomplished by nucleation, growth, and detachment. The expanding nucleation sites, decreasing bubble size, and timely detachment of bubbles from the electrode surface are key factors in bubble management. Recently, the surface engineering of electrodes has emerged as a promising strategy for bubble management in practical water splitting due to its reliability and efficiency. In this review, we start with a discussion of the bubble behavior on the electrodes during water splitting. Then we summarize recent progress in the management of bubbles from the perspective of surface physical (electrocatalytic surface morphology) and surface chemical (surface composition) considerations, focusing on the surface texture design, three-dimensional construction, wettability coating technology, and functional group modification. Finally, we present the principles of bubble management, followed by an insightful perspective and critical challenges for further development.

Three-dimensional phase-field simulation of intergranular bubble evolution in UO2 during irradiation

In this study, a three-dimensional phase-field model was developed to investigate the evolution of intergranular bubbles during irradiation. The model considered several grain structures and multiple fission gas bubbles, with the microstructure evolution being governed by improved dynamic equations. Verification of the phase-field model showed good agreement with the theory in describing the properties of grain boundary (GB). The study examined the dependency of bubble percolation on bubble shape, fission rate, and average grain size. The simulation results revealed that the changes in these three factors have significant effects on the shape of GB bubbles, the percolation rate, and the connectivity threshold of GB bubbles. The bubble connection threshold was obtained from the curve of triple junction coverage versus grain boundary coverage. The results indicate that the connectivity threshold increases with the higher fission rate and lager average grain size. This investigation into bubble percolation during irradiation provides valuable insights for further study on fission gas release in UO2 fuel within nuclear reactors.

Vibrations Analysis of Bubble Evolution in Liquids of Varying Physical Properties

In recent years there is an attempt to control the gas stirring intensity in metal-making ladles with the aid of vibration measurements. Understanding better the induced vibrations in two-phase flows can substantially improve the existing models for gas stirring control. In this work, highly sensitive accelerometers were used for the vibration measurements in a liquid metal alloy; Sn–40 wt pctBi alloy at 200 °C and water at 20 °C. The examination of the liquids was conducted in the ladle mockup integrated into the Liquid Metal Model for Steel Casting facility at Helmholtz-Zentrum Dresden Rossendorf. Single bubbles were generated in the respective liquids by controlled argon injection at low flow rates in the range of 0.01 to 0.15 NL min−1 through a single nozzle installed at the bottom of the ladle. Obtained results demonstrate differences between the induced vibrations in the examined liquids in terms of the magnitude of the root mean square values of vibration amplitude and the shape of the resulting curves with increasing flow rate. Furthermore, continuous wavelet transform reveals variations in the duration and vibrational frequency of the evolved bubble phenomena. The findings suggest that differences in the physical properties of the examined liquids result in variations in the vibrations induced during bubble evolution.

Effect of the Temperature on Interfacial Properties of CO2/H2 Mixtures Contacting with Brine and Hydrophilic Silica by Molecular Dynamics Simulations.

Underground H2 storage (UHS) is a promising technology to achieve large-scale, long-term H2 storage. Using CO2 as a cushion gas to maintain the pressure of the reservoir and withdraw stored H2 in the saline aquifer simultaneously enables the implementation of UHS and underground CO2 storage (UCS). The difference in the molecular properties of CO2 and H2 leads to distinct interfacial behavior when in contact with the brine and rock, thereby affecting the flow patterns and trapping mechanisms of gases in geological formations. Accurate prediction of the interfacial properties of CO2, H2, and the mixtures when interacting with brine and rock is crucial to minimizing the uncertainties in UHS and UCS projects. In this study, molecular dynamics (MD) simulations are performed to predict the interfacial tension, surface excess, bubble evolution, and contact angle of CO2, H2, and the mixtures at 10 MPa and 300-400 K. The MD results show that the interaction of CO2 with H2O and hydrophilic silica is considerably stronger than that of H2. The interfacial tension reduces linearly with the temperature in H2-dominated mixture systems, and the surface adsorption of H2 can diminish in a CO2-dominated system or at high-temperature conditions. The hydrophilic silica is more CO2-wet than H2-wet, and the attached CO2 bubble is more easily disconnected. Ions and the temperature play different roles in the contact angle.

Numerical simulation of Marangoni flow around a growing hydrogen bubble on a microelectrode

In this work, a numerical model for Marangoni flow around a growing hydrogen bubble on a microelectrode is developed and validated. The growth of the bubble is implemented using a body-fitted moving mesh method that moves the bubble interface based on the production of hydrogen at the cathode. The dissolved hydrogen that is produced at the cathode diffuses into the bubble. The geometry of the bubble is a spherical cap with a pinned contact line and a decreasing contact angle as the bubble grows larger. There are two cases: (1) a case in which the bubble grows on a microbubble carpet that grows during the evolution of the bubble, and (2) a case in which the carpet thickness is negligible. The purpose of this model is to simulate and investigate thermocapillary Marangoni convection during almost the entire bubble lifetime (nucleation and detachment are not included).We validate the model with data from two experiments reported in literature (Massing et al. 2019; Bashkatov et al. 2022; Bashkatov, 2022). The evolution of the bubble radius, the temperature profile around the bubble, and the profile of the Marangoni velocity, which is tangential to the bubble surface at 5 μm from the bubble, near the lower half of the bubble are accurately predicted. For the case with a bubble carpet, we investigated models with a growing carpet and fixed carpet thickness. We found that for a given bubble size the Marangoni velocity reduces with increasing carpet thickness.Stagnant bubble models are less complicated and computationally cheaper than the growing bubble model. For the case without a carpet, we compare the growing bubble model with two stagnant bubble models: one using constant current and another one using a transient current. A stagnant bubble model is recommended if the growth rate drb/dt<0.05uM where drb/dt is the instantaneous bubble growth rate and uM the average Marangoni velocity at the bubble interface. This is the case at a relatively late stage of the bubble evolution. At the early stage of the bubble evolution, the absence of bubble growth dynamics in the stagnant bubble models appears to have an effect on the velocity profile and the pressure part of the hydrodynamic force.

Influence of water layer thickness on crater volume for nanosecond pulsed laser ablation of stainless steel

Under water laser ablation is a surface texturization method used to form micrometer-sized surface structures. Plasma confinement and cavitation bubble evolution play a critical role during the ablation process and their influence on material removal is strongly tied to liquid layer thickness. To influence the effects of these processes, such that material removal is at its maximum, an optimal layer thickness was found for various laser parameters. Specifically, for nanosecond pulsed laser ablation of stainless steel, however, the relation between layer thickness and volume removal is still unknown. Here, we show the relation between water layer thickness and removed material volume for a nanosecond pulsed laser. Results reveal that volume removal is at its maximum for a 1 mm water layer and drops by a factor of 2 when the layer thickness is increased to 2 mm. A further increase of layer thickness to 3 up to 10 mm shows a negligible effect on volume removal and removed volume amounts are shown to be similar to those obtained in ambient air in this water layer thickness range. This trend echo’s results obtained for nanosecond pulsed silicon ablation. The obtained results identify processing conditions which allow for faster and therefore more cost efficient texturization of stainless steel surfaces in the future.

Effect of grain size on gas bubble evolution in nuclear fuel: Phase-field investigations

Numerous irradiation-induced gas bubbles are created in the nuclear fuel during irradiation, leading to the change of microstructure and the degradation of mechanical and thermal properties. The grain size of fuel is one of the important factors affecting bubble evolution. In current study, we first predict the thermodynamic behaviors of point defects as well as the interplay between vacancy and gas atom in both UO2 and U3Si2 according to ab initio approach. Then, we establish the irradiation-induced bubble phase-field model to investigate the formation and evolution of intra- and inter-granular gas bubbles. The effects of fission rate and temperature on the evolutions of bubble morphologies in UO2 and U3Si2 have been revealed. Especially, a comparison of porosities under different grain sizes is examined and analyzed. To understand the thermal conductivity as functions of grain size and porosity, the heat transfer capability of U3Si2 is evaluated.

Identification of Penetration and Disturbance Dynamo Electric Fields and Their Effects on the Generation of Equatorial Plasma Bubbles

AbstractA very challenging task in ionospheric studies is to determine the separate contributions of penetration and disturbance dynamo processes in the generation of equatorial plasma drift during magnetic storms. In this study, we analyze the ion drift measured by the Communications/Navigation Outage Forecasting System satellite during the magnetic storm on 15–16 July 2012. A unique feature of this storm is the exceptionally long period of continuous southward interplanetary magnetic field (IMF) for 32 hr. The storm‐induced net change of the meridional/vertical ion drift, the difference drift between the storm time and the quiet time, is derived during the storm main phase and the first 20 hr of the recovery phase with southward IMF. The difference drift during the recovery phase cannot be explained by the disturbance dynamo effect alone. A new method is used to separate the drifts caused by the penetration and disturbance dynamo processes. The penetration drift is represented by an empirical pattern of penetration electric field and depends on the IMF magnitude, and the disturbance dynamo drift is obtained by subtracting the penetration drift from the measured difference drift. The derived disturbance dynamo drift is in good agreement with previous statistical pattern. This is the first effort to identify the separate contributions of the penetration and disturbance dynamo processes to the total drift from observed data. The results have important implications in identifying storm‐time penetration and disturbance dynamo electric fields and their effects on the generation and evolution of plasma bubbles.

Magnetic Storm Effects on the Occurrence and Characteristics of Plasma Bubbles

AbstractThe Communications/Navigation Outages Forecast System satellite mission was designed to investigate the ionospheric conditions that lead to the formation of irregularities. Here, we have studied the effect of magnetic storms on the formation and evolution of plasma bubbles during the satellite's lifetime (2008–2015). During this period encompassing solar minimum and maximum conditions, many magnetic storms of varying intensity developed, producing a unique and rich data set of 248 storms (14 intense, 69 moderate, and 165 weak) that occurred during the same timeframe to examine the role of external magnetospheric drivers in the production and dynamics of equatorial plasma bubbles. We have used the Planar Langmuir Probe and Ion Velocity Meter instruments to elucidate the role of magnetic storm intensity on the bubble's depth, internal speed, width, occurrence, and lifetime. The pre‐reversal enhancement (PRE) tends to increase during the main phase and when BZ is southward. New bubbles occur during large excursions of the PRE value. The bubble lifetime extends and remains active during the main and part of the recovery phase. The plasma velocity within the bubbles increases and typically becomes over 100 m/s during significant PRE and BZ negative times. The depth of bubbles reaches values close to 100% during intense storms. In general, the intensity of the storms seems to control and augment the plasma bubbles' depth, width, and internal velocity.

Equatorial Plasma Bubbles: Zonal Thermosphere Wind Influence

There are a number of the theoretical studies pointing to the key role of the zonal thermosphere winds in the equatorial plasma bubble (EPB) generation and evolution. However, there are insufficient observational data to confirm the relationship between these phenomena. To study this relationship, the detailed comparative and correlation analysis of the LT-variations of the EPB occurrence probability and the zonal thermosphere wind velocity was carried out. The data from the EPB observations recorded aboard the ISS-b satellite (~972–1220 km) during the solstice and equinox periods were used. The data from the zonal thermosphere wind velocity observations obtained aboard the CHAMP satellite (~380–450 km) were also used. It was found that these characteristics are similar, when they are compared, and have very strong correlation (R 0.9) in summer, strong correlation (R 0.8) in winter and (R 0.79) at equinox. It was found that in all seasons the delay in the development of the EPB occurrence probability maxima with respect to the maxima of the west wind velocity is 1–3 h. It is in good agreement with the estimation of the time of the seeding perturbation development and the EPB rise to the topside ionosphere altitudes. The results can be considered a new confirmation of the theoretical conclusion (Kudeki model) about the key influence of the zonal west thermosphere winds on the equatorial plasma bubble generation.

Generation and evolution of post-sunset equatorial plasma bubbles in East and Southeast Asia during the July 2022 geomagnetic storm

In the East and Southeast Asia, equatorial plasma bubbles (EPBs) mainly occur at post-sunset in equinoctial months. On 7–8 July 2022, a moderate geomagnetic storm occurred, with minimum Dst of −85 nT. This storm could trigger unseasonal post-sunset EPBs over the East/Southeast Asian sector. Observations from GNSS TEC/scintillation receivers, VHF coherent radar and ionosonde show that on 7 July, a few clusters of EPBs, which caused strong amplitude scintillations and significant TEC fluctuations, occurred over a wide latitude range, from low to middle latitudes up to ∼ 32°N. The EPB irregularities at middle latitudes drifted westward, with velocities of about 97 m/s. Before the generation of EPBs, apparent background ionospheric disturbances were detected. An abrupt rise of ionospheric F layer, with upward drifts of about 33 m/s was observed near sunset at low latitudes. It is suggested that the moderate geomagnetic storm was the main reason for the generation of unseasonal post-sunset EPBs. The prompt penetration of eastward electric field caused by the geomagnetic storm could lead to the abrupt rise of F layer near sunset and promote the growth of Rayleigh-Taylor instability. The background ionospheric disturbance could act as a seeding role for the unseasonal EPBs.

Behaviors of cavitation bubbles driven by high-intensity ultrasound

In a multi-bubble system, the bubble behavior is modulated by the primary acoustic field and the secondary acoustic field. To explore the translational motion of bubbles in cavitation liquids containing high-concentration cavitation nuclei, evolutions of bubbles are recorded by a high-speed camera, and translational trajectories of several representative bubbles are traced. It is found that translational motion of bubbles is always accompanied by the fragmentation and coalescence of bubbles, and for bubbles smaller than 10 μm, the possibility of bubble coalescence is enhanced when the spacing of bubbles is less than 30 μm. The measured signals and their spectra show the presence of strong negative pressure, broadband noise, and various harmonics, which implies that multiple interactions of bubbles appear in the region of high-intensity cavitation. Due to the strong coupling effect, the interaction between bubbles is random. A simplified triple-bubble model is developed to explore the interaction patterns of bubbles affected by the surrounding bubbles. Patterns of bubble interaction, such as attraction, repulsion, stable spacing, and rebound of bubbles, can be predicted by the theoretical analysis, and the obtained results are in good agreement with experimental observations. Mass exchange between the liquid and bubbles as well as absorption in the cavitation nuclei also plays an important role in multi-bubble cavitation, which may account for the weakening of the radial oscillations of bubbles.

NiCo layered double hydroxides/NiFe layered double hydroxides composite (NiCo-LDH/NiFe-LDH) towards efficient oxygen evolution in different water matrices

Engineering robust non-noble metal electrocatalysts towards efficient impure water (e.g., seawater, wastewater) oxidation is a prospective approach to achieve carbon neutrality via accelerating green hydrogen energy development. Herein, a NiCo layered double hydroxides (LDH)/NiFe LDH composite (NiCo-LDH/NiFe-LDH) was developed for oxygen evolution reaction (OER) via a hydrothermal process-electrodeposition method. The optimal NiCo-LDH/NiFe-LDH-30 composite only needed an overpotential (η) of 240 mV to drive 100 mA/cm2 in alkalized freshwater, with a low Tafel slope of 16.6 mV/dec and good stability for over 90 h. Further analyses suggested that the strong interface interaction between NiCo-LDH and NiFe-LDH accelerated the oxygen gas bubble evolution and boosted interfacial charge transfer, and the formed built-in electric field and higher oxidation state species (metal oxyhydroxides) contributed to the high intrinsic catalytic activity. The NiCo-LDH/NiFe-LDH-30 composite also held excellent OER activities in different impure water environments, including alkaline 0.5 M NaCl solution (η100 = 333 mV), alkaline lake water (η100 = 345 mV), and alkaline wastewater treatment plant (WWTP) effluent (η100 = 320 mV). More importantly, the potential effects of Cl− and CO32− in impure water were revealed during the OER process. This work elaborates on the role of built-in electric field and the strong coupling interaction in composite catalysts, which pave the way for the design of cost-effective catalysts with excellent adaptability in different water environments.

The design of Pd-containing high-entropy alloys and their hardening behavior under He ion irradiation

Imparting additional chemical heterogeneities of high-entropy alloys (HEAs) is considered as a novel strategy to improve their irradiation resistance on dislocation loop growth; however, there is still a lack of knowledge about the effects of this strategy on bubble evolution and irradiation-induced hardening. In this study, a series of Pd-containing FeCrNiCo HEAs with proven potential for chemical short-range order were designed and irradiated by 400 keV He ions to the fluences of 5 × 1016 cm−2, 1 × 1017 cm−2, and 5 × 1017 cm−2 at 723 K. Microscopic characterization demonstrates a remarkably inhibited growth of dislocation loops and promoted bubble coarsening in the Pd-containing HEAs compared to the HEA without Pd alloying. The underlying mechanism is that although the alloying of Pd can suppress the evolution of dislocation loops through severe local lattice distortion, it also significantly reduces the vacancy formation energy of the entire HEA system, leading to the encouraged bubble evolution. Pd-containing HEAs exhibit improved resistance to irradiation-induced hardening at low doses, but this advantage is unexpectedly absent with increasing dose. The hardening mechanism is explained by identifying the loop-dominated hardening at low dose and the bubble-induced hardening which makes an increased contribution at high dose. In addition, the Ni segregation around bubbles and loop punching mechanism are also discussed to provide essential evidence. These insights can help to shed light on the effects of Pd alloying on bubble evolution and irradiation-induced hardening behavior in HEAs, and thus contribute to the design of irradiation-resistant HEAs.