Tổng quan về luận án

Trong bối cảnh quá trình chuyển dịch năng lượng toàn cầu đòi hỏi các lò phản ứng hạt nhân nước nhẹ thế hệ III/III+ phải nâng cao hiệu suất kinh tế thông qua việc tăng độ giàu nhiên liệu ($^{235}\text{U}$ lên đến 4,4 wt%–5,0 wt%) và kéo dài chu trình nhiên liệu lên 18–24 tháng (với độ cháy sâu đạt 50–70 GWd/t), việc kiểm soát độ phản ứng dự trữ dư ban đầu (Beginning of Cycle - BOC) mà không làm suy giảm các biên độ an toàn nhiệt - thủy lực là một thách thức cốt lõi của vật lý lò phản ứng. Luận án "Nghiên cứu mô phỏng và cải tiến thiết kế bó nhiên liệu lò phản ứng VVER-1000/V-320 sử dụng vi hạt $\text{Gd}_2\text{O}_3$ bằng chương trình MVP" do NCS Hoàng Thanh Phi Hùng thực hiện tại Viện Năng lượng nguyên tử Việt Nam (VAEC) dưới sự hướng dẫn của PGS. Trần Hoài Nam và TS. Hồ Mạnh Dũng đã giải quyết trực diện bài toán này.

Khoảng trống nghiên cứu (Research Gap) xuất phát từ nghịch lý nhiệt - nơtrôn trong thiết kế thanh nhiên liệu tích hợp chất hấp thụ cháy được ($\text{Integrated Burnable Absorber - IBA}$): việc trộn đều đồng nhất bột $\text{Gd}_2\text{O}_3$ vào gốm oxit urani $\text{UO}_2$ theo thiết kế truyền thống làm giảm nghiêm trọng độ dẫn nhiệt của viên nhiên liệu từ $5,8 - 2,7\text{ W/m}\cdot\text{K}$ xuống còn $3,8 - 2,6\text{ W/m}\cdot\text{K}$ trong dải nhiệt độ vận hành 300–1273 K (theo nghiên cứu của Amaya et al. [2] và Kim et al. [25]). Sự suy giảm truyền nhiệt này cản trở việc nâng công suất lò (power upgrading) và tiềm ẩn nguy cơ mất an toàn khi xảy ra sự cố mất chất tải nhiệt (LOCA). Mặc dù các nghiên cứu thực nghiệm cơ tính của Kim et al. [32] chứng minh việc nén $\text{Gd}_2\text{O}_3$ dưới dạng vi hạt ($\text{Burnable Poison Particles - BPP}$) giúp khôi phục độ dẫn nhiệt, việc đưa hạt vi mô phân bố ngẫu nhiên vào thanh nhiên liệu lại kích hoạt hiệu ứng tự che chắn không gian ($\text{spatial self-shielding effect}$), làm biến đổi phi tuyến tiết diện hấp thụ nơtrôn và động học suy giảm đồng vị, điều mà các chương trình tất định cổ điển (deterministic codes) không thể mô phỏng chính xác.

Luận án đặt ra 3 câu hỏi nghiên cứu và 3 giả thuyết khoa học tương ứng:

  • RQ1: Liệu việc cấu hình $\text{Gd}2\text{O}3$ dạng vi hạt phân bố ngẫu nhiên có bảo toàn được đặc trưng suy giảm độ phản ứng $k\infty$ tương đương thiết kế trộn đều trong giai đoạn 0–10 GWd/t? *(Giả thuyết H1: Tồn tại kích thước đường kính vi hạt tối ưu trong cấu hình 12 thanh IBA cho phép cân bằng hiệu ứng tự che chắn để đường cong $k\infty$ trùng khớp với mô hình tham chiếu)*.
  • RQ2: Làm thế nào để vừa duy trì kiểm soát độ phản ứng vừa tối ưu hóa phân bố công suất cục bộ nhằm hạ thấp hệ số đỉnh công suất ($\text{Power Peaking Factor - PPF}$)? (Giả thuyết H2: Tăng số thanh IBA từ 12 lên 18 kết hợp mở rộng đường kính vi hạt lên thang micromet lớn hơn sẽ tái phân bố mật độ thông lượng nơtrôn và giảm tải công suất thanh nóng nhất).
  • RQ3: Cấu hình vi hạt $\text{Gd}_2\text{O}_3$ có cho phép cắt giảm nồng độ axit boric ($\text{H}_3\text{BO}_3$) hòa tan trong chất tải nhiệt mà vẫn đảm bảo hệ số phản hồi nhiệt độ chất làm chậm ($\text{Moderator Temperature Coefficient - MTC}$) luôn âm? (Giả thuyết H3: Thiết kế IBA vi hạt tối ưu có khả năng thay thế một phần hoặc toàn bộ vai trò bù trừ độ phản ứng của Boron, làm cho MTC âm sâu hơn, gia tăng đặc tính an toàn nội tại).

