Tổng quan về luận án

Trong bối cảnh an ninh năng lượng toàn cầu đối mặt với sự cạn kiệt của nhiên liệu hóa thạch và nguy cơ cạn kiệt quặng urani ($^{235}\text{U}$ chỉ chiếm khoảng $0{,}71%$ urani tự nhiên), điện hạt nhân truyền thống đang phải giải quyết đồng thời ba thách thức cốt lõi: nguy cơ mất an toàn tới hạn (như các thảm họa Chernobyl năm 1986, Fukushima năm 2011), chi phí kinh tế khổng lồ và gánh nặng xử lý chất thải phóng xạ có đời sống dài (minor actinides). Luận án tiến sĩ vật lý của tác giả Trần Minh Tiến (2022) với đề tài "Nghiên cứu khả năng sử dụng thori làm nhiên liệu cho lò phản ứng hạt nhân điều khiển bằng máy gia tốc" đã tiên phong giải quyết các bế tắc này thông qua việc nghiên cứu hệ thống lò phản ứng hạt nhân dưới tới hạn điều khiển bằng máy gia tốc (Accelerator Driven Subcritical Reactor - ADSR), tích hợp cơ chế đột phá: sử dụng trực tiếp chì lỏng đồng thời làm bia tương tác sinh neutron $(p,n)$ và chất tải nhiệt làm mát lõi lò, kết hợp chu trình nhiên liệu Thori ($^{232}\text{Th}$).

Khoảng trống nghiên cứu (research gap) được xác định rõ rệt: hầu hết các nghiên cứu ADSR quốc tế trước đây (tiêu biểu như C. Bowman, 1992; C. Rubbia, 1995; S. Meigo, 1999) chủ yếu tập trung vào bia tương tác chì hoặc chì-bismuth ở trạng thái rắn. Bia rắn chịu ứng suất nhiệt cực lớn và hư hỏng bức xạ nghiêm trọng, đòi hỏi phải định kỳ dừng vận hành lò để thay thế bia. Một số thiết kế bia lỏng hợp kim chì-bismuth (LBE) lại sản sinh đồng vị độc hại dễ bay hơi $^{210}\text{Po}$ qua phản ứng kích hoạt neutron trên $^{209}\text{Bi}$, trong khi các nghiên cứu sử dụng chì lỏng nguyên chất vừa làm bia vừa làm chất tải nhiệt toàn phần vẫn còn thiếu vắng các tính toán vi phân vi mô và thông lượng neutron toàn diện trong cấu hình nhiên liệu thori.

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

  1. Câu hỏi 1 (RQ1): Tương tác của chùm proton năng lượng cao ($0{,}25 - 3{,}0\text{ GeV}$) trên môi trường bia chì lỏng liên tục có đặc trưng vi phân tiết diện sinh neutron và phân bố góc/năng lượng như thế nào so với bia chì rắn? Giả thuyết 1 (H1): Mặc dù mật độ khối lượng của chì lỏng ($10{,}66\text{ g/cm}^3$) thấp hơn chì rắn ($11{,}7\text{ g/cm}^3$), độ dài tương tác hiệu dụng mở rộng trên toàn đường kính lõi sẽ bù đắp hoàn toàn sự suy giảm mật độ, bảo toàn hiệu suất sinh neutron $Y_{n/p}$.
  2. Câu hỏi 2 (RQ2): Phổ phóng xạ phân rã hạt nhân và thông lượng neutron ($NF$) phân bố theo không gian (chiều cao $z$, bán kính $r$) và năng lượng bên trong lõi ADSR làm mát bằng chì lỏng biến thiên ra sao? Giả thuyết 2 (H2): Phổ neutron duy trì tỷ phần neutron nhanh và trên nhiệt đủ cao để kích hoạt quá trình hấp thụ neutron trên $^{232}\text{Th}$ chuyển hóa thành $^{233}\text{U}$.
  3. Câu hỏi 3 (RQ3): Khả năng duy trì hệ số nhân neutron hiệu dụng ($k_{eff}$) ở trạng thái dưới tới hạn an toàn ($0{,}90 \le k_{eff} < 1{,}0$) của các tổ hợp nhiên liệu thori pha tạp ($^{232}\text{Th}-^{233}\text{U}\text{O}_2$, $^{232}\text{Th}-^{235}\text{U}\text{O}_2$, $^{232}\text{Th}-^{238}\text{U}\text{O}_2$) biến thiên theo tỷ lệ phần trăm thori ($0% - 100%$) như thế nào? Giả thuyết 3 (H3): Nhiên liệu $^{232}\text{Th}-^{233}\text{U}\text{O}_2$ duy trì độ phản ứng vượt trội so với $^{232}\text{Th}-^{235}\text{U}\text{O}_2$, cho phép nâng tỷ lệ thori lên mức tối đa mà không làm mất tính ổn định dưới tới hạn.

