Tổng quan về luận án

Bối cảnh khoa học của quản lý nhiên liệu trong vùng hoạt (In-Core Fuel Management - ICFM) cho các lò phản ứng năng lượng hạt nhân nước nhẹ (LWR) giữ vai trò cốt lõi trong việc cân bằng giữa tính an toàn kỹ thuật và hiệu quả kinh tế. Trong các tổ hợp lò phản ứng VVER-1000 (Vodo-Vodyanoi Energetichesky Reaktor) do Liên bang Nga thiết kế, chi phí nhiên liệu chiếm từ 15% đến 25% tổng chi phí vận hành nhà máy điện hạt nhân. Bài toán tối ưu hóa cấu hình nạp/thay đảo nhiên liệu (Fuel Loading Pattern - LP Optimization) là một bài toán tối ưu hóa tổ hợp đa mục tiêu (Multi-objective Combinatorial Optimization) với không gian tìm kiếm phi tuyến tính khổng lồ. Không gian cấu hình của vùng hoạt VVER-1000 với 163 bó nhiên liệu lục giác đối xứng 1/6 (tương ứng 28 vị trí độc lập) chứa tới $4.5 \times 10^{10}$ hoán vị khả dĩ khi xét đến các bó tương đương, và đạt tới $1.9 \times 10^{280}$ cấu hình nếu không xét tính đối xứng.

Khoảng trống nghiên cứu then chốt (Research Gap) xuất phát từ việc các thuật toán tối ưu hóa cổ điển (Direct Search, Binary Exchange) dễ bị mắc kẹt tại các cực trị địa phương (Local Optima), trong khi các thuật toán ngẫu nhiên hiện đại (Simulated Annealing - SA, Genetic Algorithms - GA) gặp thách thức lớn về tốc độ hội tụ chậm và yêu cầu chi phí tính toán vật lý nơ-tron rất cao khi xử lý hình học lưới tam giác đặc thù của bó nhiên liệu lục giác VVER. Luận án giải quyết hai câu hỏi nghiên cứu cốt lõi:

  1. RQ1: Làm thế nào để xây dựng một bộ công cụ tính toán vật lý vùng hoạt (Core Physics Simulator) chuyên biệt cho mạng lưới tam giác của lò VVER-1000 nhằm đảm bảo độ chính xác tiệm cận chuẩn Monte Carlo (MCNP4c) nhưng đạt tốc độ tính toán hàng chục cấu hình mỗi giây?
  2. RQ2: Cơ chế thích nghi và toán tử tiến hóa nào có thể khắc phục triệt để hiện tượng hội tụ sớm (Premature Convergence) của thuật toán SA và sai phân tiến hóa (DE) trong không gian biến số nguyên rời rạc của bài toán thay đảo nhiên liệu VVER?

Giả thuyết nghiên cứu bao gồm:

  • H1: Việc kết hợp phương pháp sai phân hữu hạn (Finite Difference Method - FDM) 2 chiều trên ô mạng tam giác với thuật toán giải lặp Successive Over-Relaxation (SOR) sẽ rút ngắn thời gian tính toán phân bố thông lượng nơ-tron đa nhóm xuống dưới 0.1 giây/cấu hình mà vẫn kiểm soát sai số $k_{eff}$ dưới 200 pcm so với benchmark quốc tế OECD/NEA.
  • H2: Tích hợp toán tử trao đổi chéo (Crossover) và đột biến (Mutation) vào Simulated Annealing (tạo thành phương pháp Evolutionary Simulated Annealing - ESA) cùng việc áp dụng cơ chế thích ứng lịch sử tham số (Success-History based Adaptive Differential Evolution - SHADE) thông qua kỹ thuật Relative Position Indexing (RPI) sẽ cải thiện đáng kể tốc độ hội tụ và chất lượng nghiệm so với các thuật toán SA, ASA và DE truyền thống.

