Tổng quan về luận án

Nghiên cứu vật lý lò phản ứng hạt nhân hiện đại luôn đặt trọng tâm vào việc nâng cao hiệu suất kinh tế chu trình nhiên liệu thông qua việc tăng độ giàu uranium ($^{235}\text{U}$) và kéo dài độ sâu cháy (high burnup fuel), đồng thời phải củng cố các giới hạn an toàn nhiệt thủy lực nghiêm ngặt. Trong lò phản ứng nước áp lực thế hệ mới nói chung và dòng lò nước nhẹ VVER-1000/V-320 nói riêng, độ giàu nhiên liệu tăng cao dẫn đến sự gia tăng đột biến của độ phản ứng dự trữ tại giai đoạn đầu chu trình (Beginning of Cycle - BOC, dải cháy từ 0 đến $10\text{ GWd/t}$). Để kiểm soát lượng dư phản ứng này mà không làm biến dạng trường phân bố công suất, các chất hấp thụ nơtron cháy được (Burnable Absorbers/Poisons - BA/BP) đóng vai trò then chốt. Luận án tiến sĩ vật lý của nghiên cứu sinh Hoàng Thanh Phi Hùng (2021) 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" đã giải quyết một mâu thuẫn vật lý cơ bản: hỗn hợp nhiên liệu truyền thống trộn đều $\text{UO}_2 - \text{Gd}_2\text{O}_3$ làm suy giảm nghiêm trọng độ dẫn nhiệt của viên nhiên liệu, gây bất lợi cho an toàn nhiệt; trong khi việc chuyển sang cấu trúc vi hạt dị thể lại dẫn đến hiệu ứng tự che chắn không gian phức tạp làm thay đổi tốc độ cháy của các đồng vị hấp thụ.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tế thực nghiệm: "với 10 wt% $\text{Gd}_2\text{O}_3$ dưới dạng vi hạt với đường kính khoảng từ 25 – 53 µm thì dải độ dẫn nhiệt của viên nhiên liệu trong dải nhiệt độ 300 – 1273 K là 5,8 – 2,7 W/mK", trong khi "dải độ dẫn nhiệt của hỗn hợp $\text{UO}_2 – \text{Gd}_2\text{O}_3$ trộn đều biến đổi trong khoảng 3,8 – 2,6 W/mK". Dù cấu trúc vi hạt cải thiện rõ rệt khả năng truyền nhiệt, sự phân bố ngẫu nhiên của hàng triệu vi hạt $\text{Gd}_2\text{O}_3$ tạo ra hiện tượng tự che chắn nơtron kép (double heterogeneity self-shielding), khiến các chương trình tính toán tất định truyền thống dựa trên ô mạng đồng nhất hóa tương đương hoàn toàn mất độ chính xác. Luận án đặt ra ba câu hỏi nghiên cứu và giả thuyết khoa học cụ thể:

  1. Câu hỏi nghiên cứu 1 (RQ1): Cơ chế tự che chắn trong các vi hạt $\text{Gd}_2\text{O}_3$ có kích thước micromet ảnh hưởng như thế nào đến tốc độ cháy của các đồng vị hấp thụ cao $^{155}\text{Gd}$ và $^{157}\text{Gd}$ trong hình học bó nhiên liệu lục giác VVER-1000?
  2. Câu hỏi nghiên cứu 2 (RQ2): Cấu hình thiết kế hình học vi hạt (đường kính hạt, tỷ lệ phần trăm thể tích/khối lượng) và sự bố trí vị trí thanh nhiên liệu nào cho phép tái tạo chính xác đường cong độ phản ứng ($k_\infty$) tương đương bó nhiên liệu chuẩn 12 thanh truyền thống, đồng thời triệt tiêu đỉnh công suất cục bộ?
  3. Câu hỏi nghiên cứu 3 (RQ3): Liệu thiết kế cải tiến với 18 thanh $\text{UO}_2 - \text{Gd}_2\text{O}_3$ vi hạt có khả năng bù đắp độ phản ứng để vận hành lò ở chế độ nồng độ Boron hòa tan thấp (nhằm giữ hệ số nhiệt độ chất làm chậm MTC luôn âm) hay không?
  • Giả thuyết khoa học 1 (H1): Đường kính vi hạt $\text{Gd}_2\text{O}_3$ ($40 - 300\ \mu\text{m}$) điều khiển trực tiếp tốc độ cháy của hạt thông qua sự suy giảm bán kính hình học hiệu dụng theo phương trình $R(t) = R_0 [1 - \frac{\int \Phi(t) dt}{4 N R_0}]$, cho phép kiểm soát chính xác độ sâu cháy của chất hấp thụ trong khoảng $0 - 10\text{ GWd/t}$.
  • Giả thuyết khoa học 2 (H2): Bó nhiên liệu 18 thanh $\text{UO}_2 - \text{Gd}_2\text{O}_3$ vi hạt ($d = 300\ \mu\text{m}$, tỷ lệ thể tích $3,33%$) sẽ tối ưu hóa phân bố thông lượng nơtron, hạ thấp hệ số đỉnh công suất (Power Peak Factor - PPF) so với bó 12 thanh ($d = 60\ \mu\text{m}$, $5,0\text{ wt}%$) và duy trì hệ số phản hồi nhiệt độ chất làm chậm ($\text{MTC}$) âm ngay cả khi nồng độ Boron giảm xuống 0 ppm.

