• Journal of Mechanical Science and Technology
• /
• v.18 no.6
• /
• pp.990-1001
• /
• 2004

The retention of a molten pool vessel cooled by internal vessel reflooding and/or external vessel reactor cavity flooding has been considered as one of severe accident management strategies. The present numerical study investigates the effect of both internal and external vessel mixed cooling on an internally heated molten pool. The molten pool is confined in a hemispherical vessel with reference to the thermal behavior of the vessel wall. In this study, our numerical model used a scaled-down reactor vessel of a KSNP (Korea Standard Nuclear Power) reactor design of 1000 MWe (a Pressurized Water Reactor with a large and dry containment). Well-known temperature-dependent boiling heat transfer curves are applied to the internal and external vessel cooling boundaries. Radiative heat transfer has been considered in the case of dry internal vessel boundary condition. Computational results show that the external cooling vessel boundary conditions have better effectiveness than internal vessel cooling in the retention of the melt pool vessel failure.

• Journal of the Korean Society of Visualization
• /
• v.17 no.3
• /
• pp.24-31
• /
• 2019

Recently, an air-lift bio-reactor operated by micro bubbles has been utilized to product hydrogen fuel. To enhance the performance, characteristics of hydrodynamics inside the bio-reactor were analyzed using a numerical simulation for two-phase flow. An Eulerian model was employed for both of liquid and gas phases. The standard k-ε model was used for turbulence induced by micro bubbles. A Population Balance Model was employed to consider size distribution of bubbles. A hollow cylinder was introduced at the center of the reactor to reduce a dead area which disturbs circulation of CO bubbles. An appropriate diameter of the draft tube and hollow cylinder were optimized for better performance of the bio-reactor. The optimum model could be obtained when the cross-sectional area ratio of the hollow cylinder to the reactor, and the width ratio of the riser to the downcomer approached 0.4 and 3.5, respectively. Consequently, it is expected that the optimum model could enhance the performance of the bio-reactor with the homogeneous distribution and higher density of CO, and more effective mixing.

• Proceedings of the KSME Conference
• /
• 2000.11b
• /
• pp.228-233
• /
• 2000

The Korean Next Generation Reactor(KNGR) is a Pressurized Water Reactor adopting direct vessel injection(DVI) to optimize the performance of emergency core cooling system(ECCS). In a certain accident, however, pressurized thermal shock(PTS) of the vessel due to the sudden contact with the injected cold water is expected. In this paper, an accident of Main Steam Line Break(MSLB) has been numerically investigated with direct vessel injections and an increased volume flow rate in some cold legs. Using FLUENT code, temperature distributions of the fluid in the downcomer and of reactor vessel including the core region have been calculated, together with the distribution of convective heat transfer coefficient(CHTC) at the cladding surface of the reactor vessel. The result shows that some parts of the core region of the reactor vessel have higher temperature gradient expressing higher thermal stress.

• Nuclear Engineering and Technology
• /
• v.49 no.5
• /
• pp.1095-1099
• /
• 2017

A fission source can act as a stabilization element in coupled Monte Carlo simulations. We have observed this while studying numerical instabilities in nonlinear steady-state simulations performed by a Monte Carlo criticality solver that is coupled to a xenon feedback solver via fixed-point iteration. While fixed-point iteration is known to be numerically unstable for some problems, resulting in large spatial oscillations of the neutron flux distribution, we show that it is possible to stabilize it by reducing the number of Monte Carlo criticality cycles simulated within each iteration step. While global convergence is ensured, development of any possible numerical instability is prevented by not allowing the fission source to converge fully within a single iteration step, which is achieved by setting a small number of criticality cycles per iteration step. Moreover, under these conditions, the fission source may converge even faster than in criticality calculations with no feedback, as we demonstrate in our numerical test simulations.

• 제어로봇시스템학회:학술대회논문집
• /
• 1993.10b
• /
• pp.312-315
• /
• 1993

Optimal control theories based on the maximum principles have been evolved and applied to distributed parameter systems(DPSs) represented by partial differential equations (PDEs) and integral equations (IEs). This paper intends to show that an optimal control of a tubular reactor described by a one-dimensional partial differential equation was obtained using the distributed parameter control method for parabolic PDEs. In develping an algorithm which implements the calculation, the method of lines (MOL) was adopted through using a package called the DSS/2. For the tubular reactor system chosen for this paper, the optimal control method based on PDEs with the numerical MOL showed to be more efficient than the one based on IEs.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.781-785
• /
• 2003

