Tổng quan về luận án

Năng lượng hạt nhân giữ vai trò trụ cột trong chiến lược an ninh năng lượng toàn cầu, đặc biệt trong bối cảnh các nguồn nhiên liệu hóa thạch truyền thống đang cạn kiệt nhanh chóng và sức ép giảm phát thải carbon ngày càng cấp bách. Tuy nhiên, ngành công nghiệp điện hạt nhân thương mại hiện nay – vốn chủ yếu vận hành dựa trên các lò phản ứng nhiệt truyền thống sử dụng nhiên liệu Urani ($^{235}\text{U}$) – đang đối mặt với ba nút thắt mang tính sống còn: nguy cơ rủi ro an toàn tới hạn (thể hiện qua các thảm họa lịch sử tại Chernobyl và Fukushima), trữ lượng quặng Urani kinh tế giới hạn, và bài toán quản trị các chất thải phóng xạ có đời sống cực dài (các nguyên tố actini hiếm - Minor Actinides như $\text{Np}, \text{Am}, \text{Cm}$). Để giải quyết căn bản những thách thức này, nghiên cứu phát triển 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) kết hợp chu trình nhiên liệu Thori ($^{232}\text{Th}$) và công nghệ bia kim loại lỏng được xem là một hướng đi đột phá mang tính chuyển dịch mô thức (paradigm shift) cho ngành vật lý lò phản ứng.

Luận án tiến sĩ vật lý chuyên ngành Vật lý nguyên tử và hạt nhân 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" đã giải quyết trực diện một khoảng trống nghiên cứu cốt lõi (research gap): Trong các thiết kế ADSR truyền thống, việc sử dụng bia tán xạ phá vỡ hạt nhân (spallation target) ở trạng thái rắn (như thanh chì rắn $\text{Pb}$, vonfram $\text{W}$) thường phải đối mặt với hiện tượng tổn thương bức xạ nghiêm trọng, suy giảm cơ tính và bắt buộc phải định kỳ dừng lò để thay thế bia, gây gián đoạn chu trình vận hành liên tục. Mặt khác, các nghiên cứu trước đây về bia kim loại lỏng (như $\text{Pb}$ hoặc $\text{Pb-Bi}$) chủ yếu chỉ khoanh vùng trong các module chuyển hóa chất thải phóng xạ quy mô nhỏ, chưa thực hiện tính toán toàn diện các tham số neutron đa chiều khi tích hợp chì lỏng đồng thời làm bia bắn phá và chất tải nhiệt sơ cấp kết hợp với việc nạp nhiên liệu Thori trong toàn bộ cấu trúc lõi lò dưới tới hạn.

Để giải quyết khoảng trống này, luận án tập trung vào các câu hỏi nghiên cứu (Research Questions - RQ) và giả thuyết khoa học (Hypotheses - H) cụ thể:

  • RQ1: Tương tác hạt nhân $(p,n)$ trên môi trường bia chì lỏng dưới tác động của chùm proton năng lượng cao ($0,25\text{ GeV} - 3,0\text{ GeV}$) tạo ra các đặc trưng phân bố năng lượng, phân bố góc và vi phân bậc hai tiết diện sinh neutron như thế nào?
  • RQ2: Khi cấu hình chì lỏng vừa làm bia tương tác vừa làm chất làm mát trong lõi lò TRIGA Mark II dưới tới hạn, phân bố thông lượng neutron (theo bán kính, chiều cao và phổ năng lượng) thay đổi ra sao so với lò truyền thống?
  • RQ3: Mức độ phản ứng và hệ số nhân neutron hiệu dụng ($k_{eff}$) của các hỗn hợp nhiên liệu Thori kết hợp Urani ($\text{Th}^{233}\text{UO}_2$, $\text{Th}^{235}\text{UO}_2$, $\text{Th}^{238}\text{UO}_2$) biến thiên theo tỷ lệ thành phần như thế nào nhằm duy trì chế độ dưới tới hạn tối ưu?

Các giả thuyết nghiên cứu bao gồm:

  • H1: Bia chì lỏng với chiều dài tương tác trải rộng toàn bộ đường kính vùng chiếu xạ sẽ gia tăng hiệu suất phát neutron tổng cộng ($Y_{n/p}$) và loại bỏ hoàn toàn yêu cầu dừng lò thay thế bia cơ học.
  • H2: Quá trình chuyển hóa Thori ($^{232}\text{Th} + n \rightarrow ^{233}\text{Th} \xrightarrow{\beta^-} ^{233}\text{Pa} \xrightarrow{\beta^-} ^{233}\text{U}$) trong môi trường phổ neutron nhiệt và nhanh của ADSR chì lỏng cung cấp đủ khả năng nhân neutron để duy trì $k_{eff}$ ổn định trong dải dưới tới hạn an toàn ($0,95 \le k_{eff} < 0,98$).

Khung lý thuyết của công trình tích hợp lý thuyết phản ứng phá vỡ hạt nhân năng lượng cao (Spallation Reaction Mechanics), lý thuyết vận chuyển hạt Monte Carlo ngẫu nhiên và động học lò phản ứng dưới tới hạn. Luận án benchmark hình học dựa trên cấu trúc chuẩn của lò phản ứng nghiên cứu TRIGA Mark II, khảo sát dải năng lượng proton tới từ $250\text{ MeV}$ đến $3000\text{ MeV}$, cường độ dòng proton $25\text{ mA}$, bán kính chùm $4\text{ cm}$, lưới góc khối 18 khoảng từ $0^\circ$ đến $180^\circ$ và phổ năng lượng neutron chia thành 32 nhóm.

Literature Review và Positioning

