• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.7
• /
• pp.2315-2325
• /
• 1996

A compression wave is attenuated or distorted as it propagates in a tube. The present study investigated the propagation characteristics of the compression waves which are generated by a train in a high-speed railway tunnel. A Total Variation Diminishing (TVD) difference scheme was applied to one-dimensional, unsteady viscous compressible flow. The numerical calculation involved the effects of wall friction, heat transfer and energy loss due to the friction heat in the boundary layer behind the propagating compression wave, and compared with the measurement results of a shock tube and a real tunnel. The present results show that attenuation of the compression wave in turbulent boundary layer is stronger than in laminar boundary layer, but nonlinear effect of the compression wave is greater in the laminar boundary layer. The energy loss due to the frictional heat had not influence on attenuation and distortion of the propagating compression waves.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.5
• /
• pp.28-33
• /
• 2005

Nonlinear aeroelastic characteristics of a two dimensional typical section model with bilinear plunge spring are investigated. Doublet-point method(DPM) is used for the calculation of supersonic unsteady aerodynamic forces which are approximated by using the minimum-state approximation. For nonlinear flutter analysis structural nonlinearity is represented by an asymmetric bilinear spring and is linearized by using the describing function method. The linear and nonlinear flutter analyses indicate that the flutter characteristics are significantly dependent on the frequency ratio. From the nonlinear flutter analysis, various types of limit cycle oscillations are observed in a wide range of air speeds below or above the linear flutter boundary. The nonlinear flutter characteristics and the nonlinear aeroelastic responses are investigated.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2005.05a
• /
• pp.821-825
• /
• 2005

Tilt rotor aircraft was developed for satisfying VSTOL (vertical short take off and landing) capability and cruise performance. However the noise generated by tilt rotor system causes one of the most serious problems. In this paper, noise characteristics of tilt rotor system in hovering flight are predicted by using free wake method and Lowson's formula. The flow field of the tilt rotor is simulated by using time marching free wake method, and the free field acoustic pressure is calculated through Lowson's formula. The predicted results are compared with experimental data at various observing positions. In the near field, they show good agreement with experimental data regardless of rotating speed and collective pitch angles of 6, 8 and 10 degree, although there are some discrepancies between prediction and experiment in the far field and at the rotating axis in the near field. It seems that the reason of these discrepancies is difference of unsteady force fluctuation between experiment and calculation.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.263-270
• /
• 2010

In this study, an air cleaner is considered to improve comfort, safety, and health of automobile passengers. The performance and installation conditions of the air cleaner have been studied to investigate their effects on the air quality in the cabin room using numerical analysis. A five-passenger sedan and a seven-passenger minivan that have comparatively large indoor volume have been considered. The distributions of the local mean age and the volume averaged age of indoor air are calculated according to the variation of the placement and the air flow of the air cleaner. In addition, a decrease of contamination concentration, especially VOCs(volatile organic compounds), by the air cleaner is numerically analyzed with time-accurate unsteady calculation to quantify the effect of the air cleaner on the indoor air quality. As a result, the effective installation and operation conditions of the air cleaner for the automobile cabin room could be suggested.

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

High-speed flight vehicle have various cavities. The supersonic cavity flow is complicated due to vortices, flow separation and reattachment, shock and expansion waves. The general cavity flow phenomena include the formation and dissipation of vortices, which induce oscillation and noise. The oscillation and noise greatly affect flow control, chemical reaction, and heat transfer processes. The supersonic cavity' flow with high Reynolds number is characterized by the pressure oscillation due to turbulent shear layer, cavity geometry, and resonance phenomenon based on external flow conditions, The resonance phenomena can damage the structures around the cavity and negatively affect aerodynamic performance and stability. In the present study, we performed numerical analysis of cavities by applying the unsteady, compressible three dimensional Reynolds-Averaged Navier-Stokes(RANS) equations with the ${\kappa}-{\omega}$ turbulence model. The cavity model used for numerical calculation had a depth(D) of 15mm cavity aspect ratio(L/D) of 3, width to spanwise ratio(W/D) of 1.0 to 5.0. Based on the PSD(Power Spectral Density) and CSD(Cross Spectral Density) analysis of the pressure variation, the dominant frequency was analyized and compared with the results of Rossiter's Eq.

• Journal of computational fluids engineering
• /
• v.17 no.1
• /
• pp.94-100
• /
• 2012

