Flow Characteristics of Two-Dimensional Supersonic Under-Expanded Coanda-Reattached Jet

When there is a wall near the jet, it deflects and flows while being attached to the wall owing to the Coanda effect. The flow characteristics of the incompressible and two-imensional (2D) Coanda-reattached jets have been considerably explained. However, 2D supersonic under-expanded jets, reattached to side walls, have not been sufficiently investigated. These jets are used in gas-atomization to produce fine metal powder particles of several micrometers to several tens micrometers. In this case, the supersonic under-expanded jets are issued from an annular nozzle, which is set around a vertically installed circular nozzle for molten metal. The jet flow at the center cross-section of the annular jet resembles a 2D Coanda-reattached jet that deflects and attaches on the central axis. In this study, the flow characteristics of a supersonic under-expanded Coanda air jet from a 2D nozzle that reattaches to an offset side wall are elucidated through experiment and numerical analysis. For numerical analysis, we show how much it can express experimental results. The effects of supply pressure P0 on the flow characteristics such as the flow pattern, size of shock cell, reattachment distance, and velocity and pressure distributions, etc. are examined. The flow pattern was visualized by Schlieren method and the velocity distribution was measured using a Pitot tube. These results will be also useful in understanding the flow characteristics of a gas-atomization annular nozzle approximately.


Introduction
When there is a wall near the jet, the jet will adhere to and flow round nearby the wall. This Coanda-reattached flow is used for controlling flow in a wall attached type fluidic device, preventing flow separation, and increasing the wing lift force. The flow characteristics of incompressible two-dimensional (2D) Coanda-reattached jet flows have been elucidated considerably owing to the development of fluidics [1] [2]. For example, Bourque et al. [1] studied the reattachment distance and pressure within the separation bubble of a 2D reattached jet by a theoretical analysis using momentum theory, and compared them with experiments in terms of the offset distance and inclined angle of the side wall.
Supersonic under-expanded free jets, such as jets from jet engines and rocket engines to obtain thrust have been studied in detail. Donaldson et al. [3] classified the flow patterns of sub-and supersonic under-expanded jets into three, depending on the pressure ratio of the supply pressure to the ambient pressure.
Kojima et al. [4], Zapryagaev et al. [5], and André et al. [6] studied the structure of shock waves, the pressure and velocity distribution of supersonic under-expanded jets. Shadow et al. [7] and Zaman [8] showed the compressible spreading rates and mixing characteristics of supersonic coaxial jets. Katanoda et al. [9] [10] explained the effects of the nozzle divergence angle and Mach number on the shock cell length, and Sugawara et al. [11] showed the structure of three-dimensional microjet with a Mach disc by Mach-Zender in interferometer.
Zaman [8] examined the flow characteristics of sub-and supersonic jets from special nozzles, whereas Shakouchi [12] studied those from an orifice nozzle.
Franquet et al. [13] reviewed the studies on free under-expanded jets in a quiescent medium. However, studies regarding supersonic under-expanded Coanda-reattached jet flows and their flow characteristics are rare [14] [15].
For example, supersonic under-expanded jets are used in gas-atomization, producing fine metal particles of several micrometers to several tens of micrometers [16] [17] [18]. The molten metal from a circular nozzle is broken, subdivided, and refined by a shearing force caused by the large velocity gradient outside the supersonic under-expanded jet issued from a coaxial annular nozzle installed near the circular nozzle, and high frequency fluctuations caused by the flow (Figure 1). The atomized molten metal, which is spherical owing to surface tension, is subsequently cooled to produce fine metal powder particles. In this case, the supersonic under-expanded jet from the annular nozzle curves inward, deflects, and flows over the cone shape collectively because the pressure inside the annular jet decreases owing to jet entrainment, and the jet flow at the center section of the x-y plane of the annular jet resembles a 2D Coanda-reattached jet that deflects and attaches on the central axis. To improve the production efficiency of the gas-atomization and control the particle size and distribution, the flow characteristics of the supersonic under-expanded Coanda-reattached jet from an annular nozzle must be studied.  In this study, the flow characteristics of a supersonic under-expanded Coanda air jet from a 2D nozzle that reattaches to an offset side wall are elucidated through experimental and numerical analyses. For numerical analysis, we show how much it can express experimental results. The effects of supply pressure P 0 on the flow characteristics such as the flow pattern, size of shock cell, reattachment distance, and velocity and pressure distributions are examined. The flow pattern was visualized using the Schlieren method and the velocity distribution was measured using a Pitot tube. P 0 was changed between 0.2 MPa to 0.5 MPa.
These results will also be useful in understanding the flow characteristics of a gas-atomization annular nozzle approximately. Figure 2 shows the flow model and an illustration of how the 2D air jet issued from the nozzle of width b reattaches to the side wall with an offset distance D at the reattaching distance x R while surrounding a separation bubble region of pressure P b . In this case, the jet centerline can be approximated by a circular arc of radius R with its origin on the y-axis [1].

