1 Mở đầu Hệ số pha (tỷ phần rỗng) đóng vai trò quan trọng mô hình dòng hai pha thang tỷ lệ mô phận nhà máy điện hạt nhân (NMĐHN) Trong chương trình CTF hệ số pha tính toán giải phương trình bảo toàn Sau đó, chế độ dòng chảy xác định dựa hệ số pha Chẳng hạn, chế độ dòng xác định dựa giá trị hệ số pha mô hình tường thông thường minh họa hình 2.8 luận án Bản đồ chế độ dòng mô hình tường gia nhiệt thông thường gồm: ■ Chế độ bong bóng nhỏ xác định hệ số pha nhỏ 0.2 ■ Chuyển tiếp bong bóng nhỏ-sang túi xác định hệ số pha khoảng (0.2, 0.5) ■ Chế độ khuấy/rối xác định hệ số pha khoảng (0.5 acrit) ■ Chế độ xuyến/sương mù xác định hệ số pha lớn acrit Khi đó, với chế độ dòng riêng biệt mô hình tường gia nhiệt thông tường, diện tích mặt phân cách, lực cản mặt phân cách hệ số truyền nhiệt qua mặt phân cách xác định khác Hình 2.8 Bản đồ chế độ dòng chương trình CTF với mô hình tường thông thường (nguồn [38]) 1.1 Mục tiêu luận án Nghiên cứu hệ số pha kênh tải nhiệt vấn đề quan trọng phân tích an toàn thủy nhiệt, tiến hành với nhiều thực nghiệm chương trình mô Do dựa tư liệu đề tài cấp nhà nước (mã số ĐTĐL.2011-G/82), nội dung nghiên cứucủa luận án đặt dự đoán hệ số pha kênh nóng vùng hoạt lò phản ứng VVER-1000/V392 sử dụng phương pháp tiếp cận ước lượng tốt Trong trường hợp chương trình liên kết hai thang tỷ lệ mô RELAP-3D MARS-3D sẵn, kỹ thuật phân tích nhiều thang tỷ lệ coi cách tiếp cận ước lượng tốt nghiên cứu toán thủy nhiệt Trong luận án này, chương trình tính toán nơtron MCNP5 chương trình thủy nhiệt với thang tỷ lệ khác gồm: RELAP5, CTF CFX (Ansys CFX) sử dụng để dự đoán hệ số pha kênh nóng lò VVER-1000/V392 Luận án đặt hai nội dung nghiên cứu cụ thể sau: ■ Xây dựng quy trình phân tích nhiều thang tỷ lệ để dự đoán hệ số pha trình chuyển tiếp sở sử dụng chương trình máy tính: MCNP5, RELAP5 CTF; ■ Xem xét việc kết hợp chương trình CTF CFX nhằm nâng cao khả dự đoán hệ số pha CTF Ở đây, việc kết hợp vận dụng chương trình CTF CFX nhằm cải thiện việc dự đoán hệ số pha CTF vấn đề nghiên cứu Như thường lệ, chương trình CTF dung dự đoán hệ số pha xảy chuyển tiếp CFX sử dụng để dự đoán trạng thái dừng Có thể hy vọng chương trình CFX với thang mô nhỏ đem lại cải thiện việc dự đoán hệ số pha trạng thái dừng thời điểm cụ thể 1.1.1 Đối tượng nghiên cứu Hệ số pha kênh nóng lò VVER-1000/V392 dự đoán với thang tỷ lệ khác điều kiện chuyển tiếp 40 giây kể từ bắt đầu cố LOCAs với kích thước vỡ khác 1.1.2 Phạm vi nghiên cứu 1.2 ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ Do tính phức tạp dòng hai pha, phạm vi nghiên cứu giới hạn dòng hai pha kênh thẳng đứng tương ứng với chế độ dòng cụ thể bong bóng nhỏ, túi, khuấy xuyến Bố cục luận án Luận án gồm bốn chương phần kết luận cuối Chương nêu luận giải vấn đề dẫn đến nội dung nghiên cứu luận án với vấn đề sau: Hiện trạng phát triển điện hạt nhân giới Việt Nam; Sơ lược an toàn hạt nhân; Phân tích an toàn thủy nhiệt điều kiện chuyển tiếp; ■ Tìm hiểu công nghệ lò phản ứng VVER Việt Nam liên quan đến nghiên cứu luận án; Mục tiêu luận án; Bố cục luận án Chương hai trình bày phương pháp luận liên quan đến phân tích nhiều thang tỷ lệ mô hình vật lý chương trình với thang tỷ lệ khác nhau, chủ yếu tập trung vào mô hình chuyển pha chương trình RELAP5, CTF CFX, cụ thể sau: Tiếp cận nhiều thang tỷ lệ mô lò phản ứng nước nhẹ (LWR); Chương trình hệ thống RELAP5; Chương trình phân tích kênh CTF; Chương trình thang tỷ lệ nhỏ (meso) CFX; Kết luận Chương ba trình bày việc đánh giá mô hình chuyển pha chương trình nêu thực sở mô so sánh với kết thực nghiệm Các kết mô mức hệ thống so sánh với báo cáo phân tích an toàn (Dự án Belene) Việc đánh giá chương trình CTF thực qua mô thực nghiệm BM ENTEK, vốn toán chuẩn (ISP) dùng để khảo sát mô hình dòng sôi bó nhiên liệu lò phản ứng RBMK Liên Bang Nga Thang tỷ lệ nhỏ (meso) kiểm chứng với thực nghiệm PSBT đơn kênh, toán chuẩn (ISP) Như thế, chương ba gồm phần sau: Khái quát lò phản ứng VVER-1000/V392; ■ Kiểm chứng mô hình mô chương trình RELAP5 lò VVER-1000/V392 sở đối chứng với báo cáo phân tích an toàn; ■ Kiểm tra đánh giá mô hình chuyển pha chương trình CTF thực nghiệm BM ENTEK; ■ Kiểm tra mô hình chuyển pha Ansys CFX với thực nghiệm PSBT cho toán đơn kênh; ■ Kết luận Những mục tiêu đặt luận án chương giải chương bốn theo trình tự sau: ■ Sơ đồ tính toán; ■ Tính toán phân bố công suất chương trình MCNP5; ■ Mô cố LOCAs chương trình RELAP5; ■ Mô chuyển tiếp cố LOCAs chương trình CTF ■ Mô dừng thời điểm cụ thể cố LOCAs chương trình CFX; ■ Kết luận Tổng quan mô hình chuyến pha chương trình có thang tỷ lệ khác Chương hai trình bày nội dung sau: 2.1 Tiếp cận nhiều chương trình nhiều thang tỷ lệ mô thủy nhiệt lò PWR; 2.1.1 Các chương trình tính toán nơtron tính toán thủy nhiệt; 2.1.2 Các thang tỷ lệ khác chương trình thủy nhiệt; 2.1.3 Các phương pháp tiếp cận mô thủy nhiệt khác ; 2.2 Mô hình chuyển pha chương trình hệ thống RELAP5; 2.3 Mô hình chuyển pha chương trình phân tích kênh CTF; 2.3.1 Hóa ngưng tụ chuyển pha nhiệt ; 2.3.2 Hóa ngưng tụ pha trộn rối trượt hệ số pha hơi; 2.4 Các mô hình chuyển pha chương trình thang tỷ lệ nhỏ (meso) CFX; 2.4.1 Hóa tường; 2.4.2 Ngưng tụ khối chất lỏng; ■ 2.5 Kết luận Một cách tóm tắt, chương nêu việc vận dụng nhiều chương trình nhiều thang tỷ lệ mô lò PWR, đặc biệt việc dự đoán hệ số pha Các mô hình chuyển pha chương trình vận dụng luận án gồm RELAP5, CTF CFX trình bày cách vắn tắt Kiếm tra đánh giá mô hình chuyến pha mô số Chương ba trình bày nghiên cứu luận án kiểm tra đánh giá mô hình chuyển pha mô số với thang tỷ lệ khác nhau, cụ thể chương trình RELAP5 mô toàn hệ thống lò VVER-1000/V392, chương trình CTF mô thực nghiệm ENTEK BM chương trình CFX mô thực nghiệm PSBT toán đơn kênh Trên cở sở đó, tiến hành việc đánh giá độ xác dự đoán hệ số pha chương trình thựcnghiệm PSBT Một cách tóm tắt, kết khoa học đạt chương sau: 3.1 Kiểm tra mô hình mô RELAP5 lò Bảng 3.2 So sánh trạng thái dừng mô lò VVERVVER1000/V392 sở so sánh với báo cáo phân tích Tham số an toàn Mục đích phần nhằm kiểm tra mô hình mô chương trình RELAP5 xây dựng luận án lò VVER- 1000/V392 Do đó, kết so sánh với báo cáo phân tích an toàn kịch cố LOCAs Việc so sánh thể kết trạng thái dừng, thời điểm xảy kiện trình chuyển tiếp biểu nhiệt độ đỉnh vỏ bọc nhiên liệu 3.2.2 Kiểm tra mô trạng thái dừng Các kết cho thấy, tính toán cho trạng thái dừng phù hợp với giá trị thiết kế giá trị báo cáo phân tích an toàn [35] thể sai số giá trị nằm khoảng 4% 3.2.3 Kiểm tra thông qua kịch cố ■ Các kết so sánh dãy kiện kiện số trình bày bảng 3.4 Bảng 3.4 Dãy kiện số theo SAR [35] Thời 03 04 907 00 17 30 40 0 60 0 63 0 117 00 500 00 (hỏng máy DG DG khác tình trạng bảo dưỡng) - Dừng toàn hệ bình thường; Xuất tín hiệu dập lò áp suất phía vùng hoạt 14.70 MPa Ngưỡng khởi động bơm hệ cấp nước vùng hoạt khẩn cấp (ECCS) Xuất hiệnđạttínđến hiệu dập lò Bắt đầu phun nước vào vùng hoạt từ hệ Mở van hệ tích nước HA-2 kết nối tới vùng hoạt áp suất xuống Kết nối1.50 tới MPa hệ tải nhiệt dư thụ động Bắt đầu phun nước từ bơm hệ cấp nước vùng hoạt khẩn cấp (ECCS) Các hệ tải nhiệt dư thụ động SG 1, PHRS hoạt động đủ công suấtphun từ bể Kết thúc tích nước HA-1vào Bắt cấp nước từ vùngđầu hoạt bể KếtHA-2 thúc tính toán 08 03 90 00 15 0 30 040 0 75 0 63 011 50 khác tình trạng bảo dưỡng) - Dừng toàn hệ thống vận ành bình thường; Xuất tín hiệu dập lò áp suất phía vùng hoạt 14.70 Ngưỡng khởi động bơm hệ cấp nước vùng hoạt khẩn cấp (ECCS) đạt đến Xuất tín hiệu dập lò Bắt đầu phun nước vào vùng hoạt từ hệ tích Mở van hệ tích nước HA-2 kết nối tới vùng hoạt áp suất xuống Kết nối1.50 tới MPa hệ tải nhiệt dư thụ động (SG Bắt đầu phun nước từ bơm hệ cấp nước vùng hoạt khẩn cấp (ECCS) Các hệ tải nhiệt dư thụ động SG 1, PHRS hoạt động đủ công suất Kết thúc phun từ bể tích nước HA- 1vào vùng Bắt hoạtđầu cấp nước từ bể HA-2 Kết thúc tính toán 0.0 0.0 1 0 15 Run No 1.22 31 1.22 11 1.43 21 1.43 23 1.52 21 1.22 21 1.24 21 1.32 21 1.22 33 1.21 21 1.62 21 1.22 12 1.22 34 Nhận thấy chương trình CTF có xu hướng dự đoán thấp thực nghiệm hệ số pha nhỏ 0.2 tương ứng đồ dòng chảy bong bóng nhỏ Tại vị trí gần lối kênh, mà giá trị thực nghiệm lớn 0.2 tương ứng với miền sôi bão hòa, chương trình CTF lại có xu dự đoán cao thực nghiệm Như vậy, mô hình sôi CTF cần cải thiện cho chế độ sôi bão hòa sôi bão hòa cho hệ số pha sinh phù hợp thực nghiệm 3.4 Kiểm tra mô hình CFX thí nghiệm PSBT đơn kênh 3.4.6 Độ xác kết mô CFX so sánh với thực nghiệm PSBT đơn kênh Kết cho thấy với hiệu chỉnh số Nusselt theo Ranz Marshall, CFX cho kết dự đoán pha ứng với chế độ dòng bong bóng nhỏ chế độ truyền nhiệt nhiệt độ ngưng tốt CTF Bảng 3.24 So sánh kết tính toán CFX CTF thực nghiệm đo hệ số pha Press Mass Pow Inlet Exp CFX CTF hea D(CF ure er Tem Void t X)' Flux (kg/c (10 kg/ (kW p Mo m) ) ( C) de m2h) 0.01 0.01 0.00 sub 150.2 10.87 60.0 299 3 c 0.001 150.1 10.91 90.0 295 0.03 0.07 0.07 sub 0.037 c 100.5 5.01 59.9 209 0.04 0.00 0.00 sub c 0.037 100.5 5.03 59.9 229 0.04 0.06 0.08 sub 0.015 c 75.5 5.02 49.9 219 0.04 0.06 0.08 sub 0.015 c 150.1 10.88 69.8 299 0.04 0.03 0.00 sub c 0.009 150.2 5.02 59.8 268 0.05 0.07 0.04 sub 0.019 c 125.0 11.10 59.9 294 0.05 0.08 0.05 sub 0.032 c 150.2 10.89 59.9 309 0.06 0.07 0.02 sub 0.018 c 150.1 14.80 79.9 309 0.06 0.07 0.02 sub 0.012 5 c 50.5 5.01 50.0 189 0.07 0.04 0.10 sub c 0.026 150.1 10.88 90.0 299 0.07 0.10 0.10 sub 0.024 9 c 150.1 10.92 0.08 0.12 0.07 sub 60.1 314 0.045 c o D( CT F) 0.0 0.0 36 0.0 0.0 38 0.0 39 0.0 0.0 0.0 01 0.0 0.0 0.0 28 0.0 30 0.0 1.1221 Kết cho thấy với hiệu số Nusselt theo Warierr với 169.1 10.95 áp suất nhỏ 122 bar chế dòng chuyển tiếp sang bong 49.9 bóng lớn tương ứng với chế độ truyền nhiệt sôi bão hòa CFX 329.7 0.087 thường cho dự đoán CTF cho dự đoán thực 0.109 3.25 0.032 3.25 So sánh kết CFX CTF với thực nghiệm Table subc miền sôi bão hòa 0.022 -0.055Press Mass Run Po Inle Exp CF CT hea D(CF D(C No ure wer t X F t X)* TF) Flux (kg/c (10 kg/ (k Voi Mo Tem ) W) d0.4 0.4 0.4 de m h) p 1.52 m 75.0 5.00 50 243 nuc 0.018 0.04 22 11 29 52 b 1.44 100.2 5.03 79 248 0.5 0.6 0.5 nuc 0.101 0.08 12 04 05 84 b 1.43 100.1 5.02 60 268 0.5 0.4 0.5 nuc 0.02 26 31 83 55 b 0.048 1.43 100.2 5.03 79 248 0.5 0.4 0.5 nuc 0.02 12 66 57 88 b 0.109 1.41 99.8 10.90 69 304 0.6 0.3 0.5 nuc 22 36 88 92 b 0.248 0.04 1.52 75.6 5.03 49 263 0.6 0.5 0.6 nuc 23 47 68 37 b 0.079 0.01 1.63 50.6 1.96 20 238 0.6 0.7 0.7 nuc 0.083 0.03 12 80 63 16 b 1.43 100.1 4.96 59 289 0.6 0.5 0.6 nuc 0.00 27 88 91 89 b 0.097 1.62 50.4 5.00 79 199 0.7 0.5 0.7 nuc 0.00 12 11 29 17 b 0.182 1.62 50.5 5.03 49 239 0.7 0.6 0.7 nuc 0.00 23 18 09 b 0.109 3.4.8 Nâng cao dự đoán lệ ia CFX miền sôi bão số p hòa Bảng 3.28 hệ số pha độ nhiệt (chênh lệch nhiệt tường nhiệt bão hòa) trước sau hiệu chỉnh thể cột tương ứng “Void”, “Tsup” and “Void*”, “Tsup*”được tính toán bới CFX với hiệu chỉnh số Nusselt theo Kim Park Run Voi Tsu Exp Scalin M Delt Void Tsup D(CF d p Void g A * X) aT * * F (K) 1.43 0.42 43.2 0.531 1.4 0.6 0.54 28.5 -0.103 (K) 1.2 26 8 1.43 0.40 62.5 0.566 1.4 0.7 2.4 0.55 48.6 -0.163 12 1.41 0.41 37.2 0.636 1.4 0.9 2.9 0.58 21.5 -0.226 22 1.52 0.52 39.2 0.647 1.4 3.0 0.63 19.9 -0.121 0.6 23 1.43 0.56 43.4 0.688 1.4 0.72 17.7 -0.121 0.8 1.0 27 1.62 0.71 35.2 0.718 1.4 0.8 3.2 0.73 27.9 -0.004 23 7 *D (CFX) = (Void - Exp Void), D (CFX) = (Void* - Exp Void) Kết luận Việc kiểm tra đánh giá mô hình chuyển pha D(CF X) *0.013 -0.007 -0.052 -0.009 0.033 0.019 chương đưa kết luận sau: Chương trình hệ thống RELAP5 với khả mô toàn hệ thống gắn với việc tải nhiệt từ vùng hoạt tới môi trường tản nhiệt cuối tiến hành với lò VVER-1000/V392 hệ thống an toàn liên quan Với mục tiêu kiểm tra mô hình mô luận án chương trình RELAP5, kịch cố chất tải nhiệt vỡ lớn (LBLOCA) tài liệu đáng tin cậy báo cáo phân tích an toàn dự án điện hạt nhân Belene (Bulgaria) sử dụng để đối chứng với kết mô Dựa yếu tố: (a) độ lệch thời điểm xảy kiện chuỗi kiện trình bày Bảng 10 tối đa 15 giây; (b) Biểu nhiệt độ đỉnh vỏ bọc nhiên liệu (PCT) khoảng 300 giây nhỏ 1200 oC thời gian làm mát hạ nhiệt vỏ bọc khoảng 280 giây minh họa Hình 3.4, kết cho thấy mô luận án phù hợp tốt với kết trình bày báo cáo phân tích an toàn (SAR) Đối với dự đoán hệ số pha vùng hoạt ứng với đường kính tương đương lớn kênh chương trình hệ thống RELAP5 chưa thể công cụ phù hợp Điều xuất phát từ hai lý do: (a) RELAP5 chương trình mô chiều (1D), kênh trung bình kênh nóng vùng hoạt mô hình hóa dạng ống có thay đổi hình học dòng chảy so với thực tế dẫn đến đồ chế độ dòng thay đổi so với thực tế; (b) mô hình - - - - - chuyển pha gần tường trình bày công thức (2.6) (2.7) phụ thuộc vào nhiệt độ gần tường Các mô hình sôi CTF có xu dự đoán thấp miền sôi nhiệt độ ngưng hệ số pha nhỏ 0.2 có xu dự đoán cao miền sôi bão hòa hệ số pha lớn 0.2 Với thí nghiệm ENTEK BM, kết dự đoán phân bố hệ số pha đưa CTF so với thực nghiệm xác hầu hết trường hợp khoảng sai số phép đo (0.03) độ lệch lớn CTF thực nghiệm 0.1 lối kênh dẫn xảy vài vị trí CFX đưa dự đoán pha tốt miền sôi nhiệt độ ngưng hệ số pha nhỏ 0.2 với sai số 0.03 tương ứng với độ lệch (1ơ) Vì dải hệ số pha nhỏ 0.2 miền sôi nhiệt độ ngưng, nên sử dụng kết dự đoán hệ số pha CFX thay cho CTF Đối với miền sôi bão hòa ứng với chế dòng chyển tiếp sang bong bóng lớn, CTF thường cho dự đoán CFX thường cho dự đoán thực nghiệm Do đường dự đoán CFX CTF biên biên dọc theo kênh dự đoán hệ số pha Việc cải tiến dự đoán hệ số pha cỉa CFX miền sôi bão hòa tỷ lệ hóa đường kính bong bóng hệ số tỷ phần cực đại cho tương xem tiếp cận để cải tiến mô hình sôi RPI miền sôi bão hòa Dự đoán hệ số pha kênh nóng bó nhiên liệu lò phản ứng VVER-1000/V392 4.2 Tính toán phân bố công suất vùng hoạt chương trình MCNP5 Các kết tính toán dựa mô toàn hình học vùng hoạt với số nơtron lịch sử 2.107 (sai số tương đối keff 14x10-5) Tỷ lệ phân bố công suất bó 1/6 vùng hoạt đưa hình 4.5 Bó nhiên liệu nóng có mã 30A9P với tỷ lệ công suất tương ứng 1.72 Tỷ số phù hợp với giả định phân bố công suất nằm dải (1.6 , 1.8) đề cập tài liệu [34] Hình 4.5 Tỷ lệ phân bố công suất 1/6 vùng hoạt Hình 4.6 cho thấy phân bố công suất theo chiều dọc vùng hoạt với hệ số đỉnh 1.52 phù hợp với phân bố công suất chu trình nhiên liệu [34], Để thiết lập mô hình CTF, cần có phân bố công suất nhiên liệu bó kết tính toán đưa từ chương trình MCNP5 Hình 4.7 đưa phân bố công suất cho hệ số công suất nóng 1.374 Hình 4.7 Phân bố công suất nhiên liệu bó nóng 4.4.1 Dự đoán hệ số pha kênh nóng bó nhiên liệu lò phản ứng VVER-1000/V392 RELAP5 CTF Bảng 4.2 Các trường hợp nghiên cứu dự đoán hệ số pha Trường Mô tả Diện Đường kính Vị trí vêt vỡ hợp tích vết tương vỡ (m ) đương (m) LB010 LOCA lớn cộng với 0.11 0.374 Chân lạnh 01 SBO-1 lớn cộng với 0.095 LB010 LOCA 0.348 Chân lạnh 02 SB0100 SBO-2 LOCA nhỏ cộng với 0.07 0.298 Chân lạnh 3SB0100 SBO-1 LOCA nhỏ cộng với 0.05 0.252 Chân lạnh SBO-2 (LB01001) (LB01002) (SB01003) (SB01004) Hình 4.13 Dự đoán hệ số pha CTF RELAP5 cho LOCA+SBO ■ Độ xác cho mô hình dự đoán pha - Mô hình hóa hình học dòng chảy Mô hình dòng chảy RELAP5 1D, dòng chảy qua kênh nóng 312 nhiên liệu mô giống dòng chảy ống với đường kính tương đương 0.02538 m2 trình bày bảng 2.8 Sự biến đổi hình học dòng chảy RELAP5 nguyên nhân dẫn đến sai khác chế độ dòng so với dòng chảy thực Điều ảnh hướng tới độ xác dự đoán hệ số pha Thêm mô hình CTF mô hình dòng chảy 3D minh họa hình 4.2 với P=12.75 mm D=9.1 mm Do hình học thang tỷ lệ CTF mô xác - Xác định nhiệt độ chất lỏng gần tường Với mô hình RELAP5, nhiệt độ lấy nhiệt độ trung bình bên thể tích kiểm soát Vì vậy, nhiệt độ chất lỏng gần tường giống nhiệt độ trung tâm dòng chảy Ví dụ, với diện tích dòng chảy 0.02538 m2 khoảng cách từ tâm dòng chảy đến tường gia nhiệt cm Với mô hình CTF với dòng chảy hình học kênh hình 4.2, khoảng cách từ tâm dòng chảy đến tường 0.275 cm Có thể thấy mô hình vật lý cho biến đổi pha RELAP5 CTF trình bày công thức 2.6), (2.7) (2.8), (2.17), phụ thuộc vào enthalpy (nhiệt độ) chất lỏng gần tường, CTF đưa enthalpy chất lỏng gần tường tốt RELAP5 nên kết tính toán chuyến pha xác Bên cạnh đó, biến đối hình học RELAP5 nên chế độ dòng chảy CTF đáng tin cậy Dựa lập luận kết dự đoán hệ số pha CTF với áp suất MPa đến 7MPa hệ thực nghiệm ENTEK BM với độ xác chấp nhận [33], từ kết luận CTF đưa kết dự đoán hệ số pha tốt RELAP5 4.5 Dự đoán hệ số pha cho đơn kênh chương trình CFX Xác định trường hợp SB01003-09-37 dựa việc lấy trường hơp (SB01003) từ bảng 4.2 thời điếm 09 s từ chuyến tiếp kênh 37 hình 4.11 Hình 4.11 mặt cắt ngang 1/12 bó nhiên liệu mô hình hóa tròn CTF Dựa kết luận chương 3, kết dự đoán hệ số pha CFX thay cho CTF trường hợp hệ số pha nhỏ 0.2 Đối với miền sôi bão hòa ứng với chế độ dòng chuyến tiếp sang bong bóng lớn coi kết CTF CFX đường biên dự đoán Bảng 4.7Dự đoán hệ số pha CTF CFX lối _kênh dẫn (z = 3.48m) _ (kênh Trường hợp CTF (bó) CTF CFX SB0l003-09-37 0.429 đơn) 0.444 0.4064 SB01003-14-34 0.173 0.152 0.1964 SB01003-16-15 0.094 0.1319 0.062 SB01003-16-30 0.146 0.1371 0.101 SB01003-20-15 LB01002-15-30 LB01002-20-18 LB01002-20-20 LB01002-30-30 0.153 0.395 0.424 0.442 0.609 0.200 0.361 0.356 0.438 0.640 0.2195 0.3200 0.2979 0.3954 0.6433 Hình 4.18 Dự đoán hệ số pha dọc theo kênh CTF CFX 4.6 Dự đoán hệ số pha cho bó CFX Theo kết luận mục 3.5 chương 3, kết CTF chế độ dòng bong bóng nhỏ ứng với truyền nhiệt nhiệt độ ngưng nên thay kết CFX miền sôi bão hòa kết dự đoán CTF CFX xem đường dự đoán Để minh họa điều cần lựa chọn miền tính toán nhiệt độ ngưng bão hòa bảng 4.8 4.9 Bảng 4.8 Miền sôi nhiệt độ ngưng cho khảo sát CFX Sub cooled selected region Inl Outl Inlet Inlet Ma Po et et ss wer void Pressure (m (m) flo (k Inlet temp Table 4.9 Mien bao hda cho khao sat bang CFX CTF va nen dnoc thay bang dn doan cua CFX mien soi dudi ngnng doi vdi cac trnong hop LB01001-06, LB01002-07, SB01003- 14 va SB01004-30 Cac thi ben phai cua hinh 4.21 the hien cac dnong dn doan he so pha hoi tren va dndi cho bdi CTF va CFX Viec CTF dn doan cao hon CFX trung vdi nhan dinh d muc 3.4 va 3.5 chnong Axial channel length [m] B-0 SB01004-30-CTF-NB A-A SB01004-30-CFX-NB Axial channel length (m) Hinh 4.21 Cái thien he so pha hoi bói CFX ó thi ben trái va dudng bien duái va bien tren cho du doán he so pha hoi cho bói thi ben phái 4.7 Kết luận Tóm tắt nghiên cứu thực chương này, thấy kết nghiên cứu dự đoán hệ số pha kênh nóng lò phản ứng VVER-1000/V392 nêu sau: - Thực công mô lò phản ứng VVER-V392 chương trình tính toán hệ thống RELAP5 với phân bố công suất tính chương trình tính toán MCNP5 cho chu trình nhiên liệu kiểm chứng đề cập chương ba - Chương trình CTF CFX sử dụng để dự đoán hệ số pha vùng hoạt dựa sở tham chiếu lẫn Với hệ số pha nhỏ 0.2 ứng với mô hình truyền nhiệt CTF sôi bão hòa, CTF có xu hướng dự đoán thấp trường hợp nên sử dụng kết CFX CFX đưa dự đoán tốt với sai số khoảng 0.03 - Do CTF có xu hướng dự đoán hệ số pha cao vùng sôi bão hòa, kết tính toán hệ số pha CTF CFX sử dụng làm biên biên kết thực tế Kết luận Như đề cập chương một, đối tượng luận văn khảo sát hệ số pha kênh nóng vùng hoạt lò phản ứng VVER- 1000/V392 Luận án nêu hai vấn đề cần giải quyết: ■ Xây dựng quy trình phân tích nhiều thang tỷ lệ để dự đoán hệ số pha trình chuyển tiếp sở sử dụng chương trình máy tính: MCNP5, RELAP5 CTF; ■ Xem xét việc kết hợp chương trình CTF CFX nhằm nâng cao khả dự đoán hệ số pha CTF Có thể thấy hai nội dung nêu giải luận văn Quy trình phân tích nhiều thang tỷ lệ hệ số pha giải chương bốn sở kiểm chứng mô hình sôi chương trình với thang tỷ lệ khác thực chương ba Việc đề xuất sử dụng chương trình tính toán CTF CFX cho phân tích thủy nhiệt đặc biệt dự đoán hệ số pha vùng hoạt vấn đề Có thể kết luận dự đoán hệ số pha vùng hoạt cải thiện trình tham chiếu lẫn hai chương trình CTF CFX CTF đưa dự đoán hệ số pha suốt trình chuyển tiếp kết thời điểm xác hóa nhờ trình tham chiếu với chương trình CFX Trong vùng sôi nhiệt độ ngưng với hệ số pha nhỏ (hệ số pha nhỏ 0.2 CTF) kết CFX xác sử dụng vùng sôi bão hòa (xác định CTF với hệ số pha khoảng 0.2 đến 0.5) kết dự đoán hệ số pha CTF CFX đường biên biên dự đoán hệ số pha Có thể tóm tắt kết luận văn sau: 5.1 Phát kết đạt luận án ■ Đề xuất phương pháp tiếp cận ước lượng tốt dự đoán hệ số pha cách sử dụng nhiều chương trình tính toán khác nhiều thang tỷ lệ bao gồm MCNP5, RELAP5, CTF cho phân tích diễn biến hệ số pha kênh nóng vùng hoạt trình chuyển tiếp; ■ Trong phân tích hệ thống sử dụng chương trình RELAP5 cho lò phản ứng VVER-1000/V392, nhiệt độ gần tường xác định không xác mô bó nhiên liệu kênh nóng, mô hình chuyển pha RELAP5 không đưa kết tin cậy dự đoán hệ số pha ■ Từ trình kiểm chứng kết CTF với thực nghiệm ENTEK BM cho thấy rằng: CTF có xu hướng dự đoán hệ số pha thấp vùng sôi nhiệt độ ngưng với hệ số pha nhỏ 0.2 CTF có xu hướng dự đoán cao vùng sôi bão hòa tương ứng với chế độ dòng chuyển tiếp sang bong bóng lớn ■ Từ trình kiểm chứng với thí nghiệm PSBT đơn kênh, mô hình thiết lập cho chương trình CFX luận án đạt hội tụ với sai số tương đối (RMS) 1e-6 ổn định hệ số pha với nghiên cứu độ nhạy mô hình vật lý Đối với vùng sôi nhiệt độ ngưng ứng với chế độ dòng bong bóng nhỏ (ag < 0.2), CFX đưa dự đoán hệ số pha với sai số khoảng ±0.03 ■ Trong vùng sôi bão hòa, mô hình sôi tường RPI CFX không phân chia xác thông lượng nhiệt thành phần: đối lưu, (quenching) hóa CFX đưa dự đoán hệ số pha thấp vùng sôi bão hòa Luận án đề xuất cách hiệu chỉnh đường kính bong bóng tách thành tỷ phần cực đại cho tượng để nâng cao độ xác CFX dự đoán hệ số pha miền bão hòa ■ Luận án quy trình áp dụng chương trình CTF CFX để dự đoán hệ số pha sau: (a) miền sôi nhiệt độ ngưng tương ứng với chế độ dòng bong bóng nhỏ sử dụng kết CFX; (b) miền sôi bão hòa đường dự đoán hệ số pha từ chương trình CTF CFX dọc theo kênh sử dụng làm biên biên để dự đoán hệ số pha vùng hoạt 5.2 Đề xuất nghiên cứu Áp dụng chương trình dạng CFD để khảo sát hệ số pha vùng hoạt đến thách thức Điều xuất phát từ tính phức tạp tượng liên quan đến sôi thiếu thực nghiệm điều kiện tương tự lò phản ứng nước áp lực (PWR) để kiểm chứng hiệu chỉnh mô hình CFD Dựa nghiên cứu luận án đề xuất số vấn đề cho nghiên cứu ■ Nghiên cứu cải tiến mô hình sôi tường RPI có chương trình CFX ( FLUENT) vùng sôi bão hòa Do mô hình sôi bão hòa nhiệt độ chất lỏng vị trí chí gần tường Do đó, tượng hóa (quenching) phải đóng vai trò chủ đạo ■ Tiến hành nhiều thí nghiệm với điều kiện tương tự lò phản ứng nước áp lực PWR để có phân bố hệ số pha sở để để kiểm chứng mô hình hóa ngưng tụ phần mềm CFX ■ Nghiên cứu dự đoán hệ số pha xác chương trình CFX sở trọng vào mô hình ngưng tụ chẳng hạn mô hình hiệu chỉnh số Nusselt điều kiện sôi khác [...]... nghiên cứu được thực hiện trong chương này, có thể thấy các kết quả chính nghiên cứu dự đoán hệ số pha hơi trong kênh nóng của lò phản ứng VVER-1000/V392 được nêu ra như sau: - Thực hiện thanh công mô phỏng lò phản ứng VVER-V392 bởi chương trình tính toán hệ thống RELAP5 với phân bố công suất tính bởi chương trình tính toán MCNP5 cho chu trình nhiên liệu đầu tiên và được kiểm chứng như đề cập trong. .. đoán hệ số pha hơi trong vùng hoạt là vấn đề rất mới Có thể kết luận rằng dự đoán hệ số pha hơi trong vùng hoạt được cải thiện bởi quá trình tham chiếu lẫn nhau của hai chương trình CTF và CFX CTF đưa ra dự đoán hệ số pha hơi trong suốt quá trình chuyển tiếp và kết quả tại một thời điểm được chính xác hóa nhờ quá trình tham chiếu với chương trình CFX Trong vùng sôi dưới nhiệt độ ngưng với hệ số pha hơi... trong vùng hoạt đến nay vẫn là một thách thức Điều này xuất phát từ tính phức tạp của các hiện tượng liên quan đến sôi và thiếu các thực nghiệm trong điều kiện tương tự lò phản ứng nước áp lực (PWR) để kiểm chứng và hiệu chỉnh các mô hình của CFD Dựa trên những nghiên cứu của luận án này có thể đề xuất một số vấn đề cho các nghiên cứu tiếp theo ■ Nghiên cứu cải tiến mô hình sôi tường RPI có trong chương... kênh nóng, khi đó mô hình chuyển pha trong RELAP5 không đưa ra kết quả tin cậy dự đoán hệ số pha hơi ■ Từ quá trình kiểm chứng kết quả của CTF với thực nghiệm ENTEK BM cho thấy rằng: CTF có xu hướng dự đoán hệ số pha hơi thấp hơn trong vùng sôi dưới nhiệt độ ngưng với hệ số pha hơi nhỏ hơn 0.2 và CTF có xu hướng dự đoán cao hơn trong vùng sôi bão hòa tương ứng với chế độ dòng chuyển tiếp sang bong bóng... số pha hơi như sau: (a) tại miền sôi dưới nhiệt độ ngưng tương ứng với chế độ dòng bong bóng nhỏ sẽ sử dụng kết quả của CFX; (b) tại miền sôi bão hòa đường dự đoán hệ số pha hơi từ chương trình CTF và CFX dọc theo kênh được sử dụng làm biên trên và biên dưới để dự đoán hệ số pha hơi trong vùng hoạt 5.2 Đề xuất nghiên cứu tiếp theo Áp dụng các chương trình dạng CFD để khảo sát hệ số pha hơi trong vùng. .. lệ phân bố công suất của mỗi bó trong 1/6 vùng hoạt được đưa ra trong hình 4.5 Bó nhiên liệu nóng nhất có mã 30A9P với tỷ lệ công suất tương ứng là 1.72 Tỷ số này là phù hợp với giả định phân bố công suất nằm trong dải (1.6 , 1.8) được đề cập trong tài liệu [34] Hình 4.5 Tỷ lệ phân bố công suất trong 1/6 vùng hoạt Hình 4.6 cho thấy phân bố công suất theo chiều dọc vùng hoạt với hệ số đỉnh là 1.52 cũng... luận án ■ Đề xuất phương pháp tiếp cận ước lượng tốt nhất trong dự đoán hệ số pha hơi bằng cách sử dụng nhiều chương trình tính toán khác nhau và nhiều thang tỷ lệ bao gồm MCNP5, RELAP5, CTF cho phân tích diễn biến hệ số pha hơi tại kênh nóng vùng hoạt trong quá trình chuyển tiếp; ■ Trong phân tích hệ thống sử dụng chương trình RELAP5 cho lò phản ứng VVER-1000/V392, đã chỉ ra được rằng nhiệt độ gần tường... hệ số pha hơi trong vùng hoạt dựa trên cơ sở tham chiếu lẫn nhau Với hệ số pha hơi nhỏ hơn 0.2 ứng với mô hình truyền nhiệt trong CTF là sôi dưới bão hòa, CTF có xu hướng dự đoán thấp hơn và trong trường hợp này nên sử dụng kết quả của CFX vì CFX đưa ra dự đoán tốt hơn với sai số khoảng 0.03 - Do CTF có xu hướng dự đoán hệ số pha hơi cao hơn ở vùng sôi bão hòa, vì vậy kết quả tính toán hệ số pha hơi... biên trên và biên dưới của kết quả thực tế 5 Kết luận Như đã đề cập ở chương một, đối tượng chính của luận văn là khảo sát hệ số pha hơi trong kênh nóng vùng hoạt lò phản ứng VVER- 1000/V392 Luận án nêu hai vấn đề cần giải quyết: ■ Xây dựng quy trình phân tích nhiều thang tỷ lệ để dự đoán hệ số pha hơi trong quá trình chuyển tiếp trên cơ sở sử dụng các chương trình máy tính: MCNP5, RELAP5 và CTF; ■ Xem... đường kính bong bóng và hệ số tỷ phần cực đại cho hiện tương tôi có thể xem là tiếp cận mới để cải tiến mô hình sôi RPI đối với miền sôi bão hòa 4 Dự đoán hệ số pha hơi trong kênh nóng bó nhiên liệu lò phản ứng VVER-1000/V392 4.2 Tính toán phân bố công suất trong vùng hoạt bởi chương trình MCNP5 Các kết quả tính toán dựa trên mô phỏng toàn bộ hình học vùng hoạt với số nơtron lịch sử bằng 2.107 (sai số