Tổng quan về luận án

Nghiên cứu quản lý nhiên liệu trong vùng hoạt (In-core Fuel Management - ICFM) đóng vai trò then chốt trong kỹ thuật hạt nhân, trực tiếp quyết định tính an toàn và hiệu quả kinh tế của các lò phản ứng năng lượng. Trong cơ cấu chi phí phát điện hạt nhân, chi phí chế tạo và quản lý nhiên liệu chiếm tỷ trọng đáng kể, khoảng 15% đến 25% tổng chi phí vận hành. Luận án tiến sĩ vật lý của tác giả Trần Việt Phú, dưới sự hướng dẫn của PGS. Trần Hoài Nam và GS. Akio Yamamoto tại Viện Năng lượng Nguyên tử Việt Nam, tập trung giải quyết bài toán tối ưu hóa mô hình nạp đảo nhiên liệu (Fuel Loading Pattern - LP) cho lò phản ứng nước áp lực thế hệ VVER-1000 nạp 30% nhiên liệu oxit hỗn hợp ($MOX - Mixed\ Oxide\ Fuel$).

Bối cảnh khoa học của nghiên cứu xuất phát từ thách thức 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 vùng hoạt chuẩn VVER-1000 gồm 163 bó nhiên liệu lục giác, nếu không xét tính đối xứng, số lượng cấu hình có thể lên tới $1.9 \times 10^{280}$. Ngay cả khi áp dụng điều kiện đối xứng hình học $1/6$ và tính đến sự tương đồng của các bó nhiên liệu, không gian tìm kiếm vẫn đạt mức $4.5 \times 10^{10}$ cấu hình. Các phương pháp tối ưu hóa cổ điển như Tìm kiếm trực tiếp (Direct Search - DS) hay Trao đổi nhị phân (Binary Exchange - BE) thường nhanh chóng rơi vào bẫy cực trị địa phương (local optima), trong khi thuật toán Tôi kim mô phỏng (Simulated Annealing - SA) truyền thống có tốc độ hội tụ rất chậm.

Luận án đặt ra ba 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 công cụ tính toán vật lý vùng hoạt chuyên dụng cho hình học lục giác của lò VVER với độ chính xác tương đương các mã chuẩn quốc tế nhưng có tốc độ thực thi vượt trội nhằm đáp ứng hàng chục nghìn lượt đánh giá cấu hình?
  2. RQ2: Cải tiến toán tử sinh nghiệm của thuật toán SA như thế nào để vừa duy trì khả năng thoát cực trị địa phương vừa tăng tốc độ hội tụ toàn cục?
  3. RQ3: Làm thế nào để chuyển đổi và áp dụng hiệu quả 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 (SHADE) – vốn thiết kế cho biến liên tục – sang bài toán biến số nguyên rời rạc của cấu hình nạp nhiên liệu hạt nhân?

Tương ứng với các câu hỏi là ba giả thuyết khoa học:

  • H1: Phương pháp sai phân hữu hạn (Finite Difference Method - FDM) giải phương trình khuếch tán đa nhóm trên lưới tam giác kết hợp kỹ thuật lặp quá thư giãn liên tiếp (Successive Over-Relaxation - SOR) sẽ đảm bảo độ chính xác tiệm cận mô phỏng Monte Carlo MCNP4c nhưng rút ngắn thời gian tính toán so với mô-đun CITATION.
  • H2: Tích hợp toán tử lai ghép (crossover) và đột biến (mutation) từ thuật toán di truyền vào thuật toán tôi kim (tạo thành giải thuật Evolutionary Simulated Annealing - ESA) sẽ nâng cao rõ rệt hiệu suất tìm kiếm so với SA và Adaptive Simulated Annealing (ASA).
  • H3: Cơ chế đánh chỉ số vị trí tương đối (Relative Position Indexing - RPI) cho phép thuật toán SHADE tối ưu hóa rời rạc hiệu quả mà không làm phá vỡ cấu trúc thích ứng tham số lịch sử ($F$ và $CR$).