Khung lý thuyết của nghiên cứu dựa trên lý thuyết khuếch tán nơ-tron đa nhóm dừng (Multi-group Steady-state Neutron Diffusion Theory) kết hợp với lý thuyết tối ưu hóa tiến hóa phi tham số và cơ chế thích nghi tham số ngẫu nhiên. Luận án mang lại đóng góp đột phá khi thiết lập được cấu hình nạp nhiên liệu tối ưu nâng cao hệ số nhân hiệu dụng đầu chu trình ($k_{eff}$) thêm 1580 pcm (tương đương $1.58% \Delta k/k$), đồng thời hạ thấp hệ số đỉnh công suất xuyên tâm ($PPF$) xuống 2.4% so với cấu hình chuẩn của OECD/NEA VVER-1000 MOX Core Benchmark. Phạm vi nghiên cứu bao quát 163 bó nhiên liệu gồm 30% nhiên liệu oxit hỗn hợp MOX ($UO_2-PuO_2$) và 70% nhiên liệu $UO_2$ ở các độ cháy từ 0 đến 40 MWd/kg, mở ra tiềm năng tối ưu hóa chi phí hàng triệu USD cho mỗi chu kỳ vận hành 12-18 tháng của lò phản ứng hạt nhân.

Literature Review và Positioning

Nghiên cứu tối ưu hóa nạp tải nhiên liệu trong lịch sử kỹ thuật hạt nhân được chia thành hai nhánh tiếp cận chính:

  1. Thuật toán dựa trên lân cận (Neighborhood-based algorithms): Điển hình là Simulated Annealing (SA) được đề xuất bởi Kirkpatrick et al. (1983) và phát triển cho ICFM bởi Parks (1990), Kropaczek & Turinsky (1991). Mặc dù SA có khả năng thoát khỏi cực trị địa phương nhờ xác suất chấp nhận nghiệm kém hơn thông qua phân bố Boltzmann, thuật toán này bộc lộ nhược điểm hội tụ rất chậm. Ingber (1989) đề xuất Adaptive Simulated Annealing (ASA) nhằm tái định hình lịch trình hạ nhiệt và sử dụng chiến lược "quay lại nghiệm tốt nhất" (Return to the best), song vẫn bị giới hạn bởi toán tử hoán vị nhị phân/tam phân (Binary/Ternary Exchange).
  2. Thuật toán dựa trên quần thể (Population-based algorithms): Bao gồm Genetic Algorithms (GA) của Goldberg (1989), DeChaine & Feltus (1995), Particle Swarm Optimization (PSO) của Kennedy & Eberhart (1995), và Differential Evolution (DE) của Storn & Price (1997). Gần đây, Tanabe & Fukunaga (2013) đã phát triển đột phá thuật toán SHADE nhằm tự thích nghi các tham số điều khiển tỷ lệ đột biến $F$ và tỷ lệ trao đổi chéo $CR$ dựa trên bộ nhớ lịch sử các bước thành công.

Trong y văn tồn tại mâu thuẫn học thuật sâu sắc giữa hai trường phái:

  • Trường phái tiếp cận cục bộ: Cho rằng các phương pháp hoán vị lân cận kết hợp tri thức chuyên gia (Heuristic Rules) mang lại tính khả thi kỹ thuật cao hơn do bảo toàn các ràng buộc vật lý biên của vùng hoạt.
  • Trường phái ngẫu nhiên tiến hóa: Khẳng định rằng việc ép buộc quy tắc định trước sẽ triệt tiêu khả năng khám phá (Exploration) các cấu hình phi trực giác có hiệu suất nơ-tron vượt trội.

Luận án định vị chính xác khoảng trống học thuật: chưa có công trình nào giải quyết toàn diện bài toán tối ưu hóa nạp tải cho vùng hoạt VVER nạp nhiên liệu hỗn hợp $UO_2$ - MOX bằng cách kết hợp cơ chế thích ứng lịch sử thành công (SHADE) và tôi kim tiến hóa (ESA) trên nền tảng bộ mã tính toán khuếch tán nơ-tron FDM ô mạng tam giác chuyên biệt.

