• 한국전산유체공학회:학술대회논문집
• /
• 1997.10a
• /
• pp.10-15
• /
• 1997

Unsteady natural convection of an enclosed fluid has been one of the fundamental thermo-fluid problems, of which dynamic relevance is found in many engineering applications. Together with the inherent coupling between the boundary layers and the interior core, and strong interaction between flow and temperature fields, the unsteadiness poses serious hurdles for analytical and experimental approaches. With the recent development of computers and solution algorithms, computational fluid dynamics has become the prevailing tool to tackle the underlying problems. In this presentation, a few examples of numerical studies are introduced. The usefulness and potential of numerical simulations in investigating unsteady natural convection are elaborated.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.2 no.1
• /
• pp.63-73
• /
• 1990

A numerical analysis of initial unsteady state flow and heat transfer in an enclosed cavity has been performed by the Modified QUICK Scheme. The stable QUICK Scheme which modified the coefficient always to be positive is included in this numerical analysis. The implicit method is applied to solve the unsteady state flow; between iterations the PISO (Pressure - Implicit with Splitting of Operators) algorithm is employed to correct and update the velocity and pressure fields on a staggered grid. The accuracy of the Modified QUICK Scheme is proved by applying fewer grid systems than those which Ghia et al. and Davis applied. The initial unsteady mixed convection in an enclosed cavity is analyzed using the above numerical procedure. This study focuses on the development of the large main vortex and secondary vortex in forced convection, the effects of the Rayleigh Number in natural convection and the relative direction of the forced and natural convection.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.7
• /
• pp.919-927
• /
• 1997

The main purpose of this study is to determine bifurcation as the primary instability of a glass melting furnace. Steady-state and unsteady characteristics of natural convection in the partially open cavity as appeared in a glass melting furnace is investigated by using numerical analysis. Three types of convection, such as steady laminar, unsteady periodic or unsteady quasi-periodic convection may occur according to the temperature difference between upper two isothermal surfaces along the depth of cavity in a glass melting furnace. In the temperature difference of 150-900 K between batch and free surface, the larger the temperature difference, the weaker the convection strength and unsteadiness. Since the glass viscosity is increasing exponentially in the lower temperature, the batch freezes the thermofluidic field especially below the surface of it. If the depth of cavity is 0.5 m, the bifurcation to time-dependent natural convection may occur in the range of 60-650 K. If that is 1.0 m, it may occur in the whole range of temperature difference.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.286-288
• /
• 2003

The physical model considered here is a horizontal layer of fluid heated below and cooled above with a periodic array of evenly spaced square cylinders placed at the center of the layer, whose aspect ratio here varies from unity to six. Periodic boundary condition is employed along the horizontal direction to allow for lateral freedom for the convection cells. Two-dimensional solution for unsteady natural convection is obtained using an accurate and efficient Chebyshev spectral multi-domain methodology for a given Rayleigh numbers of $10^6$

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.12
• /
• pp.1774-1783
• /
• 1998

The transition from steady laminar to chaotic convection in a glass melting furnace specified by upper surface temperature distribution has been studied by the direct numerical analysis of the two and three-dimensional time dependent Navier-Stokes equations. The thermal instability of convection roll may take place when modified Rayleigh number($Ra_m$) is larger than $9.71{\times}10^4$. It is shown that the basic flows in a glass melting furnace are steady laminar, unsteady periodic, quasi-periodic or chaotic flow. The dimensionless time scale of unsteady period is about the viscous diffusion time, ${\tau}_d=H^2/{\nu}_0$. Through primary and secondary instability analyses the fundamental unsteady feature in a glass melting furnace is well defined as the unsteady periodic or weak chaotic flow.

• International Journal of Aeronautical and Space Sciences
• /
• v.3 no.1
• /
• pp.95-104
• /
• 2002