Khung lý thuyết của công trình tích hợp lý thuyết tương tác spallation nội hạt nhân (Intranuclear Cascade Theory), lý thuyết phân rã thống kê (Statistical Decay Model - SDM), động lực học phân tử lượng tử (Quantum Molecular Dynamics - QMD) và lý thuyết vận chuyển neutron Boltzmann. Đóng góp mang tính đột phá của nghiên cứu là việc lượng hóa toàn diện bộ thông số neutron học vi mô trên bia chì lỏng, mở rộng phạm vi trữ lượng nhiên liệu hạt nhân khi thori dồi dào gấp 4 lần urani trên vỏ trái đất (ước tính khoảng $1{,}2 \times 10^{14}\text{ tấn}$ Thori toàn cầu với tỷ lệ $6\text{ ppm}$ trong vỏ thạch quyển). Nghiên cứu được mô phỏng chi tiết trên cấu trúc hình học chuẩn của lò nghiên cứu TRIGA Mark II với dòng proton $25\text{ mA}$, bán kính chùm $4\text{ cm}$, giải quyết triệt để rào cản phải định kỳ tắt lò thay bia trong công nghệ ADSR.

Literature Review và Positioning

Nghiên cứu về hệ thống điều khiển bằng máy gia tốc khởi nguồn từ các quan sát biến đổi hạt nhân của Ernest Rutherford (1919) và phát minh cyclotron của Ernest O. Lawrence (1940), kết hợp đề xuất dùng máy gia tốc làm nguồn neutron của N. N. Semenov và khám phá $^{238}\text{Pu}$ của Glenn Seaborg. Đến thập niên 1990, C. Bowman (1992) đề xuất hệ thống ADS nhằm biến đổi chất thải phóng xạ actinide, và Carlo Rubbia (1995) công bố khái niệm Energy Amplifier sử dụng chùm proton bắn phá bia chì kết hợp chu trình $^{232}\text{Th}-^{233}\text{U}$ làm mát bằng chì lỏng ở dải công suất nhiệt $600\text{ MWth}$.

Trong dòng nghiên cứu vật lý bia tương tác spallation, các tranh luận học thuật quốc tế tập trung vào hai luồng quan điểm đối lập:

  • Luồng quan điểm thứ nhất (Bia rắn cổ điển): S. Meigo và cộng sự (1999) sử dụng MCNP4A phân tích bia chì dày với chùm proton $0{,}5\text{ GeV}$, chỉ ra rằng thông lượng neutron nhiệt vượt trội hơn neutron nhanh ở mọi góc đo; J. Goldenbaum (1999) chứng minh giới hạn bề dày bia tối ưu (chì là $55\text{ cm}$, vonfram là $30\text{ cm}$); G. Bauer (2001) cùng H. Nifenecker (2001) và A. Krasa (2009) tính toán phổ neutron ở góc $30^\circ, 90^\circ, 150^\circ$ với proton $2\text{ GeV}$, khẳng định bia rắn cho mật độ nguyên tử cao nhưng dễ bị hỏng hóc do tải nhiệt tập trung và tích tụ khí heli/hydro.
  • Luồng quan điểm thứ hai (Bia kim loại lỏng eutectic): Điển hình là thiết kế lò lai đa mục đích MYRRHA tại SCK.CEN (Bỉ) sử dụng hợp kim chì-bismuth lỏng (LBE) với proton $600\text{ MeV} - 4\text{ mA}$, công suất $65 - 100\text{ MWth}$ tạo thông lượng neutron nhanh $10^{15}\text{ n}\cdot\text{cm}^{-2}\text{s}^{-1}$; dự án C-ADS của Viện Hàn lâm Khoa học Trung Quốc (CAS: IHEP, IMP, IPP, USTC) hướng tới công suất $100\text{ MW}$ (2022) và $1\text{ GW}$ (2032); tổ hợp dưới tới hạn KIPT (Kharkov, Ukraine) hợp tác cùng Phòng thí nghiệm Quốc gia Argonne (ANL, Hoa Kỳ) dùng máy gia tốc thẳng và nhiên liệu urani oxit làm giàu thấp. Điểm yếu của hợp kim chì-bismuth là sự hình thành đồng vị phát bức xạ alpha độc tính cao $^{210}\text{Po}$ qua phản ứng $^{209}\text{Bi}(n,\gamma)^{210}\text{Bi} \xrightarrow{\beta^-} {}^{210}\text{Po}$.

Ở Việt Nam, các nghiên cứu tiên phong của nhóm PGS. Nguyễn Mộng Giao, PGS. Châu Văn Tạo và TS. Nguyễn Thị Ái Thu (2002–2015) đã phát triển mô hình màn chắn trên bia rắn (U, Pb, W, Au) ở dải năng lượng $0{,}5 - 1{,}5\text{ GeV}$. Tiếp đó, các công trình của TS. Vũ Thanh Mai cùng Cheol Ho Pyeon và Masao Yamanaka (2014–2018) tại tổ hợp tới hạn KUCA (Đại học Kyoto, Nhật Bản) khảo sát động học neutron tức thời và trễ với chùm proton $100\text{ MeV}$ trên bia Pb-Bi kết hợp nhiên liệu $^{235}\text{U}-^{232}\text{Th}$.