Numerical Analysis
The CFD program, SOLIDWORKS Flow Simulation, was used for the numerical analysis. The governing equations obtained using the k-ε turbulence model, and the Favre-averaged Navier-Stokes equation were analyzed based on the finite sedimentation and participation methods. mesh. An orthogonal coordinate system was used, and the control volume was a rectangular parallelepiped. The cells in the vicinity of the shape boundary, domain region, narrow nozzle passage, center of the jet, and the entrainment re-gion of the jet are subdivided into different levels. The cells are arbitrarily cut by the shape boundary, and one parallel six-specimen cell has a plurality of different control volumes, such as one solid and one fluid. In addition, the nozzle width, the length of the parallel section of the nozzle, and offset distance were b, L w /b = 4.0, and D/b = 5.0, respectively. The upper left of the flow channel contained an opening, for the molten metal for gas-atomization, but it was closed and calculated in this study. This was assumed not to affect the flow characteristics of the reattached jet. Calculations were performed primarily for the 2D channel and the channel with a nozzle AR = 4.0. The supply boundary entrance condition, that is, the supply air setting, was set to P 0 = 0.2 -0.5 MPa at the static pressure setting. The exit boundary condition was set to atmospheric pressure.
After investigated the effects of the number of cells on the flow, 1.13 × 10 6 cells were used.

Two-Dimensional Jet
Figure 4(a) shows the density distribution ρ/ρ 0 from numerical analysis, on the x-y plane of the supersonic under-expanded reattached jet issued from the 2D nozzle of width b. The offset distance is D/b = 5.0 and the supply pressure is P 0 = 0.4 MPa. In Figure 4(a), the length C L and width C W of the 1 st shock cell are also shown. Their details will be shown in the section 6.2. The jet from the 2D nozzle expands and compresses to form shock cells while generating shock waves and then flows away. In this case, owing to flow entrainment, the pressure on the side wall of the supersonic under-expanded jet decreases, and the jet is deflected and reattached to the side wall and flows down while forming a separation bubble. The reattaching distance x R , i.e., the value of x at the reattachment point ( Figure 2), also known as the reattaching distance, is x R /b = 16.1. Furthermore, the minimum pressure and the position on the side wall are P wmin /P 0 = −0.046, and x/b = 8.5, respectively. In addition, the jet centerline can be expressed by an arc of radius R/b = 35.0. Figure 4(b) shows the density distribution ρ/ρ 0 , from numerical analysis, of the x-y plane at the center height of the flow channel of the supersonic under-expanded reattached air jet issued from a rectangular nozzle with aspect ratio AR = H/b = 4.0 (H: nozzle height). The flow patterns in Figure 4(a) and Figure 4(b) are approximately the same, and the reattaching distance x R , minimum bubble pressure P b , and radius R are also approximately the same; x R /b = 16.1, 15.8, P wmin = −0.046, −0.046, and R/b = 35.0, 36.5, respectively. Consequently, it may be assumed that the flow state in the x-y section at the center height of the flow channel with AR = 4.0 in Figure 4   A static pressure hole, φ = 0.8 mm, was provided at the center height of the channel before the nozzle, which was used to measure the supply pressure P 0 . The supply pressure varied to a maximum of P 0 = 0.5 MPa; subsequently, flow visualization and velocity and pressure measurements were performed. A Schlieren optical system (Kato Koken, Japan, SS150) was used for flow visualization. A total or static Pitot tube with a diameter of 1.0 mm was used to  y-axis when the jet center was approximated using a circular arc. Bourque et al. [1] showed that the centerline of a 2D incompressible reattached jet can be approximated by a circular arc. It is assumed that the jet centerline of the supersonic under-expanded reattached jet can also be approximated using a circular arc.  3)] of the shock cell increases with P 0 , radius R of the circular arc, and the reattachment distance x R . Figures 6(b)(1)- (4) show the density distributions ρ/ρ 0 (from numerical analysis) for P 0 = 0.2 -0.5 MPa, respectively. Like the density distribution, the flow pattern, C L , C W , R, and x R agree with the experimental results (details will be described later). That is,

