• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.36 no.12
• /
• pp.1209-1215
• /
• 2012

Experiments with premixed flames in a tube have been conducted to investigate the transition phenomenon from a laminar flat flame to turbulent motions. To induce this phenomenon, a flat flame is formed in a tube. Then, the local velocity at the center of the flat flame surface is increased using $CO_2$ laser irradiation. The deformed flame front propagates with an increase in the total flame surface and oscillating instability. Eventually, the flame front accelerates explosively, and it shows turbulent flame motions with a strong noise. The dynamic behaviors of the flame front prior to the turbulent motions are analyzed in this study to elucidate this process. The physical model of the process is presented according to observations.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.11a
• /
• pp.19-23
• /
• 2010

Injection characteristics of a kerosene swirl injector of liquid rocket engine operating at supercritical environment have been investigated. Kerosene surrogate models are proposed to model the kerosene properties. Turbulent numerical model is based on large eddy simulation and contains Soave modification of Redlich-Kwong equation of state and Chung's model. Numerical analysis results at supercritical environment are compared with the one at transcritical condition. Differences of density and viscosity are analyzed at both liquid film and core gas in the swirl injector.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.6 no.2
• /
• pp.75-84
• /
• 2002

A numerical analysis of unsteady motion in solid rocket motors with a nozzle has been conducted. The numerical formulation including modified $\kappa$-$\varepsilon$ turbulence model treats the complete conservation equation for the gas phase and the one-dimensional equations in the radial direction for the condensed phase. A fully coupled implicit scheme based on a dual time-stepping integration algorithm has been adopted to solve the governing equations. After obtaining a steady state solution, pulse and periodic oscillations of pressure are imposed at the head-end to simulate acoustic oscillations of a travelling-wave motion in the combustion chamber. Various steady and unsteady state features in the combustion chamber of a rocket motor has been analyzed as results of numerical calculations.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.40 no.6
• /
• pp.475-484
• /
• 2012

Numerical investigation is performed to understand the effects of thickness uncertainty of a supporting airfoil due to manufacturing processes on the aerodynamic characteristics of an airfoil used for measuring data in a wind tunnel testing. This is done by comparing the coefficients of lift, drag and moment of the airfoils. In this work, the airfoil model consists of three parts, one located in the center for measuring and two outer parts used for supporting. The study is carried out with a NACA64-418 airfoil and the turbulence model of Transition SST. It is found that the effect of thickness uncertainty of the airfoils used for supporting is not significant to the performance of the test airfoil at various angles of attack and Reynolds numbers.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.48 no.5
• /
• pp.323-334
• /
• 2020

Since boundary layer transition has a significant impact on the aero-thermodynamic performance of hypersonic flight vehicles, capability of accurate prediction of transition location is essential for design and performance analysis. In this study, γ-Re_θt model is improved to predict transition of hypersonic boundary layers and validated. A coefficient in the production term of the intermittency transport equation that affects the transition onset location is constructed and applied as a function of Mach number, wall temperature, and freestream stagnation temperature based on the similarity numerical solution of compressible boundary layer. To take into account a Mach number dependency of transition onset momentum thickness Reynolds number and transition length, additional correlation equations are determined as function of Mach number and applied to Re_θc and F_length correlations of the baseline model. The suggested model is implemented to a commercial CFD code in consideration of practical use. Analysis of hypersonic flat plate and circular cone boundary layers is carried out by using the model for validation purpose. An improvement of prediction capability with respect to variation of Mach number and unit Reynolds number is identified from the comparison with experimental data.

• Tunnel and Underground Space
• /
• v.11 no.4
• /
• pp.319-327
• /
• 2001

Introduction of new design tools has been required to optimally design and operate the ventilation system of long vehicle tunnels.. The demand has led to wide spread use of the simulation technique throughout the would to analysis the dynamic relationship among the variables associated with vehicle tunnel ventilation. This paper aims at performing on-site study at local tunnels to test the applicability of NETVEN, a simulation model vehicle tunnel ventilation. The study was carried out at four urban as well as highway tunnels model of vehicle tunnel ventilation. The study was carried out at four urban as well as highway tunnels employing different ventilation systems as well as traffic methods. There were some discrepancies sound between the simulation output and measurements and the following four factors are considered to mainly cause those disagreement. (1) The real situation shows distinctive transient and retarding characteristics with respect to air flow and contaminant dispersion, while ventilation forces are not steady-state and in particular those traffic and climatic variables show significant instantaneous variation. (3) Near the exit portal, the CO levels show bigger differences. The general trend is that data with higher CO concentrations carry bigger discrepancies. Turbulent diffusion is though to be the main reason for it and also contribute to the fact hat the highest CO concentrations are found at the locations somewhat inward, not at the exit portals. (4) Higher traffic rate results in higher discrepancies of ventilation velocity. Along with the exhaust characteristics, the vehicle aerodynamic characteristics need to be studied continuously in order to reduce the velocity disagreement.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.10a
• /
• pp.15-24
• /
• 2000

