• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2011.11a
• /
• pp.178-181
• /
• 2011

Recently, micro shock tube is being extensively used in various fields of engineering applications. The flow characteristics occurring in the micro shock tube may be significantly different from that of conventional macro shock tube due to very low Reynolds number and Knudsen number effects which are, in general, manifested in such flows of rarefied gas, solid-gas two-phase, etc. In these situations, Navier-Stokes equations cannot properly predict the micro shock tube flow. In the present study, a two-dimensional CFD method has been applied to simulate the micro shock tube, with slip velocity and temperature jump boundary conditions. The effects of wall thermal conditions on the unsteady flow in the micro shock tube were also investigated. The unsteady behaviors of shock wave and contact discontinuity were, in detail, analyzed. The results obtained show much more attenuation of shock wave, compared with macro-shock tubes.

• Journal of Powder Materials
• /
• v.2 no.2
• /
• pp.120-134
• /
• 1995

The behavior of the flow about gas atomizers with a supersonic nozzle containing an under-expanded or over-expanded jet is very important with respect to performance and stability characteristics. Since detailed experiments are expensive, computational fluid mechanics have been applied recently to various relating flow field. In this study, a higher order upwind method with the 3rd order MUSCL type TVD scheme is used to solve the full Reynolds Wavier-Stokes equations. To delineate the purely exhaust jet effects, the melt flow is not considered. Comparison is made with some experimental data in terms of density fields. The influence of the exhaust-jet-to freestream pressure ratio and the effect of the protrusion length of the melt orifice are studied. The present study leads us to believe that the computational fluid mechanics should be considered as powerful tool in predicting the gas atomizer flows.

• Journal of computational fluids engineering
• /
• v.5 no.1
• /
• pp.43-59
• /
• 2000

Characteristic boundary conditions are discussed in conjunction with a flux-difference splitting formulation as modified from Roe's linearization. Details of how one can implement the characteristic boundary conditions which are made compatible with the interior point formulation are described for different types of boundaries including subsonic outflow and adiabatic wall. The validity of boundary conditions are demonstrated through computation of transonic airfoil, supersonic ogive-cylinder, hypersonic cylinder, and S-duct internal flows. The computed wall pressure distributions are compared with published experimental and computed data. Objectives of this paper are thus to give insight of formulation procedure of a flux-difference splitting method and to pave ways for other users to adopt present boundary procedure on their numerical methods.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.117-122
• /
• 1998

A numerical study is carried out to investigate the transient process of combustion phenomena associated with hypersonic propulsion devices. Reynolds averaged Navier-Stokes equations for reactive flows are used as governing equations with a detailed chemistry mechanism of hydrogen-air mixture and two-equation SST turbulence modeling. The governing equations are discretized by a high order accurate upwind scheme and solved in a fully coupled manner with a fully implicit time accurate method. At first, oscillating shock-induced combustion is analyzed and the comparison with experimental result gives the validity of present computational modeling. Secondly, the model ram accelerator experiment was simulated and the results show the detailed transient combustion mechanisms. Thirdly, the evolution of oblique detonation wave is simulated and the result shows transient and final steady state behavior at off-stability condition. Finally, shock wave/boundary layer interaction in combustible mixture is studied and the criterion of boundary layer flame and oblique detonation wave is identified.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.04a
• /
• pp.30-40
• /
• 2007

The implicit discontinuous Galerkin method for the two-dimensional Euler equations was developed on unstructured triangular meshes, which can achieve higher-order accuracy by wing hierachical basis functions based on Legendre polynomials. Numerical tests were conducted to estimate the convergence order of numerical solutions to the Ringleb flow and the supersonic vortex flow for which analytic solutions are available. And, the flows around a circle and a NACA0012 airfoil was also numerically simulated. Numerical results show that the implicit discontinuous Galerkin methods with higher-order representation of curved solid boundaries can be an efficient higher-order method to obtain very accurate numerical solutions on unstructured meshes.

• 유체기계공업학회:학술대회논문집
• /
• 2005.12a
• /
• pp.481-487
• /
• 2005

