• The KSFM Journal of Fluid Machinery
• /
• v.19 no.4
• /
• pp.19-28
• /
• 2016

Objective of this study is to compare the inherent characteristics of natural convection flow inside the canister of spent fuel dry storage system with different backfill gases by utilizing computational fluid dynamics (CFD) code. Four working fluids were selected for comparison study. Helium currently used backfill gas for canister, air, nitrogen, and argon are frequently used as coolant in many heat transfer applications. The results indicate that helium has very distinct conductive behavior and show very weak natural convective flow compared to the others. Argon showed the strongest natural convective flow but also the worst coolability. Air and nitrogen showed similar characteristics to each other. However, due to difference in Prandtl number, nitrogen showed more effective natural convective flow. These results suggest that experimental validation for the nitrogen is needed to investigate the potential coolability other than currently commercially used helium.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.20 no.4
• /
• pp.309-332
• /
• 2016

Finite element method has been applied to solve the fundamental governing equations of natural convective, electrically conducting, incompressible fluid flow past an infinite vertical plate surrounded by porous medium in presence of thermal radiation, viscous dissipation, Soret and Dufour effects. In this research work, the results of coupled partial differential equations are found numerically by applying finite element technique. The sway of significant parameters such as Soret number, Dufour number, Grashof number for heat and mass transfer, Magnetic field parameter, Thermal radiation parameter, Permeability parameter on velocity, temperature and concentration evaluations in the boundary layer region are examined in detail and the results are shown in graphically. Furthermore, the effect of these parameters on local skin friction coefficient, local Nusselt number and Sherwood numbers is also investigated. A very good agreement is noticed between the present results and previous published works in some limiting cases.

• Journal of the Korean Institute of Gas
• /
• v.16 no.2
• /
• pp.9-13
• /
• 2012

The present study has been carried out to diagnose the temperature safety of a boiler desuperheater which decreases abnormally higher temperature of superheated steam in a reheater of a power plant. The liquid water in the desuperheater stays in a closed space. It becomes heated by the high temperature superheated steam and boiling could occur. Boiling might increase internal pressure and it could destroy the desuperheater if the internal pressure exceeds the allowable pressure of the desuperheater. The present study modeled reasonably the desuperheater and four cases of heat transfer analysis are executed with the consideration of insulator and natural convective fluid flow of the inside cooling water. For the case excluded the natural convective fluid flow, the temperature exceeds the allowable temperature and pressure. On the other hand, for the real case included natural convective fluid flow and insulator, the active heat transfer from higher temperature region to lower temperature region occurs and it makes the temperature in the cooling water below the allowable temperature and pressure. From this fact, it could be thought that the desuperheater in the reheater is safe from destroy or back flow.

• Journal of Power System Engineering
• /
• v.19 no.4
• /
• pp.36-42
• /
• 2015

TMixed convective flow in a bottom heated and top cooled rectangular channel can be significantly affected by the channel aspect ratio, Prandtl number, Reynolds number, Rayleigh number and angle of inclination. In such a mixed convection, the flow pattern plays an important role in various technological processes. In this study, a numerical investigation is carried out to explore mixed convection in a three-dimensional rectangular channel with bottom heated and top cooled uniformly. The three-dimensional governing equations are discretized using the finite volume method. In the range of low Reynolds number($0{\leq}Re{\leq}9.6{\times}10^{-2}$), the effects of the aspect ratio($2{\leq}AR{\leq}12$) and Gr/Re are presented and discussed. The longitudinal roll number in the channel is increased with increasing aspect ratio, and the roll number induced, regardless of the aspect ratio number, is even in the range of aspect ratios between 2 and 12, New vortex flow structure containing inclined longitudinal rolls is found, which is affected by aspect ratio and Reynolds number. The ratio Gr/Re is used to check the relative magnitudes of forced and natural convection in the mixed convective flow of high viscous fluid.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2004.10a
• /
• pp.107-110
• /
• 2004

Application of aluminum alloy has been increasing for most of industry area because aluminum has a good mechanical properties and castability, especially automotive field for weight reduction. But, Furnace industry is sluggish. The purpose of this study is numerical analysis of aluminum holding furnace for reasonableness estimation when we design for new model of furnace. The numerical simulation involving fluid flow of inside air and heat transfer to fireproof material is presented in order to improve the understanding of aluminum furnace. First of all, we are carried out numerically for the two dimensional inside convection and surface radiation heat transfer in a square enclosure. Subsequently, we are established the analysis method of aluminum furnace considering natural convective heat transfer

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.3
• /
• pp.555-568
• /
• 1992