Khung lý thuyết của công trình tích hợp phương trình vận chuyển nơtron Boltzmann giải bằng phương pháp Monte Carlo năng lượng liên tục (Continuous-Energy Monte Carlo), lý thuyết suy giảm hạt nhân Bateman cải tiến và mô hình hình học thống kê ngẫu nhiên (Statistical Geometry Model - STGM). Phạm vi nghiên cứu bao quát bó nhiên liệu VVER-1000/V-320 gồm 331 vị trí (312 thanh nhiên liệu, 18 ống dẫn thanh điều khiển CPSAR, 1 ống đo trung tâm, assembly pitch 23,6 cm) trong dải độ sâu cháy từ 0 đến $40\text{ GWd/t}$. Tác động định lượng mang tính đột phá của luận án là xác lập thiết kế bó nhiên liệu mới vừa tăng độ dẫn nhiệt của thanh nhiên liệu chứa chất hấp thụ lên trên $35%$, vừa giảm hệ số đỉnh công suất và đảm bảo biên an toàn phản hồi nhiệt âm tuyệt đối.

Literature Review và Positioning

Các nghiên cứu tối ưu hóa chất hấp thụ cháy được trong lò phản ứng nước nhẹ (LWR) phân hóa thành hai trường phái chính:

  1. Trường phái tối ưu hóa chất hấp thụ đồng nhất và phân bố độ giàu: Hirano và cộng sự (1999) đã sử dụng phương pháp liệt kê đầy đủ (Exhaustive Enumeration) để tối ưu hóa phân bố độ giàu trong lò BWR nhằm hạ thấp hệ số đỉnh công suất. Lim và cộng sự (2001) ứng dụng kỹ thuật ma trận phi tuyến để làm phẳng phân bố công suất bó nhiên liệu. Haibach và cộng sự (2002) cùng Yilmaz và cộng sự (2006) khai thác thuật toán di truyền (Genetic Algorithm) trên lò PWR để tìm kiếm hàm lượng $\text{Gd}_2\text{O}_3$ tối ưu. Francois và cộng sự (2003) sử dụng thuật toán tìm kiếm Tabu để cực tiểu hóa độ lệch phản ứng và PPF.
  2. Trường phái thiết kế lớp phủ và chất hấp thụ lai (Hybrid/Heterogeneous Absorbers): Reda và cộng sự (2001) tăng số lượng thanh chất hấp thụ từ 12 lên 24 thanh trong lò PWR; Galahom (2016) khảo sát ảnh hưởng của các loại BA lên hệ số nhân hiệu dụng; Dandi và cộng sự (2014) kết hợp thanh WABA ($\text{Al}_2\text{O}_3/\text{B}_4\text{C}$) và Erbium ($\text{Er}_2\text{O}_3$) nhằm giảm độ phản ứng dư xuống $11% - 20%$; Choe và cộng sự (2016) đề xuất cấu hình lai $\text{UO}_2 - ^{157}\text{Gd}_2\text{O}_3$ phủ lớp $\text{ZrEr}_2$ bên ngoài viên nhiên liệu.

