• Nuclear Engineering and Technology
• /
• v.51 no.7
• /
• pp.1842-1852
• /
• 2019

The AP1000 reactor coolant pump is a vertical shielded-mixed flow pump, is the most important coolant power supply and energy exchange equipment in nuclear reactor primary circuit system, whose steadystate and transient performance affect the safety of the whole nuclear island. Moreover, safety demonstration of reactor coolant pump is the most important step to judge whether it can be practiced, among which software simulation is the first step of theoretical verification. This paper mainly introduces the fluid-solid coupling simulation method applied to reactor coolant pump, studying the feasibility of simulation results based on workbench fluid-solid coupling technology. The study found that: for the unsteady calculations of the pure liquid media, the average head of the reactor coolant pump with bidirectional fluid-solid coupling decreases to a certain extent. And the coupling result is closer to the real experimental value. The large stress and deformation of rotor under different flow conditions are mainly distributed on impeller and idler, and the stress concentration mainly occurs at the junction of front cover plate and blade outlet. Among the factors that affect the dynamic stress change of rotor, the pressure load takes a dominant position.

• Nuclear Engineering and Technology
• /
• v.51 no.7
• /
• pp.1776-1783
• /
• 2019

To accurately analyze the accidents in nuclear reactors, a thermohydraulic-neutronic coupling calculation is required to solve fluid dynamics and nuclear reactor kinetics equations in fine cells simultaneously and evaluate the local effects of neutronic and thermohydraulic parameters on each other. In the present study, a 3D thermohydraulic-neutronic coupling model is developed, validated and then applied for Isfahan MNSR (Miniature Neutron Source reactor) safety analysis. The proposed model is developed using FLUENT software and user defined functions (UDF) are applied to simulate the neutronic behavior of MNSR. The validation of the proposed model is first evaluated using 1mk reactivity insertion experiment into Isfahan MNSR core. Then, the developed coupling code is applied for a design basis accident (DBA) scenario analysis with the insertion of maximum allowed cold core reactivity of 4 mk. The results show that the proposed model is able to predict the behavior of the reactor core under normal and accident conditions with a good accuracy.

• Nuclear Engineering and Technology
• /
• v.55 no.9
• /
• pp.3213-3228
• /
• 2023

As an important part of the digital reactor, the pin-by-pin wise fine coupling calculation is a research hotspot in the field of nuclear engineering in recent years. It provides more precise and realistic simulation results for reactor design, operation and safety evaluation. CORCA-K a nodal code is redeveloped as a robust pin-by-pin wise neutronics and thermal-hydraulic coupled calculation code for pressurized water reactor (PWR) core. The nodal green's function method (NGFM) is used to solve the three-dimensional space-time neutron dynamics equation, and the single-phase single channel model and one-dimensional heat conduction model are used to solve the fluid field and fuel temperature field. The mesh scale of reactor core simulation is raised from the nodal-wise to the pin-wise. It is verified by two benchmarks: NEACRP 3D PWR and PWR MOX/UO₂. The results show that: 1) the pin-by-pin wise coupling calculation system has good accuracy and can accurately simulate the key parameters in steady-state and transient coupling conditions, which is in good agreement with the reference results; 2) Compared with the nodal-wise coupling calculation, the pin-by-pin wise coupling calculation improves the fuel peak temperature, the range of power distribution is expanded, and the lower limit is reduced more.

• Nuclear Engineering and Technology
• /
• v.55 no.4
• /
• pp.1350-1364
• /
• 2023

2D/1D coupling method is an important neutron transport calculation method due to its high accuracy and relatively low computation cost. However, 2D/1D coupling method may diverge especially in small axial mesh size. To analyze the convergence behavior of 2D/1D coupling method, a Fourier analysis for k-eigenvalue neutron transport problems is implemented. The analysis results present the divergence problem of 2D/1D coupling method in small axial mesh size. Several common attempts are made to solve the divergence problem, which are to increase the number of inner iterations of the 2D or 1D calculation, and two times 1D calculations per outer iteration. However, these attempts only could improve the convergence rate but cannot deal with the divergence problem of 2D/1D coupling method thoroughly. Moreover, the choice of axial solvers, such as DGFEM SN and traditional SN, and its effect on the convergence behavior are also discussed. The results show that the choice of axial solver is a key point for the convergence of 2D/1D method. The DGFEM SN based 2D/1D method could converge within a wide range of optical thickness region, which is superior to that of traditional SN method.

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

