• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1966-1971
• /
• 2003

Matching inviscid and boundary layer methods are developed for hypersonic flow with thick boundray layer. The new equations match all the boundary layer properties with a variation in the inviscid solution near the edge, except for the normal velocity. Computational comparison are performed for incompressible and compressible flows over a flat plate. Results from the present method are compared with Navier-Stokes solutions. The present results are in good agreement with Navier-Stokes solutions. They show that the new technique can provide improved heating rates and skin friction predictions for preliminary design of vehicles where shear layers and entropy layer swallowing are important.

• Journal of computational fluids engineering
• /
• v.13 no.2
• /
• pp.27-41
• /
• 2008

An unstructured hybrid mesh flow solver has been developed for the simulation of three-dimensional steady and unsteady incompressible flow fields. The incompressible Navier-Stokes equations with an artificial compressibility method were discretized by using a node-based finite-volume method. For the unsteady time-accurate computation, a dual-time stepping method was adopted to satisfy a divergence-free flow field at each physical time step. An implicit time integration method with local time stepping was implemented to accelerate the convergence in the pseudo-time sub-iteration procedure. The one-equation Spalart-Allmaras turbulence model has been adopted to solve high-Reynolds number flow fields. The flow solver was parallelized to minimize the CPU time and to overcome the computational overhead. This method has been applied to calculate steady and unsteady flow fields around submarine configurations and a 3-D infinite cylinder. Validations were made by comparing the predicted results with those of experiments or other numerical results. It was demonstrated that the present method is efficient and robust for the prediction of steady and unsteady incompressible flow fields.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.1 no.2
• /
• pp.41-47
• /
• 1997

The flow choking in a double orifice is expected to depend on cross-sectional area ratios of the orifices, upstream Mach number and total pressure loss between the orifices. However, no research has been reported on the problems of the compressible flow through a double orifice so far. The present study investigated analytically the choke conditions of the compressible gas through a double orifice, using a simple compressible theory. The orifice area ratio, upstream Mach number, and total pressure loss were involved to find the effects that they have on the flow choking. The results of analytical method show that for orifice area ratios below 1.0, flow choking moves from the first to the second orifice as the total pressure loss increases, however, for orifice area ratios over 1.0, it occurs only at the second orifice.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.37 no.1
• /
• pp.1-16
• /
• 2009

In this paper, we introduce two-phase versions of RoeM and AUSMPW+ schemes. Both schemes are originally developed for the gas dynamic problems, and have shown superior accuracy, efficiency and robustness. A new shock discontinuity sensing term is derived from the mixture equation of state, which is commonly used in the RoeM and AUSMPW+ schemes for the stable numerical flux calculation. The developed two-phase versions of the schemes are applied to several liquid-gas, large property discrepancy two-phase test problems, including several shock stability test problems. The results show that both schemes maintain the merits exhibited in gas dynamic problems even in two-phase flows.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2006.05a
• /
• pp.310-313
• /
• 2006

The present study aims at investigating the effectiveness of a new passive cavity flow control technique, sub-cavity. The characteristics of cavity flow oscillation with the device are compared with those with other control techniques tested previously, including a triangular bump and blowing jet. In the computation, the three-dimensional, unsteady Navier-Stokes equations governing the supersonic cavity flow are solved based on an implicit finite volume scheme spatially and multi-stage Runge-Kutta scheme temporally. Large eddy simulation (LES) is carried out to properly predict the turbulent features of cavity flow. The present results show that the pressure oscillation near the downstream edge dominates overall time-dependent cavity pressure variations, and the amplitude of the pressure oscillation can be reduced in the presence of a sub-cavity.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 1998.04a
• /
• pp.17-17
• /
• 1998

