Development and Verification of Multi-Physical Coupling Calculation Code for Typical PWR Core

Author(s):  
Zhigang Li ◽  
Ping An ◽  
Junjie Pan ◽  
Wei Liu ◽  
Wei Lu ◽  
...  

Abstract Strong feedback phenomenon between the neutronics physical and thermal-hydraulic has an important impact on the design and safety analysis of pressured water reactor (PWR). In order to accurately simulate the strong coupling effect of neutron physics and thermal hydraulic in PWR, a reactor multi-physical coupling calculation code (ARMcc) is developed, in which the three-dimensional space-time neutron dynamic equation is solved by nodal expansion method (NEM) and nodal Green's function method (NGFM), the coolant temperature and fuel temperature are solved by single channel model and the cylinder heat conduction model respectively. The 3D IAEA benchmark, the 3D Langenbuch Maurer Werner (LMW) benchmark and NEACRP-L-335 benchmark are used to verify the neutronics model and coupling calculation solution ability respectively. The results show that: 1) the NEM and NGFM have high accuracy in solving the three-dimensional space-time neutron dynamics equation; 2) the results of neutronics and thermal-hydraulic coupling steady/transient calculation such as core relative power and fuel Doppler temperature are in good agreement with those of the NEACRP-L-335 benchmark, and the calculation accuracy is equivalent to similar software such as PARCS. Four coupled neutron physics and thermal hydraulic calculation modes are used to analyze the influence of different neutron physics calculation methods and thermal hydraulic calculation methods on the key parameters of PWR transient process in this paper. The results show that the mode of NGFM + FVM can more accurately simulate the core relative power peak and Doppler temperature.

Author(s):  
Zhifeng Li ◽  
Hongchun Wu ◽  
Chenghui Wan ◽  
Tianliang Hu

In order to raise computation speed on the premise of enough numerical accuracy, the Predictor-Corrector Improved Quasi-Static (PC-IQS) method and Nodal Green’s Function Method (NGFM) were combined to solve the three-dimensional space-time neutron diffusion kinetics problems for Cartesian geometry. In addition, the improved quasi-static method and the Krylov algorithm were applied to solve the three-dimensional space-time neutron diffusion kinetics problems for cylindrical geometry. Based on the proposed model, the program of three-dimensional neutron space-time kinetic code has been tested by the two-dimensional and three-dimensional transient numerical benchmarks. The numerical results obtained by this work were in good agreement with the reference solutions.


2002 ◽  
Vol 471 ◽  
pp. 381-408 ◽  
Author(s):  
HARIS J. CATRAKIS ◽  
ROBERTO C. AGUIRRE ◽  
JESUS RUIZ-PLANCARTE ◽  
ROBERT D. THAYNE ◽  
BRENDA A. McDONALD ◽  
...  

Experiments have been conducted to investigate turbulent mixing and the dynamics of outer fluid interfaces, i.e. the interfaces between mixed fluid and pure ambient fluid. A novel six-foot-diameter octagonal-tank flow facility was developed to enable the optical imaging of fluid interfaces above the mixing transition, corresponding to fully developed turbulence. Approximately 10003 whole-field three-dimensional space– time measurements of the concentration field were recorded using laser-induced- fluorescence digital-imaging techniques in turbulent jets at a Reynolds number of Re ∼ 20 000, Schmidt number of Sc ∼ 2000, and downstream distance of ∼ 500 nozzle diameters. Multiple large-scale regions of spatially nearly uniform-concentration fluid are evident in instantaneous visualizations, in agreement with previous findings above the mixing transition. The ensemble-averaged probability density function of concentration is found to exhibit linear dependence over a wide range of concentration thresholds. This can be accounted for in terms of the dynamics of large-scale well- mixed regions. Visualization of the three-dimensional space–time concentration field indicates that molecular mixing of entrained pure ambient fluid is dynamically initiated and accomplished in the vicinity of the unsteady large scales. Examination of the outer interfaces shows that they are dynamically confined primarily near the instantaneous large-scale boundaries of the flow. This behaviour is quantified in terms of the probability density of the location of the outer interfaces relative to the flow centreline and the probability of pure ambient fluid as a function of distance from the centreline. The current measurements show that the dynamics of outer interfaces above the mixing transition is significantly different from the behaviour below the transition, where previous studies have shown that unmixed ambient fluid can extend across a wide range of transverse locations in the flow interior. The present observations of dynamical confinement of the outer interfaces to the unsteady large scales, and considerations of entrainment, suggest that the mechanism responsible for this behaviour must be the coupling of large-scale flow dynamics with the presence of small-scale structures internal to the large-scale structures, above the mixing transition. The dynamics and structure of the outer interfaces across the entire range of space–time scales are quantified in terms of a distribution of generalized level-crossing scales. The outer-interface behaviour determines the mixing efficiency of the flow, i.e. fraction of mixed fluid. The present findings indicate that the large-scale dynamics of the outer interfaces above the mixing transition provides the dominant contribution to the mixing efficiency. This suggests a new way to quantify the mixing efficiency of turbulent flows at high Reynolds numbers.


Sign in / Sign up

Export Citation Format

Share Document