• Nuclear Engineering and Technology
• /
• v.53 no.6
• /
• pp.1834-1845
• /
• 2021

Liquid lead-bismuth cooled fast reactor is one of the most promising reactor types among the fourth-generation nuclear energy systems. The flow and heat transfer characteristics of lead-bismuth eutectic (LBE) are completely different from ordinary fluids due to its special thermal properties, causing that the traditional Reynolds analogy is no longer recommended and appropriate. More accurate turbulence flow and heat transfer model for the liquid metal lead-bismuth should be developed and applied in CFD simulation. In this paper, a specific CFD solver for simulating the flow and heat transfer of liquid lead-bismuth based on the k - 𝜀 - k_𝜃 - 𝜀_𝜃 model was developed based on the open source platform OpenFOAM. Then the advantage of proposed model was demonstrated and validated against a set of experimental data. Finally, the simulation of LBE turbulent flow and heat transfer in a 7-pin wire-wrapped rod bundle with the k - 𝜀 - k_𝜃 - 𝜀_𝜃 model was carried out. The influence of wire on the flow and heat transfer characteristics and the three-dimensional distribution of key thermal hydraulic parameters such as temperature, cross-flow velocity and Nusselt number were studied and presented. Compared with the traditional SED model with a constant Pr_t = 1.5 or 2.0, the k - 𝜀 - k_𝜃 - 𝜀_𝜃 model is more accurate on predicting the turbulence flow and heat transfer of liquid lead-bismuth. The average relative error of the k - 𝜀 - k_𝜃 - 𝜀_𝜃 model is reduced by 11.1% at most under the simulation conditions in this paper. This work is meaningful for the thermal hydraulic analysis and structure design of fuel assembly in the liquid lead-bismuth cooled fast reactor.

• Journal of computational fluids engineering
• /
• v.16 no.2
• /
• pp.94-101
• /
• 2011

The fast pyrolysis characteristics of lignocellulosic biomass are investigated for a bubbling fluidized bed reactor by means of computational fluid dynamics (CFD). To simulate multiphase reacting flows for gases and solids, an Eulerian-Eulerian approach is applied. Attention is paid for the primary and secondary reactions affected by gas-solid flow field. From the result, it is scrutinized that fast pyrolysis reaction is promoted by chaotic bubbling motion of the multiphase flow enhancing the mixing of solid particles. In particular, vortical flow motions around gas bubbles play an important role for solid mixing and consequent fast pyrolysis reaction. Discussion is made for the time-averaged pyrolysis reaction rates together with time-averaged flow quantities which show peculiar characteristics according to local transverse location in a bubbling fluidized bed reactor.

• Journal of computational fluids engineering
• /
• v.15 no.4
• /
• pp.99-108
• /
• 2010

Fractional Step Methods(FSM) are popular in simulation of unsteady incompressible flow. In this study, we demonstrate that FSM, combined with a Volume-Of-Fluid method, can be further applied to simulation of multiphase flow. The interface between the fluids is constructed by the effective least squares volume-of-fluid interface reconstruction algorithm and advected by the velocity using the operator split advection algorithm. To verify our numerical methodology, our results are compared with other authors' numerical and experimental results for the benchmark problems, revealing excellent agreement. The present FSM sheds light on accurate simulation of turbulent multiphase flow which is found in many engineering applications.

• Nuclear Engineering and Technology
• /
• v.52 no.10
• /
• pp.2196-2203
• /
• 2020

A dryout mechanism model for rectangular narrow channels at high pressure conditions is developed by assuming that the Kelvin-Helmholtz instability triggered the occurrence of dryout. This model combines the advantages of theoretical analysis and empirical correlation. The unknown coefficients in the theoretical derivation are supported by the experimental data. Meanwhile, the decisive restriction of the experimental conditions on the applicability of the empirical correlation is avoided. The expression of vapor phase velocity at the time of dryout is derived, and the empirical correlation of liquid film thickness is introduced. Since the CHF value obtained from the liquid film thickness should be the same as the value obtained from the Kelvin-Helmholtz critical stability under the same condition, the convergent CHF value is obtained by iteratively calculating. Comparing with the experimental data under the pressure of 6.89-13.79 MPa, the average error of the model is -15.4% with the 95% confidence interval [-20.5%, -10.4%]. And the pressure has a decisive influence on the prediction accuracy of this model. Compared with the existing dryout code, the calculation speed of this model is faster, and the calculation accuracy is improved. This model, with great portability, could be applied to different objects and working conditions by changing the expression of the vapor phase velocity when the dryout phenomenon is triggered and the calculation formula of the liquid film.

• Nuclear Engineering and Technology
• /
• v.54 no.8
• /
• pp.3117-3129
• /
• 2022

