• 유체기계공업학회:학술대회논문집
• /
• 2002.12a
• /
• pp.337-345
• /
• 2002

The present study addresses a computational result of unsteady gas flow through a critical nozzle. The axisymmetric, unsteady, compressible, Wavier-Stokes equations are solved using a finite volume method that makes use of the second order upwind scheme for spatial derivatives and the multi-stage Runge-Kutta integral scheme for time derivatives. The steady solutions of the governing equation system are validated with the previous experimental data to ensure that the present computational method is valid to predict the critical nozzle flows. In order to simulate the effects of back pressure fluctuations on the critical nozzle flows, an excited pressure oscillation with an amplitude and frequency is assumed downstream of the exit of the critical nozzle. The results obtained show that for low Reynolds numbers, the unsteady effects of the pressure fluctuations can propagate upstream of the throat of critical nozzle, and thus giving rise to the applicable fluctuations in mass flow rate through the critical nozzle, while for high Reynolds numbers, the pressure signals occurring at the exit of the critical nozzle do not propagate upstream beyond the nozzle throat. For very low Reynolds number, it is found that the sonic line near the throat of the critical nozzle remarkably fluctuateswith time, providing an important mechanism for pressure signals to propagate upstream of the nozzle throat, even in choked flow conditions. The present study is the first investigation to clarify the unsteady effects on the critical nozzle flows.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.1
• /
• pp.54-60
• /
• 2003

Compound choking frequently occurs at a minimum area of the flow passage, where two or more streams which have different stagnation properties are merged. This phenomenon is especially important in that the flow choking may not be given by Mach number, M=1 at the nozzle throat. In order to obtain a detailed understanding of the flow characteristics involved in the compound flow choking, the two-dimensional, compressible, Wavier-Stokes equations are solved using a fully implicit finite volume method and the predicted solutions are compared with the results of the one-dimensional theoretical analysis. Stagnation pressure and temperature of each stream are changed to investigate the effects on the compound choking. The results show that stagnation pressures of each stream affect Mach number and static pressure distributions downstream of the exit of the convergent nozzle. However, the flow characteristics of the compound choking are not significantly dependent on the total temperature ratio.

• Journal of Energy Engineering
• /
• v.12 no.1
• /
• pp.29-34
• /
• 2003

Thermo-acoustic waves can be generated in a compressible fluid by rapid heating and cooling near the boundary walls. These phenomena are very important mechanism of heat transfer in the space environment in which natural convection does not exist. In this study, the generation and transmission characteristics of thermo-acoustic waves in an air filled enclosure with rapid wall heating are studied numerically. The governing equations were discretized using control volume method, and were solved using PISO algorithm and second-order upwind scheme. For the stable solution time step were considered as t=1$\times$$10^{-9}$ order, and grids are 50$\times$800. The induced thermo-acoustic wave propagates through the fluid until it decays due to viscous and heat dissipation. The wave showed sharp front shape and decreased with long tail.

• Journal of the Korean Society of Visualization
• /
• v.8 no.3
• /
• pp.56-62
• /
• 2010

The interaction patterns between shock wave and boundary layer in a rocket nozzle are mainly classified into two categories, FSS(Free Shock Separation) and RSS(Restricted Shock Separation), both of which are associated with the thrust characteristics as well as side loads of the engine. According to the previous investigations, strong side loads of the engine are produced during the period of transition from FSS to RSS or vice versa. The present work aims at investigating the unsteady behavior of the separation shock waves in a two-dimensional supersonic nozzle, using experimental method and CFD. Schlieren optical method was employed to visualize the time-mean and time-dependent shock motions in the nozzle. The unsteady, compressible N-S equations with SST K-$\omega$ turbulence closure were solved using a fully implicit finite volume scheme. The results obtained show the separation shock motions during the transition of the interaction pattern.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2005.11a
• /
• pp.117-120
• /
• 2005

The present study describes unsteady flow phenomena generated in a supersonic flow passing over a rectangular cavity and suggests a way of control of pressure oscillation, doing harm to overall performance and stable operation of aerodynamic and industrial applications. The three-dimensional, unsteady, compressible Navier-stokes equations are numerically solved based on a fully implicit finite volume scheme and large eddy simulation. The cavity flow are simulated with and without control methods, including a triangular bump and blowing jet installed near the leading edge of the cavity. The results show that the pressure oscillation is attenuated by both control techniques, especially near the trailing edge of cavity.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.14 no.6
• /
• pp.25-30
• /
• 2010

