• Nuclear Engineering and Technology
• /
• v.54 no.11
• /
• pp.4373-4391
• /
• 2022

Marine reactor systems experience platform movement, and therefore, the system thermal-hydraulic analysis code needs to reflect the motion effect on the fluid to evaluate reactor safety. A moving reactor model for MARS-KS was developed to simulate the hydrodynamic phenomena in the reactor under motion conditions; however, its applicability does not cover the MULTID component used in multidimensional flow analyses. In this study, a moving reactor model is implemented for the MULTID component to address the importance of multidimensional flow effects under dynamic motion. The concept of the volume connection is generalized to facilitate the handling of the junction of MULTID. Further, the accuracy in calculating the pressure head between volumes is enhanced to precisely evaluate the additional body force. Finally, the Coriolis force is modeled in the momentum equations in an acceleration form. The improvements are verified with conceptual problems; the modified model shows good agreement with the analytical solutions and the computational fluid dynamic (CFD) simulation results. Moreover, a simplified gravity-driven injection is simulated, and the model is validated against a ship flooding experiment. Throughout the verifications and validations, the model showed that the modification was well implemented to determine the capability of multidimensional flow analysis under ocean conditions.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.8
• /
• pp.1153-1163
• /
• 2002

A transient analysis on temperatures of fuel and oil in hydraulic and lubrication systems in an aircraft was studied using the finite difference method. Numerical calculation was performed by an explicit method with modified Dufort-Frankel scheme. Among various missions, air superiority mission was considered as a mission model with 20% hot day ambient condition in subsonic region. The ambience of the aircraft was assumed as turbulent flow. Convective heat transfer coefficient were used in calculating heat transfer between the aircraft surface and the ambience. For an aircraft on the ground, an empirical equation represented as a function of free-stream air velocity was used. And the heat transfer coefficient for flat plate turbulent flow suggested by Eckert was employed for in-flight phases. The governing equations used in this analysis are the mass and energy conservation equations on fuel and oils. Here, analysis of fuel and oil temperature in the engine was not carried out. As a result of this analysis, the ground operation phase has shown the highest temperature and the largest rate of temperature increase among overall mission phases. Also, it is shown that fuel flow rate through fuel/oil heat exchanger plays an important role in temperature change of fuel and oil. This analysis could be an important part of studies to ensure thermal stability of the aircraft and can be applicable to thermal design of the aircraft fuel system.

• Nuclear Engineering and Technology
• /
• v.41 no.7
• /
• pp.875-884
• /
• 2009

We present three nuclear/hydrogen-related R&D activities being performed at KAIST: air-ingressed LOCA analysis code development, gas turbine analysis tool development, and hydrogen-production system analysis model development. The ICE numerical technique widely used for the safety analysis of water-reactors is successfully implemented into GAMMA, with which we solve the basic equations for continuity, momentum conservation, energy conservation of the gas mixture, and mass conservation of 6 species (He, N2, O2, CO, CO2, and H2O). GAMMA has been extensively validated using data from 14 test facilities. We developed a tool to predict the characteristics of HTGR helium turbines based on the throughflow calculation with a Newton-Raphson method that overcomes the weakness of the conventional method based on the successive iteration scheme. It is found that the current method reaches stable and quick convergence even under the off-normal condition with the same degree of accuracy. The dynamic equations for the distillation column of HI process are described with 4 material components involved in the HI process: H2O, HI, I2, H2. For the HI process we improved the Neumann model based on the NRTL (Non-Random Two-Liquid) model. The improved Neumann model predicted a total pressure with 8.6% maximum relative deviation from the data and 2.5% mean relative deviation, and liquid-liquid-separation with 9.52% maximum relative deviation from the data.

• The KSFM Journal of Fluid Machinery
• /
• v.10 no.4
• /
• pp.9-14
• /
• 2007

The purpose of present work is to analyze the convective heat transfer downstream of mixing vane in subchannel of nuclear reactor with three-dimensional Navier-Stokes equations. SST model is selected as a turbulence closure by comparing the performances of two different turbulent closures. Three different shapes of mixing vane are tested. And, thermal-hydraulic performances of these vanes are discussed. The results show that twist of the vane improves the heat transfer performance far downstream of the vane.

• Nuclear Engineering and Technology
• /
• v.35 no.1
• /
• pp.45-54
• /
• 2003

A hyperbolic two-phase, two-fluid equation system developed in the previous work has been implemented in an existing nuclear safety analysis code, MARS. Although the implicit treatment of interfacial pressure force term introduced in momentum equation of the hyperbolic equation system is required to enhance the numerical stability, it is very difficult to implement in the code because it is not possible to maintain the existing numerical solution structure. As an alternative, two-step approach with stabilizer momentum equations has been selected. The results of a linear stability analysis by Von-Neumann method show the equivalent stability improvement with fully-implicit solution method. To illustrate the applicability, the new solution scheme has been implemented into the best-estimate thermal-hydraulic analysis code, MARS. This paper also includes the comparisons of the simulation results for the perturbation propagation and water faucet problems using both two-step method and the original solution scheme.

• Journal of computational fluids engineering
• /
• v.16 no.2
• /
• pp.30-48
• /
• 2011

