• 대한기계학회논문집
• /
• 제18권7호
• /
• pp.1840-1850
• /
• 1994

Enhancement of numerical algorithms for low speed compressible flow will be considered. Contemporary time-marching algorithm has been widely accepted and applied as the method of choice for transonic, supersonic and hypersonic flows. In the low Mach number regime, time-marching algorithms do not fare as well. When the velocity is small, eigenvalues of the system of compressible equations differ widely so that the system becomes very stiff and the convergence becomes very slow. This characteristic can lead to difficulties in computations of many practical engineering problems. In the present approach, the time-derivative preconditioning method will be used to control the eigenvalue stiffness and to extend computational capabilities over a wide range of flow conditions (from very low Mach number to supersonic flow). Computational capabilities of the above algorithm will be demonstrated through computation of a variety of practical engineering problems.

• 한국항공우주학회지
• /
• 제37권1호
• /
• pp.17-27
• /
• 2009

압축성 영역에서 개발된 2상 유동 RoeM과 AUSMPW+를 전 마하수 영역 해석을 위해 예조건화 하였다. 압축성 영역에서 개발된 풍상차분 수치기법들은 수치점성에서 음속을 포함하는 항들을 갖고 있는데, 이는 비압축성 영역에서 관련 수치점성 크기를 지나치게 크거나 작게 만들어 해의 정확도를 떨어뜨리거나 수치 불안정 현상을 일으킨다. 따라서 예조건화된 시스템 고유값을 수치점성에 고려하여 전 마하수 영역에서 해의 정확성을 확보하였고, 점근 분석을 통해 개발된 수치기법이 비압축 극한 영역에서 보존형 지배방정식과 일관됨을 확인하였다. 수치 예제들을 통해 예조건화된 2상 유동 RoeM과 AUSMPW+가전 마하수 영역에서 합리적인 해를 도출함을 알 수 있었다.

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

In this study, 2D computations of the Aeolian tones for some obstacles (circular cylinder, square cylinder and NACA0012 airfoil) are simulated. First of all, we calculate the flow noise generated by a uniform flow around a two-dimensional circular cylinder at Re=150 are simulated by applying the finite difference lattice Boltzmann method (FDLBM). The third-order-accurate up-wind scheme (UTOPIA) is used for the spatial derivatives, and the second-order-accurate Runge-Kutta scheme is applied for the time marching. The results show that we successively capture very small acoustic pressure fluctuation with the same frequency of the Karman vortex street compared with the pressure fluctuation around a circular cylinder. The propagation velocity of the acoustic waves shows that the points of peak pressure are biased upstream due to the Doppler effect in the uniform flow. For the downstream, on the other hand, it is faster. To investigate the effect of the lattice dependence, furthermore, simulations of the Aeolian tones at the low Reynolds number radiated by a square cylinder and a NACA0012 airfoil with a blunt trailing edge at high incidence are also investigated.

• 한국전산유체공학회지
• /
• 제17권3호
• /
• pp.53-58
• /
• 2012

Steady flow simulations through a high pressure turbine guide vanes were carried out. The main objective of the present work is to study the performance of turbulence models on the steady flow prediction from aerodynamic and aerothermal points of view. Three turbulence models were compared, namely SST, k-${\omega}$ and ${\omega}$-Reynolds stress models. The laminar results were also compared. The comparison was done with emphasis on the isentropic Mach number and heat transfer coefficient along the blade, and total pressure loss in the wake region. The calculated isentropic Mach number showed reasonable agreement with experimental data along the blade surface for all three turbulent models. For the total pressure loss in the wake region, ${\omega}$-Reynolds stress model showed the best agreement with the experimental data. However, unless using an appropriate transition model, the heat transfer coefficients of all three turbulent models showed poor agreement with experimental data.

• 대한기계학회논문집B
• /
• 제21권10호
• /
• pp.1273-1283
• /
• 1997

Experiments of shock wave/turbulent boundary layer interaction were conducted by using a supersonic wind tunnel. Nominal Mach number was varied in the range of 1.6 to 3.0 by means of different nozzles. The objective of the present study is to investigate the effects of boundary layer flow bleed on the interaction flow field in a straight tube. Two-dimensional slits were installed on the tube walls to bleed the turbulent boundary layer flows. The bleed flows were measured by an orifice. The ratio of the bleed mass flow to main mass flow was controlled within the range of 11 per cent. The wall pressures were measured by the flush mounted transducers and Schlieren optical observations were made for almost all of the experiments. The results show that the boundary layer flow bleed reduces the multiple shock waves to a strong normal shock wave. For the design Mach number of 1.6, it was found that the normal shock wave at the position of the silt was resulted from the main flow choking due to the suction of the boundary layer flow.

• 대한기계학회논문집B
• /
• 제21권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.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2008년 영문 학술대회
• /
• pp.671-679
• /
• 2008