This work presents a numerical analysis of two-phase natural circulation flow in reactor cavity under external vessel cooling. Steady, incompressible, three-dimensional Reynolds-averaged Navier-Stokes equations for multiphase flows with zero equation turbulence model are solved to predict the shear key effect on the circulation rate of cooling water and the distribution of void fraction according to the different mass flow of inlet air. Results show that shear key has a positive effect on the circulation rate of cooling water and induce a local increase of void fraction below the shear key, but not remarkably.

• Journal of the Korean Ceramic Society
• /
• v.32 no.11
• /
• pp.1246-1254
• /
• 1995

A numerical model for fabrication and deposition of ultrafine SiO2 particles were proposed in the simplified horizontal MCVD apparatus using tube furnace reactor. The model equations such as energy and mass balance equations and the 0th, 1st and 2nd moment balance equations of aerosols were considered in the reactor. The phenomena of SiCl4 chemical reaction, SiO2 particle formation and coagulation, diffusion and thermophoresis of SiO2 particles were included in the aerosol dynamic equation. The profiles of gas temperature, SiCl4 concentration and SiO2 particle volume were calculated for standard conditions. The concentrations, sizes and deposition efficiencies of SiO2 particles were calculated, changing the process conditions such as tube furnace setting temperature, total gas flow rate and inlet SiCl4 concentration.

• Nuclear Engineering and Technology
• /
• v.54 no.6
• /
• pp.2120-2134
• /
• 2022

A hypothetical core destructive accident (HCDA) has received widespread attention as one of the most serious accidents in sodium-cooled fast reactors. This study combined recent advantages in numerical methods to realize realistic modeling of the complex fluid-structure interactions during HCDAs in a full-scale sodium-cooled fast reactor. The multi-material arbitrary Lagrangian-Eulerian method is used to describe the fluid-structure interactions inside the container. Both the structural deformations and plug rises occurring during HCDAs are evaluated. Two levels of expansion energy are considered with two different reactor models. The simulation results show that the container remains intact during an accident with small deformations. The plug on the top of the container rises to an acceptable level after the sealing between the it and its support is destroyed. The methodology established in this study provides a reliable approach for evaluating the safety feature of a container design.

• Nuclear Engineering and Technology
• /
• v.19 no.4
• /
• pp.310-316
• /
• 1987

A successive iteration method to calculate the λ-modes of the diffusion equation was developed. The 2-group, 3-dimensional computer code MOGEN was developed to implement this method, The accuracy of the method was demonstrated using 2-dimensional bare homogeneous rectangular reactor. The numerical solution shows good agreement with the analytic solution in terms of eigenvalue and eigenfunction As for the standard CANDU-600 reactor, the 2-dimensional modes were generated and these represent the conventional mode characteristics well. Finally, application of theλ-mode in reactor engineering problems is described briefly.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.9
• /
• pp.709-717
• /
• 2009

Heat transfer rate is a very important factor for the performance of a steam reformer because a steam reforming reaction is an endothermic reaction. Coaxial cylindrical reactor is the reactor design which can improve the heat transfer rate. Temperature, fuel conversion and heat flux in the coaxial cylindrical steam reformer are studied in this paper using numerical method under various operating conditions. Langmuir-Hinshelwood model and pseudo-homogeneous model are incorporated for the catalytic surface reaction. Dominant chemical reactions are assumed as a Steam Reforming (SR) reaction, a Water-Gas Shift (WGS) reaction, and a Direct Steam Reforming (DSR) reaction. Although coaxial cylindrical steam reformer uses 33% less amount of catalyst than cylindrical steam reformer, its fuel conversion is increased 10 % more and its temperature is also high as about 30 degree. There is no heat transfer limitation near the inlet area at coaxial-type reactor. However, pressure drop of the coaxial cylindrical reactor is 10 times higher than that of cylindrical reactor. Operating parameters of coaxial cylindrical steam reformer are the wall temperature, the inlet temperature, and the Gas Hourly Space Velocity (GHSV). When the wall temperature is high, the temperature and the fuel conversion are increased due to the high heat transfer rate. The fuel conversion rate is increased with the high inlet temperature. However, temperature drop clearly occurs near the inlet area since an endothermic reaction is active due to the high inlet temperature. When GHSV is increased, the fuel conversion is decreased because of the heat transfer limitation and short residence time.