Numerical Investigation on Cavitation Jet in Circular Arc Curve Chamber Self-Excited Oscillation Nozzle

: In order to improve the cavitation performance of the self-excited oscillation nozzle (SEON), a novel SEON with a circular arc curve chamber was designed by changing the chamber wall proﬁle of the SEON. The performance of the circular arc curve chamber SEON was studied numerically. Taking the vapor volume distribution and the vapor volume fraction as the evaluation indexes, the inﬂuences of the chamber wall proﬁle on the cavitation performance of the circular arc curve chamber SEON were analyzed. In addition, it was compared with the broken-line chamber SEON. The numerical results show that the cavitation performance of the circular arc curve chamber SEON is ﬁrst enhanced and then weakened by increasing the circular arc radius. The circular arc curve chamber structure can form a larger central cavitation volume in the nozzle, which improves the cavitation performance of the SEON. When the circular arc radius is 2 mm, the cavitation area and the turbulent kinetic energy of the circular arc curve chamber SEON increase by 122.5% and 16.9%.


Introduction
As a ship remains in the water for a long time, a large number of aquatic organisms and algae will attach to the surface of the hull, leading to increased resistance of the ship.Cleaning the hull can not only reduce the ship's resistance and reduce greenhouse gas emissions, but also effectively prevent biological invasion and maintain the local ecological balance.Water jet cleaning is one of the common methods used in underwater hull cleaning.The self-excited oscillation cavitation jet is considered the most promising underwater hull cleaning method, due to its own pulse effect and instantaneous velocity several times higher than that of a continuous jet [1].The self-excited oscillation cavitation jet is not only outstanding in the field of underwater hull cleaning [2], but is also widely used in cutting [3], underground drilling [4,5], marine mining [6], and other industries.
Optimizing the cavitation performance of the self-excited oscillation nozzle (SEON) is the key to improving the efficiency of underwater hull cleaning.The operating conditions [7], target [8], structural parameters of the nozzle [9], and multi-field coupling [10] will affect the cavitation performance of the self-excited oscillation cavitation jet.Erosion and visualization experiments are important methods to study the characteristics of the cavitation jet [7].The erosion experiment can not only obtain the cavitation performance by calculating the mass loss before and after the erosion sample but can also obtain the surface topography of the erosion sample.The erosion process of the cavitation jet on samples is relatively complicated.Cavitation collapse affects the distribution and size of the pits on the sample surface [11,12] and is accompanied by shock wave generation [13].The visualization experiments can capture the images of the cavitation cloud, process the image data, and analyze the characteristics of the flow field.Currently, the commonly used image analysis methods include the detection of dimensionless gray intensity [14] and the proper orthogonal decomposition [15].Through the analysis of the morphology of the cavitation cloud, it can be found that the erosion ability of the cavitation jet comes from the energy released when the cavitation collapses.That is to say, the cavitation clouds with irregular boundaries or about to be broken have a better erosion effect on samples [16].It can be seen that the cavitation cloud morphology of the cavitation jet has an important influence on its erosion effect.The nozzle structure is the most important factor affecting the cavitation cloud morphology [17].Its parameters will not only affect the flow characteristics of the cavitation jet but will also influence the axial pressure of the cavitation jet [18,19].In addition, the grooves [20] and the roughness [21] on the inner surface of the nozzle will increase the energy disturbance of the jet and promote the generation of cavitation clouds.A study on the chamber structure of the SEON in the literature [22] showed that the pulse intensity of the cavitation jet was affected by changing the structural parameters of the SEON.The increase of pulse intensity in the self-excited oscillation cavitation jet is due to the transient accumulation and release of jet energy by the cavitation vortex ring in the self-excited oscillation chamber.Therefore, it is necessary to study the cavitation vortex ring in the chamber to increase the pulse intensity in the SEON.Wang et al. [23] adopted Bézier curves to optimize the chamber wall of the SEON.This structure improves the cavitation vortex ring distribution in the chamber of the SEON and enhances the cavitation performance and pulse intensity of the SEON.However, the change in the wall profile of the SEON will result in the reduction of the chamber space and the length of the collision wall.More importantly, the SEON optimized by the Bézier curve is difficult to process and apply in practice.
Therefore, in this paper, the broken line at the end of the collision wall of the chamber was designed as a circular arc curve to form a novel circular arc curve chamber SEON.The numerical investigation of the circular arc curve chamber SEON was carried out by Ansys Fluent.The cavitation performance of the SEONs with different circular arc radii was investigated.In addition, a comparative study was conducted with the broken-line chamber SEON.

