• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2000.11a
• /
• pp.20-20
• /
• 2000

2개의 경사 충격파와 하나의 수직 충격파로 초음속 유동을 압축하는 초음속 공기 흡입구의 수치적 연구를 수행하였다. 지배방정식으로는 Navier-Stokes방정식을 사용하였고 난류모델로는 SST 모델을 사용하였다. 지배방정식의 점성항 계산에는 중심차분법을 사용하였고 대류항 계산에는 풍상차분법인 Roe의 FDS기법을 MUSCL기법과 결합하여 이용하였다. 유한 체적법을 이용하여 차분된 방정식은 LU분할 기법을 이용한 완전 내재적 방법으로 2차 정확도 시간 적분으로 비정상 과정의 연구를 수행하였다. 흡입구 배압을 정해주어야 하는 어려움을 해결하기 위해 흡입구 후면에 노즐을 달고 노즐의 면적을 조절하여 배압이 형성되도록 하였다.(중략)

• Proceeding of EDISON Challenge
• /
• 2013.04a
• /
• pp.337-342
• /
• 2013

이번 연구에서는 de Laval nozzle를 이용하여 다양한 유체해석 모델과 프로그램을 비교하여 그 성능을 파악하였다. de Laval nozzle은 eigenvalue에 의해 eigenvector값이 '-'와 '+'값을 동시에 갖는 물리현상을 내포하고 있으며, 압력조건에 따라 내부에서 Normal shock이 발생하게 된다. 이러한 non-linearity를 현재 우리가 주로 사용하고 있는 상용프로그램(cfx, fluent)과 EDISON, 직접 코딩한 프로그램(Matlab이용)이 얼마나 잘 표현하는지 알아보았다. 그 결과 Van Leer Vector Splitting을 이용할 경우 물리현상을 제일 잘 표현 하였다. 또한 난류 유동(Turbulence flow)을 고려하게 될 경우, Mesh가 Boundary layer를 표현할 정도로 정밀하지 못하다면 제대로 된 해석 결과를 얻을 수 없었으며, Wall 근처에서 Non-slip condition에 의해 Vortex가 형성되고, 이 Vortex가 Back flow를 유도하여 해가 수렴하는데 방해를 하게 됨을 알 수 있었다. 이를 방지하기 위해서는 유동이 잘 표현될 수 있도록 적절한 Computational environment를 형성해 주는 것이 매우 중요하다.

• Journal of computational fluids engineering
• /
• v.7 no.4
• /
• pp.35-41
• /
• 2002

The projectile afterbodies for zero-lift drag reduction has been analyzed using the Navier-Stokes equations with the κ-εturbidence model. The numerical method of a second order upwind scheme has been used on an unstructured adaptive grid system. Base drag reduction methods that have been found effective on axisymmetric bodies are boattailing, base bleed, base combustion, locked vortex afterbodies and multistep afterbodies. In this paper, turbulence flow and pressure charateristics have been studied for geometries of multistep afterbodies. The important geometrical and flow parameters relevant to the design of such afterbodies have been identified by step number, length and height. The flow over multistep aftoerbodies or base have many kinds of compressible flow characteristics including expansion waves at the trailing edge, recompression waves, separation and recirculating flow in the base region, shear flow and wake flow. The numerical results have been compared and analyzed with the experimental data. The flow characteristics have been clearly shown.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.17 no.4
• /
• pp.43-48
• /
• 2013

A study is analyzed on the adverse pressure gradient and the transient regime of supersonic diffuser with Computational Fluid Dynamic. The flow field of supersonic diffuser is calculated using Axisymmetric two-dimensional Navier-Stokes equation with $k-{\epsilon}$ turbulence model. The transient simulation is compared in terms of mach number and static temperature of vacuum chamber according to pressure variation of rocket engine combustion chamber. Combustion gas flow into the vacuum chamber during operation of the supersonic diffuser. According to this phenomenon, the pressure and the temperature rise in the vacuum chamber were observed. Thus, the protection system will be necessary to prevent the pressure and temperature rise in the transition process during operation of the subsonic diffuser.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.2 no.1
• /
• pp.1-12
• /
• 1998