Mâu thuẫn học thuật sâu sắc nảy sinh giữa hai quan điểm:

  • Quan điểm truyền thống: Trộn đều oxit gadolinium vào ma trận $\text{UO}_2$ tạo ra môi trường đồng nhất isotropic, thuận lợi cho việc kiểm soát nơtron nhưng các nghiên cứu nhiệt lý (Bakker et al., 1995; Amaya et al., 2002; Massih et al., 2007) chứng minh hiện tượng suy thoái mạng tinh thể làm giảm mạnh độ dẫn nhiệt, gia tăng nhiệt độ tâm viên nhiên liệu (fuel centerline temperature).
  • Quan điểm vi hạt dị thể: Khẳng định việc gom chất hấp thụ thành các vi hạt phân tán sẽ bảo tồn mạng dẫn nhiệt của $\text{UO}_2$ nguyên chất, nhưng vấp phải trở ngại lớn về sự tự che chắn nơtron (spatial self-shielding effect) làm thay đổi phi tuyến tốc độ suy giảm đồng vị.

Luận án của Hoàng Thanh Phi Hùng định vị chính xác tại điểm giao thoa này. Bằng cách thiết lập mô hình tính toán hình học thống kê ngẫu nhiên (STGM) trên mã Monte Carlo MVP với thư viện JENDL-3.3, nghiên cứu đã vượt qua giới hạn của phương pháp tất định (như mã SRAC) vốn không thể mô phỏng chính xác vị trí ngẫu nhiên của các hạt micromet. So với các công trình quốc tế của Murata et al. (2002) (chỉ áp dụng vi hạt nhiên liệu phủ CFP cho lò nhiệt độ cao HTR) hay nghiên cứu của Choe et al. (2016) (trên ô mạng vuông PWR), công trình này tiên phong giải quyết trọn vẹn bài toán vi hạt $\text{Gd}_2\text{O}_3$ trong ô mạng tam giác lục giác đặc thù của lò VVER-1000.

Đó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 lý thuyết vận chuyển nơtron Boltzmann trong môi trường dị thể ngẫu nhiên bậc hai (Double Heterogeneity). Về mặt vật lý hạt nhân, tiết diện hấp thụ nơtrôn tại năng lượng nhiệt $0,0253\text{ eV}$ ($20^\circ\text{C}$) của $^{155}\text{Gd}$ đạt $60.737\text{ b}$ và $^{157}\text{Gd}$ đạt $252.912\text{ b}$ (với độ phổ biến tự nhiên lần lượt là $14,8%$ và $15,65%$). Quãng đường tự do trung bình của nơtron nhiệt trong $\text{Gd}_2\text{O}_3$ nguyên chất chỉ khoảng $8,3\ \mu\text{m}$. Khi bán kính vi hạt $R_0 \gg 8,3\ \mu\text{m}$, vi hạt đóng vai trò như một "vật hấp thụ đen" (black absorber).

Tiết diện hấp thụ vĩ mô hiệu dụng biến thiên phi tuyến theo thời gian cháy $t$ được luận án mô hình hóa qua phương trình giải tích:

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

Trong đó $n_P$ là mật độ số vi hạt, $R_0$ là bán kính ban đầu, $N$ là mật độ nguyên tử, và $\Phi(t)$ là thông lượng nơtron.

   Nơtron nhiệt tới

Công trình chứng minh bước chuyển dịch mô hình (paradigm shift): tốc độ cháy của chất hấp thụ không còn bị chi phối đơn thuần bởi mật độ nguyên tử trung bình, mà được điều tiết bởi tỷ số hình học giữa bán kính hạt $R_0$ và quãng đường tự do trung bình. Điều này thách thức giả định đồng nhất hóa truyền thống trong công thức 4 yếu tố:

$$k_\infty = \varepsilon \cdot p_r \cdot \eta \cdot f$$

Hệ số sử dụng nơtron nhiệt $f$ và xác suất tránh hấp thụ cộng hưởng $p_r$ được điều chỉnh chính xác thông qua phân bố vi mô không gian.

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

Khung phân tích của luận án là sự hợp nhất của ba trụ cột lý thuyết:

  1. Lý thuyết vận chuyển nơtron liên tục: Dựa trên phương trình tích phân mật độ va chạm Boltzmann ba chiều.
  2. Động học biến đổi đồng vị Bateman: Được giải chính xác từng bước cháy thông qua mô-đun MVP_BURN.
  3. Lý thuyết phân bố thống kê lân cận gần nhất (Nearest Neighbor Distribution - NND): Mô phỏng cấu trúc nén ngẫu nhiên của các quả cầu rắn (hard spheres packing) thông qua hàm mật độ xác suất:

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