Design of Broken-Line Chamber SEON
The structure of the broken-line chamber SEON is shown in Figure 1; it consists of an upper nozzle, a lower nozzle, and a self-excited oscillation chamber.The structural parameters of the self-excited oscillation chamber include the chamber length, chamber diameter, and collision wall angle.The length of the chamber refers to the distance from the upper nozzle outlet to the lower nozzle inlet, which affects the state of the jet inside the chamber and the cavitation characteristics of the outflow field.The chamber diameter is the diameter of the circumference of the chamber, which affects the position of the vortex ring in the chamber.The wall profile directly affected by the jet from the upper nozzle outlet in the chamber is the collision wall.The feedback process of the cavitation vortex ring to the upstream depends on the matching of the collision wall and the chamber wall.The angle between them is the collision-wall angle.The diameter of the lower nozzle is slightly larger than that of the upper nozzle, and their ratio has a great influence on the performance of the SEON.In this paper, the upper nozzle diameter of the broken-line chamber SEON is 1 mm.According to the range of dimensionless parameters of the SEON in Table 1 [24], the detailed structural parameters of the broken-line chamber SEON in this paper are shown in Table 2.

Design of Circular Arc Curve Chamber SEON
The high-speed jet enters the chamber from the upper nozzle and cavitation incepts in the chamber.Some of the cavitation bubbles leave the chamber from the lower nozzle, and the remainder of the cavitation bubbles return to the upstream part of the chamber through the collision wall to form the cavitation vortex ring.The cavitation vortex ring blocking the jet is the primary reason for the pulse generated by the SEON.The cavitation vortex ring can affect the pulse intensity.Figure 2 shows the accumulation and release process of the jet energy by the cavitation vortex ring in the SEON [2,23].The high-speed jet enters the chamber and mixes with the stationary liquid.A momentum and energy exchange occur between them, thus generating the cavitation vortex ring.There are two types of cavitation vortex rings: one is the central vortex ring formed in the shape of a ring near the central axis of the chamber, and the other is the secondary vortex ring near the collision wall, as shown in Figure 2a.The central vortex ring accumulates the jet flow at the exit of the upper nozzle.When energy accumulates to some extent, the central vortex ring breaks, and the accumulated energy is released, thus increasing the speed of the jet, as shown in Figure 2b.The cyclic occurrence of the above processes will create a pulsed jet inside the nozzle [25].As there is an angle between the collision wall and the chamber wall, a vortex will be formed when the high-speed jet flows through the collision wall returning upstream, which is called the secondary vortex ring.However, the secondary vortex ring will interact with the central vortex ring and cause energy dissipation, thus limiting the development of the central vortex ring, which is very unfavorable to the pulse intensity of the SEON.To reduce the generation of secondary vortex rings, in this paper, the wall profile of the SEON chamber was designed as a circular arc curve, as shown in Figure 3.To facilitate the flow of the jet in the nozzle, the wall surface of the chamber is tangent to the circular arc curve.However, the wall structure of the chamber will reduce  The high-speed jet enters the chamber from the upper nozzle and cavitation incepts in the chamber.Some of the cavitation bubbles leave the chamber from the lower nozzle, and the remainder of the cavitation bubbles return to the upstream part of the chamber through the collision wall to form the cavitation vortex ring.The cavitation vortex ring blocking the jet is the primary reason for the pulse generated by the SEON.The cavitation vortex ring can affect the pulse intensity.Figure 2 shows the accumulation and release process of the jet energy by the cavitation vortex ring in the SEON [2,23].The high-speed jet enters the chamber and mixes with the stationary liquid.A momentum and energy exchange occur between them, thus generating the cavitation vortex ring.There are two types of cavitation vortex rings: one is the central vortex ring formed in the shape of a ring near the central axis of the chamber, and the other is the secondary vortex ring near the collision wall, as shown in Figure 2a.The central vortex ring accumulates the jet flow at the exit of the upper nozzle.When energy accumulates to some extent, the central vortex ring breaks, and the accumulated energy is released, thus increasing the speed of the jet, as shown in Figure 2b.The cyclic occurrence of the above processes will create a pulsed jet inside the nozzle [25].As there is an angle between the collision wall and the chamber wall, a vortex will be formed when the high-speed jet flows through the collision wall returning upstream, which is called the secondary vortex ring.However, the secondary vortex ring will interact with the central vortex ring and cause energy dissipation, thus limiting the development of the central vortex ring, which is very unfavorable to the pulse intensity of the SEON.To reduce the generation of secondary vortex rings, in this paper, the wall profile of the SEON chamber was designed as a circular arc curve, as shown in Figure 3.To facilitate the flow of the jet in the nozzle, the wall surface of the chamber is tangent to the circular arc curve.However, the wall structure of the chamber will reduce the space of the chamber with the increase in the radius of the circular arc curve.Therefore, the radius of the circular arc curve needs to be limited.In this paper, three circular arc curve chamber SEONs with circular arc radii o1 mm, 2 mm, and 3 mm were designed based on the broken-line chamber SEON.
the space of the chamber with the increase in the radius of the circular arc curve.T fore, the radius of the circular arc curve needs to be limited.In this paper, three cir arc curve chamber SEONs with circular arc radii o1 mm, 2 mm, and 3 mm were des based on the broken-line chamber SEON.

