Tổng quan về luận án

Sự phát triển của nền kinh tế toàn cầu và áp lực biến đổi khí hậu đang đặt ngành năng lượng trước bài toán chuyển dịch cơ cấu cấp bách. Mặc dù năng lượng hạt nhân đóng vai trò trụ cột trong việc cung cấp nguồn điện nền phát thải carbon thấp, các công nghệ lò phản ứng hạt nhân truyền thống (chủ yếu là lò phản ứng nước nhẹ - LWR sử dụng nhiên liệu urani) đang đối mặt với ba rào cản mang tính sống còn: nguy cơ cạn kiệt tài nguyên urani tự nhiên, bài toán an toàn nghiêm ngặt sau các thảm họa Chernobyl hay Fukushima, và thách thức xử lý hàng trăm nghìn tấn chất thải phóng xạ chu kỳ bán rã dài (đặc biệt là các actini hiếm - minor actinides).

Trước bối cảnh đó, công nghệ 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 với chu trình nhiên liệu Thori ($^{232}\text{Th}$) nổi lên như một giải pháp đột phá. Hoạt động của ADSR dựa trên cơ chế ghép nối máy gia tốc hạt năng lượng cao với một tâm lò luôn duy trì ở trạng thái dưới tới hạn ($k_{eff} < 1$). Chùm proton năng lượng cao (vài trăm MeV đến vài GeV) bắn phá vào hạt nhân bia kim loại nặng, tạo ra phản ứng phá vỡ hạt nhân (spallation) sinh ra nguồn neutron bù ngoại lai dồi dào, duy trì phản ứng phân hạch mà không thể tự xảy ra phản ứng dây chuyền mất kiểm soát.

Khoảng trống nghiên cứu then chốt mà luận án tiến sĩ của NCS. Trần Minh Tiến (chuyên ngành Vật lý nguyên tử và hạt nhân, Viện Năng lượng nguyên tử Việt Nam, người hướng dẫn: PGS. Trần Quốc Dũng và PGS. Nguyễn Mộng Giao) giải quyết là: Hầu hết các công trình nghiên cứu ADSR quốc tế trước đây chủ yếu tập trung vào mô hình bia tương tác rắn (như chì rắn, vonfram, vàng) hoặc hợp kim chì - bismuth ($\text{Pb-Bi}$). Việc sử dụng bia rắn bộc lộ nhược điểm cơ học nghiêm trọng do ứng suất bức xạ cao, phá hủy cấu trúc tinh thể, bắt buộc phải dừng hoạt động của lò định kỳ để thay bia, gây gián đoạn vận hành và rủi ro an toàn phóng xạ. Trong khi đó, việc sử dụng hợp kim chì - bismuth lại sinh ra đồng vị phóng xạ cực độc $\text{Po-210}$ dễ bay hơi. Nghiên cứu của NCS. Trần Minh Tiến đề xuất mô hình mang tính tiên phong: Sử dụng trực tiếp chì lỏng đồng thời làm bia tương tác sinh neutron spallation và chất tải nhiệt làm mát lõi lò, kết hợp đánh giá toàn diện khả năng tích hợp nhiên liệu Thori tự nhiên.

Hệ thống câu hỏi nghiên cứu (Research Questions) và giả thuyết khoa học (Hypotheses) được xác lập chặt chẽ:

  • RQ1: Cơ chế phát xạ neutron từ tương tác $(p,n)$ trên bia chì lỏng thay đổi như thế nào theo năng lượng chùm proton tới ($250\text{ MeV} - 3000\text{ MeV}$) và phân bố góc khối ($0^\circ - 180^\circ$)?
    • H1: Mặc dù mật độ khối của chì lỏng ($10{,}66\text{ g/cm}^3$) thấp hơn chì rắn ($11{,}7\text{ g/cm}^3$), thể tích tương tác mở rộng toàn bộ đường kính lò sẽ bù đắp và tối ưu hóa hiệu suất sinh neutron $Y_{n/p}$, đồng thời loại bỏ nhu cầu dừng lò thay thế bia.
  • RQ2: Các đặc trưng phân rã phóng xạ của nhiên liệu Thori trong môi trường chì lỏng và phân bố phổ thông lượng neutron trong lõi lò diễn ra theo quy luật vật lý nào?
    • H2: Nhiên liệu Thori tạo ra phổ phân rã và phân bố thông lượng neutron dọc theo bán kính và chiều cao lõi lò đạt độ đồng đều cao hơn, giảm thiểu các điểm nóng cục bộ so với nhiên liệu urani truyền thống.
  • RQ3: Tỷ lệ pha trộn tối ưu giữa thori và các đồng vị urani ($\text{Th-}^{233}\text{U}$, $\text{Th-}^{235}\text{U}$, $\text{Th-}^{238}\text{U}$) trong cấu trúc hình học TRIGA Mark II là bao nhiêu để duy trì hệ số nhân neutron hiệu dụng $k_{eff}$ trong giới hạn dưới tới hạn an toàn?
    • H3: Hỗn hợp nhiên liệu $\text{Th-}^{233}\text{UO}2$ với hàm lượng thori từ $60%$ đến $80%$ vẫn đảm bảo hệ số $k{eff}$ vận hành tối ưu ($0{,}90 \le k_{eff} < 1{,}00$) khi kết hợp với nguồn neutron bù từ bia chì lỏng.

Khung lý thuyết của công trình tích hợp sâu sắc giữa lý thuyết vật lý hạt nhân năng lượng cao, lý thuyết phản ứng phá vỡ hạt nhân (nuclear spallation physics), lý thuyết chuyển hóa nhiên liệu Thori - Urani ($^{232}\text{Th}(n,\gamma)^{233}\text{Th} \xrightarrow{\beta^-} ,^{233}\text{Pa} \xrightarrow{\beta^-} ,^{233}\text{U}$), và lý thuyết vận chuyển hạt Monte Carlo. Công trình mang lại tác động định lượng đột phá: xác lập bộ dữ liệu vi phân bậc hai tiết diện sinh neutron $\frac{d^2\sigma}{dEd\Omega}$, phổ năng lượng hạt phát xạ $\alpha, \beta, \gamma, \bar{\nu}$, và phân bố thông lượng không gian 3 chiều bên trong lò phản ứng TRIGA Mark II cải tiến.

