• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.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.

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

Two-phase RoeM and AUSMPW+ schemes are preconditioned for the simulation of all Mach number flows, which are generally of interest for many gas-liquid two-phase application problems, because of large speed of sound in liquid region and low speed of sound in mixture or gas region. Conventional characteristic based schemes lose their accuracy or robustness in low Mach number flows, because their numerical dissipation terms are scaled by speed of sound, which is too large compared with local velocity magnitude in a low Mach region. All speed versions of RoeM and AUSMPW+ reflect the eigenvalues of the preconditioned governing system, which have the same order of magnitude even in low Mach number region. From the asymptotic analysis, it is observed that the discretized system by the developed schemes is consistent with the continuum system in the incompressible limit. The numerical results show the accurate and robust behavior of the proposed shcemes for all speed two-phase flows.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• 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.

• Journal of computational fluids engineering
• /
• v.17 no.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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.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.

• 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 Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• 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.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• 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}$.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 1995.11a
• /
• 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
• /
• v.18 no.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.