Luận án tiến sĩ nghiên cứu tối ưu thay đảo nhiên liệu lò phản ứng VVER-1000

Tổng quan về luận án

Bối cảnh chuyển dịch năng lượng toàn cầu và an ninh năng lượng quốc gia đặt ra yêu cầu cấp thiết về việc tối ưu hóa hiệu quả kinh tế và an toàn vận hành của các nhà máy điện hạt nhân. Trong cơ cấu chi phí phát điện hạt nhân, chi phí nhiên liệu chiếm tỷ trọng đáng kể, dao động từ 15% đến 25% tổng chi phí vận hành chu trình ("the fuel costs, accounting for about 15-25%, have a great impact on the cost of nuclear power"). Trong đó, bài toán quản lý nhiên liệu trong vùng hoạt (In-Core Fuel Management - ICFM), trọng tâm là tối ưu hóa cấu hình nạp tải và thay đảo nhiên liệu (Fuel Loading Pattern - LP Optimization), đóng vai trò then chốt nhằm tối đa hóa độ dài chu trình vận hành hoặc hệ số nhân hiệu dụng đầu chu trình ($k_{eff}$ at Beginning of Cycle - BOC), đồng thời duy trì các giới hạn an toàn nghiêm ngặt về nhiệt - thủy lực thông qua việc cực tiểu hóa hệ số đỉnh công suất xuyên tâm (Radial Power Peaking Factor - PPF).