A asic study is performed on two-dimensional natural convective flow and heat transfer in a fluid heated internally within an inclined rectangular enclosures. For Rayleign numbers from 1.0*10$^{4}$to 1.5*10$^{5}$ , aspect ratio of 1/4, 1/3 and 1/2, and inclined angle from 0deg to 90deg, the governing equations were solved numerically and the experiments were performed by MachZehnder interferometer using low salinity water as a test fluid. For aspect ratios adapted in present study, the natural convection occures the most intensive at inclined angle of 0deg. This became weak at inclined angles of 60deg and 30deg in case of aspect ratios of 1/3 and 1/2 respectively. The intensity of flow was roughly in proportion to Rayliegh numbers and in proportion to the forth power of aspect ratios.

• Journal of computational fluids engineering
• /
• v.6 no.3
• /
• pp.1-9
• /
• 2001

Mixed convection in a horizontal annulus is considered, and the effect of a forced flow on the natural convective flow is investigated. The inner cylinder is hotter than the outer cylinder, and the outer cylinder is rotating with constant angular velocity with its axis at the center of the annulus. The unsteady streamfunction-vorticity equation is solved with a finite difference method. For the fluid with Pr=0.7, there appear flows with two eddies, one eddy, or no eddy according the Rayleigh and Reynolds numbers. The rotation of the outer cylinder reduces the heat transfer rate at the wall of the annulus. The oscillatory multicellular flow of a low Prandtl number fluid with Pr=0.01 can be effectively suppressed by the forced flow.

• Journal of computational fluids engineering
• /
• v.4 no.1
• /
• pp.8-18
• /
• 1999

In devising a numerical approximation for the convective spatial transport of a fluid mechanical quantity, it is noted that the convective motion of a scalar quantity occurs in one-way, or from upstream to downstream. This consideration leads to a new scheme termed a pure upwind difference scheme (PUDS) in which an estimated value for a fluid mechanical quantity at a control surface is not influenced from downstream values. The formal accuracy of the proposed scheme is third order accurate. Two typical benchmark problems of a wall-driven fluid flow in a square cavity and a buoyancy-driven natural convection in a tall cavity are computed to evaluate performance of the proposed method. for comparison, the widely used simple upwind scheme, power-law scheme, and QUICK methods are also considered. Computation results are encouraging: the proposed PUDS sensitized to the convection direction produces the least numerical diffusion among tested convection schemes, and, notable improvements in representing recirculation of fluid stream and spatial change of a scalar. Although the formal accuracy of PUDS and QUICK are the same, the accuracy difference of approximately a single order is observed from the revealed results.

• Journal of Mechanical Science and Technology
• /
• v.16 no.5
• /
• pp.683-695
• /
• 2002

The objective of the present study is to analyze the fluid flow with moving boundary using a finite element method. The algorithm uses a fractional step approach that can be used to solve low-speed flow with large density changes due to intense temperature gradients. The explicit Lax-Wendroff scheme is applied to nonlinear convective terms in the momentum equations to prevent checkerboard pressure oscillations. The ALE (Arbitrary Lagrangian Eulerian) method is adopted for moving grids. The numerical algorithm in the present study is validated for two-dimensional unsteady flow in a driven cavity and a natural convection problem. To extend the present numerical method to engine simulations, a piston-driven intake flow with moving boundary is also simulated. The density, temperature and axial velocity profiles are calculated for the three-dimensional unsteady piston-driven intake flow with density changes due to high inlet fluid temperatures using the present algorithm. The calculated results are in good agreement with other numerical and experimental ones.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.6
• /
• pp.804-810
• /
• 2001

Transition of flows in natural convection in a horizontal cylindrical annulus is investigated for the fluid with Pr=0.2. The unsteady streamfunction-vorticity equation is solved with finite difference method. As Rayleigh number is increased, the steady crescent-shaped eddy flow bifurcates to a time-periodic flow with like-rotating eddies. After the first Hopf bifurcation, however, a reverse transition from oscillatory to a steady flow occurs by the flow pattern variation. Hysteresis phenomenon occurs between the solution branches of up-scan and down-scan stages, and dual solutions with one steady and one oscillatory flow are found. Overall Nusselt of the flows at the flows at the down-scan stage is greater than that at the up-scan stage.