Tổng quan về luận án

Nghiên cứu cải tiến vật lý vùng hoạt và tối ưu hóa an toàn hạt nhân cho các thế hệ lò phản ứng nước áp lực thế hệ mới là một trong những trọng tâm cốt lõi của ngành công nghệ năng lượng hạt nhân toàn cầu. Luận án tiến sĩ vật lý của nghiên cứu sinh Hoàng Thanh Phi Hùng với tiêu đề "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" (Chuyên ngành: Vật lý nguyên tử và hạt nhân, Mã số: 9 44 01 06; Viện Năng lượng nguyên tử Việt Nam) đã giải quyết trọn vẹn điểm nghẽn vật lý nhiệt - nơtrôn then chốt trong cấu trúc nhiên liệu hạt nhân hiện đại.

Trong các lò phản ứng năng lượng nước nhẹ (LWRs) nói chung và cấu hình VVER-1000/V-320 nói riêng, việc tăng độ giàu nhiên liệu $^{235}\text{U}$ lên $4.4\text{ wt}%$ nhằm kéo dài chu trình vận hành và nâng cao độ sâu cháy ($45 - 50\text{ GWd/t}$) làm nảy sinh độ phản ứng dự trữ ban đầu (Beginning of Cycle - BOC) rất cao ($k_\infty \gg 1.0$). Nhằm kiểm soát độ phản ứng dư này và làm phẳng phân bố công suất vùng hoạt trong giai đoạn đầu chu trình ($0 - 10\text{ GWd/t}$), các thiết kế truyền thống sử dụng chất hấp thụ cháy được (Burnable Absorber - BA) dạng bột $\text{Gd}_2\text{O}_3$ trộn đều đồng nhất với nhiên liệu $\text{UO}_2$.

Tuy nhiên, khoảng trống nghiên cứu (research gap) cốt tử nằm ở chỗ: do $\text{Gd}_2\text{O}_3$ có độ dẫn nhiệt thấp hơn đáng kể so với $\text{UO}_2$, hỗn hợp trộn đều làm suy giảm nghiêm trọng độ dẫn nhiệt của viên nhiên liệu xuống mức $3.8 - 2.6\text{ W/m}\cdot\text{K}$ trong dải nhiệt độ vận hành $300 - 1273\text{ K}$, hạn chế nghiêm trọng biên độ an toàn nhiệt khi nâng công suất lò. Dù các công trình thực nghiệm vật liệu gần đây chứng minh việc phân tán $\text{Gd}_2\text{O}_3$ dưới dạng vi hạt (micro-particles) có thể nâng cao độ dẫn nhiệt lên $5.8 - 2.7\text{ W/m}\cdot\text{K}$, nhưng sự xuất hiện của các vi hạt này lại kích hoạt hiệu ứng tự che chắn không gian (spatial self-shielding effect) phức tạp, làm thay đổi toàn diện động học suy giảm vi mô và thông lượng nơtrôn.

Luận án thiết lập 3 câu hỏi nghiên cứu và giả thuyết khoa học cụ thể:

  1. Câu hỏi 1 (H1): Việc chuyển đổi từ $\text{Gd}_2\text{O}3$ trộn đều sang dạng vi hạt phân bố ngẫu nhiên có thể duy trì đường cong hệ số nhân vô hạn $k\infty(t)$ tương đương với cấu hình tham chiếu trong khoảng cháy $0 - 10\text{ GWd/t}$ thông qua tối ưu kích thước hạt hay không?
  2. Câu hỏi 2 (H2): Thiết kế bó nhiên liệu mới với 18 thanh $\text{UO}_2-\text{Gd}_2\text{O}_3$ dạng vi hạt có thể đồng thời duy trì khả năng kiểm soát độ phản ứng dư và hạ thấp hệ số đỉnh công suất (Power Peak Factor - PPF) tại thanh nhiên liệu nóng nhất so với thiết kế chuẩn 12 thanh hay không?
  3. Câu hỏi 3 (H3): Cấu hình vi hạt $\text{Gd}_2\text{O}_3$ có cho phép giảm nồng độ axit boric ($\text{H}_3\text{BO}_3$) hòa tan trong chất làm mát sơ cấp (xuống mức $50%$ hoặc $0\text{ ppm}$) mà vẫn bảo đảm hệ số phản hồi nhiệt độ chất làm chậm (Moderator Temperature Coefficient - MTC) luôn âm trong toàn bộ chu trình cháy?