Dual throat nozzle(DTN) is recently attracting much attention as a new concept of the thrust vectoring technique. This DTN is designed with two throats, an upstream minimum and a downstream minimum at the nozzle exit, with a cavity in between the upstream throat and exit. In the present study, a computational work has been carried out to analyze the performance of a dual throat nozzle(DTN) at various mass flow rate of secondary flow and nozzle pressure ratios(NPR). Two-dimensional, steady, compressible Navier-Stokes equations were solved using a fully implicit finite volume scheme. The present computational results were validated with some experimental data available. Based upon the present results, The control effectiveness of thrust-vector is discussed in terms of the thrust coefficient and the discharge coefficient.

• 유체기계공업학회:학술대회논문집
• /
• 2006.08a
• /
• pp.227-230
• /
• 2006

The method of mass flow rate measurement using a critical nozzle is well established in the flow satisfying ideal gas law. However, in the case of measuring high-pressure gas flow, the current method shows invalid discharge coefficient because the flow does not follow ideal gas law. Therefore an appropriate equation of state considering real gas effects should be applied into the method. The present computational study has been performed to give an understanding of the physics of a critical nozzle flow for high-pressure hydrogen gas and find a way for the exact mass flow prediction. The two-dimensional, axisymmetric, compressible Navier-Stokes equations are computed using a fully implicit finite volume method. The real gas effects are considered in the calculation of discharge coefficient as well as in the computation. The computational results are compared with the previous experimental data and predict well the measured mass flow rates. It has been found that the discharge coefficient for high-pressure hydrogen gas can be corrected properly adopting the real gas effects.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.10
• /
• pp.1240-1247
• /
• 1999

Particle motion through a disk type critical orifice placed in a 3.0cm diameter chamber has been studied numerically. In the simulation, the velocity field is solved using Pantankar's SIMPLER algorithm for the compressible flow and convergence of the computation is confirmed if the mass source at each control volume is smaller than $10^{-7}$. The particle motion in the flow field is solved in Lagrangian method. The particle trajectories showed that the particles injected away from the center line are expanded rapidly. At lower pressures, this expansion phenomena are more dominant. At lower pressures, the clear difference in particle and air speed is showed all the way down to the exit plan. It was found that particles with Stokes number of ca.2.5 tend to focus close to the center line very well except the particles travelling near the wall. However, particles with Stokes number greater than ca.2.5 show a tendency to cross the center line.

• Journal of computational fluids engineering
• /
• v.15 no.4
• /
• pp.25-31
• /
• 2010

Interface tracking of two phase is significant to analyze multi-phase phenomena. The VOF(Volume of Fluid) and level set are well known interface tracking method. However, they have limitations to solve compressible flow and incompressible flow at the same time. CIP(Cubic Interpolate Propagation) method is appropriate for considering compressible and incompressible flow at once by solving the governing equation which is divided up into advection and non-advection term. In this article, we analyze the droplet impingement according to various We number using improved CIP method which treats nonlinear term once more comparison with original CIP method. Furthermore, we compare spread radius after droplet impingement on the wall with the experimental data and original CIP method. The result using improved CIP method shows the better result of the experiments, comparison with result of original CIP method, and it reduces the mass conservation error which is generated in the numerical analysis comparison with original CIP method.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.14 no.5
• /
• pp.1-7
• /
• 2010

Dual throat nozzle(DTN) is recently attracting much attention as a new concept of the thrust vectoring technique. This DTN is designed with two throats, an upstream minimum and a downstream minimum at the nozzle exit, with a cavity in between the upstream throat and exit. In the present study, a computational work has been carried out to analyze a fundamental performance of a dual throat nozzle(DTN) at various nozzle pressure ratios(NPR) and throat area ratios. Two-dimensional, axisymmetric, steady, compressible Navier-Stokes equations were solved using a fully implicit finite volume scheme. NPR was varied in the range of NPR from 2.0 to 10.0, at different throat area ratios. The present computational results were validated with some experimental data available. Based upon the present results, the performance of DTN is discussed in terms of the discharge coefficient and thrust efficiency.