• 한국항공우주학회지
• /
• 제40권9호
• /
• pp.745-751
• /
• 2012

하이브리드 로켓은 축 방향의 산화제 유동과 고체 연료의 기화로 인한 벽면 분출 유동사이의 상호 간섭에 의해 복잡한 형태의 혼합 전단층이 존재한다. 연소실 입구에 링이 설치되어 있고 질량분사가 있는 실린더 유동에 대하여 압축성 효과를 고려한 LES(Large Eddy Simulation) 기법을 적용하여 수치계산을 수행하였다. 최근의 실험에 의하면 연료 중간에 링과 같은 다이아프램이 설치된 경우, 연소율의 증가가 관찰되었다. 계산 결과에 따르면, 축방향 유동과 벽면 분출 유동이 상호 간섭하여 발생하는 벽면 와류가 국부적인 연료 표면으로의 열전달을 증가시켜 실험에서 관찰되는 딤플이 생성되는 것을 확인하였다. 또한 연소실 입구에 설치된 링에 의하여 발생되는 와류는 벽면 와류가 보다 활발하게 생성되고 열전달과 혼합을 향상시키는 역할을 하며 이 때문에 연소율이 증가되는 것으로 보인다.

• 대한기계학회논문집B
• /
• 제27권10호
• /
• pp.1464-1471
• /
• 2003