Modeling
The cavitation jet is a two-phase flow.In this paper, the mixture model was ado The governing equations are as follows.

𝜕 𝜕𝑡
+  ⋅   ⃗ = 0 the space of the chamber with the increase in the radius of the circular arc curve.Therefore, the radius of the circular arc curve needs to be limited.In this paper, three circular arc curve chamber SEONs with circular arc radii o1 mm, 2 mm, and 3 mm were designed based on the broken-line chamber SEON.

Modeling
The cavitation jet is a two-phase flow.In this paper, the mixture model was adopted.The governing equations are as follows.

Modeling
The cavitation jet is a two-phase flow.In this paper, the mixture model was adopted.The governing equations are as follows.
where → v m is the average speed for quality, m/s; ρ m is the density of the mixed phase, kg/m 3 ; µ m is the viscosity of the mixed phase, Pa • s; With the improvement in computational accuracy, the Reynolds-averaged Navier-Stokes (RANS) has been unable to reach the required accuracy in some studies.Although the large-eddy simulation (LES) can get more accurate calculation results, its calculation time is huge, several orders of magnitude higher than the RANS.The high computational cost limits the wide application of the LES.To solve the problems of these two models, the stress-blended eddy simulation (SBES) was adopted, which not only guaranteed the calculation accuracy but also saved the calculation time [26].The SBES achieves stress mixing between RANS and LES through the following equation: where τ RANS ij and τ LES ij are the stress tensors of RANS and LES, respectively.At present, the common cavitation models include the Singhal model, the Zwart-Gerber-Belamri model, and the Schnerr-Saure model.These models are derived from the Rayleigh-Plesset equation, which is the fundamental equation of cavitation dynamics.In this paper, the Zwart-Gerber-Belamri model was used as the cavitation model.The mass transfer expression of this model is as follows: When When where R B is the radius of the cavitation bubble, m; P v is the saturation pressure, Pa; P is the environmental pressure, Pa; α nuc is the volume fraction of the gas nucleus; ρ v and ρ l are the density of the gas phase and liquid phase, kg/m 3 ; F vap is the evaporation coefficient; F cond is the condensation coefficient.
For the parameter setting of the cavitation model, the radius of the cavitation bubble is set to be 10 −3 mm.Water is the first phase; the density is 998.2 kg/m 3 , and the viscosity is 1.003 × 10 −3 kg/m•s.The vapor is the second phase with a density of 0.5542 kg/m 3 and a viscosity of 1.34 × 10 −5 kg/m.The evaporation pressure of the water is set at 3450 Pa.The inlet boundary condition is the pressure inlet, and the inlet pressure of the nozzle is 20 MPa.The outlet boundary condition is the pressure outlet, and the outlet pressure of the nozzle is 0.1 MPa.The turbulent intensity was 5% and the turbulent viscosity ratio was 10.The other walls are non-slip walls.The coupled method was chosen as the calculation method.The volume fraction transport equation selects the first-order upwind equation.The continuity equation and the turbulent transport equation are analyzed by the second-order upwind equation.As cavitation is a complex unsteady phenomenon, the transient simulation was adopted in this paper.The time step length is 0.00001 s.The maximum number of iterations for each time step is 20.
To ensure the accuracy of the numerical results, a three-dimensional calculation model was used in this paper.The grid adopts a hexahedral structure grid.In addition, the grid is locally refined.The meshing of the circular arc curve chamber SEON is shown in Figure 4.In order to avoid interference with the number of grids in the calculation results and reduce the amount of computation as much as possible, grid independence verification is required.The number of grids increased from 462,813 to 799,182, and five nozzles with different numbers were calculated.The upper nozzle outlet center (point A) and the lower nozzle inlet center (point B), respectively, of the SEON were selected as reference points.The deviation was calculated by

