• 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.

• Nuclear Engineering and Technology
• /
• v.54 no.1
• /
• pp.97-109
• /
• 2022

Multi-dimensional coupling analysis is a research hot spot in nuclear reactor thermal hydraulic study and both the full-scale system transient response and local key three-dimensional thermal hydraulic phenomenon could be obtained simultaneously, which can achieve the balance between efficiency and accuracy in the numerical simulation of nuclear reactor. A one-dimensional to three-dimensional (1D-3D) coupling platform for the nuclear reactor multi-dimensional analysis is developed by XJTU-NuTheL (Nuclear Thermal-hydraulic Laboratory at Xi'an Jiaotong University) based on the CFD code Fluent and system code RELAP5 through the Dynamic Link Library (DLL) technology and Fluent user-defined functions (UDF). In this paper, the International Standard Problem (ISP) No. 43 is selected as the benchmark and the rapid boron dilution transient in the nuclear reactor is studied with the coupling code. The code validation is conducted first and the numerical simulation results show good agreement with the experimental data. The three-dimensional flow and temperature fields in the downcomer are analyzed in detail during the transient scenarios. The strong reverse flow is observed beneath the inlet cold leg, causing the de-borated water slug to mainly diffuse in the circumferential direction. The deviations between the experimental data and the transients predicted by the coupling code are also discussed.

• Journal of Ocean Engineering and Technology
• /
• v.35 no.1
• /
• pp.50-58
• /
• 2021

The demand for eco-friendly energy is expected to increase due to the recently strengthened environmental regulations. In particular, the flow inside the pipe used in a cargo handling system (CHS) or fuel gas supply system (FGSS) of hydrogen transport ships and hydrogen-powered ships exhibits a very complex pattern of multiphase-thermal flow, including the boiling phenomenon and high accuracy analysis is required concerning safety. In this study, a feasibility study applying the boiling model was conducted to analyze the multiphase-thermal flow in the pipe considering the phase change. Two types of boiling models were employed and compared to implement the subcooled boiling phenomenon in nucleate boiling numerically. One was the "Rohsenow boiling model", which is the most commonly used one among the VOF (Volume-of-Fluid) boiling models under the Eulerian-Eulerian framework. The other was the "wall boiling model", which is suitable for nucleate boiling among the Eulerian multiphase models. Moreover, a comparative study was conducted by combining the nucleate site density and bubble departure diameter model that could influence the accuracy of the wall boiling model. A comparison of the Rohsenow boiling and the wall boiling models showed that the wall boiling model relatively well represented the process of bubble formation and development, even though more computation time was consumed. Among the combination of models used in the wall boiling model, the simulation results were affected significantly by the bubble departure diameter model, which had a very close relationship with the grid size. The present results are expected to provide useful information for identifying the characteristics of various parameters of the boiling model used in CFD simulations of multiphase-thermalflow, including phase change and selecting the appropriate parameters.

• Nuclear Engineering and Technology
• /
• v.53 no.3
• /
• pp.763-775
• /
• 2021

The thermal hydraulic performances of Steam Generator (SG) under both steady and transient operation conditions are of great importance for the safety and economy in nuclear power plants. In this paper, based on our self-developed SG thermal hydraulic analysis code STAF (Steam-generator Thermalhydraulic Analysis code based on Fluent), an improved new version STAF-CT (fully Coupling and Transient) is developed and introduced. Compared with original STAF, the new version code STAF-CT has two main functional improvements including "Transient" and "Fully Three Dimensional Coupling" features. In STAF-CT, a three dimensional energy transferring module is established which can achieve energy exchange computing function at the corresponding position between two sides of SG. The STAF-CT is validated against the international benchmark experiment data and the results show great agreement. Then the U-shaped SG in AP1000 nuclear power plant is modeled and simulated using STAF-CT. The results show that three dimensional flow fields in the primary side make significant effect on the energy source distribution between two sides. The development of code STAF-CT in this paper can provide an effective method for further SG high fidelity research in the nuclear reactor system.

• Nuclear Engineering and Technology
• /
• v.53 no.12
• /
• pp.3902-3909
• /
• 2021

SACOS series of subchannel analysis codes have been developed by XJTU-NuTheL for many years and are being used for the thermal-hydraulic safety analysis of various reactor cores. To achieve fine whole core pin-level analysis, the input preprocessing and parallel capabilities of the code have been developed in this study. Preprocessing is suitable for modeling rectangular and hexagonal assemblies with less error-prone input; parallelization is established based on the domain decomposition method with the hybrid of MPI and OpenMP. For domain decomposition, a more flexible method has been proposed which can determine the appropriate task division of the core domain according to the number of processors of the server. By performing the calculation time evaluation for the several PWR assembly problems, the code parallelization has been successfully verified with different number of processors. Subsequent analysis results for rectangular- and hexagonal-assembly core imply that the code can be used to model and perform pin-level core safety analysis with acceptable computational efficiency.

• Nuclear Engineering and Technology
• /
• v.53 no.12
• /
• pp.3910-3917
• /
• 2021

Reactor fuel's power distribution plays a vital role in designing the new generation thermionic Space Reactor Power Systems (SRPS). In this paper, the 1/12th SPACE-R's full reactor core was numerically analyzed with two kinds of different axial power distribution, to identify their impacts on thermal-hydraulic and thermoelectric characteristics. In the benchmark study, the maximum error between numerical results and existing data or design values ranged from 0.2 to 2.2%. Four main conclusions were obtained in the numerical analysis: a) The axial power distribution has less impact on coolant temperature. b) Axial power distribution influenced the emitter temperature distribution a lot, when the core power was cosine distributed, the maximum temperature of the emitter was 194 K higher than that of the uniform power distribution. c) Comparing to the cosine axial power distribution, the uniform axial power distribution would make the maximum temperature in each component of the reactor core much lower, reducing the requirements for core fuel material. d) Voltage and current distribution were similar to the axial electrode temperature distribution, and the axial power distribution has little effect on the output power.

