• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1909-1914
• /
• 2004

An immersed boundary method for simulation of density-stratified flows is developed and applied to computation of viscous flows over two-dimensional obstacles in a bounded domain under stable density stratification. Density sources/sinks are introduced on the body surface. Two obstacle shapes are used, a vertical barrier and a smooth cosine-shaped hill; weak stratification, defined by $K=ND/{\pi}U{\leq}1$, where U, N, and D are the upstream velocity, buoyancy frequency, and domain height, respectively, is considered. The results are consistent with other authors' calculations, and shed light on computation of density-stratified flows in complex geometries.

• Nuclear Engineering and Technology
• /
• v.55 no.1
• /
• pp.180-189
• /
• 2023

We measured the heat load to a reactor vessel with and without the in-vessel debris bed under an IVR-ERVC condition. Mass transfer methodology was adopted based on heat and mass transfer analogy to achieve high Ra'_H of order ~10¹⁵ with compact test rigs. We postulated the in-vessel debris bed has a flat top and particulate debris was simulated as an identical diameter spheres. We conducted experiments varying the height of the debris bed and the results showed that Nusselt numbers decreased in both uppermost and curved surfaces with the increasing bed height. Once the debris bed is formed, it acts as an obstacle to the natural convective flow, which reduces the buoyancy. The reduction of driving force results in the impaired heat transfer in both upward and downward heat transfers.

• Nuclear Engineering and Technology
• /
• v.52 no.6
• /
• pp.1172-1179
• /
• 2020

As a Lagrangian particle method, Moving Particle Semi-implicit (MPS) method has great capability to capture interface/surface. In recent years, the multiphase flow simulation using MPS method has become one of the important directions of its developments. In this study, some key methods for multiphase flow have been introduced. The interface tension model in multiphase flow is modified to maintain the smooth of the interface and suitable for the three-phase flow. The mass transfer at immiscible liquid interface entrained by single bubble which could occur in Molten Core-Concrete Interaction (MCCI) has been investigated using this particle method. With the increase of bubble size, the height of entrainment column also increases, but the time of film rupture is slightly different. With the increase of density ratio between the two liquids, the height of entrained column decreases significantly due to the decreasing buoyancy of the denser liquid in the lighter liquid. In addition, the larger the interface tension coefficient is, the more rapidly the entrained denser liquid falls. This study validates that the MPS method has shown great performance for multiphase flow simulation. Besides, the influence of physical parameters on the mass transfer at immiscible interface has also been investigated in this study.

• Journal of Animal Environmental Science
• /
• v.18 no.1
• /
• pp.1-8
• /
• 2012