Với $f_p$ là tỷ lệ thể tích chiếm chỗ (packing fraction) và $r$ là khoảng cách chuẩn hóa. Khung phân tích này thiết lập điều kiện biên phản xạ tuyệt đối (specular reflection) 6 mặt trên lăng trụ lục giác, loại bỏ hoàn toàn các sai số làm tròn hình học.

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

Thiết kế nghiên cứu

Nghiên cứu đứng trên lập trường nhận thức luận thực chứng định lượng (Quantitative Positivism), sử dụng phương pháp mô phỏng ngẫu nhiên Monte Carlo năng lượng liên tục có độ chính xác cao nhất trong vật lý lò phản ứng. Thiết kế nghiên cứu đa tầng (Multi-level design) được cấu trúc chặt chẽ:

  • Tầng vi mô (Micro-scale): Mô phỏng hàng triệu hạt $\text{Gd}_2\text{O}_3$ ngẫu nhiên bên trong viên nhiên liệu $\text{UO}_2$ đường kính $7,56\text{ mm}$ (lỗ rỗng trung tâm $2,3\text{ mm}$).
  • Tầng thanh nhiên liệu (Pin-cell scale): Mô hình hóa vỏ bọc hợp kim $\text{Zr}-1%\text{Nb}$ (đường kính ngoài $9,1\text{ mm}$, trong $7,73\text{ mm}$), khe hở khí heli và nước làm chậm.
  • Tầng bó nhiên liệu (Assembly scale): Toàn bộ cấu trúc lục giác 331 vị trí với bước mạng thanh nhiên liệu (pitch) $12,75\text{ mm}$ và bước mạng bó nhiên liệu $236\text{ mm}$.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu số thông qua chuỗi công cụ chuyên dụng:

  • Hệ mã tính toán: Sử dụng chương trình Monte Carlo năng lượng liên tục MVP phiên bản II kết hợp mô-đun tính cháy 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 số liệu hạt nhân: JENDL-3.3 (Japanese Evaluated Nuclear Data Library) được xử lý thông qua hệ thống mã LICEM, cung cấp tiết diện phản ứng liên tục ở nhiệt độ vận hành chính xác của lò phản ứng ($T_{\text{fuel}} = 1027\text{ K}$, $T_{\text{mod}} = 575\text{ K}$).
  • Xử lý hiệu ứng nhiệt: Áp dụng mô hình khí tự do cho nơtron nhiệt và xử lý tán xạ nhiệt $S(\alpha, \beta)$ trong chất làm chậm $\text{H}_2\text{O}$; sử dụng phương pháp bảng xác suất (probability table method) trong vùng cộng hưởng chưa giải quyết.
  • Triangulation và kiểm định chéo: Kết quả tính toán của MVP được kiểm chứng độc lập (benchmark) với chương trình tất định đa nhóm năng lượng SRAC và số liệu chuẩn đối chuẩn quốc tế (Benchmark - BM) cho bó nhiên liệu VVER-1000 chứa $\text{Gd}_2\text{O}_3$ đồng nhất trong dải cháy $0 - 40\text{ GWd/t}$.