The effect of Pre-cooled Turbojet Engine installation and nozzle exhaust jet on Hypersonic Turbojet EXperimental aircraft(HYTEX aircraft) were investigated by three-dimensional numerical analyses to obtain aerodynamic characteristics of the aircraft during its in-flight condition. First, simulations of wind tunnel experiment using small scale model of the aircraft with and without the rectangular duct reproducing engine was performed at M=5.1 condition in order to validate the calculation code. Here, good agreements with experimental data were obtained regarding centerline wall pressures on the aircraft and aerodynamic coefficients of forces and moments acting on the aircraft. Next, full scale integrated analysis of the aircraft and the engine were conducted for flight Mach numbers of M=5.0, 4.0, 3.5, 3.0, and 2.0. Increasing the angle of attack $\alpha$ of the aircraft in M=5.0 flight increased the mass flow rate of the air captured at the intake due to pre-compression effect of the nose shockwave, also increasing the thrust obtained at the engine plug nozzle. Sufficient thrust for acceleration were obtained at $\alpha=3$ and 5 degrees. Increase of flight Mach number at $\alpha=0$ degrees resulted in decrease of mass flow rate captured at the engine intake, and thus decrease in thrust at the nozzle. The thrust was sufficient for acceleration at M=3.5 and lower cases. Lift force on the aircraft was increased by the integration of engine on the aircraft for all varying angles of attack or flight Mach numbers. However, the slope of lift increase when increasing flight Mach number showed decrease as flight Mach number reach to M=5.0, due to the separation shockwave at the upper surface of the aircraft. Pitch moment of the aircraft was not affected by the installation of the engines for all angles of attack at M=5.0 condition. In low Mach number cases at $\alpha=0$ degrees, installation of the engines increased the pitch moment compared to no engine configuration. Installation of the engines increased the frictional drag on the aircraft, and its percentage to the total drag ranged between 30-50% for varying angle of attack in M=5.0 flight.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 추계학술대회논문집B
• /
• pp.438-443
• /
• 2001

Plate impingement of the impulse wave discharged from the open end of a duct is numerically investigated using a CFD method. Harten-Yee Total Variation Diminishing method is used to solve the unsteady, compressible flow governing equations. The Mach number, the flat plate inclination and the distance between the duct exit and inclined flat plate are changed to investigate their effects on the impinging flow field. The impulse wave impingement on the inclined flat plate depends on Mach number $M_s$ and the plate inclination $\psi$. The pressure distributions on the inclined flat plate show that for a small r/D, the peak pressure at the center of an inclined flat plate decreases with an increase in the plate inclination $\psi$ in the range of $\psi$ from $45^{\circ}$ to $60^{\circ}$ but for a large r/D, the peak pressure decreases with an increase in $\psi$ in the range of $\psi$ from $75^{\circ}$ to $90^{\circ}$. It is also found that for all of r/D, the peak pressure at the center of an inclined flat plate has a maximum value in $\psi=90^{\circ}$.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 1995년도 제5회 학술강연회논문집
• /
• pp.69-83
• /
• 1995

A great number of experimental data indicating shock-induced separation(SIS) in internal or external supersonic flows were reviewed to make clear the mechanism of SIS and to present the criterion of turbulent boundary layer separation. The interesting conclusions were obtained for the considerably wide range of flow geometries that the incipient separation is almost independent of the flow geometries, and that it is relatively unaffected by changes in gas specific heat, and boundary layer Reynolds number, Furthermore, the pressure rise necessary to separate boundary layer in external flows was found to be applicable to SIS in overexpanded propulsion nozzles. This is due to the fact that the SIS phenomenon caused by the interaction between shock waves and turbulent boundary layers is processed through a supersonic deceleration. This is, the SIS in almost all of interacting flow fields is governed by the concept of free interaction, and criterion of SIS is only a Function of upstream Mach number.

• International Journal of Aeronautical and Space Sciences
• /
• 제18권1호
• /
• pp.82-90
• /
• 2017

The aim of this work is to develop a new computational program to determine the effect of using the gas of propulsion of combustion chamber at high temperature on the shape of the two-dimensional Minimum Length Nozzle giving a uniform and parallel flow at the exit section using the method of characteristics. The selected gases are $H_2$, $O_2$, $N_2$, CO, $CO_2$, $H_2O$, $NH_3$, $CH_4$ and air. All design parameters depend on the stagnation temperature, the exit Mach number and the used gas. The specific heat at constant pressure varies with the temperature and the selected gas. The gas is still considered as perfect. It is calorically imperfect and thermally perfect below the threshold of dissociation of molecules. A error calculation between the parameters of different gases with air is done in this case for purposes of comparison. Endless forms of nozzles may be found based on the choise of $T_0$, $M_E$ and the selected gas. For nozzles delivering same exit Mach number with the same stagnation temperature, we can choose the right gas for aerospace manufacturing rockets, missiles and supersonic aircraft and for supersonic blowers as needed in settings conception.