Khoảng trống nghiên cứu (Research Gap) xuất phát từ bản chất tổ hợp phi tuyến tính cao và không gian tìm kiếm khổng lồ của bài toán LP đối với lò phản ứng nước áp lực hình học lục giác kiểu Nga VVER-1000. Trong một vùng hoạt thực tế gồm 163 bó nhiên liệu với các mức độ giàu $U^{235}$ và độ cháy (burnup) khác nhau, nếu không áp dụng đối xứng, không gian tìm kiếm có thể lên tới $1.9 \times 10^{280}$ cấu hình khả dĩ; ngay cả khi áp dụng đối xứng quay bậc 6 ($1/6$ rotational symmetry), không gian vẫn chứa tới $4.5 \times 10^{10}$ phương án ("The possible number of fuel LPs would be enormous even when geometrical symmetry and constraints are considered"). Các phương pháp tìm kiếm cổ điển như Direct Search (DS) hay Binary Exchange (BE) thường bị mắc kẹt tại các điểm tối ưu cục bộ (local optima), trong khi các thuật toán ngẫu nhiên tiên tiến (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 tài nguyên tính toán khổng lồ. Đặc biệt, việc áp dụng nhiên liệu oxit hỗn hợp (Mixed Oxide Fuel - MOX) chứa Plutonium với tiết diện hấp thụ nơtron nhiệt cao làm biến dạng sâu sắc phổ nơtron và phân bố thông lượng nhiệt, đòi hỏi công cụ mô phỏng vật lý nơtron và giải thuật tối ưu thích ứng đặc thù.

Luận án xác lập hai câu hỏi nghiên cứu và hai giả thuyết cốt lõi:

1. Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để xây dựng một bộ công cụ tính toán vật lý nơtron 2 chiều chuyên biệt cho hình học ô mạng tam giác lục giác của lò VVER-1000 với tốc độ tính toán nhanh vượt bậc nhưng vẫn duy trì độ chính xác tương đương các mã chuẩn mực Monte Carlo?
2. Câu hỏi nghiên cứu 2 (RQ2): Cơ chế lai ghép tiến hóa hoặc thích ứng tham số nào có khả năng vượt qua rào cản hội tụ chậm của thuật toán tôi kim mô phỏng (SA) và chuyển đổi không gian thực - rời rạc trong thuật toán tiến hóa vi phân (DE) để tìm kiếm cấu hình nạp tải tối ưu toàn cục cho vùng hoạt VVER-1000 tải 30% MOX?

• Giả thuyết 1 (H1): Việc giải hệ phương trình khuếch tán nơtron 4 nhóm năng lượng bằng phương pháp sai phân hữu hạn (FDM) trên lưới tam giác kết hợp phương pháp siêu thư giãn liên tiếp (SOR) trong mã tính LPO-V sẽ giảm thiểu thời gian tính toán hàng chục lần so với các mã khuếch tán truyền thống mà vẫn đảm bảo độ lệch $k_{eff}$ dưới $150\text{ pcm}$ so với mô phỏng Monte Carlo MCNP4c.
• Giả thuyết 2 (H2): Việc tích hợp toán tử trao đổi chéo (crossover) vào giải thuật tôi kim mô phỏng tiến hóa (ESA) và tích hợp cơ chế chỉ số vị trí tương đối (RPI) cùng lịch sử thích ứng tham số trong giải thuật tiến hóa vi phân thích ứng rời rạc (Discrete SHADE) sẽ cải thiện đáng kể tốc độ hội tụ, tăng $k_{eff}$ BOC thêm trên $1500\text{ pcm}$ và hạ thấp hệ số đỉnh công suất $PPF$ trên $2%$ so với cấu hình chuẩn tham chiếu OECD/NEA.

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 (Multi-group Neutron Diffusion Theory) kết hợp Định luật Fick và các nguyên lý Tối ưu hóa Tổ hợp Metaheuristic. Đóng góp đột phá được lượng hóa thông qua việc xây dựng thành công mã tính toán vật lý vùng hoạt LPO-V và thiết lập cấu hình nạp tải tối ưu giúp gia tăng $k_{eff}$ đầu chu trình lên xấp xỉ $1580\text{ pcm}$, đồng thời giảm hệ số đỉnh công suất $PPF$ xuống $2.4%$ so với cấu hình tham chiếu của OECD/NEA. Phạm vi nghiên cứu thực nghiệm bao quát mô hình $1/6$ vùng hoạt lò VVER-1000 gồm 28 bó nhiên liệu ($19\text{ UO}_2$ và $9\text{ MOX}$) tại 6 trạng thái vật lý chuẩn (S1–S6), khảo sát từ 10.000 đến 100.000 cấu hình nạp tải qua 50 lần chạy độc lập để đánh giá độ tin cậy thống kê.

Literature Review và Positioning

Lịch sử nghiên cứu tối ưu hóa thay đảo nhiên liệu hạt nhân đã trải qua nhiều giai đoạn phát triển, bắt đầu từ các phương pháp xác định cổ điển như hoán đổi nhị phân (Binary Exchange) của Fagan & Sesonske (1969) và tìm kiếm trực tiếp. Tuy nhiên, tính chất đa cực trị gồ ghề của hàm mục tiêu năng lượng vùng hoạt đã thúc đẩy sự chuyển dịch sang các thuật toán ngẫu nhiên phỏng sinh học và vật lý thống kê.

Trong nhóm thuật toán dựa trên lân cận (Neighborhood-based algorithms), phương pháp Tôi kim Mô phỏng (Simulated Annealing - SA) do Kirkpatrick et al. (1983) đề xuất đã được Kropaczek & Turinsky (1991), Parks (1990) ứng dụng tiên phong vào ICFM nhờ khả năng chấp nhận các nghiệm kém hơn theo phân bố Boltzmann để thoát khỏi cực trị cục bộ. Dẫu vậy, SA thuần túy bộc lộ nhược điểm cố hữu là tốc độ làm nguội chậm, đòi hỏi số lượng hàm mục tiêu đánh giá khổng lồ. Ingber (1993) phát triển Tôi kim Mô phỏng Thích ứng (ASA) với lịch sử lưu trữ nghiệm nhằm tăng tốc độ làm nguội, song nguy cơ co hẹp không gian tìm kiếm sớm vẫn hiện hữu khi đối mặt với vùng hoạt có nạp tải nhiên liệu MOX dị bản.

Trong nhóm thuật toán dựa trên quần thể (Population-based algorithms), Thuật toán Di truyền (GA) của Holland (1975) và Goldberg (1989) được DeChaine & Feltus (1995), Chapot et al. (1999) triển khai rộng rãi cho lò phản ứng nước nhẹ hình học vuông (PWR, BWR). Tuy nhiên, GA bộc lộ hạn chế về độ phức tạp khi mã hóa nhiễm sắc thể cho các bài toán hoán vị số nguyên không lặp lại. Thuật toán Tiến hóa Vi phân (Differential Evolution - DE) do Storn & Price (1997) giới thiệu chứng minh tính vượt trội về khả năng thăm dò không gian và độ ổn định so với GA và Particle Swarm Optimization (PSO - Eberhart & Kennedy, 1995). Bước tiến lớn gần đây là thuật toán SHADE (Success-History based Adaptive Differential Evolution) do Tanabe & Fukunaga (2013) phát triển, ứng dụng bộ nhớ lịch sử thành công để tự động điều chỉnh hệ số khuếch đại đột biến $F$ và tỷ lệ trao đổi chéo $CR$.

Các tranh luận học thuật then chốt tập trung vào hai quan điểm đối lập:

1. Quan điểm trường phái Tôi kim Thống kê: Nhấn mạnh rằng các thuật toán cục bộ cải tiến (ASA, Tabu Search) có độ hội tụ tiệm cận chắc chắn và kiểm soát biên độ di chuyển hạt nhân chặt chẽ hơn.
2. Quan điểm trường phái Tiến hóa Quần thể: Khẳng định các thuật toán tiến hóa (DE, GA) sở hữu khả năng phân tán đa hướng vượt trội, tránh hiện tượng tắc nghẽn nghiệm khi không gian tìm kiếm chứa các rào cản tiềm năng sâu do sự không đồng nhất giữa các bó nhiên liệu $\text{UO}_2$ và MOX.

Vị trí của luận án được xác lập tại điểm giao thoa giữa hai trường phái: cải tiến giải thuật SA bằng cách tích hợp toán tử trao đổi chéo quần thể của GA tạo nên phương pháp lai ESA (Evolutionary Simulated Annealing), đồng thời tái cấu trúc thuật toán SHADE vốn dành cho không gian biến liên tục thành phiên bản rời rạc chuyên dụng cho bài toán hoán vị bó nhiên liệu hình học lục giác thông qua phương pháp Relative Position Indexing (RPI).

Khi đối chiếu với các nghiên cứu quốc tế về lò VVER:

• So sánh với các nghiên cứu tại Cộng hòa Séc trên lò VVER-440 Dukovany sử dụng mã tính OPTIMAL và ATHENA (Faber et al., 1994; Chladek et al., 1999) hoặc lò VVER-1000 Temelin sử dụng SCAM-W/CIGARO (Kostal et al., 2012): Các công cụ trước đây chủ yếu dựa trên trao đổi nhị phân thuần túy hoặc GA cổ điển, dẫn tới hiệu suất tìm kiếm bị giới hạn và thời gian xử lý kéo dài.
• So sánh với nghiên cứu tại Iran trên lò VVER-1000 Bushehr của Fadaei & Setayeshi (2008, 2013) sử dụng công cụ mã lai mạng nơ-ron kết hợp SA (LONSA) hoặc giải thuật đàn kiến (ACO): Các công trình này chưa giải quyết triệt để bài toán vùng hoạt nạp 30% MOX có độ bất đối xứng nơtron nhiệt nghiêm trọng và chưa chuẩn hóa kiểm định thống kê phi tham số trên số lượng mẫu chạy lặp lại đủ lớn.

Đó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 Khuếch tán Nơtron Đa nhóm trong môi trường bất đồng nhất hình học lục giác và Lý thuyết Tối ưu hóa Tổ hợp Không gian Rời rạc:

1. Mở rộng Lý thuyết Cân bằng Nơtron Vùng hoạt Phức hợp: Chứng minh rằng trong vùng hoạt tải hỗn hợp $\text{UO}2\text{-MOX}$, việc tối ưu hóa đồng thời $k{eff}$ BOC và làm phẳng phân bố thông lượng nhiệt không chỉ đơn thuần là bài toán phân bổ độ giàu, mà là sự điều tiết gradient dòng nơtron rò rỉ tại biên tiếp giáp giữa các bó MOX (có tiết diện hấp thụ nhiệt $\Sigma_a$ cao) và các bó $\text{UO}_2$ cháy thấp.
2. Xây dựng Mô hình Tối ưu hóa Lai ghép Tiến hóa - Tôi kim (Evolutionary Simulated Annealing Framework): Mở rộng lý thuyết nhiệt động học thống kê của Kirkpatrick bằng cách thay thế toán tử nhiễu loạn cục bộ (binary/ternary exchange) bằng toán tử tái tổ hợp di truyền (Crossover C1 và C2), cho phép hệ thống thực hiện các "bước nhảy lượng tử" qua các rào cản thế năng năng lượng cao trong không gian nghiệm.
3. Mở rộng Không gian Ánh xạ Tiến hóa Vi phân Rời rạc (Discrete Differential Evolution Theory): Thiết lập cơ chế ánh xạ liên tục - rời rạc thông qua giải thuật Chỉ số Vị trí Tương đối (Relative Position Indexing - RPI), bảo toàn tính toàn vẹn của các véc-tơ cá thể hoán vị không lặp lại, cho phép thuật toán SHADE tự thích ứng tham số $F$ và $CR$ trên không gian nghiệm nguyên rời rạc của cấu hình lò phản ứng.

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

Khung phân tích của luận án tích hợp 3 trụ cột lý thuyết: Lý thuyết Chuyển vận/Khuếch tán Nơtron (Boltzmann Transport/Fick's Diffusion), Lý thuyết Tối ưu hóa Tiến hóa Tự thích ứng (Self-Adaptive Evolutionary Computation), và Thống kê Phi tham số Nâng cao (Non-parametric Inferential Statistics).

Mô hình khái niệm được xác lập thông qua hai hàm thích nghi (Fitness Functions):

• Hàm thích nghi $F_1$ (Đánh giá khả năng tái tạo cấu hình tham chiếu): $$F_1 = \sqrt{\sum_{i=1}^{N_{FA}} (LP_{trial}(i) - LP_{ref}(i))^2}$$
• Hàm thích nghi $F_2$ (Tối ưu hóa đa mục tiêu năng lượng và an toàn): $$F_2 = w_k \cdot k_{eff} - w_p \cdot PPF - w_f \cdot Flatness$$ Trong đó, $w_k, w_p, w_f$ là các trọng số chuẩn hóa; $PPF = \max(P_i / \bar{P})$ là hệ số đỉnh công suất xuyên tâm; $Flatness = \sum_{i=1}^{N_{FA}} |P_i - \bar{P}| / N_{FA}$ đại diện cho độ đồng đều công suất.

Điều kiện biên (Boundary Conditions) được xác định chặt chẽ: Đối xứng quay góc $60^\circ$ ($1/6$ core symmetry), điều kiện biên phản xạ chân không tại vỏ lò phản ứng (vacuum boundary outer vessel), và các ràng buộc kỹ thuật nghiêm ngặt về giới hạn an toàn nhiệt $PPF < 1.35$.

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

Thiết kế nghiên cứu

Nghiên cứu được định vị theo trường phái Triết học Thực chứng Định lượng (Positivism / Quantitative Paradigm), dựa trên các định luật vật lý cơ bản và toán học rời rạc. Thiết kế nghiên cứu theo quy trình đa tầng (Multi-level design):

• Cấp độ vi mô (Lattice Level): Tính toán thông số tế bào mạng và tiết diện vĩ mô 4 nhóm năng lượng bằng phương pháp xác suất va chạm (PIJ module thuộc hệ mã SRAC).
• Cấp độ trung mô (Assembly Level): Mô hình hóa cấu trúc 331 thanh trong mỗi bó lục giác thành 24 ô mạng tam giác đồng nhất.
• Cấp độ vĩ mô (Core Level): Giải hệ phương trình khuếch tán 2D trên toàn bộ $1/6$ vùng hoạt thông qua mã tự phát triển LPO-V, tích hợp với các thuật toán tối ưu hóa ESA và discrete SHADE.

Quy trình nghiên cứu rigorous

Quy trình thu thập dữ liệu và mô phỏng số tuân thủ tiêu chuẩn Benchmark Quốc tế của OECD/NEA dành cho lò VVER-1000 nạp 30% MOX:

1. Chuẩn bị hằng số nhóm: Dữ liệu hạt nhân cơ bản từ thư viện ENDF/B-VII.0 với 107 nhóm năng lượng (62 nhóm nhanh, 45 nhóm nhiệt) được tính toán chi phối va chạm không gian thực 400 vùng hình học, sau đó rút gọn (collapsed) về cấu trúc 4 nhóm (3 nhóm nhanh: $10\text{ MeV} \to 0.821\text{ MeV} \to 5.53\text{ keV} \to 0.625\text{ eV}$; 1 nhóm nhiệt: $< 0.625\text{ eV}$).
2. Xây dựng và kiểm chứng mã tính LPO-V: Giải phương trình khuếch tán đa nhóm dạng ma trận: $$-\nabla \cdot (D_g \nabla \phi_g) + \Sigma_{r,g} \phi_g = \sum_{h \neq 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 không gian bằng phương pháp sai phân hữu hạn (FDM) trên lưới tam giác 2 chiều, sử dụng hệ số ghép nối dòng nơtron mặt biên: $$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}) = d_{g,i,k} (\phi_{g,i} - \phi_{g,k})$$ Hệ phương trình đại số tuyến tính lớn được giải bằng thuật toán lặp SOR với tiêu chí hội tụ thông lượng $\epsilon_\phi \le 10^{-6}$ và giá trị riêng $\epsilon_k \le 10^{-5}$.
3. Phát triển thuật toán tối ưu hóa:
   • ESA: Khởi tạo từ cấu hình cơ sở, sinh cấu hình thử nghiệm bằng toán tử Crossover C1 (hoán đổi ngẫu nhiên 2 bó giữa hai cha mẹ) hoặc C2 (hoán đổi cả một khối bó phụ) kết hợp toán tử đột biến đảo vị trí ngẫu nhiên, chấp nhận nghiệm mới theo tiêu chuẩn Metropolis: $P = \exp(-\Delta F / T)$.
   • Discrete SHADE: Sử dụng cơ chế nhớ lịch sử $H = 6$, kích thước quần thể $N_P = 20$. Tại mỗi thế hệ, tham số $F_i$ và $CR_i$ được sinh ngẫu nhiên từ phân bố Cauchy $C(M_F, 0.1)$ và phân bố chuẩn $N(M_{CR}, 0.1)$, sau đó cập nhật lại bộ nhớ thành công theo trung bình Lehmer ($mean_L$) và trung bình trọng số số học. Véc-tơ liên tục sau đột biến và lai ghép được chuyển sang véc-tơ hoán vị nguyên bằng giải thuật RPI.

+-------------------------------------------------------------------------+
|                  QUY TRÌNH TỐI ƯU HÓA TỔ HỢP LPO-V                      |
+-------------------------------------------------------------------------+
|  [Thư viện ENDF/B-VII.0] -> [SRAC/PIJ (107 nhóm)] -> [Rút gọn 4 nhóm]  |
|                                     │                                   |
|                                     ▼                                   |
|  [Khởi tạo Quần thể LP / Cấu hình Cơ sở (1/6 Đối xứng - 28 Bó VVER)]   |
|                                     │                                   |
|      ┌──────────────────────────────┴─────────────────────────────┐     |
|      ▼                                                            ▼     |
|  [Thuật toán ESA]                                       [Thuật toán SHADE] 
|  - Crossover C1 / C2                                    - Đột biến thích ứng (F)
|  - Đột biến hoán vị                                     - Trao đổi chéo (CR)
|  - Làm nguội Boltzmann                                  - Ánh xạ RPI rời rạc
|      └──────────────────────────────┬─────────────────────────────┘     |
|                                     │                                   |
|                                     ▼                                   |
|             [MÃ MÔ PHỎNG VẬT LÝ VÙNG HOẠT LPO-V]                        |
|             - Sai phân hữu hạn (FDM) Lưới tam giác                      |
|             - Bộ giải ma trận Siêu thư giãn (SOR)                       |
|             - Tính toán: k_eff, PPF, Phân bố công suất                  |
|                                     │                                   |
|                                     ▼                                   |
|             [Đánh giá Hàm Mục tiêu Đa Biến (F_2)]                       |
|                                     │                                   |
|                                     ▼                                   |
|        [Kiểm định Thống kê Phi tham số Mann-Whitney U Test]            |
+-------------------------------------------------------------------------+

Data và phân tích

Đặc tính mẫu phân tích bao gồm cấu hình $1/6$ vùng hoạt đối xứng với 28 vị trí bó nhiên liệu, phân bổ thành 7 nhóm nhiên liệu có mức độ cháy và thành phần đồng vị khác biệt:

• Nhóm $U_0$ (5 bó $\text{UO}2$ mới, độ cháy $0\text{ MWd/kg}$, $k{\infty} = 1.30604$)
• Nhóm $U_{15}$ (4 bó $\text{UO}2$, độ cháy $15\text{ MWd/kg}$, $k{\infty} = 1.13524$)
• Nhóm $U_{32}$ (5 bó $\text{UO}2$, độ cháy $32\text{ MWd/kg}$, $k{\infty} = 1.00293$)
• Nhóm $U_{40}$ (5 bó $\text{UO}2$, độ cháy $40\text{ MWd/kg}$, $k{\infty} = 0.94726$)
• Nhóm $M_0$ (3 bó MOX mới, độ cháy $0\text{ MWd/kg}$, $k_{\infty} = 1.18956$)
• Nhóm $M_{17}$ (4 bó MOX, độ cháy $17\text{ MWd/kg}$, $k_{\infty} = 1.05436$)
• Nhóm $M_{33}$ (2 bó MOX, độ cháy $33\text{ MWd/kg}$, $k_{\infty} = 0.94314$)

Phân tích thống kê suy luận sử dụng kiểm định phi tham số Mann-Whitney U Test trên 50 lần chạy độc lập để so sánh các phân phối giá trị thích nghi cực đại giữa các thuật toán. Mức ý nghĩa thống kê được thiết lập tại $\alpha = 0.05$ (với $p$-value $< 0.05$ biểu thị sự vượt trội có ý nghĩa thống kê của giải thuật đề xuất).

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

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

1. Hiệu năng và độ chính xác vượt trội của mã tính LPO-V: Kiểm chứng trên 5 trạng thái vận hành chuẩn (S1–S5) của VVER-1000 MOX Benchmark chứng minh sai số $k_{eff}$ giữa LPO-V và mã chuẩn Monte Carlo MCNP4c là cực kỳ nhỏ: Trạng thái làm việc S1 đạt $k_{eff} = 1.04169$ (LPO-V) so với $1.04018 \pm 0.00030$ (MCNP4c), độ lệch chỉ $151\text{ pcm}$; trạng thái không Bo S4 đạt $k_{eff} = 1.10304$ so với $1.10170 \pm 0.00030$, độ lệch $134\text{ pcm}$. Đặc biệt, LPO-V cho tốc độ tính toán trung bình đạt 10 LPs/giây, nhanh hơn đáng kể so với mô-đun CITATION khi giải trên cùng hệ lưới tam giác.
2. Đột phá về nâng cao tính kinh tế và an toàn của cấu hình nạp tải tối ưu: Cấu hình nạp tải tối ưu tìm được bởi cả hai thuật toán ESA và discrete SHADE đạt kết quả hội tụ đồng nhất, cải thiện vượt bậc các chỉ số vận hành so với cấu hình tham chiếu của OECD/NEA:
   • Hệ số nhân hiệu dụng đầu chu trình ($k_{eff}$) tăng thêm $1580\text{ pcm}$ ($0.01580\ \Delta k$), tương đương với việc kéo dài đáng kể thời gian vận hành chu trình công suất định mức ("The results show that the keff of the optimal LP is greater than that of the reference core by about 1580 pcm").
   • Hệ số đỉnh công suất xuyên tâm ($PPF$) giảm $2.4%$, từ $1.346$ ở cấu hình tham chiếu xuống $1.314$ ở cấu hình tối ưu ("Whereas, the radial power peaking factor (PPF) of the optimal LP is about 2.4% smaller than that of the reference core"), gia tăng biên dự trữ an toàn chống quá nhiệt cục bộ (Thermal Margin).
3. Hiện tượng dịch chuyển phân bố nhiên liệu MOX và bó nhiên liệu cháy: Phát hiện quy luật bố trí vật lý phi trực giác: Để hạn chế rò rỉ nơtron và kiềm chế đỉnh công suất do hiệu ứng biên của các bó MOX, cấu hình tối ưu đã bố trí các bó MOX độ cháy cao ($M_{33}$) xen kẽ với các bó $\text{UO}_2$ mới ($U_0$) tại khu vực vành trong, đồng thời đẩy các bó $\text{UO}2$ độ cháy sâu ($U{40}$) ra sát biên ngoài giáp vách phản xạ (low-leakage loading pattern strategy).
4. Bằng chứng thống kê về tính ưu việt của ESA và discrete SHADE: Kết quả kiểm định Mann-Whitney U Test khẳng định sự vượt trội có ý nghĩa thống kê ($p < 0.05$) của ESA-C1 và SHADE so với SA nguyên bản và ASA. Cụ thể, trong bài toán tái tạo cấu hình tham chiếu ($F_1$), SA chỉ đạt tỷ lệ tái tạo thành công $12%$, ASA đạt $44%$, trong khi ESA-C1 đạt $96%$. Đối với hàm mục tiêu đa biến $F_2$, SHADE và ESA đạt giá trị fitness trung bình cao hơn và độ lệch chuẩn nhỏ hơn một bậc độ lớn so với SA và ASA.

Implications đa chiều

• Về mặt lý thuyết: Cung cấp cơ sở lý luận vững chắc về sự tương tác giữa phổ nơtron nhiệt của nhiên liệu MOX và hình học lưới tam giác VVER, khẳng định tính khả thi của việc tối ưu hóa đa mục tiêu xung đột mà không cần hy sinh biên độ an toàn nhiệt.
• Về mặt phương pháp luận: Khẳng định thuật toán discrete SHADE kết hợp RPI là công cụ metaheuristic tổng quát, có khả năng chuyển giao và ứng dụng trực tiếp cho các bài toán tối ưu hóa tổ hợp hoán vị số nguyên phức tạp khác trong kỹ thuật lò phản ứng.
• Về mặt thực tiễn và chính sách: 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 cắt giảm chi phí chu trình nhiên liệu hàng triệu USD mỗi năm thông qua việc nâng cao độ sâu cháy và kéo dài chu trình nhiên liệu, đồng thời hỗ trợ các cơ quan pháp quy hạt nhân thiết lập quy trình thẩm định an toàn vùng hoạt độc lập.

Limitations và Future Research

Nghiên cứu thừa nhận một số giới hạn phương pháp luận và phạm vi ứng dụng:

1. Mô hình hóa 2 chiều ở trạng thái đầu chu trình (2D BOC): Luận án tập trung tính toán tối ưu tại thời điểm BOC mà chưa tích hợp đầy đủ quá trình cháy nhiên liệu theo thời gian 3 chiều (3D depletion calculation) kết hợp phản hồi nhiệt - thủy lực liên kết (Thermal-Hydraulic Feedback - $T_{fuel}, T_{mod}$ distribution).
2. Phụ thuộc vào bộ hằng số nhóm rút gọn tĩnh: Tiết diện vĩ mô 4 nhóm được sinh ra từ các tính toán tế bào mạng đơn lẻ tại trạng thái biên đẳng hướng, chưa xét đến hiệu ứng phổ giao thoa giữa các bó lân cận (cross-assembly spectral history effects).
3. Giới hạn hình học đối xứng quay $1/6$: Giả định đối xứng hoàn hảo loại trừ khả năng khảo sát các cấu hình nạp tải bất đối xứng hoàn toàn trong các trường hợp có bó nhiên liệu bị khuyết tật hoặc kiểm soát thanh điều khiển cục bộ.

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

• Mở rộng mã tính LPO-V sang mô hình 3D đa nhóm kết hợp tính toán biến đổi độ cháy từng bước thời gian (burnup steps).
• Ghép nối mã LPO-V với mô-đun nhiệt - thủy lực kênh phụ (sub-channel thermal-hydraulics code) để đánh giá tỷ số truyền nhiệt tới hạn (DNBR) trong thời gian thực của quá trình tìm kiếm.
• Phát triển các thuật toán tiến hóa đa mục tiêu Pareto đa cá thể (Multi-Objective SHADE / MO-DE) để cung cấp tập nghiệm Pareto tường minh giữa $k_{eff}$ và $PPF$.
• Ứng dụng công nghệ tính toán song song trên nền tảng GPU để nâng tốc độ đánh giá lên hàng trăm LPs/giây.

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

• Ảnh hưởng học thuật: Thiết lập chuẩn mực nghiên cứu mới trong lĩnh vực vật lý lò phản ứng VVER tại Việt Nam; công bố các bài báo khoa học chất lượng cao trên các tạp chí chuyên ngành quốc tế ISI/Scopus (như Annals of Nuclear Energy, Nuclear Engineering and Design).
• Chuyển đổi công nghiệp: Cung cấp bộ công cụ tính toán nạp đảo nhiên liệu độc lập có thể chuyển giao cho các trung tâm nghiên cứu hạt nhân, cơ sở đào tạo chuyên sâu và các đơn vị thiết kế lò phản ứng.
• Tác động quốc tế: Đóng góp dữ liệu số và giải pháp tối ưu cho cộng đồng nghiên cứu quốc tế thuộc Cơ quan Năng lượng Nguyên tử Quốc tế (IAEA) và Tổ chức Hợp tác và Phát triển Kinh tế (OECD/NEA) trong chương trình quản lý nhiên liệu MOX nâng cao.

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

• Nghiên cứu sinh và học viên sau đại học: Tiếp cận phương pháp luận mô hình hóa vật lý nơtron lưới tam giác FDM và kỹ thuật lập trình giải thuật metaheuristic lai ghép tiên tiến.
• Các nhà nghiên cứu và kỹ sư vật lý lò: Sử dụng mã tính LPO-V như một công cụ đánh giá nhanh, hiệu năng cao để sơ loại hàng vạn cấu hình nạp tải trước khi chạy kiểm chứng chi tiết bằng các mã Monte Carlo nặng.
• Cơ quan quản lý an toàn bức xạ và hạt nhân: Nắm vững công cụ độc lập để đánh giá và thẩm tra các báo cáo phân tích an toàn vùng hoạt do nhà thầu cung cấp.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào? Đóng góp độc đáo nhất là việc phát triển khung tích hợp tối ưu hóa tiến hóa vi phân thích ứng rời rạc (Discrete SHADE Framework), mở rộng Lý thuyết Tiến hóa Vi phân từ không gian liên tục sang không gian tổ hợp rời rạc thông qua thuật toán RPI, đồng thời giải quyết thành công bài toán cân bằng nơtron đa mục tiêu trong vùng hoạt VVER-1000 tải 30% MOX.
2. Cải tiến phương pháp luận thể hiện thế nào khi so sánh với các nghiên cứu trước? So với phương pháp Tôi kim cổ điển (SA) của Kirkpatrick (1983) và Tôi kim Thích ứng (ASA) của Ingber (1993) vốn chỉ dùng hoán đổi nhị phân, phương pháp ESA của luận án tích hợp toán tử trao đổi chéo C1/C2 giúp tỷ lệ tìm thấy cấu hình tối ưu tăng từ $12%$ lên $96%$. So với GA truyền thống (Chapot et al., 1999), discrete SHADE loại bỏ hoàn toàn việc phải hiệu chỉnh thủ công các tham số điều khiển ($F, CR$) nhờ cơ chế thích ứng lịch sử thành công.
3. Phát hiện nào gây bất ngờ nhất và dữ liệu nào chứng minh điều đó? Phát hiện bất ngờ nhất là việc đạt được mức tăng $k_{eff}$ rất lớn ($+1580\text{ pcm}$) nhưng không hề làm tăng $PPF$ mà ngược lại còn kéo giảm $PPF$ xuống $2.4%$ ($1.314$ so với $1.346$). Dữ liệu phân bố công suất 2D chứng minh rằng sự kết hợp vi vị trí giữa các bó MOX giàu Plutonium và các bó $\text{UO}_2$ cháy vừa đã tạo ra vùng đệm thông lượng nhiệt tự nhiên làm phẳng toàn bộ biên dạng công suất xuyên tâm.
4. Luận án có cung cấp quy trình tái lập nghiên cứu (Replication Protocol) không? Luận án cung cấp đầy đủ thông số hình học, thành phần đồng vị chi tiết từng vùng tế bào mạng (Phụ lục A), bộ tiết diện vĩ mô 4 nhóm năng lượng (Phụ lục B), cấu trúc thuật toán chi tiết và lưu đồ thực thi của mã LPO-V, bảo đảm khả năng tái lập hoàn toàn trên bất kỳ hệ thống tính toán nào.
5. Kế hoạch nghiên cứu 10 năm được phác thảo như thế nào? Lộ trình 10 năm tập trung vào việc nâng cấp mã LPO-V thành hệ thống phần mềm mô phỏng 3D thời gian thực tích hợp ghép nối đa vật lý (Neutronics - Thermal Hydraulics - Fuel Performance) phục vụ thiết kế các thế hệ lò VVER-1200, VVER-TOI và lò phản ứng thế hệ IV làm mát bằng nước siêu tới hạn (VVER-SCWR).

Kết luận

1. Xây dựng và hoàn thiện thành công mã mô phỏng vật lý vùng hoạt LPO-V giải hệ phương trình khuếch tán nơtron 4 nhóm trên lưới tam giác FDM-SOR, đạt độ chính xác cao (sai số $k_{eff} < 150\text{ pcm}$ so với MCNP4c) và tốc độ tính toán vượt trội phục vụ bài toán tối ưu hóa.
2. Phát triển thành công hai giải thuật tối ưu hóa tiên tiến: Thuật toán Tôi kim mô phỏng tiến hóa (ESA) với toán tử trao đổi chéo di truyền và Thuật toán Tiến hóa vi phân thích ứng dựa trên lịch sử thành công rời rạc (Discrete SHADE) sử dụng kỹ thuật RPI.
3. Xác lập cấu hình nạp tải tối ưu cho lò VVER-1000 nạp 30% MOX, nâng cao $k_{eff}$ đầu chu trình thêm $1580\text{ pcm}$ và làm giảm hệ số đỉnh công suất xuyên tâm $PPF$ $2.4%$ so với cấu hình tham chiếu quốc tế OECD/NEA.
4. Chứng minh bằng thực nghiệm số và kiểm định phi tham số Mann-Whitney U Test ($p < 0.05$) về tính ưu việt vượt trội của ESA và discrete SHADE so với các thuật toán truyền thống SA và ASA.
5. Mở ra 3 hướng nghiên cứu mới: Mô phỏng cháy 3D ghép nối nhiệt - thủy lực liên kết, Tối ưu hóa đa mục tiêu Pareto đa cá thể, và Ứng dụng thuật toán tối ưu rời rạc cho các thế hệ lò nước nhẹ tiên tiến và lò phản ứng nghiên cứu hạt nhân đa mục tiêu.
6. Thiết lập di sản khoa học thực chứng với bộ công cụ tính toán mã nguồn độc lập, nâng cao năng lực nghiên cứu làm chủ công nghệ tính toán an toàn hạt nhân chuyên sâu.