Tổng quan nghiên cứu

Ngành năng lượng hạt nhân toàn cầu bước vào giai đoạn phát triển thế hệ mới với hơn 440 lò phản ứng thương mại vận hành, cung cấp khoảng 10% sản lượng điện năng toàn thế giới. Đối với các thiết kế lò phản ứng nước áp lực tiên tiến thế hệ III+, việc tối ưu hóa hiệu suất và nâng cao biên độ an toàn là điều kiện tiên quyết. Tại Việt Nam, nghiên cứu chuyên sâu về công nghệ hạt nhân phục vụ mục đích hòa bình, đặc biệt là các phân tích kỹ thuật gắn liền với công nghệ VVER, giữ vai trò cốt lõi trong xây dựng năng lực đánh giá an toàn quốc gia.

Vấn đề nghiên cứu trọng tâm của luận văn là đánh giá toàn diện hành vi cơ - nhiệt, thủy nhiệt động và bức xạ của thanh nhiên liệu hạt nhân TVS-2006 trong điều kiện vận hành ổn định của lò phản ứng VVER AES-2006 (công suất nhiệt 3.200 MWt, công suất điện 1.198,8 MWe). Mục tiêu cụ thể bao gồm việc làm rõ các cơ chế biến đổi vi cấu trúc, phân bố nhiệt độ, sự tương tác cơ học - hóa học giữa viên gốm nhiên liệu và vỏ bọc hợp kim zirconi, từ đó kiểm chứng các tiêu chuẩn kỹ thuật an toàn nhiệt vật lý.

Phạm vi nghiên cứu tập trung vào vùng hoạt lò VVER-1200 với chu trình nhiên liệu 4 năm và chu kỳ thay đảo 12 tháng, mô phỏng 163 bó nhiên liệu lục lăng tại các mức độ cháy mở rộng từ 60 đến 70 MWd/kgU. Kết quả của luận văn mang ý nghĩa thực tiễn to lớn: cung cấp bộ thông số định lượng chuẩn xác, hỗ trợ nâng cao biên dự trữ an toàn vận hành thêm khoảng 15%, giảm thiểu nguy cơ rò rỉ phóng xạ và tối ưu hóa quản lý nhiên liệu dài hạn.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu ứng dụng hệ thống lý thuyết đa lĩnh vực bao gồm cơ học rạn nứt vật liệu, truyền nhiệt phi tuyến trong môi trường xốp và nhiệt thủy lực hai pha. Mô hình truyền nhiệt dẫn suất dừng xuyên tâm từ tâm viên gốm qua khe hở chứa khí heli đến vỏ bọc zirconi E110 (Zr-1%Nb) và chất tải nhiệt nước nhẹ được thiết lập dựa trên phương trình cân bằng năng lượng dừng kết hợp hiệu ứng thay đổi độ dẫn nhiệt theo độ cháy.

Các khái niệm chuyên ngành cốt lõi gồm có:

  1. Hiện tượng tái cấu trúc vi mô viên gốm UO2 dưới gradient nhiệt độ lớn từ 2.000 đến 3.000 độ C/cm, tạo nên các vùng cấu trúc hạt đẳng trục và hạt dạng cột.
  2. Tương tác cơ nhiệt giữa viên gốm và vỏ bọc (PCI/PCMI) do giãn nở nhiệt và trương nở vật liệu.
  3. Cơ chế oxy hóa và hydro hóa vỏ bọc zirconi dẫn đến hiện tượng giòn hóa và ăn mòn ứng suất (SCC).
  4. Thông lượng nhiệt tới hạn (CHF) và hiện tượng rời khỏi chế độ sôi nhân (DNB) trong điều kiện áp suất 16,2 MPa.

Phương pháp nghiên cứu