𝑥
Table 3 shows the results of the grid independence verification.When the grid number gradually increases from 462,813 to 799,182, the velocity deviation and the pressure deviation of the reference points gradually decrease.When the grid number increases from 641,547 to 799,182, both the velocity deviation and the pressure deviation are less than 0.1%, which is negligible.Therefore, the number of grids was determined to be about 641,547 in this paper.

Effect of Circular Arc Radius on Cavitation Performance of SEON
Vapor volume distribution, vorticity, and turbulent kinetic energy (TKE) are important evaluation indexes of the cavitation performance of the nozzle.The effect of the circular arc radius on the cavitation performance of the SEON was investigated when the circular arc radii were 1 mm, 2 mm, and 3 mm, respectively.
The vapor volume distribution and the vapor volume fraction can most intuitively reflect the intensity of the cavitation vortex ring.Figure 5 shows the vapor volume distribution of the nozzle chamber with different circular arc radii.When the circular arc radius increases from 1 mm to 3 mm, the vapor phase region increases first and then decreases.Table 3 shows the results of the grid independence verification.When the grid number gradually increases from 462,813 to 799,182, the velocity deviation and the pressure deviation of the reference points gradually decrease.When the grid number increases from 641,547 to 799,182, both the velocity deviation and the pressure deviation are less than 0.1%, which is negligible.Therefore, the number of grids was determined to be about 641,547 in this paper.Vapor volume distribution, vorticity, and turbulent kinetic energy (TKE) are important evaluation indexes of the cavitation performance of the nozzle.The effect of the circular arc radius on the cavitation performance of the SEON was investigated when the circular arc radii were 1 mm, 2 mm, and 3 mm, respectively.
The vapor volume distribution and the vapor volume fraction can most intuitively reflect the intensity of the cavitation vortex ring.Figure 5 shows the vapor volume distribution of the nozzle chamber with different circular arc radii.When the circular arc radius increases from 1 mm to 3 mm, the vapor phase region increases first and then decreases.With the increase in the arc radius, the chamber volume will decrease, which will affect the development of the central vortex ring.Although the circular arc radius can enhance the cavitation vortex ring in the chamber, it can also shorten the length of the collision wall.This will reduce the buffering effect of the collision wall on the cavitation jet at the upper nozzle outlet.The vapor volume fraction can represent the cavitation intensity.Figure 6 shows the variation of the vapor volume fraction in the nozzle chamber with time under different circular arc radii.When the circular arc radius increases from 1 mm to 2 mm, the vapor volume fraction increases.When the radius of the circular arc continues to increase to 3 mm, the curvature of the nozzle chamber becomes smaller, resulting in a sharp reduction of the chamber volume.This variation is very unfavorable to the formation of the central vortex ring.At this time, the vapor volume fraction of the nozzle decreases, and the cavitation performance weakens.
the development of the central vortex ring.Although the circular arc radius can enhance the cavitation vortex ring in the chamber, it can also shorten the length of the collision wall.This will reduce the buffering effect of the collision wall on the cavitation jet at the upper nozzle outlet.The vapor volume fraction can represent the cavitation intensity.Figure 6 shows the variation of the vapor volume fraction in the nozzle chamber with time under different circular arc radii.When the circular arc radius increases from 1 mm to 2 mm, the vapor volume fraction increases.When the radius of the circular arc continues to increase to 3 mm, the curvature of the nozzle chamber becomes smaller, resulting in a sharp reduction of the chamber volume.This variation is very unfavorable to the formation of the central vortex ring.At this time, the vapor volume fraction of the nozzle decreases, and the cavitation performance weakens.Vorticity is one of the most important physical parameters in describing vortex motion.Under the action of shear stress, a vortex is formed in the shear layer of the outflow field jet, and a low-pressure region is formed in the center of the vortex.Cavitation will occur when the pressure is below the saturated vapor pressure.Therefore, the more vortices formed in the shear layer of the outflow field jet, the easier cavitation is to form. the cavitation vortex ring in the chamber, it can also shorten the length of the collision wall.This will reduce the buffering effect of the collision wall on the cavitation jet at the upper nozzle outlet.The vapor volume fraction can represent the cavitation intensity.Figure 6 shows the variation of the vapor volume fraction in the nozzle chamber with time under different circular arc radii.When the circular arc radius increases from 1 mm to 2 mm, the vapor volume fraction increases.When the radius of the circular arc continues to increase to 3 mm, the curvature of the nozzle chamber becomes smaller, resulting in a sharp reduction of the chamber volume.This variation is very unfavorable to the formation of the central vortex ring.At this time, the vapor volume fraction of the nozzle decreases, and the cavitation performance weakens.Vorticity is one of the most important physical parameters in describing vortex motion.Under the action of shear stress, a vortex is formed in the shear layer of the outflow field jet, and a low-pressure region is formed in the center of the vortex.Cavitation will occur when the pressure is below the saturated vapor pressure.Therefore, the more vortices formed in the shear layer of the outflow field jet, the easier cavitation is to form.Vorticity is one of the most important physical parameters in describing vortex motion.Under the action of shear stress, a vortex is formed in the shear layer of the outflow field jet, and a low-pressure region is formed in the center of the vortex.Cavitation will occur when the pressure is below the saturated vapor pressure.Therefore, the more vortices formed in the shear layer of the outflow field jet, the easier cavitation is to form. Figure 7 shows the vorticity distribution in the chamber with different circular arc radii.When the circular arc radius is 1 mm, the vorticity is mostly concentrated in the lower nozzle, and the vorticity is small.When the circular arc radius increases to 2 mm, the vorticity in the chamber increases, and the vorticity distribution moves to the upstream part of the chamber.A large number of vortex motions occur at the upper nozzle outlet.This indicates that the cavitation vortex rings increase at the upper nozzle outlet.However, the vorticity distribution is more dispersed when the arc radius is 3 mm.In addition, there is no large vorticity distribution at the upper nozzle outlet, and the vorticity decreases.Figure 8 shows the isosurface of vorticity in the chamber with different circular arc radii, which can reflect the detailed morphological characteristics of the cavitation vortex rings.When the circular arc radius is 1 mm, the vortex rings are evenly distributed along the jet direction, and large-scale vortex rings are generated near the collision wall.In addition, the distribution area of the vorticity in the chamber increases, and the size of the vortex ring also increases obviously when the circular arc radius is 2 mm.This indicates that there is a violent vortex movement in the chamber at this time.When the circular arc radius continues to increase, the vorticity area in the chamber decreases, and the size of the vortex ring becomes smaller.Especially at the outlet of the upper nozzle, the variation is most obvious.The reason for this phenomenon is that the length of the collision wall decreases and the impact of the jet on the collision wall increases with the increase of the circular arc radius.This will cause an energy disturbance in the chamber, and a relatively stable large-scale cavitation vortex ring cannot be formed.
ticity in the chamber increases, and the vorticity distribution moves to the upstream part of the chamber.A large number of vortex motions occur at the upper nozzle outlet.This indicates that the cavitation vortex rings increase at the upper nozzle outlet.However, the vorticity distribution is more dispersed when the arc radius is 3 mm.In addition, there is no large vorticity distribution at the upper nozzle outlet, and the vorticity decreases.Figure 8 shows the isosurface of vorticity in the chamber with different circular arc radii, which can reflect the detailed morphological characteristics of the cavitation vortex rings.When the circular arc radius is 1 mm, the vortex rings are evenly distributed along the jet direction, and large-scale vortex rings are generated near the collision wall.In addition, the distribution area of the vorticity in the chamber increases, and the size of the vortex ring also increases obviously when the circular arc radius is 2 mm.This indicates that there is a violent vortex movement in the chamber at this time.When the circular arc radius continues to increase, the vorticity area in the chamber decreases, and the size of the vortex ring becomes smaller.Especially at the outlet of the upper nozzle, the variation is most obvious.The reason for this phenomenon is that the length of the collision wall decreases and the impact of the jet on the collision wall increases with the increase of the circular arc radius.This will cause an energy disturbance in the chamber, and a relatively stable large-scale cavitation vortex ring cannot be formed.When the circular arc radius is 1 mm, the vortex rings are evenly distributed along the jet direction, and large-scale vortex rings are generated near the collision wall.In addition, the distribution area of the vorticity in the chamber increases, and the size of the vortex ring also increases obviously when the circular arc radius is 2 mm.This indicates that there is a violent vortex movement in the chamber at this time.When the circular arc radius continues to increase, the vorticity area in the chamber decreases, and the size of the vortex ring becomes smaller.Especially at the outlet of the upper nozzle, the variation is most obvious.The reason for this phenomenon is that the length of the collision wall decreases and the impact of the jet on the collision wall increases with the increase of the circular arc radius.This will cause an energy disturbance in the chamber, and a relatively stable large-scale cavitation vortex ring cannot be formed.Due to the phase change, the cavitation bubbles absorb a large amount of energy when they incept and release the same amount of energy when they collapse.The inception and collapse of the cavitation bubble are accompanied by the variation of the TKE in the flow field.Therefore, the variation of the TKE can reflect the intensity of the cavitation.Figure 9 shows the variations of the TKE in the chamber with time under different circular arc radii.When the circular arc radius is 1 mm, the TKE changes dramatically with time.This is because the cavitation vortex ring in the chamber has little influence on the pulse intensity of the nozzle at the initial stage of cavitation generation.With the development of the central vortex ring, the pulse intensity of the nozzle is enhanced.When the circular arc radius increases to 2 mm, the TKE increases, and the vortex motion becomes more intense.
When the circular arc radius increases from 2 mm to 3 mm, the cavitation intensity and the TKE will decrease due to the smaller chamber space.
tion and collapse of the cavitation bubble are accompanied by the variation of the TKE in the flow field.Therefore, the variation of the TKE can reflect the intensity of the cavitation.Figure 9 shows the variations of the TKE in the chamber with time under different circular arc radii.When the circular arc radius is 1 mm, the TKE changes dramatically with time.This is because the cavitation vortex ring in the chamber has little influence on the pulse intensity of the nozzle at the initial stage of cavitation generation.With the development of the central vortex ring, the pulse intensity of the nozzle is enhanced.When the circular arc radius increases to 2 mm, the TKE increases, and the vortex motion becomes more intense.When the circular arc radius increases from 2 mm to 3 mm, the cavitation intensity and the TKE will decrease due to the smaller chamber space.