Khung lý thuyết của luận án tích hợp chặt chẽ giữa Lý thuyết khuếch tán nơ-tron đa nhóm (Multi-group Neutron Diffusion Theory) và Lý thuyết siêu nghiệm phỏng sinh học thích ứng (Adaptive Metaheuristics). Đột phá định lượng của nghiên cứu thể hiện qua việc phát triển mã tính LPO-V và tìm ra cấu hình nạp tải tối ưu giúp tăng hệ số nhân hiệu dụng ($k_{eff}$) thêm 1580 pcm, đồng thời giảm hệ số đỉnh công suất hướng tâm ($PPF - Power\ Peaking\ Factor$) khoảng 2.4% so với cấu hình chuẩn của OECD/NEA. Phạm vi nghiên cứu bao quát bài toán chuẩn VVER-1000 MOX với 28 vị trí bó nhiên liệu trong góc đối xứng $1/6$, trải qua 5 trạng thái vận hành từ trạng thái làm việc có boron đến trạng thái nguội nồng độ boron cao.

Literature Review và Positioning

Lịch sử tối ưu hóa nạp đảo nhiên liệu lò phản ứng đã chứng kiến sự phát triển qua nhiều giai đoạn. Dòng nghiên cứu dựa trên không gian lân cận (Neighborhood-based) mở đầu với thuật toán SA của Kirkpatrick et al. (1983), được Parks (1996) và các cộng sự ứng dụng rộng rãi trong thiết kế vùng hoạt lò PWR nhờ khả năng chấp nhận nghiệm xấu hơn theo phân bố Boltzmann nhằm vượt qua rào cản thế năng. Tuy nhiên, Ingber (1993) chỉ ra rằng SA cổ điển đòi hỏi lịch trình hạ nhiệt quá chậm, dẫn đến chi phí tính toán khổng lồ. Để khắc phục, các biến thể như ASA (Szu & Hartley, 1987) đã đưa vào cơ chế tái sử dụng nghiệm tốt nhất ("return to the best") và ma trận xác suất chuyển trạng thái. Song song đó, kỹ thuật Tìm kiếm Tabu (Tabu Search - TS) của Glover (1986) được triển khai để ngăn ngừa vòng lặp lân cận.

Dòng nghiên cứu dựa trên quần thể (Population-based) phát triển mạnh với Giải thuật Di truyền (GA) của Goldberg (1989), được DeChaine & Feltus (1995) đưa vào ICFM. Sau đó, các phương pháp Tối ưu hóa bầy đàn (PSO) của Kennedy & Eberhart (1995), Tối ưu hóa đàn kiến (ACO) và Tiến hóa vi phân (Differential Evolution - DE) của Storn & Price (1997) lần lượt được thử nghiệm. Năm 2013, Tanabe & Fukunaga đề xuất thuật toán SHADE, tạo nên bước nhảy vọt trong tối ưu hóa liên tục nhờ cơ chế tự điều chỉnh hệ số đột biến $F$ và tỷ lệ lai ghép $CR$ dựa trên bộ nhớ lịch sử thành công $H$.

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 tìm kiếm lân cận (SA/ASA): Ưu tiên tính bảo toàn cấu trúc tổ hợp của bó nhiên liệu, không làm phát sinh nghiệm vi phạm ràng buộc số lượng bó, nhưng đánh đổi bằng tốc độ quét không gian giải pháp rất chậm.
  • Trường phái tiến hóa quần thể (GA/DE): Có khả năng đa dạng hóa tìm kiếm (exploration) xuất sắc nhưng gặp trở ngại lớn khi xử lý biến rời rạc. Các toán tử đột biến liên tục trong DE gốc thường phá vỡ tính nguyên vẹn của danh mục nhiên liệu đầu vào.

Luận án của tác giả Trần Việt Phú định vị chính xác tại giao điểm này, giải quyết triệt để hạn chế của cả hai trường phái: cải tiến SA bằng toán tử lai ghép tiến hóa để tạo ra ESA, đồng thời làm chủ cơ chế RPI để đưa thuật toán tiên tiến SHADE vào không gian rời rạc của lò VVER.

