
==== Front
Ultrason Sonochem
Ultrason Sonochem
Ultrasonics Sonochemistry
1350-4177
1873-2828
Elsevier

S1350-4177(24)00312-2
10.1016/j.ultsonch.2024.107064
107064
Original Research Article
Bubble shape instability of acoustic cavitation in molten metal used in ultrasonic casting
Yamamoto Takuya takuya.yamamoto@omu.ac.jp

Department of Chemical Engineering, Graduate School of Engineering, Osaka Metropolitan University, 1-1, Gakuen-cho, Naka-ku, Sakai, Osaka 599-8531, Japan
13 9 2024
12 2024
13 9 2024
111 10706430 6 2024
28 8 2024
9 9 2024
© 2024 The Author(s)
2024
https://creativecommons.org/licenses/by/4.0/ This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).
In this study, we estimated the equilibrium bubble size of acoustic cavitation in a molten metal, which is basic information in ultrasonic casting. For this, the bubble shape instability of acoustic cavitation in the melt was numerically investigated by solving the Keller–Miksis equation and dynamic equation of the distortion amplitude. The acoustic cavitation bubbles are more stable in aluminum and magnesium melts than in water, and the parametric instability mainly determines the bubble stability at 20–160 kHz in molten metals. However, the afterbounce instability does not significantly affect the bubble stability in molten metals owing to the small number of bubble oscillations after the first rapid compression, and the distortion amplitude cannot grow significantly after the first compression. The bubbles in the melt become more unstable with an increase in the ultrasonic frequency owing to the corresponding increase in the bubble wall velocity. Through this stability analysis, we can estimate that the stable bubble size in the aluminum or magnesium melt is approximately three or four times larger than that in water at the same ultrasonic pressure amplitude.

Keywords

Acoustic cavitation
Keller-Miksis equation
Shape instability
Aluminum melt
Magnesium melt
==== Body
pmc1 Introduction

When an ultrasound is irradiated into a liquid, acoustic cavitation occurs causing nonlinear bubble oscillations [1], [2]. During these oscillations, the bubble is significantly compressed with a significant increase in the temperature inside the bubble. In addition, the acoustic cavitation bubbles cause shock waves and micro-jet formation [3], [4], [5], [6], [7]. In particular, the water vapor in a bubble is thermally decomposed into radicals, causing many chemical reactions [8], [9]. These phenomena have been used to develop many applications, such as ultrasonic emulsification [10], [11], [12], [13], [14], [15], [16], [17], atomization [18], [19], [20], [21], [22], nano-particle fabrication [23], [24], [25], organic chemical decomposition [26], [27], [28], and ultrasonic casting [29], [30], [31], [32]. In this study, we focused on phenomena that occur during ultrasonic casting.

Ultrasonic casting is an ultrasonic application that has been investigated for a long time [32]. It is expected to be used especially for metals with comparably low-temperature melting points, such as light metals [30], [32]. This is because the lifetime of sonotrodes is significantly short for metals with high-temperature melting points [33]. The solidified microstructures of grains or intermetallic compounds are refined by the ultrasound irradiation [34], [35], [36]. The homogenized size and morphology of these microstructures improve the mechanical characteristics of the cast metal. In addition, grain refiners, which are used to refine the grains of the microstructure, are activated by ultrasound irradiation to remove impurities attached to them [37]. Research has been widely conducted, particularly on aluminum and magnesium alloy melts, to improve the material properties. Attempts have been made to elucidate the phenomena occurring during casting by measuring the sound field distribution in the molten metals [38], [39], [40], acoustic streaming [41], and the direct observation of bubble motion and size [42], [43], [44]. Although many phenomena have been clarified in these studies, some remain unclear because they are difficult to observe. For instance, although the bubble size of acoustic cavitation in an Al-Cu alloy melt was measured using synchrotron observation, the bubble size was measured at a location far from the sonotrode tip owing to the limitation of synchrotron observation [43], [45], and the measured bubble size largely varied in the experiment [46]. In addition, the nanoparticles in the melt increased the bubble size in the melt [47]. Under ultrasound conditions, many bubbles of different sizes exist in the liquid. The cavitation bubble is known to grow owing to the bubble coalescence and rectified diffusion [48], [49]. The sufficiently grown bubbles cannot oscillate when the bubble radius exceeds the resonance radius. Therefore, two types of bubbles exist in a liquid: fine oscillating and large unoscillating. Fine oscillating bubbles emit shock waves and form microjets. Therefore, they are important in ultrasonic casting for the fragmentation of intermetallic compounds and primary crystals. However, active and inactive bubbles cannot be identified separately by the direct observation. Researchers have attempted to measure the time variation of the bubble radius at considerably high speeds, even in the melt. Although bubble motion in a melt with a low-temperature melting point has been measured directly using a high-speed camera [50], measuring the equilibrium bubble size is difficult. Hence, whether the measured bubble size corresponds to the equilibrium bubble size of the melt is unclear. In numerical simulation studies, the equilibrium bubble size in an aluminum alloy melt was set to be almost identical to that in water [51], [52], [53], [54], [55], [56], [57]. In addition to ultrasonic casting in an aluminum melt, ultrasound has been used during the casting of magnesium alloy melt [58], [59], [60]. However, the bubble size of acoustic cavitation in a magnesium alloy melt has not been measured at all, and the equilibrium bubble radius is still unclear.

