• Water for future
• /
• v.31 no.4
• /
• pp.291-297
• /
• 1998

The two dimensional time dependent flow past a circular cylinder is analyzed numerically. In the analysis, equations of conservation of mass and momentum are transformed to equations of stream function-vorticity and vorticity transport, and nondimensionalized by nondimensional parameters representing flow characteristics, The resulting stream function-vorticity equation and vorticity transport equation are solved by successive over relaxation scheme and alternating direct implicit scheme. Numerical experiments are performed for the flow in the range of Reynolds number 125 to 275. The time dependent streamlines, vorticities, pressure on cylinder surface, separation angle, and drag and lift coefficients are calculated, and the method for estimation of pressure on cylinder surface and the outer boundary limit are developed.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.7 no.1
• /
• pp.40-48
• /
• 2003

Large eddy simulation is used to investigate the compressible flow over a open cavity, The sub-grid scale stresses are modeled using the dynamic model. The compressible Navier-Stokes equations are solved with the sixth order accurate compact finite difference scheme in the space and the 4th order Runge-Kutta scheme in the time. The results show a typical flow pattern of the shear layer mode of oscillation over the cavity. The votical disturbances, the roll-up of vorticity, and impingement and scattering of vorticity at the downstream cavity edge can be seen in the shear layer. Predicted acoustic resonant frequency is in good agreement with that of the empirical formula. The mean flow streamlines are nearly horizontal along the mouth of the cavity. The pressure has its minimum value in the vortex core inside the cavity.

• Journal of Mechanical Science and Technology
• /
• v.17 no.1
• /
• pp.129-138
• /
• 2003

A 1 : 20 laboratory scale test rig of a 200 MW tangentially fired boiler is built up with completely simulated structures such as platen heaters and burners. Iso-thermal turbulent flow in the boiler is mapped by 3-D PDA (Particle Dynamic Analyzer). The 3-D numerical models for the same case are proposed based on the solution of к-$\varepsilon$ model closed RANS (Reynolds time-Averaged Navier-Stokes) equations, which are written in the framework of general coordinates and discretized in the corresponding body-fitted meshes. Not only are the grid lines arranged to fit the inner/outer boundaries. but also to align with the streamlines to the best possibility in order to reduce the NDE (numerical diffusion errors). Extensive comparisons of profiles of mean velocities are carried out between experiment and calculation. Predicted velocities in burner region were quantitatively similar with measured ones, while those in other area have same tendency with experimental counterpart.

• Journal of computational fluids engineering
• /
• v.20 no.1
• /
• pp.16-25
• /
• 2015

The present study investigates the sloshing phenomena in a two-dimensional rectangular tank at a low filling level by using a level set method based on finite volume method. The code validations are performed by comparing between the present results and previous numerical and experimental results, which gives a good agreement. Various excitation frequencies and excitation amplitude of the 30% filling height tank have been considered in order to observe the dependence of the sloshing behavior on the excitation frequency and amplitude. Regardless of excitation amplitude, the maximum value of wall pressure occurs when the excitation frequency reaches the natural frequency. The time sequence of free surface and corresponding streamlines for excitation frequencies have been presented to analysis the variation of wall pressure according to time, which contributes to explain the double peaks in the time variation of wall pressure.

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

This paper describes the shape optimization of a stator blade in a single-stage transonic axial compressor. The blade optimization has been performed using response surface method and three-dimensional Navier-Stokes analysis. Two shape variables of the stator blade, which are used to define a stacking line, are introduced to increase an adiabatic efficiency. Data points for response evaluations have been selected by D-optimal design, and linear programming method has been used for an optimization on a response surface. Throughout the shape optimization of a stator blade, the adiabatic efficiency is increased to 5.8 percent compared to that of the reference shape of the stator. The increase of the efficiency is mainly caused by the pressure enhancement in the stator blade. Flow separation on the blade suction surface of the stator is also improved by optimizing the stator blade. It is noted that the optimization of the stator blade is also useful method to increase the adiabatic efficiency in the axial compressor as well as the optimization of a rotor blade, which is widely used now.

• Journal of the Society of Naval Architects of Korea
• /
• v.50 no.3
• /
• pp.167-174
• /
• 2013