Comparative Study on Circular Arc Curve Chamber SEON and Broken-Line Chamber SEON
Taking the vapor volume distribution, the vorticity distribution, and the TKE as the comparison parameters, the circular arc curve chamber SEON with the circular arc radius of 2 mm was compared with the broken-line chamber SEON.The isosurface diagram of the vapor volume distribution can reflect not only the size of the vapor phase region but also the morphological characteristics of the cavitation vortex ring.Figures 10 and 11 are the isosurface diagrams of the vapor volume distribution of the broken-line chamber and the circular arc curve chamber, respectively.It can be seen that the central vortex ring forms a circle in the chamber.The position of the central vortex ring in the nozzle of the broken-line chamber is closer to the collision wall.As time goes by, the central vortex ring moves to the upstream of the chamber.In Figure 10d, the collapse of the central vortex ring has basically been completed, and the jet flow completes the accumulation of the energy.The central vortex ring has a long round-trip distance from the downstream to the upstream, and the energy is dissipated in the process.In Figure 11, the position of the central vortex ring in the circular arc curve chamber is closer to the upper nozzle outlet, and the scale of the central vortex ring is larger than that of the broken-line chamber.In addition, the circular arc curve chamber reduces the space of the chamber and the roundtrip distance of the vortex ring.Figure 12 shows the comparison of the vapor volume fraction between the broken-line chamber and the circular arc curve chamber.In the figure, the vapor volume fraction of the broken-line chamber ranges from 0.017 to 0.055, and its