In the case of acoustic cavitation in water, the equilibrium bubble size has been measured accurately using laser techniques [61], [62], [63], [64]. However, this technique cannot be used for molten metals because of their opacity. In another method, the bubble equilibrium size in water was estimated using stability analysis. Many researchers estimated the stable bubble size by solving both the Rayleigh-Plesset equation and dynamic equation of distortion of spherical harmonics [65], [49], [66], [67], [68], [69], [70], [71], [72], [73] and, also, by solving the Keller–Miksis equation [74], [75], [76], [77], [78], [79] instead of the Rayleigh-Plesset equation. This method was also used to investigate the effects of the surfactants [77], liquid viscosity [74], and ultrasonic frequency [79] in water. In addition, the stabilities of acoustic cavitation bubbles in a glycerin solution [78] and sulfuric mixture [73] were investigated. However, the stability of acoustic bubbles in a molten metal has not been investigated, and the equilibrium bubble radius in the molten metal remains unknown.

The information regarding the bubble stability of acoustic cavitation is significantly important for understanding the phenomena occurring in ultrasonic casting, which is in high demand in the material manufacturing industry. In this study, we used the Keller–Miksis equation and the dynamic equation of distortion to investigate bubble shape instability in aluminum and magnesium melts, which are expected to be used in ultrasonic casting, through stability analysis.

2 Numerical analysis

In this study, we solved the Keller–Miksis equation [80] and dynamic equation of distortion of spherical harmonics to simulate the time evolution of the bubble dynamics behavior and small distortions of the bubble shape, respectively. The Keller–Miksis equation was derived from the continuity and Euler equations by including the effect of liquid compressibility, and the dynamic distortion equation was derived from the radial oscillation model equation. The detailed models are presented below.

2.1 The Keller-Miksis equation

The bubble radial dynamics were modeled using the following Keller–Miksis equation [1], [80]:(1) 1-R˙c∞RR¨+32R˙21-R˙c∞=1ρL,∞1+R˙c∞pB-pst-p0+Rc∞ρL,∞dpBdt

where R is the bubble radius, c∞ is the sound speed, ρL∞ is the liquid density, pB is the liquid pressure at the bubble wall, ps is the sinusoidal ultrasonic pressure oscillation, p0 is the atmosphere pressure, and t is time. The liquid pressure at the bubble wall is expressed as(2) pB=pg+pv-2σR-4μR˙R

where pg is the pressure inside the bubble, pv is the vapor pressure, σ is the surface tension, and μ is the dynamic viscosity. The pressure inside the bubble is modeled as(3) pg=p0+2σR0-pvR0R3γ

where γ is the specific heat ratio. The sinusoidal ultrasonic wave is given as(4) pst=-Pasinωt

where Pa is pressure amplitude, and ω is angular frequency.

In this simulation, to avoid divergence in simulation at high pressure amplitude, the bubble wall velocity (R˙) was replaced by the following sound speed in liquid at the bubble wall (cLB) when the bubble wall velocity exceeded the sound speed (c∞) [81]. The sound speed in the liquid at the bubble wall was modeled as(5) cLB=7.15pB+BρLi

where B is the pressure parameter, and ρLi is the liquid density at the bubble wall, which were modeled as in the literature [1], [82]. The sound speed was derived from the modified Tait equation [83].

The physical properties and used parameters are shown in Table 1.Table 1 Physical properties and parameters used in this stability analysis [84], [85].

Properties	Value	Unit	
Aluminum	Magnesium	
Liquid density, ρL,∞	2.35 × 103	1.59 × 103	kgm−3	
External pressure, P0	1.00 × 105	1.00 × 105	Pa	
Liquid sound velocity, c∞	4.65 × 103	4.00 × 103	ms−1	
Liquid kinematic viscosity, ν	4.42 × 10−7	7.86 × 10−6	m2s−1	
Surface tension, σ	8.57 × 10−1	5.59 × 10−1	Nm−1	
Ratio of specific heat, γ	1.40	−	
Sound pressure amplitude, Pa	7.00–15.0 × 104	Pa	
Equilibrium bubble radius, R0	5.00–50.0 × 10−6	m	
Angular frequency, ω	20000, 40000, 80000, 160,000
× 2π	rads−1	

2.2 Dynamic equation for small distortion

As introduced by Plesset [65] and Eller et al. [49], [66], the fluctuating bubble radius is considered as a linear combination of the bubble radius and distortion of spherical harmonics. The dynamic equation for small distortions of the spherical harmonics can be derived by introducing a fluctuating radius into the dynamic equation of the bubble radius. The dynamic equation for a small distortion is described as(6) a¨n+Bnta˙n-Antan=0

where An and Bn are expressed as [70], [71](7) Ant=n-1R¨R-βnσρR3-2νR˙R3n-1n+2+2nn+2n-1δR

(8) Bnt=3R˙R+2νR2n+22n+1-2nn+22δR

where βn is modeled as(9) βn=n-1n+1n+2,

where n is the degree of spherical harmonics. The boundary layer thickness δ was modeled using Brenner’s model [70], [71] as(10) δ=minνω,R2n

The boundary layer thickness was modeled as a diffusion length scale, and a cutoff length, R/2n, was imposed. The time variations of the bubble radius and distortion amplitude were calculated by numerically solving the above Keller–Miksis equation and dynamics equation for small distortions.

After solving these equations, the bubble stability was evaluated using two factors: the magnitude of the maximal eigenvalue of the Floquet transition matrix and the maximal value of the distortion amplitude. The Floquet transition matrix Fn (T) is defined as(11) ant+Ta˙n(t+T)=Fnt+Tanta˙nt

where T is an ultrasonic cycle period, which is defined by(12) T=2πω

The distortion amplitude decays when the magnitude of the maximal eigenvalue of the Floquet transition matrix is less than unity. The parametric instability was evaluated using this index.

The other factor is the maximal value of the distortion amplitude. During the ultrasonic cycle, the bubble fragments when the distortion amplitude exceeds the bubble radius. In particular, after the first bubble compression, the distortion amplitude begins to increase. This type of instability is called the afterbounce instability. The threshold of this instability is calculated as:(13) maxt′|t<t′<t+Tant′Rt′≥1