In this paper, the cavitating flows around a hydrofoil have been numerically investigated by using a 2-d multi-phase RANS flow solver based on pseudo-compressibility and a homogeneous mixture model on unstructured meshes. For this purpose, a vertex-centered finite-volume method was utilized in conjunction with 2nd-order Roe's FDS to discretize the inviscid fluxes. The viscous fluxes were computed based on central differencing. The Spalart-Allmaras one equation model was employed for the closure of turbulence. A dual-time stepping method and the Gauss-Seidel iteration were used for unsteady time integration. The phase change rate between the liquid and vapor phases was determined by Merkle's cavitation model based on the difference between local and vapor pressure. Steady state calculations were made for the modified NACA66 hydrofoil at several flow conditions. Good agreements were obtained between the present results and the experiment for the pressure coefficient on a hydrofoil surface. Additional calculation was made for cloud cavitation around the hydrofoil. The observation of the vapor structure, such as cavity size and shape, was made, and the flow characteristics around the cavity were analyzed. Good agreements were obtained between the present results and the experiment for the frequency and the Strouhal number of cavity oscillation.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.365-368
• /
• 2002

The design of industrial arc plasma systems is still largely based on trial and error although the situation is rapidly improving because of the available computational power at a cost which is still fast coming down. The desire to predict the behavior of arc plasma system, thus reducing the development cost, has been the motivation of arc research. To interrupt fault current, the most enormous duty of a circuit breaker, is achieved by separating two contacts in a interruption medium, $SF_{6}$ gas or air etc., and arc plasma is inevitably established between the contacts. The arc must be controlled and interrupted at an appropriate current zero. In order to analyze arc behavior in $SF_{6}$ gas circuit breakers, a numerical calculation method combined with flow field and electromagnetic field has been developed. The method has been applied to model arc generated in the Aachen nozzle and compared the results with the experimental results. Next, we have simulated the unsteady flow characteristics to be induced by arcing of AC cycle, and conformed that the method can predict arc behavior in account of thermal transport to $SF_{6}$ gas around the arc, such as increase of arc voltage near current zero and dependency of arc radius on arc current to maintain constant arc current density.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.5 s.236
• /
• pp.609-616
• /
• 2005

While the SIMPLE algorithm is most widely used for the simulations of flow phenomena that take place in the industrial equipments or the manufacturing processes, it is less adopted for the simulations of the free surface flow. Though the SIMPLE algorithm is free from the limitation of time step, the free surface behavior imposes the restriction on the time step. As a result, the explicit schemes are faster than the implicit scheme in terms of computation time when the same time step is applied to, since the implicit scheme includes the numerical method to solve the simultaneous equations in its procedure. If the computation time of SIMPLE algorithm can be reduced when it is applied to the unsteady free surface flow problems, the calculation can be carried out in the more stable way and, in the design process, the process variables can be controlled based on the more accurate data base. In this study, a modified SIMPLE algorithm is presented fur the free surface flow. The broken water column problem is adopted for the validation of the modified algorithm (MoSIMPLE) and for comparison to the conventional SIMPLE algorithm.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.11
• /
• pp.1422-1436
• /
• 1997

Numerical investigations were conducted to study the convective phenomena of an initially stably stratified salt water solution with lateral heating in a uniformly rotating annulus. The method of investigation is the finite difference analysis of the basic conservation equation for an axisymmetric, unsteady, double-diffusive convection and calculation is made for R $a_{\eta}$=2*10$^{5}$ and Ta=10$^{7}$ ~ 2.5*10$^{8}$ . Formation of layered flow structure, merging process of layers, the corresponding temperature and concentration distributions, Nu variation with time are examined. Numerical results show that in each layer, the temperature profile looks 'S'-shaped and the concentration profile is uniform due to the convective mixing. At the interface between adjacent layers, the temperature changes smoothly but the concentration changes rapidly. As the effect of the rotation increases, the generation of rolls at hot wall, the formation and merging of layers are delayed. The average Nu shows the trend of conduction heat transferees the effect of the rotation increases.n increases.

• Solar Energy
• /
• v.17 no.3
• /
• pp.3-11
• /
• 1997

Recently the construction of atrium buildings has increased but along with it many problems in thermal environment have arised. since the exterior wall of glass, indoor temperature is greatly influenced by weather conditions and since the space volume is very large, the vertical air temperature is not uniform. So, in this study, a Vertical Temperature Distribution Model was developed to predict the vertical air temperature of an atrium and evaluate the effects of the design parameters on the air temperature distribution of an atrium. To consider the characteristics of the vertical air temperature distribution in an atrium, the Satosh Togari's Macroscopic Model was used basically for the calculation of the vertical air temperature distribution in large space and the solar radiation analysis model and natural ventilation analysis model in atrium. And to calculate the unsteady-state inside wall surface temperature(boundary condition), the finite difference method was used. For the verification of the developed temperature distribution program, numerical evaluation of air flow by the ${\kappa}-{\varepsilon}$ turbulence model and in-situ test was conducted in parallel. The results of this study, the developed temperature distribution program was seen to predict the thermal condition of the atrium very accurately.