• Journal of the Korean Society of Safety
• /
• v.35 no.3
• /
• pp.96-104
• /
• 2020

The numerical simulations were conducted to investigate the thermal-fluid phenomena occurred inside the experimental apparatus during a PCCS, used to remove heat released in accidents from a containment of light water nuclear power plant, operation. Numerical simulations of the flow and heat transfer caused by wall condensation inside the containment simulation vessel (CSV), which equipped with 18 vertical heat exchanger tubes, were conducted using the commercial computational fluid dynamics (CFD) software ANSYS-CFX. Shear stress transport (SST) and the wall condensation model were used for turbulence closure and wall condensation, respectively. The simulation using the actual size of the apparatus. However, rather than simulating the whole experimental apparatus in consideration of the experimental cases, calculation resources, and calculation time, the simulation model was prepared only in CSV. Selective simulation was conducted to verify the effects of non-condensable gas(NC gas) concentration, CSV internal pressure, and wall sub-cooling conditions. First, as a result of the internal flow of CSV, it was observed that downward flow due to condensation occurred surface of the vertical tube and upward flow occurred in the distant place. Natural convection occurred actively around the heat exchanger tube. Due to this rising and falling internal flow, natural circulation occurred actively around the heat exchanger tubes. Next, in order to check the performance of built-in condensation model using according to the non-condensable gas concentration, CSV internal flow and wall sub-cooling, the heat flux values were compared with the experimental results. On average, the results were underestimated with and error of about 25%. In addition, the influence of CSV internal pressure and wall sub-cooling was small, but when the condensate was highly generated due to the low non-condensable gas concentration, the error was large compared to the experimental values. This is considered to be due to the nature of the condensation model of the CFX code. However, in spite of the limitations of CFD, it is valid to use the built-in condensation model of CFD for PCCS performance prediction from a conservative perspective.

• 한국신재생에너지학회:학술대회논문집
• /
• 2008.05a
• /
• pp.33-36
• /
• 2008

The transient behavior of a passive air breathing direct methanol fuel cell (DMFC) operated on vapor-feeding mode is studied in this paper. It generally takes 30 minutes after starting for the cell response to come to its steady-state and the response is sometimes unstable. A mathematical dynamic one-dimensional model for simulating transient response of the DMFC is presented. In this model a DMFC is decomposed into its subsystems using lumped model and divided into five layers, namely the anodic diffusion layer, the anodic catalyst layer, the proton exchange membrane (PEM), the cathodic catalyst layer and the cathodic diffusion layer. All layers are considered to have finite thickness, and within every one of them a set of differential-algebraic governing equations are given to represent multi-components mass balance, such as methanol, water, oxygen and carbon dioxide, charge balance, the electrochemical reaction and mass transport phenomena. A one-dimensional, isothermal and mass transport model is developed that captures the coupling between water generation and transport, oxygen consumption and natural convection. The single cell is supplied by pure methanol vapor from a methanol reservoir at the anode, and the oxygen is supplied via natural air-breathing at the cathode. The water is not supplied from external source because the cell uses the water created at the cathode using water back diffusion through nafion membrane. As a result of simulation strong effects of water transport were found out. The model analysis provides several conclusions. The performance drop after peak point is caused by insufficiency of water at the anode. The excess water at the cathode makes performance recovery impossible. The undesired crossover of the reactant methanol through the PEM causes overpotential at the cathode and limits the feeding methanol concentration.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.4
• /
• pp.1019-1030
• /
• 1994

Natural convection heat transfer from a horizontal ice cylinder immersed in quiescent cold pure water was studied experimentally. The experiment was conducted for the ambient water temperatures ranging from $2.0^{\cric}C$ to $10.0^{\circ}C$. The flow fields around an ice cylinder and its melting shapes were visualized and local Nusselt numbers obtained. Especially, its attention was focused on the density maximum effects and stagnation point Nusselt number. From the visualized photographs of flow fields, three distinct flow patterns were observed with the ambient water temperature variation. The melting shapes of ice cylinder are various in shape with flow patterns. Steady state upflow was occured at the range of $2.0^{\circ}C \leq T_{\infty} \leq 4.6^{\circ}C$ and steady state downflow was occured at $T_{\infty} \geq 6.0^{\circ}C$. In the range of $4.7^{\circ}C < T_{\infty} < 6.0^{\circ}C$, three-dimensional unsteady state flow was observed. Especially, the melting shapes of ice cylinder have formed the several spiral flutes for the temperatures ranging from $5.5^{\circ}C$ to $5.8^{\circ}C$. For upflow regime, the maximum stagnation point Nusselt number exists at $T_{\infty} = 2.5^{\circ}C$ and as the ambient water temperature increases the Nusselt number decreases. At ambient water temperature of about $5.7^{\circ}C$, Nusselt number shows its minimum value.