So với các công trình quốc tế tiêu biểu:

  1. Mã OPTIMAL và ATHENA (Cộng hòa Séc): Áp dụng phương pháp hoán đổi nhị phân và ngẫu nhiên cho lò VVER-440 Dukovany (Faber et al., 1994), nhưng dễ dừng lại ở cực trị cục bộ do thiếu cơ chế tiến hóa thích ứng.
  2. Hệ thống SCAM-W / CIGARO (Temelin NPP, Séc) và LONSA (Bushehr NPP, Iran): Tích hợp mạng nơ-ron và GA/SA (Sadeghi et al., 2013) nhưng phụ thuộc nặng nề vào việc tinh chỉnh siêu tham số thủ công. Discrete SHADE trong luận án này vượt trội hơn nhờ khả năng tự thích ứng tham số hoàn toàn tự động thông qua bộ nhớ lịch sử.

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

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

Luận án đóng góp quan trọng vào lý thuyết tính toán vật lý hạt nhân và lý thuyết tối ưu hóa tổ hợp thông qua việc mở rộng các nguyên lý thích ứng tham số trong môi trường rời rạc. Nghiên cứu đã hình thành hai mệnh đề lý thuyết cốt lõi:

  • Proposition 1 (Về tính bảo toàn và đa dạng hóa trong ESA): Việc tích hợp toán tử trao đổi chéo ngẫu nhiên (Crossover C1 - hoán đổi 2 vị trí bất kỳ) và trao đổi khối (Crossover C2 - hoán đổi cụm vị trí liên kết) giữa hai cấu hình bố trí cơ sở ($LP_{base}$) bảo toàn tuyệt đối tính khả thi của danh mục nhiên liệu, đồng thời tạo ra bước nhảy trạng thái đủ lớn để phá vỡ các bẫy thế năng cục bộ mà các phép hoán đổi đơn điểm cổ điển không thể vượt qua.
  • Proposition 2 (Về tính bất biến thứ tự trong discrete SHADE): Ánh xạ không gian vector liên tục $\mathbb{R}^D$ sang không gian hoán vị nguyên $\mathbb{Z}^D$ thông qua hàm thứ tự vị trí tương đối (RPI) bảo toàn nguyên vẹn động lực học thích ứng của các phân bố Cauchy (cho hệ số đột biến $F$) và phân bố Chuẩn (cho tỷ lệ lai ghép $CR$), duy trì trọn vẹn sức mạnh tìm kiếm của SHADE trên miền không gian cấu hình vùng hoạt.

Nghiên cứu tạo ra sự chuyển dịch mô thức từ việc tối ưu hóa dựa trên các quy tắc kinh nghiệm cố định (fixed heuristics) sang tối ưu hóa tự học thích ứng theo lịch sử tiến hóa (history-based adaptive search) trong bài toán quản lý nhiên liệu hạt nhân.

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

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

  1. Lý thuyết trường nơ-tron tĩnh (Steady-state Neutron Transport & Multi-group Diffusion): Mô hình hóa phân bố thông lượng nơ-tron và công suất thông qua hệ phương trình đạo hàm riêng elliptic giải trên lưới tam giác.
  2. Lý thuyết siêu nghiệm tiến hóa phỏng sinh học: Điều khiển quá trình tìm kiếm cấu hình nạp tải thông qua hàm mục tiêu đa biến (Fitness Function).
  3. Lý thuyết kiểm định thống kê phi tham số (Non-parametric Statistics): Ứng dụng kiểm định Mann-Whitney U nhằm lượng hóa sự khác biệt có ý nghĩa thống kê giữa các thuật toán tối ưu ngẫu nhiên qua 50 lần chạy độc lập.