Shock Cell Size
The jet issued from the nozzle formed a shock cell after expansion and compression; thereafter, while expansion and compression continued, it curved and reattached to the side wall and flowed downstream. Figure 7 shows the size, length C L , and maximum width C W of the first shock cell (Figure 4(a)). The abscissa is P 0 and P 0a /P aa , where P 0a and P aa are supply and ambient pressures expressed in the absolute pressure, respectively. The C L was obtained as the distance from the nozzle exit to the position where the measured jet centerline velocity is initially minimal (Figure 11). The C W was obtained from the visualized flow pattern and the density distribution. At that time, the jet boundary was set at the position where the shade suddenly changed in the jet width direction.
Both C L and C W increase almost linearly with P 0 . The numerical analysis is consistent with the experimental results. In addition, the sizes of the second and third shock cells were similar to the first one.

Centerline of the Jet and Reattached Distance
As mentioned in Section 5.1, the jet centerline can be represented by a circular arc.  Figure 9 shows the pressure distribution P w /P 0 on the side wall for P 0 = 0.2 -0.5

Pressure Distribution on Side Wall
MPa. The jet issued from the nozzle entrains the surrounding fluid owing to the large velocity gradient and the fluid viscosity. When a solid wall exists near the jet, it entrains and removes the fluid from the finite area between the jet and the wall. Consequently, the pressure decreases, and the jet deflects and reattaches, enclosing the circulating separation bubble to the side wall. The pressure P w /P 0 on the side wall has a negative value corresponding to the separation bubble; it decreases with increasing x/b and reaches a minimum and approaches the ambient pressure. The minimum of the negative value and maximum of P w /P 0 increase with P 0 , and are almost the same at P 0 ≥ 0.4 MPa. This is because increasing the velocity owing to the increase in P 0 hinders the jet from entraining the surrounding fluid [12], jet flows penetrates. The numerical analysis reflected the experimental well, particularly when P 0 ≥ 0.3 MPa, that is, when P 0 is large.   Figure 10 shows the minimum pressure P wmin /P 0 on the side wall. P wmin /P 0 decreases with increasing P 0 and approaches a constant at P 0 ≥ 0.4 MPa, and the numerical analysis reflects the experiment effectively. In addition, according to Bourque et al. [1], the position x of the minimum pressure on the side wall of a 2D incompressible reattached jet corresponds to the position x of the minimum bubble pressure. This is reflected in the numerical analysis of the supersonic under-expanded reattached jet investigated in this study.

