• Proceedings of the Korean Nuclear Society Conference
• /
• 1998.05a
• /
• pp.149-154
• /
• 1998

A subcritical reactor driven by a linear proton accelerator has been considered as a nuclear waste incinerator at Korea Atomic Energy Research Institute(KAERI). Since the multiplication factor of a subcritical reactor is less than unity, to compensate exponentially decreasing fission neutrons from spallation reactions are essentially required for operating the reactor in its steady state. furthermore, the profile of accelerator beam currents is very important in controlling a subcritical reactor, because the reactor power varies in accordance of the profile of external neutrons. We have developed a code system to find numerical solutions of reactor kinetics equations, which are the simplest dynamic model for controlling reactors. In a due course of our previous numerical study of point kinetics equations for critical reactors, however, we learned that the same code system can be used in studying dynamic behavior of the subcritical reactor. Our major motivation of this paper is to investigate responses of subcritical reactors for small changes in thermal hydraulic parameters. Building a thermal hydraulic model for the subcritical reactor dynamics, we performed numerical simulations for dynamic responses of the reactor based on point kinetics equations with a source term. Linearizing a set of coupled differential equations for reactor responses, we focus our research interest on dynamic responses of the reactor to variations of the thermal hydraulic parameters in transient phases.

• Nuclear Engineering and Technology
• /
• v.45 no.2
• /
• pp.191-202
• /
• 2013

A numerical analysis of thermal stratification in the upper plenum of the MONJU fast breeder reactor was performed. Calculations were performed for a 1/6 simplified model of the MONJU reactor using the commercial code, CFX-13. To better resolve the geometrically complex upper core structure of the MONJU reactor, the porous media approach was adopted for the simulation. First, a steady state solution was obtained and the transient solutions were then obtained for the turbine trip test conducted in December 1995. The time dependent inlet conditions for the mass flow rate and temperature were provided by JAEA. Good agreement with the experimental data was observed for steady state solution. The numerical solution of the transient analysis shows the formation of thermal stratification within the upper plenum of the reactor vessel during the turbine trip test. The temporal variations of temperature were predicted accurately by the present method in the initial rapid coastdown period (~300 seconds). However, transient numerical solutions show a faster thermal mixing than that observed in the experiment after the initial coastdown period. A nearly homogenization of the temperature field in the upper plenum is predicted after about 900 seconds, which is a much shorter-term thermal stratification than the experimental data indicates. This discrepancy may be due to the shortcoming of the turbulence models available in the CFX-13 code for a natural convection flow with thermal stratification.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1997.05a
• /
• pp.102-107
• /
• 1997

A stiffness in a dynamical system is numerically studied to investigate a sensitivity of a reactor to the delayed neutron spectra with the Doppler feedback. To test numerical procedure, we adopted a case of a reactivity accident in a point reactor model. We found that the stiffness is sensitive to a reactivity insertion rate and the delayed neutron spectra in the Doppler feedback phase. Our numerical results show that global reactor characteristics are not very sensitive to the delayed neutron spectra even though their instantaneous ones are sensitive. We present the time evolution of each precursor group explicitly.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.19 no.4
• /
• pp.285-290
• /
• 2007

This study has been numerically conducted to investigate the removal efficiency of Volatile Organic Compound (VOC) in photocatalytic micro-reactor. This study has placed emphasis on the improvements of the working condition of photocatalytic micro-reactor. The micro-reactor consists of 19 microchannels with a rectangular cross-section. For the validation of the current numerical study, a computation has been carried out to simulate an existing experimental study on the cylindrical reactor, which has shown a good agreement. The degradation characteristics with different inlet concentrations and velocities have been obtained. The current results can be used for the design of advanced photocatalytic micro-reactor.

• Korean Chemical Engineering Research
• /
• v.59 no.1
• /
• pp.85-93
• /
• 2021

For accurate and reliable process design for phenol oxidation in a plug flow reactor with supercritical water, modeling can be very insightful. Here, the velocity and density distribution along the reactor have been predicted by a numerical model and variations of temperature and phenol mass fraction are calculated under various flow conditions. The numerical model shows that as we proceed along the length of the reactor the temperature falls from above 430 ℃ to approximately 380 ℃. This is because the generated heat from the exothermic reaction is less that the amount lost through the walls of the reactor. Also, along the length, the linear velocity falls to less than one-third of the initial value while the density more than doubles. This is due to the fall in temperature which results in higher density which in turn demands a lower velocity to satisfy the continuity equation. Having a higher oxygen concentration at the reactor inlet leads to much faster phenol destruction; this leads to lower capital costs (shorter reactor will be required); however, the operational expenditures will increase for supplying the needed oxygen. The phenol destruction depends heavily on the kinetic parameters and can be as high as 99.9%. Using different kinetic parameters is shown to significantly influence the predicted distributions inside the reactor and final phenol conversion. These results demonstrate the importance of selecting kinetic parameters carefully particularly when these predictions are used for reactor design.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.11 no.6
• /
• pp.28-35
• /
• 2012

