• Nuclear Engineering and Technology
• /
• v.55 no.5
• /
• pp.1802-1813
• /
• 2023

Conjugate heat transfer between liquid metal and solid is a common phenomenon in a liquid-metal-cooled fast reactor's fuel assembly and heat exchanger, dramatically affecting the reactor's safety and economy. Therefore, comprehensively studying the sophisticated conjugate heat transfer in a liquid-metal-cooled fast reactor is profound. However, it has been evidenced that the traditional Simple Gradient Diffusion Hypothesis (SGDH), assuming a constant turbulent Prandtl number (Pr_t,, usually 0.85 - 1.0), is inappropriate in the Computational Fluid Dynamics (CFD) simulations of liquid metal. In recent decades, numerous studies have been performed on the four-equation model, which is expected to improve the precision of liquid metal's CFD simulations but has not been introduced into the conjugate heat transfer calculation between liquid metal and solid. Consequently, a four-equation model, consisting of the Abe k - ε turbulence model and the Manservisi k_𝜃 - ε_𝜃 heat transfer model, is applied to study the conjugate heat transfer concerning liquid metal in the present work. To verify the numerical validity of the four-equation model used in the conjugate heat transfer simulations, we reproduce Johnson's experiments of the liquid lead-bismuth-cooled turbulent pipe flow using the four-equation model and the traditional SGDH model. The simulation results obtained with different models are compared with the available experimental data, revealing that the relative errors of the local Nusselt number and mean heat transfer coefficient obtained with the four-equation model are considerably reduced compared with the SGDH model. Then, the thermal-hydraulic characteristics of liquid metal turbulent pipe flow obtained with the four-equation model are analyzed. Moreover, the impact of the turbulence model used in the four-equation model on overall simulation performance is investigated. At last, the effectiveness of the four-equation model in the CFD simulations of liquid sodium conjugate heat transfer is assessed. This paper mainly proves that it is feasible to use the four-equation model in the study of liquid metal conjugate heat transfer and provides a reference for the research of conjugate heat transfer in a liquid-metal-cooled fast reactor.

• Journal of the Korean Wood Science and Technology
• /
• v.17 no.3
• /
• pp.67-81
• /
• 1989

This experiment was carried out to know the mothod of changing the step of moisture content schedule with time in conventional kiln drying. For the purpose of this object. we made drying model by applying the moisture diffusion model by J.FSiau(1984) to average moisture content equation by J.Crank(1956) derived it from Fick's second law. And to verify this method of drying model. 2.5cm-thick boards and 5.0cm-thick dimension lumbers of Pinus densiflora were kiln-dried with the schedule of T11-C3 and T10-C4, respectively. And then the drying rates were investigated and compared with those calculated from drying model. The results obtained were as follows 1. Average drying rate and total drying time of board to dry to 6.5% moisture content were 0.64%/hr and 109hr., and those of dimension lumber to dry to 8.3% moisture content were 0.4%/hr. and 162hr., respectively. 2. The moisture content of shell and core decreased by equalizing treatment and increased by conditioning treatment both on board and dimension lumber. But the moisture gradient was lower after conditioning than after equalizing. 3. As the drying was proceeded, the transverse bound water diffusion coefficient all but linearly decreased, the water vapor diffusion coefficient abruptly curvilinearly increased, while the transverse diffusion coefficient curvilinearly decreased both on board and dimension lumber. But each of diffusion coefficients on board was larger than that on dimension lumber. 4. Compared to experimential drying rate of board. theoretical drying rate was larger at 30.0%-21.8% moisture content range and was similiar at 21.8%-5.4% moisture content. And in case of dimension lumber, the drying rate was similiar at 30.0%-16.1% moisture content range but theoretical drying rate was much lower at 16.1%-8.3% moisture content range. 5. The possibility of adapting this drying model to changing the moisture content schedule step with time was in the range of 21.8%-5.4% moisture content on board. And in the case of dimension lumber that was in the range of 30.0%-16.1% moisture content.

• Journal of Korea Water Resources Association
• /
• v.45 no.4
• /
• pp.349-360
• /
• 2012

A horizontally curvilinear non-hydrostatic free surface model that was applicable to three-dimensional viscous flows was developed. The proposed model employed a top-layer equation to close kinematic free-surface boundary condition, and an isotropic k-${\varepsilon}$ model to close turbulence viscosity in the Reynolds averaged Navier-Stokes equation. The model solved the governing equations with a fractional step method, which solved intermediate velocities in the advection-diffusion step, and corrects these provisional velocities by accounting for source terms including pressure gradient and gravity acceleration. Numerical applications were implemented to the wind-driven currents in a two-dimensional closed basin, the flow in a steep-sided trench, and the flow in a strongly-curved channel accounting for secondary current by the centrifugal force. Through the numerical simulations, the model showed its capability that were in good agreement with experimental data with respect to free surface elevation, velocity, and turbulence characteristics.

• Journal of Mechanical Science and Technology
• /
• v.17 no.3
• /
• pp.440-448
• /
• 2003

In the present paper, numerical simulations of buoyant plume dispersion in a neutral and stable atmospheric boundary layer have been carride out. A Lagrangian Stochastic Model (LSM) with a Non-Linear Eddy Viscosity Model (NLEVM) for turbulence is used to generate a Reynolds stress field as an input condition of dispersion simulation. A modified plume-rise equation is included in dispersion simulation in order to consider momentum effect in an initial stage of plume rise resulting in an improved prediction by comparing with the experimental data. The LSM is validated by comparing with the prediction of an Eulerian Dispersion Model (EDM) and by the measured results of vertical profiles of mean concentration in the downstream of an elevated source in an atmospheric boundary layer. The LSM predicts accurate results especially in the vicinity of the source where the EDM underestimates the peak concentration by 40% due to inherent limitations of gradient diffusion theory. As a verification study, the LSM simulation of buoyant plume dispersions under a neutral and stable atmospheric condition is compared with a wind-tunnel experiment, which shows good qualitative agreements.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.19 no.2
• /
• pp.205-214
• /
• 2021