Luận án của Trần Minh Tiến đã định vị chính xác khoảng trống học thuật: chuyển đổi hoàn toàn từ bia rắn và hợp kim Pb-Bi sang chì lỏng nguyên chất hai chức năng (vừa làm bia spallation, vừa làm chất làm mát) trong cấu hình lõi TRIGA Mark II nạp nhiên liệu Thori oxit. So với mô hình ADTR ($85%\text{ Th} + 15%\text{ Pu}$) của Victoria B. Ashley (2014) hay mô hình ADS muối nóng chảy $60\text{NaCl}-34\text{ThCl}_4-6{}^{233}\text{UCl}4$ của Nhật Bản ($800\text{ MW}$, $Y{n/p} = 27\text{ n/p}$ tại $1{,}5\text{ GeV}$), nghiên cứu này cung cấp bộ dữ liệu vi phân bậc hai chính xác cao từ thư viện JENDL-HE-2007, chứng minh tính khả thi vật lý vượt trội mà không cần dùng đến muối nóng chảy ăn mòn hay hợp kim chứa Bismuth nguy hiểm.

Đó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 phát triển năng lượng phân hạch dưới tới hạn của Carlo Rubbia và lý thuyết chu trình nhiên liệu Thori ($^{232}\text{Th}-^{233}\text{U}$ breeding cycle):

  1. Cơ chế tái tạo neutron liên tục trên môi trường lỏng: Mở rộng lý thuyết tương tác spallation trên môi trường đối lưu lỏng liên tục, chứng minh rằng tương tác $(p,n)$ không bị giới hạn bởi hiện tượng suy thoái cấu trúc mạng tinh thể (lattice damage) như ở vật liệu bia rắn.
  2. Mô hình hóa chuỗi phân rã phóng xạ trong môi trường tải nhiệt kim loại nặng: Định lượng hóa chính xác phổ năng lượng các bức xạ $\alpha, \beta, \gamma, \bar{\nu}_e$ và động năng của các hạt nhân con trong môi trường chì lỏng, cung cấp bằng chứng lý thuyết về khả năng tự che chắn bức xạ và ổn định phổ của chì.
  3. Quy luật diễn tiến độ phản ứng ($k_{eff}$) theo tỷ lệ pha tạp Thori: Làm sáng tỏ cơ chế vật lý hạt nhân giải thích tại sao hỗn hợp $^{232}\text{Th}-^{233}\text{U}\text{O}2$ duy trì giá trị $k{eff}$ ổn định vượt trội hơn hẳn so với $^{232}\text{Th}-^{235}\text{U}\text{O}_2$ và $^{232}\text{Th}-^{238}\text{U}\text{O}_2$ khi nâng tỷ lệ thori từ $0%$ lên đến $100%$, củng cố nền tảng lý thuyết cho các lò phản ứng thế hệ IV tăng cường an toàn nội tại (inherent safety).

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

Khung phân tích của luận án kết hợp 3 trụ cột lý thuyết hạt nhân: (1) Lý thuyết phản ứng spallation năng lượng cao (ENDF/JENDL-HE); (2) Lý thuyết vận chuyển bức xạ ngẫu nhiên Monte Carlo; (3) Lý thuyết động học lò phản ứng dưới tới hạn ($k_{eff} < 1$).

Toàn bộ quá trình sinh neutron từ phản ứng $(p,n)$ trên bia chì lỏng được mô hình hóa toán học thông qua phương trình vi phân bậc hai của tiết diện sinh neutron theo năng lượng và theo góc khối (neutron production double-differential cross section): $$\frac{d^2\sigma(\mu, E_p, E_n)}{dE d\Omega} = \sigma(E_p) \cdot \frac{\gamma(E_p)}{4\pi} \cdot f(\mu, E_p, E_n)$$

Tiết diện vi phân theo góc khối được rời rạc hóa trên lưới 19 điểm giá trị góc ($\mu = [\mu_1, \mu_2, \dots, \mu_{19}]$ tương ứng 18 khoảng phân chia từ $0^\circ$ đến $180^\circ$) và 32 mức năng lượng neutron rời rạc ($E_n = [E_1, E_2, \dots, E_{32}]$): $$\frac{d\sigma}{d\Omega} = \sigma(E_p) \sum_{i=1}^{19} \frac{f_{i+1}(\mu, E_p, E_n) + f_i(\mu, E_p, E_n)}{2}$$