Ý tưởng sử dụng máy gia tốc hạt năng lượng cao để kích hoạt và duy trì phản ứng phân hạch hạt nhân bắt nguồn từ những quan sát nền tảng của Ernest Rutherford (1919) khi bắn phá hạt nhân $^{14}\text{N}$ bằng hạt $\alpha$, tiếp nối bởi các sáng chế máy gia tốc cyclotron của Ernest O. Lawrence (1930) và các đề xuất độc lập của Lawrence cùng Semenov (1940). Trong giai đoạn thập niên 1990, khái niệm ADSR thực sự trở thành tâm điểm khoa học toàn cầu thông qua các công trình tiên phong của C. Bowman cùng cộng sự (1992) [3] tại Phòng thí nghiệm Quốc gia Los Alamos (LANL) và công trình đột phá về "Bộ khuếch đại năng lượng" (Energy Amplifier) của Carlo Rubbia cùng cộng sự (1995) [1, 12] tại CERN. Rubbia đã chứng minh khả năng vận hành một hệ thống neutron nhanh công suất nhiệt $600\text{ MWth}$ làm mát bằng chì, sử dụng nhiên liệu Thori tự nhiên trong trạng thái dưới tới hạn nội tại.

Trong dòng nghiên cứu tính toán đặc trưng neutron trên bia bắn phá, S. Meigo và cộng sự (1999) [34] đã sử dụng mã MCNP4A để mô phỏng phân bố thông lượng neutron từ bia chì dày dưới chùm proton $0,5\text{ GeV}$, chỉ ra rằng phổ neutron sinh ra trải rộng từ dải nhiệt đến dải nhanh. Boudard và cộng sự (2000) [33] xác định mối quan hệ nghịch đảo giữa khoảng cách hình học và năng lượng neutron trung bình ở chùm proton $0,6\text{ GeV}$. Letourneau, Goldenbaum và cộng sự (2000) [35] đã thiết lập các giới hạn độ dày tới hạn đối với bia tán xạ hạt nhân rắn: $30\text{ cm}$ đối với Vonfram ($\text{W}$) và $55\text{ cm}$ đối với Chì ($\text{Pb}$). G. Bauer (2001) [36] phân tích chi tiết đặc trưng vật lý của chùm proton $2\text{ GeV}$ bắn phá bia chì dày $20\text{ cm}$, phát hiện năng lượng neutron trung bình tại các góc $30^\circ, 90^\circ, 150^\circ$ lần lượt đạt giá trị khoảng $21\text{ MeV}$, và dải góc từ $10^\circ - 90^\circ$ sinh ra neutron có năng lượng vượt trội so với dải $90^\circ - 150^\circ$. David Sangcheol Lee (2018) [39] đã hoàn thiện các phép đo và mô phỏng hiệu suất phát neutron, phân bố thông lượng trên bia chì và chì-bismuth ở năng lượng $1,5\text{ GeV}$.

Tại Việt Nam, trường phái nghiên cứu ADSR được khởi xướng bởi Nguyễn Mộng Giao và các cộng sự từ năm 2002. Đáng chú ý, nghiên cứu của Nguyễn Thị Ái Thu, Châu Văn Tạo và Nguyễn Mộng Giao [16-19] đã xây dựng thành công "mô hình màn chắn trên bia" để hiệu chỉnh sai số giữa lý thuyết và thực nghiệm cho các tương tác $(p,n)$ trên các bia nặng ($\text{U}, \text{Pb}, \text{W}, \text{Au}$) trong dải năng lượng $0,5 - 1,5\text{ GeV}$. Tiếp đó, các công trình của Vũ Thanh Mai cùng nhóm nghiên cứu Cheol Ho Pyeon, Masao Yamanaka tại Tổ hợp tới hạn Đại học Kyoto (KUCA) và JAEA (Nhật Bản) [20-25] đã đóng góp các dữ liệu thực nghiệm chuẩn mực về động học neutron trễ và neutron tức thời khi sử dụng máy gia tốc proton $100\text{ MeV}$ trên bia kim loại nặng kết hợp nhiên liệu $^{235}\text{U}$ và $^{232}\text{Th}$.

                 TIẾN TRÌNH PHÁT TRIỂN & ĐỊNH VỊ NGHIÊN CỨU
===================================================================================
1940-1990: Khởi nguồn ý tưởng      ---> Lawrence, Semenov, Seaborg, C. Bowman (1992)
1995: Mô hình Energy Amplifier    ---> Carlo Rubbia (CERN) - ADS dùng Thori & Chì
2000-2018: Nghiên cứu Bia rắn      ---> Meigo, Bauer, Letourneau, D. S. Lee (MCNP)
2002-2020: Nghiên cứu tại VN/QT    ---> Nguyễn Mộng Giao, Nguyễn Thị Ái Thu, Vũ Thanh Mai (KUCA)
2022: LUẬN ÁN TRẦN MINH TIẾN       ---> ĐỘT PHÁ: BIA CHÌ LỎNG ĐỒNG THỜI LÀM CHẤT TẢI NHIỆT
                                         + NHIÊN LIỆU THORI ĐA DẠNG (Th-U) TRONG TRIGA MARK II
===================================================================================