Phạm vi nghiên cứu bao quát dải năng lượng proton từ $250\text{ MeV}$ đến $3\text{ GeV}$ với dòng chùm $25\text{ mA}$, bán kính chùm $4\text{ cm}$, khai thác thư viện dữ liệu hạt nhân năng lượng cao JENDL-HE-2007 (chia 19 phân khoảng góc và 32 nhóm năng lượng), kết hợp với các chương trình mô phỏng chuẩn quốc tế MCNPX, MCNP5 và GEANT4.

Literature Review và Positioning

Lịch sử phát triển của hệ thống lò phản ứng điều khiển bằng máy gia tốc bắt nguồn từ những quan sát sơ khai của Ernest Rutherford (1919) và việc Ernest Lawrence phát minh máy gia tốc cyclotron (1930s). Đến năm 1940, Lawrence và Semenov đã độc lập đề xuất việc sử dụng máy gia tốc hạt làm nguồn phát sinh neutron ngoại lai. Khái niệm ADSR hiện đại thực sự bùng nổ khi nhóm nghiên cứu của C. Bowman (1992, 2001) và người đoạt giải Nobel Carlo Rubbia (1995) công bố thiết kế "Bộ khuếch đại năng lượng" (Energy Amplifier), mở ra hướng đi kết hợp máy gia tốc chùm proton với lò phản ứng dưới tới hạn để thiêu hủy chất thải hạt nhân và sản xuất điện năng.

Tổng quan các luồng nghiên cứu quốc tế chỉ ra ba trục chính:

  1. Nghiên cứu tương tác spallation và phân bố neutron trên bia nặng: S. Meigo và cộng sự (1999) đã sử dụng MCNP4A tính toán phân bố thông lượng neutron sinh ra từ bia chì dày với chùm proton $0{,}5\text{ GeV}$, xác định phổ neutron trải dài từ dải nhiệt đến dải nhanh. A. Goldenbaum và cộng sự (2000) khảo sát hiệu suất neutron trên các bia dày vonfram, thủy ngân, chì ở dải năng lượng $0{,}5 - 2{,}5\text{ GeV}$, đưa ra kết luận về chiều dài giới hạn của bia (chì là $55\text{ cm}$, vonfram là $30\text{ cm}$). G. Bauer (2001) làm sáng tỏ phân bố góc và năng lượng của neutron phân hạch ở các góc $30^\circ, 90^\circ, 150^\circ$ khi bắn phá proton $2\text{ GeV}$ lên bia chì dày $20\text{ cm}$. A. Krasa và cộng sự (2008), David Sangcheol Lee (2018) tiếp tục làm sâu sắc đặc trưng phổ spallation và phổ phân hạch trên các bia $\text{Pb}$ và $\text{Pb-Bi}$.
  2. Nghiên cứu chu trình nhiên liệu Thori: Từ thí nghiệm lò phản ứng muối nóng chảy MSRE tại Oak Ridge (Alvin Weinberg, Ralph Moir, 1969) đến chu trình nhiên liệu Thori của Radkowsky và lò phản ứng nước nặng CANDU, tiềm năng của $^{232}\text{Th}$ đã được thừa nhận rộng rãi nhờ trữ lượng tự nhiên dồi dào gấp 4 lần Urani (ước tính khoảng $1{,}2 \times 10^{14}$ tấn trên vỏ Trái Đất, chiếm khoảng 6 ppm trong đất). Các đề xuất ADTR (Accelerator Driven Thorium Reactor) của Victoria B. Ashley và cộng sự, cùng mô hình ADSR muối nóng chảy tại Nhật Bản ($60\text{NaCl} - 34\text{ThCl}_4 - 6^{233}\text{UCl}_4$) đã khẳng định tính khả thi của việc sử dụng thori trong phổ neutron nhanh và nhiệt.
  3. Các nghiên cứu tiền phong tại Việt Nam: Nhóm nghiên cứu của Nguyễn Mộng Giao, Châu Văn Tạo và Nguyễn Thị Ái Thu (2002–2015) đã phát triển mô hình màn chắn trên bia dày để hiệu chỉnh sai số lý thuyết - thực nghiệm trong tương tác $(p,n)$ trên các bia $\text{U, Pb, W, Au}$ ở năng lượng $0{,}5 - 1{,}5\text{ GeV}$. Tiếp đó, Vũ Thanh Mai và các cộng sự (Đại học Kyoto) đã thực hiện hàng loạt thực nghiệm trên tổ hợp tới hạn KUCA với chùm proton $100\text{ MeV}$ trên bia kim loại nặng kết hợp nhiên liệu hỗn hợp $^{235}\text{U}$ và $^{232}\text{Th}$.
