• Journal of computational fluids engineering
• /
• v.13 no.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.03a
• /
• 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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.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
• /
• v.54 no.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
• /
• v.45 no.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
• /
• v.34 no.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.

• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• 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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.38 no.11
• /
• pp.935-942
• /
• 2014

A numerical analysis of the effect of the Prandtl number on the natural convection in a cold outer tilted square enclosure with an inner hot circular cylinder is presented. Several Prandtl numbers (Pr=0.1, 0.7, 7) are considered, with different angles($0^{\circ}$, $15^{\circ}$, $30^{\circ}$, $45^{\circ}$) for the enclosure and Rayleigh numbers ($Ra=10^3$, $10^4$, $10^5$). The effect of the Prandtl number on the natural convection is analyzed using isotherms and streamline and surface-averaged Nusselt numbers. The flow and heat transfer characteristics are found to be dependent on the time for $Ra=10^5$ and Pr=0.1 at angles of $0^{\circ}$ and $45^{\circ}$. However, in the other cases, the flow and heat transfer characteristics are independent of the time.The surfaceaveraged Nusselt number increases with an increase in the Prandtl number. As the Prandtl number increases, the Nusselt number becomes larger regardless of the angle for $Ra=10^5$. In particular, the Nusselt number steeply increases when the angle is $45^{\circ}$ for $Ra=10^5$ and Pr=0.1.

• Korean Chemical Engineering Research
• /
• v.57 no.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
• /
• v.5 no.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.