• Nuclear Engineering and Technology
• /
• v.52 no.4
• /
• pp.721-733
• /
• 2020

In nuclear engineering, the occurrence of critical heat flux (CHF) is complicated for rod bundle, and it is much more difficult to predict the CHF when it is in natural circulation under motion condition. In this paper, the dryout-type CHF is investigated for the rod bundle in a natural circulation loop under rolling motion condition based on the coupled analysis of subchannel method, a one-dimensional system analysis method and a CHF mechanism model, namely the three-fluid model for annular flow. In order to consider the rolling effect of the natural circulation loop, the subchannel model is connected to the one-dimensional system code at the inlet and outlet of the rod bundle. The subchannel analysis provides the local thermal hydraulic parameters as input for the CHF mechanism model to calculate the occurrence of CHF. The rolling motion is modeled by additional motion forces in the momentum equation. First, the calculation methods of the natural circulation and CHF are validated by a published natural circulation experiment data and a CHF empirical correlation, respectively. Then, the CHF of the rod bundle in a natural circulation loop under both the stationary and rolling motion condition is predicted and analyzed. According to the calculation results, CHF under stationary condition is smaller than that under rolling motion condition. Besides, the CHF decreases with the increase of the rolling period and angular acceleration amplitude within the range of inlet subcooling and mass flux adopted in the current research. This paper can provide useful information for the prediction of CHF in natural circulation under motion condition, which is important for the nuclear reactor design improvement and safety analysis.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.19 no.4
• /
• pp.61-68
• /
• 2015

Reactive sprays of two propellant combinations(GOx/kerosene and $GN_2O$/ethanol) were observed and compared with each other as a basic research of visualizing supercritical combustion. A shadowgraph imaging method was used to visualize the reactive sprays, and shadowgraph images were converted to density gradient magnitude images to analyse the structure of reactive sprays. The gas-liquid interface of GOx/kerosene spray showed rougher boundary and steeper density gradient near the injector face than the $N_2O$/ethanol at similar combustion chamber pressure. Spray core length was calculated from averaged density gradient magnitude images and it was revealed that spray core length of GOx/kerosene was shorter than that of $GN_2O$/ethanol, although momentum flux ratio of GOx/kerosene propellant combination was lower.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.21 no.2
• /
• pp.41-47
• /
• 2017

By using liquid and gaseous nitrogen as simulants, effects of pintle opening distance and simulant supplying condition on spray from pintle injector were investigated qualitatively. The experiment was performed by fixing the liquid supply pressure and varying the opening of the pintle and the supply pressure of gaseous nitrogen. Shadowgraph method with CCD camera was used for image visualization, and the differences of liquid-gas spray were compared by image post processing. In case of liquid nitrogen single injection, as the opening distance was increased, the injection speed was decreased. Therefore, the end of liquid sheet were converged to the center of spray. In the case of liquid nitrogen/Gaseous nitrogen spray, although shadowgraph images look similar to each other, It is shown that it can be analysed by image post processing.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.1
• /
• pp.97-103
• /
• 2003

Numerical research has been done for the transverse jet behind a rearward-facings step in turbulent supersonic flow without chemical reaction. Purpose of transverse jet is to enhance mixing of the fuel in the combustor. Two-dimensional unsteady flowfields generated by slot injection into supersonic flow are numerically simulated with the Navier-Stokes equations with two-equation k-$\varepsilon$ turbulence model. Numerical method is used high-order upwind TVD scheme. Eight cases are computed for different slot momentum flux ratios and slot position at downstream of the step. The flow is very similar to the cavity flow, because the jet acts as an obstacle. The numerical results thus show the periodic phenomenon.

• Journal of Korea Water Resources Association
• /
• v.41 no.8
• /
• pp.761-772
• /
• 2008

The object of this study is to develop the model that solves the numerically difficult problems in hydraulic engineering and to demonstrate the applicability of this model by means of various test examples, such as, verification in the gradually varied unsteady condition, three steady flow problems with the change of bottom slope with exact solution, and frictional bed with analytical solution. The governing equation of this model is the integral form of the Saint-Venant equation satisfying the conservation laws, and finite volume method with the Riemann solver is used. The evaluation of the mass and momentum flux with the HLL Riemann approximate solver is executed. MUSCL-Hancock scheme is used to achieve the second order accuracy in space and time. This study introduce the new and simple technique to discretize the source terms of gravity and hydrostatic pressure force due to longitudinal width variation for the balance of quantity between nonlinear flux and source terms. The results show that the developed model's implementation is accurate, robust and highly stable in various flow conditions with source terms, and this model is reliable for one-dimensional applications in hydraulic engineering.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2014.02a
• /
• pp.221-221
• /
• 2014

The present work proposes an improved numerical simulator for design and modification of large area capacitively coupled plasma (CCP) processing chamber. CCP, as notoriously well-known, demands the tremendously huge computational cost for carrying out transient analyses in realistic multi-dimensional models, because electron dissociations take place in a much smaller time scale (${\Delta}t{\approx}10-8{\sim}10-10$) than time scale of those happened between neutrals (${\Delta}t{\approx}10-1{\sim}10-3$), due to the rf drive frequencies of external electric field. And also, for spatial discretization of electron flux (Je), exponential scheme such as Scharfetter-Gummel method needs to be used in order to alleviate the numerical stiffness and resolve exponential change of spatial distribution of electron temperature (Te) and electron number density (Ne) in the vicinity of electrodes. Due to such computational intractability, it is prohibited to simulate CCP deposition in a three-dimension within acceptable calculation runtimes (<24 h). Under the situation where process conditions require thickness non-uniformity below 5%, however, detailed flow features of reactive gases induced from three-dimensional geometric effects such as gas distribution through the perforated plates (showerhead) should be considered. Without considering plasma chemistry, we therefore simulated flow, temperature and species fields in three-dimensional geometry first, and then, based on that data, boundary conditions of two-dimensional plasma discharge model are set. In the particular case of SiH4-NH3-N2-He CCP discharge to produce deposition of SiNxHy thin film, a cylindrical showerhead electrode reactor was studied by numerical modeling of mass, momentum and energy transports for charged particles in an axi-symmetric geometry. By solving transport equations of electron and radicals simultaneously, we observed that the way how source gases are consumed in the non-isothermal flow field and such consequences on active species production were outlined as playing the leading parts in the processes. As an example of application of the model for the prediction of the deposited thickness uniformity in a 300 mm wafer plasma processing chamber, the results were compared with the experimentally measured deposition profiles along the radius of the wafer varying inter-electrode gap. The simulation results were in good agreement with experimental data.