Khung phân tích thiết lập hai hàm mục tiêu chuẩn hóa:

  • Hàm mục tiêu $F_1$ (Kiểm tra khả năng tái lập cấu hình tham chiếu): $$F_1 = \sum_{i=1}^{N_{ass}} |k_{\infty, i}^{calc} - k_{\infty, i}^{ref}|$$
  • Hàm mục tiêu $F_2$ (Tối ưu hóa kinh tế và an toàn vận hành): $$F_2 = w_k \cdot (k_{eff} - 1.0) - w_p \cdot \max(0, PPF - PPF_{limit}) - w_f \cdot Flatness$$ Trong đó, các điều kiện biên được xác lập nghiêm ngặt: đối xứng $60^\circ$ góc quay, dòng nơ-tron triệt tiêu tại bề mặt ngoài vỏ lò phản ứng (vacuum boundary condition), bảo toàn số lượng 19 bó $UO_2$ và 9 bó MOX với 7 cấp độ cháy khác nhau.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ lập trường nhận thức luận thực chứng (Positivism), kết hợp mô hình hóa giải tích số học tất định với các thuật toán tối ưu hóa ngẫu nhiên (stochastic optimization). Thiết kế nghiên cứu đa tầng bao gồm ba cấp độ tích hợp:

  • Cấp độ mạng vi mô (Lattice level): Tính toán tiết diện vi mô và vĩ mô 107 nhóm năng lượng bằng mô-đun PIJ của mã chuẩn SRAC (sử dụng phương pháp xác suất va chạm - Collision Probability Method) trên cơ sở thư viện dữ liệu hạt nhân đánh giá ENDF/B-VII.0. Cấu trúc 107 nhóm (62 nhóm nhanh, 45 nhóm nhiệt) sau đó được rút gọn (collapsed) về 4 nhóm năng lượng (3 nhóm nhanh, 1 nhóm nhiệt).
  • Cấp độ vùng hoạt (Core level): Phát triển mã nguồn LPO-V giải hệ phương trình khuếch tán 4 nhóm năng lượng ở chế độ dừng trên lưới tam giác 2D: $$-\nabla \cdot (D_g \nabla \phi_g) + \Sigma_{r,g} \phi_g = S_g$$ Với nguồn nơ-tron $S_g$ bao gồm nơ-tron tán xạ từ các nhóm khác và nơ-tron phân hạch chia cho hệ số $k_{eff}$.
  • Cấp độ tối ưu hóa (Optimization level): Ghép nối trực tiếp mã LPO-V với hai thuật toán mới ESA và discrete SHADE.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu triển khai theo các bước chuẩn hóa nghiêm ngặt:

  1. Không gian hóa hình học: Mỗi bó nhiên liệu lục giác được chia nhỏ thành 24 phần tử lưới tam giác đồng nhất.
  2. Rời rạc hóa không gian: Ứng dụng FDM thiết lập dòng nơ-tron qua mặt phân cách $A_{i,k}$ giữa hai lưới tam giác lân cận $i$ và $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}) = d_{g,i,k} (\phi_{g,i} - \phi_{g,k})$$
  3. Giải hệ đại số tuyến tính: Ma trận thưa kích thước lớn được giải bằng thuật toán SOR với tiêu chí hội tụ thông lượng và nguồn phân hạch đạt sai số $\epsilon \le 10^{-6}$.
  4. Kiểm chuẩn mã tính (Code Verification): LPO-V được kiểm chứng chéo đối chuẩn với mã Monte Carlo MCNP4c và mã khuếch tán CITATION qua 5 trạng thái chuẩn từ S1 đến S5 của OECD/NEA.

Data và phân tích

Dữ liệu đầu vào sử dụng bộ thông số chuẩn của bài toán VVER-1000 MOX Core Benchmark:

  • 163 bó nhiên liệu, chiều cao vùng hoạt 3.55 m, công suất nhiệt 3000 MWt.
  • Góc đối xứng $1/6$ gồm 28 vị trí: Bó $UO_2$ chia thành 4 mức độ cháy (0, 15, 32, 40 MWd/kg) với các loại làm giàu $3.7%$ và $4.2%$ U-235 cùng thanh hấp thụ $Gd_2O_3$; Bó MOX chia thành 3 mức độ cháy (0, 17, 33 MWd/kg) với hàm lượng Pu phân hạch $2.4%$, $2.7%$ và $3.6%$.

