• 한국전산유체공학회지
• /
• 제13권3호
• /
• pp.62-68
• /
• 2008

This paper reports briefly on the computational results of a turbulent Rayleigh-Benard convection with the elliptic-blending second-moment closure (EBM). The primary emphasis of the study is placed on an investigation of accuracy and numerical stability of the elliptic-blending second-moment closure for the turbulent Rayleigh-Benard convection. The turbulent heat fluxes in this study are treated by the algebraic flux model with the temperature variance and molecular dissipation rate of turbulent heat flux. The model is applied to the prediction of the turbulent Rayleigh-Benard convection for Rayleigh numbers ranging from Ra=$2{\times}10^6$ to Ra=$10^9$ and the computed results are compared with the previous experimental correlations, T-RANS and LES results. The predicted cell-averaged Nusselt number follows the correlation by Peng et al.(2006) (Nu=$0.162Ra^{0.286}$) in the 'soft' convective turbulence region ($2{\times}10^6{\leq}Ra{\leq}4{\times}10^7$) and it follows the experimental correlation by Niemela et al. (2000) (N=$0.124Ra^{0.309}$) in the 'hard' convective turbulence region ($10^8{\leq}Ra{\leq}10^9$) within 5% accuracy. This results show that the elliptic-blending second-moment closure with an algebraic flux model predicts very accurately the Rayleigh-Benard convection.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 학술대회
• /
• pp.111-117
• /
• 2008

This paper reports briefly on the computational results of a turbulent Rayleigh-Benard convection with the elliptic-blending second-moment closure (EBM). The primary emphasis of the study is placed on an investigation of accuracy and numerical stability of the elliptic-blending second-moment closure for the turbulent Rayleigh-Benard convection. The turbulent heat fluxes in this study are treated by the algebraic flux model with the temperature variance and molecular dissipation rate of turbulent heat flux. The model is applied to the prediction of the turbulent Rayleigh-Benard convection for Rayleigh numbers ranging from $Ra=2{\times}10^6$ to $Ra=10^9$, and the computed results are compared with the previous experimental correlations, T-RANS and LES results. The predicted cell-averaged Nusselt number follows the correlation by Peng et al.(2006) ($Nu=0.162Ra^{0.286}$) in the 'soft' convective turbulence region ($2{\times}10^6{\leq}Ra{\leq}4{\times}10^7$) and it follows the experimental correlation by Niemela et al. (2000) ($Nu=0.124Ra^{0.309}$) in the 'hard' convective tubulence region ($10^8{\leq}Ra{\leq}10^9$) within 5% accuracy. This results show that the elliptic-blending second-moment closure with an algebraic flux model predicts very accurately the Rayleigh Benard convection.

• 대한기계학회논문집B
• /
• 제33권6호
• /
• pp.435-442
• /
• 2009

Three-dimensional heat transfer characteristics for natural convection flows are numerically investigated in the doubly-inclined cubical-cavity according to the variation of a newly defined orientation angle �� of the hot wall surface from horizontal plane at moderate Rayleigh numbers. Numerical simulations of laminar flows are conducted in the range of Rayleigh numbers($10^4{\leq}Ra{\leq}10^5$) and $0^{\circ}{\leq}{\alpha}90^{circ}$ with a solution code(PowerCFD) employing unstructured cell-centered method. Comparisons of the average Nusselt number at the cold face are made with benchmark solutions and experimental results found in the literature. It is found that the average Nusselt number at the cold wall has a maximum value around the specified orientation ${\alpha}$ at each Rayleigh number. Special attention is also paid to three-dimensional thermal characteristics in natural convection according to new orientation angles at Ra��= $1{\times}10^5$, in order to investigate a new additional heat transfer characteristic found in the range of above Ra = $6{\times}10^4$.

• Nuclear Engineering and Technology
• /
• 제54권9호
• /
• pp.3540-3550
• /
• 2022