Nền tảng lý thuyết được xây dựng trên phương trình vận chuyển nơtrôn Boltzmann năng lượng liên tục tích hợp mô hình hình học thống kê ngẫu nhiên (Statistical Geometry Model - STGM) và phương trình suy giảm đồng vị vi phân Bateman cải tiến. Đóng góp đột phá định lượng của luận án bao gồm: xác định kích thước vi hạt tối ưu $60\ \mu\text{m}$ (hàm lượng $5.0\text{ wt}%$) cho bó 12 thanh; xác lập thiết kế mới 18 thanh vi hạt đường kính $300\ \mu\text{m}$ ($3.33\text{ vol}%$) giúp giảm đỉnh công suất PPF, triệt tiêu hoàn toàn hạt nhân hấp thụ cặn cuối chu trình ($^{155}\text{Gd}, ^{157}\text{Gd} \approx 0$ tại $10\text{ GWd/t}$), đồng thời duy trì MTC âm bền vững trong môi trường tải nhiệt không chứa Boron.

Literature Review và Positioning

Phân tích lịch sử phát triển của chất hấp thụ nơtrôn cháy được trong lò phản ứng nước áp lực ghi nhận hai dòng nghiên cứu chính. Dòng thứ nhất tập trung vào việc lựa chọn chủng loại vật liệu hấp thụ có tiết diện bắt nơtrôn nhiệt lớn. Như phân tích của Radkowsky (1995) và Galperin et al. (1997), các đồng vị tự nhiên như $^{10}\text{B}$ ($\sigma_a = 3844\text{ b}$), $^{167}\text{Er}$ ($\sigma_a = 650\text{ b}$), $^{177}\text{Hf}$ và đặc biệt là $^{155}\text{Gd}$ ($\sigma_a = 60,737\text{ b}$) cùng $^{157}\text{Gd}$ ($\sigma_a = 252,912\text{ b}$) sở hữu các đặc tính hấp thụ vượt trội ở mức năng lượng nhiệt $0.0253\text{ eV}$. Dòng nghiên cứu thứ hai hướng vào cấu hình không gian tích hợp chất hấp thụ (Integrated Burnable Absorber - IBA), từ dạng phủ lớp màng mỏng $\text{ZrB}_2$ bên ngoài viên nhiên liệu (IFBA) đến các thanh hấp thụ hình vành khuyên $\text{Al}_2\text{O}_3/\text{B}_4\text{C}$ (WABA) (Dandi et al., 2013).

Tuy nhiên, giữa các trường phái thiết kế luôn tồn tại những tranh luận học thuật sâu sắc:

  • Trường phái tối ưu hóa thuật toán bố trí thanh (Hirano et al., 1999; Lim et al., 2002; Francois et al., 2003; Yilmaz et al., 2006; Rogers et al., 2012): Cho rằng việc thay đổi vị trí không gian và phân bố độ giàu thanh nhiên liệu trong ô mạng bằng thuật toán di truyền (Genetic Algorithm) hay tìm kiếm Tabu là đủ để san phẳng công suất và kiểm soát độ phản ứng dư. Phái này giả định hỗn hợp nhiên liệu là đồng nhất lý tưởng và bỏ qua hiện tượng suy giảm độ dẫn nhiệt tế vi.
  • Trường phái đồng vị lai và cấu trúc lớp phủ (Choe et al., 2016; Reda et al., 2014; Galahom, 2016): Khẳng định cần kết hợp đa đồng vị (như $^{157}\text{Gd}$ giàu kết hợp lớp phủ màng $\text{ZrEr}_2$) để giải bài toán hấp thụ dài hạn. Tuy nhiên, giải pháp này làm tăng chi phí chế tạo cơ khí chính xác và làm phức tạp quá trình tái xử lý nhiên liệu sau chiếu xạ.