Các cuộc tranh luận học thuật quốc tế hiện nay tập trung vào hai đối cực:

  1. Lựa chọn chất làm mát và vật liệu bia: Tranh luận giữa sử dụng hợp kim Chì-Bismuth ($\text{Pb-Bi}$, nhiệt độ nóng chảy thấp $123,5^\circ\text{C}$ nhưng sinh ra đồng vị độc hại $^{210}\text{Po}$ dễ bay hơi và gây ăn mòn kim loại) so với Chì nguyên chất ($\text{Pb}$, nhiệt độ nóng chảy $327^\circ\text{C}$, mật độ khối chất lỏng $10,66\text{ g}\cdot\text{cm}^{-3}$, an toàn phóng xạ cao hơn nhưng đòi hỏi hệ thống duy trì nhiệt độ khắt khe).
  2. Chu trình nhiên liệu Thori trong lò phản ứng nhanh so với lò nhiệt: Trong khi các nghiên cứu thế hệ IV (như SFR, LFR) đánh giá chu trình $^{232}\text{Th}-^{233}\text{U}$ gặp trở ngại do chu kỳ bán rã dài của $^{233}\text{Pa}$ ($T_{1/2} \approx 27\text{ ngày}$) dẫn đến hiện tượng bắt neutron ký sinh làm giảm tốc độ sinh sản, thì hệ thống ADSR lại khắc phục triệt để nhờ nguồn neutron ngoại lai bù đắp liên tục từ phản ứng phá vỡ hạt nhân.

So sánh với các đại dự án quốc tế:

  • Dự án MYRRHA (SCK.CEN, Bỉ) [29]: Sử dụng máy gia tốc $600\text{ MeV} - 4\text{ mA}$, bia $\text{Pb-Bi}$, công suất $70\text{ MWth}$, thông lượng neutron nhanh $10^{15}\text{ n}\cdot\text{cm}^{-2}\text{s}^{-1}$ chủ yếu phục vụ chiếu xạ vật liệu và chuyển hóa chất thải.
  • Dự án C-ADS (Viện Hàn lâm Khoa học Trung Quốc - CAS) [31]: Hợp tác giữa IHEP, IMP, IPP và USTC với lộ trình $100\text{ MW}$ (2022) tiến tới $1\text{ GW}$ (2032).
  • Hệ thống KIPT (Kharkov, Ukraine hợp tác với ANL, Hoa Kỳ) [32]: Máy gia tốc tuyến tính kết hợp tổ hợp dưới tới hạn sử dụng Urani độ giàu thấp (LEU) làm mát bằng nước và vành phản xạ $\text{Be-C}$.

Luận án của Trần Minh Tiến định vị độc đáo khi chọn giải pháp thiết kế tích hợp: dùng chì lỏng nguyên chất vừa làm bia bắn phá $(p,n)$ vừa làm chất tải nhiệt đối lưu cho toàn bộ lõi lò TRIGA Mark II dưới tới hạn nạp nhiên liệu Thori.

Đó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 sinh và truyền tải neutron năng lượng cao trong môi trường vật chất lỏng tỷ trọng lớn. Về mặt toán lý hạt nhân, tiết diện sinh neutron vi phân bậc hai theo năng lượng và góc khối trong phản ứng phá vỡ hạt nhân $(p,n)$ được mô hình hóa chính xác theo công thức:

$$\frac{d^2 \sigma (\mu, E_p, E_n)}{dE_n , d\Omega} = \sigma (E_p) \cdot \frac{d\gamma (E_p)}{d\Omega} \cdot f(\mu, E_p, E_n)$$

Trong đó:

  • $E_p$ là năng lượng của hạt proton tới ($\text{eV}$).
  • $E_n$ là năng lượng của neutron phát ra ($\text{eV}$).
  • $\sigma(E_p)$ là tiết diện toàn phần của phản ứng ($\text{barn}$).
  • $\gamma(E_p)$ là hiệu suất phát neutron tích hợp.
  • $f(\mu, E_p, E_n)$ là hàm phân bố xác suất hai biến biểu diễn tương quan giữa cosin góc phát $\mu = \cos\theta$ và năng lượng neutron $E_n$.

Mật độ dòng neutron sinh ra theo phổ năng lượng được chuẩn hóa qua biểu thức:

$$\frac{dN_n}{dE_n} = N_p \cdot \frac{d\sigma}{dE_n}$$

Tổng số neutron sinh ra $N_n$ từ chùm hạt proton $N_p$ được tích phân gián đoạn trên không gian pha gồm 32 nhóm năng lượng $E_j$ và 19 điểm nút góc khối $\mu_i$:

$$N_n = N_p \sum_{j=1}^{32} \left[ \frac{d\sigma (\mu, E_p, E_n)}{dE_n} \right]_j \Delta E_j$$

Mô hình lý thuyết này giải thích rõ ràng sự dịch chuyển phổ từ các neutron spallation năng lượng cao định hướng hẹp phía trước ($0^\circ - 30^\circ$) sang các neutron bay hơi (evaporation neutrons) phân bố đẳng hướng ở năng lượng thấp hơn, cung cấp cơ sở giải tích chuẩn mực cho các hệ hạt nhân dưới tới hạn.

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

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

  1. Lý thuyết tương tác nội hạt nhân (Intra-Nuclear Cascade - INC) và mô hình phân rã thống kê (Statistical Decay Model - SDM): Mô tả giai đoạn proton kích thích hạt nhân chì giải phóng hạt nhanh và giai đoạn hạt nhân kích thích nguội đi bằng cách bay hơi neutron.
  2. Lý thuyết vận chuyển hạt Monte Carlo ngẫu nhiên: Xử lý quá trình làm chậm, khuếch tán, hấp thụ và phân hạch neutron qua ma trận nhiên liệu và chất làm mát.
  3. Động học biến đổi chu trình nhiên liệu Thori: Mô phỏng sự tích tụ $^{233}\text{U}$ và phát xạ phóng xạ thứ cấp ($\alpha, \beta^-, \gamma, \bar{\nu}$).
                     KHUNG PHÂN TÍCH TÍCH HỢP TOÀN DIỆN
