• Journal of Power System Engineering
• /
• v.13 no.5
• /
• pp.11-17
• /
• 2009

A numerical study of unsteady mixed convection in a cavity with high viscous fluid is presented. Finite volume method was employed for the discretization and PISO algorithm was used for calculating pressure term. The parameters governing the problem are the Rayleigh number ($10^3\;{\leq}\;Ra\;{\leq}\;10^5$), the Reynolds number (0 < Re $\leq$ 1), and the aspect ratio (0.5 $\leq$ AR $\leq$ 2). The fluid used is silicon oil, a high prandtl number fluid, Pr = 909.1. The results show velocity vectors and temperature distributions. It is found that the periodic flows in a cavity are observed at very low Reynolds numbers, and the period of periodic flow decreases with increasing Reynolds and Rayleigh numbers, and increases with increasing aspect ratio. Also, the Reynolds number range of periodic flow increases with increasing Rayleigh numbers and aspect ratio.

• Journal of the Korean Society of Manufacturing Technology Engineers
• /
• v.21 no.2
• /
• pp.196-201
• /
• 2012

This study represents numerical study on the thermal and fluid flow characteristics of exhaust gas in a motorcycle muffler. The engine generates 125cc of displacement. Numerical analysis with computational fluid dynamics(CFD) was carried out to investigate the exhaust gas that flow into a motorcycle muffler. The STAR-CD S/W used to analyze three dimensional steady state and transient in a muffler. The Navier-Stokes Equation is solved with SAMPLE and PISO method in cartesian coordinates system. As the numerical result, it could be confirmed pulsating pressure generated from inlet of muffler become closer to the atmospheric pressure through baffle. Also the numerical result was almost identical to experimental result and, through the comparison, the reliability of numerical result was confirmed.

• Proceedings of the Korean Society of Marine Engineers Conference
• /
• 2005.06a
• /
• pp.210-215
• /
• 2005

A three-dimensional numerical simulation is performed to investigate on a low Reynolds number mixed convection in a horizontal rectangular channel with the upper part cooled and the lower part heated uniformly. The three-dimensional governing equations are solved using a finite volume method. For convective term, the central differencing scheme is used and for the pressure correction, the PISO algorithm is used. Solutions are obtained for A=4, Pr=0.72, 10, 909, the Reynolds number ranging from $2.1{\times}10^{-2}$ to $1.2{\times}10^{-1}$, the Rayleigh number is $3.5{\times}10^4$. It is found that vortex roll structures of mixed convection in horizontal rectangular channel can be classified into three roll structures which affected by Prandtl number and Reynolds number.

• Journal of Power System Engineering
• /
• v.11 no.4
• /
• pp.56-64
• /
• 2007

A numerical simulation is carried out mixed convection in horizontal channel with a heat source from below of rectangular cavity. Finite volume method was employed for the discretization and PISO algorithm was used for calculating pressure term. The parameters governing the problem are the Reynolds number ($10^{-2}{\leq}Re{\leq}50$), the Rayleigh number ($10^3{\leq}Ra{\leq}2.06{\times}10^5$), the Prandtl number ($0.72{\leq}Pr{\leq}909$), the aspect ratio ($0.5{\leq}AR=W/H{\leq}2$) and the angle of inclination ($0^{\circ}{\theta}60^{\circ}$). Mean Nusselt number distributions were obtained and effect of Reynolds number, Rayleigh number and Prandtl number on mixed convection in the horizontal channel with rectangular cavity were investigated.

• Journal of computational fluids engineering
• /
• v.16 no.3
• /
• pp.22-28
• /
• 2011

Laminar Flow over an equilateral triangular cylinder is studied for several inflow angles. Under an uniform flow of $Re_d$=50,75,100,125,150, the triangular cylinder is rotated by ${\theta}$=$0^{\circ}$,$15^{\circ}$,$30^{\circ}$,$45^{\circ}$,$60^{\circ}$,$75^{\circ}$,$90^{\circ}$,$105^{\circ}$. The governing equations are solved by the PISO algorithm based on the finite volume method of the unstructured grid system. The effects of the inflow angle on the vortex-shedding flows are investigated. The Strouhal number shows a minimum at ${\theta}$=$60^{\circ}$. It is closely related to the variation of pressure and flow structure induced by the movement of separation points.

• Korean Journal of Computational Design and Engineering
• /
• v.1 no.1
• /
• pp.56-64
• /
• 1996

A quadrilateral-triangular mixed grid method for the solution of incompressible viscous flow is presented. The solution domain near the body surface is meshed using elliptic grid geneator to acculately simulate the viscous flow. On the other hand, we used unstructured triangular grid system generated by advancing front technique of a simple automatic grid generation algorithm in the rest of the computational domain. The present method thus is capable of not only handling complex geometries but providing accurate solutions near body surface. The numerical technique adopted here is PISO type finite element method which was developed by the present author. Investigations have been made of two-dimensional unsteady flow of Re=550 past a circular cylinder. In the case of use of the unstructured grid only, there exists a considerable amount of difference with the existing results in drag coefficient and vorticity at the cylinder surface; this may be because of the lack of the grid clustering to the surface that is a inevitable requirement to resolve the viscous flow. However, numerical results on the mixed grid show good agreements with the earlier computations and experimental data.

• International Journal of Air-Conditioning and Refrigeration
• /
• v.10 no.3
• /
• pp.117-128
• /
• 2002

This paper describes the smoke filling process of a fire field model based on a self-deve-loped SMEP (Smoke Movement Estimating Program) code to the simulation of fire induced flows in the two types of atrium space containing a ceiling heat flux. The SMEP using PISO algorithm solves conservation equations for mass, momentum, energy and species, together with those for the modified k-$\varepsilon$ turbulence model with buoyancy production term. Also it solves the radiation equation using the discrete ordinates method. Compressibility is assumed and the perfect gas law is used. Comparison of the calculated upper-layer average tempera-ture and smoke layer clear height with the zone models has shown reasonable agreement. The zone models used are the CFAST and the NBTC one-room. For atrium fires with ceiling glass the ceiling heat flux by solar heat causes a high smoke temperature near the ceiling. However, it has no effect on the smoke movement such as the smoke layer clear heights that are important in fire safety. In conclusion, the smoke layer clear heights that are important in evacuation activity except the early of a fire were not as sensitive as the smoke layer tem-perature to the nature of ceiling heat flux condition. Thus, a fire sensor in atrium with ceiling glass has to consider these phenomena.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.14 no.11
• /
• pp.914-922
• /
• 2002

Thermo-acoustic waves can be thermally generated in a compressible flow field by rapid heating and cooling, and chemical reaction near the boundary walls. This mechanism is very important in the space environment in which natural convection does not exist. Also this may be a significant factor for heat transfer when the fluids are close to the thermodynamic critical point. In this study, the generation and transmission characteristics of thermo-acoustic waves in an air-filled confined domain with two-step pulsed heating are studied numerically. The governing equations are discretized using control volume method, and are solved using PISO algorithm and second-order upwind scheme. For the purpose of stable solution, time step was set to the order of $1\times10_-9s,\;and\;grids\;are\;50\times2000$. Results show that temperature and pressure distributions of fluid near the boundary wall subjected to a rapid heating are increased abruptly, and the induced thermo-acoustic wave propagates through the fluid until it decays due to viscous and heat dissipation. Pressure waves have sharp front shape and decay with a long tail in the case of step heating, but these waves have sharp pin shape in the case of pulsed heating.