Data và phân tích

  • Quy mô thống kê Monte Carlo: Mỗi điểm cháy tính toán với 2.000.000 lịch sử nơtron (100 batches, mỗi batch 20.000 hạt, loại bỏ 20 batches đầu để nguồn nơtron hội tụ ổn định), đảm bảo sai số thống kê $1\sigma$ đối với $k_\infty$ luôn nhỏ hơn $0,0005$ ($\pm 50\text{ pcm}$).
  • Mô hình chuỗi cháy: Chuỗi cháy tiêu chuẩn u4cm6fp50bp16T tích hợp 4 actinides ($\text{UO}_2$ đến $\text{Cm}$), 50 sản phẩm phân hạch chính (FPs) và 16 đồng vị hấp thụ nơtron cháy được.
  • Môi trường tính toán: Thực thi song song hóa bộ nhớ phân tán thông qua giao thức MPI (Message Passing Interface) và máy tính ảo PVM (Parallel Virtual Machine).

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

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

  1. Xác định đường kính vi hạt tối ưu cho cấu hình 12 thanh $\text{UO}_2 - \text{Gd}_2\text{O}_3$: Khảo sát dải đường kính hạt từ 40 đến $100\ \mu\text{m}$ (với hàm lượng $5,0\text{ wt}%$) cho thấy đường kính $d = 60\ \mu\text{m}$ tái tạo hoàn hảo đường cong độ phản ứng $k_\infty$ của bó nhiên liệu tham chiếu (trộn đều đồng nhất) trong suốt giai đoạn $0 - 10\text{ GWd/t}$. Sai lệch $k_\infty$ giữa cấu hình vi hạt $60\ \mu\text{m}$ và cấu hình đồng nhất tại $0\text{ GWd/t}$ chỉ là $0,0012$ ($0,12%\ \Delta k$).
  2. Quy luật suy giảm đồng vị $^{155}\text{Gd}$ và $^{157}\text{Gd}$: Ở đường kính vi hạt $60\ \mu\text{m}$, mật độ nguyên tử của hai đồng vị bắt nơtron mạnh nhất giảm nhanh và cháy kiệt gần như hoàn toàn tại bước cháy $10\text{ GWd/t}$, không để lại lượng độc nơtron dư thừa (residual poison penalty) tại giai đoạn cuối chu trình (EOC), đảm bảo $k_\infty$ ở vùng $10 - 40\text{ GWd/t}$ bám sát bó nhiên liệu không chứa chất hấp thụ.
  3. Phát hiện giới hạn kích thước hạt: Khi đường kính vi hạt vượt quá $100\ \mu\text{m}$ trong cấu hình 12 thanh, hiệu ứng tự che chắn quá mạnh khiến tốc độ cháy bị kìm hãm quá mức, kéo dài thời gian tồn tại của gadolinium qua ngưỡng $10\text{ GWd/t}$, làm giảm độ phản ứng hữu ích ở cuối chu trình.
  4. Đột phá thiết kế bó nhiên liệu mới với 18 thanh $\text{UO}_2 - \text{Gd}_2\text{O}_3$: Bằng việc phân tán chất hấp thụ ra 18 thanh (thay vì 12 thanh) và tối ưu hóa đường kính vi hạt lên $d = 300\ \mu\text{m}$ với tỷ lệ thể tích $3,33%$, công trình đã hạ thấp đáng kể hệ số đỉnh công suất cục bộ (PPF). Tại thời điểm $0\text{ GWd/t}$, hệ số PPF giảm từ $1,342$ (ở bó 12 thanh) xuống dưới $1,285$ (ở bó 18 thanh), loại bỏ hoàn toàn các điểm tập trung nhiệt gây nguy hiểm.
  5. Vận hành an toàn ở nồng độ Boron thấp và 0 ppm Boron: Thiết kế bó nhiên liệu 18 thanh vi hạt cho phép giảm nồng độ axit boric hòa tan trong nước làm mát xuống $50%$ và $0\text{ ppm}$ mà vẫn duy trì hệ số phản hồi nhiệt độ chất làm chậm ($\text{MTC}$) mang giá trị âm sâu ($\text{MTC} < 0$) trong toàn bộ chu trình cháy. Điều này triệt tiêu hoàn toàn rủi ro xuất hiện $\text{MTC}$ dương – một mối đe dọa an toàn nghiêm trọng đối với các lò PWR/VVER truyền thống khi vận hành với nồng độ boron cao tại đầu chu trình (BOC).
Thông số kỹ thuật Thiết kế truyền thống (Tham chiếu) Cải tiến 1: Bó 12 thanh vi hạt Cải tiến 2: Bó 18 thanh vi hạt (Đột phá)
Dạng phân bố $\text{Gd}_2\text{O}_3$ Trộn đều đồng nhất Vi hạt ngẫu nhiên (STGM) Vi hạt ngẫu nhiên (STGM)
Đường kính vi hạt ($d$) Không hạt (Hỗn hợp bột) $60\ \mu\text{m}$ $300\ \mu\text{m}$
Hàm lượng $\text{Gd}_2\text{O}_3$ $5,0\text{ wt}%$ $5,0\text{ wt}%$ $3,33\text{ vol}%$
Độ dẫn nhiệt viên nhiên liệu $3,8 - 2,6\text{ W/mK}$ $5,8 - 2,7\text{ W/mK}$ $5,8 - 2,7\text{ W/mK}$
Hệ số đỉnh công suất PPF (BOC) $1,348$ $1,342$ $1,285$ (Giảm mạnh)
Khả năng giảm nồng độ Boron Bắt buộc Boron cao ($>1000\text{ ppm}$) Hạn chế Hoạt động tối ưu ở 50% & 0 ppm Boron
Hệ số MTC tại BOC Nguy cơ MTC kém âm / dương Âm Âm sâu ($\text{MTC} \ll 0$)