Nghiên cứu / Dự án Quốc gia / Tổ chức Loại bia tương tác Hệ nhiên liệu Năng lượng hạt & Công suất Hạn chế / Điểm nghẽn kỹ thuật
MYRRHA (SCK.CEN) Bỉ / Châu Âu $\text{Pb-Bi}$ lỏng/rắn MOX (Pu-U oxit) $600\text{ MeV}$ ($4\text{ mA}$), $70\text{ MWth}$ Nguy cơ tạo thành $\text{Po-210}$ độc hại; ăn mòn vật liệu vỏ ở nhiệt độ cao
C-ADS (CAS) Trung Quốc Bia hạt rắn / $\text{Pb-Bi}$ $\text{UO}_2$ / Nhiên liệu hạt $100\text{ MWth} \rightarrow 1\text{ GW}$ Phức tạp trong cơ cấu tản nhiệt bia rắn và độ bền chùm gia tốc
KIPT Ukraine / ANL (Mỹ) Bia kim loại LEU ($\text{UO}_2$ làm giàu thấp) Linac electron/proton, làm mát bằng $\text{H}_2\text{O}$ Hiệu suất sinh neutron hạn chế, sử dụng chu trình Urani truyền thống
Luận án NCS. Trần Minh Tiến Việt Nam (VAEC) Chì nguyên chất lỏng ($\text{Pb}$ lỏng) Hỗn hợp $\text{ThO}_2 - \text{UO}_2$ $250\text{ MeV} - 3\text{ GeV}$, TRIGA Mark II Giải quyết triệt để sự suy thoái bia rắn; loại bỏ $\text{Po-210}$; tận dụng tối đa chu trình $^{232}\text{Th}$

Các tranh luận khoa học lớn trong y văn xoay quanh hai trường phái:

  • Trường phái bia rắn vs. Bia lỏng: Các nghiên cứu của Nifenecker và Letourneau ưu tiên bia rắn vì mật độ hạt nhân cao ($11{,}7\text{ g/cm}^3$), cho hiệu suất sinh neutron tức thời trên một đơn vị thể tích lớn hơn. Ngược lại, quan điểm của Rubbia và công trình của NCS. Trần Minh Tiến chứng minh rằng bia lỏng tự đối lưu loại bỏ hoàn toàn hiện tượng phá hủy cấu trúc do sốc nhiệt và chiếu xạ, kéo dài tuổi thọ vận hành liên tục của lò phản ứng mà không suy giảm tổng thông lượng neutron.
  • Trường phái phổ nhanh vs. Phổ nhiệt cho Thori: Trong khi các lò LFR thế hệ IV thiên về phổ nhanh để thiêu hủy actinide, thì việc ứng dụng cấu trúc TRIGA Mark II với chất làm mát chì lỏng trong luận án này tạo ra một phổ neutron hỗn hợp tối ưu: vừa đủ cứng ở vùng gần bia để kích hoạt phản ứng phân hạch và spallation, vừa có thành phần neutron nhiệt hóa để tối đa hóa tiết diện bắt neutron chuyển đổi $^{232}\text{Th}$ thành $^{233}\text{U}$.

Đó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 các khung lý thuyết vật lý hạt nhân then chốt:

graph TD
    A["Chùm Proton Năng Lượng Cao (250 MeV - 3 GeV)"] --> B["Bia Chì Lỏng (Đồng thời là Chất Tải Nhiệt)"]
    B --> C["Tương tác Spallation (p,n) sinh Neutron Ngoại Lai"]
    C --> D["Phổ Neutron Hỗn Hợp (Nhanh & Nhiệt Hóa)"]
    D --> E["Bắt Neutron: 232Th + n -> 233Th -> 233Pa -> 233U"]
    D --> F["Phân Hạch Duy Trì Dưới Tới Hạn (keff = 0.90 - 0.98)"]
    E --> F
    F --> G["Truyền Nhiệt Đối Lưu Liên Tục qua Chì Lỏng (327°C)"]
    G --> H["An Toàn Nội Tại: Tắt Máy Gia Tốc -> Lò Tự Ngừng Hoạt Động"]
  1. Mở rộng lý thuyết tương tác spallation trên môi trường kim loại lỏng liên tục: Trước đây, các mô hình phân rã thống kê (Statistical Decay Model - SDM) và mô hình động lực học phân tử lượng tử (Quantum Molecular Dynamics - QMD) chủ yếu được kiểm chứng trên các bia rắn tĩnh. Luận án đã chứng minh tính tương thích và mở rộng phạm vi áp dụng của các mô hình này lên môi trường chì lỏng đồng nhất có mật độ $10{,}66\text{ g/cm}^3$, xác lập mối quan hệ hàm số giữa năng lượng chùm proton $E_p$, góc phát xạ $\theta$, và tiết diện vi phân sinh neutron.
  2. Hệ thống hóa các định đề lý thuyết (Theoretical Propositions):
    • Định đề P1: Hiệu suất phát neutron spallation $Y_{n/p}$ trên bia chì lỏng tỉ lệ thuận phi tuyến với năng lượng proton tới $E_p$ trong dải $0{,}25 - 3{,}0\text{ GeV}$, và đạt trạng thái tăng tuyến tính tiệm cận khi $E_p > 1{,}5\text{ GeV}$.
    • Định đề P2: Vi phân bậc hai tiết diện sinh neutron $\frac{d^2\sigma}{dEd\Omega}$ trong môi trường chì lỏng có tính bất đối xứng góc rõ rệt ở dải năng lượng neutron cao ($E_n > 20\text{ MeV}$ tập trung ở góc phát xạ nhỏ $\theta < 30^\circ$) và chuyển dịch sang phân bố đẳng hướng ở dải neutron nhiệt ($E_n < 1\text{ eV}$).
    • Định đề P3: Năng lượng phân rã phóng xạ của hệ nhiên liệu thori trong chì lỏng phân bố ưu thế qua kênh phân rã alpha ($\alpha$) của các hạt nhân con trong dãy phóng xạ, tạo ra nguồn nhiệt bổ sung ổn định mà không gây nhiễu loạn thông lượng neutron điều khiển.
    • Định đề P4: Hệ số nhân neutron hiệu dụng $k_{eff}$ của lõi ADSR tỷ lệ nghịch với tỷ lệ thori hóa trong hỗn hợp nhiên liệu oxit, nhưng độ dốc suy giảm của $\text{Th-}^{233}\text{U}$ thấp hơn đáng kể so với $\text{Th-}^{235}\text{U}$ và $\text{Th-}^{238}\text{U}$, cho phép duy trì trạng thái dưới tới hạn sâu nhưng ổn định.
  3. Chuyển dịch mô thức (Paradigm Shift): Nghiên cứu tạo ra bước chuyển dịch từ tư duy "kiểm soát độ phản ứng dựa trên thanh điều khiển cơ học trong lò tới hạn" sang "kiểm soát công suất tức thời thông qua cường độ dòng máy gia tốc trong hệ thống dưới tới hạn an toàn nội tại" (inherently safe subcritical systems).

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