Phân tích tối ưu hóa tiến hành khảo sát đa thông số:

  • Khảo sát trọng số hàm mục tiêu xác định bộ thông số tối ưu: $w_p = 2.5$, $w_f = 0.5$.
  • Khảo sát siêu tham số thuật toán SHADE: xác định kích thước quần thể $N_P = 20$, dung lượng bộ nhớ lịch sử $H = 6$.
  • Phân tích thống kê: Thực hiện 50 lần chạy độc lập cho mỗi cấu hình thuật toán (SA, ASA, ESA, DE, SHADE). Kiểm định Mann-Whitney U test được áp dụng ở mức ý nghĩa $\alpha = 0.05$ để chứng minh tính vượt trội về mặt thống kê ($p < 0.05$).

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

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

  1. Hiệu năng vượt trội của mã LPO-V: Kết quả tính toán $k_{eff}$ và phân bố công suất của LPO-V sai lệch dưới $0.2%$ so với kết quả tham chiếu MCNP4c trên toàn bộ các trạng thái S1–S5, trong khi tốc độ tính toán nhanh hơn đáng kể so với mô-đun CITATION, cho phép đánh giá từ $10^4$ đến $10^5$ cấu hình trong thời gian khả thi.
  2. Nâng cao vượt bậc hệ số nhân hiệu dụng: Cấu hình tối ưu do ESA và discrete SHADE tìm được đạt mức tăng $k_{eff}$ xấp xỉ 1580 pcm so với cấu hình tham chiếu ban đầu. Điều này đồng nghĩa với việc kéo dài đáng kể độ dài chu kỳ nạp nhiên liệu (cycle length).
  3. Cải thiện biên độ an toàn nhiệt thủy lực: Hệ số đỉnh công suất hướng tâm ($PPF$) của cấu hình tối ưu đạt mức giảm $2.4%$ so với lõi tham chiếu ($PPF$ giảm từ giá trị tham chiếu xuống mức an toàn tuyệt đối dưới ngưỡng giới hạn 1.35), giúp phân bố công suất vùng hoạt phẳng hơn đáng kể.
  4. Sự tương đương tối ưu giữa ESA và discrete SHADE: Cả hai phương pháp phát triển mới đều tìm ra cấu hình nạp tải tối ưu đồng nhất, vượt qua hoàn toàn các thuật toán SA và ASA cổ điển về cả tốc độ hội tụ lẫn giá trị hàm thích nghi tốt nhất.
  5. Ý nghĩa thống kê được kiểm chứng: Phép kiểm định Mann-Whitney U test với giá trị $p < 0.05$ khẳng định sự vượt trội của ESA và discrete SHADE so với SA và ASA không phải do ngẫu nhiên mà xuất phát từ ưu thế cấu trúc giải thuật.

Implications đa chiều

  • Về mặt học thuật: Khẳng định tính khả thi và hiệu quả vượt trội của việc số hóa thuật toán vi phân thích ứng (SHADE) cho các bài toán tối ưu tổ hợp rời rạc trong vật lý hạt nhân; làm phong phú kho tàng phương pháp luận ICFM thế giới.
  • Về mặt phương pháp: Bộ công cụ LPO-V cung cấp một nền tảng tính toán mã nguồn mở, độc lập, có khả năng mở rộng sang các hệ thống hình học phức tạp khác có cấu trúc thanh nhiên liệu tam giác.
  • Về mặt thực tiễn công nghiệp: Tăng $k_{eff}$ thêm 1580 pcm cho phép các nhà máy điện hạt nhân VVER-1000 tiết kiệm hàng triệu USD chi phí nhiên liệu mỗi chu kỳ thông qua việc khai thác tối đa độ cháy của các bó nhiên liệu cũ và giảm khối lượng nhiên liệu làm giàu nạp mới.
  • Về mặt an toàn và chính sách: Phân bố công suất phẳng hơn giúp giảm thiểu rủi ro quá nhiệt cục bộ tại các thanh nhiên liệu, gia tăng giới hạn an toàn khi vận hành ở công suất định mức, hỗ trợ các cơ quan pháp quy hạt nhân xây dựng tiêu chuẩn thẩm định an toàn vùng hoạt.