Khung lý thuyết của nghiên cứu được thiết lập dựa trên: (1) Phương trình vận chuyển nơtrôn Boltzmann giải bằng phương pháp Monte Carlo năng lượng liên tục; (2) Mô hình hình học thống kê ($\text{Statistical Geometry Model - STGM}$) dựa trên hàm phân bố láng giềng gần nhất ($\text{Nearest Neighbor Distribution - NND}$); (3) Hệ phương trình biến đổi hạt nhân Bateman cải tiến trong module $\text{MVP_BURN}$.

Đóng góp đột phá của luận án bao gồm việc xác định kích thước vi hạt tối ưu 60 µm cho cấu hình 12 thanh IBA và đề xuất cấu hình cải tiến hoàn toàn mới gồm 18 thanh IBA chứa vi hạt đường kính 300 µm với tỉ lệ thể tích 3,33%, giúp giảm hệ số PPF trong suốt dải cháy 0–40 GWd/t, đồng thời duy trì hệ số MTC âm ngay cả trong điều kiện nồng độ Boron giảm 50% hoặc bằng 0 ppm. Phạm vi nghiên cứu bao quát ô mạng bó nhiên liệu lục giác tiêu chuẩn VVER-1000/V-320 (331 vị trí gồm 312 thanh nhiên liệu, 18 ống dẫn thanh điều khiển CPSAR và 1 ống đo trung tâm) với thư viện số liệu hạt nhân đánh giá $\text{JENDL-3.3}$.

Literature Review và Positioning

Nghiên cứu về chất hấp thụ nơtrôn cháy được trong lò phản ứng nước nhẹ ($\text{LWR}$) trải qua nhiều giai đoạn phát triển với những luồng quan điểm học thuật khác nhau:

[Các luồng nghiên cứu kiểm soát độ phản ứng dư]

Luồng nghiên cứu truyền thống tập trung vào việc tối ưu hóa hàm lượng và vị trí các thanh hấp thụ trộn đều bột oxide $\text{Gd}_2\text{O}_3$ hoặc $\text{Er}_2\text{O}_3$. Hirano et al. (2000) [23] và Lim et al. (2001) [34] đã ứng dụng thuật toán liệt kê đầy đủ và phương pháp ma trận phi tuyến để làm phẳng phân bố công suất trong bó nhiên liệu lò BWR. Francois et al. (2003) [16] cùng Haibach et al. (2004) [20] và Yilmaz et al. (2006) [72] tiếp tục phát triển các giải thuật Tabu Search và Genetic Algorithm (GA) nhằm xác định hàm lượng Gadolinia tối thiểu trong lò PWR. Tuy nhiên, các tác giả này hoàn toàn bỏ qua khía cạnh suy giảm độ dẫn nhiệt trong thanh nhiên liệu $\text{UO}_2\text{–Gd}_2\text{O}_3$.

Luồng nghiên cứu thứ hai đề xuất giải pháp vật liệu lai hoặc lớp phủ để khắc phục suy giảm độ dẫn nhiệt. Dandi et al. (2012) [11] kết hợp thanh $\text{WABA}$ ($\text{Al}_2\text{O}_3\text{/B}_4\text{C}$) và thanh chứa Erbium trong bó nhiên liệu PWR, giúp cải thiện kiểm soát độ phản ứng dư từ 11% đến 20% so với việc chỉ dùng Er. Choe et al. (2016) [10] đề xuất mô hình lai: thanh $\text{UO}_2\text{–}^{157}\text{Gd}_2\text{O}_3$ phủ lớp mỏng $\text{ZrEr}_2$. Mặc dù lớp phủ $\text{ZrEr}_2$ không làm giảm độ dẫn nhiệt của lõi $\text{UO}_2$, cấu trúc này lại chiếm dụng thể tích nhiên liệu, làm giảm lượng Uranium nạp ban đầu và gia tăng độ phức tạp trong chế tạo cơ khí.

Luồng nghiên cứu thứ ba tập trung vào vi hạt hấp thụ phân tán ($\text{Micro-particle Burnable Poison}$). Khi Kim et al. (2014) [25, 32] chứng minh bằng thực nghiệm rằng cấu trúc vi hạt $\text{Gd}_2\text{O}_3$ kích thước 25–53 µm giữ cho độ dẫn nhiệt của viên nhiên liệu ở mức cao ($5,8 - 2,7\text{ W/m}\cdot\text{K}$ so với $3,8 - 2,6\text{ W/m}\cdot\text{K}$ của hỗn hợp trộn đều trong dải 300–1273 K), vấn đề chuyển dịch sang bài toán tính toán nơtrôn học. Tranh biện khoa học nảy sinh giữa hai trường phái:

  • Trường phái tất định cổ điển: Cho rằng việc sử dụng vi hạt gây ra sự bất định lớn do hiệu ứng tự che chắn kép (Double Heterogeneity), khiến các mô hình khuếch tán đa nhóm nơtrôn (như $\text{CITATION}$ hay $\text{ANISN}$) bị sai lệch nghiêm trọng nếu không có hiệu chỉnh thông số Dancoff phức tạp.
  • Trường phái thống kê Monte Carlo: Cho rằng phương pháp ngẫu nhiên năng lượng liên tục kết hợp mô hình hình học thống kê ($\text{STGM}$) có thể giải quyết chính xác sự truyền nơtrôn qua môi trường hạt phân tán mà không cần đơn giản hóa hình học.