Khung phân tích của luận án tích hợp đồng thời ba trụ cột lý thuyết:

  • Lý thuyết vận chuyển hạt Boltzmann đa nhóm năng lượng: Mô tả hành vi chuyển động, tán xạ đàn hồi, tán xạ không đàn hồi, hấp thụ và phân hạch của neutron trong môi trường không đồng nhất.
  • Lý thuyết phản ứng hạt nhân năng lượng cao: Dựa trên mô hình tầng va chạm nội hạt nhân (Intranuclear Cascade - INC) kết hợp mô hình bốc hơi hạt nhân (Evaporation Model).
  • Lý thuyết nhiệt thủy động học chất lỏng kim loại nặng: Phân tích cơ chế trao đổi nhiệt đối lưu của chì lỏng ở nhiệt độ vận hành trên điểm nóng chảy $327^\circ\text{C}$.

Phương pháp tiếp cận mới mẻ này xác lập các điều kiện biên tường minh: năng lượng chùm hạt proton $250\text{ MeV} \le E_p \le 3000\text{ MeV}$, cường độ dòng $I_p = 25\text{ mA}$, bán kính chùm $r_{beam} = 4\text{ cm}$, vùng hoạt lò phản ứng TRIGA Mark II với cấu trúc thanh nhiên liệu chuẩn, và môi trường làm mát chì lỏng bao bọc toàn bộ thể tích lõi.

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 tính toán (computational positivism) kết hợp phương pháp mô phỏng ngẫu nhiên Monte Carlo và phân tích số trị tất định. Thiết kế nghiên cứu đa tầng (multi-level design) được cấu trúc thành ba cấp độ tương tác:

  • Cấp độ vi mô: Tương tác $(p,n)$ giữa proton đơn lẻ và hạt nhân chì ($\text{Pb}$), xác định tiết diện vi phân sinh neutron.
  • Cấp độ trung mô: Quá trình phân rã phóng xạ của Thori và các hạt nhân con trong môi trường chì lỏng; phổ năng lượng của các hạt $\alpha, \beta, \gamma$ và neutrino ($\bar{\nu}$).
  • Cấp độ vĩ mô: Sự vận chuyển neutron, phân bố thông lượng $\Phi(r, z, E)$ và hệ số nhân neutron hiệu dụng $k_{eff}$ trong toàn bộ thể tích hình học lò phản ứng TRIGA Mark II.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu được thực hiện qua các giao thức nghiêm ngặt:

graph LR
    subgraph Giai đoạn 1: Dữ liệu hạt nhân
        A["Thư viện JENDL-HE-2007 (ENDF Format)"] --> B["Trích xuất Tiết diện: 19 Phân khoảng góc & 32 Nhóm năng lượng"]
    end
    subgraph Giai đoạn 2: Tính toán Spallation
        B --> C["Chương trình MCNPX / GEANT4"]
        C --> D["Tính d²σ/dEdΩ, Y_n/p, Phổ Năng lượng Neutron"]
    end
    subgraph Giai đoạn 3: Mô hình Lõi lò TRIGA Mark II
        D --> E["Mô phỏng 3D: Bia Chì Lỏng + Chất làm mát Pb"]
        E --> F["Nạp Nhiên liệu: Th-233UO2 / Th-235UO2 / Th-238UO2"]
    end
    subgraph Giai đoạn 4: Đánh giá Tham số Neutron
        F --> G["Phân bố Thông lượng Φ(r, z, E)"]
        F --> H["Tính toán Hệ số keff & Kiểm tra An toàn Dưới tới hạn"]
    end

  1. Giao thức trích xuất dữ liệu: Sử dụng thư viện dữ liệu hạt nhân năng lượng cao JENDL-HE-2007 định dạng ENDF của Cơ quan Năng lượng Nguyên tử Nhật Bản (JAEA), chứa dữ liệu tương tác hạt nhân của 406 đồng vị lên đến năng lượng $3\text{ GeV}$.

  2. Thuật toán xử lý ma trận: Vi phân bậc hai của tiết diện sinh neutron theo năng lượng và góc khối được tính toán thông qua công thức giải tích chuẩn: $$\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)$$ Từ đó, tiết diện vi phân theo góc khối được tích phân: $$\frac{d\sigma (\mu, E_p, E_n)}{d\Omega} = \int \frac{d^2 \sigma (\mu, E_p, E_n)}{dE d\Omega} dE = \frac{\sigma(E_p) \cdot \gamma(E_p)}{4\pi} \sum_{i=1}^{32} f(\mu, E_p, E_i) \Delta E_i$$ Trong đó:

    • $E_p$: Năng lượng của proton tới ($\text{eV}$).
    • $E_n$: Năng lượng của neutron phát ra ($\text{eV}$).
    • $\sigma(E_p)$: Tiết diện phản ứng toàn phần ($\text{barn}$).
    • $\gamma(E_p)$: Hiệu suất phát neutron toàn phần.
    • $f(\mu, E_p, E_n)$: Hàm phân bố mật độ xác suất chuẩn hóa (ma trận $32 \times 19$).
    • $\mu = \cos\theta$: Cosine của góc phát xạ neutron, chia thành 19 điểm nút ($\mu_1, \mu_2, \dots, \mu_{19}$) tương ứng 18 khoảng góc từ $0^\circ$ đến $180^\circ$.
    • $E_n$: Chia thành 32 nhóm năng lượng ($E_1, E_2, \dots, E_{32}$).

    Tổng số neutron sinh ra $N_n$ trên bia được xác định qua biểu thức: $$N_n = N_p \sum_{j=1}^{32} N \frac{d\sigma (\mu, E_p, E_{n,j})}{dE} d$$ với $N_p$ là số hạt proton tới, $N$ là mật độ nguyên tử bia chì lỏng, và $d$ là quãng đường tương tác.

  3. Kiểm định độ tin cậy và giá trị (Validity & Reliability): Kết quả tính toán vi phân bậc hai tiết diện sinh neutron $\frac{d^2\sigma}{dEd\Omega}$ trên bia chì lỏng được đối chiếu chéo (cross-validation) trực tiếp với kết quả thực nghiệm trên bia chì rắn của nhóm tác giả X. Ledoux và cộng sự ở các góc $45^\circ - 55^\circ$ và $60^\circ$ tại các mức năng lượng $250\text{ MeV}, 500\text{ MeV}, 1000\text{ MeV}, 2000\text{ MeV}$, cho thấy sự trùng khớp chặt chẽ về dạng phổ và độ lớn, khẳng định độ tin cậy tuyệt đối của mô hình mô phỏng.

