27Bảng 3: Kích thước vỏ thanh nhiên liệu và ng dố ẫn thanh điều khi n, thanh trung tâmể .... 15Hình 4 : Mô hình hóa trong chương trình tính toán COBRA-EN.. Các kết quả tính toán kết hợp
BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI TẠ DUY LONG XÂY DỰNG MƠ HÌNH TÍNH TỐN KẾT HỢP VẬT LÝ – THỦY NHIỆT CHO BĨ NHIÊN LIỆU LÕ PHẢN ỨNG VVER SỬ DỤNG CHƯƠNG TRÌNH MCNP5 VÀ COBRA-EN Chuyên ngành : KỸ THUẬT HẠT NHÂN LUẬN VĂN THẠC SĨ KHOA HỌC CHUYÊN NGÀNH KỸ THUẬT HẠT NHÂN NGƯỜI HƯỚNG DẪN KHOA HỌC: TS Trần Hoài Nam Hà Nội – Năm 2016 1708330017536df18d912-0bd3-4b88-a375-d17d2c58372a 170833001753698656bde-3537-44b9-9463-992607a4871f 170833001753606c232f7-748e-4cc2-8178-3839679e9dd3 LỜI CAM ĐOAN COBRA- T Duy Long LỜI CẢM ƠN qua MỤC LỤC L LI C DANH MC BNG BIU DANH MC HÌNH V DANH MC CÁC T VIT TT M u 1.1 Lý lựa chọn đề tài……………………………………………………………… 1.2 Mục tiêu nghiên cứu………………………………………………………………….9 1.3 Nội dung nghiên cứu……………………………………………………………… 10 lý thuyt 12 2.1 Tổng quan lò phản ứng VVER-1000………………………………………….12 2.2 Tổng quan hiệu ứng phản hồi nhiệt độ………………………………….14 2.3 Ảnh hưởng mật độ chất tải nhiệt lên tính tốn vùng hoạt lị phản ứng 16 Xây dng mơ hình tính tốn kt hp Vt lý Thy nhit 21 3.1 Mơ tả chương trình tính tốn……………………………………………… 21 3.1.1 21 3.1.2 -EN22 3.2 Mơ hình bó nhiên liệu lị phản ứng VVER-1000……………………………… 25 3.3 Mơ hình hóa bó nhiên liệu lị phản ứng VVER-1000 với chương trình MCNP5 COBRA-EN…………………………………………………………………………… 28 3.3.1 Mơ hình hóa bó nhiên liu vng trình MCNP5.28 3.3.2 Mơ hình COBRA-EN 32 3.4 Chương trình tính tốn kết hợp………………………………………………… 36 Các kt qu tính tốn 41 4.1 Tính tốn benchmark với chương trình MCNP5……………………………….41 4.2 Các kết tính tốn kết hợp Vật lý – Thủy nhiệt…………………………… 43 4.2.1 Phân b công sut bó nhiên liu.43 4.2.2 Phân b nhi m cht ti nhit46 4.2.3 Phân b nhi nhiên liu v bc nhiên liu 50 Kt lun 52 TÀI LIU THAM KHO 54 PH LC 55 Input mô hình hóa bó nhiên liệu với chương trình MCNP5………………… 55 Input mơ hình hóa bó nhiên liệu với chương trình COBRA-EN…………… 64 Chương trình cobra2mc (Tạo input cho MCNP5 từ output COBRA- EN) ……………………………………………………………………………………… 67 DANH MỤC BẢNG BIỂU Bảng 1: Một số đặc trưng nhiên liệu kênh chất tải nhiệt 24 Bảng 2: Mật độ hạt nhân đồng vị thành phần nhiên liệu vật liệu cấu trúc 27 Bảng 3: Kích thước vỏ nhiên liệu ống dẫn điều khiển, trung tâm 28 Bảng 4: Kết hệ số nhân vô hạn trạng thái S2 (Hot Full Power) [1] 41 Bảng 5: Phân bố cơng suất theo chiều cao bó nhiên liệu 45 Bảng 6: Phân bố nhiệt độ chất tải nhiệt theo chiều cao bó nhiên liệu (Đơn vị: K) 47 Bảng 7: Hệ số nhân vơ hạn theo vịng lặp tính tốn 48 Bảng : Mức sai khác vịng lặp mơ hình tính tốn kết hợp 48 DANH MỤC HÌNH VẼ Hình 1: Cấu hình vùng hoạt lị phản ứng VVER-1000 (Bushehr) 12 Hình 2: Bó nhiên liệu lị phản ứng VVER-1000 13 Hình 3: Mối tương quan nhiệt độ mật độ nước áp suất 15.7 Mpa 15 Hình : Mơ hình hóa chương trình tính tốn COBRA-EN 24 Hình 5: Tiết diện mặt cắt ngang bó nhiên liệu lị VVER-1000 26 Hình 6: Cấu hình mặt cắt ngang mạng nhiên liệu ống dẫn 26 Hình 7: Mơ hình bó nhiên liệu lị phản ứng VVER-1000 với chương trình MCNP5 31 Hình 8: Mơ hình hóa 1/6 bó nhiên liệu lị phản ứng VVER-1000 với chương trình COBRA-EN 32 Hình 9: Lưu đồ thuật tốn tính tốn kết hợp vật lý – thủy nhiệt 37 Hình 10: Phân bố cơng suất bó nhiên liệu lị phản ứng VVER-1000 42 Hình 11: Phân bố cơng suất theo chiều bán kính bó nhiên liệu lị VVER-1000 44 Hình 12: Phân bố cơng suất theo chiều cao bó nhiên liệu 46 Hình 13: Phân bố nhiệt độ chất tải nhiệt theo chiều cao bó nhiên liệu 47 Hình 14: Hệ số nhân vơ hạn k-inf theo bước lặp 49 Hình 15: Mức sai khác vịng lặp mơ hình tính tốn kết hợp 49 Hình 16: Phân bố nhiệt độ trung bình nhiên liệu theo chiều cao 50 Hình 17: Phân bố nhiệt độ nhiên liệu theo chiều bán kính 51 DANH MỤC CÁC TỪ VIẾT TẮT Các từ viết tắt Ý nghĩa VVER MCNP5 -Carlo N-Particle) COBRA-EN, COBRA-TF PERL MVP-GMVP -EN, COBRA-TF Practical Extraction and Reporting Language) - pháp Monte-Carlo FSAR Báo cáo phân tích an tồn (final safety analysis report) SRAC PARCS v2.7 ENDF Th BOC NJOY nhân PWR Chương Mở đầu 1.1 Lý lựa chọn đề tài lý H - (Neutronics-Thermal hydraulics) - àn (safety margins) t vài [2] - Zwermann und W Bernnat [3]; A Abarca, R Miró, T Barrachina and G Verdú [4] -TF (CTF) cho tính toán vùng Tuy nhiên, tính to lý 1.2 Mục tiêu nghiên cứu thay