+-------------------------------------------------------------------------+
|                  CHÙM PROTON NĂNG LƯỢNG CAO (Máy gia tốc)               |
|                  Ep = 250 - 3000 MeV | I = 25 mA | R = 4 cm             |
+------------------------------------+------------------------------------+
                                     |
                                     v
+-------------------------------------------------------------------------+
|                BIA CHÌ LỎNG (Vừa làm Bia vừa làm Tải nhiệt)             |
|   Tương tác phá vỡ hạt nhân (p,n) -> Vi phân bậc hai: d²σ/(dE·dΩ)       |
|   Dữ liệu: JENDL-HE-2007 (ENDF Format) | Mật độ Pb = 10,66 g/cm³        |
+------------------------------------+------------------------------------+
                                     |
                                     v
+-------------------------------------------------------------------------+
|              LÕI LÒ DƯỚI TỚI HẠN (Mô hình TRIGA Mark II)                |
|   Nhiên liệu: ThUO₂ hỗn hợp (Th-233U, Th-235U, Th-238U)                 |
|   Mô phỏng MCNPX: Thông lượng Φ(r, z, E), Hiệu suất Yn/p, Hệ số keff    |
+------------------------------------+------------------------------------+
                                     |
                                     v
+-------------------------------------------------------------------------+
|                     ĐÁNH GIÁ ĐẶC TRƯNG VẬT LÝ                           |
|   - An toàn tuyệt đối: keff < 1 (Lò tắt tức thì khi ngắt máy gia tốc)   |
|   - Vận hành liên tục: Không cần dừng lò thay thế bia rắn               |
|   - Không phát sinh Plutoni vũ khí | Tận dụng tài nguyên Thori vô tận   |
+-------------------------------------------------------------------------+

Điều kiện biên của mô hình được xác lập chặt chẽ: Vùng hoạt động duy trì nghiêm ngặt dưới tới hạn ($k_{eff} < 1$), chùm proton công suất cao được dẫn qua ống chân không chiếu thẳng đứng vào tâm bể chì lỏng, vành phản xạ neutron bằng than chì và berylli bọc ngoài nhằm giảm thiểu tổn thất neutron rò rỉ.

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ủ triết lý nghiên cứu thực chứng (positivism) dựa trên phương pháp luận mô phỏng số học ngẫu nhiên (stochastic simulation) kết hợp đánh giá dữ liệu hạt nhân thực nghiệm chuẩn quốc tế. Cấu trúc mô phỏng được phân tầng đa cấp:

  • Cấp độ 1 (Cấp vi mô hạt nhân): Tính toán tiết diện phản ứng $(p,n)$, phân bố góc và phổ năng lượng neutron sinh ra từ tương tác proton đơn năng trên hạt nhân chì.
  • Cấp độ 2 (Cấp tế bào nhiên liệu và thanh bia): Phân tích phổ phân rã phóng xạ của Thori trong môi trường chì lỏng và quá trình hấp thụ neutron chuyển hóa thành $^{233}\text{U}$.
  • Cấp độ 3 (Cấp toàn lõi lò phản ứng - Core level): Tích hợp cấu trúc hình học lò TRIGA Mark II, phân bố các bó thanh nhiên liệu, chất làm mát chì lỏng và vành phản xạ để xác định thông lượng neutron 3 chiều $\Phi(r, z, E)$ và hệ số nhân neutron $k_{eff}$.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu sử dụng hệ thống công cụ và thư viện dữ liệu tiêu chuẩn cao nhất trong ngành vật lý hạt nhân:

  1. Thư viện dữ liệu JENDL-HE-2007 (Nhật Bản): Đây là thư viện dữ liệu hạt nhân năng lượng cao chuẩn quốc tế của JAEA, lưu trữ tiết diện phản ứng của các đồng vị hạt nhân tương tác với neutron và proton lên đến $3\text{ GeV}$ theo định dạng chuẩn ENDF-6.
  2. Chương trình mô phỏng MCNPX và MCNP5: Sử dụng kỹ thuật Monte Carlo theo dõi lịch sử chuyển vận của hàng chục triệu hạt ($NPS \ge 10^7$) để đảm bảo sai số thống kê tương đối luôn dưới $1 - 2%$.
  3. Mã GEANT4: Dùng để kiểm chứng chéo phổ phân bố các bức xạ hạt nhân con và hạt thứ cấp.
                           QUY TRÌNH MÔ PHỎNG VẬT LÝ
[Dữ liệu JENDL-HE-2007] ---> [Xử lý ma trận tiết diện σ(Ep), f(μ, Ep, En)]
                                          |
                                          v
[Thiết lập hình học TRIGA Mark II] -> [Mô phỏng MCNPX/MCNP5 (NPS > 10^7)]
                                          |
                                          v
    +-------------------------------------+-----------------------------------+
    |                                     |                                   |
    v                                     v                                   v
[Phân bố thông lượng Φ(r,z,E)]   [Hiệu suất phát Yn/p]           [Hệ số keff theo % Th]

Phương pháp tam giác đạc lý thuyết và mô hình (model triangulation) được thực hiện bằng cách đối sánh kết quả tính toán trên bia chì lỏng với các mô hình hạt nhân kinh điển: Mô hình phân rã thống kê (SDM - Statistical Decay Model) và Mô hình động lực học phân tử lượng tử (QMD - Quantum Molecular Dynamics) trên đồng vị $^{208}\text{Pb}$.

Data và phân tích

