• 유체기계공업학회:학술대회논문집
• /
• 유체기계공업학회 1999년도 유체기계 연구개발 발표회 논문집
• /
• pp.179-186
• /
• 1999

Three-dimensional flow analysis is implemented to investigate the flow through transonic axial-flow compressor rotor(NASA R67), and to evaluate the performances of k-$\epsilon$ and Baldwin-Lomax turbulence models. A finite volume method is used for spatial discretization. And, the equations are solved implicitly in time with the use of approximate factorization. Upwind difference scheme is used for inviscid terms, but viscous terms are centrally differenced. The flux-difference-splitting of Roe is used to obtain fluxes at the cell faces. Numerical analysis is performed near peak efficiency and near stall. And, the results are compared with the experimental data for NASA R67 rotor. Blade-to-Blade Mach number distributions are compared to confirm the accuracy of the code. From the results, we conclude that k-$\epsilon$ model is better for the calculation of flow rate and efficiency than Baldwin-Lomax model. But, the predictions for Mach number and shock structure are almost same.

• 유체기계공업학회:학술대회논문집
• /
• 유체기계공업학회 1999년도 유체기계 연구개발 발표회 논문집
• /
• pp.27-32
• /
• 1999

A numerical study based on the three-dimensional thin-layer Navier-Stokes solver is carried out to analyze the flowfield through a single stage transonic compressor. Explicit four-step Runge-Kutta scheme with spatially variable time step and implicit residual smoothing is used. The governing equations are discretized with exploit finite difference method. Mixed-out average method is used at the interface between rotor and stator. And, an artificial dissipation model is used to assure the stability of solution. The results with k-$\omega$ turbulence model were compared to the results with Baldwin-Lomax model, and physical phenomena of transonic compressor are presented. The two turbulence models give the results that show reasonably good agreements with experimental data.

• 대한기계학회논문집B
• /
• 제38권3호
• /
• pp.219-225
• /
• 2014

본 연구에서는 NACA0012/14/15 천음속 에어포일 유동에서 비평형 응축이 충격파 진동에 미치는 영향을 TVD 수치해석을 통하여 연구하였다. 주류 마하수 0.81-0.90에 대해, 정체점 상대습도 및 에어포일의 기하학적 형상이 유동 특성에 미치는 영향이 구명되었다. 받음각 ${\alpha}=0^{\circ}$, 정체점 온도(288K) 및 주류 마하수가 0.87인 경우, 정체점 상대습도의 증가는 Terminating Shock의 충격파 강도를 약화시킨다. 정체점 상대습도가 30%인 경우 961Hz이던 충격파의 진동수가 60%일 때는 912Hz로 약 5% 감소한다. 정체점 상대습도가 동일한 경우는 주류 마하수가 클수록 충격파의 진동수 및 이동거리는 크게 된다. 또, 진동의 한 주기에 대해 항력계수의 변화도 구명되었다. 정체점 상대습도가 높을수록 최대 항력 계수는 작고, 항력계수의 변화폭 또한 감소한다. 한편 에어포일의 최대 두께가 두꺼울수록 초음속 영역의 크기는 증가하며 충격파의 진동수 및 이동거리도 증가한다.

• 한국항공우주학회지
• /
• 제36권3호
• /
• pp.211-219
• /
• 2008

본 연구에서는 삼차원 비정렬 비점성 유동 해석 코드를 이용하여 전투기 형상의 외부 장착물이 꼬리 날개에 미치는 공력 간섭효과에 대한 연구를 수행하였다. 수치적 기법으로는 격자점 중심(vertex-centered)에 기초한 유한체적법과 시간 적분을 위한 내재적인 point Gauss-Seidel 반복 계산법을 사용하였다. 해석 코드의 검증을 위해 전투기 형상에 대한 정상 유동 계산을 수행하였으며 결과를 실험 데이터와 비교하였다. 외부 장착물의 후류(wake)를 정확히 포착하기 위해 예상되는 후류 영역에 대해 국부적인 격자 조밀화를 수행하였으며 천음속 영역에서의 비정상 유동 해석을 통해 외부 장착물에서 발생하는 후류가 수평꼬리날개에 미치는 간섭효과를 확인하였다. 공력 간섭효과를 감소시키기 위한 대안으로는 뒷전 플랩 꺽임각(trailing edge flap deflection)을 고려하였으며 이에 대한 정량적인 감소효과를 제시하였다.

• 한국유체기계학회 논문집
• /
• 제15권2호
• /
• pp.12-19
• /
• 2012