Deeply understanding the phase change of thin liquid sodium film inside wick pore is very important for further studying high-temperature sodium heat pipe's heat transfer. For the first time, the evaporation and condensation of thin liquid sodium film are investigated by the Lennard-Jones potential of molecular dynamics. Based on the startup and normal operation of the sodium heat pipe, three different cases are simulated. First, the equilibrium is achieved and the Mass Accommodation Coefficients of the three cases are 0.3886, 0.2119, 0.2615 respectively. Secondly, the non-equilibrium is built. The change of liquid film thickness, the number of gas atoms, the net evaporation flux (J_net), the heat transfer coefficient (h) at the liquid-gas interface are acquired. Results indicate that the magnitude of the J_net and the h increase with the basic equilibrium temperature. In 520-600 K (the startup of the heat pipe), the h has approached 5-6 W m^-2 K^-1 while liquid film thickness is in 11-13 nm. The fact shows that during the initial startup of the sodium heat pipe, the thermal resistance at the liquid-gas interface can't be negligible. This work is the complement and extension for macroscopic investigation of heat transfer inside sodium heat pipe. It can provide a reference for further numerical simulation and optimal design of the sodium heat pipe in the future.

• Journal of computational fluids engineering
• /
• v.9 no.1
• /
• pp.18-24
• /
• 2004

The unified simulation for the multiphase flow by predictor-corrector scheme based on CIP method is introduced. In this algorithm, the interface between different phases is identified by a density function and tracked by solving an advection equation. Solid body motion is modeled by the translation and angular motion. The mathematical formulation and numerical results are also described. To verify the efficiency, accuracy and capability of proposed algorithm, two dimensional incompressible cavity flow, the motion of a floating ball into water and a single rising bubble by buoyancy force are numerically simulated by the present scheme. As results, it is confirmed that the present scheme gives an efficient, stable and reasonable solution in the multiphase flow problem.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.17 no.3
• /
• pp.209-219
• /
• 2013

The present paper deals with the extension of AUSMPW+ scheme into two-fluid model for multiphase flow. AUSMPW+ scheme is the improvement of a single-phase AUSM+ scheme by designing pressure-based weighting functions to prevent oscillations near a wall and shock instability after a strong shock. Recently, Kitamura and Liou assessed a family of AUSM-type schemes with two-fluid model governing equations [K. Kitamura and M.-S. Liou, Comparative study of AUSM-Family schemes in compressible multi-phase flow simulations, ICCFD7-3702 (2012)]. It was observed that the direct application of the single-phase AUSMPW+ did not provide satisfactory results for most of numerical test cases, which motivates the current study. It turns out that, by designing pressure-based weighting functions, which play a key role in controlling numerical diffusion for two-fluid model, problems reported in can be overcome. Various numerical experiments validate the proposed modification of AUSMPW+ scheme is accurate and robust to solve multiphase flow within the framework of two-fluid model.

• Nuclear Engineering and Technology
• /
• v.52 no.5
• /
• pp.937-946
• /
• 2020

The aim in the present work is to simulate accident scenarios of AP1000 during the small-break loss-of-coolant accident (SBLOCA) and investigate the performance and behavior of automatic depressurization system (ADS) during accidents by using MIDAC (The Module In-vessel Degradation severe accident Analysis Code). Four types of accidents with different hypothetical conditions were analyzed in this study. The impact on the thermal-hydraulic of the reactor coolant system (RCS), the passive core cooling system and core degradation was researched by comparing these types. The results show that the RCS depressurization becomes faster, the core makeup tanks (CMT) and accumulators (ACC) are activated earlier and the effect of gravity water injection is more obvious along with more ADS valves open. The open of the only ADS1-3 can't stop the core degradation on the basis of the first type of the accident. The open of ADS1-3 has a great impact on the injection time of ACC and CMT. The core can remain intact for a long time and the core degradation can be prevent by the open of ADS-4. The all results are significant and meaningful to understand the performance and behavior of the ADS during the typical SBLOCA.

• Particle and aerosol research
• /
• v.12 no.1
• /
• pp.21-26
• /
• 2016

In order to realize In-line and convenient measurement for solid-gas two phase flows, Light Transmission Fluctuation (LTF) based on the random variation of transmitted light intensity, light scattering theory and cross-correlation method was presented for online measurement of particle size, concentration and velocity. The statistical relationship among transmitted light intensity, particle size and particle number in measurement zone was described by Beer-Lambert Law. Accordingly, the particle size and concentration were determined from the fluctuation signal of transmitted light intensity. Simultaneously, the particle velocity was calculated by cross-correlation analysis of two neighboring light beams. By considering the influence of concentration variation in industrial applications, the improved algorithm based on spectral analysis of transmitted light intensity was proposed to improve measurement accuracy and stability. Therefore, the online measurement system based on LTF was developed and applied to measure pulverized coal in power station and raw material in cement plant. The particle size, concentration and velocity of powder were monitored in real-time. It can provide important references for optimal control, energy saving and emission reduction of energy-intensive industries.

• Proceedings of the Korea Water Resources Association Conference
• /
• 1995.05a
• /
• pp.238-244
• /
• 1995