Dữ liệu đầu vào và các tham số kỹ thuật được định lượng chi tiết:

  • Thông số chùm hạt proton: Năng lượng $E_p = 250, 500, 1000, 1500, 2000, 3000\text{ MeV}$; dòng điện chùm $I_p = 25\text{ mA}$; bán kính chùm proton $R_{beam} = 4\text{ cm}$.
  • Thông số vật liệu bia và tải nhiệt: Mật độ khối chì lỏng $\rho_{\text{Pb_liquid}} = 10,66\text{ g}\cdot\text{cm}^{-3}$ (so với chì rắn $\rho_{\text{Pb_solid}} = 11,7\text{ g}\cdot\text{cm}^{-3}$); nhiệt độ nóng chảy $327^\circ\text{C}$.
  • Cấu hình nhiên liệu khảo sát:
    • Nhiên liệu truyền thống: $\text{UZrH}$ làm mát bằng nước nhẹ hoặc chì lỏng.
    • Nhiên liệu oxit urani: $\text{UO}_2$.
    • Nhiên liệu hỗn hợp Thori-Urani: $\text{Th}^{233}\text{UO}_2$, $\text{Th}^{235}\text{UO}_2$, $\text{Th}^{238}\text{UO}_2$ với các tỷ lệ Thori thay đổi từ $0%, 20%, 40%, 60%, 80%$ đến $100%$.
  • Lưới phân tích không gian: Chiều cao lõi lò khảo sát từ $z = 5\text{ cm}$ đến $z = 65\text{ cm}$; bán kính lõi lò khảo sát từ $r = 2\text{ cm}$ đến $r = 26,25\text{ cm}$ (tại các vị trí $8,75; 12,25; 15,75; 19,25; 22,75; 26,25\text{ cm}$); lưới góc phát xạ phân tích tại $20^\circ, 25^\circ, 35^\circ, 45^\circ, 55^\circ, 60^\circ, 65^\circ, 70^\circ, 75^\circ, 80^\circ, 85^\circ$.

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

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

Thứ nhất, nghiên cứu đã chứng minh tính ưu việt vượt trội của bia chì lỏng so với bia rắn về mặt vật lý hạt nhân và kỹ thuật nhiệt. Dù chì lỏng có mật độ khối thấp hơn chì rắn ($10,66\text{ g}\cdot\text{cm}^{-3}$ so với $11,7\text{ g}\cdot\text{cm}^{-3}$), nhưng do toàn bộ thể tích chì đối lưu dọc theo đường đi của chùm proton đóng vai trò làm bia tương tác, chiều dài tương tác thực tế tăng lên đáng kể, dẫn đến hiệu suất phát neutron $Y_{n/p}$ tăng tỷ lệ thuận với năng lượng chùm proton từ $500\text{ MeV}$ đến $3000\text{ MeV}$. Kết quả so sánh vi phân bậc hai tiết diện sinh neutron $\frac{d^2\sigma}{dEd\Omega}$ tại góc $60^\circ$ hoàn toàn phù hợp với tính toán thực nghiệm trên bia chì rắn của X. Ledoux và cộng sự (ở góc $45^\circ - 55^\circ$).

Thứ hai, phổ phân rã phóng xạ của Thori trong môi trường chì lỏng đã được định lượng hóa toàn diện. Luận án xác định chính xác phổ năng lượng của 4 loại bức xạ chính ($\alpha, \beta, \gamma, \bar{\nu}$) cùng các mức năng lượng cực tiểu, trung bình và cực đại của các hạt nhân con sinh ra trong chuỗi phân rã phóng xạ Thori, cung cấp bộ hằng số phóng xạ thiết yếu cho việc tính toán che chắn an toàn sinh học.

Thứ ba, phân bố thông lượng neutron trong lõi ADSR sử dụng chì lỏng và nhiên liệu Thori thể hiện tính đồng nhất và tối ưu hóa không gian.

  • Dọc theo chiều cao lò ($z$), thông lượng neutron đạt giá trị cực đại tại vùng tâm lò ($z = 30 - 40\text{ cm}$) và giảm dần đối xứng về hai đầu biên theo dạng hàm cosin chuẩn.
  • Dọc theo bán kính ($r$), thông lượng neutron đạt đỉnh tại tâm bia chì lỏng ($r = 2 - 4\text{ cm}$) và suy giảm theo hàm Bessel khi tiến ra vùng ngoài ($r = 26,25\text{ cm}$).
  • Phổ năng lượng neutron bên trong lò duy trì cấu trúc hai vùng rõ rệt: đỉnh neutron nhiệt phục vụ chuyển hóa $^{232}\text{Th} \rightarrow ^{233}\text{U}$ và đuôi neutron nhanh ($E_n > 1\text{ MeV}$) trực tiếp kích thích phản ứng phân hạch, tương thích cao hơn hẳn so với cấu hình làm mát bằng nước nhẹ truyền thống.

Thứ tư, quy luật biến thiên hệ số nhân neutron hiệu dụng ($k_{eff}$) theo tỷ lệ Thori đã được xác lập:

  • Đối với hỗn hợp $\text{Th}^{233}\text{UO}2$: Do $^{233}\text{U}$ có tiết diện phân hạch cao và số neutron sinh ra trên mỗi phân hạch lớn ($\eta \approx 2,3$), $k{eff}$ duy trì ở mức cao và chỉ giảm nhẹ khi tăng tỷ lệ Thori. Ngay cả ở tỷ lệ Thori cao ($60% - 80%$), hệ thống vẫn duy trì hoàn hảo trạng thái dưới tới hạn mong muốn ($k_{eff} \approx 0,95 - 0,98$).
  • Đối với hỗn hợp $\text{Th}^{235}\text{UO}2$: $k{eff}$ suy giảm nhanh hơn khi tăng hàm lượng $^{232}\text{Th}$, đòi hỏi tỷ lệ bù neutron lớn hơn từ máy gia tốc.
  • Đối với hỗn hợp $\text{Th}^{238}\text{UO}2$: $k{eff}$ suy giảm mạnh nhất do cả hai đồng vị $^{232}\text{Th}$ và $^{238}\text{U}$ đều là vật liệu màu (fertile material), không thể tự duy trì phản ứng dây chuyền nếu không có độ giàu $^{235}\text{U}$ ban đầu thích hợp.
                           BẢNG SO SÁNH ĐẶC TRƯNG HỆ THỐNG
