• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.110-119
• /
• 2006

An accurate and cost efficient method PSE is used for the stability analysis of 2D or 3D compressible boundary layers. A highly accurate finite difference PSE code has been developed at a general curvilinear coordinate system using an implicit marching procedure to deal with a broad range of transition predictions problems. Evolution of disturbances in compressible flat plate boundary layers are studied for free-stream Mach numbers ranging from 0 to 1.5. The effect of mean-flow nonparallelism is found to be weak on two dimensional waves and strong on three dimensional waves. The maximum amplification rate increases monotonically with Mach number. The present PSE solutions are compared with previous numerical investigations and experimental results and are found to be in good agreement.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.737-742
• /
• 2008

Methods for the solution of linear stability for compressible boundary layers are developed. Both the global and local methods for stability analysis are used. Both methods are use in solution of Coutte shear flow and the results are analysis and compare. Some well-known conclusions of Coutte flow are proved by these methods again.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.4
• /
• pp.509-517
• /
• 1997

Results of flat plate compressible boundary layer calculation, based on discrete formulation of DSMC method, are presented in low Mach number and low Knudsen number range. The free stream is a uniform flow of pure nitrogen at various Mach numbers in low pressures (i.e. rarefied gas). Complete thermal accommodation and diffuse molecular reflections are used as the wall boundary condition, replacing unreal no-slip condition used in continuum calculations. In the discrete formulation of DSMC method, there is no need to use ad hoc assumptions on transport properties like viscosity and thermal conductivity, instead viscosity is calculated from values of other field variables (velocity and shear stress). Also the results are compared with existing self-similar continuum solutions. In all Mach number cases computed, velocity slip is most pronounced in regions near the leading edge where continuum formulation renders the solution singular. As the boundary layer develops further downstream, velocity slips asymptote to values that are between 10 to 20% of the magnitude of free stream velocity. When the free stream number density is reduced, so the gas more rarefied, the velocity slip increases as expected.

• 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 the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.5
• /
• pp.407-419
• /
• 2008

In this study, the transition Reynolds number of compressible axi-symmetric sharp cone boundary layer is predicted by using a linear stability theory and the -method. The compressible linear stability equation for sharp cone boundary layer was derived from the governing equations on the body-intrinsic axi-symmetric coordinate system. The numerical analysis code for the stability equation was developed based on a second-order accurate finite-difference method. Stability characteristics and amplification rate of two-dimensional second mode disturbance for the sharp cone boundary layer were calculated from the analysis code and the numerical code was validated by comparing the results with experimental data. Transition prediction was performed by application of the -method with N=10. From comparison with wind tunnel experiments and flight tests data, capability of the transition prediction of this study is confirmed for the sharp cone boundary layers which have an edge Mach number between 4 and 8. In addition, effect of wall cooling on the stability of disturbance in the boundary layer and transition position is investigated.

• 유체기계공업학회:학술대회논문집
• /
• 2006.08a
• /
• pp.235-238
• /
• 2006

This paper is about the code that realizes the $e^N$-Method for boundary-layer transition prediction. The $e^N$-Method based on the linear stability theory is applied to predicting boundary-layer transition frequently. This paper deals with the construction of code, stability analysis and the calculation of N-factor. The results of transition prediction using the $e^N$-Method for flat plate/cone compressible boundary-layers are presented.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.37 no.11
• /
• pp.1055-1065
• /
• 2009

In this study, stability analysis of boundary layers on airfoils is performed by using parabolized stability equations(PSE). Boundary layer edge conditions are obtained by compressible inviscid flow calculations. Mean velocity and temperature profiles of the laminar boundary layer are obtained by solving compressible boundary layer equations in generalized curvilinear coordinates with fourth order accuracy in the wall normal direction. Laminar mean flow profiles are used as input data for PSE to investigate growth rates of disturbances and stability characteristics. For the cases of boundary layer on NACA0012 and HSNLF(1)-0213 airfoils at Mach number 0.5, growth rates with respect to disturbance frequencies and profiles of disturbance amplitude are investigated. The effect of angle of attack on stability characteristics are examined at both upper and lower surfaces. The neutral stability curves, effect of Mach number and effect of airfoil section shapes are also analyzed.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2000.06a
• /
• pp.104-109
• /
• 2000

The modes of oscillation and radiated acoustic fields of compressible flows over open cavities are investigated computationally. The compressible Navier-Stokes equations are solved for two-dimensional cavities with laminar boundary layers upstream. The high-order and high-resolution numerical schemes are used for the evaluation of spatial derivatives and the time integration. Physically correct numerical boundary conditions are implemented to produce time-accurate solutions in the whole computation domain. The computational domain is large enough to directly resolve a portion of the radiated acoutic field. The results show a transition from a shear layer mode, for shorter cavities and lower Mach numbers, to a wake mode for longer cavities and higher Mach numbers. The shear layer mode is well characterized by Rossiter modes and these oscillations lead to intense upstream acoustic radiation dominated by a single frequency. The wake mode is characterized instead by a large-scale vortex shedding. Acoustic radiation is more intense, with multiple frequencies present.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.90-95
• /
• 1998

The two-dimensional incompressible and compressible Navier-Stokes codes are developed for the computation of the viscous turbulent flow over high-lift airfoils. Incompressible code using pseudo-compressibility and dual-time stepping method involves a conventional upwind differencing scheme for the convective terms and LU-SGS scheme for time integration. Compressible code also adopts an FDS scheme and LU-SGS scheme. Several two-equation turbulence models (the standard $k-{\varepsilon}$ model, the $k-{\omega}$ model. and $k-{\omega}$ SST model) are evaluated by computing the flow over single and multi-element airfoils. The compressible and incompressible codes are validated by computing the flow around the transonic RAE2822 airfoil and the NACA4412 airfoil, respectively. Both the results show a good agreement with experimental surface pressure coefficients and velocity profiles in the boundary layers. Also, the GA(W)-1 single airfoil and the NLR7301 airfoil with a flap are computed using the two-equation turbulence models. The grid systems around two- and three-element airfoil are efficiently generated using Chimera grid scheme, one of the overlapping grid generation methods.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.7 no.3 s.18
• /
• pp.165-173
• /
• 2004

Matching inviscid and boundary layer methods are developed for analysis of hypersonic flow with thick boundary 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 peformed for incompressible and compressible flows over a flat plate. Results from the present method are compared with Wavier-Stokes solutions. The present results are in good agreement with Wavier-Stokes solutions. They show that the new technique can provide improved predictions of heating rates and skin friction predictions for preliminary design of vehicles where shear layers and entropy layer swallowing are importantfor for preliminary design.