This study presents unsteady and unstable characteristics of three stage transonic axial compressor, developed by Korea Aerospace Research Institute. As approaching to the unstable operating region at the 103% design speed of the compressor, a modal type stall precursor appears in front of highly loaded 3rd rotor row at first, and it propagates to the upstream. On the contrary, actual stall cell initiates from the stall precursor in front of the 1st rotor row, and it propagates to the downstream of the compressor. After the stall region reached the 3rd stage and stall cell rotates circumferentially about 360 deg, it develops to one dimensional compressor surge mode. It shows a mild surge behaviour with 3~4 Hz frequency. From the test data, it can be suggested that there is a priority to give an optimum blade loading distributions to construct a multi stage transonic axial compressor stages either to secure more stable compressor operating ranges, or to maximize the compressor efficiency.

• 한국항공우주학회지
• /
• 제39권10호
• /
• pp.912-918
• /
• 2011

본 연구의 목적은 풍동실험을 통해 회전발사체에서 발생하는 동적 롤댐핑 특성을 실험적으로 측정하는데 있었으며, 이를 위해 약 12,000 rpm으로 회전하는 회전발사체 실험모형에 작용하는 롤댐핑 특성의 측정을 위한 고속풍동실험을 국방과학연구소 삼중음속풍동에서 수행하였다. 실험시의 마하수는 0.7~1.05까지의 천음속 영역이었으며 이때의 받음각 구간은 -4~+10 deg이었다. 풍동실험 측정기법의 유효성 평가를 위해 동일형상 모형에 대해 기 수행하였던 롤댐핑 측정결과와의 비교검토를 수행하였다.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2008년도 추계학술대회논문집
• /
• pp.167-174
• /
• 2008

In this study, flutter analyses for supercritical airfoil have been conducted in transonic region. Advanced computational analysis system based on computational fluid dynamics (CFD) and computational structural dynamics (CSD) has been developed in order to investigate detailed static and dynamic responses of supercritical airfoil. Reynolds-averaged Navier-Stokes equations with Spalart-Allmaras (S-A) and SST ${\kappa}-{\omega}$ turbulence models are solved for unsteady flow problems. A fully implicit time marching scheme based on the Newmark direct integration method is used for computing the coupled aeroelastic governing equations of cascades for fluid-structure interaction (FSI) problems. Also, flow-induced vibration (FIV) analyses for various supercritical airfoil models have been conducted. Detailed flutter responses for supercritical are presented to show the physical performance and vibration characteristics in various angle of attack.

• 한국전산유체공학회지
• /
• 제20권1호
• /
• pp.47-56
• /
• 2015

In this paper, we study the effect of various turbulence models by comparing the aerodynamic characteristics and the flow patterns computed for aircraft models. An in-house CFD solver, MSAPv, that solves the three dimensional RANS equations with the turbulence model equations is used. The turbulence models used in this study are the Spalart-Allmaras model, Menter's $k-{\omega}$ SST model, Coakley's $q-{\omega}$ model, and Huang and Coakley's $k-{\varepsilon}$ model. DLR-F6 WB and WBNP configurations are selected for the study. We concentrate on the separated flow pattern variations with the turbulence models at the wing-body junction and the wing-pylon junction as well as drag polar curves.

• 대한기계학회논문집B
• /
• 제26권7호
• /
• pp.997-1004
• /
• 2002

The current study describes experimental and computational works on the passive control of the steady and unsteady condensation shock waves, which are generated in a transonic nozzle. The bleed slots are installed on the contoured wall of the transonic nozzle in order 10 control the magnitude of the condensation shock wave and its oscillations. For computations, a droplet growth equation is copuled with two-dimensional Navier-Stokes equation systems. Computations are carried out using a third-order MUSCL type TVD finite-difference scheme with a second-order fractional time step. Baldwin-Lomax turbulence model is employed to close the governing equations. An experiment using an indrafi transonic wind tunnel is made to validate the computational results. The current computations represented well the experimental flows. From both the experimental and computational results it is found that the magnitude of the condensation shock wave in the bleed slotted nozzle is signi ficantly reduced, compared with no passive control of solid wall. The oscillations of the condensation shock wave are successfully suppressed by a bleed slot system.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 춘계 학술대회논문집
• /
• pp.168-173
• /
• 1998

A numerical analysis of axisymetric backward facing step main nozzle flow in air jet loom has been accomplished. To obtain basic design data for an optimum main nozzle for an air-jet loom and to predict the transonic/supersonic flow, a characteristic based upwind flux difference splitting compressible Navier-Stokes method has been used. The wall static pressure of the main nozzle and the flow velocity changes in the nozzle tube were analyzed by changing air tank pressures and acceleration tube lengths. The flow inside the nozzle experiences double choking one at the needle tip and the other at the acceleration tube exit at tank pressures over $4kg_f/cm^2$. The tank pressure $P_t$ leading to the critical condition depends on the acceleration tube length; i.e, $P_t$ is higher for longer acceleration tubes. The $P_t$ value required to bring the acceleration tube exit to the critical condition is nearly constant regardless of acceleration tube length. The round needle tip shape might lead to less total pressure loss when compared with step shape.