Cường độ dòng neutron $N_n$ sinh ra từ chùm proton $N_p$ được xác định qua tổng tích chập tiết diện: $$N_n = N_p \sum_{j=1}^{32} N \cdot d_j \frac{d\sigma(\mu, E_p, E_n)}{dE}$$

Điều kiện biên xác định: Phổ proton tới $E_p \in [0{,}25; 3{,}0]\text{ GeV}$, chùm hạt đơn năng bán kính $r = 4\text{ cm}$, dòng điện $I_p = 25\text{ mA}$; môi trường chì lỏng ở nhiệt độ vận hành trên điểm nóng chảy $327^\circ\text{C}$ với mật độ chuẩn $\rho = 10{,}66\text{ g/cm}^3$.

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

Thiết kế nghiên cứu

Luận án tuân thủ chặt chẽ thế giới quan bản thể luận thực chứng (positivism) và nhận thức luận định lượng (quantitative epistemology) với độ tin cậy mô phỏng đa tầng. Thiết kế nghiên cứu tích hợp phương pháp Monte Carlo stochastics thông qua hệ mã vận chuyển hạt tiên tiến MCNPX (phiên bản mở rộng đa hạt của Los Alamos National Laboratory), phối hợp với MCNP5 và GEANT4 để chéo hóa kết quả.

Cấu trúc hình học mô phỏng được xây dựng chuẩn xác dựa trên lõi lò phản ứng TRIGA Mark II — cấu hình lò nghiên cứu tiêu chuẩn quốc tế đã được kiểm chứng thực nghiệm rộng rãi tại nhiều phòng thí nghiệm trên thế giới. Kênh dẫn chùm proton được bố trí tại vị trí trung tâm lõi, đưa chùm proton $25\text{ mA}$ tương tác trực tiếp với thể tích chì lỏng.

       +-------------------------------------------------------+
       |             Chùm Proton Máy Gia Tốc                    |
       |         (Ep = 0.25 - 3.0 GeV, Ip = 25 mA)             |
       +---------------------------+---------------------------+
                                   |
                                   v
       +-------------------------------------------------------+
       |         BIA CHÌ LỎNG & CHẤT LÀM MÁT CHÍNH             |
       |        (T > 327°C, Mật độ rho = 10.66 g/cm3)          |
       |       - Tương tác Spallation (p,n) liên tục           |
       |       - Không cần dừng lò để thay thế bia vật lý       |
       +---------------------------+---------------------------+
                                   |
                +------------------+------------------+
                |                                     |
                v                                     v
+-------------------------------+   +----------------------------------+
|   VÙNG NHIÊN LIỆU PHÂN HẠCH   |   |   HỆ THỐNG PHÂN TÍCH VẬT LÝ      |
|    DƯỚI TỚI HẠN (keff < 1)    |   | - Thư viện JENDL-HE-2007 (<=3GeV)|
|  - TRIGA Mark II Geometry     |   | - Mã mô phỏng: MCNPX, GEANT4     |
|  - Nhiên liệu ThO2 pha tạp:   |   | - Rời rạc hóa: 19 góc x 32 năng  |
|    * Th-232 + U-233           |   |   lượng; Phổ vi phân d2sigma     |
|    * Th-232 + U-235           |   | - Tính toán thông lượng NF (r, z)|
|    * Th-232 + U-238           |   | - Đánh giá độ an toàn keff       |
+-------------------------------+   +----------------------------------+

Quy trình nghiên cứu rigorous

Quy trình tính toán và kiểm định diễn ra theo 4 giai đoạn chuẩn mực hạt nhân:

  1. Trích xuất và chuẩn hóa dữ liệu hạt nhân: Khai thác tập tin dữ liệu tương tác proton-hạt nhân và neutron-hạt nhân định dạng ENDF từ thư viện JENDL-HE-2007 (chứa dữ liệu của 406 hạt nhân với mức năng lượng hạt lên đến $3\text{ GeV}$).
  2. Rời rạc hóa ma trận tương tác: Xây dựng mảng ma trận phân bố $f(\mu, E_p, E_n)$ kích thước $32 \times 19$ tương ứng với 32 mức năng lượng và 19 điểm góc trong không gian bán cầu $0^\circ - 180^\circ$.
  3. Mô phỏng vận chuyển hạt và phân rã phóng xạ: Thiết lập card vật lý trong MCNPX để mô phỏng tương tác spallation, bốc hơi hạt nhân, quá trình làm chậm neutron và kích hoạt phóng xạ của chuỗi phân rã thori trong môi trường chì lỏng.
  4. Kiểm định độ giá trị và độ tin cậy (triangulation): So sánh đối chiếu kết quả vi phân bậc hai $d^2\sigma/dE d\Omega$ ở các góc $45^\circ - 60^\circ$ với số liệu thực nghiệm và tính toán của X. Ledoux, S. Meigo (1999) và A. Letourneau; kiểm chứng tính toán phân bố góc với mô hình SDM và QMD trên bia $^{208}\text{Pb}$.