Limitations và Future Research

Nghiên cứu thừa nhận một số giới hạn khoa học mang tính điều kiện biên:

  1. Giới hạn chiều không gian: Mô hình tính toán LPO-V hiện tại được xây dựng trên không gian 2D, chưa phản ánh đầy đủ hiệu ứng bất đối xứng phân bố công suất theo phương trục đứng (axial power distribution) và hiện tượng biến dạng thông lượng 3D khi có sự chuyển dịch của các chùm thanh điều khiển.
  2. Đơn giản hóa quá trình cháy: Đánh giá tối ưu hóa hiện tập trung tại thời điểm bắt đầu chu kỳ (Beginning of Cycle - BOC) mà chưa ghép nối tính toán độ cháy phụ thuộc thời gian (burnup calculation) trong từng bước lặp tối ưu nhằm tiết kiệm thời gian tính toán.
  3. Phản hồi nhiệt thủy lực: Mô hình chưa xét đến sự phản hồi phi tuyến tức thời của nhiệt độ nhiên liệu (hiệu ứng Doppler) và mật độ chất làm chậm trong quá trình tính toán phân bố thông lượng.

Chương trình nghiên cứu tương lai định hướng 5 mũi nhọn:

  • Nâng cấp mã LPO-V lên mô hình 3D đa nhóm giải tích số mạng tam giác - lục giác.
  • Ghép nối tính toán song song đa vật lý (Neutronics - Thermal Hydraulics coupled calculations).
  • Tích hợp mô hình tính toán suy giảm đồng vị (depletion chain kinetics) trực tiếp vào vòng lặp tối ưu hóa toàn chu kỳ.
  • Ứng dụng các biến thể tiên tiến hơn của DE như L-SHADE hoặc SHADE-cnEpSin.
  • Mở rộng ứng dụng công cụ tối ưu hóa cho các thế hệ lò phản ứng mới như VVER-1200 (AES-2006), VVER-TOI và lò phản ứng làm nguội bằng nước siêu tới hạn VVER-SCWR (Thế hệ IV).

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

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

  • Tác động học thuật: Các bài báo khoa học công bố từ luận án trên các tạp chí quốc tế chuyên ngành hạt nhân uy tín (thuộc danh mục ISI/Scopus) mở ra hướng tiếp cận mới trong việc ứng dụng thuật toán tiến hóa thích ứng cho lò phản ứng nước áp lực kiểu Nga.
  • Chuyển giao công nghệ và công nghiệp hạt nhân: Cung cấp cho các kỹ sư vận hành vùng hoạt lò VVER một công cụ tính toán độc lập, tin cậy để thẩm định các kế hoạch đảo nhiên liệu do nhà cung cấp thiết kế.
  • An ninh năng lượng và môi trường: Nâng cao hiệu suất cháy của nhiên liệu MOX góp phần quan trọng vào việc tái chế Plutonium từ nhiên liệu đã qua sử dụng, giảm thiểu chất thải phóng xạ mức độ cao và thúc đẩy phát triển năng lượng hạt nhân bền vững.

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

  • Nghiên cứu sinh và học giả chuyên ngành Vật lý Hạt nhân: Tiếp cận phương pháp luận chặt chẽ trong việc giải số phương trình khuếch tán trên lưới tam giác và kỹ thuật số hóa thuật toán tiến hóa.
  • Chuyên gia thiết kế và quản lý vùng hoạt (R&D Nuclear Engineers): Sử dụng trực tiếp cấu trúc thuật toán ESA và discrete SHADE để tối ưu hóa chi phí vận hành cho các lò phản ứng thương mại.
  • Cơ quan Pháp quy và An toàn Hạt nhân: Có thêm công cụ tính toán độc lập để kiểm chứng các thông số an toàn giới hạn ($PPF$, $k_{eff}$) trong hồ sơ xin cấp phép nạp đảo nhiên liệ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ì? Đóng góp độc đáo nhất là việc phát triển và hoàn thiện thành công cơ chế chuyển đổi không gian biến liên tục sang rời rạc thông qua Relative Position Indexing (RPI) kết hợp trong thuật toán SHADE, chứng minh rằng cơ chế thích ứng tham số lịch sử ($F$ và $CR$) hoàn toàn hoạt động tối ưu trên không gian hoán vị tổ hợp của bài toán vật lý lò phản ứng.