Comparative Study on Circular Arc Curve Chamber SEON and Broken-Line Chamber SEON
Taking the vapor volume distribution, the vorticity distribution, and the TKE as the comparison parameters, the circular arc curve chamber SEON with the circular arc radius of 2 mm was compared with the broken-line chamber SEON.The isosurface diagram of the vapor volume distribution can reflect not only the size of the vapor phase region but also the morphological characteristics of the cavitation vortex ring.Figures 10 and 11 are the isosurface diagrams of the vapor volume distribution of the broken-line chamber and the circular arc curve chamber, respectively.It can be seen that the central vortex ring forms a circle in the chamber.The position of the central vortex ring in the nozzle of the broken-line chamber is closer to the collision wall.As time goes by, the central vortex ring moves to the upstream of the chamber.In Figure 10d, the collapse of the central vortex ring has basically been completed, and the jet flow completes the accumulation of the energy.The central vortex ring has a long round-trip distance from the downstream to the upstream, and the energy is dissipated in the process.In Figure 11, the position of the central vortex ring in the circular arc curve chamber is closer to the upper nozzle outlet, and the scale of the central vortex ring is larger than that of the broken-line chamber.In addition, the circular arc curve chamber reduces the space of the chamber and the round-trip distance of the vortex ring.Figure 12 shows the comparison of the vapor volume fraction between the broken-line chamber and the circular arc curve chamber.In the figure, the vapor volume fraction of the broken-line chamber ranges from 0.017 to 0.055, and its average value is 0.040.However, the vapor volume fraction of the circular arc curve chamber ranges from 0.030 to 0.112, with an average value of 0.089.The above data indicate that the cavitation area of the circular arc curve chamber SEON is increased by 122.5% compared with that of the broken-line chamber SEON.
average value is 0.040.However, the vapor volume fraction of the circular arc curve chamber ranges from 0.030 to 0.112, with an average value of 0.089.The above data indicate that the cavitation area of the circular arc curve chamber SEON is increased by 122.5% compared with that of the broken-line chamber SEON.    Figure 13 shows the distribution of vortex rings in the broken-line chamber and the circular arc curve chamber.As shown in Figure 13a, the central vortex ring is generated in the chamber near the upper nozzle outlet, which can squeeze and block the jet flow at the exit of the upper nozzle.The angle between the collision wall and the chamber wall is altered, and the vortex is formed when the jet flows back upstream through the collision wall.A secondary vortex ring on a smaller scale is generated near the collision wall.The secondary vortex ring is formed in the process of energy dissipation of the jet flow.When the scale is large enough, it will offset the energy of the central vortex ring and weaken the development of the central vortex ring.The distribution of vortex rings inside the circular arc curve chamber nozzle is shown in Figure 13b.It can be seen that the secondary vortex rings at the collision wall are significantly weakened, and their interference with the central vortex ring decreases.This makes the central vortex ring at the upper nozzle outlet effectively enhance the pulse intensity of the SEON.The vorticity distribution of the nozzle outflow field is an intuitive representation of the nozzle vortex motion.In order to make a further comparative study on circular arc curve chamber SEON and broken-line chamber SEON, the same modeling described above was adopted to investigate the outflow field of the two nozzles.Figure 14 is the isosurface diagram of the vorticity distribution of the nozzle outflow field.The energy of the jet flow from the broken-line chamber SEON is limited.The vortex sheds substantially at the end of the jet flow, as shown in Figure 14a.However, the vorticity distribution of the circular arc curve chamber SEON is much higher than that of the broken-line chamber SEON, especially at the end of