• Nuclear Engineering and Technology
• /
• v.52 no.1
• /
• pp.59-68
• /
• 2020

ECCMIX component was introduced in RELAP5/MOD3 for calculating the interfacial condensation. Compared to other existing components in RELAP5, user experience of ECCMIX component is restricted to developmental assessment applications. To evaluate the capability of the ECCMIX component, ECCS experiment was conducted which included single-phase and two-phase thermal mixing. The experiment was carried out with test sections containing a main pipe (70 mm inner diameter) and a branch pipe (21 mm inner diameter) under the atmospheric pressure. The steam mass flow in the main pipe ranged from 0 to 0.0347 kg/s, and the subcooled water mass flow in the branch pipe ranged from 0.0278 to 0.1389 kg/s. The comparison of the experimental data with the calculation results illuminated that although the ECCMIX component was more difficult to converge than Branch component, it was a more appropriate manner to simulate interfacial condensation under two-phase thermal mixing circumstance, while the two components had no differences under single-phase circumstance.

• 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.

• 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.

• Nuclear Engineering and Technology
• /
• v.53 no.5
• /
• pp.1540-1555
• /
• 2021

The plate-type fuel assembly adopted in nuclear research reactor suffers from complicated effect induced by non-uniform irradiation, which might affect its stress conditions, mechanical behavior and thermal-hydraulic performance. A reliable numerical method is of great importance to reveal the complex evolution of mechanical deformation, flow redistribution and temperature field for the plate-type fuel assembly under non-uniform irradiation. This paper is the first part of a two-part study developing the numerical methodology for the thermal-fluid-structure coupling behaviors of plate-type fuel assembly under irradiation. In this paper, the thermal-fluid-structure coupling methodology has been developed for plate-type fuel assembly under non-uniform irradiation condition by exchanging thermal-hydraulic and mechanical deformation parameters between Finite Element Model (FEM) software and Computational Fluid Dynamic (CFD) software with Mesh-based parallel Code Coupling Interface (MpCCI), which has been validated with experimental results. Based on the established methodology, the effects of non-uniform irradiation and fluid were discussed, which demonstrated that the maximum mechanical deformation with irradiation was dozens of times larger than that without irradiation and the hydraulic load on fuel plates due to differential pressure played a dominant role in the mechanical deformation.