When the distortion amplitude exceeds the bubble radius, the bubbles are considered unstable.

2.3 Numerical methods

The Keller–Miksis equation and distortion equations are discretized using 4th order Runge-Kutta method. The time step for this simulation (Δt) was set to T/500,000. At frequencies of 20 and 160 kHz, the time steps were 0.10 and 0.0125 ns, respectively. The initial distortion amplitude was set to 1.0 nm. It should be noted that we previously investigated the effect of the initial distortion amplitude on the stability diagram in our previous study [79] and found that the initial distortion amplitude does not change the stability diagram. Bubble stability was investigated with different bubble radii (5.0–50.0 μm at intervals of 0.10 μm) and pressure amplitudes (7.0 × 104–1.5 × 105 Pa at intervals of 200 Pa). The number of sampled data points in one stability diagram was 450×400 for the bubble radius and pressure amplitude. In our previous study, we validated and verified this program [79].

3 Numerical results

3.1 Instability in aluminum melt

Fig. 1 shows the stability diagram of the bubble shape of acoustic cavitation in an aluminum melt under a 20 kHz ultrasound. In this graph, the black and white zones indicate the stable and unstable zones for bubble shape instability, respectively. As explained in Section 2.3, the number of data points was 400 and 450 in the horizontal and vertical directions, respectively. The investigated parameter ranges were divided into 400 × 450 squares. When the bubble shape became stable, the square domain was marked in black. Therefore, the spatial resolution of each graph was 400 × 450 pixels, which was sufficiently large. The resolution of the figure might look coarse and unclear. However, this plot of the stable domain originated from the characteristics of the stable zone. The thin stable and unstable band-shaped domains overlapped alternately. Therefore, a locally stable domain was formed, as shown in this figure. The stability diagram of parametric instability is completely different from that of afterbounce instability. The main stable zone of the parametric instability is located at the bottom-left corner of Fig. 1 (a). This implies a spherical bubble shape corresponding to a small pressure amplitude of ultrasound and small bubble radius. Meanwhile, the afterbounce instability was hard to occur in this parameter range. Therefore, the bubble became aspherical owing to parametric instability. Afterbounce instability does not occur easily in any of the other parameters. Hence, we focused only on parametric instability.Fig. 1 Stability diagram of bubble shape of acoustic cavitation in the aluminum melt at 20 kHz ultrasound when the degree of spherical harmonics is 2: (a) parametric instability and (b) afterbounce instability.

Next, we evaluated the influence of the spherical harmonic degree on the stability diagram. Fig. 2 shows the stability diagram of the parametric instability in an aluminum melt at 20 kHz ultrasound with different degrees of spherical harmonics. Similar to the acoustic cavitation in water [71], [79], the most unstable mode was n = 2. The stable zone expanded slightly with an increase in the degree of spherical harmonics, although a sharp stable band-shaped domain was formed near the boundary between the stable and unstable zones. For example, the boundary at a pressure amplitude of 70000 Pa was located at approximately 30, 35, and 50 μm for n = 2, 4, and 6, respectively. Similar to acoustic cavitation in water, Mathieu toughs are formed near the border of stable and unstable zones. Notably, the tough length was much longer, and the tough shape was much sharper in the aluminum melt than in water. In this section, we have focused on the most unstable mode (n = 2).Fig. 2 Stability diagram of parametric instability in the aluminum melt at 20 kHz ultrasound with different degree of spherical harmonics, n = (a) 2, (b) 3, (c) 4, (d) 5, and (e) 6.

The influence of ultrasonic frequency on bubble shape instability was investigated for an aluminum melt. Fig. 3 shows the stability diagram of the parametric instability in the aluminum melt ultrasound with n = 2 at different ultrasonic frequencies. Similar to parametric instability in water, the stable zone in the stability diagram shrank with an increase in ultrasonic frequency. Surprisingly, the stable zone aligned horizontally at higher ultrasonic frequencies. This result indicates that the bubble instability is not affected by the pressure amplitude, but by the bubble radius. When the ultrasonic frequency increased, the bubbles became unstable, and the slope of the border between the stable and unstable zones became gentle. This gentle slope was more remarkable in the case of the aluminum melt, causing a horizontally aligned stable zone.Fig. 3 Stability diagram of parametric instability in the aluminum melt ultrasound with n = 2 degree with different ultrasonic frequencies: (a) 20, (b) 40, (c) 80, and (d) 160 kHz.

3.2 Instability in magnesium melt

Fig. 4 shows the stability diagram of the bubble shape of the acoustic cavitation in the magnesium melt at 20 kHz ultrasound when the degree of spherical harmonics was 2. The tendency of the unstable zone in the magnesium melt was the same as that in the aluminum melt. The stable zone of the parametric instability for the magnesium melt is completely different from that of the afterbounce instability. The stable zone expands widely in the case of afterbounce instability. Hence, the most unstable mode is determined by the parametric instability. Hereafter, we focused only on the parametric instability. Smaller bubbles and pressure amplitudes were required to maintain spherical bubbles. Fig. 1, Fig. 4 were compared, and the stable zone was found narrower in the case of the magnesium melt for parametric instability. Hence, the acoustic cavitation bubbles in the magnesium melt are more unstable than those in the aluminum melt.Fig. 4 Stability diagram of bubble shape of acoustic cavitation in the magnesium melt at 20 kHz ultrasound when the degree of spherical harmonics is 2: (a) parametric instability and (b) afterbounce instability.

Next, we investigated the influence of the spherical harmonic degree on the stability diagram. Fig. 5 shows a stability diagram of the parametric instability in the magnesium melt at 20 kHz ultrasound with different degrees of spherical harmonics. Similar to aluminum melt, as shown in Fig. 2, the stable zone expanded with an increase in the degree of the spherical harmonics. The most unstable mode was n = 2. Therefore, we focused on the parametric instability with n = 2.Fig. 5 Stability diagram of parametric instability in the magnesium melt at 20 kHz ultrasound with different degree of spherical harmonics, n = (a) 2, (b) 3, (c) 4, (d) 5, and (e) 6.