A turbopump turbine consists of rotational part including a rotor and stationary part including nozzles and exit guide vanes, of which shape and relative position affect turbine performance owing to supersonic flows with prevailing unsteadiness. In this study, numerical 3-D flow calculations of the turbine with the different number of exit guide vanes and different relative position of each component are conducted and the effect of flow passage variations on turbine performance is analyzed.

• Journal of Mechanical Science and Technology
• /
• v.15 no.10
• /
• pp.1423-1433
• /
• 2001

A flow modeling method has been developed to analyze the flow in the annular base (rear- facing surface) of a circular engine nacelle flying at subsonic speed but with a supersonic exhaust jet. Real values of exhaust gas properties and temperature at an altitude of 30, 000 feet are employed. Potential flows of the air and gas streams are computed for the flow past a separated wake. Then a viscous jet mixing is superimposed on this inviscid solution. Conserva- tion of mass, momentum and energy for the wake flow field is achieved by multiple iterations with modest computer requirements.

• 한국전산유체공학회:학술대회논문집
• /
• 2002.05a
• /
• pp.29-34
• /
• 2002

TIn this study, the counter-jet flows which designed for improvement of aerodynamic performance of the blunt body vehicle have been analyzed. The variations of the drag force and jet penetration depth due to changes in the stagnation properties of counter jet new such as total pressure, mach number, and total temperature. The counter jet flow, which is injected toward incoming supersonic freestream at stagnation region of blunt cone-cylinder vehicle, have been studied by using upwind flux difference splitting navier-stokes method. The changes in the stagnation pressure and Mach number resulted in large effects on the wall pressure and drag force, on the other hand tile total temperature changes did not.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.3068-3073
• /
• 2007

Aerodynamic forces and moments have been used to control rocket propelled vehicles. If control is required at very low speed, Those systems only provide a limited capability because aerodynamic control force is proportional to the air density and low dynamic pressure. But thrust vector control(TVC) can overcome the disadvantages. TVC is the method which generates the side force and roll moment by controlling exhausted gas directly in a rocket nozzle. TVC is classified by mechanical and fluid dynamic methods. Mechanical methods can change the flow direction by several objects installed in a rocket nozzle exhaust such as tapered ramp tabs and jet vane. Fluid dynamic methods control the flight direction with the injection of secondary gaseous flows into the rocket nozzle. The tapered ramp tabs of mechanical methods are used in this paper. They installed at the rear in the rocket nozzle could be freely moved along axial and radial direction on the mounting ring to provide the mass flow rate which is injected from the rocket nozzle. In this paper, the conceptual design and the performance study on the tapered ramp tabs of the thurst vector control has been carried out using the supersonic cold flow system and shadow graph. Numerical simulation was also performed to study flow characteristics and interactions between ramp tabs. This paper provides to analyze the location of normal shock wave and distribution of surface pressure on the region enclosed by the tapered ramp tabs.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.2
• /
• pp.316-328
• /
• 1997

There were appreciable progresses on the study of shock wave / boundary layer interaction control in the transonic flow without nonequilibrium condensation. But in general, the actual flows associated with those of the airfoil of high speed flight body, the cascade of steam turbine and so on accompany the nonequilibrium condensation, and under a certain circumstance condensation shock wave occurs. Condensation shock wave / boundary layer interaction control is quite different from that of case without condensation, because the droplets generated by the result of nonequilibrium condensation may clog the holes of the porous wall for passive control and the flow interaction mechanism between the droplets and the porous system is concerned in the flow with nonequilibrium condensation. In these connections, it is necessary to study the condensation shock wave / boundary layer interaction control by passive cavity in the flow accompanying nonequilibrium condensation with condensation shock wave. In the present study, experiments were made on a roof mounted half circular arc in an indraft type supersonic wind tunnel to evaluate the effects of the porosity, the porous wall area and the depth of cavity on the pressure distribution around condensation shock wave. It was found that the porosity of 12% which was larger than the case of without nonequilibrium condensation produced the largest reduction of pressure fluctuations in the vicinity of condensation shock wave. The results also showed that wider porous area, deeper cavity for the same porosity of 12% are more favourable "passive" effect than the cases of its opposite. opposite.