• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.112-122
• /
• 1996

A finite element scheme using the concept of PISO method has been developed to solve the Navier-Stokes viscous flows in all speed range. This scheme includes development of new pressure equation that retains both the hyperbolic term related with the density variation and the elliptic term reflecting the incompressibility constraint. The present method is applied to the incompressible two-dimensional driven cavity flow problems(Re=100, 400 and 1,000). For compressible flows, the Carter plate problem(M=3 and Re=1,000) is computed. Finally, we have simulated the shock-boundary layer interaction(M=2 and Re=2.96×10/sup 5/), a more difficult problem, and compared its results with the experiment to demonstrate the shock capturing capability of the present solution algorithm.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.41 no.9
• /
• pp.615-622
• /
• 2017

A normal shock wave is initially formed in the shock tube that migrates towards the closed end of the tube, which, in turn, leads to the reflection of shock. Due to the interaction of the reflected shock with the boundary layer, bifurcation of the shock wave takes place. A shock train will be generated after the bifurcated shock wave approaches the contact surface. Until now, only a few studies have been conducted to investigate this shock train phenomenon inside the shock tube. For the present study, a CFD analysis has been performed on a two dimensional axisymmetric model of a shock tube using unsteady, compressible Navier-Stokes equations. In order to investigate the detailed characteristics of the shock train phenomenon, quantitative studies have been performed by varying shock tube length, diameter under fixed diaphragm, and pressure ratio inside a shock tube.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.40 no.5
• /
• pp.375-384
• /
• 2012

This paper presents numerical results comparing the performance of the 2008 Wilcox $\mathcal{k}-{\omega}$ turbulence model to the one of the 1988 Wilcox $\mathcal{k}-{\omega}$ model for supersonic flows. A comparison with experimental data is offered for a shock wave/turbulent boundary layer interaction case and two ramp injector mixing cases. Furthermore, a comparison is performed with empirical correlations on the basis of skin friction for flow over a flat plate and shear layer growth for a free shear layer. It is found that the maximum injectant mass fraction of some ramp injector cases is better predicted using the 1988 Wilcox model. On the other hand, the 2008 model performs better in simulating shock-boundary layer cases.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.34 no.12
• /
• pp.18-24
• /
• 2006

Passive control of the shock wave/turbulent boundary-layer interaction utilizing slotted plates and a porous plate over a cavity has been carried out. Effect of various slot configurations on the characteristics of the interaction has been observed. Pitot/wall surface pressure distributions and flow visualizations including Schlieren images, kerosene-lampblack tracings and interference fringe patterns over a thin oil-film have been obtained at the downstream of the shock interactions. For the streamwise-slot configuration, a local higher pitot pressure was noticed at the downstream of the interaction as compared with the case of no control, however, not much improvement in pitot pressure was observed for the spanwise-slot configuration.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.75-78
• /
• 2004

Various jet engines (Turbine engine family and RAM Jet engine) have been developed for high speed aircrafts. but their application to hypersonic flight is restricted by principle problems such as increase of total pressure loss and thermal stress. Therefore, the development of next generation propulsion system for hypersonic aircraft is a very important subject in the aerospace engineering field, SCRAM Jet engine based on a key technology, Supersonic Combustion. is supposed as the best choice for the hypersonic flight. Since Supersonic Combustion requires both rapid ignition and stable flame holding within supersonic air stream, much attention have to be given on the mixing state between air stream and fuel flow. However. the wider diffusion of fuel is expected with less total pressure loss in the supersonic air stream. So. in this study the direction of fuel injection is inclined 30 degree to downstream and the total pressure of jet is controlled for lower penetration height than thickness of boundary layer. Under these flow configuration both streams, fuel and supersonic air stream, would not mix enough. To spread fuel wider into supersonic air an aerodynamic force, baroclinic torque, is adopted. Baroclinic torque is generated by a spatial misalignment between pressure gradient (shock wave plane) and density gradient (mixing layer). A wedge is installed in downstream of injector orifice to induce an oblique shock. The schlieren optical visualization from side transparent wall and the total pressure measurement at exit cross section of combustor estimate how mixing is enhanced by the incidence of shock wave into supersonic boundary layer composed by fuel and air. In this study non-combustionable helium gas is injected with total pressure 0.66㎫ instead of flammable fuel to clarify mixing process. Mach number 1.8. total pressure O.5㎫, total temperature 288K are set up for supersonic air stream.