Currently, the interim storage pools of spent fuels in South Korea are expected to become saturated from 2024. It is required to prepare an operation plan of a domestic dry storage facility during a long-term period, with the researches on safety evaluation methods. This study modified the FRAPCON code to predict the spent fuel integrity evaluation such as the axial cladding temperature, the hoop stress and hydrogen distribution in dry storage. The cladding temperature in dry storage was calculated using the COBRA-SFS code with the burnup information which was calculated using the FRAPCON code. The hoop stress was calculated using the ideal gas equation with spent fuel information such as rod internal pressure. Numerical analysis method was used to calculate the degree of hydrogen diffusion according to the hydrogen concentration and temperature distribution during a dry storage period. Before 50 years of dry storage, the cladding temperature and hoop stress decreased rapidly. However, after 50 years, they decreased gradually and the cladding temperature was below 400 K. The initial temperature distribution and hydrogen concentration showed a parabolic line, but hydrogen was transferred by the hydrogen concentration and temperature gradient over time.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.1B
• /
• pp.73-81
• /
• 2009

This paper presents a numerical investigation of the impact of the secondary currents on solute transport in open-channel flows. The RANS model with Reynolds stress model is used for flow modeling, and the GGDH(generalized gradient diffusion hypothesis) model is used to close the scalar transport equation. Using the developed model, the impact of secondary currents on solute transport in open channel flows over smooth-rough strip is investigated. Through numerical experiments, the secondary currents are found to affect the solute spreading, leading a movement of the position of the peak concentration and a skewed distribution of solute concentration. Due to the lateral flow of secondary currents near the free surface, the concentration at the rough strip is found to be larger than that at the smooth strip bed. The solute at the rough strip is more rapidly transported than smooth bed. A magnitude analysis of the solute transport rate in scalar transport equation is also carried out to investigate the effect of secondary currents and scalar flux on the concentration distribution.

• International Journal of Highway Engineering
• /
• v.7 no.1 s.23
• /
• pp.31-38
• /
• 2005

Concrete has higher vapor pressure than its surrounding ambient air immediately after placement. Moisture at concrete surface evaporates to the ambient air to adjust equilibrium of the vapor pressure between them. The moisture inside the concrete moves to the surface because the evaporation at the surface causes gradient of vapor pressure inside the concrete. Plastic cracking, degree of hydration, strength development, and others caused by velocity of the moisture movement significantly influences quality of concrete. In this paper, the moisture diffusivity of early-age concrete was back-calculated using governing equation of the moisture diffusion, and temperature and relative humidity of concrete measured in a laboratory. The moisture diffusivity of the concrete was modeled using the back-calculated moisture diffusivity. The relative humidity of the concrete calculated by finite element method (FEM) using the modeled moisture diffusivity as Input data coincided with the measured relative humidity well.

• Journal of Advanced Marine Engineering and Technology
• /
• v.35 no.5
• /
• pp.616-624
• /
• 2011

The flow characteristics and heat transfer augment on the periodically arranged semi-circular ribs in a rectangular channel for turbulent flow has been investigated numerically. The aspect ratio of the rectangular channel was AR=5, the rib height to hydraulic diameter ratio were 0.07 and rib height to channel height ratio was set as e/H=0.117 for various PR(rib pitch-to-rib height rate) between 8~14, respectively. The SST k-${\omega}$ turbulence model and v2-f turbulence model were used to find out the heat transfer and the flow characteristics of near the wall which are suited to obtain realistic phenomena. The numerical analysis results show turbulent flow characteristics, heat transfer enhancement and friction factor as observed experimentally. The results predict that turbulent kinetic energy(k) is closely relative to the diffusion of recirculation flow. and v2-f turbulence model simulation results have a good agreement with experimental values.

• Journal of the Korean Vacuum Society
• /
• v.18 no.1
• /
• pp.66-72
• /
• 2009

The detection methods are required to monitor and diagnose the abnormality on the insulation condition inside a gas-insulated switchgear (GIS). Due to a good sensitivity to the products decomposed by partial discharges (PDs) in $SF_6$ gas, the development of a SWNT gas sensor is actively in progress. However, a few numerical studies on the diffusion mechanism of the $SF_6$ decomposition products by PD have been reported. In this study, we modeled $SF_6$ decomposition process in a chamber by calculating temperature, pressure and concentration of the decomposition products by using a commercial CFD program in conjunction with experimental data. It was assumed that the mass production rate and the generation temperature of the decomposition products were $5.04{\times}10^{-10}$ [g/s] and over 773 K respectively. To calculate the concentration equation, the Schmidt number was specified to get the diffusion coefficient functioned by viscosity and density of $SF_6$ gas instead rather than setting it directly. The results showed that the drive potential is governed mainly by the gradient of the decomposition concentration. A lower concentration of the decomposition products was observed as the sensors were placed more away from the discharge region. Also, the concentration increased by increasing the discharge time. By installing multiple sensors the location of PD is expected to be identified by monitoring the response time of the sensors, and the information should be very useful for the diagnosis and maintenance of GIS.