Fig. 6 shows the stability diagram of the parametric instability in the magnesium melt ultrasound with n = 2 at different ultrasonic frequencies. The tendency of the stable zone in the magnesium melt was the same as that in the aluminum melt. When the ultrasonic frequency increased, the stable zone shrank, and the bubble became unstable. Similar to aluminum melt, the stable zone became horizontally aligned at higher ultrasonic frequencies. In the next section, the mechanism of the above-mentioned phenomena is discussed, and the bubble stability in different liquids is compared.Fig. 6 Stability diagram of parametric instability in the magnesium melt ultrasound with n = 2 with different ultrasonic frequencies: (a) 20, (b) 40, (c) 80, and (d) 160 kHz.

3.3 Discussion

Here, the generation mechanism of bubble shape instability is discussed separately, and bubble instabilities are compared in different liquids.

3.3.1 Generation mechanism of bubble shape instability

In this subsection, the bubble shape instability is discussed in terms of the following three points:• The reason for bubble instability with an increase in ultrasonic frequency.

• The reason for the difficulty of afterbounce instability occurrence in molten metals, although the afterbounce instability is most unstable mode in water at 20 kHz [79].

• The reason for the local stable zone to appear at high pressure amplitudes and large bubble radii.

The phenomenon was investigated in our previous study [79]. The moving velocity and acceleration of bubble wall increase with ultrasonic frequency. Hence, the distortion amplitude increases significantly, and the bubbles become unstable at higher ultrasonic frequencies. For further details, the readers are requested to refer to our previous study [79].

In the previous study, parametric and afterbounce instability occurred simultaneously at 20 kHz in water [79], and the parametric instability became dominant at higher ultrasonic frequencies. Meanwhile, afterbounce instability did not occur in molten metals, even at low ultrasonic frequencies, as shown in Fig. 3, Fig. 6. Second, we have explained the reason for difficulty of afterbounce instability occurrence in molten metal at 20 kHz below.

Fig. 7 shows the time variation of the bubble radius and distortion amplitude with n = 2 at 20 kHz ultrasound with a pressure amplitude of 1.1 × 105 Pa. The bubble radius was 40 μm in the aluminum and magnesium melt. The distortion amplitude increased after the first collapse of the bubble. After collapse, the bubble was bound one or two times before the start of the next ultrasonic cycle. Compared to the bubble oscillations in water at 20 kHz, the number of bounces after the first collapse was small. For reference, the time variations of the bubble radius and distortion amplitude in water are shown in Fig. 8. In this simulation, the pressure amplitude was 9.0 × 104 Pa and the bubble radius was 10 μm. Contrary to the molten metals, small oscillations with high frequency were observed in water after the first collapse of the bubble. During these small bubble oscillations, the distortion amplitude increased significantly and afterbounce instability occurred. However, the number of small oscillations after the first collapse is small in the molten metal. Therefore, afterbounce instability occurrence is difficult in the molten metal. Conversely, the distortion amplitude increased with the number of ultrasonic cycles in the molten metal, as shown in Fig. 7(a). Therefore, parametric instability can occur in molten metals.Fig. 7 Time variation of (a) bubble radius and (b) distortion amplitude with n = 2 at 20 kHz ultrasound with the pressure amplitude of 1.1 × 105 Pa; the bubble radius is 40 μm in the aluminum and the magnesium melts.

Fig. 8 Time variation of (a) bubble radius and (b) distortion amplitude with n = 2 at 20 kHz ultrasound with the pressure amplitude of 9.0 × 104 Pa; the bubble radius is 10 μm in water.

Third, we have discussed the reason for a local stable zone to appear at high pressure amplitudes and large bubble radii below. This phenomenon can be simulated by direct numerical simulations using the volume of the fluid method [86]. Here, we focused on the stability diagram of the aluminum melt under 80 kHz ultrasound, as shown in Fig. 3(c). A long and wide stable domain was observed at approximately R0 = 30–40 μm. Therefore, the time variations of the bubble radius and distortion amplitude were evaluated for the equilibrium radii R0 = 25, 35, and 40 μm at a pressure amplitude of 8.0 × 104 Pa. The bubbles were unstable at R0 = 25 and 40 μm but stable at R0 = 35 μm. Fig. 9 shows the time variation of the (a) bubble radius and (b) distortion amplitude with n = 2 at 80 kHz ultrasound with a pressure amplitude of 8.0 × 104 Pa for the bubble radii of 25, 35, and 40 μm in the aluminum melt. The time at which the bubble was the most compressed shifted later as the bubble size increased. Accordingly, the periods of the bubble and ultrasonic oscillations did not match. Owing to this mismatch, the distortion amplitude did not increase, especially at approximately R0 = 35 μm.Fig. 9 Time variation of (a) bubble radius and (b) distortion amplitude with n = 2 at 80 kHz ultrasound with the pressure amplitude of 8.0 × 104 Pa; the bubble radii are 25, 35 and 40 μm in the aluminum melt.

3.3.2 Comparison of bubble instability between different liquids