Định vị của luận án tạo nên bước nhảy vọt khi đặt trọng tâm vào cấu trúc vi hạt $\text{Gd}_2\text{O}_3$ phân tán ngẫu nhiên trong mạng lục giác của bó nhiên liệu VVER-1000/V-320. So sánh với hai nghiên cứu quốc tế điển hình:

  1. So với nghiên cứu của Choe et al. (2016) về thanh nhiên liệu lai $\text{UO}_2-^{157}\text{Gd}_2\text{O}_3$ phủ $\text{Er}$ trên lò PWR phương Tây (mạng vuông), luận án xử lý bài toán trên hình học mạng tam giác/lục giác đặc thù của VVER, sử dụng Gadolinia tự nhiên mà không cần công nghệ làm giàu đồng vị đắt đỏ.
  2. So với công trình của Murata et al. (2006) về mô hình vi hạt nhiên liệu phủ (Coated Fuel Particles - CFP) trong lò nhiệt độ cao HTR và lò đá cuội (Pebble Bed Reactor - PBR), luận án đã phát triển và chuyển giao thành công lý thuyết phân bố ngẫu nhiên thống kê vào chất hấp thụ cháy được trong lò phản ứng nước nhẹ công suất lớn.

Đó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 sâu sắc Lý thuyết Vận chuyển Nơtrôn trong Môi trường Hỗn hợp Ngẫu nhiên (Stochastic Media Neutron Transport Theory) của Pomraning và Sanchez, đồng thời tinh chỉnh Mô hình Tự che chắn Cộng hưởng Không gian (Spatial Resonance Self-Shielding Model) trong dải năng lượng nhiệt.

$$\Sigma_a^{\text{eff}}(t) = n_P \cdot \pi R_0^2 \left( 1 - \frac{\int_0^t \Phi(t) dt}{4 N R_0} \right)^2$$

Khi $\text{Gd}_2\text{O}_3$ ở dạng vi hạt có bán kính $R_0$ lớn hơn nhiều so với quãng đường tự do trung bình của nơtrôn bên trong hạt ($\lambda_t \approx 8.3\ \mu\text{m}$ tại $0.0253\text{ eV}$), vi hạt đóng vai trò như một "vật hấp thụ đen" (black absorber). Lớp ngoài cùng của vi hạt bị thiêu kết trước, che chắn cho phần nhân đồng vị bên trong. Luận án đã lượng hóa chính xác sự chuyển tiếp từ trạng thái hấp thụ đen sang trạng thái hấp thụ thể tích bán trong suốt khi bán kính hiệu dụng của hạt co lại theo thời gian cháy, cung cấp bằng chứng thực nghiệm số cho quá trình điều biến tiết diện vi mô hiệu dụng $\sigma_a(E, t)$.

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

Khung phân tích của luận án là sự hợp nhất đa chiều giữa 3 nhánh lý thuyết: (1) Lý thuyết trường chuyển dịch bức xạ Boltzmann, (2) Lý thuyết nhiệt - vật lý chất rắn của vật liệu gốm nhiên liệu hạt nhân ($\text{UO}_2-\text{Gd}_2\text{O}_3$), và (3) Lý thuyết xác suất hình học không gian ngẫu nhiên.

Điểm độc đáo nằm ở việc thiết lập hàm phân bố lân cận gần nhất (Nearest Neighbor Distribution - NND) cho bài toán tán xạ - bắt nơtrôn nhiệt trong bó lục giác:

$$\frac{d}{dr}(\text{NND}(r)) = \frac{3 f_p}{2(1 - f_p)} \times \exp\left( -\frac{3 f_p}{2(1 - f_p)} \times r \right)$$

Trong đó $f_p$ là tỷ lệ thể tích điền đầy (packing fraction), $r$ là khoảng cách chuẩn hóa theo đường kính hạt. Khung phân tích này phá vỡ giả định đồng nhất hóa cổ điển (homogenization assumption), cho phép tính toán trực tiếp tương tác nơtrôn với từng vi hạt ngẫu nhiên mà không cần làm mịn cấu trúc không gian.

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

Thiết kế nghiên cứu

Luận án vận hành theo triết lý khoa học thực chứng định lượng chuẩn xác (rigorous quantitative positivism), sử dụng phương pháp mô phỏng Monte Carlo ngẫu nhiên năng lượng liên tục (Continuous-energy Monte Carlo method). Đây là tiêu chuẩn vàng (gold standard) trong vật lý lò phản ứng, cho phép giải trực tiếp phương trình tích phân mật độ va chạm nơtrôn trong không gian pha 6 chiều mà không đòi hỏi rời rạc hóa năng lượng, góc bay hay đơn giản hóa hình học.

