Structure Optimization and Numerical Simulation of Nozzle for High Pressure Water Jetting

Three kinds of nozzles normally used in industrial production are numerically simulated, and the structure of nozzle with the best jetting performance out of the three nozzles is optimized. The R90 nozzle displays the most optimal jetting properties, including the smooth transition of the nozzle’s inner surface. Simulation results of all sample nozzles in this study show that the helix nozzle ultimately displays the best jetting performance. Jetting velocity magnitude along Y and Z coordinates is not symmetrical for the helix nozzle. Compared to simply changing the jetting angle, revolving the jet issued from the helix nozzle creates a grinding wheel on the cleaning surface, which makes not only an impact effect but also a shearing action on the cleaning object. This particular shearing action improves the cleaning process overall and forms awider, effective cleaning range, thus obtaining a broader jet width.


Introduction
In the process of material surface cleaning by high pressure water jet, the enhancement of cleaning effectiveness and reduction of energy consumption have been significant areas of research for quite some time.Nurick (1976) studied the cavitation characteristics of sharp-edged orifices including circular orifices and rectangular orifices and their effect on spray mixing [1].Hashish and du Plessis (1978) developed a theory based on a control volume analysis, which presented a basic method for building cutting equations to evaluate hydrodynamic forces [2].Soyama et al. (1996) focused on peculiar phenomena occurring around high speed submerged water jets, clarified the effectiveness of the nozzle structure and injection pressure therein, and then discussed the observed characteristics and configuration of the water jet [3].A clear relationship between parameters (standoff distance and installation angle) and cleaning results was drawn by Gao and Chen (2006) [4].
Estimation in nozzle cleaning performance is of paramount importance in high axial velocity and radial velocity.We hold the opinion based on experiment and analysis that different nozzle configurations determined different axial velocities.Considering shearing action will function in cleaning process, we have made a contribution to increasing axial and radial velocity by optimizing nozzle configuration.This study, which focuses on nozzle configuration, posits a new nozzle and verifies its highly effective jet velocity by numerical simulation.Three kinds of nozzles commonly used in industrial production, the Cylindrical Contracting (CC) nozzle, Taper Contracting (TC) nozzle, and Stepped nozzle, are simulated by FLUENT to obtain exact jet velocity comparison data.The most applicable of these is then sequentially optimized by simulating tentative structure design models.For the sequential optimization, a 3D model is used to simulate the helix nozzle outlet structure.Comparison results from the simulation of different nozzle inner surface configurations indicate the most favorable, best performing nozzle structure.determines the cleaning process.The Cylindrical Contracting (CC) nozzle, Taper Contracting (TC) nozzle, and Stepped nozzle are those generally used in the industry currently.Figure 1 details their respective structures.The type of water jet that issues from all these different orifice nozzles has similar jet structures in air.Yanaida and Ohashi first detailed a geometrical description of the water jet in 1980 [5].