• Journal of Energy Engineering
• /
• v.7 no.1
• /
• pp.113-121
• /
• 1998

In natural convection flows, the fluid velocities are highly dependent on the thermal field and property variations can have a strong effect on both the flow and thermal fields. To examine the effect of property variations, at first, numerical analyses covering wide range of the Prandtl number under the same Rayleigh numbers have been carried out. Next, we have modeled the viscosity and thermal conductivity as parabolic functions of temperature and a comprehensive set of numerical solutions have been obtained to understand the effect. The Prandtl number dependence of Nusselt number is fairly strong even though the effect is still weak compared to the Rayleigh number dependence. When thermophysical properties are dependent on temperature, the flow field showed a fairly weak variation except near boundaries, whereas the temperature field is strongly affected, especially by the temperature dependent thermal conductivity.

• Solar Energy
• /
• v.19 no.3
• /
• pp.29-41
• /
• 1999

Three-dimensional natural convective heat transfer in a vertical channel with a protruding single module was investigated experimentally. The particular interest was in the removal of thermal energy from the module by convective heat transfer. Hence radiative and conductive heat losses were estimated by using thermocouples and heat flux sensor respectively. The flow fields in the channel were visualized by means of a smoke-method. Also, local temperatures were measured by thermocouples inside the channel, along the vertical wall and module surface. It is found that convective heat transfer was promoted at the lower comer of the module and was decreased at the upper comer due to a recirculation zone. A general correlation of the critical channel ratios was found as a function of Rayleigh number. For the range of $8.28{\times}10^3<Ra^*_c<3.48{\times}10^6$, a useful correlation for the mean Nusselt number was proposed as a function of modified channel Rayleigh number.

• Nuclear Engineering and Technology
• /
• v.25 no.1
• /
• pp.8-19
• /
• 1993

The natural convection model of the consolidated system has been developed to make sure the removal of decay heat generated in the spent fuel for the loss of forced cooling accident. The numerical technique employed was based on the ADI scheme. The calculation of heat generation rate in the spent fuel was peformed by the ANS-79 decay heat model, and the nonuniform surface heat flux is assumed with a chopped sine curve for the conservative decay heat generation input. The sensitivity study was performed to examine the possibility of the pool bulk boiling by varying the various parameters, i.e. inter-fuel spacing ratio, heat generation power, and radius of the fuel rod. The application results of this model show that the natural circulation flow through compacted spent fuel bundles enables the pool temperature to control in a safe and effective manner, after the required cooling time. The corresponding acceptance criteria of the cooling time for rearranging the spent fuel rods were also found.

• Journal of the Korean Magnetics Society
• /
• v.23 no.2
• /
• pp.77-81
• /
• 2013

The objective of this study is experimentally to investigate natural convective heat transfer characteristics of the ferrofluid for a concentric annuli under rotating magnetic field with variations of the revolution and the magnetic field strength. The rotating magnetic field was provided by induction motor with 6 poles and 3 phases and the revolution and the magnetic field strength were controlled by an inverter driver and a voltage meter, respectively. Temperatures of the inner pipe and the outer pipe in the tested concentric annuli were maintained at $30^{\circ}C$ and $25^{\circ}C$, respectively, during the test and the direction of the rotating magnetic field was a counterclockwise. As a result, the natural convective heat transfer characteristics of the ferrofluid for a concentric annuli were increased with the rise of the revolution and magnetic field strength due to the increased heat dissipation between hot side and cold side of the concentric annuli.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.37 no.11
• /
• pp.1423-1428
• /
• 2013