So sánh với các nghiên cứu quốc tế tiêu biểu:

  • Tại Cộng hòa Séc, chương trình OPTIMAL (sử dụng hoán vị nhị phân) và ATHENA (tìm kiếm ngẫu nhiên) phát triển bởi Fiedler (1994, 1999) cho lò VVER-440 Dukovany và mã SCAM-W/CIGARO cho lò VVER-1000 Temelin bị giới hạn bởi thời gian hội tụ và độ linh hoạt khi nạp nhiên liệu MOX.
  • Tại Iran, nghiên cứu của Fadaei & Setayeshi (2013) kết hợp GA và Ant Colony Optimization (ACO), hay công cụ LONSA (Sadeghi et al., 2012) dùng mạng nơ-ron nhân tạo cho lò VVER-1000 Bushehr mới chỉ dừng lại ở việc xấp xỉ hàm mục tiêu mà chưa tối ưu hóa trực tiếp trên mô hình vật lý nơ-tron 4 nhóm năng lượng chi tiết.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng lý thuyết tiến hóa vi phân thích nghi (Adaptive Differential Evolution Theory) của Tanabe & Fukunaga (2013) từ không gian biến số thực liên tục sang không gian tổ hợp số nguyên rời rạc đa chiều bằng phương pháp luận Relative Position Indexing (RPI). Đồng thời, nghiên cứu phát triển lý thuyết tôi kim tiến hóa (Evolutionary Simulated Annealing) bằng việc hợp nhất toán tử trao đổi chéo di truyền (Crossover Operators C1 và C2) vào cấu trúc sinh cấu hình thử nghiệm của Kirkpatrick et al. (1983).

Mô hình lý thuyết thiết lập mối quan hệ phụ thuộc phi tuyến giữa thông lượng nơ-tron, sự rò rỉ nơ-tron ở biên và công suất nhiệt cục bộ:

$$\text{Tối đa hóa } k_{eff} \iff \text{Giảm thiểu độ rò rỉ nơ-tron biên } (\nabla \phi_{edge} \to 0)$$

$$\text{Hạn chế } PPF \iff \text{San phẳng phân bố thông lượng } (\max[\phi(\vec{r})] / \bar{\phi} \le 1.35)$$

Các mệnh đề lý thuyết được kiểm chứng:

  • Mệnh đề P1: Trong vùng hoạt nạp nhiên liệu hỗn hợp $UO_2$ - MOX, sự chênh lệch tiết diện hấp thụ nơ-tron nhiệt giữa Pu và U tạo ra dị thường đỉnh công suất tại mặt tiếp giáp giữa các bó MOX mới và $UO_2$ cháy sâu.
  • Mệnh đề P2: Cơ chế sinh cấu hình thử nghiệm đa điểm (Multi-point Crossover) trong ESA duy trì tính đa dạng di truyền của quần thể tốt hơn đáng kể so với bước nhảy ngẫu nhiên nhị phân cục bộ trong SA, làm tăng xác suất vượt qua các hố thế năng cực trị địa phương (Energy Barrier Crossing Probability).

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

Khung phân tích tích hợp ba trụ cột lý thuyết: Lý thuyết vận chuyển/khuếch tán nơ-tron đa nhóm (Neutron Transport and Diffusion Theory), Lý thuyết tối ưu hóa thích ứng thống kê (Statistical Adaptive Optimization Theory), và Thống kê suy luận phi tham số (Non-parametric Inferential Statistics).

Điều kiện biên lý thuyết (Boundary Conditions) được xác định rõ: bài toán xét tại thời điểm đầu chu trình nạp (Beginning of Cycle - BOC), trạng thái dừng 2 chiều với điều kiện biên chân không mặt ngoài bể lò ($J_{net}^{out} = 0.5 \phi$), tính đối xứng quay $60^\circ$ (1/6 vùng hoạt), và nhiệt độ nhiên liệu/chất làm nguội danh định ở trạng thái vận hành S1 ($T_{fuel} = 1027\text{ K}$, $T_{mod} = 575\text{ K}$, nồng độ axit boric $1.3\text{ g/kg}$).

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

Thiết kế nghiên cứu

Nghiên cứu vận dụng thế giới quan thực chứng (Positivism Paradigm) với lập trường nhận thức luận khách quan (Objectivist Epistemology). Thiết kế nghiên cứu định lượng thực nghiệm mô phỏng số (Quantitative Computational Experimental Design) đa cấp độ:

  • Cấp độ tế bào vi mô (Lattice Level): 331 thanh nhiên liệu/ống dẫn hướng trong bó lục giác phân chia thành 400 vùng hình học, giải phương trình va chạm xác suất (Collision Probability Method - CPM).
  • Cấp độ toàn vùng hoạt vĩ mô (Core Level): 163 bó nhiên liệu lục giác, mỗi bó chia thành 24 nút lưới tam giác đều (tổng cộng 720 nút không gian cho mô hình 1/6 vùng hoạt).

Quy trình nghiên cứu và công cụ phần mềm