• Journal of computational fluids engineering
• /
• v.20 no.4
• /
• pp.36-43
• /
• 2015

Accurate prediction of supersonic transition is required for the heat transfer estimation over supersonic double compression ramp flows. Correlation-based transition models were assessed for a supersonic double ramp problem. Numerical results were compared with experimental data from RWTH Aachen University. A parametric study on the nose bluntness was performed using a selected transition model. As the nose bluntness increases, the boundary layer thickness is increased and the triple point of shock interactions moves downstream. The peak magnitude of the heat transfer is consequently decreased with the nose bluntness.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.117-122
• /
• 1998

A numerical study is carried out to investigate the transient process of combustion phenomena associated with hypersonic propulsion devices. Reynolds averaged Navier-Stokes equations for reactive flows are used as governing equations with a detailed chemistry mechanism of hydrogen-air mixture and two-equation SST turbulence modeling. The governing equations are discretized by a high order accurate upwind scheme and solved in a fully coupled manner with a fully implicit time accurate method. At first, oscillating shock-induced combustion is analyzed and the comparison with experimental result gives the validity of present computational modeling. Secondly, the model ram accelerator experiment was simulated and the results show the detailed transient combustion mechanisms. Thirdly, the evolution of oblique detonation wave is simulated and the result shows transient and final steady state behavior at off-stability condition. Finally, shock wave/boundary layer interaction in combustible mixture is studied and the criterion of boundary layer flame and oblique detonation wave is identified.

• Journal of the Korean Society of Visualization
• /
• v.13 no.3
• /
• pp.29-34
• /
• 2015

The interaction between a normal shock wave and a boundary layer along a wall surface in internal compressible flows causes a very complicated flow. This interaction region containing shock train and mixing region is called as pseudo-shock waves. Pseudo-shock waves in the divergent part of a rectangular nozzle have been investigated by using large-eddy simulation (LES). LES studies have been done for the complex flow phenomena of three-dimensional pseudo-shock waves. The LES results have been validated against experimental wall-pressure measurements. The LES results are in good agreement with experimental results. Pseudo-shock length and corner separation have been studied in three-dimensional LES model. Comparison of centerline pressure measurement and 3D visualization measurement has been discussed for the corner separation position. It has been concluded that the pseudo-shock length should be measured by using 3D visualization measurement.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.6 no.4
• /
• pp.124-131
• /
• 2003

When a secondary gaseous flow is injected vertically into a supersonic flow through circular nozzle, a complicated structure of flow field is produced around the injection area. The interaction between the two streams produces a strong bow shock wane on the upstream side of the side-jet. The results show that bow shock wave and turbulent boundary layer interaction induces the boundary layer separation in front of the side-jet. This study is to analyze the structure of flow fields and distribution of surface pressure on the flat plate according to total pressure ratio using a supersonic cold-flow system and also to study the control force of affected side-jet. The nozzle of main flow was designed to have Mach 2.88 at the exit. The injector has a sonic nozzle with 4mm diameter at the exit of the side-jet. In experiments, The oil flow visualization using a silicone oil and ink was conducted in order to analyze the structure of flow fields around the side-jet. The flow fields are visualized using the schlieren method. In this study, a computational fluid dynamic solution is also compared with experimental results.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.88-97
• /
• 1996

This paper treats an adaptive finite-element method for the viscous compressible flow governed by Navier-Stokes equations in two dimensions. The numerical algorithm is the two-step Taylor-Galerkin mettled using unstructured triangular grids. To increase accuracy and stability, combined moving node method and grid refinement method have been used for grid adaption. Validation of the present algorithm has been made by comparing the present computational results with the existing experimental data and other numerical solutions. Four benchmark problems are solved for demonstration of the present numerical approach. They include a subsonic flow over a flat plate, the Carter flat plate problem, a laminar shock-boundary layer interaction. and finally a laminar flow around NACA0012 airfoil at zero angle of attack and free stream Mach number of 0.85. The results indicates that the present adaptive triangular grid method is accurate and useful for laminar viscous flow calculations.