Simulation of Ultrasonic Induced Cavitation and Acoustic Streaming in Liquid and Solidifying Aluminum

: Ultrasonic treatment (UST), more precisely, cavitation and acoustic streaming, of liquid light metal alloys is a very promising technology for achieving grain and structure reﬁnement, and therefore, better mechanical properties. The possibility of predicting these process phenomena is an important requirement for understanding, implementing, and scaling this technology in the foundry industry. Using an established (casting) computational ﬂuid dynamics (CFD)-simulation tool, we studied the ability of this software to calculate the onset and expansion of cavitation and acoustic streaming for the aluminum alloy A356, partly depending on different radiator geometries. A key aspect was a holistic approach toward pressure distribution, cavitation, and acoustic streaming prediction, and the possibility of two- and (more importantly) three-dimensional result outputs. Our feasibility analysis showed that the simulation tool is able to predict the mentioned effects and that the results obtained are in good agreement with the results and descriptions of previous investigations. Finally, capabilities and limitations as well as future challenges for further developments are discussed.

Another important phenomenon is the presence of acoustic streaming and convective flows, which support and increase mass and heat transfer in the bulk of the liquid, and influence the solidification mechanisms of the melt. Thus, acoustic streaming distributes by shock waves newly created Al nucleation sites or other particles throughout the bulk of the melt and ensures a more homogeneous distribution of contents [1,7,9,11,[28][29][30]. Together, cavitation and acoustic streaming have high potential for the effective isothermal and solidification-accompanying treatment of an alloy and allow for the improvement of the grain structure and final mechanical properties of casting products [1,7,[31][32][33].
So far, one obstacle for the comprehensive entrance of this technology into the foundry industry has been the missing potential for simulating the processing and influence of ultrasonic melt treatments [10,34]. In recent years, many researchers have investigated acoustic streaming using (for instance) particle image velocimetry (PIV) analyses to predict the propagation of acoustic waves, as well as the development of acoustic streaming and its influence on the fluid using simulation tools [8,17,21,[35][36][37][38][39][40]. An overview of (some) research publications on the subject of simulation of ultrasonic treatment so far can be taken from [41] and Table 1.
However, to the best of our knowledge, almost none of these solutions offers a comprehensive approach. Some simulations concentrate mainly on acoustic streaming [36,38], but simplify the flow behaviour by neglecting the turbulences or treat the radiator as fluid inlet to evoke acoustic streaming and neglect cavitation to reduce the necessary computational power [21,35,42]. Others focus on the calculation of the occurring acoustic pressure distribution [8,17], sometimes in combination with acoustic streaming [43], or either solely on cavitation [44][45][46]. Nearly all simulation approaches are based on two-dimensional modeling. Therefore, the aim of this investigation was to create a three-dimensional simulation model with the possibility of two-dimensional sectional views using the commercial computational fluid dynamics (CFD) code (CCC) simulation tool FLOW-3D for aluminum alloy A356. This way, it could be possible to predict the beginning and development of cavitation and acoustic streaming, allow for overlap with simulations of the solidification process (as in [47]), or enable preinvestigations for process design. Furthermore, Lebon et al. [41] noted rightly that the modeling behind the results in [48] was not further described. Therefore, a further goal is the more detailed description of the modeling and results behind this work. The overriding aim is to move towards the industrial commercialization of this technology in the foundry industry. Table 1. List of publications on the simulation of ultrasonic treatment, listed according to year of publication and calculated ultrasonic effects (AP/WP-Acoustic pressure/wave propagation, C-Cavitation, AS-Acoustic Streaming, S-Solidification). Note that most of the studies listed below are using the respective software for incorporation of individual ultrasonic-related codes.

General
The simulation model was created and run using the CCC-simulation software FLOW-3D v11.2, while the results were analyzed with FlowSight v11.2. Both software tools were developed by Flow Science, Inc. (Santa Fe, NM, USA). To enable a straightforward application and adaption of the presented model, the following model description follows the workflow within FLOW-3D.
The simulated time for wave propagation was set to 0.0001 s, which conforms to two full radiator oscillations. Cavitation and acoustic streaming were investigated for the first 0.01 s and 5 s of UST, respectively, to investigate both developments over a proper period of time. Differing times will be mentioned.

Fluid
Within the scope of the model, UST was simulated for aluminum silicon alloy A356 at 973.15 K because this alloy has significant industrial applications. The main parameters are presented in Table 2. In the case of FLOW-3D, a linear system of equations to be solved is derived from the continuity equation, which maintains mass conservation of fluid. The general mass continuity equation is with V F as the fractional volume open to flow, ρ as the fluid density, R DIF as a turbulent diffusion term, R SOR as a mass source, and R and ξ as mesh-dependent coefficients. A x is the fractional area open to flow in the x-direction, A y and A z are similar area fractions for flow in the yand z-directions, respectively. The computation of acoustic pressure waves in fluids that would otherwise be treated as incompressible is done by an approximation of the density time derivative: ∂ρ ∂t where ρ and c are the respective fluid density and adiabatic speed of sound whereas p corresponds to the pressure. Integration of sound velocity is realized by calculation and inclusion of the compressibility factor β according to From Equation (2), a modified continuity equation results as follows: The complete description of the model can be taken from [70,71].

Physics
The appropriate choice of physical models directly influences the quality of the simulation model. Besides gravitation and the activation of the general moving object (GMO) model, for the later movement of the radiator, the bubble and phase change model, cavitation model, surface tension model, and viscosity and turbulence model are of primary importance. The respective specific parameters and values used are listed in Table 3. The influence of the resulting pressure changes is calculated by Equations (2) and (4). Cavitation is solved by a transport equation which calculates advection, production, and dissipation of cavitation volume fraction, according to the following equations [70]: where V cav is the computed cavitation volume fraction, C p is the cavitation production coefficient, C e the evaporation coefficient, C c is the condensation coefficient, E turb is the turbulent kinetic energy (alternatively, 10% of the total kinetic energy if no turbulence model selected), σ is the surface tension, P cav is the specified cavitation pressure, P is the local fluid pressure, and f cav is the mass fraction of cavitation in the cell with ρ l and ρ v as densities of liquid and vapor, respectively. The coefficients used for cavitation production and dissipation are default values (0.02 and 0.01, respectively) that can be adjusted. The gas species inside the cavitation bubbles was assumed to be hydrogen gas at 973.15 K. Gas density was calculated by the general gas equation where p is the gas pressure, V is the gas volume, n is the amount of substance of gas, R m is the universal gas constant, and T is the thermodynamic temperature of the gas [72]. Using the data in Table 4 with the molar mass of H 2 , the gas density of cavitation bubbles amounts to 0.025 kg/m 3 . Regarding the bubbles, the relationship between pressure, volume, and temperature will follow an adiabatic law. Pressure p in each void region is evaluated on the void region volume V using the isentropic model of expansion or compression in which pV γ is constant, where γ is the isentropic exponent. Break up, coalesce, receive, and loss of mass and energy at the mesh boundaries are approximated by the bubble pressure model. The pressure in each bubble will be inversely proportional to its volume to the power of the isentropic exponent γ. If a bubble neither breaks up nor coalesces, the pV γ relation determines the new pressure. For an accurate solution, the bubbles must be resolved by a minimum of three cells across the diameter [70]. Since this criterion is not met to ensure acceptable computing times, bubble opening and pulsation are neglected. For the same reason, the effect of nonequilibrium evaporation and condensation at the bubble wall was also not taken into account. This simplification is appropriate as, in the case of aluminum and hydrogen, evaporation and condensation of gas within the bubbles are negligible at the involved pressures and temperatures. Hence, the model still allows taking the surface tension and density of cavitation bubbles into account. Regarding the turbulence model, a laminar flow model was chosen, because no differences were observed in comparisons of different models (e.g., k − or RNG) with the chosen cell size. The streaming is provoked by the radiator movement and the resulting pressure changes that were calculated by the GMO-model and controlled by a sinusoidal translational velocity component in the z-direction, corresponding to an ultrasonic system with 20 kHz and a peak-to-peak amplitude of 35 µm. The sinusoidal translational velocity is calculated by with where v(t) corresponds to the angular frequency, ω is the angular frequency with f as the frequency, s 0 is the amplitude (17.5 µm), and φ 0 is the phase angle. The GMO-model calculates the kinetic energy, which is caused by the radiator movement, and its transfer on the fluid by momentum and mass conservation. For moving object/fluid coupling, the explicit numerical method was chosen, in which the fluid and GMO motions-i.e., the radiator-of each time step are calculated using force and velocity data from the previous time step.
To investigate the influence of acoustic streaming on distribution of potential oxide particles and dendrite fragments, additionally, a particle model was activated. Two kinds of particles were defined: Al 2 O 3 with a particle size of 0.1 µm and a density of 3950 kg/m 3 ; as well as pure Al with a particle size of 20 µm and a density of 2700 kg/m 3 . Such considerations with regard to particle distribution could, for example, be helpful in the simulation-based investigation of UST-supported production of aluminum matrix composites (AMC).

Geometry and Movement
To reduce the required calculation time, only essential geometric elements were simulated, without neglecting the basic conditions of a real system. Therefore, the vessel wall was defined as a boundary without heat transfer, to avoid additional cells for the mesh and calculation of heat transfer. At the mesh boundaries, full reflection of sound waves/without any losses was calculated. The fluid volume was 100 × 100 × 100 mm (xyz), for a total volume of 1000 cm 3 (less volume of immersed radiator) with temperature for aluminum of 973.15 K. The modeled radiator consists of a ceramic material that has a length of 60 mm, a diameter of 22 mm and, unless otherwise mentioned, the same temperature as the fluid. The radiator was positioned in the middle of the vessel at an immersion depth of 30 mm. The atmospheric pressure and temperature were set to 101,325 Pa and 293.15 K, respectively. The particle source was placed directly below the radiator tip and had a size of 10 × 10 × 2 mm. A total of 800 particles per second (400 Al 2 O 3 and 400 Al) were produced with a uniform generation rate. Fluid motion affects the particles, but not vice versa. To obtain adequate simulation results, a minimum cell size of 1 mm was necessary. The overall three-dimensional system, discretized with 1,100,000 cells, can be seen in Figure 1. The simulation setup for the solidification-accompanying UST can be taken from Figure 2.
To achieve a precise simulation, the cooling curve was calculated in advance using the ThermoCalc software and was afterwards implemented in FLOW-3D. The simulation started at a melt temperature of 890 K (liquidus temperature). The mold wall was defined as a boundary with heat transfer and a starting temperature of 673 K. Similar to the isothermal model, at the mesh boundaries, full reflection of sound waves/without any losses was calculated. According to the experimental setup in [48], the radiator had a frequency of 20 kHz, a peak-to-peak amplitude of 35 µm, a temperature of 473 K and was placed concentrically directly above the mold cavity. The same simulation was performed without UST to use the results as a reference. As a result of the additional effort for the calculation of solidification and due to the symmetrical setup, only one-fourth of the model was discretized with a cell size of 2 mm to show the qualitative influence and simulation potential.

Time Step Definition
In view of the calculation time steps, it is important to also take the frequency into account. For calculation of the correct sinusoidal movement of the radiator, the defined time steps have to hit characteristic, positive and negative equally distributed spots on the curve, depending on the frequency of the radiator, but at least the peak points of the curves. The reduction of the time step is largely unproblematic. The smaller the time step, the higher the accuracy of simulation but the calculation periods of the model as well, whereas the maximum time step is limited by frequency, especially in the scope of larger simulation durations. For this reason, with the chosen frequency (20 kHz), a minimum of 80,000 calculation steps per second is necessary to adequately calculate the radiator's movement. Table 5 contains an overview of the chosen time steps for the following simulations.  Figure 3 shows pressure initiation and distribution through aluminum for three different stages as well as the related min/max pressures. At t = 3 µs (Figure 3a), negative pressure increases by several MPa immediately below the radiator face, then propagates in the -z-direction. Subsequently, the positive pressure amplitude develops and starts propagating in -z-direction (Figure 3b). Due to the spherical radiator geometry, propagation in the xy-direction is also clearly discernible. The highest values arise in the center of the radiator faces with 4.5 MPa and −4.2 MPa for the spherical tip. From the tip, the pressure intensity decreases with increasing distance from the radiator and shows exponential intensity loss. The measured simulated half-wavelength is about 11.5 cm. The propagating waves are reflected by the sidewalls, leading to a slightly slower pressure decrease in these areas because of interfering waves (depending on the simulation boundary conditions) (Figure 3c).

Cavitation Development
The results for cavitation development for the spherical tip and additionally, for reasons of comparison, for a flat radiator tip, are demonstrated in Figure 4. Furthermore, the radiator-depending collapse activity is demonstrated, since FLOW-3D allows tracking of collapsed bubbles by void particles. Cavitation mainly develops in the area immediately below the radiator, where the highest pressure peaks occur and, consequently, the best conditions for cavitation exist. Note that the cavitation gas volume is bigger in the case of the spherical radiator, but even if the cavitation zone is smaller in the case of the flat radiator, collapse activity is much higher compared to the spherical radiator. Figure 5a demonstrates that in the case of the spherical tip, the total cavitation gas volume fraction at the end of the calculated time period is more than twice as high compared to the flat radiator. However, Figure 5b reveals the much higher collapsing activity in the case of the flat radiator.   The comparison of simulation results with and without the activated cavitation model (apart from that identical) reveals an influence of the cavitation area on wave propagation, as demonstrated in Figure 6 for t = 0.005 s. Within the cavitation zone, a uniform pressure ambient seems to exist (Figure 6a). Without the activated cavitation model, such a zone is not recognizable (Figure 6b). This is supported by Figure 7, where the cavitation activity ties to the current pressure conditions, for the first 7-mm is constant. The pressure waves (and thus, the pressure peaks) are damped by the cavitation cloud, so for the current example, the negative pressure peak differs by a factor of almost four.   Figure 8 shows the development, pattern, and propagation of acoustic streaming. First, the streaming front propagates in the -z-direction, with a streaming center velocity of 1.5 to 1.7 m/s. The streaming reaches the bottom after approximately 0.1 s and decelerates because of redirection in the xy-plane and loses the majority of its velocity. Additionally, a feedback takes place which leads to a streaming deceleration from the radiator to the bottom. The streaming velocity in this area henceforth is 0.6 to 0.7 m/s. The continuing streaming forces the horizontal circular propagation of the fluid, and leads to perpendicular ascension on the sidewalls. After reaching the walls, the stream is redirected once more. It decelerates once again and climbs up the vessel walls. The velocities at the bottom and walls are approximately 0.2 to 0.4 m/s and 0.05 to 0.2 m/s, respectively. Reaching a specific height, the sidewall-stream falls perpendicular to the wall in the middle of the vessel near the mainstream and completes a circular flow, which enables a continuing mixing of the fluid. After 5 s, the whole fluid volume is in motion. Figure 9 illustrates the corresponding centerline-velocities in -z-direction for 1, 2, 3, 4, and 5 s. Figure 10 shows the acoustic streaming velocity pattern for four different stages of acoustic streaming development and reveals the formation of eddies along the main stream (radiator to the bottom). At t = 3.950 s, the whole fluid is in motion and turbulences in the fluid are recognizable.

Cavitation and Mass Transport
The influence of acoustic streaming on cavitation is shown in Figure 11. After 0.05 s, the cavitation zone around the tip reaches a steady-state and gets penetrated from the downward-directed streaming front. Cavitation bubbles are entrained by acoustic streaming and transported in the bulk liquid. As a result, the size of the zone decreases considerably. It can be seen that the density of cavitation bubbles drastically decreases with growing distance from the radiation surface (e.g., t = 0.2 s). After approximately 40 mm, the majority of the created cavitation bubbles disappeared because of the decreasing pressure values in the growing distance to the radiator. A certain bubble proportion seems to be stable and gets transported and distributed through the bulk of the fluid volume. The remaining bubbles follow the flow of acoustic streaming on the bottom and side walls (t = 1.5 s) and are distributed throughout the whole fluid volume (t = 5 s).  Figure 12 shows the influence of acoustic streaming on the distribution of particles within the fluid. As described above, two kinds of particles were defined: Al 2 O 3 and pure Al. In Figure 13a, the Al 2 O 3 -and Al-particles are the dark and bright dots, respectively. In comparison with Figure 11, it is recognizable that all particles, independent of their size and density within the model, are transported by the stream and follow its direction. In a matter of 5 s, the particles are distributed in the bulk of the melt volume. In Figure 13b, the particles are colored as a function of appearance. The course of the particles, particularly the early ones, reveals the characteristic pattern of the circular flow, since some of the first particles (dark) are within the mainstream again.

Solidification and Heat Transport
The superimposition of UST (primarily acoustic streaming) on a cooling/solidification process is demonstrated in Figure 13, where the radiator is brought into the melt via the feeder, which is a possible way for solidification accompanying treatment of real castings. At t = 10 s (Figure 13a) and t = 20 s (Figure 13b), it is recognizable that acoustic streaming continuously transports warmer melt from the center area immediately under the radiator to the more rapidly cooling boundary areas. This process seems to counteract classical exogenous, directional solidification, since a widely homogeneous temperature regime exists until the final phase of the solidification (Figure 13c Although Figure 13 shows that even with UST an exogenous solidification takes place, the analysis of the temperature regime in Figure 15a shows that acoustic streaming leads to a significant reduction of temperature gradient compared to solidification without UST. In the present case, the average temperature gradient for the solidification process with and without UST is about 119.5 K and 180.3 K, respectively. As a result of the lower temperature gradient and thus the lower average temperature (Figure 15b), the cast is faster solidified compared to the cast without UST.
Furthermore, a chilling effect is recognizable resulting from the radiator stem, that is demonstrated in Figure 16. The lower temperature of the radiator compared to the melt leads to a temperature reduction of the melt of about 10 K in its immediate vicinity.

Acoustic Pressure Wave Propagation
The simulation allows us to estimate the influence of radiator and vessel geometry on the creation and propagation of pressure waves, establishing a basis for the calculation and prediction of the cavitation zone. Regarding the pressure propagation, the round tip has only one point in the center of it that moves strictly perpendicular on the fluid. For this reason, fluid can stream past around the spherical tip which makes the tip move smoother and produce lower pressure peaks compared to the usage of a flat radiator. For the same reason, pressure propagates spherically and in more volume around the tip cavitation where gas volume is produced. The pressure decrease is calculated corresponding to where α is the loss coefficient, I is the intensity, and x is the distance to radiator [1]-and thus shows correct calculation regarding the theoretical description of the pressure/sound wave propagation.
The results show high qualitative agreement with the results of other authors, e.g., [8,17].

Cavitation Development
The comparison between the different radiator geometries confirms that the radial wave propagation of the round tip seems to cause the necessary pressure conditions for cavitation in a higher fluid volume around the tip. As a consequence, cavitation occurs in a higher volume of the fluid as well compared to the flat radiator. In both cases, the cavitation cloud exhibits a cavitation transition-layer-zone which allows a clear distinction of the cavitation cloud compared to the rest of the fluid. In both cases of radiator geometry, it seems to exist a sector of almost six to seven mm of thickness with steady cavitation activity. From there, the intensity drops exponentially. The whole cavitation cloud in -z-direction is about 40 mm for both tips. This seems logical, since the lowest point of both tip geometries are at the same level, so the wave propagation proceeding from this point is the same.
An interesting aspect is the consideration of the shielding effect caused by the cavitation zone. Propagating pressure waves first must pass this area and therefore lose a large portion of the pressure wave energy to the cavitation cloud. Generally, this result coincides with the experimental investigation in [73,74] in terms of the shielding of acoustic waves by acoustic cavitation. The subsequent pressure amplitude is less than that compared to the case without cavitation zone. Figure 7 shows that the cavitation zone prevents the development and propagation of the pressure wave along the z-axis, demonstrates that the simulation takes the shielding effect into account, and considers it in the further calculation of the sound wave propagation. However, due to the coarse mesh compared to the size of the cavitation bubbles (void regions are just opened where the cavitation volume fraction is predicted to exceed 1% of a cell's volume [70]), a homogeneous distributions of cavitation gas, i.e., hydrogen, and corresponding pressure ambient within the cavitation zone was calculated, which in this form does not correspond to experimental observations. Since no speed of sound for hydrogen was/could be defined, a constant pressure and sound wave propagation occurred. Investigations on a smaller scale and with a finer mesh did confirm that higher resolution results in a more detailed calculation of the pressure distribution within the cavitation zone could be obtained.

Acoustic Streaming
The obtained results regarding the acoustic streaming velocity and streaming pattern show high-accordance with the results of [38,39,42] and especially with the PIV-analysis for a comparable ultrasonic system in [40], as demonstrated in Figure 17. In these studies, the prediction and measurements are of the same order of magnitude. Furthermore, the calculated flow pattern shows high-consensus with measurements in the PIV-study. Presumably, minor deviations derive from the different radiator geometries and the PIV measurements of a certain sectional plane requiring looking through the fluid volume, which can lead to distortions. In contrast, the evaluation of the simulation results allows a clear, uninfluenced view of the selected sectional plane. The differing centerline-velocities illustrated in Figure 9 emphasize the turbulent character of the acoustic streaming. The entire flow profile could presumably not be formed with the given modeling, but a suitable adjustment of the model is possible without problems. A great advantage is the possibility to simulate the interaction of acoustic streaming and the development of cavitation. With the distance-dependent decreasing pressure amplitude, cavitation gas volume fraction decreases as well during the transport of the cavitation bubbles from the active cavitation area through the acoustic streaming. In view of the influence of acoustic streaming on the solidification, a more homogeneous temperature distribution was recognizable, which together with cavitation finally influences the material's structure, as shown in [48]. Therefore, two streaming-related effects appear: (1) a more homogeneous solidification and structure as a result of continuous fluid mixing by acoustic streaming and (2) a faster solidification of the total cast, confirming the descriptions of [12,47]. The results emphasize the fact that acoustic streaming can influence the solidification process by enlarging the solidification front. The results concerning the cooler radiator reveal a cooling effect proceeding from the radiator. The effect increases with decreasing distance to the fluid surface. The combined chilling effect of atmosphere and radiator seems possible and, on a theoretical base, conforms to the considerations and results of [12,47].

FLOW-3D
The presented results demonstrate the ability of the investigated simulation tool to calculate and predict the most important phenomena connected with UST (i.e., cavitation, acoustic streaming, and their interaction as well as their influence on solidification). Most of the obtained results are in good agreement with the outcomes of other studies. The great advantage in the present case is the usage of a well-established CFD simulation tool and its ability to produce both two-and three-dimensional results. Additionally, it provides a foundation with which additional models and ideas can be integrated. In the current state, the tool could be used for fundamental investigations into UST or for process design, e.g., ultrasonic melt flow treatment, or isothermal and solidification-accompanying treatment. For mixing processes, the usage of the particle model allows an estimation of the time necessary for homogeneous distribution of contents, e.g., in case of AMC-production. Thus, experimental process times can be saved. Depending on cell size, the investigated time frame, and computing power, the calculation times for wave propagation and cavitation zone can take anywhere from a few minutes to several hours. Options for reducing lengthy calculation durations are two-dimensional simulations and the use of symmetric models. That way, in simple cases, only half to a fourth of a model must be calculated for qualitative results.
Although FLOW-3D provides interesting and promising results, some aspects of the method require adjustments and further developments for the specific calculation of UST. Equations (5)- (7) used for the calculation of cavitation are designed for high-velocity flows, not for the special application case of UST. A specification of the equations or the implementation of special ones for UST-which, for example, include the sound speed of gas within cavitation bubbles-might achieve a more precise result. In this vein, accurate sound wave propagation in the presence of formed cavitation zones with the corresponding impact on acoustic streaming could be calculated.
Within the scope of simulation evaluation described above, a validation of grid independency was not described in detail. The calculation time for the isothermal model with a grid resolution of 1 mm takes about five days (Intel Xeon Gold 6146 CPU 3.20 GHz, 2 × 12 cores). Investigations of two-dimensional discretizations showed that a mesh-refinement eventuates in a more detailed visual presentation of the results, while quantitative values remained mostly unchanged. For that reason, a refinement of grid resolution from 1 mm to 0.5 mm was not performed as it would result in an eight-fold increase in predicted calculation time. A coarser mesh reveals a simplified, more laminar flow calculation with largely similar velocities of acoustic streaming. The observation that the application of other turbulence models does not seem to influence the calculated flow behaviour could not be clarified so far. In this context, a correlation between calculated ultrasonic effects and mesh size is probable.

Summary
In this work, the potential for an established CFD fluid and casting simulation tool was analyzed, considering two-and three-dimensional pressure, cavitation, and acoustic streaming prediction. The most important results are as follows: • The simulation tool can predict the development and distribution of pressure waves in a fluid, induced through a modeled radiator.

•
Based on this, the development of the cavitation zone below and around the radiator tip can be calculated. Furthermore, information about cavitation intensity is provided. The results are consistent with the theoretical descriptions and demonstrated experimental behaviour of cavitation of previous investigations. With growing distance from the radiator surface, the pressure oscillations fall below the cavitation threshold, and most cavities are not able to survive in these regions. Furthermore, a shielding effect on acoustic field, proceeding from the cavitation zone, is measurable. However, for more accurate results, the used model may require adjustment for the application case of UST.

•
Rapidly changing pressure conditions lead to calculated acoustic streaming, consistent with the qualitative descriptions and circumstances leading to acoustic streaming. Furthermore, the results are in good agreement and have the same order of magnitude as the results of other studies and experimental PIV-measurements.

•
The influence of acoustic streaming on the cavitation can be calculated as well as particle transport and, for example, heat and mass transfer during solidification.

•
The simulation tool is able to create a three-dimensional prediction of cavitation and acoustic streaming. The possibility of producing two-and three-dimensional results is a great advantage for analyzing, developing, and adjusting upcoming ultrasonic systems (e.g., different radiator geometries and properties).

•
All parameters for the fluid and vessel geometries (as well as the ultrasonic system) are easily changeable and thus allow for preinvestigations into the influence of frequency, amplitude, etc., on the process and different fluids.

Conflicts of Interest:
The authors declare no conflicts of interest.

Abbreviations
The following abbreviations are used in this manuscript: AMC Aluminum matrix composites CCC Commercial CFD-Code CFD Computational fluid dynamics GMO General moving object PIV Particle image velocimetry UST Ultrasonic treatment