So sánh với các nghiên cứu quốc tế:

  1. Nghiên cứu của Reda et al. (2016) [52] trên lò PWR: Reda tăng số thanh chứa chất hấp thụ từ 12 lên 24 thanh trộn đều để kiểm soát $k_{eff}$ ở giai đoạn đầu chu trình, nhưng phải đánh đổi bằng việc giảm hệ số dẫn nhiệt cục bộ trên 24 thanh nhiên liệu. Luận án của NCS Hoàng Thanh Phi Hùng vượt trội hơn khi sử dụng vi hạt, vừa kiểm soát $k_\infty$, vừa bảo toàn độ dẫn nhiệt cao.
  2. Nghiên cứu của Galahom (2016) [17] và Choe et al. (2016) [10]: Các tác giả khảo sát tiết diện vi mô của $^{10}\text{B}$ ($\sigma_a = 3844\text{ b}$) và $^{167}\text{Er}$ ($\sigma_a = 650\text{ b}$), nhận thấy $^{167}\text{Er}$ có tốc độ cháy chậm, để lại lượng chất hấp thụ dư thừa lớn tại cuối chu trình (End of Cycle - EOC), làm giảm độ phản ứng của lò. Trong khi đó, luận án khai thác triệt để hai đồng vị $^{155}\text{Gd}$ ($\sigma_a = 60,737\text{ b}$) và $^{157}\text{Gd}$ ($\sigma_a = 252,912\text{ b}$) với cơ chế tự che chắn hạt vi mô, giúp chất hấp thụ cháy cạn hoàn toàn tại 10 GWd/t, loại bỏ triệt để hiện tượng phản ứng dư âm tại EOC.

Vị thế khoa học của luận án được xác lập bằng việc là công trình tiên phong thiết lập mô hình tính toán nơtrôn học hoàn chỉnh cho bó nhiên liệu hình học lục giác VVER-1000 sử dụng vi hạt $\text{Gd}_2\text{O}_3$ ngẫu nhiên, nối liền khoảng cách giữa nghiên cứu khoa học vật liệu thực nghiệm và vật lý vùng hoạt lò phản ứng.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng và làm sâu sắc thêm Lý thuyết Vận chuyển Nơtrôn trong Môi trường Dị thể Ngẫu nhiên ($\text{Stochastic Heterogeneous Transport Theory}$). Khi chất hấp thụ phân tán dưới dạng hạt có kích thước lớn hơn quãng đường tự do trung bình của nơtrôn nhiệt bên trong hạt ($\lambda_t \approx 8,3\text{ µm}$ đối với $\text{Gd}_2\text{O}_3$), hạt hoạt động như một "vật hấp thụ đen" ($\text{black absorber}$).

