 Open Access
 Total Downloads : 0
 Authors : Manoj Kumar, A. K. Gupta, P. Kumar, R. K. Sahoo
 Paper ID : IJERTCONV6IS16013
 Volume & Issue : RDME – 2018 (Volume 06 – Issue 16)
 Published (First Online): 05012019
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Numerical Analysis to Predict the Fluid Flow Pattern Through Convergent Nozzle
Numerical Analysis to Predict the Fluid Flow Pattern Through Convergent Nozzle
Manoj Kumar*, A. K. Gupta, P. Kumar
Ph.D. Scholar, Cryogenics Engineering Laboratory, Department of Mechanical Engineering
National Institute of Technology, Rourkela Odisha, India, 769008
R. K. Sahoo
Professor, Cryogenics Engineering Laboratory, Department of Mechanical Engineering National Institute of Technology, Rourkela Odisha, India
AbstractCurrent paper reports the fluid flow characteristics inside an airfoil convergent nozzle used in a turboexpander. The coordinates of upper (5th order) and lower (3rd order) surface of the nozzle is obtained using curve fitting method in MatlabÂ© after that the computational domain has been generated in SolidWorksÂ©. Numerical analysis has been done to visualize the fluid flow and thermal characteristics inside the nozzle at a pressure ratio of two. Shear stress transport (SST) turbulence model has been used to solve the computational domain. The key feature of this implementation is to obtain subsonic velocity at the nozzle exit. The computing results show that the reduction in temperature for air is approximately 14.2 K with Mach number of 0.88 at the outlet.
KeywordsFluid flow; Convergent nozzle; Ideal gas

INTRODUCTION
Cryogenic fluids like liquid helium, nitrogen, oxygen, hydrogen, etc. are used due to its variety of applications in the fields such as rocket propulsion and aerospace appliances, superconducting equipment, industrial applications, etc. The design of an efficient nozzle plays a vital role in this process. In this framework, researchers are interested to design an efficient and optimized nozzle profile, which is desirable to minimize the losses and compact size of a turboexpander [1, 10].
The effect of geometrical parameters on the performance of the supersonic micronozzle which reduces the thrust losses by increasing the divergent angle above 300 and discussed the thermal, viscous and effect of rarefaction on the performance of nozzle using continuum approach [2]. Researchers suggest that shock wave generation and boundary layer separation in overexpanded supersonic nozzles have unsteadiness during the flow, which causes unsymmetrical flow separation [3]. The experimental analysis of double divergent nozzle to find the shock wave behavior under overexpanded regimes and concluded that the flow field was similar to that of dualbell nozzle [4]. Geron et al. [5] investigated the threedimensional flow field inside a plug nozzle by partitioning it into different modules to enhance the performance of rocket engines. The selection of crosssection at the exit plays a vital role to reduce the thrust loss. Researchers suggest that rectangular cross section possess less thrust loss as compare to circular one whereas density profiles through the circular nozzle have good agreement with the experimental results as compare to that of for the rectangular or square nozzles [6]. The numerical work was carried on to visualize the oscillatory flow pattern and hysteresis phenomenon for different pressure ratios between
two shock structures (FSS and RSS) and its effect on annulus height, wall pressure and shear stress in axisymmetric and truncated contour nozzle [7]. The computational work to visualize the effect of nozzle pressure ratio on flow structure, shockinduced boundary layer separation inside a non axisymmetric supersonic convergentdivergent nozzle was performed by Hasan [8]. Mousavi and Roohi [9] investigated the shock train in a threedimensional convergentdivergent nozzle for compressible and turbulent fluid flow using the Reynolds stress turbulence model (RSM) which is validated with the experimental data.
In this paper, the convergent type nonaxisymmetric airfoil nozzle has been designed which can convert the pressure and thermal energy of the fluid into kinetic energy to achieve subsonic velocity at the outlet, which is essential for an effective design of turboexpander used in gas liquefaction applications. The design method is based on 3rd and 5th order curve fitting under a specified boundary condition which is suitable for the meanline design and provides uniform flow at the inlet of the turbine. The flow behavior of cryogenic fluid air is predicted inside the nozzle using commercially available software ANSYS CFXÂ©. This type of nozzle is essential for mixedflow type having radial inlet and an axial outlet. The rotor is of impulse type with a small amount of reaction. The expansion ratio in the nozzle is approximately two. For exact analysis point of view, experiments provided the exact results, but it is expensive. Accurate prediction of flow field behavior, flow separation, pressure, velocity, Mach number, temperature, TKE, etc. inside the domain is necessary to design a nozzle for high performance of the turboexpander.

MATHEMATICAL MODEL
The physical model considered for the present study is nonaxisymmetric convergent nozzle designed from the curve fitting approach. The oblique shocks occurred inside the nozzle can be minimized by obtaining the continuous contour profile of a nozzle. Due to this, the curvefitting method has been used. The slope of the nozzle wall is continuously decreasing (convergent type) from the inlet to the outlet. After the outlet of the nozzle, it is assumed that the fluid strikes the rotor tangentially. From the manufacturing point of view, the clearance of 1 mm has been taken in the present case.

Turbulence model and discretization scheme
The ShearStress Turbulence (SST) model is widely used due to it possesses combined advantages of k near walls and k in wakes and freeshear regions in the outer boundary layer. It can also control the eddyviscosity by restricting the turbulent shear stress, which improves the model performance in adverse pressure gradients and flows separation cases in particular. The governing equations, mass, momentum, and energy along with pressure based solver transient formulation are adopted in this case. To achieve the accurate simulation results, the selection of the numerical scheme is equally important as the turbulence model. Due to this reason, in the current study, the secondorder upwind scheme is selected for discretization of momentum and energy equations.