+-------------------------+-------------------------+-------------------------+
|     THAM SỐ KHẢO SÁT    |    LÒ TRADITIONAL (LWR) |    MÔ HÌNH ADSR (LUẬN ÁN)|
+-------------------------+-------------------------+-------------------------+
| Nhiên liệu chính        | Urani giàu (U-235)      | Thori (Th-232) + U-233  |
| Trạng thái tới hạn      | Tới hạn (keff = 1,0)    | Dưới tới hạn (keff < 1) |
| Nguy cơ mất kiểm soát   | Có (Chernobyl/Fukushima)| Hoàn toàn triệt tiêu    |
| Bản chất Bia bắn phá    | Không có (Tự duy trì)   | Chì lỏng (Tải nhiệt kép)|
| Thời gian dừng thay bia | Không áp dụng           | 0 ngày (Bia lỏng vĩnh cửu)|
| Sản sinh Plutoni vũ khí | Lượng lớn (Pu-239)      | Hầu như không có        |
| Xử lý Minor Actinides   | Rất hạn chế             | Khả năng chuyển hóa cao |
+-------------------------+-------------------------+-------------------------+

Implications đa chiều

  • Ý nghĩa lý thuyết: Cung cấp bộ dữ liệu giải tích và số học chuẩn mực về tương tác $(p,n)$ trên bia chì lỏng năng lượng cao ($0,25 - 3\text{ GeV}$) bổ sung vào cơ sở dữ liệu vật lý lò phản ứng hạt nhân thế hệ mới.
  • Đổi mới phương pháp luận: Thiết lập quy trình mô phỏng liên kết hoàn chỉnh từ dữ liệu vi mô (JENDL-HE-2007) qua công cụ trung gian (MCNPX/MCNP5) đến phân tích lõi lò TRIGA Mark II phức tạp.
  • Ứng dụng thực tiễn: Mở đường cho thiết kế các lò phản ứng nghiên cứu và module phát điện hạt nhân quy mô nhỏ (SMR-ADS) an toàn tuyệt đối, loại bỏ cơ chế nổ tới hạn và giải quyết bài toán cạn kiệt quặng Urani bằng nguồn tài nguyên Thori phong phú ($1,2 \times 10^{14}\text{ tấn}$ trên toàn cầu).
  • Khuyến nghị chính sách: Cung cấp luận cứ khoa học vững chắc để các cơ quan năng lượng nguyên tử quốc gia xây dựng chiến lược phát triển điện hạt nhân thế hệ mới, định hướng đầu tư nghiên cứu các tổ hợp máy gia tốc chùm proton công suất trung bình kết hợp lò hạt nhân nghiên cứu.

Limitations và Future Research