Conclusions
In order to investigate the effect of the chamber wall profile on the cavitation performance of the SEON, a novel circular arc curve chamber SEON was designed.The influence of the circular arc radius on the cavitation performance of the circular arc curve chamber SEON was analyzed by a numerical method.Additionally, a comparative study was con-

Conclusions
In order to investigate the effect of the chamber wall profile on the cavitation performance of the SEON, a novel circular arc curve chamber SEON was designed.The influence of the circular arc radius on the cavitation performance of the circular arc curve chamber SEON was analyzed by a numerical method.Additionally, a comparative study was conducted between the broken-line chamber SEON and the circular arc curve chamber SEON.The following conclusions were obtained: (1) With the increase in the circular arc radius, the cavitation performance of the circular arc curve chamber SEON is first enhanced and then weakened.When the circular arc radius is 2 mm, the vortex ring in the chamber reaches the maximum, and the cavitation performance of the nozzle is the best.
(2) Compared to the broken-line chamber SEON, the cavitation performance of the circular arc curve chamber SEON is better.The numerical results show that the cavitation area of the circular arc curve chamber SEON increases by 122.5% and the TKE increases by 16.9%.

Figure 1 .
Figure 1.Structural diagram of the SEON.

Figure 1 .
Figure 1.Structural diagram of the SEON.