Implications đa chiều

  • Tiến bộ lý thuyết: Cung cấp dữ liệu thực chứng chuẩn mực về sự biến thiên tiết diện vi mô trong môi trường hạt ngẫu nhiên, hoàn thiện lý thuyết tự che chắn nơtron trong ô mạng tam giác lục giác.
  • Đổi mới phương pháp: Chuẩn hóa quy trình mô phỏng liên hợp STGM - MVP_BURN - JENDL-3.3 cho các bài toán phân tích dị thể cấp độ cao trong vật lý lò hạt nhân.
  • Ứng dụng thực tiễn: Nâng cao giới hạn biên an toàn nhiệt (Thermal Margin), cho phép các nhà máy VVER-1000 nâng công suất phát điện (Power Uprating) lên $104% - 107%$ định mức mà không cần thay đổi kích thước hình học thùng lò.
  • Chính sách & An toàn hạt nhân: Cung cấp cơ sở khoa học để các cơ quan pháp quy hạt nhân phê duyệt các thiết kế nhiên liệu kéo dài chu trình vận hành (18 - 24 tháng) với độ an toàn nội tại nâng cao.

Limitations và Future Research

  1. Giới hạn mô hình nhiệt thủy lực: Các tính toán Monte Carlo được thực hiện ở trạng thái nhiệt độ nhiên liệu và chất làm chậm đẳng nhiệt/phân tầng cố định, chưa ghép cặp động học tức thời (online coupled TH) với các chương trình nhiệt thủy lực đa chiều như RELAP5 hoặc COBRA-IV.
  2. Giả định hình học vi hạt đơn phân tán (Monodisperse): Mô hình STGM giả định tất cả vi hạt $\text{Gd}_2\text{O}_3$ có cùng một kích thước đường kính hình học duy nhất ($60\ \mu\text{m}$ hoặc $300\ \mu\text{m}$), trong khi thực tế chế tạo bột thiêu kết luôn tồn tại phân bố kích thước hạt ngẫu nhiên dạng Gauss (Polydisperse).
  3. Phạm vi mô phỏng cấp độ bó nhiên liệu: Nghiên cứu tập trung vào bó nhiên liệu 2D/3D với điều kiện biên phản xạ; chưa mở rộng mô phỏng toàn vùng hoạt (full-core 3D) 163 bó với các sơ đồ đảo đảo nhiên liệu phức tạp (fuel reloading schemes).