A three-dimensional thermal-hydraulic code, CUPID, has been developed for the analysis of transient two-phase flows at component scale. The CUPID code adopts a two-fluid three-field model for two-phase flows. A semi-implicit two-step numerical method was developed to obtain numerical solutions on unstructured grids. This paper presents an overview of the CUPID code development and assessment strategy. The governing equations, physical models, numerical methods and their improvements, and the systematic verification and validation processes are discussed. The code couplings with a system code, MARS, and, a three-dimensional reactor kinetics code, MASTER, are also presented.

• Nuclear Engineering and Technology
• /
• v.26 no.1
• /
• pp.19-31
• /
• 1994

The objective of this study is to conduct a thermal-hydraulic analysis on the spent fuel pool and to evaluate a parametric effect for the thermal-hydraulic analysis of spent fuel pool. The selected parameters are the Reynolds Number and the gap flow through the oater gap between fuel cell and fuel bundle. The simplified flow network for a path of fuel cells is used to analyze the natural circulation phenomenon. In the flow network analysis, the pressure drop for each assembly from the entrance of the fuel rack to the exit of the fuel assembly is balanced by the driving head due to the density difference between the pool fluid and the average fluid in each spent fuel assembly. The governing equations ore developed using this relation. But, since the parameters(flow rate, pressure loss coefficient, decay heat, density)are coupled each other, iteration method is used to obtain the solution. For the analysis of the YGN 3&4 spent fuel rack, 12 channels are considered and the inputs such as decay heat and pressure loss coefficient are determined conservatively. The results show the thermal-hydraulic characteristics(void fraction, density, boiling height)of the YGN 3&4 spent fuel rack. There occurs small amount of boiling in the cells. Fuel cladding temperature is lower than 343.3$^{\circ}C$. The evaluation of parametric effect indicates that flow resistances by geometric effect are very sensitive to Reynolds number in the transition region and the gap flow is negligible because of the larger flow resistance in the gap flow path than in the fuel bundle.

• Transactions of the Korean hydrogen and new energy society
• /
• v.30 no.5
• /
• pp.448-454
• /
• 2019

This paper presents a high-efficiency design technique for developing the serialized models of a single-channel pump based on the diameter, flow rate and head as the main performance parameters. The variation in pump performance by changing of the single-channel pump geometry was predicted based on computational fluid dynamics (CFD). Numerical analysis was conducted by solving three-dimensional steady Reynolds-averaged Navier-Stokes equations with the shear stress transport (SST) turbulence model. The tendencies of the hydraulic performance depending on the pump geometry scale were analyzed with the fixed rotational speed. These performances were expressed and evaluated as the functionalization for designing the serialized models of a single-channel pump in this work.

• Nuclear Engineering and Technology
• /
• v.41 no.10
• /
• pp.1347-1360
• /
• 2009

A multi-dimensional component for the thermal-hydraulic system analysis code, MARS, was developed for a more realistic three-dimensional analysis of nuclear systems. A three-dimensional and two-fluid model for a two-phase flow in Cartesian and cylindrical coordinates was employed. The governing equations and physical constitutive relationships were extended from those of a one-dimensional version. The numerical solution method adopted a semi-implicit and finite-difference method based on a staggered-grid mesh and a donor-cell scheme. The relevant length scale was very coarse compared to commercial computational fluid dynamics tools. Thus a simple Prandtl's mixing length turbulence model was applied to interpret the turbulent induced momentum and energy diffusivity. Non drag interfacial forces were not considered as in the general nuclear system codes. Several conceptual cases with analytic solutions were chosen and analyzed to assess the fundamental terms. RPI air-water and UPTF 7 tests were simulated and compared to the experimental data. The simulation results for the RPI air-water two-phase flow experiment showed good agreement with the measured void fraction. The simulation results for the UPTF downcomer test 7 were compared to the experiment data and the results from other multi-dimensional system codes for the ECC delivery flow.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.11 no.1
• /
• pp.18-30
• /
• 1999

A heat transfer analysis for the two-dimensional (2-D) steady state using finite difference method (FDM) is performed to predict the thermal behavior of the primary first-wall (FW) system of fusion reactor under various geometric and thermo-hydraulic conditions, such as the beryllium (Be) armor thickness, pitch of cooling tube, and coolant velocity. The FW consists of authentic steel (type 316 stainless steel solution annealed) for cooling tubes, Cu for cooling tubes embedding material, and Be for a protective armor, based on the International Thermonuclear Experiment Reactor (ITER) report. The present 2-D analysis, the control volume discretized with hybrid grid (rectangular grid and polar grid) and Gauss-Seidel iteration method are adapted to solve the governing equations. In the present study, geometric and thermo-hydraulic parameters are optimized with consideration of several limitations. Consequently, it is suggested that the adequate pitch of cooling tube is 22-32mm, the beryllium armor thickness is 10-12mm, and that the coolant velocity is 4.5m/s-6m/s for $100^{\circ}C$ of inlet coolant temperature. The cooling tube should locate near beryllium armor. But, it would be better for locating the center of Cu wall, considering problems of material and manufacturing. Also, 2-D analysis neglecting the axial temperature distribution of cooling tube is appropriate, regarding the discretization error in axial direction.