In this study, a heat dissipation measurement system is developed to analyze a 1.2-kW BLDC motor. It is important to check the temperature of the motor because an increase in temperature causes problems in the motor insulations, which in turn influences the motor life. A generator for a vehicle is installed to set up a load. We changed the load from 165 to 495 W. While the rpm varies from 2000 to 4000 under various load conditions, the changes in temperature were measured for the operating period by using a thermocouple. The results of experiments conducted under natural convection conditions suggest that the temperature was not stationary with the rpm, load, and coil of the motor and it kept increasing over $120^{\circ}C$. However, under forced convection conditions, the temperature stationarily reached $84^{\circ}C$ after 4000 s. The difference between the maximum and the minimum temperatures was $10-26^{\circ}C$ with an increase in the rpm and load. The orders of high temperature were as follows: motor coil (Ch#1), side of motor surface (Ch#5), inside of motor cap (Ch#2), upper side of motor surface (Ch#4), and inner wall of the motor (Ch#3).

• Journal of the Korean Wood Science and Technology
• /
• v.44 no.3
• /
• pp.370-380
• /
• 2016

The black pine logs damaged by pine wilt disease in Jeju-do were heat-treated to extend the utilization of domestic trees damaged by pine wilt disease. The heat-treatment of wood requires wood to be heated to $56^{\circ}C$ for 30 min at the core. The average moisture content and top-diameter of the black pine logs were ranged from 46% to 141% and from 180 mm to 500 mm, respectively. And the basic specific gravity and oven-dry specific gravity of the black pine logs were 0.47 and 0.52, respectively. The time required for heat-treatment at $105^{\circ}C$ temperature was ranged from 7.7 h to 44.2 h, depending on moisture content and top-diameter. The temperature distribution was used to predict the time required for heat-treatment of black pine log with various moisture contents and top-diameters using finite difference method. The thermal properties of wood including the thermal conductivity and specific heat in accordance with moisture content were calculated. Heat transfer coefficient for mixed convection in form of adding natural convection and forced convection was used for heat transfer analysis. The error between the measured and predicted values ranged from 3% to 45%. The predicted times required for heat-treatment of black pine log with 50% moisture content and 200 mm, 300 mm, and 400 mm top-diameter were 10.9 h, 18.3 h, and 27.0 h, respectively. If the initial moisture content of black pine log is 75%, heat treatment times of 13.6 h, 22.5 h, and 32.8 h were predicted in accordance with top-diameter. And if the initial moisture content of black pine log is 100%, heat treatment times of 16.2 h, 26.5 h, and 38.2 h were predicted in accordance with top-diameter. When the physical properties of logs damaged by pine wilt disease are presented, these results can be applicable to the heat-treatment of red pine and Korean pine logs as well.

• Korean Journal of Ecology and Environment
• /
• v.42 no.1
• /
• pp.85-94
• /
• 2009

This study examines the operational effectiveness of a Convection Current Aeration System (CCAS) in reservoir. CCAS was run from June, 2008 when the thermocline begun forming in the reservoir. This paper reviews the influence of stratification, dissolved oxygen dynamics and temperature in the lake's natural state from June to October 2008. The survey was done on a week basis. Upwelling flow effects a radius of $7{\sim}10m$ at a surface directly and was irrelevant to the strength of thermocline. On the other hand, it was affected the number of working days, and strength of thermocline at vertical profiles of the reservoir. Longer CCAS run, the deeper was the vertical direct flow area. However it didn't break the thermocline during summer season of 2008. The operating efficiency of the CCAS in the reservoir depends on hydraulics and meteological conditions. Computational Fluid Dynamics (CFD) is a very useful tool for evaluating the operating efficiency of fluid dynamics. The geometry for CFD simulation consists of a cylindrical vessel 25 m radius and 40 m height. The CCAS is located in center of domain. The non-uniform tetrahedral meshes had a bulk of the geometry. The meshes ranged from the coarse to the very fine. This is attributed to the cold water flowing into the downcomer and rising, creating a horizontal flow to the top of the CCAS. The result of CFD demonstrate a closer agreement with surveyed data for temperature and flow velocity. Theoretical dispersion volume were calculated at 8m depth, 120 m diameter working for 30 days and 10 m depth, 130 m diameter working for 50 days.