• 한국안전학회지
• /
• 제35권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.

• 대한기계학회논문집B
• /
• 제41권1호
• /
• pp.37-46
• /
• 2017

가열된 표면에서의 액적은 일반적으로 내부에 마랑고니 유동이 발생하고, 이는 불균일한 침전물 패턴 형상을 구성하게 된다. 본 연구는 마랑고니 유동을 가시화하고, 수직진동을 사용해서 이를 제어하는 것을 주 목적으로 한다. 액적이 증발하는 동안, 액적의 접촉각 변화와 부피변화를 실험적으로 알 수 있었고, PIV(Particle Image Velocimetry) 실험 장치를 이용하여, 평판 온도별 마랑고니 유동의 내부 유동의 흐름을 가시화하였다. 그리고 평판에 각 주파수별 수직진동을 가해주는 실험을 진행하여, 그 결과 마랑고니 유동의 유동 방향과 수직진동의 유동 방향이 서로 반대인 것을 확인할 수 있었다. 마지막으로 증발하는 액적에 수직진동을 가해줌으로써, 액적의 하단부분에서 내부유동의 흐름변화를 관찰하였다. 마랑고니유동에 의해 발생하는 내부유동 방향과 수직진동으로 발생하는 내부유동의 방향이 서로 반대 방향이므로 가열된 평판에 진동을 가해주었을 때 액적 내부유동의 흐름이 변화가 발생하였고, 이는 곧 불균일한 침전물 패턴이 억제된 것을 증발 후 침전물의 패턴형상을 통해 확인할 수 있었다.

• 설비공학논문집
• /
• 제6권3호
• /
• pp.302-314
• /
• 1994

The effects of the inclination of enclosure and partition on natural convective flow and heat transfer were investigated numerically. The enclosure was composed of the lower hot and the upper cold horizontal walls and the adiabatic vertical walls, and a partition was positioned perpendicularly at the mid-height of one vertical insulated wall. The governing equations are solved by using the finite element method with Galerkin method. The computations were performed with the variations of the partition length and Rayleigh number based on the temperature difference between two horizontal walls and the enclosure height with water(Pr=4.95). The effects of the inclination angle of enclosure and partition on the heat transfer within an enclosure were also studied. As the results, the increase of the inclination angle of enclosure rapidly raised the heat transfer rate, while the inclination angle for the maximum Nusselt number was retarded with the increase of the partition length and the decrease of the heat transfer rate became larger in proportion to the increase of the partition length. The Nusselt number obtained by the inclination of partition was smaller than that of the inclination of enclosure. However, the difference of the heat transfer rates was considerably decreased at the longer partition lengths and the trends for the variation of the average Nusselt number were more similar with that of the inclination of enclosure. The upward oriented partition increases the convective heat transfer distinctly in contrast to that of the inclination of enclosure as the partition length increases.

• Nuclear Engineering and Technology
• /
• 제44권3호
• /
• pp.261-272
• /
• 2012

Molten Core Concrete Interaction (MCCI) is a complex process characterized by concrete ablation and volatile generation; Thermal and solutal convection in a bubble-agitated melt; Physico-chemical evolution of the corium pool with a wide solidification range (of the order of 1000 K). Twelve experiments have been carried out in the VULCANO facility with prototypic corium and sustained heating. The dry oxidic corium tests have contributed to show that silica-rich concrete experience an anisotropic ablation. This unexpected ablation pattern is quite reproducible and can be recalculated, provided an empirical anisotropy factor is assumed. Dry tests with oxide and metal liquid phases have also yielded unexpected results: a larger than expected steel oxidation and unexpected topology of the metallic phase (at the bottom of the cavity and also on the vertical concrete walls). Finally, VULCANO has proved its interest for the study of mitigation solutions such as the COMET bottom flooding core catcher.

• 대한기계학회논문집
• /
• 제19권12호
• /
• pp.3329-3343
• /
• 1995

The interaction of mixed convection and surface radiation in a printed circuit board(PCB) is investigated numerically. The electronic equipment is modeled by a two-dimensional channel with three hot blocks. In order to calculate the turbulent flow characteristics, the low Reynolds number k-.epsilon. model which is proposed by Launder and Sharma is applied. The S-4 approximation is used to solve the radiative transfer equation. The effects of the Reynolds number and geometric configuration variation of PCB on the flow and heat transfer characteristics are analyzed. As the results of this study, it is found that the thermal boundary layer occured at adiabatic wall in case with thermal radiation included, and the effect of radiation is also found to be insignificant for high Reynolds numbers. It is found, as well, that the heat transfer increases as the Reynolds number and block space increase and the channel height decreases and the heat transfer of vertical channel is greater than that of horizontal channel.

• 대한기계학회논문집B
• /
• 제30권8호
• /
• pp.764-771
• /
• 2006