2. Đổi mới phương pháp luận so với các nghiên cứu tiền nhiệm? So với các nghiên cứu trước đây vốn chỉ sử dụng SA đơn thuần (như LONSA) hoặc GA cổ điển (như CIGARO), luận án đã tạo ra bước đột phá kép: phát triển mã LPO-V chuyên biệt giải FDM lưới tam giác nhanh hơn CITATION, đồng thời phát triển ESA (tích hợp toán tử lai ghép vào tôi kim) và discrete SHADE (tự thích ứng tham số), vượt trội hoàn toàn về tốc độ hội tụ và chất lượng nghiệm.

3. Phát hiện gây bất ngờ nhất về mặt số liệu thực nghiệm? Phát hiện ấn tượng nhất là cấu hình tối ưu đồng thời thỏa mãn hai mục tiêu đối nghịch nhau: tăng $k_{eff}$ thêm tới 1580 pcm (tương đương gia tăng kinh tế lớn) trong khi vẫn kéo giảm hệ số đỉnh công suất $PPF$ xuống $2.4%$ (tăng cường biên độ an toàn), phá vỡ định kiến cho rằng tăng độ dài chu kỳ luôn đi kèm với việc làm xấu phân bố công suất vùng hoạt.

4. Quy trình tái lập (Replication Protocol) có được cung cấp đầy đủ không? Luận án cung cấp chi tiết toàn bộ thông số hình học, thành phần đồng vị 6 trạng thái S1–S5, thư viện tiết diện 4 nhóm năng lượng, cấu trúc lưới 24 tam giác/bó, sơ đồ thuật toán FDM-SOR, mã giả (pseudocode) của ESA và discrete SHADE, đảm bảo khả năng tái lập kết quả tính toán $100%$.

5. Lộ trình nghiên cứu 10 năm tới được vạch ra như thế nào? Lộ trình hướng tới việc hoàn thiện nền tảng tính toán 3D đa vật lý nơ-tron - nhiệt thủy lực thời gian thực, mở rộng sang tối ưu hóa chu kỳ đa bước (multi-cycle optimization) cho các dòng lò VVER-1200 thương mại và các thiết kế lò thế hệ IV làm nguội bằng nước siêu tới hạn (VVER-SCWR).

Kết luận

Luận án tiến sĩ của tác giả 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 5 đóng góp cụ thể:

  1. Xây dựng và kiểm chứng thành công mã tính vật lý vùng hoạt LPO-V trên hình học lưới tam giác lò VVER-1000 với độ chính xác tiệm cận MCNP4c và hiệu năng vượt trội CITATION.
  2. Phát triển thành công thuật toán Evolutionary Simulated Annealing (ESA) kết hợp toán tử lai ghép tiến hóa vào cấu trúc tôi kim lân cận.
  3. Phát triển thành công thuật toán discrete SHADE ứng dụng kỹ thuật đánh chỉ số vị trí tương đối RPI cho bài toán tối ưu hóa tổ hợp rời rạc.
  4. Tìm ra cấu hình nạp tải tối ưu cho vùng hoạt VVER-1000 MOX giúp tăng $k_{eff}$ thêm 1580 pcm và giảm PPF 2.4%.
  5. Xác lập khung đánh giá thống kê phi tham số chuẩn tắc bằng Mann-Whitney U test cho các thuật toán tối ưu hóa vùng hoạt hạt nhân.

Công trình đánh dấu bước tiến vững chắc trong việc làm chủ công nghệ mô phỏng và tối ưu hóa lò phản ứng hạt nhân công suất cao, tạo tiền đề khoa học vững chắc cho các chương trình nghiên cứu ứng dụng năng lượng nguyên tử hòa bình trong tương lai.