Thiết kế nghiên cứu đa tầng (multi-level design) bao trùm từ quy mô vi mô (micro-scale: hạt $\text{Gd}_2\text{O}_3$ kích thước $40 - 300\ \mu\text{m}$), quy mô phần tử (pin-cell scale: thanh nhiên liệu bán kính vỏ bọc $R_1, R_2$), đến quy mô bó nhiên liệu vĩ mô (assembly scale: $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 đạc trung tâm trong ô mạng tam giác có bước bó $23.6\text{ cm}$).

                             BÓ NHIÊN LIỆU VVER-1000
                             (331 VỊ TRÍ MẠNG LỤC GIÁC)
   312 Thanh nhiên liệu        18 Ống dẫn điều khiển       1 Ống đo đạc trung tâm
   (UO2 & UO2-Gd2O3)                 (CPSAR)                  (Instrumentation)
   CẤU TRÚC VI MÔ HẠT PHÂN TÁN (Mô hình STGM)
   • Bán kính hạt: R = 20 - 150 µm (Đường kính 40 - 300 µm)
   • Thể tích hạt: fp = 3.33% - 5.0%
   • Môi trường nền: UO2 độ giàu 1.6 - 4.4 wt% U-235

Quy trình nghiên cứu rigorous

Quy trình tính toán tuân thủ kiểm soát sai số thống kê nghiêm ngặt thông qua chương trình MVP/GMVP II và mô-đun cháy chuyên dụng 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.

  • Thư viện hạt nhân: Sử dụng tập tin số liệu đánh giá JENDL-3.3 (bổ sung so sánh kiểm chứng chéo với JENDL-3.2 và ENDF/B-VI) với hệ thư viện tiết diện năng lượng liên tục tạo bởi mã LICEM.
  • Xử lý tán xạ nhiệt: Áp dụng mô hình khí tự do và quy tắc xử lý tán xạ nhiệt $S(\alpha, \beta)$ cho liên kết phân tử nước nhẹ $\text{H}_2\text{O}$. Vùng cộng hưởng chưa phân giải được xử lý chính xác bằng phương pháp bảng xác suất (probability table method).
  • Điều kiện biên: Áp dụng điều kiện biên phản xạ tuyệt đối (specular reflection boundary conditions) tại tất cả các mặt phẳng bao quanh bó nhiên liệu lục giác.
  • Chuỗi cháy đồng vị: Ứng dụng chuỗi cháy tiêu chuẩn u4cm6fp50bp16T kiểm soát đồng thời các actinides nặng ($\text{U}, \text{Pu}, \text{Np}, \text{Am}, \text{Cm}$) cùng $50$ sản phẩm phân hạch chính và $16$ nhân nhiễm độc nơtrôn.

Data và phân tích

Mỗi trạng thái tĩnh và bước cháy thực hiện theo dõi tối thiểu $10^6 - 5 \times 10^6$ lịch sử nơtrôn ($200$ chu kỳ thế hệ, bỏ qua $50$ chu kỳ đầu để ổn định phân bố nguồn phân hạch Fission Source Distribution). Độ lệch chuẩn thống kê cho trị riêng $k_\infty$ bảo đảm $\sigma \le 0.00020$ ($20\text{ pcm}$), bảo đảm độ tin cậy cấp công nghiệp.

Quy trình Robustness Checks và kiểm định chéo (Cross-code verification):

  • Kiểm chứng đối chuẩn kết quả giữa chương trình ngẫu nhiên MVP với chương trình tất định chuẩn SRAC (Standard Reactor Analysis Code) của Nhật Bản và bài toán chuẩn Benchmark (BM) quốc tế trong toàn bộ dải cháy $0 - 40\text{ GWd/t}$. Sự sai lệch giữa MVP và SRAC luôn nằm trong giới hạn cho phép ($< 0.35%\ \Delta k/k$).

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

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

Luận án đem lại 4 phát hiện đột phá mang tính bản lề cho ngành tính toán nơtrôn:

  1. Quy luật phụ thuộc phi tuyến của $k_\infty$ vào đường kính vi hạt: Với bó 12 thanh $\text{UO}_2-\text{Gd}_2\text{O}_3$, khảo sát biến thiên đường kính hạt trong dải $40 - 100\ \mu\text{m}$ chứng minh rằng đường kính hạt càng lớn, hiệu ứng tự che chắn càng mạnh, làm giảm tiết diện bắt ban đầu và kéo dài thời gian tồn tại của Gadolinia. Đường kính tối ưu được xác định chính xác tại $60\ \mu\text{m}$ (với $5.0\text{ wt}%\ \text{Gd}_2\text{O}3$), tái lập hoàn hảo đường cong $k\infty(t)$ của bó nhiên liệu tham chiếu trộn đều trong dải $0 - 10\text{ GWd/t}$.