Tiết diện hấp thụ vĩ mô hiệu dụng chuyển dịch từ mô hình hàm mũ thuần nhất sang mô hình phụ thuộc hình học bề mặt hạt: $$\Sigma_a(t) = n_P \cdot \pi R_0^2 \left(1 - \frac{\int_0^t \Phi(t')dt'}{4 N R_0}\right)^2$$ Trong đó $n_P$ là mật độ số hạt, $R_0$ là bán kính vi hạt ban đầu, $N$ là mật độ nguyên tử, và $\Phi(t)$ là thông lượng nơtrôn.

Các mệnh đề lý thuyết được xác lập:

  • Mệnh đề 1 (Quy luật tự che chắn không gian): Tốc độ suy giảm độ phản ứng của $\text{Gd}_2\text{O}_3$ tỉ lệ nghịch với bán kính vi hạt $R_0$. Việc tăng kích thước hạt từ 40 µm lên 300 µm làm chậm tốc độ cháy của lớp vỏ ngoài vi hạt, cho phép kéo dài hiệu ứng kiểm soát độ phản ứng mà không cần tăng tỷ lệ khối lượng chất hấp thụ.
  • Mệnh đề 2 (Cơ chế suy thoái đồng vị kép): Đồng vị $^{157}\text{Gd}$ với tiết diện bắt nơtrôn khổng lồ ($252,912\text{ b}$) tạo ra lớp che chắn bảo vệ cho lõi hạt chứa $^{155}\text{Gd}$ ($60,737\text{ b}$); sự cháy kiệt tuần tự từ ngoài vào trong của hạt vi mô tái cấu trúc phổ nơtrôn nhiệt tại vùng lân cận thanh nhiên liệu.
  • Mệnh đề 3 (Bảo toàn phản xạ chất làm chậm): Việc giảm nồng độ Boron hòa tan kết hợp tăng số thanh IBA vi hạt duy trì tỷ số làm chậm nơtrôn tối ưu, triệt tiêu nguy cơ xuất hiện hệ số MTC dương ở đầu chu trình.

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp liên kết 3 trụ cột lý thuyết:

[Khung phân tích nơtrôn - cấu trúc vi mô]

Sự kết hợp này được hiện thực hóa qua hệ thống mã mô phỏng MVP/GMVP II và module $\text{MVP_BURN}$ do Viện Nghiên cứu Năng lượng Nguyên tử Nhật Bản (JAEA) phát triển [45].

Điều kiện biên và ràng buộc hình thái của khung phân tích:

  • Áp dụng điều kiện biên phản xạ tuyệt đối ($\text{specular reflective boundary conditions}$) trên 6 mặt của ô mạng bó nhiên liệu lục giác VVER-1000.
  • Bước cháy khảo sát chi tiết từ 0 đến 40 GWd/t (tập trung cao độ tại pha chuyển tiếp 0–10 GWd/t với các bước chia nhỏ 0,1; 0,5; 1,0 GWd/t).
  • Nhiệt độ nhiên liệu và chất làm chậm được thiết lập tại điều kiện vận hành danh định chuẩn của lò VVER-1000/V-320 (nhiệt độ lối ra chất tải nhiệt 320 °C, áp suất vòng sơ cấp 15,7 MPa).

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ thế giới quan thực chứng khách quan ($\text{Positivism}$) trong vật lý tính toán hạt nhân, sử dụng phương pháp mô phỏng ngẫu nhiên Monte Carlo năng lượng liên tục kết hợp giải tích ma trận chuỗi cháy vi mô.

Thiết kế đa cấp độ ($\text{Multi-level hierarchical design}$):

  1. Cấp độ vi mô (Micro-scale): Mô hình hóa các vi hạt $\text{Gd}_2\text{O}_3$ hình cầu đơn lẻ có đường kính từ 40 µm, 60 µm, 80 µm, 100 µm đến 300 µm phân bố ngẫu nhiên trong chất nền $\text{UO}_2$.
  2. Cấp độ thanh nhiên liệu (Pin-cell scale): Viên nhiên liệu bán kính 3,8 mm, lỗ rỗng trung tâm bán kính 1,15 mm, khe hở khí Heli, vỏ bọc hợp kim Zr-1%Nb ($\text{E110}$) bán kính ngoài 4,55 mm, đặt trong ô mạng tam giác bước mạng 12,75 mm.
  3. Cấp độ bó nhiên liệu (Assembly scale): Bó nhiên liệu lục giác 331 vị trí, khoảng cách tâm hai bó liền kề 23,6 cm, gồm 312 thanh nhiên liệu (độ giàu $^{235}\text{U}$ đồng nhất 3,7 wt% hoặc 4,4 wt%), 18 ống dẫn thanh điều khiển CPSAR và 1 ống đo trung tâm.

Quy trình nghiên cứu rigorous

Quy trình tính toán vận hành theo chuỗi liên hoàn:

  • Dữ liệu tiết diện phản ứng vi mô liên tục được trích xuất từ thư viện $\text{JENDL-3.3}$ thông qua hệ mã xử lý $\text{LICEM}$.
  • Xử lý vùng cộng hưởng chưa giải tỏa bằng phương pháp bảng xác suất ($\text{probability table method}$).
  • Mô hình hình học ngẫu nhiên ($\text{STGM}$) trong mã MVP xử lý phân bố vi hạt theo hàm phân bố lân cận gần nhất $\text{NND}(r)$ được tạo bởi chương trình $\text{MCRDF}$: $$\frac{d}{dr}(\text{NND}(r)) = \frac{3 f_p}{2 (1 - f_p)} \exp\left(-\frac{3 f_p}{2 (1 - f_p)} r\right)$$ với $f_p$ là tỷ lệ điền đầy thể tích vi hạt ($3,33% - 5,0%$).
  • Quá trình cháy hạt nhân được giải bằng $\text{MVP_BURN}$ với mô hình chuỗi cháy chuẩn u4cm6fp50bp16T gồm các actinides chính và 50 sản phẩm phân hạch ($\text{FPs}$), 16 chất hấp thụ nơtrôn độc hại.

Độ tin cậy và kiểm chứng chéo ($\text{Triangulation & Validation}$): Kết quả mô phỏng từ chương trình MVP được đối chứng độc lập với mã tính toán chuẩn công nghiệp $\text{SRAC}$ (Standard Reactor Analysis Code) và bộ dữ liệu đối chuẩn quốc tế ($\text{Benchmark - BM}$) của OECD/NEA [13]. Sai số thống kê Monte Carlo của hệ số nhân vô hạn $\sigma(k_\infty)$ được kiểm soát nghiêm ngặt ở mức dưới $\pm 0,05%$ ($< 50\text{ pcm}$) với số lượng lịch sử hạt $10^6 - 10^7$ nơtrôn trên mỗi điểm cháy.

Data và phân tích

Thông số vật lý hạt nhân Giá trị định lượng Ý nghĩa vật lý / Vai trò trong mô hình
Tiết diện hấp thụ $^{157}\text{Gd}$ ($\sigma_a$ tại 0,0253 eV) 252.912 barn Đồng vị hấp thụ nơtrôn nhiệt mạnh nhất, chi phối giai đoạn 0–5 GWd/t
Tiết diện hấp thụ $^{155}\text{Gd}$ ($\sigma_a$ tại 0,0253 eV) 60.737 barn Đồng vị hấp thụ mạnh thứ hai, duy trì kiểm soát độ phản ứng đến 10 GWd/t
Quãng đường tự do trung bình $\lambda_t$ của nơtrôn trong $\text{Gd}_2\text{O}_3$ 8,3 µm Cơ sở định cỡ hạt: hạt $> 8,3\text{ µm}$ xuất hiện hiệu ứng tự che chắn mạnh
Độ dẫn nhiệt hạt vi mô $\text{UO}_2\text{–Gd}_2\text{O}_3$ (300–1273 K) 5,8 – 2,7 W/m·K Cao hơn vượt trội so với hỗn hợp trộn đều (3,8 – 2,6 W/m·K) [25, 32]
Điểm nóng chảy của $\text{Gd}_2\text{O}_3$ 2350 °C Đảm bảo tính toàn vẹn cấu trúc trong điều kiện vận hành nhiệt độ cao
Số vị trí trong bó nhiên liệu VVER-1000 331 vị trí 312 thanh nhiên liệu, 18 ống CPSAR, 1 ống đo trung tâm
Cấu hình nghiên cứu 1 12 thanh IBA, hạt 60 µm Tỷ lệ thể tích $\text{Gd}_2\text{O}_3$ là 5,0%, so sánh với mô hình trộn đều chuẩn
Cấu hình nghiên cứu 2 (Đột phá) 18 thanh IBA, hạt 300 µm Tỷ lệ thể tích $\text{Gd}_2\text{O}_3$ là 3,33%, bố trí đối xứng tối ưu hóa PPF

Các kiểm tra độ vững chắc ($\text{Robustness Checks}$) được thực hiện qua các phép tính phân nhánh nhiệt độ và biến thiên nồng độ Boron tải nhiệt từ 0 ppm, 50% Boron danh định đến 100% Boron danh định.

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Khả năng thay thế tương đương của cấu hình 12 thanh IBA vi hạt 60 µm: Mô phỏng chỉ ra rằng vi hạt $\text{Gd}2\text{O}3$ đường kính 60 µm (hàm lượng thể tích 5,0%) tạo ra đường cong biến thiên hệ số nhân vô hạn $k\infty$ theo độ cháy sâu gần như trùng khớp hoàn toàn với thiết kế trộn đều truyền thống trong dải 0–10 GWd/t. Độ lệch $\Delta k\infty$ giữa hai mô hình tại mọi bước cháy đều nhỏ hơn 150 pcm. Hai đồng vị $^{155}\text{Gd}$ và $^{157}\text{Gd}$ cháy kiệt hoàn toàn sau 10 GWd/t, trả lại đặc tính nơtrôn học bình thường cho bó nhiên liệu mà không để lại phản ứng dư âm.

  2. Cấu hình cải tiến 18 thanh IBA vi hạt 300 µm làm phẳng phân bố công suất: Khi phân tán 18 thanh nhiên liệu chứa $\text{Gd}_2\text{O}_3$ vi hạt kích thước 300 µm (tỷ lệ thể tích 3,33% - Mô hình 1 và Mô hình 2), hiệu ứng tự che chắn mạnh ở hạt 300 µm làm giảm tiết diện hấp thụ hiệu dụng ban đầu trên từng thanh, giúp loại bỏ hiện tượng lõm sâu cục bộ của công suất tại các thanh lân cận thanh IBA.

  3. Cắt giảm đột phá hệ số đỉnh công suất (PPF): So với bó nhiên liệu tham chiếu (12 thanh trộn đều có $\text{PPF}_{\text{max}} = 1,32$ tại BOC), cấu hình 18 thanh vi hạt 300 µm hạ thấp $\text{PPF}$ xuống dưới 1,18 tại 0 GWd/t và duy trì $\text{PPF} < 1,15$ trong toàn bộ dải cháy từ 0 đến 40 GWd/t. Điều này làm tăng đáng kể biên độ an toàn nhiệt (DNBR margin).

  4. Vận hành an toàn ở nồng độ Boron thấp và duy trì MTC âm sâu: Trong điều kiện nồng độ Boron hòa tan giảm 50% hoặc hạ về 0 ppm, bó nhiên liệu 18 thanh IBA vi hạt 300 µm vẫn kiểm soát hoàn hảo độ phản ứng dự trữ. Đặc biệt, hệ số phản hồi nhiệt độ chất làm chậm ($\text{MTC}$) luôn duy trì giá trị âm sâu (từ $-35\text{ pcm/}^\circ\text{C}$ đến $-55\text{ pcm/}^\circ\text{C}$) trong suốt chu trình cháy, triệt tiêu hoàn toàn rủi ro xuất hiện MTC dương ở đầu chu trình.

  5. Bảo tồn trọn vẹn đặc tính truyền nhiệt gốm nhiên liệu: Việc ứng dụng vi hạt $\text{Gd}_2\text{O}_3$ giải phóng chất nền $\text{UO}_2$ khỏi tình trạng hòa tan mạng tinh thể của ion $\text{Gd}^{3+}$, bảo toàn độ dẫn nhiệt của thanh nhiên liệu ở mức $5,8 - 2,7\text{ W/m}\cdot\text{K}$, giảm nhiệt độ tâm thanh nhiên liệu từ 100–150 °C ở công suất đỉnh.

[So sánh hệ số đỉnh công suất PPF theo độ cháy sâu]
PPF
1.40 |  x (Bó 12 thanh trộn đều tham chiếu: PPF_max = 1.32)
1.30 |  \
1.20 |   \-----x (Bó 12 thanh hạt 60µm: PPF ~ 1.25)
1.10 |        \========================o (Bó 18 thanh hạt 300µm cải tiến: PPF < 1.18)
1.00 +--------------------------------------------------> Độ cháy (GWd/t)
     0        5        10       20       30       40

Implications đa chiều

  • Lý thuyết: Đặt nền móng cho lý thuyết thiết kế nhiên liệu hạt nhân cấu trúc vi mô dị thể ($\text{Micro-heterogeneous Nuclear Fuel Design}$), chứng minh rằng việc thao túng cấu trúc hình học ở thang micromet có thể điều khiển chính xác các đặc trưng nơtrôn học vĩ mô ở thang mét của vùng hoạt.
  • Phương pháp luận: Thiết lập quy trình tính toán đối chuẩn tin cậy cao kết hợp MVP/GMVP II, STGM và MVP_BURN, có thể chuyển giao áp dụng cho các cấu hình lò phản ứng nước áp lực phương Tây (PWR AP1000, EPR) và lò phản ứng module nhỏ (SMR).
  • Thực tiễn công nghiệp: Cung cấp giải pháp kỹ thuật khả thi cho Tập đoàn Năng lượng Nguyên tử Quốc gia Nga (Rosatom) và Công ty Nhiên liệu TVEL để nâng cấp các bó nhiên liệu $\text{TVS-2M}$ cho đội tàu lò VVER-1000 đang vận hành trên toàn cầu.
  • Chính sách an toàn: Cung cấp bằng chứng khoa học thực nghiệm tính toán cho các cơ quan pháp quy hạt nhân (như Cục An toàn Bức xạ và Hạt nhân - VARANS) trong việc thẩm định an toàn thiết kế nhiên liệu hạt nhân dung sai sự cố ($\text{Accident Tolerant Fuel - ATF}$).

Limitations và Future Research

Luận án thẳng thắn thừa nhận 4 giới hạn nghiên cứu cụ thể:

  1. Mô phỏng ở cấp độ bó nhiên liệu đơn lẻ: Nghiên cứu áp dụng điều kiện biên phản xạ trên một bó nhiên liệu 2D/3D vô hạn, chưa tính đến tương tác phổ nơtrôn giữa các bó nhiên liệu có độ giàu và độ cháy khác nhau trong sơ đồ đảo nhiên liệu thực tế của toàn vùng hoạt (Full-core loading pattern).
  2. Chưa ghép tầng nhiệt thủy lực động lực học (Thermal-Hydraulic Feedback Coupling): Các thông số nhiệt độ chất làm chậm và nhiên liệu được cố định theo trạng thái chuẩn, chưa tích hợp mã nhiệt thủy lực đa chiều (như RELAP5 hoặc COBRA-IV) để phản hồi biến thiên tỷ trọng nước tải nhiệt theo chiều cao thanh nhiên liệu.
  3. Giả định phân bố vi hạt đẳng hướng lý tưởng: Mô hình STGM giả định các hạt phân tán ngẫu nhiên đồng đều theo phân bố Poisson/NND, chưa xét tới hiện tượng kết tụ hạt ($\text{particle agglomeration}$) có thể xảy ra trong quá trình thiêu kết gốm $\text{UO}_2$ thực tế.
  4. Thiếu dữ liệu chiếu xạ thực tế trong lò phản ứng nghiên cứu: Kết quả dựa trên mô phỏng số độ chính xác cao; cần các thử nghiệm chiếu xạ trong thanh nhiên liệu thực nghiệm (Lead Test Assemblies) để kiểm chứng mức độ giải phóng khí phân hạch (Fission Gas Release).

Chương trình nghiên cứu tiếp nối khuyến nghị 4 hướng đi:

  • Phát triển mô hình ghép tầng đa vật lý ($\text{Multi-physics Coupling}$) giữa MVP-BURN và mã nhiệt thủy lực mã nguồn mở OpenFOAM.
  • Mở rộng khảo sát sang cấu hình bó nhiên liệu lò VVER-1200 thế hệ III+ với độ giàu 4,95 wt% và chu trình 24 tháng.
  • Nghiên cứu cơ tính và ứng suất nhiệt vi mô tại ranh giới hạt $\text{UO}_2\text{–Gd}_2\text{O}_3$ dưới hiệu ứng chiếu xạ nơtrôn nhanh.
  • Thiết kế thử nghiệm chiếu xạ mẫu nhiên liệu vi hạt tại Lò phản ứng nghiên cứu hạt nhân Đà Lạt hoặc Trung tâm Nghiên cứu Khoa học Công nghệ Hạt nhân (CNST) trong tương lai.

Tác động và ảnh hưởng

  • Ảnh hưởng học thuật: Mở ra một nhánh nghiên cứu chuyên sâu về ứng dụng vi hạt hấp thụ trong vật lý lò phản ứng, với tiềm năng tạo ra các công bố khoa học trên các tạp chí quốc tế đầu ngành nhóm Q1 (như Annals of Nuclear Energy, Progress in Nuclear Energy, Nuclear Engineering and Design).
  • Chuyển đổi công nghiệp: Tăng cường hệ số an toàn nhiệt của các nhà máy điện hạt nhân chạy lò VVER-1000, cho phép các nhà vận hành tăng công suất phát điện định mức thêm 5%–10% mà không cần tái thiết kế cấu trúc thùng lò.
  • Tác động chính sách: Đóng góp luận cứ khoa học quan trọng cho Chiến lược phát triển ứng dụng năng lượng nguyên tử vì mục đích hòa bình, hỗ trợ quá trình đào tạo nguồn nhân lực chất lượng cao trong ngành Vật lý hạt nhân và Kỹ thuật hạt nhân tại Việt Nam.
  • Lợi ích kinh tế - xã hội: Giảm thiểu phát sinh chất thải phóng xạ thông qua việc tăng độ cháy sâu của nhiên liệu lên 50–60 GWd/t, tối ưu hóa chi phí chu trình nhiên liệu hạt nhân hàng triệu USD mỗi năm cho mỗi tổ máy 1000 MWe.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và học viên sau đại học: Tiếp cận phương pháp luận mô phỏng Monte Carlo tiên tiến với mã MVP và kỹ thuật xử lý hình học thống kê ngẫu nhiên STGM.
  • Các nhà khoa học vật lý vùng hoạt: Thu nhận bộ thông số thiết kế tối ưu định lượng (kích thước 60 µm và 300 µm, tỷ lệ thể tích $3,33% - 5,0%$) làm mốc đối chuẩn cho các nghiên cứu cải tiến nhiên liệu LWR.
  • Kỹ sư thiết kế nhiên liệu tại các tập đoàn hạt nhân: Ứng dụng trực tiếp bản thiết kế phân bố 18 thanh IBA vi hạt để cải tiến các thế hệ bó nhiên liệu thương mại cho lò VVER.
  • Cơ quan quản lý an toàn bức xạ và hạt nhân: Sử dụng các dữ liệu về hệ số phản hồi MTC và hệ số đỉnh công suất PPF để xây dựng tiêu chí đánh giá giấy phép cho các thiết kế vùng hoạt sử dụng nhiên liệu tiên tiến.

Câu hỏi chuyên sâu

  1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào? Luận án đã mở rộng Lý thuyết Tự che chắn Không gian trong Môi trường Hạt Dị thể Ngẫu nhiên bằng cách giải quyết thành công phương trình vận chuyển nơtrôn Boltzmann năng lượng liên tục qua cấu trúc hạt vi mô $\text{Gd}_2\text{O}_3$ phân tán. Luận án chứng minh rằng thông số bán kính hạt vi mô $R_0$ đóng vai trò như một bậc tự do độc lập cho phép điều khiển tốc độ cháy của chất hấp thụ mà không làm thay đổi hàm lượng nạp chất hấp thụ ban đầu.

  2. Cải tiến phương pháp luận so với các nghiên cứu trước đây thể hiện như thế nào? So với nghiên cứu của Reda et al. [52] (chỉ tăng số thanh trộn đều, làm giảm truyền nhiệt) và Choe et al. [10] (phủ lớp $\text{ZrEr}_2$, làm giảm lượng Uranium nạp), luận án ứng dụng mô hình hình học thống kê STGM kết hợp thuật toán láng giềng gần nhất NND trong mã MVP. Phương pháp này cho phép mô phỏng chính xác sự chồng lấn không gian của hàng triệu vi hạt ngẫu nhiên mà không cần đơn giản hóa ô mạng, điều mà các mã tất định truyền thống như SRAC hay CITATION không thể thực hiện được.

  3. Phát hiện bất ngờ và phản trực giác nhất trong dữ liệu là gì? Phát hiện phản trực giác nhất là việc tăng kích thước vi hạt $\text{Gd}_2\text{O}_3$ lên 300 µm trong cấu hình 18 thanh lại giúp hệ số đỉnh công suất PPF giảm mạnh xuống dưới 1,18 (thay vì tăng lên do lo ngại tính bất đồng nhất nơtrôn). Hiệu ứng tự che chắn cực mạnh của hạt 300 µm đã kìm hãm sự sụt giảm thông lượng nơtrôn nhiệt tại các thanh nhiên liệu lân cận, tạo nên một phân bố mật độ phân hạch đồng đều hơn cả thiết kế trộn đều truyền thống.

  4. Nghiên cứu có cung cấp giao thức tái lập (Replication Protocol) không? Có. Luận án mô tả chi tiết toàn bộ thông số hình học ô mạng VVER-1000 (bước mạng 12,75 mm, bán kính nhiên liệu 3,8 mm, vỏ hợp kim E110), thành phần đồng vị từ thư viện JENDL-3.3, quy chuẩn chuỗi cháy u4cm6fp50bp16T, phương pháp sinh số ngẫu nhiên Monte Carlo và các cấu hình tọa độ thanh IBA trong bó nhiên liệu, cho phép các nhóm nghiên cứu độc lập tái lập hoàn toàn các kết quả mô phỏng.

  5. Chương trình nghiên cứu 10 năm được phác thảo như thế nào? Chương trình 10 năm bao gồm: Giai đoạn 1 (Năm 1–3): Ghép tầng đa vật lý nơtrôn - nhiệt thủy lực 3D toàn vùng hoạt; Giai đoạn 2 (Năm 4–6): Nghiên cứu cơ tính và công nghệ thiêu kết vi hạt $\text{UO}_2\text{–Gd}_2\text{O}_3$ trong điều kiện phòng thí nghiệm; Giai đoạn 3 (Năm 7–10): Chế tạo thanh nhiên liệu mẫu và thực hiện thử nghiệm chiếu xạ in-pile tại lò phản ứng nghiên cứu quốc gia.

Kết luận

Công trình nghiên cứu của NCS Hoàng Thanh Phi Hùng đã hoàn thành xuất sắc các mục tiêu khoa học đề ra với các đóng góp cốt lõi:

  1. Xác lập tính khả thi nơtrôn học của việc sử dụng vi hạt $\text{Gd}_2\text{O}_3$ trong thanh nhiên liệu $\text{UO}_2$ bó nhiên liệu VVER-1000/V-320, giải quyết dứt điểm mâu thuẫn giữa hiệu suất hấp thụ nơtrôn và độ suy giảm truyền nhiệt gốm oxit.
  2. Xác định kích thước vi hạt tối ưu 60 µm cho thiết kế 12 thanh IBA truyền thống, tạo ra đường cong $k_\infty$ trùng khớp với mô hình trộn đều chuẩn trong dải cháy 0–10 GWd/t.
  3. Sáng tạo thiết kế cải tiến mới gồm 18 thanh IBA chứa vi hạt $\text{Gd}_2\text{O}_3$ đường kính 300 µm (tỷ lệ thể tích 3,33%), giúp cắt giảm hệ số đỉnh công suất PPF từ 1,32 xuống dưới 1,18 trong toàn bộ chu trình cháy.
  4. Chứng minh khả năng vận hành an toàn của bó nhiên liệu cải tiến trong môi trường chất tải nhiệt nồng độ Boron thấp (50% và 0 ppm), đảm bảo hệ số nhiệt độ chất làm chậm MTC luôn âm sâu.
  5. Thiết lập quy trình tính toán ngẫu nhiên Monte Carlo năng lượng liên tục chuẩn mực với mã MVP/GMVP II và mô hình STGM cho các bài toán nhiên liệu hạt nhân cấu trúc vi mô phân tán dị thể.
  6. Mở ra 3 hướng nghiên cứu mới về nhiên liệu hạt nhân dung sai sự cố (ATF), mô phỏng ghép tầng đa vật lý nơtrôn - nhiệt thủy lực, và tối ưu hóa vùng hoạt các thế hệ lò phản ứng tiên tiến VVER-1200 và SMR.