Data và phân tích

Không gian pha tính toán bao phủ toàn diện:

  • Dải năng lượng hạt proton: $E_p = 250\text{ MeV}, 500\text{ MeV}, 1000\text{ MeV}, 1500\text{ MeV}, 2000\text{ MeV}, 3000\text{ MeV}$.
  • Thông số cấu trúc lõi lò TRIGA Mark II: Bán kính khảo sát $r$ từ $2\text{ cm}$ đến $24\text{ cm}$ (chia nhỏ các vị trí $8{,}75\text{ cm}, 12{,}25\text{ cm}, 15{,}75\text{ cm}, 19{,}25\text{ cm}, 22{,}75\text{ cm}, 26{,}25\text{ cm}$); chiều cao lõi lò $z$ từ $5\text{ cm}$ đến $65\text{ cm}$.
  • Khảo sát góc tán xạ: Phân tích chi tiết tại các góc $20^\circ, 25^\circ, 35^\circ, 45^\circ, 55^\circ, 60^\circ, 65^\circ, 70^\circ, 75^\circ, 80^\circ, 85^\circ$.
  • Tổ hợp nhiên liệu thử nghiệm: Nhiên liệu chuẩn $\text{UZrH}$ (làm mát bằng nước nhẹ đối chứng với làm mát bằng chì lỏng); $\text{UO}_2$ nguyên chất; và các hỗn hợp oxit chứa thori $\text{Th}^{233}\text{UO}_2$, $\text{Th}^{235}\text{UO}_2$, $\text{Th}^{238}\text{UO}_2$ tại các cấp tỷ lệ nạp Thori: $0%, 20%, 40%, 60%, 80%, 100%$.

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

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

Luận án cung cấp 5 phát hiện cốt lõi với đầy đủ bằng chứng định lượng:

  1. Hiệu suất phát neutron ($Y_{n/p}$) trên bia chì lỏng tương đương bia rắn: Dòng proton năng lượng từ $0{,}5\text{ GeV}$ đến $3{,}0\text{ GeV}$ tương tác trên thể tích chì lỏng sản sinh mật độ neutron tăng tuyến tính theo năng lượng chùm hạt. Tại $E_p = 1{,}5\text{ GeV}$, hiệu suất đạt xấp xỉ $27 - 30\text{ n/p}$. Kết quả vi phân bậc hai $d^2\sigma/dE d\Omega$ trên bia chì lỏng khớp chính xác với tính toán của nhóm tác giả X. Ledoux và cộng sự ở các góc $45^\circ - 55^\circ$ trên bia chì rắn, chứng minh rằng sự sụt giảm mật độ vật lý ($10{,}66\text{ g/cm}^3$ so với $11{,}7\text{ g/cm}^3$) hoàn toàn không làm suy giảm hiệu suất tương tác nhờ thể tích bia chất lỏng bao phủ toàn bộ đường đi chùm tia.
  2. Phân bố góc và phổ năng lượng bất đối xứng: Phổ năng lượng neutron sinh ra trải rộng từ dải neutron nhiệt ($< 0{,}625\text{ eV}$) đến neutron nhanh ($> 1\text{ MeV}$). Ở các góc tán xạ nhỏ ($0^\circ - 30^\circ$), thông lượng neutron mang động năng cao chiếm ưu thế do cơ chế va chạm trực tiếp (direct cascade); ở các góc rộng ($90^\circ - 150^\circ$), neutron bốc hơi hạt nhân (evaporation neutrons) chiếm ưu thế với năng lượng trung bình giảm rõ rệt, dao động từ $21\text{ MeV}$ xuống dưới $5\text{ MeV}$.
  3. Đặc tính phân bố thông lượng neutron 3 chiều ($r, z, E$): Thông lượng neutron cực đại tập trung tại tâm hình học của vùng bia ($z = 30 - 35\text{ cm}, r \le 5\text{ cm}$) và suy giảm dần theo hàm mũ đối xứng ra biên lò. Việc thay thế chất làm mát nước nhẹ bằng chì lỏng giúp làm cứng phổ neutron (tăng tỷ phần neutron nhanh FNF và neutron trên nhiệt), tạo điều kiện tối ưu cho phản ứng bắt neutron trên $^{232}\text{Th}$ để biến đổi thành $^{233}\text{U}$.
  4. Bản đồ phổ bức xạ phân rã Thori trong chì lỏng: Luận án đã xác định định lượng phổ năng lượng của 4 thành phần bức xạ chính phát ra từ chuỗi phóng xạ thori ($\alpha, \beta, \gamma, \bar{\nu}_e$) cùng năng lượng giật lùi (recoil energy) cực tiểu, trung bình và cực đại của các hạt nhân con trong môi trường chì lỏng, chứng minh môi trường chì đóng vai trò lớp che chắn bức xạ tự nhiên vô cùng hiệu quả.
  5. Hành vi hệ số nhân neutron hiệu dụng $k_{eff}$ giữa các hỗn hợp nhiên liệu:
    • Với hỗn hợp $\text{Th}^{233}\text{UO}2$: $k{eff}$ duy trì giá trị cao và suy giảm rất chậm khi tăng tỷ lệ Thori từ $0%$ lên $60%$ và $80%$, luôn bảo đảm trạng thái hoạt động dưới tới hạn lý tưởng ($0{,}92 \le k_{eff} \le 0{,}98$).
    • Với hỗn hợp $\text{Th}^{235}\text{UO}_2$ và $\text{Th}^{238}\text{UO}2$: $k{eff}$ suy giảm nhanh chóng khi tỷ lệ Thori vượt quá $40%$, đòi hỏi tăng mạnh cường độ dòng proton máy gia tốc để bù neutron.