Luận văn kết hợp chặt chẽ giữa phương pháp hồi cứu tư liệu pháp quy chuẩn quốc tế (IAEA, US.NRC, tiêu chuẩn pháp quy Liên bang Nga) và phương pháp mô phỏng số trực quan. Công cụ phân tích chính là phần mềm mô phỏng nhiệt - cơ thanh nhiên liệu FRAPCON-3.5. Lý do lựa chọn mã tính toán này là khả năng giải đồng thời hệ phương trình vi phân truyền nhiệt và cơ học biến dạng đàn - dẻo một cách chính xác dưới điều kiện chiếu xạ neutron cao.

Mô hình tính toán thực hiện trên tập mẫu đại diện gồm toàn bộ 163 bó thanh nhiên liệu trong vùng hoạt, phân tích chi tiết tại các vị trí thanh nhiên liệu chịu tải nhiệt cao nhất với chiều cao cột nhiên liệu 3,73 m và độ làm giàu tối đa 4,95% U-235. Dữ liệu đầu vào sử dụng thông số vận hành thực tế tại nhiệt độ làm mát lối vào 298,2 độ C, lối ra 328,9 độ C và lưu lượng dòng chảy 83.900 m3/h. Timeline nghiên cứu được thực hiện xuyên suốt quá trình tính toán mô phỏng 4 chu kỳ vận hành tương đương thời gian cháy tích lũy lên tới 49.000 EFPD (Effective Full Power Days).

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình mô phỏng và phân tích bằng công cụ FRAPCON-3.5 đã đưa ra những phát hiện định lượng quan trọng:

  1. Phân bố nhiệt độ xuyên tâm: Ở 100% công suất danh định, nhiệt độ tâm cực đại của viên gốm nhiên liệu đạt mức xấp xỉ 1.600 độ C, trong khi nhiệt độ bề mặt viên gốm là khoảng 470 độ C. Chênh lệch nhiệt độ qua bán kính viên gốm 3,8 mm lên tới hơn 1.130 độ C. Chênh lệch nhiệt độ qua khe hở khí đạt 100 độ C và qua bề dày vỏ bọc E110 là 23 độ C, giữ nhiệt độ mặt ngoài vỏ ở mức an toàn khoảng 350 độ C.
  2. Động học biến dạng hình học: Do hiện tượng ép nóng và tái cấu trúc mạng tinh thể trong 200 giờ đầu vận hành, mật độ viên gốm tăng lên đạt 98% đến 99% mật độ lý thuyết, gây co ngót cột nhiên liệu từ 3,5 cm đến 7,5 cm trước khi quá trình phồng nở do tích tụ sản phẩm phân hạch chiếm ưu thế ở độ cháy trên 40.000 MWd/tU.
  3. Ứng suất tiếp tuyến và rão vật liệu vỏ: Sự tiếp xúc cơ học giữa viên gốm và vỏ bọc xuất hiện khi độ cháy tăng cao, làm tăng ứng suất tiếp tuyến nhưng biến dạng dẻo tích lũy của hợp kim E110 vẫn duy trì dưới ngưỡng giới hạn 1,0%, thấp hơn 35% so với giới hạn an toàn quy định.
  4. Oxy hóa và tích tụ hydro: Chiều dày lớp oxit ZrO2 bề mặt ngoài sau 4 chu kỳ vận hành duy trì ở mức dưới 40 micromet, lượng hydro tích lũy trong vỏ bọc duy trì dưới ngưỡng 250 ppm, đảm bảo độ dẻo dai và ngăn ngừa nguy cơ gãy giòn.

Thảo luận kết quả

Các kết quả phân tích chỉ ra rằng thiết kế nhiên liệu TVS-2006 của lò VVER-AES2006 đáp ứng hoàn toàn các tiêu chuẩn kỹ thuật vận hành an toàn nghiêm ngặt nhất. Sự bổ sung chất hấp thụ neutron cháy được Gadolini (Gd2O3) trong hỗn hợp UO2-Gd2O3 đã kiểm soát hiệu quả độ phản ứng dự trữ ban đầu và làm phẳng mật độ công suất vùng hoạt, giảm hiện tượng đỉnh công suất cục bộ khoảng 12% so với các cấu hình truyền thống.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua các đường cong phân bố công suất dọc trục, đồ thị biến thiên nhiệt độ tâm viên gốm theo 4 chu kỳ vận hành và bảng tương quan giữa độ cháy tích lũy với áp suất khí phân hạch (xenon, krypton) thoát vào khoảng trống thanh nhiên liệu. So sánh với các nghiên cứu tương đương trên thanh nhiên liệu Zircaloy-4 của lò PWR phương Tây, hợp kim E110 (Zr-1%Nb) của Nga thể hiện tốc độ ăn mòn đồng đều thấp hơn 20% và độ bền rão vượt trội trong môi trường chiếu xạ nhiệt độ cao.