Preliminary Optimization of Nozzle
Over the past 40 years, a recognizable characteristic pattern has emerged.Jet structure, based on theory and observation, is shown in Figure 2 [6].The jet velocity is equally distributed at the nozzle but produces a considerable velocity difference as soon as it leaves the nozzle, compared to the boundary layer formed by the environment medium.The process of exchanging momentum and mass of jet medium and environment medium is the same as the process of diffusing jet fluid.The jet in the potential core region is applied to cutting, and jetting in the diffused droplet region is applied to dustlaying and aspirating, whereas the jet in the main region, the focus of this study, is applied to cleaning and surface finish.Nozzles shown in Figure 1 are widely used in engineering projects; however, higher velocity, higher dynamic pressure, lower pump pressure, and lower cost are necessary considerations.The optimization of nozzle structure offers a simple solution to these.The most suitable nozzle is selected for the study from the three kinds of nozzles listed above.
FLUENT was used as the simulation tool.Two-dimensional simulation of the jet formation of the CC nozzle, TC nozzle, and Stepped nozzle was performed first.The computational area and boundary conditions were as shown in Figure 3.The structured grid system of the nozzles was created in GAMBIT.The inlet boundary condition was a pressure inlet, and the outlet boundary condition was a pressure outlet.Water-liquid, the subject of flow in axisymmetric space, was the uncompressed ideal fluid and followed a steady flow pattern.
Figure 4 shows the main nozzle structural parameters.According to the theoretical analysis and experimental data, geometric parameters of the CC nozzle are defined as  = 13 ∘ , / = 3, and / = 4. Boundary conditions were considered the inlet pressure = 100 MPa and outlet pressure of normal atmospheric pressure.The fluent fluid material was water-liquid, and the standard - model is appropriate in our case.Nozzle outlet diameter was 0.8 mm, inlet diameter was 3.2 mm, and computational area size was 288 × 25 mm 2 to ensure that the water jet was completely divergent.The numerical methodology expressed by Tu et al. (2007) was adopted in our flow model [7].The flow velocity at the inlet and outlet could be preliminarily roughly calculated based on the Bernoulli equation: where  1 ,  2 are pressures at inlet and outlet cross sections along the horizontal pipeline of the nozzle, V 1 , V 2 are flow velocities at the two cross sections, and  denotes flow specific weight.
The equation for counting volume flux in the pipeline is then deduced: According to the definition  = V, (1) is then written as follows: The term  0 is a theoretical value of volume flux.Considering friction loss, more accurate volume flux  is obtained from  0 multiplying by a friction coefficient (set to 1.2). 1 ,  2 are areas of the two cross sections.Certain concrete values must be used to analyze our flow model, among which  1 is inlet cross section area and  2 is outlet cross section area, defined as follows: By substituting (3) into (2) and adding the actual values to these equations, we obtain Theoretical values V 1 and V 2 , slightly different compared to the values, are used for quantitative analysis in distinguishing flow type.One important use of the Reynolds number (Re) obtained by V 1 , regarded an essential nondimensional parameter that describes the flow characteristics, is to indicate whether the flow is turbulent or laminar.A channel flow will become turbulent if the critical Reynolds number is above 1400.Re is defined as the ratio of the inertia force over the friction force: It is observed that the equation encapsulates three variables, density (), dynamic viscosity (), and inlet velocity (V 1 ).(= 2) is the diameter of the nozzle inlet.Combining (1) and ( 3) into (5), we obtain Re = 9.0 × 10 4 ≫ 1400.
Thus, we confirm that the flow model is turbulent.Our case is a constant property fluid flow, indicating that flow density is constant and body force in the equation can be neglected.
The momentum equations employed are given as follows [8]: where   is the volume fraction;   is the density of water; and V  is the axial velocity of water droplets.The term ṁ →  is the mass transfer from  (water) phase.In order to simulate the air entrainment process, we set ṁ →  and   as zero and where Ṁ is known from the semiempirical mass model of Guha et al. [9]: where  0 and  0 are the volume fraction and axial velocity of water droplets at the nozzle exit, respectively, and  is the nozzle radius.The variable   is defined as where   is the radial width of the continuous flow region and is given by Outside of this region is the droplet flow region, the radial width of which varies as follows: 1 and  are spread coefficients.The value of  is usually in between 0.02 and 0.06. 1 is estimated as follows: where const is usually (10 −1 ).Schematic description of   and  0 has been provided in Figure 2.
The standard - model equations, expressed in nonconservation form, cause an increase in demand.The transport equations are as follows: where   is the mixture density and ⃗ V  is the mixture velocity.The turbulent viscosity ( , ) and the production of turbulent kinetic energy ( , ) are calculated as follows: The adjustable constants   ,   ,   ,  1 , and  2 are obtained by sufficient data statistics (Launder and Spalding, 1974) [10]: (17)

Optimization Structure of Helix Nozzle.
To better understand the principles at work in this study, let us consider the firing of a bullet.Propelled by gunpowder, the bullet accelerates and fires through air in a rotational motion.The pressure of powder gases expels the bullet, and it continues along a predictable trajectory.As shown in Figure 8, by observing the inside of the gun barrel, the bullet expelled does not move in a straight line but a spiral at the angles which correspond to the centerline of the barrel.It functions similarly to a top-heavy gyroscope, in the fact that it requires a rapid spinning motion to continue its trajectory.As it flies off the barrel, the bullet must overcome air resistance.A tight spiral trajectory ensures high velocity and straight, forward shooting direction.
If we think of the flow element as a bullet, as shown in Figure 5, we can effectively structure the nozzle as a gun barrel in order to simulate rotational movement.The formation of a spiral not only guides the flow, but also speeds the flow and keeps it steady.As shown in Figure 6, for precise simulation of the helix nozzle, the 2D model must be transformed to 3D, as 2D cannot properly represent the complexity of a spiral.The 3D is built while keeping nozzle geometry parameters and simulating boundary conditions invariable (Table 1).A structured grid system is created in the 3D model, distinguished from 2D model by two independent mesh bodies (the nozzle and computational domain) which are connected by an interface.

Results and Discussion
Nozzle structure design models with better performance were chosen to make further optimization.Table 2 shows optimization procedures of nozzle structure design models.
To compare the effect of optimization design models, average axial and radial velocities of structure models in target distance range of  = 0.08 ∼ 0.09 m are compared in Table 3.