Knowledge of ground thermal properties is most important for the proper design of large BHE(borehole heat exchanger) systems. Thermal response tests with mobile measurement devices were first introduced in Sweden and USA in 1995. Thermal response tests have so far been used primarily for in insitu determination of design data for BHE systems, but also for evaluation of grout material, heat exchanger types and ground water effects. The main purpose has been to determine insitu values of effective ground thermal conductivity, including the effect of ground-water flow and natural convection in the boreholes. Test rig is set up on a small trailer, and contains a circulation pump, a heater, temperature sensors and a data logger for recording the temperature data. A constant heat power is injected into the borehole through the pipe system of test rig and the resulting temperature change in the borehole is recorded. The recorded temperature data are analysed with a line-source model, which gives the effective insitu values of rock thermal conductivity and borehole thermal resistance.

• 설비공학논문집
• /
• 제9권3호
• /
• pp.312-322
• /
• 1997

The natural convective flow in a two-dimensional square enclosure with horizontal partitions is investigated numerically. The enclosure was composed of the lower hot and the upper cold horizontal walls and the adiabatic vertical walls, and two identical partitions were positioned perpendicularly at the mid-height of the right and left walls, respectively. The governing equations are solved by using the finite element method with Galerkin method. Calculations are made for different partition lengths, partition conductivites, and Rayleigh numbers based on the temperature difference between two horizontal walls and the enclosure height with water(Pr=4.95). An oscillatory motion of the natural convective flow is affected significantly by the variation of the gap width and Rayleigh number. When the gap width is comparatively short, the heat transfer rate is raised with the increase of the thermal conductivity of partitions. However, for sufficiently large gap widths at higher Rayleigh numbers, the average Nusselt numbers of the conductive partitions are smaller than those of the adiabatic partitions.

• 한국결정성장학회지
• /
• 제19권5호
• /
• pp.215-222
• /
• 2009

Many concepts of external magnetic field applications in crystal growth processes have been developed to control melt convection, impurity content and growing interface shape. Especially, travelling magnetic fields (TMF) are of certain advantages. However, strong shielding effects appear when the TMF coils are placed outside the growth vessel. To achieve a solution of industrial relevance within the framework of the $KRISTMAG^{(R)}$ project inner heater-magnet modules(HMM) for simultaneous generation of temperature and magnetic field have been developed. At the same time, as the temperature is controlled as usual, e.g. by DC, the characteristics of the magnetic field can be adjusted via frequency, phase shift of the alternating current (AC) and by changing the amplitude via the AC/DC ratio. Global modelling and dummy measurements were used to optimize and validate the HMM configuration and process parameters. GaAs and Ge single crystals with improved parameters were grown in HMM-equipped industrial liquid encapsulated Czochralski (LEC) puller and commercial vertical gradient freeze (VGF) furnace, respectively. The vapour pressure controlled Czochralski (VCz) variant without boric oxide encapsulation was used to study the movement of floating particles by the TMF-driven vortices.

• Journal of Radiation Protection and Research
• /
• 제45권3호
• /
• pp.95-100
• /
• 2020

Background: Compared with reported data of radon concentration, data of radon progeny concentration is limited in general, especially in outdoor environment. Materials and Methods: To know both the level and the variation of radon progeny concentration in outdoor environment in Beijing area, one-year continuous measurement with a cycle of 60 minutes was carried out by a step-advanced filter (SAF) monitor for radon progeny measurement. The observation site was located in a park in Eastern Beijing area, and the observation period was from October 17, 2018 to September 29, 2019. Results and Discussion: The equivalent equilibrium concentration (EEC) of radon progeny varies from 0.7 to 19.1 Bq·m^-3, with an annual average of 4.9 ± 2.7 Bq·m^-3. A clear diurnal variation of EEC, higher in the early morning and lower in the late afternoon, is observed due to the high sensitivity of the SAF monitor. Conclusion: Vertical convection of atmospheric boundary layer is thought to be the main reason of this phenomenon. For annual variation, the lowest monthly average EEC appeared in April, while the highest appeared in November, which might attribute to the atmospheric stability in different seasons.

• 대한기계학회논문집B
• /
• 제20권11호
• /
• pp.3589-3597
• /
• 1996

The Reynolds-averaged Navier-Stokes equations with the equations of the k-.epsilon. turbulence model are solved numerically in a general curvilinear system for a three-dimensional turbulent flow around an Ahmed body. The simulation is especially aimed at the evaluation of three finite differencing schemes for the convection term, which include the upwind differencing scheme(UDS), the second order upwind differencing scheme(SOU scheme) and the QUICK scheme. The drag coefficient, the velocity and pressure fields are found to be changed considerably with the adopted finite differencing schemes. It is clearly demonstrated that the large difference between computation and experiment in the drag coefficient is due to relatively high predicted values of pressure drag from both front part and vertical rear end base. The results also show that the simulation with the QUICK or SOU scheme predicts fairly well the flow field and gives more accurate drag coefficient than other finite differencing scheme.