Data và phân tích

  • Thông số cấu trúc hình học: Lõi lò TRIGA Mark II được mô phỏng chi tiết trong MCNPX với các bó thanh nhiên liệu hình trụ, bố trí theo các vành tròn đồng tâm, bao quanh là vành phản xạ neutron và bể chứa chì lỏng.
  • Môi trường vật liệu: Chì lỏng nguyên chất đóng vai trò kép: bia spallation trung tâm và môi chất truyền nhiệt bao quanh các thanh nhiên liệu với khối lượng riêng $10{,}66\text{ g/cm}^3$ tại trạng thái nóng chảy ($> 327^\circ\text{C}$).
  • Hệ nhiên liệu khảo sát: Khảo sát ba nhóm hỗn hợp nhiên liệu gốm oxit với các dải tỷ lệ phần trăm Thori biến thiên từ $0%$ đến $100%$ (bước nhảy $20%$):
    1. $\text{ThO}_2 - ,^{233}\text{UO}_2$
    2. $\text{ThO}_2 - ,^{235}\text{UO}_2$
    3. $\text{ThO}_2 - ,^{238}\text{UO}_2$ kết hợp đối chứng với nhiên liệu chuẩn $\text{UZrH}$ làm mát bằng nước nhẹ và làm mát bằng chì lỏng.

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

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

  1. Đặc trưng hiệu suất phát neutron $Y_{n/p}$ vượt trội trên bia chì lỏng: Kết quả tính toán chứng minh hiệu suất phát sinh neutron $Y_{n/p}$ tăng đơn điệu theo năng lượng chùm proton tới. Ở mức năng lượng $E_p = 500\text{ MeV}$, hiệu suất đạt xấp xỉ $10 - 12\text{ neutron/proton}$; khi năng lượng tăng lên $E_p = 1500\text{ MeV}$ ($1{,}5\text{ GeV}$), hiệu suất đạt tới $27 - 30\text{ neutron/proton}$; và tại $E_p = 3000\text{ MeV}$ ($3{,}0\text{ GeV}$), hiệu suất vượt ngưỡng $55\text{ neutron/proton}$. Mặc dù mật độ chì lỏng giảm khoảng $8{,}89%$ so với chì rắn, nhưng do chùm proton tương tác dọc toàn bộ trục chiều dài lõi lò chì lỏng, tổng số neutron sinh ra thực tế không hề suy giảm mà phân bố đồng đều hơn vào không gian buồng phản ứng.
  2. Quy luật phân bố góc và năng lượng của neutron spallation: Phổ vi phân bậc hai tiết diện sinh neutron $\frac{d^2\sigma}{dEd\Omega}$ chỉ ra rằng ở các góc phát xạ hẹp phía trước ($\theta < 30^\circ$), thông lượng neutron mang năng lượng cao ($E_n > 50\text{ MeV}$) chiếm tỷ trọng áp đảo. Tuy nhiên, ở các góc tán xạ rộng ($\theta > 60^\circ$ đến $150^\circ$), phổ năng lượng neutron dịch chuyển mạnh về dải năng lượng dưới $20\text{ MeV}$ và dải nhiệt ($E_n < 1\text{ eV}$). Đây là điều kiện lý tưởng để kích hoạt phản ứng bắt neutron của $^{232}\text{Th}$ ở các vùng nhiên liệu ngoại vi.
  3. Bản chất phổ phóng xạ tự nhiên của Thori trong môi trường chì lỏng: Tính toán chi tiết phổ năng lượng phân rã cho thấy các tia alpha ($\alpha$) sinh ra từ chuỗi phân rã Thori mang năng lượng tập trung trong dải $4{,}0 - 8{,}8\text{ MeV}$, các tia beta ($\beta^-$) có phổ liên tục đạt đỉnh cực đại ở $1{,}5 - 2{,}2\text{ MeV}$, và tia gamma ($\gamma$) có các đỉnh năng lượng đặc trưng tại $0{,}238\text{ MeV}$ (từ $^{212}\text{Pb}$) và $2{,}614\text{ MeV}$ (từ $^{208}\text{Tl}$). Toàn bộ năng lượng bức xạ này bị hấp thụ hoàn toàn trong khối chì lỏng xung quanh, chuyển hóa thành nhiệt năng hữu ích và được lớp chì lỏng che chắn an toàn sinh học tự nhiên.
  4. Tái lập cấu trúc phân bố thông lượng neutron $\Phi(r, z)$ trong lõi TRIGA Mark II: Khi thay thế chất làm mát nước nhẹ bằng chì lỏng và nạp nhiên liệu Thori, thông lượng neutron cực đại $\Phi_{max}$ dọc theo chiều cao lõi ($z$) đạt giá trị đỉnh tại vị trí trung tâm lò ($z = 30 - 35\text{ cm}$) và phân bố đối xứng dạng cosine cụt. Theo phương bán kính ($r$), thông lượng neutron giảm dần từ tâm bia chì lỏng ($r = 0 - 4\text{ cm}$) ra vành phản xạ ($r = 24 - 26\text{ cm}$), nhưng độ dốc giảm thoải hơn nhiều so với lò làm mát bằng nước nhẹ, giúp độ sâu cháy (burnup) của các thanh nhiên liệu thori ở các vành ngoài diễn ra đồng đều.
  5. Hành vi hệ số nhân neutron hiệu dụng $k_{eff}$ theo tỷ lệ cấp phối Thori:
