• 대한기계학회논문집B
• /
• 제24권9호
• /
• pp.1210-1218
• /
• 2000

Thermal convection in a horizontal annulus is considered, and the bifurcation phenomena of flows from time-periodic to chaotic convection are numerically investigated. The unsteady two-dimensional streamfunction-vorticity equation is solved with finite difference method. As Rayleigh number is increased, the steady flow bifurcates to a time-periodic flow with a fundamental frequency, and afterwards a period-tripling bifurcation occurs with further increase of the Rayleigh number. Chaotic convection is established after a period-doubling bifurcation. A periodic convection with period 4 appears after the first chaotic convection. At still higher Rayleigh numbers, chaotic flows reappear.

• 설비공학논문집
• /
• 제13권2호
• /
• pp.88-95
• /
• 2001

Transition to chaotic convection is investigated for natural convection of a fluid with Pr=0.1 in a wide-gap horizontal annuls. The unsteady two-dimensional stream-function-vorticity equation is solved with finite difference method. As the Rayleigh number is increased, the steady 'downward flow' bifurcates to a time-periodic flow with a fundamental frequency, and afterwards a period-doubling bifurcation occurs. As the Rayleigh number is increased further, the chaotic flow regime is reached after a sequence of successive Hopf bifurcation to quasi-periodic and chaotic flow regimes. The route to chaos shows the Ruelle-Takens-Newhouse scenario. The flow of chaotic regime displays complex coalescence and separation of eddies in the side and lower region of the annulus.

• 대한기계학회논문집B
• /
• 제22권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.

• 대한기계학회논문집B
• /
• 제25권3호
• /
• pp.433-441
• /
• 2001

Natural convection of a fluid with intermediate Prand시 number of Pr=0.2 in a horizontal annulus is considered, and the bifurcation phenomena and chaotic flows are numerically investigated. The unsteady two-dimensional streamfunction-vorticity equation is solved with finite difference method. The steady downward flow with two counter-rotating eddies bifurcates to a simple periodic flow with a fundamental frequency. And afterwards, second Hopf bifurcation occurs, and a quasi-periodic flow with two incommensurable frequencies appears. However, a new time-periodic flow is established after experiencing quasi-periodic states. As Rayleigh number is increased further, the chaotic flow regime is reached after a sequence of successive Hopf bifurcation to quasi-periodic and chaotic flow regimes. A scenario similar to the Ruelle-Takens-Newhouse scenario of the onset of chaos is observed.

• 한국전산유체공학회지
• /
• 제16권3호
• /
• pp.52-58
• /
• 2011

This study investigates the bifurcation sequence to chaos in a horizontal annulus with a constant heat flux wall. After the first Hopf bifurcation from a steady to a simple time-periodic flow with a fundamental frequency, quasi-periodic flows with two or three incommensurable frequencies appear. A reverse transition from a quasi-periodic flow to a simple periodic flow is observed with increase of Rayleigh number. And finally, chaotic convection is established after appearance of three incommensurable frequencies at a high Rayleigh number. Simple periodic flows exist between quasi periodic flows. The transition route to chaos of the present simulations follows the Ruelle-Takens route.

• 대한기계학회논문집B
• /
• 제22권12호
• /
• pp.1784-1795
• /
• 1998

Natural convection of low Prandtl number fluids with $Pr{\leq}0.2$ in a narrow horizontal annulus is numerically investigated. For $Pr{\leq}0.2$, hydrodynamic instability induces oscillatory multicellular flows consisting of multiple like-rotating cells. For a fluid with $Pr{\approx}0$, the region in which instability of conduction regime first forms is near the vertical section of annulus, and the multiple cells are distributed uniformly in the lower and upper regions of annulus. As Pr increases, however, the cells are shifted upwards. The like-rotating cells drift downward, as time goes on, and the speed of travel increases with increase of Pr. For a fluid with Pr=0.1, a flow with period-4 solution is observed between chaotic states.

• 대한기계학회논문집B
• /
• 제22권3호
• /
• pp.360-372
• /
• 1998

This study deals with the interaction between buoyancy-induced convection and externally imposed excitation in the form of harmonic rocking and the effect of the interaction upon heat transfer in low-Pr liquids. A wide array of system responses are discussed using the spectral collocation numerical technique. The superposition of buoyancy and Coriolis forces leads to complex fluid flow and heat transfer. The transition to chaotic convection is accelerated, and heat transfer rates are reduced as the enclosure is excited at the fundamental frequency of oscillation associated with the pure buoyancy-driven case. Average heat transfer rates are correlated for Pr=0.02 and 0.03. The heat transfer is affected more in the Pr=0.03 liquid than the case of Pr=0.02.

• 한국전산유체공학회지
• /
• 제11권2호
• /
• pp.49-55
• /
• 2006

This study investigates the oscillatory thermal convection of a fluid with Pr=0.02 in a wide-gap horizontal annulus with constant heat flux inner wall. When Pr=0.02, dual steady-state flows are not found. After the first Hopf bifurcation from a steady to a time-periodic flow, five successive period-doubling bifurcations are recorded before chaos. The power spectrum shows the $period-2^4\;and\;2^5$ flows clearly, and a window of period $3{\times}2^3$ flow is found in the chaotic regime. The approximate value of the Feigenbaum number for the last three period-doubling bifurcations is 4.76. The transition route to chaos of the present simulations is consistent with the period-doubling route of Feigenbaum.

• 설비공학논문집
• /
• 제6권1호
• /
• pp.20-28
• /
• 1994

A numerical study has been performed on enhancement of heat transfer in a forced convection of the modified driven-cavity flow which was previously found by the author to give a regular or chaotic stirring depending on the parameter value. It is found that for the present case wherein heat is transmitted between fluid and the surrounding walls, the chaotic stirring enhances the heat transfer at high Peclet numbers. The optimal condition of the flow modulation for the best heat transfer can be predicted by purely investigating the hydrodynamic facet, i.e. the stirring effect.

• 대한기계학회논문집B
• /
• 제21권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.