Preliminary Optimization.
Compared to the Stepped nozzle in Figure 7, all CC nozzles and TC nozzles demonstrate better jetting performance in the range of effective target distance.Despite the fact that the highest velocity TC nozzle performed at the outlet, its velocity characteristic in the main region is similar to the CC nozzle.Evolved from the TC nozzle, the CC nozzle is a solid flow nozzle whose jetting behavior displays great concentrated force and standoff distance, where its shape proves beneficial in high precision micromachining.To this effect, the CC nozzle was chosen as the study object in this paper.Its configuration optimization design is processed in two stages.
As is commonly known, a sharp transition from the convergent section and the two cylindrical sections in the CC nozzle body would increase local resistance loss.As a solution, smooth transition of the inner surface can ensure fluid streamlining and form jetting that includes a high velocity coherent core.
Figure 8 shows the preliminary optimization structures, in which (a) is CC nozzle and (b)-(f) are optimized nozzles with different radii of circular arcs.These nozzles were then simulated in order to select the one with the best jetting performance.According to the different nozzle structures first attempted, these nozzles were successively simulated.Velocity magnitude was chosen for simulation, as a series of parameters such as velocity magnitude, axial velocity, dynamic pressure, and shooting force, all of which have positive correlation, are able to effectively measure jetting performance.
As shown in Figure 9, (b) is a local display of (a) at the range of 0.078 m to 0.10 m in the  direction.Compared with the CC nozzle, the jetting velocity magnitude of the R90 nozzle is higher than the others if the target distance was set in the range of  = 0.08 ∼ 0.09 m.A gripping tool or manipulator installed in the high pressure water jet equipment enables the nozzle to reach any range of target distance; therefore, the R90 nozzle is the proposed model.

Optimization of Helix Nozzle.
Observing the velocity vectors as shown in Figure 10, flow velocity magnitude apparently increases along the spiral line, and horizontal velocity almost matches the axial velocity of the same cross section.
Figure 11 shows that jetting velocity magnitude along  and  coordinates is not symmetrical, indicating that the high speed jet performs a shearing action on the cleaning object.In the high pressure water jet cleaning process, the jet generally loses some kinetic energy while vertically shooting at the cleaning surface.However, the jet issuing from a helix nozzle actually avoids this problem.Compared to simply changing the jetting angle, a revolving jet acts like a grinding wheel on the cleaning surface, displaying not only an impact effect but also a shearing action on the cleaning object.
In order to compare the velocity magnitude of the helix nozzle with the R90 nozzle, high pressure water jetting simulation data was loaded in MATLAB, and Lagrange Interpolation Method was used to fit the velocity magnitude curves.As shown in Table 3, R90 nozzle performs the highest average axial velocity but very low radial velocity; helix nozzle performs both comparatively high average axial velocity and radial velocity.The removal of surface materials partly depends on a considerable shearing action of water jets issuing from the helix nozzle.Furthermore, the spiral line added to the helix nozzle's inner surface had no distinctly opposite influence on axial velocity.
Cleaning effectiveness was obvious only when the velocity value was above 150 m/s.The curve denoted R90 nozzle shows an effective cleaning range of about 0.078-0.12m (the length of nozzle in our case was 0.078 m), while the effective Advances in Materials Science and Engineering cleaning range of the helix nozzle was about 0.078-0.143m.In a specific limit, a wider effective cleaning range produces a broader jet width, resulting in more effective cleaning.

Conclusion
In this paper, an effective optimization of nozzle configuration for high pressure water jetting has been performed.was chosen for subsequent experimentation.The relevant advantages and feasibility of adding a helix to the nozzle at the outlet were confirmed by applying rifle theory.
In future study, it is necessary to comprehensively analyze the quantity and basic parameters of a spiral line and its corresponding jetting divergency phenomenon in air, considering factors of cavitation and gravity.Further research on jetting performance by surface stress analysis of different cleaning objects will be discussed in the future.

Figure 7 :
Figure 7: Three different nozzles structure comparison from simulation of jet velocity magnitude.

Figure 11 :
Figure 11: Comparison of velocity vectors of R90 nozzle and helix nozzle.

Table 1 :
Values of nozzle structure parameters and boundary conditions.

Table 2 :
Optimization procedures of nozzle structure design models.

Table 3 :
Average axial and radial velocities of structure models in target distance range of  = 0.08∼0.09m.
[7]et al. (2007)y demonstrable superiority and reliability of the designed object, which shows better manifestation of dynamic pressure, cohesion, and low local resistance loss, are also provided.As previous studies show, flow around the variable-area channel can dramatically change, where particles collide and produce a whirlpool, causing local resistance loss.To remedy this, we construct a highly effective, pliable structure in the nozzle's inner surface which creates an optimized structure, based on the principle of a helical accelerator.The governing equations of turbulence proposed byTu et al. (2007)[7]were used to simulate our models.At the first step of optimization, nozzle R90 (Figure8(d)) Figure 10: Velocity vectors of helix nozzle.