Chương trình nghiên cứu tiếp theo (Future Research Agenda):

  • Phát triển hệ mã ghép cặp nơtron - nhiệt thủy lực độ chính xác cao (Continuous-Energy Monte Carlo MVP / Thermal-Hydraulics Code Coupling).
  • Khảo sát ảnh hưởng của phân bố kích thước vi hạt đa phân tán (Polydisperse size distribution) đến hiệu ứng tự che chắn.
  • Mở rộng thiết kế vi hạt $\text{Gd}_2\text{O}_3$ sang cấu hình bó nhiên liệu của lò phản ứng VVER-1200 (thế hệ III+) với độ cháy mục tiêu $70\text{ MWd/kgU}$.
  • Đánh giá thực nghiệm cơ tính và độ bền bức xạ của viên nén nhiên liệu $\text{UO}_2 - \text{Gd}_2\text{O}_3$ vi hạt sau chiếu xạ (Post-Irradiation Examination - PIE).

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

  • Tác động học thuật: Mở ra hướng nghiên cứu chuyên sâu về vi cấu trúc hạt hấp thụ trong lò phản ứng năng lượng; dự kiến tạo nguồn trích dẫn quan trọng cho các nghiên cứu mô phỏng nơtron dị thể trên các tạp chí quốc tế hàng đầu (như Annals of Nuclear Energy, Nuclear Engineering and Design, Progress in Nuclear Energy).
  • Chuyển đổi công nghiệp: Cung cấp giải pháp kỹ thuật trực tiếp cho các tập đoàn chế tạo nhiên liệu hạt nhân (như TVEL/Rosatom) cải tiến quy trình nén ép viên nhiên liệu, nâng cao độ dẫn nhiệt và kéo dài tuổi thọ thanh nhiên liệu.
  • Chính sách năng lượng: Hỗ trợ luận cứ kỹ thuật cho các quốc gia đang vận hành hoặc xây dựng kế hoạch phát triển lò VVER trong việc thẩm định an toàn thiết kế nhiên liệu mới nồng độ Boron thấp.
  • Lợi ích kinh tế - xã hội: Giảm chi phí chu trình nhiên liệu thông qua việc nâng cao hiệu suất cháy uranium, hạn chế lượng chất thải phóng xạ sinh ra trên mỗi megawatt-giờ điện, đồng thời loại trừ nguy cơ sự cố mất nước làm mát (LOCA) do thanh nhiên liệu có biên an toàn nhiệt cao hơn.

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

  • Nghiên cứu sinh & Nhà khoa học hạt nhân: Thừa hưởng phương pháp luận mô hình hóa hình học thống kê ngẫu nhiên (STGM) trên mã MVP để áp dụng cho các mô hình nhiên liệu tiên tiến khác (như ATF - Accident Tolerant Fuel, nhiên liệu hạt TRISO).
  • Kỹ sư thiết kế vùng hoạt & An toàn lò phản ứng: Nhận được bộ thông số thiết kế tối ưu ($d = 60\ \mu\text{m}$ cho bó 12 thanh; $d = 300\ \mu\text{m}$ cho bó 18 thanh) để ứng dụng trực tiếp vào tối ưu hóa nạp đảo nhiên liệu lò VVER-1000.
  • Nhà sản xuất nhiên liệu hạt nhân: Nắm bắt quy trình công nghệ chế tạo viên nén vi hạt $\text{Gd}_2\text{O}_3$ giúp nâng cao độ dẫn nhiệt thực tế lên $5,8\text{ W/mK}$.
  • Cơ quan quản lý an toàn bức xạ và hạt nhân: Có cơ sở dữ liệu đối chuẩn độc lập để đánh giá các báo cáo phân tích an toàn (SAR) khi các nhà máy điện hạt nhân xin cấp phép nâng công suất hoặc chuyển đổi loại nhiên liệu.

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ở rộng lý thuyết vận chuyển nơtron Boltzmann trong môi trường dị thể ngẫu nhiên bậc hai (Double Heterogeneity) trên cấu trúc ô mạng lục giác VVER. Luận án đã giải quyết bài toán tự che chắn không gian nơtron nhiệt của các vi hạt hấp thụ đen $\text{Gd}_2\text{O}_3$ bằng mô hình giải tích tiết diện hiệu dụng $\Sigma_a(t)$ phụ thuộc bán kính hạt biến thiên, chứng minh rằng sự kết hợp giữa mô hình hình học thống kê (STGM) và mã Monte Carlo liên tục MVP có thể lượng hóa chính xác hiệu ứng tự che chắn mà không cần qua các bước đồng nhất hóa trung gian gây sai số.

2. Đổi mới phương pháp luận nghiên cứu thể hiện như thế nào khi so sánh với các nghiên cứu trước đây?

So với các nghiên cứu của Hirano et al. (1999) hay Francois et al. (2003) vốn sử dụng các chương trình tất định đa nhóm với xấp xỉ đồng nhất hóa mạng nơtron, luận án đã tiên phong ứng dụng mã Monte Carlo năng lượng liên tục MVP II kết hợp mô hình phân bố lân cận gần nhất (NND) để mô tả ngẫu nhiên cấu trúc hạt micromet. So với mã SRAC truyền thống, phương pháp của luận án phản ánh trung thực tương tác nơtron trong từng vi hạt mà không gặp hiện tượng đánh giá thấp (underestimation) hoặc đánh giá cao (overestimation) tốc độ cháy đồng vị $^{155}\text{Gd}/^{157}\text{Gd}$.

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 mức rất lớn ($d = 300\ \mu\text{m}$) khi phân tán vào 18 thanh nhiên liệu không những không làm suy giảm khả năng kiểm soát độ phản ứng, mà còn tạo ra sự phân bố công suất đồng đều vượt trội (PPF hạ xuống $1,285$) và cho phép lò vận hành ổn định tại nồng độ Boron $0\text{ ppm}$ với hệ số $\text{MTC}$ hoàn toàn âm. Thông thường, kích thước hạt lớn bị xem là gây bất lợi do hiện tượng che chắn quá mức, nhưng trong cấu hình 18 thanh, nó đã biến thành lợi thế điều hòa độ phản ứng phẳng và mượt mà.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) không?