In order to investigate the flows with shock wave in branch, 108$^{\circ}$ elbow and T-junction of the IRWST system of standard Korean nuclear reactor, detail time dependent behaviors of unsteady flow with shock wave, vortex and so on are obtained by numerical method using compressible three-dimensional Navier-Stokes equations. At first, the complex flow including the incident and reflected shock waves, vortex and expansion waves which are generated at the corner of T-junction is calculated by the commercial code of FLUENT6 and is compared with the experimental result to obtain the validation of numerical method. Then the flow fields in above mentioned units are analyzed by numerical method of [mite volume method. In numerical analysis, the distributions of flow properties with the moving of shock wave and the forces acting on the wall of each unit which can be used to calculate the size of supporting structure in future are calculated specially. It is found that the initial shock wave of normal type is re-established its type from an oblique one having the same strength of the initial shock wave at the 4 times hydraulic diameters of downstream from the branch point of each unit. Finally, it is turned out that the maximum force acting on the pipe wall becomes in order of the T-junction, 108$^{\circ}$ elbow and branch in magnitude, respectively.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2003년도 The Fifth Asian Computational Fluid Dynamics Conference
• /
• pp.294-298
• /
• 2003

The complex, unsteady, self-sustained pressure oscillations induced by supersonic flow past a rectangular cavity is investigated using numerical simulations. The present numerical study is performed using a parallel, multiblock solver for the two-dimensional, compressible Navier­Stokes equations. Open cavities with length-to-depth (L / D) ratio in the range 0.5 - 3.3 are considered. This paper sheds light on the cavity physics, cavity oscillatory mechanism, and the organisation of vortical structures inside the cavity. The vortex shedding phenomenon, the shear layer impingement event at the aft wall and the movement of the acoustic/compression wave within the cavity are well predicted. The vortical structures· and the source of the acoustic disturbances are found to be located near the aft wall of the cavity. With the increase in the cavity length, strong recompression of the flow near the aft wall leading to a sudden jump in the cavity form drag is observed. The estimated cavity tones are in good agreement with the available semi­empirical relation. Multiple peaks are noticed in deep and long cavities. For the present free­stream Mach number 1.71, it is observed that around L/D=2.0, the cavity oscillatory mechanism changes from the transverse to longitudinal oscillatory mode. The effects of this transition on various fluid dynamics and acoustic properties are also discussed.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• pp.151-154
• /
• 2002

The twin impulse wave leads to very complicated flow fields, such as Mach stem, spherical waves, and vortex ring. The twin impulse wave discharged from the exits of the two tubes placed in parallel is investigated to understand detailed flow physics associated with the twin impulse wave, compared with those in a single impulse wave. In the current study, the merging phenomena and propagation characteristics of the impulse waves are investigated using a shock tube experiment and by numerical computations. The Harten-Yee's total variation diminishing (TVD) scheme is used to solve the unsteady, two-dimensional, compressible, Euler equations. The Mach number $M_{s}$, of incident shock wave is changed below 1.5 and the distance between two-parallel tubes, L/d, is changed from 1.2 to 4.0. In the shock tube experiment, the twin impulse waves are visualized by a Schlieren optical system for the purpose of validation of computational work. The results obtained show that on the symmetric axis between two parallel tubes, the peak pressure produced by the twin-impulse waves and its location strongly depend upon the distance between two parallel tubes, L/d and the incident shock Mach number, $M_{s}$. The predicted Schlieren images represent the measured twin-impulse wave with a good accuracy.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2004년도 춘계학술대회
• /
• pp.1643-1648
• /
• 2004

When a shock wave propagates into a Helmholtz resonator, very complicated wave phenomena are formed both inside and outside the resonator tube. Shock wave reflection, shock focusing phenomena and shock-vortex interactions cause strong pressure fluctuations inside the resonator, consequently leading to powerful sound emission. In the present study, the wave phenomena inside and outside the Helmholtz resonator are, in detail, investigated with a help of CFD. The Mach number of the incident shock wave is varied below 2.0 and several types of resonators are tested to investigate the influence of resonator geometry on the wave phenomena. A TVD scheme is employed to solve the axisymmetric, compressible, Euler equations. The results obtained show that the configuration of the Helmholtz resonator significantly affects the peak pressure of shock wave focusing, its location, the amplitude of the discharged wave and resonance frequency.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2009년도 학술발표회 초록집
• /
• pp.516-523
• /
• 2009

The artificial compressibility (AC) method for the incompressible Navier-Stokes equations in the generalized curvilinear coordinates using the primitive form is implemented. The main advantage of the AC approach is that the resulting system of equations resembles the system of compressible N-S equations and can thus be integrated in time using standard, well-established time-marching methods. The errors, which are the odd-even oscillation, for pressure field in using the artificial compressibility can be eliminated by using the $4^{th}$ order artificial dissipation term which is explicitly included. Even though this paper focuses exclusively on 2D laminar flows to validate and assess the performance of this solver, this numerical method is general enough so that it can be readily extended to carry out 3D URANS simulation of engineering flows. This algorithm yields practically identical velocity profiles and primary vortex and secondary vortices that are in excellent overall agreement with the results of the vorticity-stream function formulation (Ghia et al., 1982). However, the grid resolution have to be required to be large enough to express the various vortices.

• 한국추진공학회지
• /
• 제14권3호
• /
• pp.9-15
• /
• 2010

과팽창 로켓노즐에서 발생하는 충격파 박리패턴의 천이 유동장을 예측하기 위해 축대칭 수치해석적연구를 수행하였다. 비정상, 압축성 N-S 방정식에 k-$\omega$ SST 난류모델을 적용하여, 유한 체적법으로 계산하였다. 종래의 실험적 연구 결과와 비교하였으며, 계산된 결과와 정성적으로 잘 일치하였다. 본 연구의 결과로부터 RSS에서 FSS로 천이할 때 가장 큰 횡력이 발생하며, 이는 비점성 제트 코어 영역에서 발생하는 Vortex ring의 발달로 기인됨을 예측하였다. 또한 엔진 시동과정과 정지과정에서 발생하는 히스테리시스 현상을 잘 모사하였다.

• 한국해양공학회지
• /
• 제18권5호
• /
• pp.15-21
• /
• 2004

An edge tone is the discrete tone or narrow-band sound produced by an oscillating free shear layer, impinging on a rigid surface. In this paper, we present a 2-D edge tone to predict the frequency characteristics of the discrete oscillations of a jet-edge feedback cycle, using the finite difference lattice Boltzmann method (FDLBM). We use a modified version of the lattice BGK compressible fluid model, adding an additional term and allowing for longer time increments, compared to a conventional FDLBM, and also use a boundary fitted coordinates system. The jet is chosen long enough in order to guarantee the parabolic velocity profile of the jet at the outlet, and the edge consists of a wedge with an angle of ${\alpha}$ = 23. At a stand-off distance, the edge is inserted along the centerline of the jet, and a sinuous instability wave, with real frequency, is assumed to be created in the vicinity of the nozzle and propagates towards the downstream. We have succeeded in capturing very small pressure fluctuations, resulting from periodical oscillations of a jet around the edge. The pressure fluctuations propagate with the speed of sound. Its interaction with the wedge produces an non-rotational feedback field, which, near the nozzle exit, is a periodic transverse flow, producing the singularities at the nozzle lips.

• 한국전산유체공학회지
• /
• 제18권1호
• /
• pp.1-6
• /
• 2013

Missile and fighter aircraft have been challenged by low restoring nose-down pitching moment at high angle of attach. The consequence of weak nose-down pitching moment can be resulting in a deep stall condition. Especially, the pressure oscillation has a huge effect on noise generation, structure damage, aerodynamic performance and safety, because the flow has strong unsteadiness at high angle of attack. In this paper, the unsteady aerodynamics coefficients were analyzed at high angle of attack up to 50 degrees around two dimensional NACA0012 airfoil. The two dimensional unsteady compressible Navier-Stokes equation with a LES turbulent model was calculated by OHOC (Optimized High-Order Compact) scheme. The flow conditions are Mach number of 0.3 and Reynolds number of $10^5$. The lift, drag, pressure, entropy distribution, etc. are analyzed according to the angle of attack. The results of average lift coefficients are compared with other results according to the angle of attack. From a certain high angle of attack, the strong vortex formed by the leading edge are flowing downstream as like Karman vortex around a circular cylinder. The primary and secondary oscillating frequencies are analyzed by the effects of these unsteady aerodynamic characteristics.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2003년도 추계학술대회
• /
• pp.502-507
• /
• 2003

Spiral jet is characterized by a wide region of the free vortex flow with a steep axial velocity gradient, while swirl jet is largely governed by the forced vortex flow and has a very low axial velocity at the jet axis. However, detailed generation mechanism of spiral flow components is not well understood, although the spiral jet is extensively applied in a variety of industrial field. In general, it is known that spiral jet is generated by the radial flow injection through an annular slit which is installed at the inlet of convergent nozzle. The objective of the present study is to understand the flow characteristics of the spiral jet, using a computational method. A finite volume scheme is used to solve 3-dimensional Navier-Stokes equations with RNG ${\kappa}-{\varepsilon}$ turbulent model. The computational results are validated by the previous experimental data. It is found that the spiral jet is generated by coanda effect at the inlet of the convergent nozzle and its fundamental features are dependent the pressure ratio of the radial flow through the annular slit and the coanda wall curvature.