• Nuclear Engineering and Technology
• /
• v.54 no.2
• /
• pp.579-588
• /
• 2022

In order to meet the needs of domestic reactor severe accident analysis program, a MIDAC (Module Invessel Degraded severe accident Analysis Code) is developed and maintained by Xi'an Jiaotong University. As the accuracy of the calculation results of the analysis program is of great significance for the formulation of severe accident mitigation measures, the article select three experiments to evaluate the updated severe accident models of MIDAC. Among them, QUENCH-06 is the international standard No.45, QUENCH-16 is a test for the analysis of air oxidation, and FROMA is an out-of-pile fuel rod melting experiment recently carried out by Xi'an Jiaotong University. The heating and melting model with lumped parameter method and the steam oxidation model with Cathcart-Pawel and Volchek-Zvonarev correlations combination in MIDAC could better meet the needs of severe accident analysis. Although the influence of nitrogen still need to be further improved, the air oxidation model with NUREG still has the ability to provide guiding significance for engineering practice.

• Smart Structures and Systems
• /
• v.22 no.3
• /
• pp.249-257
• /
• 2018

This study contributes to evaluate multiphase topology optimization design of plate-like elastic structures with constant thickness and Reissner-Mindlin plate theory. Stiffness and adjoint sensitivity formulations linked to Reissner-Mindlin plate potential energy of bending and shear are derived in terms of multiphase design variables. Multiphase optimization problem is solved through alternative active-phase algorithm with Gauss-Seidel version as an optimization model of optimality criteria. Numerical examples verify efficiency and diversity of the present topology optimization method of Reissner-Mindlin elastic plates depending on multiphase and Poisson's ratio.

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

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.115-121
• /
• 2010

Cavitating flow is widely shown in many engineering systems, such as marine propellers, pump impellers, nozzles, injectors, torpedoes, etc. The present work focuses on the numerical analysis of the multiphase flow around the underwater vehicle which was launched from a submarine. The governing equation is the Navier-Stokes equation with a homogeneous mixture mode. The multiphase flow solver uses an implicit preconditioning scheme in curvilinear coordinate. For the code validation, the results from the present work are compared with the existing experimental and numerical results, and a reasonably good agrements are obtained. The multiphase flow around an underwater vehicle is simulated which includes submarine wake effects.

• Steel and Composite Structures
• /
• v.27 no.2
• /
• pp.177-184
• /
• 2018

This study contributes to evaluate multiphase topology optimization design of plate-like structures on elastic foundations by using classic plate theory. Multi-material optimal topology and shape are produced as an alternative to provide reasonable material assignments based on stress distributions. Multi-material topology optimization problem is solved through an alternative active-phase algorithm with Gauss-Seidel version as an optimization model of optimality criteria. Stiffness and adjoint sensitivity formulations linked to thin plate potential strain energy are derived in terms of multiphase design variables and Winkler-Pasternak parameters considering elastic foundation to apply to the current topology optimization. Numerical examples verify efficiency and diversity of the present topology optimization method of elastic thin plates depending on multiple materials and Winkler-Pasternak parameters with the same amount of volume fraction and total structural volume.

• Nuclear Engineering and Technology
• /
• v.55 no.10
• /
• pp.3754-3767
• /
• 2023

Turbulent jet discharges in subcooled water pools are essential for safety systems in nuclear power plants, specifically in the pressure suppression pool of boiling water reactors and In-containment Refueling Water Storage Tank of advanced pressurized water reactors. The gas and liquid flow in these systems is investigated using multiphase flow analysis. This field has been extensively examined using a combination of experiments, theoretical models, and Computational Fluid Dynamics (CFD) simulations. ANSYS CFX offers two approaches to model multiphase flow behavior. The non-homogeneous Eulerian-Eulerian Model has been used in this work; it computes global information and is more convenient to study interpenetrated fluids. This study utilized the Large Eddy Simulation Model as the turbulence model, as it is better suited for non-stationary and buoyant flows. The CFD results of this study were validated with experimental data and theoretical results previously obtained. The figures of merit dimensionless penetration length and the dimensionless buoyancy length show good agreement with the experimental measurements. Correlations for these variables were obtained as a function of dimensionless numbers to give generality using only initial boundary conditions. CFD numerical model developed in this research has the capability to simulate the behavior of non-condensable gases discharged in water.

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

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

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