As shown in Fig. 3, Fig. 6, the stable zone in the aluminum melt was slightly larger than that in the magnesium melt. These results were compared with the stability diagram of water reported in the literature [70], [71], [72], [74], [77], [79], [87]: the stable zone was much wider, and the bubbles were more stable in aluminum and magnesium melts. For example, when the ultrasonic frequency was 20 kHz and the pressure amplitude was 1.0 × 105 Pa, the bubbles became unstable at ∼6 μm in water and at ∼20–25 μm in the aluminum and magnesium melts. These results indicate that the bubble radius of acoustic cavitation was larger in aluminum and magnesium melts than in water. This result is in good quantitative agreement with experimentally observed bubbles [43]. Fig. 10 shows the time variation of the bubble radius (=10 μm) with different liquids at 20 kHz with a pressure amplitude of 1.0 × 105 Pa. The bubble oscillation amplitude was considerable small in the case of molten metals. This is because the surface tension of molten metals is one order of magnitude higher than that of water. In addition, the surface tension of the aluminum melt is higher than that of the magnesium melt as shown in Table 1. Therefore, the acoustic cavitation bubbles are more stable in the aluminum melt. The distortion cannot grow owing to the small bubble oscillations at small pressure amplitudes and bubbles in the case of molten metals. The descending order of the bubble oscillation amplitude, as shown in this figure, was the same as that of the bubble instability.Fig. 10 Time variation of bubble radius with different liquids at 20 kHz with the pressure amplitude of 1.0 × 105 Pa, when the bubble radius is 10 μm.

The bubble size and stability were compared with the experimental results [43]. According to the in-situ observation of the Al-Cu alloy system at 30 kHz ultrasound, the bubble modal radius in the aluminum alloy melt was ∼10–20 μm, which was much larger than that in water. The measured bubble modal radius corresponded to the simulated radius when the pressure amplitude was ∼1.0 × 105 Pa. From these results, we can conclude that the stability analysis predicted the stable bubble diameter. The measured bubble size was slightly larger than those in other studies [46], [47]. For example, the most frequent radii are ∼30–50 μm [46] and ∼25 μm [47] in Al-Cu alloy melts. As shown in the stability diagrams, the bubble stability was largely affected by the pressure amplitude, and inactive bubbles were counted in the experiments. Therefore, the variation in the measured values was not significant in this range. No experimentally measured values are available for magnesium alloy melts. However, the bubble size in a magnesium melt can be predicted to be almost the same or slightly smaller than that in an aluminum melt. This information is helpful for understanding the phenomena occurring in ultrasonic casting of magnesium melts.

Previously, the radial dynamics of cavitation bubbles have been solved in aluminum alloy melts [51], [52], [53], [54], [55], [56], [57], [88] and magnesium alloy melts [89], [90]. To investigate ultrasonic processing in aluminum alloy melts, many researchers set the equilibrium bubble radius to 5 μm [51], [52], [53], [54], [55], [56], [57] and to 4.5 μm in magnesium ultrasonic casting [89], [90]. These values were similar to the equilibrium bubble size of water. However, as found in this study, the equilibrium bubble size was much larger than these values. The aforementioned studies used a smaller equilibrium bubble size. So, the influence of the bubble size on the phenomena occurring in ultrasonic casting should be carefully investigated in the future study.

4 Conclusion

In this study, a linear stability analysis was conducted to investigate the bubble shape stability and stable bubble size in molten aluminum and magnesium, which are expected to be used in ultrasonic casting applications. The main findings of this study are summarized as follows.• Acoustic cavitation bubbles in aluminum or magnesium melts are more stable than those in water at the same pressure amplitude. This is because the larger surface tension in the molten metal suppresses bubble oscillation, leading to smaller distortion growth.

• The acoustic cavitation bubble in the aluminum melt is slightly more stable than that in the magnesium melt. This is because the surface tension of the aluminum melt is larger than that of magnesium.

• Afterbounce instability does not occur in aluminum and magnesium melts, although it can occur in water at 20 kHz. This is because the number of small bubble oscillations after the first collapse is smaller in the molten metal than in water.

• An acoustic cavitation bubble in an aluminum or magnesium melt is destabilized with an increase in the ultrasonic frequency. This is because the moving velocity and acceleration of the bubble wall increase with an increase in the ultrasonic frequency.

• Stable zones in the stability diagrams of aluminum and magnesium melts exist even with a larger bubble radius. This is because the periods of bubble oscillations and ultrasonic oscillations do not match, leading to less distortion growth.

In a previous study on the ultrasonic casting of aluminum using numerical simulations, the equilibrium radius was set to be much smaller than that derived in this study. Previous studies used the equilibrium radius of water instead of that of a molten metal. In future studies, researchers should reconsider the influence of the equilibrium bubble radius on the numerical results of ultrasonic casting of molten metals.

CRediT authorship contribution statement

Takuya Yamamoto: Writing – review & editing, Writing – original draft, Visualization, Validation, Software, Methodology, Investigation, Funding acquisition, Formal analysis, Data curation, Conceptualization.

Declaration of competing interest

The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Takuya Yamamoto reports financial support was provided by Japan Science and Technology Agency. If there are other authors, they declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgements

This research was supported in part by JST 10.13039/501100009023 PRESTO (Grant Number JPMJPR22OA ), Japan.
==== Refs
References