Mặc dù luận án đã đạt được những kết quả đột phá mang tính hệ thống, một số giới hạn nghiên cứu (boundary conditions) cần được nhìn nhận khách quan:

  1. Giới hạn ghép nối thủy động lực học nhiệt (Thermal-Hydraulic Coupling): Mô phỏng MCNPX chủ yếu tập trung vào vật lý hạt nhân (neutronics) ở trạng thái dừng, chưa ghép nối đồng thời (coupled simulation) với các phương trình thủy động lực học chất lưu (CFD) để khảo sát trường vận tốc, hiện tượng đối lưu nhiệt tự nhiên và sự phân bố nhiệt độ không đồng đều của chì lỏng trong điều kiện công suất cao.
  2. Hiện tượng ăn mòn vật liệu cấu trúc: Chì lỏng ở nhiệt độ trên $327^\circ\text{C}$ có tính ăn mòn hóa học và cơ học rất mạnh đối với vỏ bọc thanh nhiên liệu và thành thùng lò kim loại; luận án chưa đi sâu vào động học ăn mòn và độ bền vật liệu dài hạn.
  3. Mô phỏng cháy nhiên liệu (Burnup Analysis) theo thời gian thực: Luận án khảo sát các trạng thái tĩnh của hỗn hợp nhiên liệu với các tỷ lệ khác nhau, chưa mô phỏng liên tục chu trình cháy nhiên liệu trong 1000 - 1500 ngày hiệu dụng để đánh giá chi tiết sự tích tụ các sản phẩm phân hạch gây ngộ độc neutron (như $^{135}\text{Xe}, ^{149}\text{Sm}$) và động thái chuyển hóa trung gian của $^{233}\text{Pa}$ ($T_{1/2} = 27\text{ ngày}$).
  4. Kiểm chứng thực nghiệm: Do điều kiện cơ sở hạ tầng trong nước chưa có máy gia tốc proton năng lượng $\text{GeV}$ kết hợp tổ hợp tới hạn, các kết quả mới dừng lại ở mô phỏng số và so sánh chéo với các công bố quốc tế.

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Phát triển mô hình mô phỏng đa vật lý (Multiphysics) ghép nối mã neutronics (MCNPX/OpenMC) với mã nhiệt thủy lực CFD (ANSYS Fluent/OpenFOAM).
  • Nghiên cứu cơ chế chiết tách trực tiếp Protactini ($^{233}\text{Pa}$) ra khỏi vùng chiếu xạ để tối ưu hóa hiệu suất chuyển hóa thành $^{233}\text{U}$ tinh khiết.
  • Nghiên cứu động học chuyển tiếp khi chùm máy gia tốc bị ngắt đột ngột (beam trip transients) và ảnh hưởng đến ứng suất nhiệt của cấu trúc lò.
  • Đánh giá khả năng chuyển hóa các đồng vị actinide hiếm ($\text{Am}, \text{Cm}$) khi pha trộn trực tiếp vào ma trận nhiên liệu Thori trong ADSR chì lỏng.

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

  • Ảnh hưởng học thuật: Công trình đóng góp trực tiếp vào kho tàng tri thức vật lý hạt nhân Việt Nam và quốc tế, tạo tiền đề trích dẫn cho các nghiên cứu về hệ thống năng lượng dưới tới hạn, đồng thời khẳng định năng lực tính toán mô phỏng vật lý lò phản ứng tiên tiến của Viện Năng lượng Nguyên tử Việt Nam.
  • Chuyển đổi công nghiệp năng lượng: Cung cấp thiết kế ý niệm hoàn chỉnh cho ngành công nghiệp hạt nhân tương lai, nơi các lò ADSR nhiên liệu Thori có thể vừa sản xuất điện năng vừa đóng vai trò như các "nhà máy thiêu hủy chất thải phóng xạ" độc hại từ các lò phản ứng truyền thống.
  • Tác động chính sách: Định hình tầm nhìn chiến lược cho các nhà hoạch định chính sách năng lượng về việc tiếp cận công nghệ hạt nhân an toàn nội tại (inherent passive safety), xua tan lo ngại xã hội sau sự cố Fukushima.
  • Lợi ích xã hội và môi trường: Giảm thiểu hàng nghìn lần khối lượng chất thải phóng xạ có thời gian phân rã hàng trăm nghìn năm, rút ngắn chu kỳ quản lý địa chất rác thải hạt nhân xuống dưới 300 năm, bảo vệ môi trường sinh thái bền vững.

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

  • Nghiên cứu sinh và Giảng viên Vật lý Hạt nhân: Tiếp cận được quy trình mô phỏng chuẩn tắc kết hợp dữ liệu JENDL-HE-2007 và mã MCNPX, làm tài liệu tham khảo chuyên sâu cho các đề tài nghiên cứu tiếp theo.
  • Các Nhà Vật lý và Kỹ sư Thiết kế Lò phản ứng: Sử dụng các dữ liệu về phân bố thông lượng $\Phi(r,z,E)$, vi phân tiết diện sinh neutron và hệ số $k_{eff}$ để thiết kế chi tiết cấu trúc lõi lò, hệ thống che chắn sinh học và mạch tải nhiệt chì lỏng.
  • Các Viện Nghiên cứu R&D và Tập đoàn Năng lượng: Nắm bắt cơ sở khoa học để phát triển các dự án thử nghiệm máy gia tốc công suất lớn kết hợp module hạt nhân dưới tới hạn.
  • Cơ quan Quản lý An toàn Bức xạ và Hạt nhân: Có cơ sở lý thuyết định lượng để xây dựng quy chuẩn an toàn đối với các công nghệ điện hạt nhân thế hệ mới sử dụng nhiên liệu Thori và kim loại lỏng.

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 lý thuyết độc đáo nhất là việc mô hình hóa và xác lập bản chất vật lý của phản ứng phá vỡ hạt nhân $(p,n)$ trong môi trường bia chì lỏng thể tích lớn đóng vai trò kép (vừa làm bia sinh neutron vừa làm chất truyền nhiệt). Luận án đã mở rộng lý thuyết va chạm nội hạt nhân và phân rã thống kê (INC-SDM) khi áp dụng vào môi trường kim loại lỏng đồng nhất, chỉ ra rằng hiệu suất phát neutron không chỉ phụ thuộc vào tiết diện vi phân vi mô $\frac{d^2\sigma}{dEd\Omega}$ mà còn được tối ưu hóa nhờ chiều dài tương tác thực tế mở rộng dọc theo đường kính vùng chiếu xạ, khắc phục giới hạn suy giảm neutron trong các cấu hình bia rắn truyền thống.

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

So với nghiên cứu của S. Meigo và cộng sự (1999) [34] (chỉ khảo sát bia chì rắn ở năng lượng cố định $0,5\text{ GeV}$) và nghiên cứu của G. Bauer (2001) [36] (khảo sát hình học bia chì rắn dày $20\text{ cm}$ ở $2\text{ GeV}$), luận án của Trần Minh Tiến đã đổi mới vượt bậc về phương pháp luận khi:

  1. Mở rộng dải năng lượng khảo sát từ $0,25\text{ GeV}$ đến $3,0\text{ GeV}$ kết hợp khai thác thư viện dữ liệu năng lượng cao nhất JENDL-HE-2007 với định dạng ENDF-6.
  2. Tích hợp trực tiếp dữ liệu tương tác $(p,n)$ trên bia chì lỏng vào mô hình hình học toàn phần 3D của lõi lò TRIGA Mark II dưới tới hạn nạp nhiên liệu Thori đa dạng ($\text{Th}^{233}\text{UO}_2, \text{Th}^{235}\text{UO}_2, \text{Th}^{238}\text{UO}_2$), thay vì chỉ mô phỏng bia đơn lẻ tách rời khỏi lõi lò như các công trình tiền nhiệm.

3. Phát hiện bất ngờ hoặc có tính phản trực giác nhất từ dữ liệu nghiên cứu là gì?