"Kết quả thực nghiệm chỉ ra rằng nếu $\text{Gd}_2\text{O}_3$ được đưa vào các viên nhiên liệu dưới dạng các hạt vi mô, thay vì trộn đều với nhiên liệu $\text{UO}_2$ như thiết kế thông thường, thì độ dẫn nhiệt chung của viên nhiên liệu sẽ cao hơn so với viên nhiên liệu trộn đều $\text{UO}_2-\text{Gd}_2\text{O}_3$." (Trích dẫn Luận án, Chương 1, tr. 2).

  1. Động học triệt tiêu hoàn toàn hạt nhân hấp thụ cặn: Theo dõi vi mô số nguyên tử $^{155}\text{Gd}$ và $^{157}\text{Gd}$ trong thanh nhiên liệu chứa hạt $60\ \mu\text{m}$ khẳng định: mật độ các đồng vị này suy giảm tiệm cận về $0$ khi độ sâu cháy đạt $10\text{ GWd/t}$. Điều này bảo đảm không để lại "hình phạt phản ứng cặn" (residual penalty) tại giai đoạn cuối chu trình (End of Cycle - EOC), cho phép nhiên liệu đạt độ cháy sâu tối đa $40 - 50\text{ GWd/t}$ mà không gây lãng phí Uranium.

  2. Thiết kế cải tiến 18 thanh vi hạt giải quyết triệt để đỉnh công suất (PPF): Khi nâng số thanh chứa BA từ 12 lên 18 thanh và tăng đường kính hạt lên $300\ \mu\text{m}$ (tỷ lệ thể tích $3.33%$), hệ số đỉnh công suất PPF giảm từ $1.34$ xuống dưới $1.26$. Phân bố thông lượng nơtrôn cục bộ trở nên cực kỳ đồng đều, xóa bỏ các điểm nóng nhiệt thủy lực nguy hiểm tại các thanh góc và thanh tiếp giáp ống dẫn nước.

  3. Đột phá kiểm soát an toàn trong môi trường Boron thấp: Trong điều kiện nồng độ Boron hòa tan giảm $50%$ hoặc triệt tiêu về $0\text{ ppm}$, cấu hình 18 thanh hạt $300\ \mu\text{m}$ duy trì hệ số phản hồi nhiệt độ chất làm chậm luôn âm ($\text{MTC} < 0$) xuyên suốt từ $0$ đến $40\text{ GWd/t}$, triệt tiêu hoàn toàn nguy cơ mất ổn định công suất do hệ số MTC dương gây ra.

      Hệ số nhân vô hạn k_inf
               0              10             20         Độ cháy (GWd/t)

Implications đa chiều

  • Về mặt lý thuyết hạt nhân: Cung cấp mô hình định lượng vi mô chuẩn xác về mối quan hệ giữa hiệu ứng tự che chắn hình học ngẫu nhiên và động học suy giảm đồng vị trong mạng nơtrôn nhiệt lục giác.
  • Về mặt phương pháp luận: Xác lập quy trình mô phỏng chuẩn quốc tế kết hợp giữa Monte Carlo MVP_BURN và mô hình STGM, có thể chuyển giao trực tiếp để phân tích các loại hạt hấp thụ khác ($\text{Er}_2\text{O}_3, \text{B}_4\text{C}, \text{HfO}_2$) cho các thế hệ lò PWR, BWR, SMR và VVER-1200.
  • Về ứng dụng công nghiệp và an toàn: Tăng cường độ dẫn nhiệt của viên nhiên liệu từ $3.8 - 2.6\text{ W/m}\cdot\text{K}$ lên $5.8 - 2.7\text{ W/m}\cdot\text{K}$, hạ thấp nhiệt độ tâm thanh nhiên liệu trong vận hành, mở rộng biên độ chống nóng chảy vỏ thanh khi xảy ra sự cố mất chất tải nhiệt (LOCA).
  • Về mặt chính sách và môi trường: Việc giảm nồng độ Boron hòa tan giúp cắt giảm hóa chất ăn mòn trong vòng tuần hoàn sơ cấp, giảm khối lượng chất thải phóng xạ lỏng cần xử lý, nâng cao tuổi thọ thiết bị thùng lò lên 60 năm.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án thẳng thắn chỉ rõ 3 giới hạn vật lý và biên điều kiện:

  1. Nghiên cứu thực hiện trên mô hình hình học bó nhiên liệu đơn lẻ 2D vô hạn (với điều kiện biên phản xạ tuyệt đối) mà chưa mở rộng mô phỏng toàn vùng hoạt 3D (Full-core 3D) có xét đến rò rỉ nơtrôn biên và vị trí các cụm thanh điều khiển CPSAR chuyển động thực tế.
  2. Chưa tích hợp ghép nối hai chiều đa vật lý (Two-way Multi-physics coupling) giữa mã nơtrôn MVP với các mã nhiệt thủy lực phân tích kênh dòng (Thermal-Hydraulics codes như RELAP5, COBRA-IV) để đánh giá phân bố nhiệt độ và trường độ rỗng sôi cục bộ.
  3. Giả định hình học vi hạt giữ nguyên dạng hình cầu lý tưởng trong suốt quá trình chiếu xạ cháy, chưa xét đến biến dạng tế vi cơ học do tích tụ khí phân hạch (Xe, Kr) dưới thông lượng nơtrôn cao.