Luận án cung cấp chi tiết toàn bộ thông số hình học, thành phần vật liệu, mật độ nguyên tử, cấu hình bó nhiên liệu VVER-1000/V-320 (331 vị trí, bán kính viên, vỏ bọc $\text{Zr}-1%\text{Nb}$), bộ thư viện tiết diện hạt nhân JENDL-3.3 qua mã LICEM, chuỗi cháy u4cm6fp50bp16T, cùng các điều kiện biên phản xạ tuyệt đối. Mọi nhà nghiên cứu có mã MVP/GMVP II đều có thể tái lập chính xác $100%$ các kết quả tính toán.

5. Khung chương trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Chương trình 10 năm tập trung vào 4 giai đoạn chiến lược:

  1. (Năm 1 - 2): Xây dựng nền tảng ghép cặp trực tiếp mã nơtron MVP với mã nhiệt thủy lực 3D Subchannel (MVP-COBRA).
  2. (Năm 3 - 5): Mô phỏng toàn vùng hoạt 3D (Full-core) cho lò VVER-1000 và VVER-1200 với các chu trình nạp đảo nhiên liệu đa dạng.
  3. (Năm 6 - 8): Thử nghiệm công nghệ chế tạo viên nén nhiên liệu $\text{UO}_2 - \text{Gd}_2\text{O}_3$ vi hạt đa phân tán trong điều kiện thiêu kết công nghiệp.
  4. (Năm 9 - 10): Tiến hành chiếu xạ thử nghiệm trong lò nghiên cứu và phân tích kiểm tra sau chiếu xạ (PIE) để chuẩn hóa đưa vào ứng dụng thương mại.

Kết luận

  1. Giải quyết triệt để mâu thuẫn nhiệt - nơtron: Công trình đã chứng minh tính khả thi tuyệt đối của việc thay thế bột $\text{Gd}_2\text{O}_3$ trộn đều bằng vi hạt $\text{Gd}_2\text{O}_3$ ngẫu nhiên, giúp nâng độ dẫn nhiệt của viên nhiên liệu lên $5,8 - 2,7\text{ W/mK}$ (tăng trên $35%$) trong dải nhiệt độ $300 - 1273\text{ K}$, cải thiện vượt bậc an toàn nhiệt lò phản ứng.
  2. Xác lập chuẩn thiết kế bó 12 thanh vi hạt: Tìm ra kích thước vi hạt tối ưu $d = 60\ \mu\text{m}$ ($5,0\text{ wt}%$) cho bó 12 thanh $\text{UO}_2 - \text{Gd}_2\text{O}3$, tái lập chính xác đường cong độ phản ứng $k\infty$ của thiết kế truyền thống trong dải cháy $0 - 10\text{ GWd/t}$ và cháy kiệt hoàn toàn ở cuối chu trình.
  3. Đột phá thiết kế bó 18 thanh vi hạt: Phát triển thành công cấu hình bó nhiên liệu mới 18 thanh $\text{UO}_2 - \text{Gd}_2\text{O}_3$ ($d = 300\ \mu\text{m}$, $3,33\text{ vol}%$) giúp hạ thấp hệ số đỉnh công suất PPF xuống mức $1,285$, triệt tiêu các điểm quá nhiệt cục bộ trong vùng hoạt.
  4. Hiện thực hóa chế độ vận hành Boron thấp: Thiết kế cải tiến 18 thanh mở ra khả năng vận hành lò VVER-1000 ở chế độ giảm $50%$ hoặc $0\text{ ppm}$ Boron hòa tan, đảm bảo hệ số phản hồi nhiệt độ chất làm chậm ($\text{MTC}$) luôn âm sâu, nâng cao độ an toàn nội tại của nhà máy điện hạt nhân.
  5. Đóng góp phương pháp luận chuẩn xác cao: Xác lập quy trình mô phỏng nơtron vi cấu trúc dị thể tiên tiến bằng mã Monte Carlo MVP II trên thư viện JENDL-3.3, đặt nền móng vững chắc cho các nghiên cứu thiết kế nhiên liệu hạt nhân thế hệ mới tại Việt Nam và quốc tế.