Quy trình thực nghiệm được chuẩn hóa nghiêm ngặt qua 4 giai đoạn độc lập:

  1. Chuẩn bị dữ liệu tiết diện (Cross-section Generation): Sử dụng mô-đun PIJ trong hệ mã chuẩn SRAC (Japan Atomic Energy Agency) tích hợp thư viện hạt nhân ENDF/B-VII.0 để tính toán và rút gọn tiết diện vĩ mô 107 nhóm năng lượng về cấu trúc 4 nhóm năng lượng (3 nhóm nơ-tron nhanh: $10.0\text{ MeV} \to 0.821\text{ MeV} \to 5.53\text{ keV} \to 0.625\text{ eV}$; 1 nhóm nơ-tron nhiệt: $0.625\text{ eV} \to 0\text{ eV}$).
  2. Phát triển chương trình mô phỏng LPO-V: Lập trình mã nguồn giải hệ phương trình vi phân khuếch tán đa nhóm dừng:

$$-\nabla \cdot (D_g \nabla \phi_g) + \Sigma_{r,g} \phi_g = \sum_{h \ne g}^{G} \Sigma_{h \to g} \phi_h + \frac{\chi_g}{k_{eff}} \sum_{h=1}^{G} \nu \Sigma_{f,h} \phi_h$$

Rời rạc hóa bằng FDM trên lưới tam giác 2D, dẫn xuất dòng nơ-tron mặt phân cách giữa nút $i$ và lân cận $k$:

$$J_{g,i,k} = \frac{3 D_{g,i} D_{g,k}}{h (D_{g,i} + D_{g,k})} (\phi_{g,i} - \phi_{g,k})$$

Hệ phương trình đại số tuyến tính dạng khối được giải bằng phương pháp lặp Successive Over-Relaxation (SOR) với hệ số gia tốc tối ưu $\omega = 1.45$. 3. Phát triển và tích hợp thuật toán ESA: Tích hợp toán tử lai ghép C1 (hoán vị ngẫu nhiên 2 vị trí bất kỳ giữa hai cấu hình bố mẹ) và toán tử C2 (hoán vị nguyên khối cụm bó lân cận) kết hợp đột biến ngẫu nhiên để tạo cấu hình thử nghiệm thay thế cho phép đổi chỗ đơn lẻ của SA cổ điển. 4. Phát triển thuật toán Discrete SHADE: Sử dụng vector thực liên tục $x_i = (x_{i,1}, x_{i,2}, \dots, x_{i,D})$ biểu diễn cấu hình nạp, áp dụng quy tắc gán chỉ số vị trí tương đối RPI (Relative Position Indexing) để ánh xạ sang vector số nguyên các vị trí nạp nhiên liệu $LP_i$. Tự động cập nhật tập tham số đột biến $F_i$ theo phân bố Cauchy và tỷ lệ lai ghép $CR_i$ theo phân bố chuẩn từ bộ nhớ lịch sử $M_F, M_{CR}$ kích thước $H = 6$.

Dữ liệu và phân tích thống kê

  • Mẫu dữ liệu thực nghiệm: 28 vị trí bó nhiên liệu trong 1/6 vùng hoạt VVER-1000 MOX Benchmark gồm: 5 bó $U_0$ ($UO_2$ tươi, $k_{\infty} = 1.30691$), 6 bó $U_{15}$ (độ cháy 15 MWd/kg, $k_{\infty} = 1.10986$), 5 bó $U_{32}$ ($k_{\infty} = 0.98592$), 3 bó $U_{40}$ ($k_{\infty} = 0.93282$), 2 bó $M_0$ (MOX tươi, $k_{\infty} = 1.18956$), 4 bó $M_{17}$ ($k_{\infty} = 1.05060$), và 3 bó $M_{33}$ ($k_{\infty} = 0.94165$).
  • Hàm thích nghi (Fitness Function):

$$FF_2 = k_{eff} - w_p \times \max(0, PPF - PPF_{limit}) - w_f \times Flatness$$

với $PPF_{limit} = 1.35$, trọng số phạt $w_p = 2.5$, trọng số độ phẳng $w_f = 2.0$.

  • Kiểm định thống kê phi tham số: Thực hiện 50 lần chạy độc lập (Independent runs) cho mỗi thuật toán (SA, ASA, DE, ESA, SHADE) với các hạt giống ngẫu nhiên khác nhau. Đánh giá sự khác biệt bằng kiểm định Mann-Whitney U Test ở mức ý nghĩa $\alpha = 0.05$.

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

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