Pressure Difference between Both Sides of Reattached Jet and Jet Circular Motion
Bourque et al. [1] reported that when a 2D incompressible jet flow reattaches to a side wall set parallel to the nozzle axis at an offset distance, the jet centerline can be approximated to a circular arc of radius R (Figure 2). The relationship between the pressure difference ΔP on both sides of the jet and the centrifugal force exerted on the jet can be expressed as follows, where the inside pressure of the reattached jet is the minimum pressure of the separation bubble. .
where J is the jet momentum.
Equation (1) can be rewritten as where Q and ρ 0 are the flow rate and density of air at standard conditions, respectively. This relationship was investigated for a supersonic under-expanded reattached jet of P 0 = 0.4 MPa, assuming that the jet momentum was conserved, as in Bourque et al. [1]. Here, because the minimum pressure on the side wall, flow rate, and arc radius were P wmin = 1.58 × 10 4 [Pa], Q = 58.0 [Nm 3 /h], and R = 86.5 mm, respectively. Inserting the jet momentum derived from Q and R into Equation (1) yields ΔP = 1.55 × 10 4 [Pa]. This is almost the same as P wmin , which is reflected in the case of the supersonic reattached jet. Therefore, the relationship shown in Equation (1) holds well. Figure 10. Minimum surface pressure P wmin /P 0 on the side wall.
T. Ohmura et al. Figure 11 shows the numerical calculations of the centerline velocity u c of the supersonic reattached jet. The abscissa is the angle θ (angle in the circumferential direction with nozzle exit set at θ = 0˚, Figure 2.) when the jet centerline is approximated by a circular arc. The case of experiment for P 0 = 0.4 MPa is also shown in Figure 11. Based on the numerical results for P 0 = 0.4 MPa, the jet issues from the nozzle initially expanding and the u c reaches a maximum at θ = 3.25˚ and then decreases by compression to a minimum at θ = 6.0˚. It is assumed that the region from θ = 0˚ to 6.0˚ corresponds to the first shock cell. Thereafter, the u c decays gradually while repeating expansion and compression. The amplitude of the increase/decrease in u c reduces, but the amplitude of u c is moderate.    (Figure 12), respectively. The y' on the vertical axis in Figure 13 is the distance in the radius direction from the centerline of the reattached jet ( Figure 2). The cross-sectional velocity profile of a supersonic under-expanded jet differs depending on the expansion or compression regions. Kojima, T. et al. [4] showed that the velocity profile in the shock cell of a supersonic under-expanded free jet was concave initially, convex in the middle, and nipple-shaped at the end. The velocity profile at θ = 3.25˚ was concave and asymmetric, with a minimum at the jet center, which differs from Kojima's classification of convex. This could because this jet was a reattached jet, and the pressure on the separation bubble side was negative.

Equi-Velocity Distribution and Velocity Distribution at Cross-Section
The profiles at θ = 6.5˚, 12.5˚, and 18.25˚ are convex at the most compressed the first shock cell is convex, that at θ = 6.5˚ ("1st Com.") is concave, and those at θ˚ = 12.5˚, 18.25˚ ("2nd and 3rd Com.") are convex. As shown by the jet centerline velocity in point of the local minimum. Additionally, numerical calculations are   shown in Figure 13. The jet width and range of the numerical results are smaller than those of the experimental results, however, the numerical calculation reflects the approximate distribution of the experiment.

Conclusions
In this study, flow characteristics, such as shock cell formations, moving trajectory of reattached centerline approximated by a circular arc, reattachment dis- 2) The centerline of a 2D supersonic under-expanded reattached jet can be approximated by a circular arc of radius R, in which R and the reattachment distance x R increased with supply pressure P 0 . Furthermore, the numerical calculations for supersonic reattached jet represent the experimental results fairly well.
The relationship ΔP = J/R, obtained by assuming that the force owing to the pressure difference ΔP applied to both sides of the jet and the centrifugal force applied to the jet element were balanced, was applicable.
3) The length and width of the shock cell of the supersonic reattached jet increased with P 0 , and the minimum pressure on the side wall decreased with increasing P 0 . The numerical calculations represented the experimental outcomes well. 5) The above results for the 2D reattached jet will also be useful in understanding the flow characteristics from a gas-atomization nozzle approximately.