1 / 19

DiFis 2.0 – 3D Finite Element Neutron Kinetic Code

17 th 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.

tejano
Download Presentation

DiFis 2.0 – 3D Finite Element Neutron Kinetic Code

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

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

  2. Objective • Existing 1D kinetic codes too conservative • SNRCU intend to require 3D analysis for accidents with reactivity insertion (Control Rod Eject Accident and others)

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

  4. Nuclear Model 2 group approximation

  5. Thermal Hydraulic Model Average Fuel Rod (one per assembly) for fuel, gap, cladding for coolant

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

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

  8. Finite Element Technique - 2 Weighted Residual Equations or

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

  10. Finite Element Technique - 4 Fuel Rod Normally Fuel – 7 zones Gap – 2 zones Clad – 5 zones Axial – 24 zones

  11. Code capabilities • Core Power vs. time • Spatial and time power distribution • Spatial and time temperature distribution in fuel, cladding and coolant

  12. Example of Steady-State Calculations EOC Relative FA Power ANC DiFis-ANC Max Err = 4.0% Rms Err = 1.6%

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

  14. Example of Transient Calculations - 2 Rod Cluster Worth 0.18%

  15. Example of Transient Calculations - 3 Doppler only FA under ejected RCCA

  16. Example of Transient Calculations - 4 All Banks drop FA under ejected RCCA

  17. Example of Transient Calculations - 5 Doppler only FA under ejected RCCA

  18. Example of Transient Calculations - 6 All Banks drop FA under ejected RCCA

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

More Related