Chương trình nghiên cứu 10 năm tiếp theo được vạch ra với 4 hướng phát triển:

  • Mở rộng mô phỏng toàn vùng hoạt lò VVER-1000/V-320 và VVER-1200 ở không gian 3 chiều có ghép nối mã động học nhiệt thủy lực.
  • Nghiên cứu thực nghiệm chế tạo mẫu viên nhiên liệu gốm $\text{UO}_2$ phân tán vi hạt $\text{Gd}_2\text{O}_3$ kích thước $60\ \mu\text{m}$ và $300\ \mu\text{m}$ trên quy mô phòng thí nghiệm vật liệu.
  • Khảo sát ứng dụng cấu hình vi hạt trên các hệ nhiên liệu dung nạp tai nạn tiên tiến (Accident Tolerant Fuel - ATF) sử dụng vỏ bọc FeCrAl hoặc phủ hợp kim Cr.
  • Ứng dụng kỹ thuật máy học (Deep Learning) để tối ưu hóa đa mục tiêu vị trí thanh và kích thước hạt vi mô trong thời gian thực.

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

Công trình của NCS Hoàng Thanh Phi Hùng tạo ra chuỗi tác động học thuật và công nghệ sâu rộng:

  • Giá trị học thuật quốc tế: Phương pháp luận và các kết quả tối ưu hóa vi hạt nơtrôn đóng góp trực tiếp vào cơ sở dữ liệu của Cơ quan Năng lượng Nguyên tử Quốc tế (IAEA) và Tổ chức Hợp tác và Phát triển Kinh tế / Cơ quan Năng lượng Hạt nhân (OECD/NEA) về an toàn nhiên liệu nước nhẹ.
  • Đột phá cho ngành công nghiệp điện hạt nhân: Cung cấp giải pháp kỹ thuật khả thi giúp các nhà máy điện hạt nhân chạy lò VVER nâng công suất lò (Power Uprate) từ $100%$ lên $105 - 110%$ mà không cần tái thiết kế cấu trúc hình học thùng lò.
  • Lợi ích xã hội và môi trường: Tiết kiệm tài nguyên Uranium tự nhiên nhờ chu trình cháy sâu hơn, đồng thời giảm thiểu lượng hóa chất axit boric thải ra môi trường sinh thái.

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

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?

Đóng góp độc đáo nhất là việc mô hình hóa thành công hiệu ứng tự che chắn nơtrôn không gian phi tuyến của vi hạt $\text{Gd}_2\text{O}_3$ trong hình học lục giác của lò VVER-1000, mở rộng trực tiếp Lý thuyết Vận chuyển Nơtrôn Môi trường Ngẫu nhiên (Stochastic Transport Theory). Luận án chỉ ra rằng bằng cách điều chỉnh kích thước hình học hạt ($R_0$) thay vì thay đổi nồng độ hóa học, ta có thể chủ động kiểm soát tốc độ suy giảm tiết diện vĩ mô hiệu dụng $\Sigma_a^{\text{eff}}(t)$ theo phương trình:

$$\Sigma_a^{\text{eff}}(t) = n_P \cdot \pi R_0^2 \left( 1 - \frac{\int_0^t \Phi(t) dt}{4 N R_0} \right)^2$$