1 Yasui K. Acoustic Cavitation and Bubble Dynamics 2018 Springer
2 Lauterborn W. Kurz T. Physics of bubble oscillations Rep. Prog. Phys. 73 2010 106501
3 Philipp A. Lauterborn W. Cavitation erosion by single laser-produced bubbles J. Fluid Mech. 361 1998 75 116
4 Zhang Z.Y. Zhang H.S. Surface tension effects on the behavior of a cavity growing, collapsing, and rebounding near a rigid wall Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 70 2004 056310
5 Mifsud J. Lockerby D. Chung Y. Gordon J. Numerical simulation of a confined cavitating gas bubble driven by ultrasound Phys. Fluids 33 2021 122114
6 Huang J. Wang J. Huang J. Lv P. Li H. Wang Y. Effects of wall wettability on vortex flows induced by collapses of cavitation bubbles: a numerical study Phys. Fluids 35 2023 087122
7 Yamamoto T. Komarov S.V. Dynamic behavior of acoustic cavitation bubble originated from heterogeneous nucleation J. Appl. Phys. 128 2020 044702
8 Suslick K.S. Sonochemistry Science 247 1990 1439 1445 17791211
9 Suslick K.S. Didenko Y. Fang M.M. Hyeon T. Kolbeck K.J. McNamara W.B. III Mdleleni M.M. Wong M. Acoustic cavitation and its chemical consequences Philos. Trans. Math. Phys. Eng. Sci. 357 1999 335 353
10 Canselier J.P. Delmas H. Wilhelm A.M. Abismaïl B. Ultrasound emulsification—An overview J. Dispers. Sci. Technol. 23 2002 333 349
11 Modarres-Gheisari S.M.M. Gavagsaz-Ghoachani R. Malaki M. Safarpour P. Zandi M. Ultrasonic nano-emulsification - A review Ultrason. Sonochem. 52 2019 88 105 30482437
12 Raman K.A. Rosselló J.M. Ohl C.-D. Cavitation induced oil-in-water emulsification pathways using a single laser-induced bubble Appl. Phys. Lett. 121 2022 194103
13 Stepisnik Perdih T. Zupanc M. Dular M. Revision of the mechanisms behind oil-water (O/W) emulsion preparation by ultrasound and cavitation Ultrason. Sonochem. 51 2019 298 304 30327174
14 Yamamoto T. Matsutaka R. Komarov S.V. High-speed imaging of ultrasonic emulsification using a water-gallium system Ultrason. Sonochem. 71 2021 105387
15 Yamamoto T. Komarov S.V. Liquid jet directionality and droplet behavior during emulsification of two liquids due to acoustic cavitation Ultrason. Sonochem. 62 2020 104874
16 Ren Z. Han H. Zeng H. Sun C. Tagawa Y. Zuo Z. Liu S. Interactions of a collapsing laser-induced cavitation bubble with a hemispherical droplet attached to a rigid boundary J. Fluid Mech. 976 2023 A11
17 Komarov S. Yamamoto T. Sun J. Fabrication of Al-Bi frozen emulsion alloys due to high-intense ultrasound irradiation J. Alloy. Compd. 859 2021 158231
18 Zhang Y. Yuan S. Wang L. Investigation of capillary wave, cavitation and droplet diameter distribution during ultrasonic atomization Exp. Therm. Fluid Sci. 120 2021 110219
19 Zhang Y. Yuan S. Gao Y. Spatial distribution and transient evolution of sub-droplet velocity and size in ultrasonic atomization Exp. Therm. Fluid Sci. 140 2023 110761
20 Yuan S. Zhang Y. Gao Y. Faraday wave instability characteristics of a single droplet in ultrasonic atomization and the sub-droplet generation mechanism Exp. Therm. Fluid Sci. 134 2022 110618
21 Kudo T. Sekiguchi K. Sankoda K. Namiki N. Nii S. Effect of ultrasonic frequency on size distributions of nanosized mist generated by ultrasonic atomization Ultrason. Sonochem. 37 2017 16 22 28427620
22 Monti C. Turani M. Papis K. Bambach M. A new Al-Cu alloy for LPBF developed via ultrasonic atomization Mater. Des. 229 2023 111907
23 Yamanaka T. Hayashi Y. Takizawa H. Near-room temperature synthesis of Zn<sup>2+</sup>-doped γ-Ga<sub>2</sub>O<sub>3</sub> nanoparticles via direct oxidation of Zn–Ga alloy by ultrasound J. Ceram. Soc. Jpn. 131 2023 100 105
24 Fuentes-Garcia J.A. Santoyo-Salzar J. Rangel-Cortes E. Goya G.F. Cardozo-Mata V. Pescador-Rojas J.A. Effect of ultrasonic irradiation power on sonochemical synthesis of gold nanoparticles Ultrason. Sonochem. 70 2021 105274
25 Zonarsaghar A. Mousavi-Kamazani M. Zinatloo-Ajabshir S. Sonochemical synthesis of CeVO4 nanoparticles for electrochemical hydrogen storage Int. J. Hydrogen Energy 47 2022 5403 5417
26 Ilić N. Andalib A. Lippert T. Knoop O. Franke M. Bräutigam P. Drewes J.E. Hübner U. Ultrasonic degradation of GenX (HFPO-DA) – Performance comparison to PFOA and PFOS at high frequencies Chem. Eng. J. 472 2023 144630
27 Kewalramani J.A. Bezerra de Souza B. Marsh R.W. Meegoda J.N. Contributions of reactor geometry and ultrasound frequency on the efficiency of sonochemical reactor Ultrason. Sonochem. 98 2023 106529
28 Komarov S. Yamamoto T. Fang Y. Hariu D. Combined effect of acoustic cavitation and pulsed discharge plasma on wastewater treatment efficiency in a circulating reactor: a case study of Rhodamine B Ultrason. Sonochem. 68 2020 105236
29 Eskin D.G. Tzanakis I. Wang F. Lebon G.S.B. Subroto T. Pericleous K. Mi J. Fundamental studies of ultrasonic melt processing Ultrason. Sonochem. 52 2019 455 467 30594518
30 Eskin G.I. Influence of cavitation treatment of melts on the processes of nucleation and growth of crystals during solidification of ingots and castings from light alloys Ultrason. Sonochem. 1 1994 S59 S63
31 Eskin G.I. Broad prospects for commercial application of the ultrasonic (cavitation) melt treatment of light alloys Ultrason. Sonochem. 8 2001 319 325 11441617
32 Eskin G.I. Ultrasonic Treatment of Light Alloy Metals 1998 CRC Press London
33 Komarov S. Yamamoto T. Development and application of large-sized sonotrode systems for ultrasonic treatment of molten aluminum alloys Miner. Met. Mater. Ser. 2019 2019 1597 1604
34 Pragathi P. Elansezhian R. Mechanical and microstructure behaviour of aluminum nanocomposite fabricated by squeeze casting and ultrasonic aided squeeze casting: a comparative study J. Alloy. Compd. 956 2023
35 Yusof N.S.M. Anandan S. Sivashanmugam P. Flores E.M.M. Ashokkumar M. A correlation between cavitation bubble temperature, sonoluminescence and interfacial chemistry - A minireview Ultrason. Sonochem. 85 2022 105988
36 Jaime R.F. Puga H. Prokic M. Söderhjelm C. Apelian D. Fundamentals of ultrasonic treatment of aluminum alloys Int. J. Met. 2024
37 Balasubramani N. Venezuela J. StJohn D. Wang G. Dargusch M. A review of the origin of equiaxed grains during solidification under mechanical stirring, vibration, electromagnetic, electric-current, and ultrasonic treatments J. Mater. Sci. Technol. 144 2023 243 265
38 Tzanakis I. Lebon G.S.B. Eskin D.G. Pericleous K.A. Characterisation of the ultrasonic acoustic spectrum and pressure field in aluminium melt with an advanced cavitometer J. Mater. Process. Technol. 229 2016 582 586
39 Tzanakis I. Lebon G.S.B. Eskin D.G. Pericleous K. Investigation of the factors influencing cavitation intensity during the ultrasonic treatment of molten aluminium Mater. Des. 90 2016 979 983
40 Tzanakis I. Lebon G.S. Eskin D.G. Pericleous K.A. Characterizing the cavitation development and acoustic spectrum in various liquids Ultrason. Sonochem. 34 2017 651 662 27773292
41 Yamamoto T. Kubo K. Komarov S.V. Characterization of acoustic streaming in water and aluminum melt during ultrasonic irradiation Ultrason. Sonochem. 71 2021 105381
42 Tzanakis I. Xu W.W. Eskin D.G. Lee P.D. Kotsovinos N. In situ observation and analysis of ultrasonic capillary effect in molten aluminium Ultrason. Sonochem. 27 2015 72 80 26186822
43 Xu W.W. Tzanakis I. Srirangam P. Mirihanage W.U. Eskin D.G. Bodey A.J. Lee P.D. Synchrotron quantification of ultrasound cavitation and bubble dynamics in Al-10Cu melts Ultrason. Sonochem. 31 2016 355 361 26964960
44 Wang F. Tzanakis I. Eskin D. Mi J. Connolley T. In situ observation of ultrasonic cavitation-induced fragmentation of the primary crystals formed in Al alloys Ultrason. Sonochem. 39 2017 66 76 28732991
45 Tzanakis I. Xu W.W. Lebon G.S.B. Eskin D.G. Pericleous K. Lee P.D. In situ synchrotron radiography and spectrum analysis of transient cavitation bubbles in molten aluminium alloy Phys. Procedia 70 2015 841 845
46 Huang H. Shu D. Fu Y. Wang J. Sun B. Synchrotron radiation X-ray imaging of cavitation bubbles in Al-Cu alloy melt Ultrason. Sonochem. 21 2014 1275 1278 24433976
47 Mirihanage W. Xu W. Tamayo-Ariztondo J. Eskin D. Garcia-Fernandez M. Srirangam P. Lee P. Synchrotron radiographic studies of ultrasonic melt processing of metal matrix nano composites Mater. Lett. 164 2016 484 487
48 Leighton T.G. The Acoustic Bubble 1994 Academic Press Ashford
49 Eller A. Flynn H.G. Rectified diffusion during nonlinear pulsations of cavitation bubbles J. Acoust. Soc. Am. 37 1965 493 503
50 Wang B. Tan D. Lee T.L. Khong J.C. Wang F. Eskin D. Connolley T. Fezzaa K. Mi J. Ultrafast synchrotron X-ray imaging studies of microstructure fragmentation in solidification under ultrasound Acta Mater. 144 2018 505 515
51 Lebon B. Tzanakis I. Pericleous K. Eskin D. Numerical modelling of the ultrasonic treatment of aluminium melts: an overview of recent advances Materials (Basel) 12 2019
52 Lebon G.S.B. Tzanakis I. Djambazov G. Pericleous K. Eskin D.G. Numerical modelling of ultrasonic waves in a bubbly Newtonian liquid using a high-order acoustic cavitation model Ultrason. Sonochem. 37 2017 660 668 28427680
53 Lebon G.S.B. Tzanakis I. Pericleous K. Eskin D. Experimental and numerical investigation of acoustic pressures in different liquids Ultrason. Sonochem. 42 2018 411 421 29429686
54 Lebon G.S.B. Salloum-Abou-Jaoude G. Eskin D. Tzanakis I. Pericleous K. Jarry P. Numerical modelling of acoustic streaming during the ultrasonic melt treatment of direct-chill (DC) casting Ultrason. Sonochem. 54 2019 171 182 30755390
55 Subroto T. Lebon G.S.B. Eskin D.G. Skalicky I. Roberts D. Tzanakis I. Pericleous K. Numerical modelling and experimental validation of the effect of ultrasonic melt treatment in a direct-chill cast AA6008 alloy billet J. Mater. Res. Technol. 12 2021 1582 1596
56 Tonry C.E.H. Bojarevics V. Djambazov G. Pericleous K. Contactless ultrasonic treatment in direct chill casting JOM 72 2020 4082 4091
57 Beckwith C. Djambazov G. Pericleous K. Subroto T. Eskin D.G. Roberts D. Skalicky I. Tzanakis I. Multiphysics modelling of ultrasonic melt treatment in the hot-top and launder during direct-chill casting: path to indirect microstructure simulation Metals 11 2021 674
58 Chen Y.-J. Hsu W.-N. Shih J.-R. The effect of ultrasonic treatment on microstructural and mechanical properties of cast magnesium alloys Mater. Trans. 50 2009 401 408
59 Khosro Aghayani M. Niroumand B. Effects of ultrasonic treatment on microstructure and tensile strength of AZ91 magnesium alloy J. Alloy. Compd. 509 2011 114 122
60 Emadi P. Ravindran C. The influence of high temperature ultrasonic processing time on the microstructure and mechanical properties AZ91E magnesium alloy J. Mater. Eng. Perform. 30 2021 1188 1199
61 Burdin F. Tsochatzidis N.A. Guiraud P. Wilhelm A.M. Delmas H. Characterisation of the acoustic cavitation cloud by two laser techniques Ultrason. Sonochem. 6 1999 43 51 11233937
62 Tsochatzidis N.A. Guiraud P. Wilhelm A.M. Delmas H. Determination of velocity, size and concentration of ultrasonic cavitation bubbles by the phase-Doppler technique Chem. Eng. Sci. 56 2001 1831 1840
63 Brotchie A. Grieser F. Ashokkumar M. Effect of power and frequency on bubble-size distributions in acoustic cavitation Phys. Rev. Lett. 102 2009 084302
64 Iida Y. Ashokkumar M. Tuziuti T. Kozuka T. Yasui K. Towata A. Lee J. Bubble population phenomena in sonochemical reactor: I estimation of bubble size distribution and its number density with pulsed sonication - laser diffraction method Ultrason. Sonochem. 17 2010 473 479 19811943
65 Plesset M.S. On the stability of fluid flows with spherical symmetry J. Appl. Phys. 25 1954 96 98
66 Eller A.I. Crum L.A. Instability of the motion of a pulsating bubble in a sound field J. Acoust. Soc. Am. 47 1970 762 767
67 Hsieh D.Y. On thresholds for surface waves and subharmonics of an oscillating bubble J. Acoust. Soc. Am. 56 1974 392 393
68 Prosperetti A. Bubble phenomena in sound fields: part two Ultrasonics 22 1984 115 124
69 Mei C.C. Zhou X. Parametric resonance of a spherical bubble J. Fluid Mech. 229 1991 29 50
70 Brenner M.P. Lohse D. Dupont T.F. Bubble shape oscillations and the onset of sonoluminescence Phys. Rev. Lett. 75 1995 954 957 10060160
71 Hilgenfeldt S. Lohse D. Brenner M.P. Phase diagrams for sonoluminescing bubbles Phys. Fluids 8 1996 2808 2826
72 Brenner M.P. Hilgenfeldt S. Lohse D. Single-bubble sonoluminescence Rev. Mod. Phys. 74 2002 425 484
73 Puente G.F. Garcia-Martinez P. Bonetto F.J. Single-bubble sonoluminescence in sulfuric acid and water: bubble dynamics, stability, and continuous spectra Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 75 2007 016314
74 Hao Y. Prosperetti A. The effect of viscosity on the spherical stability of oscillating gas bubbles Phys. Fluids 11 1999 1309 1317
75 Holzfuss J. Surface-wave instabilities, period doubling, and an approximate universal boundary of bubble stability at the upper threshold of sonoluminescence Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 77 2008 066309
76 Godinez F.A. Navarrete M. Influence of liquid density on the parametric shape instability of sonoluminescence bubbles in water and sulfuric acid Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 84 2011 016312
77 Leong T. Yasui K. Kato K. Harvie D. Ashokkumar M. Kentish S. Effect of surfactants on single bubble sonoluminescence behavior and bubble surface stability Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 89 2014 043007
78 Klapcsik K. Hegedus F. Study of non-spherical bubble oscillations under acoustic irradiation in viscous liquid Ultrason. Sonochem. 54 2019 256 273 30718178
79 Yamamoto T. Linear stability analysis for bubble shape of acoustic cavitation with different ultrasonic frequency Phys. Fluids 36 2024 in Press
80 Keller J. Miksis M. Bubble oscillations of large amplitude J. Acoust. Soc. Am. 68 1980 628 633
81 Yasui K. Effect of liquid temperature on sonoluminescence Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 64 2001 016310
82 Yasui K. Variation of liquid temperature at bubble wall near the sonoluminescence threshold J. Phys. Soc. Jpn. 65 1996 2830 2840
83 Storey B.D. Szeri A.J. Mixture segregation within sonoluminescence bubbles J. Fluid Mech. 396 1999 203 221
84 Leitner M. Leitner T. Schmon A. Aziz K. Pottlacher G. Thermophysical properties of liquid aluminum Metall. Mater. Trans. A 48 2017 3036 3045
85 Cao X. Jahazi M. Immarigeon J.P. Wallace W. A review of laser welding techniques for magnesium alloys J. Mater. Process. Technol. 171 2006 188 204
86 Yamamoto T. Effect of ultrasonic frequency on mass transfer of acoustic cavitation bubble Chem. Eng. Sci. 300 2024 120654
87 Holt R.G. Gaitan D.F. Observation of stability boundaries in the parameter space of single bubble sonoluminescence Phys. Rev. Lett. 77 1996 3791 3794 10062309
88 Sun J. Yamamoto T. Komarov S. Behavior of Al-Zr intermetallic compound particles under high-amplitude ultrasound irradiation into molten aluminum J. Alloy. Compd. 912 2022 165128
89 Shao Z.-W. Le Q.-C. Zhang Z.-Q. Cui J.-Z. Numerical simulation of acoustic pressure field for ultrasonic grain refinement of AZ80 magnesium alloy Trans. Nonferrous Met. Soc. Chin. 21 2011 2476 2483
90 Song S. Zhou X. Li L. Ma W. Numerical simulation and experimental validation of SiC nanoparticle distribution in magnesium melts during ultrasonic cavitation based processing of magnesium matrix nanocomposites Ultrason. Sonochem. 24 2015 43 54 25559849