연소반응을 가지는 후류(wake)는 가스터빈 연소실의 flame holder 등에서 발생한다. 후류유동의 안정성 혹은 불안정성은 이러한 유동에 있어서 많은 영향을 끼치므로 상당히 중요하다. 본 연구는 위와 같이 연소 반응에 의한 밀도구배를 가지는 후류유동에 대하여 불안정성 해석을 수행하였다. 교란에 대한 지배방정식은 Navier-Stokes 방정식에서 점성항을 제외한 Euler 방정식을 고려하였다. 충류유동의 압력은 일정하다고 가정하였다. 교란 방정식을 유도하기 위하여 충류 유동이 평행하여 유동 방향에 수직한 방향의 구배만이 존재한다고 가정하였다. 모든 변수들은 충류 유동의 값과 움직이는 파장의 형태를 가지는 작은 교란의 합으로 생각하여 압력에 대한 교란방정식을 구하며, shooting법과 Newton-Raphson법에 근거를 두는 반복법을 사용하여 풀었다. 불안정성 해석을 수행하는 기본 유동의 속도장 및 온도장은 불안정성 해석을 수행하는 기본 유동의 속도장 및 온도장은 비압축성의 경우 우선 Gaussian Profile 가정함과 동시에 연소반응을 포함하는 유동장을 보다 정확히 구하기 위하여 Navier-Stokes 방정식으로부터 구한 결과를 사용하였다. 연소반응을 포함하는 유동장을 구할 때에는 해석상 편의를 위해 예혼합물질은 이상기체로, 화학반응은 1단계의 비가역반응으로 가정하였으며, 반응열로 인한 부력의 효과는 무시하였다. 위와 같은 유동장을 가지고 불안정성 해석을 수행한 결과 후류유동은 두 개의 변곡점을 가지며 sinuous 모드와 varicose모드의 두 개의 불안정성 모드가 있음을 보였다. 밀도 변화가 있는 경우나 밀도 변화가 없는 경우 모두 sinuous 모드의 가장 불안정한 모드가 varicose 모드의 가장 불안정한 모드보다 더 불안정함을 보여주어 후류 유동은 자유 유동에 가까운 위상 속도를 가지는 sinuous 모드에 의해 지배될 것임을 예측할 수 있다. 연소반응이 완전연소에 가까울수록 그리고 압축성 효과가 클수록 유동내부의 온도가 증가하고 점성 또한 증가하여 후류유동은 안정됨을 알 수 있었다 유동변수들의 contour로부터 유동의 특성을 예측한 결과 baroclinic 항이 dilatational 항보다 상대적으로 크며, 중심선 상하에 생기는 vortex를 더욱 성장시킬 것으로 생각된다.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.41 no.10
• /
• pp.649-657
• /
• 2017

A rupture in a high-pressure pipe causes the fluid in the pipe to be discharged in the atmosphere at a high speed resulting in a supersonic jet that generates the compressible flow. This supersonic jet may display complicated and unsteady behavior in general. In this study, Computational Fluid Dynamics (CFD) analysis was performed to investigate the compressible flow generated by a supersonic jet ejected from a high-pressure pipe. A Shear Stress Transport (SST) turbulence model was selected to analyze the unsteady nature of the flow, which depends upon the various gases as well as the diameter of the pipe. In the CFD analysis, the basic boundary conditions were assumed to be as follows: pipe of diameter 10 cm, jet pressure ratio of 5, and an inlet gas temperature of 300 K. During the analysis, the behavior of the shockwave generated by a supersonic jet was observed and it was found that the blast wave was generated indirectly. The pressure wave characteristics of hydrogen gas, which possesses the smallest molecular mass, showed the shortest distance to the safety zone. There were no significant difference observed for nitrogen gas, air, and oxygen gas, which have similar molecular mass. In addition, an increase in the diameter of the pipe resulted in the ejected impact caused by the increased flow rate to become larger and the zone of jet influence to extend further.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2005.11a
• /
• pp.117-120
• /
• 2005

The present study describes unsteady flow phenomena generated in a supersonic flow passing over a rectangular cavity and suggests a way of control of pressure oscillation, doing harm to overall performance and stable operation of aerodynamic and industrial applications. The three-dimensional, unsteady, compressible Navier-stokes equations are numerically solved based on a fully implicit finite volume scheme and large eddy simulation. The cavity flow are simulated with and without control methods, including a triangular bump and blowing jet installed near the leading edge of the cavity. The results show that the pressure oscillation is attenuated by both control techniques, especially near the trailing edge of cavity.

• Proceedings of the KSME Conference
• /
• 2005.11a
• /
• pp.42.2-42.2
• /
• 2005

• Proceedings of the KSME Conference
• /
• 2005.11a
• /
• pp.82.1-82.1
• /
• 2005