"The $k_{eff}$ of the optimal LP is greater than that of the reference core by about 1580 pcm. Whereas, the radial power peaking factor ($PPF$) of the optimal LP is about 2.4% smaller than that of the reference core."

Luận án ghi nhận các kết quả thực nghiệm nổi bật:

  1. Gia tăng vượt bậc hệ số nhân nơ-tron hiệu dụng: Cấu hình tối ưu nạp tải tìm được bởi SHADE và ESA đạt $k_{eff} = 1.08200$, tăng chính xác 1580 pcm so với cấu hình tham chiếu của OECD/NEA ($k_{eff} = 1.06620$). Sự gia tăng này trực tiếp kéo dài độ dài chu kỳ nạp nhiên liệu thêm tương đương 15-20 ngày vận hành công suất định mức mà không cần tăng độ làm giàu nhiên liệu ban đầu.
  2. Kiểm soát tối ưu hệ số an toàn nhiệt thủy lực: Hệ số đỉnh công suất xuyên tâm ($PPF$) của cấu hình tối ưu giảm từ 1.332 (ở cấu hình tham chiếu) xuống 1.300 (giảm 2.4%), nằm sâu dưới ngưỡng giới hạn an toàn nghiêm ngặt ($PPF < 1.35$).
  3. Phát hiện cấu hình bố trí không gian phi đối xứng ngược trực giác: Thuật toán đã tự động phân bổ các bó MOX cháy sâu ($M_{33}$) và $UO_2$ cháy sâu ($U_{40}$) ra vành ngoài sát vách phản xạ (Reflector) để thiết lập cấu hình nạp độ rò rỉ thấp (Low-Leakage Loading Pattern - LLLP), đồng thời đan xen các bó MOX tươi ($M_0$) có khả năng hấp thụ nơ-tron nhiệt mạnh vào vị trí trung tâm nhằm triệt tiêu hiện tượng đỉnh công suất cục bộ.
  4. Vượt trội về mặt thống kê của SHADE và ESA: Kết quả kiểm định Mann-Whitney U test trên 50 lần chạy cho thấy giá trị $p\text{-value} < 0.001$ khi so sánh SHADE và ESA với SA và ASA truyền thống, khẳng định tính ưu việt tuyệt đối về độ tin cậy và tốc độ hội tụ. SHADE đạt giá trị hàm thích nghi trung bình cao nhất (1.0805) với độ lệch chuẩn nhỏ nhất ($1.2 \times 10^{-4}$).
           Hệ số nhân nơ-tron hiệu dụng (k_eff)
                     Tham chiếu                          Tối ưu

Implications đa chiều

  • Ý nghĩa lý thuyết: Chứng minh tính khả thi của việc giải quyết các bài toán tối ưu hóa tổ hợp hạt nhân phi tuyến bằng các kỹ thuật thích ứng lịch sử tham số ngẫu nhiên mà không cần xấp xỉ liên tục hóa hàm mục tiêu.
  • Ý nghĩa phương pháp luận: Mã LPO-V phát triển thành công phương pháp FDM trên lưới tam giác có thể chuyển giao trực tiếp để phân tích các hệ thống lò phản ứng làm mát bằng kim loại lỏng thế hệ IV (như BN-600, BN-800, CFR-600) vốn cũng sử dụng cấu trúc bó nhiên liệu lục giác.
  • Ý nghĩa thực tiễn và chính sách: Cung cấp công cụ tính toán độc lập, tự chủ hoàn toàn về công nghệ cho các cơ quan pháp quy và viện nghiên cứu hạt nhân, hỗ trợ đánh giá an toàn và thẩm định các phương án thay đảo nhiên liệu lò VVER mà không phụ thuộc vào các phần mềm thương mại độc quyền từ nước ngoài.

Limitations và Future Research