Đề xuất và khuyến nghị

Dựa trên kết quả nghiên cứu, 4 nhóm giải pháp chiến lược được đề xuất nhằm nâng cao hiệu quả quản trị công nghệ và an toàn hạt nhân:

  1. Thiết lập quy trình giám sát trực tuyến thông lượng nhiệt: Áp dụng hệ thống cảm biến đo đạc tại chỗ nhằm kiểm soát chặt chẽ biên độ dời khỏi sôi nhân (DNBR), mục tiêu duy trì biên an toàn nhiệt trên 1,30 trong suốt chu kỳ, thực hiện liên tục trong giai đoạn vận hành thương mại bởi đội ngũ kỹ sư vận hành nhà máy.
  2. Chuẩn hóa quy chuẩn kỹ thuật tiếp nhận nhiên liệu: Ban hành bộ tiêu chuẩn đánh giá khuyết tật mạng tinh thể và độ dày lớp phủ ngoài của thanh TVS-2006, giảm tỷ lệ hư hại thanh nhiên liệu xuống dưới 0,001%, hoàn thành trước quý IV năm 2026 bởi cơ quan pháp quy và viện nghiên cứu chuyên ngành.
  3. Tối ưu hóa chu trình nạp đảo thanh nhiên liệu: Triển khai mô hình nạp nhiên liệu phân tán kết hợp các thanh UO2-Gd2O3 để kéo dài chu kỳ nạp đảo từ 12 tháng lên 18 tháng, nâng hiệu suất kinh tế lò phản ứng thêm 4,5%, thực hiện trong lộ trình 2026 - 2028 bởi đơn vị quản lý vận hành.
  4. Đầu tư nâng cấp năng lực tính toán mô phỏng chuyên sâu: Ứng dụng tích hợp mã nguồn FRAPCON cùng các module tính toán thủy nhiệt 3D tiên tiến, nâng độ chính xác dự báo hành vi cơ lý lên 95%, thực hiện định kỳ hàng năm do các nhóm nghiên cứu vật lý lò chủ trì.

Đối tượng nên tham khảo luận văn

  1. Kỹ sư vận hành và an toàn nhà máy điện hạt nhân: Cung cấp tài liệu tra cứu chuyên sâu về trường nhiệt độ, giới hạn ứng suất và cơ chế hư hỏng nhiên liệu trong các kịch bản vận hành thực tế.
  2. Cán bộ quản lý pháp quy hạt nhân: Hỗ trợ cơ sở khoa học định lượng để thẩm định báo cáo phân tích an toàn (SAR), xây dựng tiêu chí cấp phép cho các tổ tổ hợp nhiên liệu thế hệ mới.
  3. Nghiên cứu sinh và học viên cao học ngành Vật lý hạt nhân/Kỹ thuật năng lượng: Cung cấp phương pháp luận chuẩn mực về mô hình hóa cơ - nhiệt, kỹ năng ứng dụng mã tính toán chuyên ngành FRAPCON.
  4. Các chuyên gia tư vấn công nghệ năng lượng: Làm tài liệu tham chiếu tin cậy khi so sánh các giải pháp công nghệ lò phản ứng nước áp lực VVER so với PWR và BWR trên thị trường quốc tế.

Câu hỏi thường gặp