Implications đa chiều

  • Về mặt lý thuyết hạt nhân: Cung cấp bộ hằng số vi phân tương tác $(p,n)$ trên chì lỏng cho cộng đồng vật lý năng lượng cao; làm sáng tỏ động học neutron dưới tới hạn trong môi trường kim loại lỏng.
  • Về mặt phương pháp luận: Thiết lập quy trình tính toán liên hợp giữa thư viện dữ liệu hạt nhân năng lượng cao JENDL-HE-2007 và mã mô phỏng MCNPX áp dụng cho cấu hình lò TRIGA Mark II cải tiến.
  • Về mặt ứng dụng công nghệ: Giải quyết triệt để vấn đề quá nhiệt và phá hủy bia rắn trong ADSR, loại bỏ hoàn toàn các chu kỳ tắt máy tốn kém để thay thế bia, nâng cao hệ số khả dụng của nhà máy điện hạt nhân lên mức tiệm cận $90 - 95%$.
  • Về mặt chính sách và môi trường: Khai thác hiệu quả nguồn tài nguyên Thori dồi dào, triệt tiêu nguy cơ phổ biến vũ khí hạt nhân (do chu trình $^{232}\text{Th}-^{233}\text{U}$ không tạo ra plutoni cấp độ vũ khí $^{239}\text{Pu}$ và luôn đi kèm đồng vị phát gamma năng lượng cao $^{232}\text{U}$), đồng thời mở ra lộ trình thiêu hủy rác thải phóng xạ actinides một cách bền vững.

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn kỹ thuật cần tiếp tục hoàn thiện:

  1. Giới hạn về mô phỏng thủy nhiệt liên hợp (coupled neutronics-thermal hydraulics): Các tính toán thông lượng neutron và $k_{eff}$ hiện được thực hiện ở trạng thái tĩnh với mật độ chì lỏng đồng nhất ($10{,}66\text{ g/cm}^3$), chưa xét đến biến thiên nhiệt độ cục bộ gây giãn nở thể tích và dòng chảy đối lưu nhiệt thực tế trong quá trình tải nhiệt $25\text{ mA}$ proton.
  2. Vấn đề ăn mòn vật liệu chì lỏng ở nhiệt độ cao: Chì lỏng ở nhiệt độ trên $400 - 500^\circ\text{C}$ có tính ăn mòn mạnh đối với vỏ bọc thanh nhiên liệu thép không gỉ, đòi hỏi công nghệ kiểm soát nồng độ oxy hòa tan nghiêm ngặt chưa được mô hình hóa trong phạm vi luận án.
  3. Hiện tượng bắt neutron ký sinh của $^{233}\text{Pa}$: Hạt nhân trung gian Protactini ($^{233}\text{Pa}$) có chu kỳ bán rã tương đối dài ($T_{1/2} = 26{,}97\text{ ngày}$); trong điều kiện thông lượng neutron rất cao, $^{233}\text{Pa}$ có thể bắt tiếp 1 neutron biến thành $^{234}\text{Pa}$ thay vì phân rã $\beta^-$ thành $^{233}\text{U}$, làm giảm nhẹ hiệu suất sinh nhiên liệu.

Chương trình nghiên cứu 5–10 năm tiếp theo được định hình:

  • Hướng 1: Phát triển mô hình mô phỏng ghép cặp đa vật lý (Multiphysics coupling) giữa MCNPX/OpenMC và hệ mã động lực học dòng chảy CFD (ANSYS Fluent/OpenFOAM) để khảo sát trường nhiệt - thủy lực động của chì lỏng.
  • Hướng 2: Nghiên cứu thực nghiệm kiểm chứng trên mô hình bia chì lỏng quy mô nhỏ kết hợp máy gia tốc proton tại các trung tâm nghiên cứu hạt nhân quốc tế.
  • Hướng 3: Đánh giá chi tiết cơ chế cháy sâu nhiên liệu (burnup credit) và động học tích tụ các sản phẩm phân hạch độc neutron ($^{135}\text{Xe}, {}^{149}\text{Sm}$) trong chu trình dài hạn $1000\text{ ngày}$.
  • Hướng 4: Nghiên cứu tối ưu hóa lớp vỏ bọc chống ăn mòn bằng vật liệu gốm hoặc hợp kim phủ nano Si/Al chịu nhiệt độ cao trong môi trường chì lỏng.

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