This study is on SCR reactor, NOx reduction system in Marine that has been an issue nowadays. Especially design data was obtained by numerical on flow uniformity that is one of the design factor in SCR reactor. Also pressure drop on catalyst size inserted into SCR reactor was compared by experiment and numerical analysis. S/W, numerical analysis used for this study was confirmed that the result of numerical analysis used STAR-CCM+, common use CFD code, pressure drop on catalyst is not big different from the result of numerical analysis. In addition, degree of uniformity of liquid on SCR reactor was over 0.9. Whereas it was assured that degree of uniformity of liquid was changed depends on the shape of pipe at the entrance of SCR.

• Nuclear Engineering and Technology
• /
• v.56 no.3
• /
• pp.1002-1012
• /
• 2024

To remove insoluble fission products, which could possibly cause reactor instability and significantly reduce heat transfer efficiency from primary system of molten salt reactor, a helium bubbling method is employed into a passive molten salt fast reactor. In this regard, two-phase flow behavior of molten salt and helium bubbles was investigated experimentally because the helium bubbles highly affect the circulation performance of working fluid owing to an additional drag force. As the helium flow rate is controlled, the change of key thermal-hydraulic parameters was analyzed through a two-phase experiment. Simultaneously, to assess the applicability of numerical model for the analysis of two-phase flow behavior, the numerical calculation was performed using the OpenFOAM 9.0 code. The accuracy of the numerical analysis code was evaluated by comparing it with the experimental data. Generally, numerical results showed a good agreement with the experiment. However, at the high helium injection rates, the prediction capability for void fraction of helium bubbles was relatively low. This study suggests that the multiphaseEulerFoam solver in OpenFOAM code is effective for predicting the helium bubbling but there exists a room for further improvement by incorporating the appropriate drag flux model and the population balance equation.

• Journal of Hydrogen and New Energy
• /
• v.26 no.3
• /
• pp.260-270
• /
• 2015

One-dimensional packed bed reactor model accounting for interfacial and intra-particle gradients was developed and based on it numerical analyses were performed to investigate the dynamic behavior of a commercial scale methanation reactor. Methanation reaction was almost complete near the reactor inlet and gases with equilibrated composition were discharged from the reactor. Both the intra-particle temperature gradient and differential surface temperature rise were found to be severe near the reactor inlet. To reduce the possible degradation or fracture of catalyst particles and prevent local overheating on the catalyst, addition of inert material can be an effective way.

• Journal of computational fluids engineering
• /
• v.21 no.3
• /
• pp.8-14
• /
• 2016

A numerical analysis is conducted to estimate the core expansion and the energy behaviors induced by a core disruptive accident in a sodium-cooled fast reactor. The numerical formulation based on underwater explosion theory is carried out to simulate the core explosion inside the reactor vessel. The transient pressure, temperature and expansion of the core are examined by solving the equation of state and nonlinear governing equation of momentum conservation in one-dimensional spherical coordinates. The energy balance inside the computation domain is examined during the core expansion process. Heat transfer between the core and the sodium coolant, and the bubble rise during the expansion process are briefly investigated.

• Journal of Korean Society for Atmospheric Environment
• /
• v.23 no.2
• /
• pp.203-213
• /
• 2007

Numerical analysis has been performed to understand flow characteristics in the reactor with bag filters in an integrated adsorption/catalytic process which can treat dioxin and $NO_{x}$ together. Computational fluid dynamics technique was employed with Euler-Lagrangian model to consider flue gas and activated carbon particles simultaneously, so that residence time of flue gas and activated carbon particle could be obtained from the numerical analysis directly. The numerical analysis has been performed with different three particle sizes and compared each flow characteristics with particle's size. Fundamental flow patterns of flue gas and activated carbon particles, pressure distribution, residence time of flue gas and activated carbon particles, and distribution of activated carbon have been obtained from the numerical analysis. Flow patterns of flue gas and activated carbon particles in the reactor were very complicated and they moved along very various paths. Therefore, their residence time in the reactor was also various. The results obtained would be effectively used to estimate the removal efficiency in the reactor once the residence time is combined with the reaction equation.