The transient response of an airfoil to a rapidly deploying spoiler is numerically investigated using the turbulent compressible Navier-Stokes equations in two dimensions. Algebraic Baldwin-Lomax model, Wilcox $\kappa-\omega$ model, and SST $\kappa-\omega$ turbulence model are used to calculate the unsteady separated flow due to the rapid spoiler deployment. The spoiler motion relative to a stationary airfoil is treated by an overset grid hounded by a Dynamic Domain-Dividing Line which has been devised by the authors. The adverse effects of the spoiler influenced by the spoiler location and the hinge gap are expounded. The numerical results are in reasonably good agreement with the existing experimental data.

• Transactions of The Korea Fluid Power Systems Society
• /
• v.8 no.2
• /
• pp.36-41
• /
• 2011

Analysis of two-dimensional flow in an involute gear pump has been done by using FLUENT. Analysis extended to the turbulent flow includes the mass flow rate with functions of pressure difference between inlet and outlet, rotational velocities of involute gear, and clearances between tip of gear and housing. In general mass flow rate decreases with decreasing rotational velocity, and with increasing clearance and pressure difference. The flow rate efficiency of gear pump, which is defined with the theoretical flow rate, has been presented in terms of the above parameters.

• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.75-75
• /
• 2011

본 논문에서는 진동중인 교량 거더에 작용하는 풍하중을 산정하고 그에 따른 플러터 발생풍속을 예측하기 위하여 분산형 전산환경을 활용한 수치해석 연구를 수행하였다. 분산형 전산환경은 웹 포탈을 기반으로 수치해석 환경을 제공하는 수치풍동 시스템으로서, 전산유체역학(CFD : Computational Fluid Dynamics)에 대한 전문지식이 부족한 사용자들도 격자생성, 수치해석자를 이용한 계산, 가시화 등의 전 과정을 편리하게 수행할 수 있는 차세대 토목분야 연구 환경이다. 본 시스템은 그리드스피어(GfidSphere)를 기반으로 구성되었으며, 기본적으로 사용자 관리, 세션 관리, 그룹 관리, 레이아웃 관리 등을 제공하여 사용자가 포탈을 통해서 다양한 서비스를 쉽게 사용할 수 있는 환경을 구축하도록 도와준다. 수치해석을 위한 유체 지배방정식은 2차원 비정상 비압축성 RANS(Reynolds-Averaged Navier-Stokes) 방정식이며, pseudo compressibility 방법을 적용하였다. 비정상 유동장을 해석하기 위하여 이중시간 전진법(dual time stepping)을 사용하였으며, 수렴가속화를 위해 Multi-grid 기법을 적용하였다. 또한 난류 유동장 해석을 위해서 $k-{\omega}$ SST 난류 모델을 사용하였으며, 난류 천이 과정에서의 유동을 모사하기 위하여 Total stress limitation 방법을 적용하였다. 교량 거더의 연직과 회전방향의 2자유도 움직임을 모사하기 위하여 동적격자 기법을 도입하였다. 교량 거더 주변의 비정상 유동해석 결과를 통해, 거더 표면에서 떨어져나가는 크고 작은 와류의 영향으로 양력 및 모멘트 계수 그래프가 중첩된 진폭과 주기를 갖고 주기적으로 나타나는 것을 확인할 수 있었다. 또한 계산된 비정상 공기력을 적용한 2자유도 플러터 방정식을 통하여 플러터 발생풍속을 산정하였다. 최종적으로 본 연구에서 계산된 결과의 타당성을 검증하기 위하여 수치적으로 구한 플러터 발생풍속과 기존의 실험 및 수치해석 결과를 비교하였으며, 결과는 잘 일치하였다.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2017.05a
• /
• pp.679-681
• /
• 2017

Numerical analysis was carried out by applying the concept of E-D(expansion-deflection) nozzle to dual bell nozzle. We used the CEA code to calculate the chemical composition of the nozzle and to analyze the freezing flow of 8 species. The turbulence model was chosen as the $k-{\omega}$ SST. We applied the concept of E-D nozzle to the dual bell nozzle and performed the calculated transition altitude and performance. As a result of the interpretation, the application of the E-D nozzle concept led to the formation of over-expansion conditions, which resulted in an increase in the transition altitude.