Công trình mang lại tác động sâu rộng trên nhiều bình diện:

  • Tác động học thuật: Đóng góp bộ dữ liệu vi phân sinh neutron trên chì lỏng vào kho tàng vật lý hạt nhân Việt Nam và quốc tế; dự kiến thu hút trích dẫn cao trong các nghiên cứu thiết kế lò phản ứng hạt nhân thế hệ IV (Lead-cooled Fast Reactor - LFR và ADS).
  • Chuyển dịch công nghiệp năng lượng: Cung cấp luận cứ khoa học vững chắc cho việc thiết kế các lò phản ứng module nhỏ an toàn nội tại (SMR-ADS), mở ra giải pháp công nghệ khả thi cho ngành công nghiệp điện hạt nhân tương lai không phụ thuộc vào quặng Urani làm giàu.
  • Định hình chính sách an toàn hạt nhân: Hỗ trợ các cơ quan quản lý nhà nước (Bộ Khoa học & Công nghệ, Cục An toàn Bức xạ và Hạt nhân) trong việc xây dựng tiêu chuẩn kỹ thuật an toàn cho các cơ sở nghiên cứu và ứng dụng công nghệ lò phản ứng tiên tiến.
  • Lợi ích xã hội và toàn cầu: Giảm thiểu triệt để rủi ro thảm họa hạt nhân dạng mất kiểm soát phản ứng dây chuyền (nhờ cơ chế tự tắt tức thì khi ngắt dòng proton máy gia tốc), giải quyết dứt điểm bài toán rác thải phóng xạ xuyên thế hệ.

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

  • Nghiên cứu sinh và học giả chuyên ngành Vật lý hạt nhân: Tiếp cận bộ dữ liệu vi phân $d^2\sigma/dE d\Omega$, phương pháp mô phỏng MCNPX trên bia lỏng và quy trình tính toán thông lượng neutron đa chiều.
  • Các kỹ sư thiết kế lò phản ứng hạt nhân: Ứng dụng mô hình bia chì lỏng hai chức năng để thiết kế hệ thống làm mát và bia spallation liên tục, loại bỏ cơ cấu cơ khí thay bia phức tạp.
  • Các tổ chức R&D công nghệ hạt nhân (VINATOM, SCK.CEN, JAEA, CAS): Kế thừa kết quả tính toán để hoàn thiện các dự án tiền khả thi về lò ADSR thương mại.
  • Nhà hoạch định chính sách năng lượng: Có cơ sở dữ liệu định lượng về tiềm năng chu trình nhiên liệu Thori an toàn, không phổ biến vũ khí hạt nhân để xây dựng chiến lược năng lượng quốc gia dài hạ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 mở rộng học thuyết nào?
Luận án mở rộng trực tiếp lý thuyết chu trình nhiên liệu Thori và khái niệm Energy Amplifier của Carlo Rubbia (1995) thông qua việc thiết lập mô hình giải tích định lượng cho tương tác spallation $(p,n)$ trên môi trường chì lỏng đối lưu liên tục, chứng minh bằng toán học và mô phỏng rằng môi trường chì lỏng vừa đảm bảo vai trò bia sinh neutron hiệu suất cao ($Y_{n/p} \approx 27-30\text{ n/p}$ tại $1{,}5\text{ GeV}$), vừa duy trì trường neutron trên nhiệt lý tưởng để biến đổi $^{232}\text{Th}$ thành $^{233}\text{U}$.

2. Đột phá phương pháp luận so với ít nhất 2 công trình quốc tế trước đây là gì?
So với nghiên cứu của S. Meigo và cộng sự (1999) và G. Bauer (2001) vốn chỉ giới hạn trên mô hình bia chì rắn tĩnh chịu giới hạn hư hỏng cơ học, luận án đã tích hợp thư viện JENDL-HE-2007 (năng lượng đến $3\text{ GeV}$) để rời rạc hóa ma trận vi phân bậc hai $32 \times 19$ trên bia chì lỏng. So với thiết kế MYRRHA của SCK.CEN (Bỉ) dùng hợp kim chì-bismuth tạo chất độc $^{210}\text{Po}$, luận án sử dụng chì lỏng nguyên chất kết hợp cấu hình lõi TRIGA Mark II, loại bỏ triệt để nguy cơ nhiễm độc bismuth mà vẫn bảo toàn thông lượng neutron mong muốn.