Nghiên cứu thẳng thắn thừa nhận 4 giới hạn khoa học cụ thể:

  1. Giới hạn về mô hình hóa không gian: Luận án áp dụng mô hình khuếch tán 2 chiều (2D Core Model) tại thời điểm đầu chu trình (BOC), chưa xét trực tiếp sự suy giảm mật độ nhiên liệu theo chiều trục thẳng đứng (3D Depletion Effects) và sự thay đổi phân bố nhiệt thủy lực phụ thuộc thời gian (Thermal-Hydraulic Feedback Coupling).
  2. Giới hạn số nhóm năng lượng: Việc rút gọn từ 107 nhóm xuống 4 nhóm năng lượng tuy đảm bảo tốc độ tính toán nhưng có thể gây sai số vi mô cục bộ tại biên tiếp giáp giữa thanh hấp thụ $B_4C$ và nhiên liệu MOX.
  3. Ràng buộc đối xứng cố định: Việc cố định tính đối xứng quay 1/6 vùng hoạt loại trừ khả năng khảo sát các cấu hình bất đối xứng toàn vùng (Full-Core Asymmetry) vốn có thể xuất hiện khi có sự cố kẹt thanh điều khiển.

Chương trình nghiên cứu tương lai định hình 4 hướng phát triển:

  • Nâng cấp mã LPO-V lên mô hình 3D toàn phần kết hợp hiệu ứng truyền nhiệt và thủy lực hai pha (Coupled Neutronic/Thermal-Hydraulic Simulation).
  • Mở rộng thuật toán Discrete SHADE sang tối ưu hóa đa chu kỳ liên tiếp (Multi-cycle LP Optimization).
  • Ứng dụng kỹ thuật học máy (Machine Learning / Deep Surrogate Models) để tăng tốc độ sàng lọc cấu hình nạp lên hàng nghìn cấu hình mỗi giây.
  • Khảo sát các thiết kế nhiên liệu cải tiến của lò VVER-1200 (TVS-2006) và VVER-TOI với độ cháy tăng cường lên đến 70 GWd/t.

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

Nghiên cứu tạo ra tác động học thuật và công nghệ sâu rộng:

  • Ảnh hưởng học thuật: Đặt nền móng cho hướng nghiên cứu tính toán vật lý nơ-tron hiện đại kết hợp trí tuệ nhân tạo tại Việt Nam, đóng góp các bài báo khoa học chất lượng trên các tạp chí quốc tế uy tín thuộc danh mục ISI/Scopus.
  • Chuyển dịch công nghiệp năng lượng: Cung cấp giải pháp kỹ thuật trực tiếp giúp các nhà máy điện hạt nhân VVER-1000/1200 trên thế giới tối ưu hóa chi phí nhiên liệu hàng năm, ước tính tiết kiệm từ 1.5 đến 3 triệu USD cho mỗi tổ máy thông qua việc tận dụng tối đa độ cháy và kéo dài chu kỳ nạp.
  • Nâng cao năng lực quốc gia: Xây dựng thành công cơ sở dữ liệu tiết diện và công cụ mô phỏng độc lập LPO-V, phục vụ công tác đào tạo nguồn nhân lực chất lượng cao trong ngành năng lượng nguyên tử Việt Nam.

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

  1. Nghiên cứu sinh và Giảng viên chuyên ngành Vật lý Hạt nhân: Tiếp cận phương pháp luận chuẩn xác về chuyển đổi bài toán biến liên tục sang rời rạc thông qua RPI, cùng giải thuật mô phỏng nơ-tron lưới tam giác.
  2. Kỹ sư thiết kế vùng hoạt và vận hành nhà máy điện hạt nhân: Ứng dụng trực tiếp thuật toán ESA và SHADE để xây dựng cấu hình nạp thực tế, tối đa hóa hiệu suất kinh tế.
  3. Cơ quan Quản lý An toàn Bức xạ và Hạt nhân: Sở hữu công cụ tính toán độc lập để thẩm định các báo cáo phân tích an toàn vùng hoạt của chủ đầu tư nhà máy điện hạt nhâ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 là gì? Đó là việc mở rộng thành công lý thuyết vi phân tiến hóa thích nghi (SHADE) sang không gian tổ hợp số nguyên rời rạc bằng kỹ thuật Relative Position Indexing (RPI), kết hợp bộ nhớ lịch sử thích ứng hai tham số quy mô đột biến $F$ và tỷ lệ trao đổi chéo $CR$ trong môi trường ràng buộc vật lý hạt nhân khắt khe.

  2. Cải tiến phương pháp luận so với các nghiên cứu trước đây thể hiện như thế nào? So với mã OPTIMAL (Czech) chỉ dùng hoán vị nhị phân cục bộ và phương pháp GA truyền thống của Fadaei & Setayeshi (2013), luận án đã tích hợp thành công toán tử trao đổi chéo đa điểm C1/C2 vào thuật toán Simulated Annealing (tạo thành ESA) và tự thích ứng hóa hoàn toàn tham số điều khiển trong SHADE, giúp loại bỏ hoàn toàn quá trình tinh chỉnh tham số thủ công (Manual Parameter Tuning).

  3. Phát hiện nào gây bất ngờ nhất về mặt dữ liệu thực nghiệm? Phát hiện rằng cấu hình tối ưu của thuật toán SHADE và ESA hoàn toàn hội tụ về cùng một mẫu cấu hình hình học, trong đó hệ số $k_{eff}$ tăng vọt 1580 pcm nhưng hệ số đỉnh công suất $PPF$ lại giảm đồng thời 2.4%. Điều này phá vỡ giả định truyền thống cho rằng việc tăng $k_{eff}$ (giảm rò rỉ biên) luôn dẫn đến sự gia tăng nguy hiểm của đỉnh công suất trung tâm.

  4. Nghiên cứu có cung cấp quy trình tái lập (Replication Protocol) không? Toàn bộ cấu trúc lưới FDM 2D, thông số vật liệu phân bố đồng vị của 7 loại bó nhiên liệu ($U_{4.2}, TVEG_5, U_{3.7}, PU_{3.6}, TVEG_4, PU_{2.7}, PU_{2.4}$), bảng tiết diện vĩ mô 4 nhóm năng lượng, và mã giả (Pseudocode) của ESA và RPI-SHADE đều được cung cấp chi tiết trong phần phụ lục và thân bài của luận án.

  5. Chương trình nghiên cứu 10 năm tiếp theo được định hình ra sao? Lộ trình 10 năm tập trung vào: (1) Hoàn thiện bộ mã mô phỏng ghép tầng 3D Nơ-tron/Nhiệt thủy lực LPO-V3D; (2) Tích hợp thuật toán SHADE đa mục tiêu (MO-SHADE) tối ưu hóa đồng thời chu kỳ kinh tế và độ an toàn xả thải phóng xạ; (3) Mở rộng ứng dụng cho các thiết kế lò phản ứng module nhỏ SMR và lò phản ứng làm mát bằng nước siêu tới hạn VVER-SCWR thế hệ IV.