Tỷ lệ Thori (% khối lượng) $k_{eff}$ với $\text{Th-}^{233}\text{UO}_2$ $k_{eff}$ với $\text{Th-}^{235}\text{UO}_2$ $k_{eff}$ với $\text{Th-}^{238}\text{UO}_2$ Trạng thái vật lý
0% (100% U) $1{,}185 \pm 0{,}002$ $1{,}120 \pm 0{,}002$ $0{,}320 \pm 0{,}001$ Trên tới hạn (với $^{233}\text{U}, ,^{235}\text{U}$)
20% Th $1{,}112 \pm 0{,}002$ $1{,}045 \pm 0{,}002$ $0{,}285 \pm 0{,}001$ Trên tới hạn
40% Th $1{,}038 \pm 0{,}001$ $0{,}972 \pm 0{,}001$ $0{,}248 \pm 0{,}001$ Cận tới hạn / Dưới tới hạn
60% Th $0{,}965 \pm 0{,}001$ $0{,}895 \pm 0{,}001$ $0{,}210 \pm 0{,}001$ Dưới tới hạn tối ưu cho ADSR
80% Th $0{,}882 \pm 0{,}001$ $0{,}810 \pm 0{,}001$ $0{,}175 \pm 0{,}001$ Dưới tới hạn an toàn cao
100% Th $0{,}125 \pm 0{,}001$ $0{,}125 \pm 0{,}001$ $0{,}125 \pm 0{,}001$ Không thể duy trì nếu thiếu U phân hạch

Kết quả khẳng định: Hỗn hợp $\text{Th-}^{233}\text{UO}_2$ ở tỷ lệ thori $60% - 80%$ và $\text{Th-}^{235}\text{UO}2$ ở tỷ lệ thori $40% - 60%$ là cấu hình vàng, duy trì hệ số $k{eff}$ trong khoảng $0{,}90 - 0{,}98$, bảo đảm lò phản ứng hoạt động hoàn toàn dưới tới hạn, triệt tiêu mọi khả năng xảy ra sự cố quá nhiệt mất kiểm soát.

Implications đa chiều

  • Ý nghĩa lý thuyết: Xác lập mô hình toán - lý hoàn chỉnh mô tả tương tác spallation và động học neutron trong hệ thống lai kim loại lỏng - thori, bổ sung dữ liệu thực nghiệm mô phỏng chuẩn xác cho cơ sở dữ liệu vật lý hạt nhân năng lượng cao.
  • Đổi mới phương pháp luận: Thiết lập quy trình mô phỏng liên hoàn từ thư viện hạt nhân JENDL-HE qua mã MCNPX để giải quyết đồng thời bài toán nguồn spallation và bài toán tới hạn của tâm lò.
  • Ứng dụng thực tiễn trong công nghệ hạt nhân: Cung cấp thiết kế cơ sở cho các dự án phát triển lò phản ứng nghiên cứu và lò năng lượng thế hệ mới, giải quyết triệt để vấn đề thay bia cơ học và loại trừ nguy cơ tích tụ chất độc phóng xạ $\text{Po-210}$.
  • Khuyến nghị chính sách và giải trừ vũ khí hạt nhân: Chu trình nhiên liệu Thori - ADSR không tạo ra đồng vị $\text{Pu-239}$ cấp độ vũ khí (weapon-grade plutonium), loại bỏ nguy cơ phổ biến vũ khí hạt nhân, tạo tiền đề vững chắc cho việc phê chuẩn các dự án năng lượng nguyên tử hòa bình.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn chỉ ra 4 giới hạn khoa học mang tính khách quan:

  1. Giới hạn về động học tích tụ Protactinium ($^{233}\text{Pa}$): Trong chu trình Thori, đồng vị trung gian $^{233}\text{Pa}$ có chu kỳ bán rã tương đối dài ($T_{1/2} = 27\text{ ngày}$) và có tiết diện bắt neutron ký sinh đáng kể. Mô hình mô phỏng tĩnh của luận án chưa tính toán chi tiết sự suy giảm hiệu suất chuyển đổi do hiệu ứng bắt neutron của $^{233}\text{Pa}$ theo thời gian vận hành dài hạn (burnup time-dependent kinetics).
  2. Giới hạn tương tác nhiệt - thủy động lực học ghép nối (Coupled Neutronic-CFD): Luận án giả định mật độ chì lỏng là đại lượng không đổi ($10{,}66\text{ g/cm}^3$) trên toàn bộ lõi lò. Trong điều kiện vận hành thực tế, sự chênh lệch nhiệt độ giữa tâm lò và ngoại vi sẽ gây ra gradient mật độ và dòng chảy đối lưu nhiệt phức tạp, có thể làm thay đổi nhẹ phân bố thông lượng neutron cục bộ.
  3. Thách thức ăn mòn vật liệu ở nhiệt độ cao: Chì lỏng ở nhiệt độ $> 400^\circ\text{C}$ có tính ăn mòn hóa học rất mạnh đối với các loại thép hợp kim và vỏ thanh nhiên liệu. Luận án chỉ giới hạn ở việc mô phỏng vật lý neutron mà chưa khảo sát cơ chế suy thoái cơ - hóa của vật liệu cấu trúc.
  4. Hiện tượng ứng suất màng cửa sổ chùm proton (Beam Window Stress): Mô phỏng giả định chùm proton đi trực tiếp vào khối chì lỏng mà chưa mô hình hóa chi tiết sự suy giảm thông lượng hạt qua màng ngăn chân không - kim loại lỏng.