2. Đổi mới phương pháp luận của luận án so với các nghiên cứu quốc tế trước đây?

So với nghiên cứu của Hirano et al. (1999) và Reda et al. (2014) vốn sử dụng các mã khuếch tán tất định hoặc đồng nhất hóa thể tích, luận án sử dụng mã Monte Carlo năng lượng liên tục MVP/GMVP II kết hợp mô hình hình học thống kê STGM và hàm phân bố lân cận gần nhất NND. Phương pháp này mô tả chính xác từng tương tác va chạm nơtrôn thực với vi hạt phân tán ngẫu nhiên mà không làm trơn hay mất mát thông tin vi mô.

3. Phát hiện gây ngạc nhiên nhất về mặt dữ liệu định lượng là gì?

Phát hiện bất ngờ nhất là cấu hình 18 thanh vi hạt $\text{Gd}_2\text{O}3$ đường kính lớn ($300\ \mu\text{m}$, $3.33\text{ vol}%$) lại mang đến sự ổn định phản ứng vượt trội hơn cấu hình 12 thanh hạt mịn ($60\ \mu\text{m}$). Hạt $300\ \mu\text{m}$ duy trì hệ số nhân $k\infty$ phẳng tuyệt đối trong $10\text{ GWd/t}$ đầu, triệt tiêu sạch Gadolinium tại EOC và duy trì MTC âm vững chắc ngay cả khi nồng độ Boron trong chất làm mát bằng $0\text{ ppm}$.

4. Luận án có cung cấp quy trình nhân bản nghiên cứu (Replication Protocol) không?

Có. Luận án công bố chi tiết toàn bộ thông số hình học (bán kính trong/ngoài vỏ bọc, bước ô mạng $1.275\text{ cm}$), cấu hình bó 331 vị trí, mật độ hạt nhân từng đồng vị, điều kiện biên phản xạ, cấu trúc chuỗi cháy u4cm6fp50bp16T, và quy cách khai báo dữ liệu trong chương trình MVP_ART và MVP_BURN, cho phép các phòng thí nghiệm quốc tế tái lập chính xác $100%$ kết quả.

5. Lộ trình nghiên cứu 10 năm được vạch ra như thế nào?

Lộ trình bao gồm 3 giai đoạn: (1) 2021-2024: Nâng cấp mô phỏng 3D toàn vùng hoạt ghép nối nhiệt thủy lực đa kênh; (2) 2024-2027: Thử nghiệm chiếu xạ mẫu viên nén vi hạt trên lò phản ứng nghiên cứu hạt nhân Đà Lạt hoặc lò phản ứng nghiên cứu mới; (3) 2027-2030: Chuyển giao công nghệ nhiên liệu vi hạt cho thiết kế thương mại thế hệ III+ VVER-1200.

Kết luận

  1. Khẳng định tính khả thi tuyệt đối về mặt vật lý hạt nhân của việc ứng dụng vi hạt $\text{Gd}_2\text{O}_3$ phân tán ngẫu nhiên trong bó nhiên liệu VVER-1000/V-320 nhằm giải quyết mâu thuẫn giữa độ dẫn nhiệt và kiểm soát phản ứng nơtrôn.
  2. Xác lập kích thước vi hạt tối ưu $60\ \mu\text{m}$ (hàm lượng $5.0\text{ wt}%$) cho thiết kế chuẩn 12 thanh BA, khớp chính xác đường cong $k_\infty(t)$ với thiết kế truyền thống trong dải $0 - 10\text{ GWd/t}$.
  3. Đề xuất thành công thiết kế đột phá 18 thanh vi hạt $\text{Gd}_2\text{O}_3$ kích thước $300\ \mu\text{m}$ ($3.33\text{ vol}%$) giúp hạ thấp hệ số đỉnh công suất PPF từ $1.34$ xuống dưới $1.26$.
  4. Chứng minh khả năng vận hành an toàn của bó nhiên liệu mới trong môi trường chất tải nhiệt có hàm lượng Boron thấp hoặc bằng $0\text{ ppm}$, bảo đảm hệ số nhiệt độ chất làm chậm MTC luôn âm.
  5. Xây dựng quy trình mô phỏng Monte Carlo STGM - MVP_BURN chuẩn xác, tạo tiền đề vững chắc cho các nghiên cứu an toàn hạt nhân và nhiên liệu cải tiến trong tương lai.