Industrial ejector system is a facility to transport, to compress or to pump out a low pressure secondary flow by using a high pressure primary flow. An advantage of the ejector system is in its geometrical simplicity, not having any moving part, compared with other fluid machinery. Most of the previous works have been performed experimentally and analytically. The obtained data. are too insufficient to improve our current understanding on the detailed flow field inside the ejector. In order to provide more comprehensive data on this ejector flow field, two-dimensional computations using Reynolds-averaged Navier-Stokes equations were performed for a very wide range of operating pressure ratio of the supersonic ejector with a secondary throat. The current results showed that the supersonic ejector system has an optimum pressure ratio for the secondary flow total pressure to be minimized. The numerical results clearly revealed the shock system, shock/boundary layer interaction, and secondary flow entrainment inside the supersonic ejector.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.17 no.1
• /
• pp.61-69
• /
• 2013

Numerical simulations have been performed for assessment of compressibility correction of two-equation turbulent models suitable for large scale separation flows perturbed by a pintle strokes. Two-equation turbulence models, the low Reynolds k-${\varepsilon}$ and the k-${\omega}$ SST models with or without compressibility correction proposed by Wilcox and Sarkar are evaluated. The detail flow structures are observed and static pressures along nozzle wall are compared with experimental results. Mach disk location and pressure recovery profiles in flow separation region are noticeably distinct between turbulent models of k-${\varepsilon}$ and k-${\omega}$ SST. The compressible effect corrections to those models improve resolving of separation flow behaviors. The compressibility corrections to k-${\varepsilon}$ model have provided very comparable results with test data.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.37 no.10
• /
• pp.955-966
• /
• 2009

DES method is applied to an axisymmetric base flow at supersonic mainstream. The model is based on the Spalart-Allmaras (S-A) turbulence model in the RANS mode, and is based on the subgrid scale model in the Large-eddy simulation (LES) mode. Accurate predictions of the base flowfield and base pressure are successfully achieved by using the DES methodology which is less expensive than LES. Flow properties at the edge of base, such as boundary layer thickness, momentum thickness and skin fraction are compared with Dutton et al [experimental data to proper prediction of base flowfiled. From the present results, The DES accurately resolves the physics of unsteady turbulent motions, such as shear layer rollup, large-eddy motions in the downstream region and small eddy motions inside the recirculating region. Moreover, The present results of using an empirical constant $C_{DES}$ of 1.2 shows good agreement with experimental data than conventional empirical constant $C_{DES}$ of 0.65.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.35 no.2
• /
• pp.94-100
• /
• 2007

An axisymmetric supersonic jet screech in the Mach number range from 1.07 to 1.2 is numerically simulated. The axisymmetric mode is the dominant screech mode for an axisymmetric jet. The Reynolds-averaged Navier-Stokes equations in the conjunction with a modified Spalart-Allmaras turbulence model are employed. A high resolution finite volume essentially non-oscillatory(ENO) schemes are used along with nonreflecting characteristic boundary conditions that are crucial to screech tone computations to accurately capture the sound waves, shock-cell structures and large-scale instability waves.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2007.04a
• /
• pp.329-334
• /
• 2007

An axisymmetric supersonic jet screech in the Mach number range from 1.07 to 1.2 is numerically simulated. The axisymmetric mode is the dominant screech mode for an axisymmetric jet. The Reynolds-averaged Navier-Stokes equations in the conjunction with modified Spalart-Allmaras turbulence model are employed. A high resolution finite volume essentially non-oscillatory(ENO) schemes are used along with nonreflecting characteristic boundary conditions that are crucial to screech tone computations to accurately capture the sound waves, shock-cell structures, unsteady shock motions and large-scale instability waves.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• v.y2005m4
• /
• pp.399-402
• /
• 2005

The present numerical study describes the flow physics on the interaction between the supersonic freestream and jet plume. The compressible flow past a simplified afterbody model with a sonic nozzle is investigated using mass-averaged Navier-Stokes equations, discretized by a fully implicit finite volume scheme, and the standard $k-{\omega}$ turbulence model. The results obtained through the present study are discussed specifically regarding the effect of the plume pressure ratio, freestream Mach number and base dimensions on the location of the plume-induced shock wave generated on the afterbody by the underexpansion of the jet plume.