Chương trình nghiên cứu tiếp theo (Future Research Agenda) được đề xuất gồm 5 định hướng cụ thể:

  • Hướng 1: Ghép nối đa vật lý (Multiphysics Coupling) giữa bộ mã neutron Monte Carlo (MCNPX/Serpent) với phần mềm thủy động lực học tính toán (CFD OpenFOAM/ANSYS Fluent) để khảo sát trường nhiệt - thủy lực thực tế của chì lỏng.
  • Hướng 2: Nghiên cứu động thái cháy nhiên liệu chi tiết (Long-term Burnup and Depletion Analysis) theo chu kỳ 1000 - 3000 ngày vận hành, theo dõi chính xác nồng độ các đồng vị $^{233}\text{Pa}, ,^{233}\text{U}$, và các sản phẩm phân hạch ngộ độc neutron như $^{135}\text{Xe}, ,^{149}\text{Sm}$.
  • Hướng 3: Khảo sát công nghệ vật liệu vỏ bọc chống ăn mòn chì lỏng (như hợp kim FeCrAl phủ lớp oxit $\text{Al}_2\text{O}_3$ tự phục hồi) trong điều kiện chiếu xạ neutron cường độ cao.
  • Hướng 4: Mở rộng bài toán thiêu hủy actini hiếm ($\text{Np, Am, Cm}$) bằng cách đưa chất thải phóng xạ vào thành phần nhiên liệu Thori - Urani trong ADSR chì lỏng.
  • Hướng 5: Thiết kế tối ưu hóa hình học cửa sổ chùm tia (beam windowless concept) để chùm proton tiếp xúc trực tiếp bề mặt tự do của xoáy chì lỏng.

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

  • Tác động học thuật: Công trình đóng góp bộ tham số vật lý và phương pháp luận tính toán chuẩn xác cho cộng đồng vật lý hạt nhân quốc tế. Bộ dữ liệu vi phân bậc hai tiết diện sinh neutron $\frac{d^2\sigma}{dEd\Omega}$ trên chì lỏng là nguồn tham chiếu quý giá cho các nghiên cứu mô phỏng ADSR toàn cầu, với tiềm năng trích dẫn cao trong các tạp chí chuyên ngành hàng đầu như Annals of Nuclear Energy, Nuclear Engineering and Design, Progress in Nuclear Energy.
  • Tác động chuyển dịch công nghiệp: Mở ra lộ trình kỹ thuật khả thi cho ngành công nghiệp chế tạo lò phản ứng hạt nhân thế hệ IV (Gen-IV Lead-cooled Fast Reactors - LFR và ADSR). Giải pháp loại bỏ bia rắn giúp các nhà vận hành nhà máy điện hạt nhân tương lai tiết kiệm hàng trăm triệu USD chi phí bảo dưỡng định kỳ và xử lý bia thải nhiễm xạ.
  • Ảnh hưởng chính sách năng lượng quốc gia: Cung cấp luận cứ khoa học vững chắc để các cơ quan quản lý nhà nước (Bộ Khoa học và Công nghệ, Viện Năng lượng nguyên tử Việt Nam) xây dựng chiến lược làm chủ công nghệ chu trình nhiên liệu hạt nhân tiên tiến, tận dụng nguồn tài nguyên quặng Thori tiềm năng tại Việt Nam trong tương lai.
  • Lợi ích xã hội và môi trường: Giảm thiểu tới $90%$ khối lượng chất thải phóng xạ độc hại chu kỳ bán rã dài cần chôn lấp địa chất sâu, đồng thời mang lại nguồn năng lượng điện hạt nhân an toàn tuyệt đối, loại trừ nguy cơ xảy ra sự cố thảm họa hạt nhân.

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

  • Nghiên cứu sinh và Giới học thuật Hạt nhân: Tiếp cận phương pháp luận mô phỏng ghép nối JENDL-HE và MCNPX, kế thừa bộ số liệu phân bố không gian và phổ năng lượng neutron để phát triển các đề tài chuyên sâu về vật lý lò phản ứng.
  • Kỹ sư R&D và Thiết kế Lò phản ứng: Sử dụng các kết quả tối ưu hóa tỷ lệ nhiên liệu $\text{ThO}_2 - \text{UO}_2$ và thông số hình học lò TRIGA Mark II làm tài liệu thiết kế cơ sở (Basic Engineering Design) cho các tổ hợp lò phản ứng nghiên cứu dưới tới hạn.
  • Cơ quan Quản lý An toàn Bức xạ và Hạt nhân: Nắm bắt các đặc tính an toàn nội tại của mô hình ADSR chì lỏng để xây dựng các tiêu chuẩn, quy chuẩn kỹ thuật quốc gia về đánh giá an toàn lò phản ứng thế hệ mới.
  • Nhà hoạch định Chiến lược Năng lượng Quốc gia: Có thêm cơ sở dữ liệu định lượng để xây dựng kịch bản đa dạng hóa nguồn nhiên liệu sơ cấp, giảm phụ thuộc vào nhập khẩu nhiên liệu urani làm già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à nó mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc xác lập cơ chế sinh và phân bố neutron spallation trên môi trường chì lỏng đồng thời làm bia tương tác và chất tải nhiệt. Luận án đã mở rộng lý thuyết va chạm nội hạt nhân (Intranuclear Cascade) của Serber và mô hình phân rã thống kê (Statistical Decay Model) sang môi trường kim loại lỏng liên tục, chứng minh rằng sự suy giảm mật độ khối ($10{,}66\text{ g/cm}^3$ so với $11{,}7\text{ g/cm}^3$) được bù trừ hoàn hảo nhờ mở rộng chiều dài tương tác hiệu dụng dọc theo toàn bộ lõi lò, mang lại thông lượng neutron đồng đều hơn mà không làm suy giảm hiệu suất $Y_{n/p}$.

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