Figure 2 .
Figure 2. Schematic diagram of the accumulation and release process of the jet energy by th tation vortex ring: (a) accumulation process; (b) release process.

Figure 3 .
Figure 3. Profile of inner wall of the SEON.

Figure 2 .
Figure 2. Schematic diagram of the accumulation and release process of the jet energy by the cavitation vortex ring: (a) accumulation process; (b) release process.

Figure 2 .
Figure 2. Schematic diagram of the accumulation and release process of the jet energy by the cavitation vortex ring: (a) accumulation process; (b) release process.

Figure 3 .
Figure 3. Profile of inner wall of the SEON.

Figure 3 .
Figure 3. Profile of inner wall of the SEON.

Figure 4 .
Figure 4.The meshing of the circular arc curve chamber SEONs.

Figure 4 .
Figure 4.The meshing of the circular arc curve chamber SEONs.

Figure 6 .
Figure 6.Variations of vapor volume fraction in nozzle chamber with time.

Figure 6 .
Figure 6.Variations of vapor volume fraction in nozzle chamber with time.

Figure 9 .
Figure 9. Variations of TKE in nozzle chamber with time.

Figure 9 .
Figure 9. Variations of TKE in nozzle chamber with time.

Figure 13 Figure 13 .
Figure13shows the distribution of vortex rings in the broken-line chamber and the circular arc curve chamber.As shown in Figure13a, the central vortex ring is generated in the chamber near the upper nozzle outlet, which can squeeze and block the jet flow at the exit of the upper nozzle.The angle between the collision wall and the chamber wall is altered, and the vortex is formed when the jet flows back upstream through the collision wall.A secondary vortex ring on a smaller scale is generated near the collision wall.The secondary vortex ring is formed in the process of energy dissipation of the jet flow.When the scale is large enough, it will offset the energy of the central vortex ring and weaken the development of the central vortex ring.The distribution of vortex rings inside the circular arc curve chamber nozzle is shown in Figure13b.It can be seen that the secondary vortex rings at the collision wall are significantly weakened, and their interference with the central vortex ring decreases.This makes the central vortex ring at the upper nozzle outlet effectively enhance the pulse intensity of the SEON.The vorticity distribution of the nozzle outflow field is an intuitive representation of the nozzle vortex motion.

Figure 12 .
Figure 12.Comparison of vapor volume fraction of two nozzles.

Figure 13
Figure13shows the distribution of vortex rings in the broken-line chamber and the circular arc curve chamber.As shown in Figure13a, the central vortex ring is generated in the chamber near the upper nozzle outlet, which can squeeze and block the jet flow at the exit of the upper nozzle.The angle between the collision wall and the chamber wall is altered, and the vortex is formed when the jet flows back upstream through the collision wall.A secondary vortex ring on a smaller scale is generated near the collision wall.The secondary vortex ring is formed in the process of energy dissipation of the jet flow.When the scale is large enough, it will offset the energy of the central vortex ring and weaken the development of the central vortex ring.The distribution of vortex rings inside the circular arc curve chamber nozzle is shown in Figure13b.It can be seen that the secondary vortex rings at the collision wall are significantly weakened, and their interference with the central vortex ring decreases.This makes the central vortex ring at the upper nozzle outlet effectively enhance the pulse intensity of the SEON.The vorticity distribution of the nozzle outflow field is an intuitive representation of the nozzle vortex motion.

Figure 15 .
Figure 15.Comparison of TKE of two nozzles.

Figure 15 .
Figure 15.Comparison of TKE of two nozzles.

Table 1 .
Range of dimensionless parameters of the SEON.

Table 1 .
Range of dimensionless parameters of the SEON.

Table 2 .
Structural parameters of the SEON.

Table 2 .
Structural parameters of the SEON.

. Results and Discussion 4
.1.Effect of Circular Arc Radius on Cavitation Performance of SEON