Governing equations
Governing equations, which is used to solve the computational domain, are as follows:
Continuity Equation:
. U 0
(a)
t
Momentum Equations:
U
. U U p . SM
t
(1)
(2)
Where the stress tensor is related to the strain rate as:
U U T 2 .U
3
(3)
Total Energy Equation:
htot p .Uh
.T
(b)
t t
tot
(4)
Figure1. (a) Computational mesh used in this study (b) Grid
independence test
.U. U.SM SE
Where htot is the total enthalpy, which is related to static enthalpy:

Boundary conditions
Time integration uses a secondorder implicit scheme for the advancement of the solution with a specified linear solver. The upwind and highresolution methods are selected for a
htot
h 1 U 2
2
(5)
different scheme. The CFD models and their boundary and initial conditions are shown in Table 1. The convergence criteria (RMS) and conservation target are set to be 106 and
0.001 respectively. For dynamic model control, automatc pressure level information, temperature damping, and velocity pressure coupling with RhieChow fourth order model are used.

Grid independence test
To compute the flow field inside the nozzle, a grid independence analysis is done for air with the four different mesh resolutions. Figure 1 represents the velocity at the outlet of the nozzle for four different mesh. Further refinement is not realistic due to very small changes in results but takes more computational times. The computational mesh inside the domain contains 71,972 nodes which are sufficient to predict the flow and thermal characteristics as shown in Figure 1 (a).
P (bar)
T (K)
Mass flow rate (kg/s)
Inlet
16
100
–
Outlet
–
–
0.01
Wall
Noslip
Adiabatic
Table 1. Boundary conditions


RESULTS AND DISCUSSION
The flow and thermal behavior of air, assuming as an ideal fluid, has been predicted inside the nozzle. Generally, ideal nozzles provide uniform and parallel flow at the exit. Figure 2
(a) represents the variation of pressure along the axial distance at different planes for the pressure ratio of two. It clearly shows that the pressure variation inside the domain is from 16 bar (inlet) to 8 bar (outlet).
(a)
(b)
Figure 2. Contour of (a) Pressure (b) Temperature at different locations inside the domain
Figure 2 (b) represents the temperature contours of at different crosssections. It shows that the temperature of the air at the outlet of the nozzle is 86.5 K which shows the decrement of 14.5 K as compared to that of an inlet. This decrease in temperature is very much significant at ultralow temperature because it happens inside the nozzle which reduces the turbine work and hence increases the efficiency of turboexpander. It is desirable for the liquefaction of gases.
Figure 3 represents the velocity and Mach number profile along the axial distance. The velocity at axial distance x = 2.5 mm is approximately 35.20 m/s and continuously increases as one move towards the outlet where it is approximately 166.00 m/s with Mach number of 0.88. It is also observed that the velocity increases rapidly after x = 5 mm and becomes stable at the time of impact with the turbine (outlet of the nozzle) due to which a uniform subsonic flow has occurred. The Mach number at the outlet for helium is lowest for this nozzle.
(a)
(b)
Figure 3. Contour of (a) Velocity (b) Mach number at different locations inside the domain

CONCLUSIONS
The present work reports the numerical investigation to visualize the fluid flow properties such as velocity, Mach number and thermal properties temperature, pressure, etc. inside the nonaxisymmetric convergent nozzle. Shear stress transport (SST) turbulence model is applied in the current simulation. Threedimensional CFD model for turbulent fluid flow inside the convergent nozzle is simulated. The interesting point of this study is the reduction in temperature for air is around 14.2 K with subsonic Mach number of 0.88 at the outlet of the nozzle. The investigation shows the productive results for the design and numerical analysis of nozzles used in turbomachinery applications.
REFERENCES
[1] Kumar, M., Behera, S. K., Sahoo, R. K. Design and Numerical Analysis of Fluid Flow Characteristics in a Nonaxisymmetric Convergent Nozzle. 5th International Conference on Computational Methods for Thermal Problems. 2018 [2] Louisos, William, and Darren Hitt. "Optimal expansion angle for viscous supersonic flow in 2D micronozzles." 35th AIAA Fluid dynamics conference and exhibit. 2005. [3] Hadjadj, A., Y. Perrot, and S. Verma. "Numerical study of shock/boundary layer interaction in supersonic overexpanded nozzles." Aerospace science and technology 42 (2015): 158168. [4] Arora, Rajat, and Aravind Vaidyanathan. "Experimental investigation of flow through planar double divergent nozzles." Acta Astronautica 112 (2015): 200216. [5] Geron, Marco, et al. "Flowfield analysis of a linear clustered plug nozzle with roundtosquare modules." Aerospace Science and Technology 11.23 (2007): 110118. [6] Nguyen, T. V., R. G. Spencer, and A. Siebenhaar. "Aerodynamic performance of a roundtosquare nozzle." (1997). [7] Lijo, Vincent, et al. "Numerical simulation of transient flows in a rocket propulsion nozzle." International Journal of Heat and Fluid Flow 31.3 (2010): 409417. [8] Hasan, ABM Toufique. "Characteristics of overexpanded nozzle flows in imposed oscillating condition." International Journal of Heat and Fluid Flow 46 (2014): 7083. [9] Mousavi, Seyed Mahmood, and Ehsan Roohi. "Three dimensional investigation of the shock train structure in a convergentdivergent nozzle." Acta Astronautica 105.1 (2014): 117127. [10] Kumar, M., Behera, S. K., Sahoo, R. K. Design and numerical investigation to visualize the fluid flow and thermal characteristics of nonaxisymmetric convergent nozzle. Engineering Science and Technology, an International Journal (2018) (Accepted)