• Nuclear Engineering and Technology
• /
• v.55 no.1
• /
• pp.310-323
• /
• 2023

The k-eigenvalue neutronics/thermal-hydraulics coupling calculation is a key issue for reactor design and analysis. Jacobian-free Newton-Krylov (JFNK) method, featured with super-linear convergence rate and high efficiency, has been attracting more and more attention to solve the multi-physics coupling problem. However, it may converge to the high-order eigenmode because of the multiple solutions nature of the k-eigenvalue form of multi-physics coupling issue. Based on our previous work, a modified JFNK with a line search method is proposed in this work, which can find the fundamental eigenmode together with thermal-hydraulics feedback in a wide range of initial values. In detail, the existing modified JFNK method is combined with the line search strategy, so that the intermediate iterative solution can avoid a sudden divergence and be adjusted into a convergence basin smoothly. Two simplified 2-D homogeneous reactor models, a PWR model, and an HTR model, are utilized to evaluate the performance of the newly proposed JFNK method. The results show that the performance of this proposed JFNK is more robust than the existing JFNK-based methods.

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

• Nuclear Engineering and Technology
• /
• v.55 no.6
• /
• pp.1974-1987
• /
• 2023

The load-following operation of small modular reactors (SMRs) requires accurate prediction of transient behaviors that can occur in the balance of plants (BOP) and the nuclear steam supply system (NSSS). However, 1-D thermal-hydraulics analysis codes developed for safety and performance analysis have conventionally excluded the BOP from the simulation by assuming ideal boundary conditions for the main steam and feed water (MS/FW) systems, i.e., an open loop. In this study, we introduced a lumped model of BOP fluid system and coupled it with NSSS without any ideal boundary conditions, i.e., in a closed loop. Various methods for coupling boundary conditions at MS/FW were tested to validate their combination in terms of minimizing numerical instability, which mainly arises from the coupled boundaries. The method exhibiting the best performance was selected and applied to a transient simulation of an integrated NSSS and BOP system of a SMART. For a transient event with core power change of 100-20-100%, the simulation exhibited numerical stability throughout the system without any significant perturbation of thermal-hydraulic parameters. Thus, the introduced boundary-condition coupling method and BOP fluid system model can expectedly be employed for the transient simulation and performance analysis of SMRs requiring daily load-following operations.

• Journal of Mechanical Science and Technology
• /
• v.18 no.12
• /
• pp.2125-2136
• /
• 2004

In this paper, a numerical application algorithm for applying the CFAM (Consistent Fluid Added Mass) matrix for a core seismic analysis is developed and applied to the 7-ducts core system to investigate the fluid effects on the dynamic characteristics and the seismic time history responses. To this end, three cases such as the in-air condition, the in-water condition without the fluid coupling terms, and the in-water condition with the fluid coupling terms are considered in this paper. From modal analysis, the core duct assemblies revealed strongly coupled out-of-phase vibration modes unlike the other cases with the fluid coupling terms considered. From the results of the seismic time history analysis, it was also verified that the fluid coupling terms in the CFAM matrix can significantly affect the impact responses and the seismic displacement responses of the ducts.

• Nuclear Engineering and Technology
• /
• v.53 no.5
• /
• pp.1556-1568
• /
• 2021

The plate-type fuel assembly adopted in nuclear research reactor suffers from complicated effect induced by non-uniform irradiation, which might affect stress conditions, mechanical behaviors and thermal-hydraulic performance of the fuel assembly. This paper is the Part II work of a two-part study devoted to analyzing the complex unique mechanical deformation and thermal-hydraulic characteristics for the typical plate-type fuel assembly under irradiation effect, which is on the basis of developed and verified numerical thermal-fluid-structure coupling methodology under irradiation in Part I of this work. The mechanical deformation, thermal-hydraulic performance and Mises stress have been analyzed for the typical plate-type fuel assembly consisting of support plates under non-uniform irradiation. It was interesting to observe that: the plate-type fuel assembly including the fuel plates and support plates tended to bend towards the location with maximum fission rate; the hot spots in the fuel foil appeared at the location with maximum thickness increment; the maximum Mises stress of fuel foil was located at the adjacent location with the maximum plate thickness increment et al.