Model ice forming process in ice tank needs several steps of seeding, freezing, tempering. In those process, one of the most important factors to affect the accuracy of experiment is the homogeneity of the ice thickness and the temperature. This paper investigated a computational and statistical method to assess the uniformity of the model ice. In addition, the different configurations of freezing systems were considered to improve the uniformity. Qualitative assessment using streamlines from the cooling units was carried out by computational fluid dynamics (CFD) and the quantitative evaluations of the homogeneity were compared using the temperature distribution on the ice surface. In addition, multi species transport analysis is introduced to understand the circulation efficiency of cold air from the cooling units. As the results, optimized configurations were determined by adjusting the angles of vane in the cooling units.

• Journal of the Korean Society of Visualization
• /
• v.12 no.2
• /
• pp.18-22
• /
• 2014

Caenorhabditis elegans (C. elegans) is an undulatory nematode which exhibits two distinct locomotion types of swimming and crawling. Although in its natural habitat C. elegans lives in a non-Newtonian fluidic environment, our current understanding has been limited to the behavior of C. elegans in a simple Newtonian fluid. Here, we present some experimental results on the penetrating behavior of C. elegans at the interface from liquid to solid environment. Once C. elegans, which otherwise swims freely in a liquid, makes a contact to the solid gel boundary, it begins to penetrate vertically to the surface by changing its stroke motion characterized by a stiffer body shape and a slow stroke frequency. The particle image velocimetry (PIV) analysis reveals the flow streamlines produced by the stroke of worm. For the worm that crawls on a solid surface, we utilize a technique of traction force microscopy (TFM) to find that the crawling nematode forms localized force islands along the body where makes direct contacts to the gel surface.

• Journal of the Korean Society of Safety
• /
• v.12 no.1
• /
• pp.77-93
• /
• 1997

The natural convection and combined heat transfer induced by fire in a rectangular enclosure is numerically studied. The model for this numerical analysis is partially opened, it is divided by a vertical baffle projecting from ceiling. The solution procedure Includes the standard k- $\varepsilon$ model for turbulent flow and the discrete ordinates method (DOM ) is used for the calculation of radiative heat transfer equation. In this study, numerical simulation on the combined naturnal convection and radiation is carried out in a partial enclosure filled with absorbed-emitted gray media, but is not considered scattering problem. The velocity vectors, streamlines, and isothermal lines are compared the results of pure convection with those of the combined convection-radiation, the combined heat transfer. Comparing the results of pure convection with those of the combined convection-radiation, the combined heat transfer analysis shows the stronger circulation than those of the pure convection. Three different locations of heat source are considered to observe the effect of heat source location on the heat transfer phenomena. As the results, the circulation and the heat transfer In the left region from heating block are much more influenced than those in the right region. It is also founded that the radiation effect cannot be neglected in analyzing the building in fire.

• Journal of computational fluids engineering
• /
• v.21 no.3
• /
• pp.1-7
• /
• 2016

In the present study, the homogeneous model is used to simulate the natural convection heat transfer of the CuO-water nanofluid in a concentric annular enclosure. Simulations have been carried while the Rayleigh number ranges from $10^3$ to $10^6$, solid volume fraction ranges from 0.01 to 0.04 and the radius ratio varies between 0.1 and 0.7. Results are presented in the form of streamlines, isotherm patterns and averaged Nusselt numbers for different values of solid volume fraction, radius ratio of the annulus and Rayleigh numbers. The results show that by decreasing the radius ratio and/or increasing the Rayleigh number, the averaged Nusselt number increases. Also the heat transfer rate increases as increased solid volume fractions.

• Journal of computational fluids engineering
• /
• v.21 no.3
• /
• pp.48-54
• /
• 2016

In the present work, numerical simulations were conducted on the scaled model of the BWB type UCAV in the subsonic region using ANSYS FLUENT V15. The prediction method was validated through comparison with experimental results and the effect of the twisted wing was investigated. To consider the transitional flow phenomenon, ${\gamma}$ transition model based on SST model was adopted. The coefficients of lift, drag and pitching moment were compared with experimental results and the pressure distribution and streamlines were investigated. The twisted wing decreases the lift force but increases lift-to-drag ratio through delay of stall and leading edge vortex's movement to the front, also the non-linearity of the pitching moment is decreased.