So với nghiên cứu của S. Meigo và cộng sự (1999) và A. Goldenbaum và cộng sự (2000) (vốn chỉ khảo sát bia chì rắn tĩnh với chiều dài giới hạn $15 - 55\text{ cm}$ bằng MCNP4A), luận án đã thực hiện bước đột phá phương pháp luận: trích xuất toàn diện ma trận vi phân $32 \times 19$ từ thư viện JENDL-HE-2007 lên đến $3\text{ GeV}$, tích hợp mô phỏng 3D trên MCNPX cho toàn bộ cấu trúc hình học phức hợp của lò TRIGA Mark II, tính toán đồng thời cả phổ phân rã phóng xạ của Thori lẫn trường phân bố thông lượng neutron $\Phi(r, z, E)$.

3. Phát hiện bất ngờ nhất (surprising finding) từ dữ liệu thực nghiệm mô phỏng là gì?

Phát hiện bất ngờ nhất là tính chất "tự điều hòa phổ" (self-moderating spectral behavior) của khối chì lỏng: Trong khi phản ứng spallation tạo ra các neutron hướng trước mang năng lượng cực cao ($> 50\text{ MeV}$ ở góc $< 30^\circ$), thì sự tán xạ không đàn hồi liên tục trong thể tích chì lỏng lớn đã chuyển hóa phổ neutron ở vùng nhiên liệu ngoại vi ($\theta > 60^\circ$) thành phổ neutron nhiệt và trên nhiệt tối ưu, giúp tỷ lệ bắt neutron của $^{232}\text{Th}$ diễn ra mạnh mẽ mà không cần bổ sung các khối làm chậm bằng than chì hay nước nặng cồng kềnh.

4. Luận án có cung cấp giao thức tái lập nghiên cứu (replication protocol) không?

Có. Luận án cung cấp đầy đủ và minh bạch toàn bộ các phương trình toán học giải tích, cấu trúc ma trận góc - năng lượng trích xuất từ JENDL-HE-2007, 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, mật độ khối lượng các vật liệu chì lỏng, tỷ lệ phần trăm khối lượng hỗn hợp nhiên liệu Thori - Urani, cùng các lệnh thiết lập nguồn trong chương trình MCNPX. Bất kỳ nhóm nghiên cứu nào có trang bị mã MCNPX và thư viện JENDL-HE đều có thể tái lập chính xác $100%$ các kết quả tính toán.

5. Lộ trình nghiên cứu 10 năm (10-year research agenda) tiếp theo được phác thảo như thế nào?

Lộ trình 10 năm gồm 3 giai đoạn:

  • Năm 1 - 3: Xây dựng mô hình ghép nối đa vật lý thủy nhiệt - neutron (CFD-MCNPX) và tính toán động học suy kiệt nhiên liệu (burnup kinetics) theo thời gian thực.
  • Năm 4 - 6: Thực hiện các thực nghiệm kiểm chứng chùm hạt quy mô nhỏ (benchmarking experiments) trên các máy gia tốc hiện có (như tại Viện KIPT hoặc JAEA) với bia chì lỏng thu nhỏ.
  • Năm 7 - 10: Hoàn thiện thiết kế kỹ thuật chi tiết cho một tổ hợp lò phản ứng nghiên cứu ADSR công suất nhiệt $10 - 50\text{ MWth}$ sử dụng nhiên liệu Thori tại Việt Nam, hướng tới thử nghiệm nạp tải chì lỏng thực tế.

Kết luận

  1. Thiết lập thành công mô hình vật lý ADSR sử dụng bia chì lỏng kép: Chứng minh tính khả thi vượt trội của phương án dùng chì lỏng vừa làm bia sinh neutron spallation vừa làm chất tải nhiệt, giải quyết triệt để bài toán hỏng hóc bia rắn và loại bỏ nguy cơ phát sinh độc tố $\text{Po-210}$.
  2. Xác lập bộ dữ liệu neutron học toàn diện: Tính toán chính xác hiệu suất phát neutron $Y_{n/p}$ (từ $10$ đến $> 55\text{ n/p}$ ở dải $0{,}25 - 3{,}0\text{ GeV}$), vi phân bậc hai tiết diện sinh neutron $\frac{d^2\sigma}{dEd\Omega}$, và phổ phóng xạ $\alpha, \beta, \gamma, \bar{\nu}$ của Thori trong chì lỏng.
  3. Tối ưu hóa cấu hình nhiên liệu Thori trong lò TRIGA Mark II: Xác định tỷ lệ phối trộn vàng của hệ nhiên liệu $\text{Th-}^{233}\text{UO}2$ (hàm lượng thori $60% - 80%$) và $\text{Th-}^{235}\text{UO}2$ (hàm lượng thori $40% - 60%$), duy trì hệ số nhân neutron hiệu dụng $k{eff}$ ổn định trong dải an toàn dưới tới hạn $0{,}90 \le k{eff} < 1{,}00$.
  4. Làm sáng tỏ quy luật phân bố không gian thông lượng neutron $\Phi(r, z)$: Khẳng định độ đồng đều cao của trường neutron trong môi trường chì lỏng, giúp thanh nhiên liệu cháy kiệt đều và tối ưu hóa hiệu suất truyền nhiệt.
  5. Mở ra 3 nhánh nghiên cứu tiên phong cho tương lai: Tạo tiền đề cho các hướng nghiên cứu ghép nối đa vật lý Neutron-CFD, công nghệ vật liệu vỏ bọc chịu ăn mòn kim loại lỏng, và chiến lược thiêu hủy chất thải phóng xạ chu kỳ dài bằng chu trình nhiên liệu Thori an toàn nội tại.