190 likes | 318 Vues
The 17th AER Symposium highlights the DiFis 2.0 code, a breakthrough in 3D Finite Element Neutron Kinetic simulations essential for nuclear reactor safety. Developed by A.I. Zhukov and A.M. Abdullayev, this code overcomes the conservative nature of existing 1D kinetic models, addressing SNRCU's requirement for 3D analyses during critical reactivity insertion accidents, such as Control Rod Ejection. DiFis 2.0 ensures enhanced accuracy with around 4% deviation compared to the established ANC-H code, making it a vital tool for reactor physics research.
E N D
17th AER Symposium on VVER Reactor Physics and Reactor Safety DiFis 2.0 – 3D Finite Element Neutron Kinetic Code A.I. Zhukov and A.M. Abdullayev NSC Kharkov Institute of Physics and Technology September 24 - 28, 2007 Yalta, Ukraine
Objective • Existing 1D kinetic codes too conservative • SNRCU intend to require 3D analysis for accidents with reactivity insertion (Control Rod Eject Accident and others)
Background • Westinghouse PHOENIX-H code (XS) • Westinghouse Advanced Nodal Code (ANC-H) • Westinghouse 1D-3D neutron kinetic code TWINKLE (for PWR square FA lattice) • DiFis 1.0 – 3D FE steady-state code • FE methodology
Nuclear Model 2 group approximation
Thermal Hydraulic Model Average Fuel Rod (one per assembly) for fuel, gap, cladding for coolant
Feedback Models • PHOENIX-H provides all XS with feedback • Moderator density correction includes: • Diffusion coefficient correction • Moderator absorption correction • Boron concentration correction • Spectrum correction • Fuel temperature correction • Xe, Sm, Pm correction • Control Rod correction
6 5 4 3 2 1 Finite Element Technique - 1 Expansion of fluxes into series over linear functions Set of simplex-functions Set of linear functions
Finite Element Technique - 2 Weighted Residual Equations or
12 11 14 13 10 9 34 33 15 2 1 8 36 35 32 31 16 3 6 7 38 37 12 11 30 29 17 4 5 24 39 14 13 10 9 28 18 19 22 23 40 15 2 1 8 27 20 21 41 16 3 6 7 26 42 17 4 5 24 25 43 18 19 22 23 54 44 45 20 21 52 53 46 47 50 51 48 49 6 2 1 5 3 6 4 3 4 5 2 1 Finite Element Technique - 3
Finite Element Technique - 4 Fuel Rod Normally Fuel – 7 zones Gap – 2 zones Clad – 5 zones Axial – 24 zones
Code capabilities • Core Power vs. time • Spatial and time power distribution • Spatial and time temperature distribution in fuel, cladding and coolant
Example of Steady-State Calculations EOC Relative FA Power ANC DiFis-ANC Max Err = 4.0% Rms Err = 1.6%
Example of Transient Calculations - 1 Control Rod Eject Accident Timeline start finish all banks start to drop finish 0 s 0.1 s 0.4 s 3.0 s
Example of Transient Calculations - 2 Rod Cluster Worth 0.18%
Example of Transient Calculations - 3 Doppler only FA under ejected RCCA
Example of Transient Calculations - 4 All Banks drop FA under ejected RCCA
Example of Transient Calculations - 5 Doppler only FA under ejected RCCA
Example of Transient Calculations - 6 All Banks drop FA under ejected RCCA
Summary DiFis 2.0 compatible with well-known codes such as PHOENIX-H and ANC-H DiFis 2.0 provides accuracy ~4% (in compare with ANC) for steady-state calculations DiFis 2.0 justification and benchmarking to be continued