Phát hiện bất ngờ nhất là mặc dù mật độ khối của chì lỏng ($\rho = 10,66\text{ g}\cdot\text{cm}^{-3}$) thấp hơn chì rắn ($\rho = 11,7\text{ g}\cdot\text{cm}^{-3}$) khoảng $8,9%$, nhưng thông lượng neutron sinh ra và hiệu suất phát neutron tổng cộng $Y_{n/p}$ trong lõi lò không hề suy giảm mà ngược lại còn có xu hướng phân bố đồng đều và bao quát hơn. Điều này được giải thích bởi hiệu ứng mở rộng thể tích tương tác: chùm proton thâm nhập sâu hơn vào bể chì lỏng mà không bị giới hạn bởi bề mặt phân cách cấu trúc cơ học của bia rắn, tạo ra vùng sinh neutron spallation trải dài, giúp tăng cường xác suất bắt neutron của các thanh nhiên liệu Thori xung quanh.

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 quy trình tái lập nghiên cứu vô cùng minh bạch và chuẩn xác:

  • Đặc tả hình học: Cung cấp chi tiết thông số kích thước cấu trúc lõi lò TRIGA Mark II, cấu tạo thanh nhiên liệu, bán kính ống dẫn chùm proton ($4\text{ cm}$).
  • Nguồn dữ liệu: Xác định cụ thể thư viện JENDL-HE-2007 của JAEA và định dạng ENDF.
  • Tham số chùm hạt: Năng lượng từ $250 - 3000\text{ MeV}$, dòng điện $25\text{ mA}$.
  • Hệ tọa độ phân tích: Lưới phân chia 18 khoảng góc khối, 32 nhóm năng lượng neutron, các tọa độ bán kính ($r = 2 - 26,25\text{ cm}$) và chiều cao ($z = 5 - 65\text{ cm}$) được số hóa cụ thể trong các biểu bảng và đồ thị của luận án.

5. Chương trình nghiên cứu 10 năm được vạch ra từ kết quả luận án bao gồm những gì?

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

  • Năm 1 - 3: Xây dựng nền tảng mô phỏng đa vật lý ghép nối Neutronics - CFD - Burnup thời gian thực trên các hệ máy tính hiệu năng cao (HPC).
  • Năm 4 - 6: Nghiên cứu vật liệu vỏ bọc chống ăn mòn chì lỏng ở nhiệt độ cao ($400 - 550^\circ\text{C}$) và phát triển công nghệ cảm biến đo đạc neutron trực tuyến trong môi trường kim loại lỏng.
  • Năm 7 - 8: Hợp tác quốc tế triển khai các phép đo thực nghiệm kiểm chứng chéo trên các tổ hợp máy gia tốc chùm proton kết hợp tổ hợp dưới tới hạn (như KUCA tại Nhật Bản hoặc các cơ sở tại Châu Âu).
  • Năm 9 - 10: Hoàn thiện thiết kế kỹ thuật chi tiết cho một lò phản ứng nghiên cứu dưới tới hạn quy mô pilot ($10 - 20\text{ MWth}$) sử dụng nhiên liệu Thori và chì lỏng tại Việt Nam.

Kết luận

  1. Luận án đã đề xuất và xây dựng thành công mô hình vật lý đột phá của hệ thống lò phản ứng hạt nhân dưới tới hạn (ADSR) sử dụng chì lỏng đồng thời làm bia bắn phá $(p,n)$ và chất tải nhiệt đối lưu, giải quyết triệt để bài toán suy giảm tuổi thọ và gián đoạn vận hành của bia rắn truyền thống.
  2. Thiết lập bộ dữ liệu toàn diện về vi phân bậc hai tiết diện sinh neutron $\frac{d^2\sigma}{dEd\Omega}$, hiệu suất phát neutron $Y_{n/p}$ và phân bố góc-năng lượng của chùm proton $0,25 - 3,0\text{ GeV}$ trên bia chì lỏng dựa trên thư viện JENDL-HE-2007.
  3. Định lượng hóa đầy đủ phổ phát xạ phóng xạ ($\alpha, \beta, \gamma, \bar{\nu}$) và động học các hạt nhân con trong quá trình phân rã và chuyển hóa của Thori trong môi trường chì lỏng.
  4. Mô phỏng không gian 3 chiều phân bố thông lượng neutron $\Phi(r,z,E)$ trong lõi lò TRIGA Mark II dưới tới hạn, chứng minh sự hình thành phổ neutron kép (nhiệt và nhanh) lý tưởng cho việc chuyển hóa $^{232}\text{Th}$ thành $^{233}\text{U}$ và kích hoạt phân hạch hiệu quả.
  5. Xác định chính xác quy luật biến thiên của hệ số nhân neutron hiệu dụng $k_{eff}$ trên các ma trận nhiên liệu hỗn hợp $\text{Th}^{233}\text{UO}_2$, $\text{Th}^{235}\text{UO}_2$, $\text{Th}^{238}\text{UO}_2$, khẳng định tính khả thi tuyệt đối của việc nạp nhiên liệu Thori hàm lượng cao mà vẫn duy trì biên độ an toàn dưới tới hạn hoàn hảo.
  6. Mở ra 3 hướng nghiên cứu mới mang tầm quốc tế: mô phỏng đa vật lý ghép nối nhiệt thủy lực chất lưu kim loại lỏng; công nghệ chu trình nhiên liệu Thori sạch không sinh Plutoni vũ khí; và giải pháp thiêu hủy chất thải phóng xạ actinide hiếm bảo vệ môi trường toàn cầu. Công trình là một dấu mốc học thuật quan trọng, khẳng định vị thế và năng lực hội nhập đỉnh cao của ngành vật lý hạt nhân Việt Nam trên bản đồ năng lượng nguyên tử thế giới.