Numerical solution of the full Navier-Stokes equation as well as the energy equation has been obtained for the unsteady natural convection in a rectangular enclosure. One side wall was maintained at very high temperature simulating fires. Especially the effect of surface radiation was taken into account. While the enclosed air was assumed to be transparent, the internal walls directly interacted one another through the surface radiation. Due to a significant temperature difference in the flow field, the equation of state was used instead of the Boussinesq approximation. It was found that the rapid heating of the adiabatic ceiling and floor by the incoming radiation from the hot wall made the evolution at thermo-fluid field highly unstable in the initial period. Therefore, the secondary cells brought about at the floor region greatly affected the heat transfer mechanism inside the enclosure. The heat transfer rate was augmented by the radiation, resulting in requiring less time for the flow to reach the steady state. At the steady state neglecting radiation two internal hydraulic jumps were clearly observed in upper/left as well as in lower/right comer. However, the hydraulic jump in the lower/right comer could not be observed for the case including radiation due to its high momentum flow over the bottom wall. Radiation resulted in a faster establishment of the steady state phenomena.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.302-310
• /
• 2011

A fluid in an enclosure can be heated by electric heating, chemical reaction, or fission heat. In order to remove the volumetric heat of the fluid, the walls surrounding the enclosure must be cooled. In this case, a natural convection occurs in the pool of the fluid, and it has a dominant role in heat transfer to the surrounding walls. It can augment the heat transfer rates tens to hundreds times larger than conductive heat transfer. The heat transfer by a natural convection in a regular shape such as a square cavity or semi-circular pool has been studied experimentally and numerically for many years. A pool of an inverted triangular shape with 10 degree inclined bottom walls has a good cooling performance because of enhanced boiling critical heat flux (CHF) compared to horizontal downward surface. The coolability of the pool is determined by comparing the thermal load from the pool and the maximum heat flux removable by cooling mechanism such as radiative or boiling heat transfer on the pool boundaries. In order to evaluate the pool coolability, it is important to correctly expect the thermal load by a natural convection heat transfer of the pool. In this study, turbulence models with modifications for buoyancy effect were validated for unsteady natural convections by volumetric heating. And natural convection in the triangular pool was evaluated by using the models.

• Journal of computational fluids engineering
• /
• v.16 no.3
• /
• pp.66-74
• /
• 2011

A fluid in an enclosure can be heated by electric heating, chemical reaction, or fission heat. In order to remove the volumetric heat of the fluid, the walls surrounding the enclosure must be cooled. In this case, a natural convection occurs in the pool of the fluid, and it has a dominant role in heat transfer to the surrounding walls. It can augment the heat transfer rates tens to hundreds times larger than conductive heat transfer. The heat transfer by a natural convection in a regular shape such as a square cavity or semi-circular pool has been studied experimentally and numerically for many years. A pool of an inverted triangular shape with 10 degree inclined bottom walls has a good cooling performance because of enhanced boiling critical heat flux (CHF) compared to horizontal downward surface. The coolability of the pool is determined by comparing the thermal load from the pool and the maximum heat flux removable by cooling mechanism such as radiative or boiling heat transfer on the pool boundaries. In order to evaluate the pool coolability, it is important to correctly expect the thermal load by a natural convection heat transfer of the pool. In this study, turbulence models with modifications for buoyancy effect were validated for unsteady natural convections by volumetric heating. And natural convection in the triangular pool was evaluated by using the models.

• Journal of computational fluids engineering
• /
• v.5 no.3
• /
• pp.40-46
• /
• 2000

Natural convection in a wide-gap horizontal annulus is considered, and the transition of flows from steady to oscillatory convection is investigated for the fluid with Pr=0.1. The unsteady streamfunction-vorticity equation is solved with finite difference method. As Rayleigh number is increased, the steady crescent-shaped flow bifurcates to a time-periodic flow with like-rotating eddies. And afterwards, a transition to an oscillatory multicellular flow with a counter-rotating eddy on the top of the annulus occurs. A transition from steady to an oscillatory flow occurs, but dual solutions and hysteresis phenomena are not observed.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.3 s.258
• /
• pp.273-282
• /
• 2007

Numerical calculations are carried out for the natural convection induced by temperature difference between a cold outer square cylinder and a hot inner circular cylinder. A two-dimensional solution for unsteady natural convection is obtained, using the immersed boundary method (IBM) to model an inner circular cylinder based on finite volume method, for different Rayleigh numbers varying over the range of $10^4\;to\;10^6$. The study goes further to investigate the effect of an inner cylinder location on the heat transfer and fluid flow. The location of inner circular cylinder is changed vertically along the center-line of square enclosure. The number, size and formation of cell strongly depend on Rayleigh number and the position of inner circular cylinder. The changes in heat transfer quantities have been presented.