The supercritical carbon dioxide (S-CO₂) Brayton cycle is an important energy conversion technology for the fourth generation of nuclear energy. Since the printed circuit heat exchanger (PCHE) used in the S-CO₂ Brayton cycle has narrow channels, Rayleigh-Bénard (RB) convection is likely to exist in the tiny channels. However, there are very few studies on RB convection in supercritical fluids. Current research on RB convection mainly focuses on conventional fluids such as water and air that meet the Boussinesq assumption. It is necessary to study non-Boussinesq fluids. PRB convection refers to RB convection that is affected by horizontal incoming flow. In this paper, the computational fluid dynamics simulation method is used to study the PRB convection phenomenon of non-Boussinesq fluid-supercritical carbon dioxide. The result shows that the inlet Reynolds number (Re) of the horizontal incoming flow significantly affects the PRB convection. When the inlet Re remains unchanged, with the increase of Rayleigh number (Ra), the steady-state convective pattern of the fluid layer is shown in order: horizontal flow, local traveling wave, traveling wave convection. If Ra remains unchanged, as the inlet Re increases, three convection patterns of traveling wave convection, local traveling wave, and horizontal flow will appear in sequence. To characterize the relationship between traveling wave convection and horizontal incoming flow, this paper proposes the relationship between critical Reynolds number and relative Rayleigh number (r).

• Nuclear Engineering and Technology
• /
• 제45권7호
• /
• pp.969-978
• /
• 2013

A full-sized model for the horizontally oriented metal cask containing 21 spent fuel assemblies has been considered to evaluate the internal natural convection behavior within a dry shield canister (DSC) filled with helium as a working fluid. A variety of two-dimensional CFD numerical investigations using a turbulent model have been performed to evaluate the heat transfer characteristics and the velocity distribution of natural convection inside the canister. The present numerical solutions for a range of Rayleigh number values ($3{\times}10^6{\sim}3{\times}10^7$) and a working fluid of air are further validated by comparing with the experimental data from previous work, and they agreed well with the experimental results. The predicted temperature field has indicated that the peak temperature is located in the second basket from the top along the vertical center line by effects of the natural convection. As the Rayleigh number increases, the convective heat transfer is dominant and the heat transfer due to the local circulation becomes stronger. The heat transfer characteristics show that the Nusselt numbers corresponding to $1.5{\times}10^6$ < Ra < $1.0{\times}10^7$ are proportional to 0.5 power of the Rayleigh number, while the Nusselt numbers for $1.0{\times}10^7$ < Ra < $8.0{\times}10^7$ are proportional to 0.27 power of the Rayleigh number. These results agreed well with the trends of the experimental data for Ra > $1.0{\times}10^7$.

• Journal of Advanced Marine Engineering and Technology
• /
• 제34권1호
• /
• pp.53-61
• /
• 2010

In the present study, a laminar natural convection flow of $H_2O$-Cu nanofluid in a two dimensional enclosure has been investigated using a thermal lattice Boltzmann approach with the Bhatnagar-Gross-Krook (BGK) model. The effect of suspended nanoparticles on the fluid flow and heat transfer process have been studied for different controlling parameters such as particle volume fraction ($\Phi$), Rayleigh number (Ra). For this investigation the Rayleigh number changes from 104 to 106 and volume fraction varied from 0 to 10% with three different particle diameters (dp), say 10 nm, 20 nm and 40 nm. It is shown that increasing the Rayleigh number (Ra) and the volume fraction of nanofluid causes an increase of the effective heat transfer rate in terms of average Nusselt number (Nu) as well as the thermal conductivity of nanofluid. On the other hand, increasing the particle diameter causes the decrease of the heat transfer rate and thermal conductivity. The result of the analysis are compared with experimental and numerical data both for pure and nanofluids and it is seen a relatively good agreement.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2008년도 추계학술대회B
• /
• pp.2744-2749
• /
• 2008

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 for Rayleigh number of $Ra=10^7$. This study investigates the effect of the inner cylinder location on the heat transfer and fluid flow. The location of inner circular cylinder ($\delta$) is changed vertically along the center-line of square enclosure. The natural convection bifurcates from unsteady to steady state according to $\delta$. Two critical positions of ${\delta}_{C,L}$ and ${\delta}_{C,U}$ as a lower bound and an upper bound are ${\delta}_{C,L}=0.05$ and ${\delta}_{C,U}=0.18$, respectively. Within the defined bounds, the thermal and flow fields are steady state. When the inner cylinder locates at ${\delta}{\geq}{\delta}_{C,U}$, the space between the upper surface of inner cylinder and the top surface of the enclosure forms a relatively shallow layer where the natural convection characterized as the pure Rayleigh-Benard convection forms alternately the upwelling and downwelling plums, as a result that a series of cells known as Benard cells is derived.