The study was conducted to investigate the effects of climate on indoor environment of a swine house with natural. This study was tested in the beef swine stall at Young-in, Kyung-ki do. The test was experimented for the effect of interior environment by the outdoor environment and the interior-pan. The results are as follows. 1. In test 1 ($T_{out}$ : $25.7^{\circ}C$, without fan), an indoor air flow pattern was showed that entered from sidewall winch-curtain to went out of a indoor by the ridge winch-curtain. And the velocity of a section of the center was measured two times as large as the velocity of the floor. It is the acceleration of the velocity by thermal buoyancy. And, the entered air was rapidly dissipated by flow energy. So that in the swain livestock with sidewall winch-curtain is effected by thermal buoyancy. And the air temperature of the indoor was distributed more higher as compared with the outdoor temperature. This result is caused by the sensible heat from swine and the ventilation is restricted. 2. In test 2 (($T_{out}$ : $25.7^{\circ}C$, with fan), the velocity of a section of the center was measured more higher as compared with the test 1. And the variance of air velocity was distributed higher as compared with the test 1. This result is showed dead region of air flow with a fan operation. And, the variance of gas density was distributed lower as compared with the test 1.

• Nuclear Engineering and Technology
• /
• v.54 no.7
• /
• pp.2728-2735
• /
• 2022

The research on the flow and heat transfer characteristics of lead bismuth(LBE) is significant for the thermal-hydraulic calculation, safety analysis and practical application of lead-based fast reactors(LFR). In this paper, a new CFD model is proposed to solve the thermal-hydraulic analysis of LBE. The model includes two parts: turbulent model and turbulent Prandtl, which are the important factors for LBE. In order to find the best model, the experiment data and design of 19-pin hexagonal rod bundle with spacer grid, undertaken at the Karlsruhe Liquid Metal Laboratory (KALLA) are used for CFD calculation. Furthermore, the turbulent model includes SST k - 𝜔 and k - 𝜀; the turbulent Prandtl includes Cheng-Tak and constant (Pr_t =1.5,2.0,2.5,3.0). Among them, the combination between SST k - 𝜔 and Cheng-Tak is more suitable for the experiment. But in the low Pe region, the deviation between the experiment data and CFD result is too much. The reason may be the inlet-effect and when Pe is in a low level, the number of molecular thermal diffusion occupies an absolute advantage, and the buoyancy will enhance. In order to test and verify versatility of the model, the NCCL performed by the Nuclear Thermal-hydraulic Laboratory (Nuthel) of Xi'an Jiao tong University is used for CFD to calculate. This paper provides two verification examples for the new universal model.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.36 no.2
• /
• pp.189-196
• /
• 2012

The heavy moderator acts as the ultimate heat-sink in an operating CANDU reactor. HUKINS has been developed to investigate moderator flow patterns. HUKINS consists of a 38.4-mm-thick cylindrical shell with a 0.95 m inner diameter and 88 sus-tubes that produce a total heat of 10 kW. A chemical visualization method was selected to estimate the occurrence of typical moderator flow patterns. Momentum-dominated flow, mixed flow, and buoyancy-dominated flow are detected under conditions of a heat load of 7.7 kW and input mass flow rates of 4, 7, and 11 L/min. The experimental results are similar to the results of a CFD simulation that consisted of approximately 1.9 million grids and was conducted using the k-${\varepsilon}$ turbulence model. Therefore, both the present experiments and simulations using HUKINS, a 1/8-scale model, represent all three important flow patterns expected in the real CANDU6 reference reactor. Thus, it has been demonstrated that HUKINS could be useful in the study of CANDU6 moderator circulation.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.42 no.2
• /
• pp.86-96
• /
• 2006

This research aims at establishing the fundamental characteristics of the kite through the analysis of the flow field around various types of kites. The approach of this study were adopted for the analysis; visualization by PIV(particle image velocimetry). Also, the lift and drag tests of kites had been performed in our previous finding(Bae et al., 2004a; Bae et al., 2004b). For this situation, models of canvas kite were deployed in the circulating water channel for the PIV test using the same conditions as in the lift and drag tests. The results obtained from the above approach are summarized as follows: Given the rectangular and triangular kites when attack angle is $20^{\circ}$, vortex by the boundary layer separation was seen in the leading edge and the flow towards the trailing edge was more turbulent. But, the inverted triangular type kite was seen to be stable without any boundary layer separation or turbulence. The increase of the attack angle resulted in the eddy in order of the rectangular, triangular and inverted triangular type. The magnitude of the eddy followed the same order. The effect of edge-eddy was biggest in the triangular type followed by the rectangular and then the inverted triangular type. The kite as the buoyancy device or the opening device will be very useful when the appropriate applications and the stability are met.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.3
• /
• pp.551-560
• /
• 1988

Gas-particle two phase flow and heat transfer in a cavity receiving thermal radiation through selectively transparent walls have been analyzed by a finite difference method. Particles injected from the upper hole of the cavity are accelerated downward by gravity and exit through the lower hole while they absorb, emit and scatter the incident thermal radiation. Gas phase is heated through convection heat transfer from particles, and consequently buoyancy induced flow field is formed. Two-equation model with two-way coupling is adopted and interaction terms are treated as sources by PSI-Cell method. For the particulate phase, Lagrangian method is employed to describe velocities and temperatures of particles. As thermal radiation is incident upon horizontally, radiative heat transfer in the vertical direction is assumed negligible and two-flux model is used for the solution of radiative heat flus. Gas phase velocity and temperature distributions, and particle trajectories, velocities and temperatures are presented. The effects of particle inlet condition, particle size, injection velocity and particle mass rate are mainly investigated.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.3
• /
• pp.702-715
• /
• 1990

This paper describes behaviors of two-phase open channel flow inside the flash chamber of a horizontal Multi-Stage-Flash evaporator numerically along with the experimental observations. Bubble trajectories and the velocity and temperature distributions of the liquid phase were predicted by using the particle-source-in-cell(PSI-Cell) method with the appropriate bubble motion/growth equations. Size and number of bubble nuclei embedded in the incoming liquid(brine) were taken into account as important parameters in addition to the conventional ones such as the velocity, degree of inlet superheat, inlet opening height, and the liquid level. Bubble motions, which are unsteady, appeared to be mostly determined by the buoyancy and the drag forces. The calculations, though a number of simplifying assumptions were made, reasonably simulated the hydrodynamic behaviors of the two-phase horizontal stream observed in the experiments. The simulated temperature distributions also agreed fairly well with the other's measurements. Non-equilibrium allownaces, evaluated from the simulated temperature distributions, were within the range of those obtained from the existing correlations, and reduced with the increases of the number and size of incoming bubble nuclei due to vigorous flashing.

• Nuclear Engineering and Technology
• /
• v.51 no.6
• /
• pp.1487-1503
• /
• 2019

The methods and performance of a pin-level nuclear reactor core thermal-hydraulics (T/H) code ESCOT employing the drift-flux model are presented. This code aims at providing an accurate yet fast core thermal-hydraulics solution capability to high-fidelity multiphysics core analysis systems targeting massively parallel computing platforms. The four equation drift-flux model is adopted for two-phase calculations, and numerical solutions are obtained by applying the Finite Volume Method (FVM) and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE)-like algorithm in a staggered grid system. Constitutive models involving turbulent mixing, pressure drop, and vapor generation are employed to simulate key phenomena in subchannel-scale analyses. ESCOT is parallelized by a domain decomposition scheme that involves both radial and axial decomposition to enable highly parallelized execution. The ESCOT solutions are validated through the applications to various experiments which include CNEN $4{\times}4$, Weiss et al. two assemblies, PNNL $2{\times}6$, RPI $2{\times}2$ air-water, and PSBT covering single/two-phase and unheated/heated conditions. The parameters of interest for validation include various flow characteristics such as turbulent mixing, spacer grid pressure drop, cross-flow, reverse flow, buoyancy effect, void drift, and bubble generation. For all the validation tests, ESCOT shows good agreements with measured data in the extent comparable to those of other subchannel-scale codes: COBRA-TF, MATRA and/or CUPID. The execution performance is examined with a mini-sized whole core consisting of 89 fuel assemblies and for an OPR1000 core. It turns out that it is about 1.5 times faster than a subchannel code based on the two-fluid three field model and the axial domain decomposition scheme works as well as the radial one yielding a steady-state solution for the OPR1000 core within 30 s with 104 processors.