Điểm khác biệt cốt lõi giữa thanh nhiên liệu VVER-AES2006 và PWR phương Tây là gì?
Thanh nhiên liệu TVS-2006 được sắp xếp trong bó nhiên liệu hình lục lăng với mạng kênh tam giác, sử dụng vỏ bọc hợp kim Zirconi E110 (Zr-1%Nb). Thiết kế này mang lại khả năng phân bố dòng làm mát đồng đều hơn và nâng cao tính chịu ăn mòn so với cấu trúc bó vuông dùng Zircaloy-4 của lò PWR.

Tại sao nhiệt độ tâm viên gốm nhiên liệu lại đạt mức cao tới 1.600 độ C?
Mức nhiệt độ này bắt nguồn từ độ dẫn nhiệt thấp của gốm UO2 kết hợp với mật độ sinh nhiệt thể tích rất cao do phản ứng phân hạch hạt nhân. Khoảng cách truyền nhiệt 3,8 mm tạo ra gradient nhiệt độ lớn hơn 2.000 độ C/cm truyền nhiệt ra chất làm mát ở 329 độ C.

Mã tính toán FRAPCON-3.5 đóng vai trò gì trong đánh giá thanh nhiên liệu?
FRAPCON-3.5 là công cụ tính toán mô phỏng trạng thái dừng chuyên dụng, cho phép tính toán kết hợp truyền nhiệt, phát tán khí phân hạch, phồng nở viên gốm, oxy hóa và biến dạng cơ học của vỏ bọc trong suốt tuổi thọ làm việc của thanh nhiên liệu hạt nhân.

Hợp kim Zirconi E110 có ưu điểm gì trong môi trường lò phản ứng?
Hợp kim E110 với thành phần Zirconi bổ sung 1% Niobi có tiết diện hấp thụ neutron nhiệt cực thấp, khả năng chống ăn mòn đồng đều vượt trội và giảm thiểu tối đa sự hấp thụ hydro, ngăn chặn hiện tượng rão giòn dưới thông lượng bức xạ neutron cao.

Chất hấp thụ Gd2O3 có vai trò gì trong việc kiểm soát vùng hoạt?
Gadolini là đồng vị có tiết diện bắt giữ neutron rất lớn nhưng tự tiêu hao trong quá trình cháy. Việc pha trộn Gd2O3 vào UO2 giúp kìm hãm độ phản ứng dư ban đầu khi độ làm giàu đạt xấp xỉ 5%, đồng thời hỗ trợ san phẳng trường công suất trong lò phản ứng.

Kết luận

  • Luận văn đã giải quyết toàn diện bài toán phân tích các đặc trưng cơ - lý - nhiệt - hóa và bức xạ của thanh nhiên liệu TVS-2006 trong lò phản ứng hạt nhân thế hệ III+ VVER-AES2006.
  • Ứng dụng thành công mã số tính toán FRAPCON-3.5 để mô hình hóa chính xác sự biến thiên nhiệt độ xuyên tâm, ứng suất cơ học và mức độ giải phóng khí phân hạch qua 4 chu kỳ vận hành.
  • Khẳng định các chỉ số nhiệt - cơ của thanh nhiên liệu TVS-2006 hoàn toàn thỏa mãn các tiêu chuẩn chấp thuận an toàn nghiêm ngặt của IAEA và Liên bang Nga.
  • Cung cấp cơ sở dữ liệu khoa học tin cậy phục vụ công tác thẩm định an toàn, đào tạo nguồn nhân lực chất lượng cao trong ngành vật lý lò phản ứng.
  • Đề xuất lộ trình nghiên cứu mở rộng giai đoạn 2026 - 2030 tập trung vào việc mô phỏng thanh nhiên liệu trong các điều kiện chuyển tiếp và kịch bản sự cố mất chất tải nhiệt (LOCA).

Để tiếp cận toàn bộ dữ liệu tính toán chi tiết, biểu đồ phân bố và hệ phương trình mô hình hóa chuyên sâu, độc giả có thể liên hệ trực tiếp đơn vị lưu trữ học thuật để khai thác toàn văn công trình nghiên cứu.