• 대한기계학회논문집B
• /
• 제38권11호
• /
• pp.935-942
• /
• 2014

본 연구에서는 내부에 고온의 원형 실린더가 위치한 저온의 기울어진 사각 밀폐계에서 Prandtl 수 변화에 따른 밀폐계 내부 자연대류 현상에 대한 수치해석을 수행하였다. Rayleigh 수가 $10^3$, $10^4$, $10^5$ 그리고 밀폐계의 각도가 $0^{\circ}$, $15^{\circ}$, $30^{\circ}$, $45^{\circ}$ 인 경우에 대해서 Prandtl 수를 0.1, 0.7, 7 로 변화시키며 Prandtl 수의 영향을 분석하였다. Prandtl 수의 변화에 따른 밀폐계 내부의 자연대류 현상은 실린더 표면과 밀폐계 벽면의 평균 Nusselt 수, 밀폐계 내부의 등온선 및 유선을 바탕으로 분석되었다. Prandtl 수가 0.7, 7 인 경우에는 각도 및 Rayleigh 수에 관계없이 열유동은 정상상태의 특성을 보였지만 Prandtl 수가 0.1 인 경우에는 Rayleigh 수가 $10^5$ 이고 밀폐계의 각도가 $0^{\circ}$ 와 $45^{\circ}$ 일 때 시간에 따라 열유동장이 변하였다. 실린더 표면과 밀폐계 벽면의 평균 Nusselt 수는 Rayleigh 수가 $10^5$ 인 경우, 각도에 관계없이 Prandtl 수가 클수록 증가하였다. 특히 Rayleigh 수가 $10^5$ 이고 Prandtl 수가 0.1 일 때 밀폐계의 각도가 $45^{\circ}$ 인 경우 실린더 표면과 밀폐계 벽면의 평균 Nusselt 수가 크게 증가하였다.

• Korean Chemical Engineering Research
• /
• 제57권1호
• /
• pp.142-147
• /
• 2019

The onset of buoyancy-driven convection in an initially isothermal and quiescent horizontal fluid layer was analyzed theoretically. It is well-known that at the critical Rayleigh number $Ra_c=669$ convective motion sets in with a constant-heat-flux cooling through the upper boundary. Here, based on the momentary instability concept, the dimensionless critical time ${\tau}_m$ to mark the onset of convective motion for Ra > 669 was analyzed theoretically. The energy method under the momentary stability concept was used to find the critical conditions as a function of the Rayleigh number Ra and the Prandtl number Pr. The predicted critical conditions were compared with the previous theoretical and experimental results. The momentary stability criterion gives more reasonable wavenumber than the conventional energy method.

• Advances in aircraft and spacecraft science
• /
• 제5권4호
• /
• pp.411-430
• /
• 2018

This computational study examines the augmentation of classic 2-D Rayleigh-$B{\acute{e}}nard$ convection by the addition of periodically-spaced transverse fins. The fins are attached to the heated base of the cavity and serve to partition the cavity into 'units' with different aspect ratios. The respective impacts upon heat transfer of the fin configuration parameters - including spacing, height, thickness and thermal conductivity - are systematically examined through numerical simulations for a range of laminar Rayleigh numbers (0 < Ra < $2{\times}10^5$) and reported in terms of an average Nusselt number. The selection of the low Rayleigh number regime is linked to likely scenarios within aerospace applications (e.g. avionics cooling) where the cavity length scale and/or gravitational acceleration is small. The net heat transfer augmentation is found to result from a combination of competing fin effects, most of which are hydrodynamic in nature. Heat transfer enhancement of up to $1.2{\times}$ that for a Rayleigh-$B{\acute{e}}nard$ cavity without fins was found to occur under favorable fin configurations. Such configurations are generally characterized by short, thin fins with half-spacings somewhat less than the convection cell diameter from classic Rayleigh-$B{\acute{e}}nard$ theory. In contrast, for unfavorable configurations, it is found that the introduction of fins can result in a significant reduction in the heat transfer performance.