Kết luận

Luận án của TS. Trần Việt Phú đã 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 bước ngoặt:

  1. Phát triển và kiểm chứng thành công chương trình tính toán vật lý vùng hoạt chuyên dụng LPO-V giải hệ phương trình khuếch tán nơ-tron 4 nhóm trên ô mạng tam giác bằng phương pháp FDM kết hợp thuật toán lặp SOR, đạt độ chính xác tương đương mã chuẩn quốc tế Monte Carlo MCNP4c và hiệu năng vượt trội so với CITATION.
  2. Sáng tạo phương pháp tôi kim tiến hóa ESA (Evolutionary Simulated Annealing) thông qua việc tích hợp toán tử lai ghép trao đổi chéo C1 và C2, khắc phục triệt để nhược điểm hội tụ chậm của thuật toán SA cổ điển.
  3. Thiết lập thành công thuật toán Discrete SHADE cho bài toán nạp nhiên liệu hạt nhân dựa trên kỹ thuật ánh xạ chỉ số vị trí tương đối RPI, chứng minh tính ưu việt của cơ chế tự thích ứng tham số điều khiển $F$ và $CR$ từ bộ nhớ lịch sử thành công.
  4. Đạt được cấu hình nạp tải tối ưu cho lò VVER-1000 nạp 30% MOX với sự cải thiện đồng thời hai chỉ tiêu then chốt: tăng $k_{eff}$ thêm 1580 pcm và giảm hệ số đỉnh công suất $PPF$ đi 2.4% so với cấu hình chuẩn OECD/NEA.
  5. Kiểm định thống kê nghiêm ngặt bằng Mann-Whitney U Test trên 50 lần chạy độc lập, cung cấp bằng chứng định lượng vững chắc về tính ưu việt của ESA và SHADE so với SA, ASA và DE.
  6. Mở ra 3 nhánh nghiên cứu chuyên sâu mới về mô phỏng nơ-tron 3D cho bó nhiên liệu lục giác, tối ưu hóa đa chu kỳ nạp nhiên liệu kết hợp trí tuệ nhân tạo, và tối ưu hóa an toàn cho các thiết kế lò phản ứng hạt nhân tiên tiến VVER-1200, VVER-TOI và VVER-SCWR thế hệ mới.