3. Phát hiện bất ngờ nhất được chứng minh bằng dữ liệu là gì?
Mặc dù chì lỏng có mật độ thấp hơn chì rắn khoảng $8{,}9%$ ($10{,}66\text{ g/cm}^3$ so với $11{,}7\text{ g/cm}^3$), nhưng hiệu suất sinh neutron và vi phân tiết diện $d^2\sigma/dE d\Omega$ ở các góc trung gian ($45^\circ - 60^\circ$) không hề suy giảm so với bia chì rắn của X. Ledoux. Thể tích chì lỏng trải dài trên toàn đường kính lò đã tăng quãng đường tự do tương tác của proton, bù đắp trọn vẹn sự suy giảm mật độ hạt nhân.

4. Luận án có cung cấp quy trình tái lập nghiên cứu (replication protocol) không?
Có. Luận án cung cấp đầy đủ các phương trình giải tích vi phân (2.1–2.4), thông số chùm hạt ($I_p = 25\text{ mA}, r = 4\text{ cm}, E_p = 0{,}25 - 3{,}0\text{ GeV}$), lưới rời rạc 19 giá trị góc $\mu$ và 32 kênh năng lượng $E_n$, định danh thư viện JENDL-HE-2007 định dạng ENDF, kích thước hình học chi tiết của thanh nhiên liệu và lõi lò TRIGA Mark II, cùng bảng thành phần đồng vị oxit thori-urani để các nhà nghiên cứu độc lập có thể tái lập hoàn toàn trên MCNPX/GEANT4.

5. Chương trình nghị sự nghiên cứu 10 năm tới được phác thảo như thế nào?
Lộ trình 10 năm được xác định gồm 3 giai đoạn: (1) 2022–2025: Mô phỏng thủy nhiệt liên hợp CFD-MCNP đa vật lý và tối ưu hóa hình học buồng bia chì lỏng; (2) 2025–2028: Thử nghiệm thực nghiệm tương tác proton trên bia chì lỏng quy mô phòng thí nghiệm và nghiên cứu vật liệu vỏ thanh nhiên liệu chống ăn mòn; (3) 2028–2032: Thiết kế kỹ thuật chi tiết tổ hợp lò phản ứng nghiên cứu dưới tới hạn ADSR công suất nhiệt $10 - 50\text{ MWth}$ sử dụng nhiên liệu Thori tại Việt Nam.

Kết luận

Luận án tiến sĩ của tác giả Trần Minh Tiến đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 6 đóng góp mang tính đột phá:

  1. Xây dựng thành công mô hình bia tương tác chì lỏng hai chức năng vừa đóng vai trò bia sinh neutron spallation liên tục, vừa làm môi trường tải nhiệt làm mát lõi lò, khắc phục triệt để nhược điểm phải ngừng lò thay thế bia của các ADSR bia rắn truyền thống.
  2. Tính toán và kiểm chứng chính xác bộ tham số vi mô của tương tác $(p,n)$ trên chì lỏng ở dải năng lượng proton $0{,}25 - 3{,}0\text{ GeV}$, bao gồm phân bố năng lượng, phân bố góc, hiệu suất phát neutron $Y_{n/p}$ và vi phân bậc hai tiết diện $d^2\sigma/dE d\Omega$ dựa trên thư viện JENDL-HE-2007.
  3. Định lượng hóa toàn diện đặc trưng phân rã phóng xạ của Thori trong chì lỏng, xác định chi tiết phổ năng lượng các bức xạ $\alpha, \beta, \gamma, \bar{\nu}_e$ và động năng giật lùi của các hạt nhân con.
  4. Xác lập bản đồ phân bố thông lượng neutron 3 chiều ($r, z, E$) trong lõi ADSR dựa trên cấu trúc TRIGA Mark II, chứng minh môi trường làm mát chì lỏng tạo ra phổ neutron tối ưu cho quá trình chuyển hóa $^{232}\text{Th}$ thành $^{233}\text{U}$.
  5. Định lượng chính xác độ phản ứng dưới tới hạn ($k_{eff}$) của các tổ hợp nhiên liệu $\text{Th}^{233}\text{UO}_2$, $\text{Th}^{235}\text{UO}_2$, $\text{Th}^{238}\text{UO}2$, khẳng định ưu thế vượt trội của hỗn hợp $^{232}\text{Th}-^{233}\text{U}$ trong việc duy trì trạng thái an toàn $k{eff} \approx 0{,}95 - 0{,}98$ ở tỷ lệ nạp Thori cao.
  6. Mở ra 3 hướng nghiên cứu ứng dụng mới: Công nghệ lò SMR-ADS làm mát bằng chì lỏng; kỹ thuật thiêu hủy rác thải phóng xạ actinide năng lượng cao; và lộ trình tự chủ nhiên liệu hạt nhân bền vững từ nguồn quặng Thori trong nước. Công trình là một bước tiến quan trọng đóng góp vào sự phát triển chung của ngành vật lý hạt nhân và công nghệ năng lượng nguyên tử hiện đại.