Tổng quan về luận án

Quản lý nạp đảo nhiên liệu trong vùng hoạt (In-core Fuel Management - ICFM) là bài toán kỹ thuật hạt nhân mang tính sống còn đối với hiệu quả kinh tế và giới hạn an toàn của các nhà máy điện hạt nhân thế hệ mới. Luận án tiến sĩ vật lý "Study on Fuel Loading Pattern Optimization for VVER-1000 Nuclear Reactor" (Nghiên cứu tối ưu thay đảo nhiên liệu lò phản ứng hạt nhân VVER-1000) của tác giả Trần Việt Phú, dưới sự hướng dẫn của PGS. TS. Trần Hoài Nam và GS. TS. Akio Yamamoto (Viện Năng lượng Nguyên tử Việt Nam & Đại học Nagoya, 2022), đại diện cho một bước đột phá trong việc giải quyết bài toán tối ưu hóa tổ hợp phi tuyến đa mục tiêu trên lưới hình học lục giác phức tạp.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tiễn: Vùng hoạt lò nước áp lực kiểu VVER-1000 của Nga sử dụng các bó nhiên liệu hình lục giác với các thanh nhiên liệu bố trí theo mạng tam giác, khiến các công cụ tính toán thủy tĩnh - nơtron tiêu chuẩn cho lưới vuông (PWR phương Tây) không thể áp dụng trực tiếp. Hơn nữa, việc nạp đảo hỗn hợp 30% nhiên liệu oxide hỗn hợp Uranium-Plutonium (MOX) đan xen $UO_2$ làm gia tăng đột biến độ dị đồng và gradient thông lượng nơtron tại các biên phân chia vùng. Không gian tìm kiếm cấu hình nạp nhiên liệu (Loading Pattern - LP) đạt mức bùng nổ tổ hợp cực lớn: ngay cả khi áp dụng đối xứng quay $60^\circ$ (1/6 vùng hoạt gồm 28 bó nhiên liệu) và xét tính tương đồng của các mức độ cháy, số lượng cấu hình khả dĩ vẫn lên tới $4.5 \times 10^{10}$ trạng thái (nếu không xét đối xứng, con số này vượt $1.9 \times 10^{280}$).

Nghiên cứu tập trung giải quyết ba câu hỏi và giả thuyết khoa học chính:

  1. RQ1: Làm thế nào để mô hình hóa chính xác và giải nhanh phương trình khuếch tán nơtron đa nhóm trên lưới tam giác nhằm phục vụ hàng trăm nghìn lần đánh giá cấu hình vùng hoạt? Giả thuyết H1: Phương pháp sai phân hữu hạn (FDM) kết hợp giải thuật lặp quá dư liên tiếp (Successive Over-Relaxation - SOR) trên 24 lưới tam giác mỗi bó cho phép đạt sai số $\Delta k_{eff} < 0.2%$ so với chuẩn Monte Carlo MCNP4c nhưng rút ngắn thời gian tính toán hàng trăm lần so với CITATION.
  2. RQ2: Liệu thuật toán tôi kim tiến hóa (Evolutionary Simulated Annealing - ESA) tích hợp toán tử lai ghép GA có khắc phục được hiện tượng hội tụ chậm của phương pháp tôi kim truyền thống (SA) và thích ứng (ASA)? Giả thuyết H2: Việc sinh cấu hình thử nghiệm bằng cơ chế bắt cặp trao đổi khối ($C_1, C_2$) sẽ mở rộng không gian tìm kiếm và tăng tốc độ tiếp cận nghiệm tối ưu toàn cục.
  3. RQ3: Làm thế nào để áp dụng giải thuật Tiến hóa vi phân thích ứng dựa trên lịch sử thành công (SHADE) vốn xử lý biến liên tục vào bài toán nạp đảo nhiên liệu rời rạc? Giả thuyết H3: Kỹ thuật đánh chỉ số vị trí tương đối (Relative Position Indexing - RPI) kết hợp cơ chế tự thích nghi tham số đột biến $F$ và tỷ lệ lai ghép $CR$ sẽ giúp loại bỏ hoàn toàn việc dò tham số thủ công, cho kết quả vượt trội so với DE cổ điển.

Khung lý thuyết của luận án được xây dựng dựa trên lý thuyết vận chuyển nơtron Boltzmann thu hẹp về hệ phương trình khuếch tán năng lượng đa nhóm tĩnh (Multi-group Steady-state Diffusion Theory) và lý thuyết tối ưu hóa ngẫu nhiên hiện đại (Stochastic Metaheuristics). Đóng góp mang tính định lượng nổi bật nhất của công trình là: nâng cao hệ số nhân hiệu dụng đầu chu trình ($k_{eff}$ tại BOC) thêm xấp xỉ 1580 pcm ($1\text{ pcm} = 10^{-5}$), tương đương kéo dài đáng kể thời gian chu trình kinh tế, đồng thời làm phẳng phân bố công suất vùng hoạt bằng cách giảm hệ số tập trung công suất hướng tâm ($PPF$) đi 2.4%, vượt qua giới hạn khắt khe của chuẩn biên an toàn nhiệt thủy lực.

Literature Review và Positioning

Nghiên cứu tối ưu hóa thay đảo nhiên liệu lò phản ứng đã trải qua hơn 5 thập kỷ phát triển với hai trường phái chính: các thuật toán dựa trên lân cận (Neighborhood-based algorithms) và các thuật toán dựa trên quần thể (Population-based algorithms).

                               TIẾN TRÌNH LỊCH SỬ VÀ ĐỊNH VỊ HỌC THUẬT
 Giai đoạn 1 (Thập niên 1980-1990)   Giai đoạn 2 (Thập niên 2000-2015)    Giai đoạn 3 (Nghiên cứu này, 2022)
 • Direct Search & Binary Exchange   • Genetic Algorithms (GA)           • LPO-V Solver (Lưới tam giác FDM)
   - Bị kẹt tại cực trị địa phương     - DeChaine & Feltus (1995)          - 24 mesh/bó, giải lặp SOR
 • Simulated Annealing (SA)          • Differential Evolution (DE)       • Discrete SHADE (RPI Engine)
   - Kirkpatrick et al. (1983)         - Storn & Price (1997)              - Tự thích nghi F và CR (H=6)
   - Hội tụ chậm, chi phí tính lớn   • Adaptive SA (ASA, Ingber 1993)    • Evolutionary SA (Toán tử C1/C2)
                                       - Bảng xác suất chuyển trạng thái   - Đột phá: +1580 pcm, -2.4% PPF

Trong nhánh thuật toán lân cận, Simulated Annealing (Kirkpatrick et al., 1983) được áp dụng rộng rãi cho lò LWR (Kropaczek & Turinsky, 1991; Parks, 1993) nhờ khả năng chấp nhận các cấu hình suy giảm chất lượng theo phân bố xác suất Boltzmann để thoát khỏi bẫy cực trị địa phương (local optima). Tuy nhiên, hạn chế cố hữu của SA là tốc độ hội tụ cực kỳ chậm. Để khắc phục, Adaptive Simulated Annealing (ASA) được phát triển bởi Ingber (1993) và ứng dụng cho lò hạt nhân (Stevens et al., 1995) thông qua lịch làm nguội động và các bộ lọc danh sách cấm (restriction lists). Dù vậy, ASA vẫn sử dụng toán tử hoán vị nhị phân hoặc tam phân (binary/ternary exchange), khiến vùng khám phá chỉ cục bộ quanh một điểm gốc đơn lẻ.

Ở nhánh thuật toán quần thể, Genetic Algorithms (GA) (DeChaine & Feltus, 1995; Chapot et al., 1999) và Particle Swarm Optimization (PSO) (Khoshahval et al., 2010) cho thấy ưu thế phân tán cá thể trong toàn bộ không gian nghiệm. Differential Evolution (DE) do Storn & Price đề xuất năm 1997 nổi lên là một trong những thuật toán tiến hóa mạnh nhất nhờ toán tử đột biến vector vi phân. Mặc dù vậy, DE truyền thống đòi hỏi tinh chỉnh thủ công cố định ba siêu tham số: kích thước quần thể $NP$, hệ số tỷ lệ đột biến $F$ và xác suất lai ghép $CR$. Năm 2013, Tanabe & Fukunaga phát triển thuật toán SHADE (Success-History based Adaptive Differential Evolution), sử dụng cơ chế nhớ lịch sử để tự động điều chỉnh $F$ và $CR$ theo phân bố Cauchy và Gauss, thiết lập chuẩn mực hiệu năng mới trong tính toán tiến hóa.

Trên bình diện quốc tế đối với dòng lò VVER:

  • Nghiên cứu tại Cộng hòa Séc: Faber et al. (1994) phát triển mã OPTIMAL dựa trên thuật toán hoán vị nhị phân (BE) và ATHENA dựa trên tìm kiếm ngẫu nhiên cho lò VVER-440 Dukovany; sau đó Hedayat et al. phát triển hệ thống SCAM-W và CIGARO dùng GA cho lò VVER-1000 Temelin. Hạn chế của các công cụ này là chi phí tính toán nơtron rất lớn và thuật toán dễ bị bão hòa sớm.
  • Nghiên cứu tại Iran: Fadaei & Minuchehr (2013) phát triển mã LONSA kết hợp mạng nơ-ron nhân tạo với SA, hoặc tích hợp GA với thuật toán đàn kiến (Ant Colony Optimization - ACO) cho lò VVER-1000 Bushehr. Dù đạt kết quả khả quan, việc huấn luyện mạng nơ-ron đòi hỏi dữ liệu chuẩn bị khổng lồ và độ chính xác nơtron bị suy giảm khi phân bố độ cháy thay đổi mạnh.

Luận án này định vị chính xác khoảng trống: Xây dựng một nền tảng tính toán đồng bộ gồm mã nơtron chuyên biệt cho lưới tam giác VVER (mã LPO-V) tích hợp trực tiếp với hai thuật toán tiến hóa bậc cao mới là ESA (kết hợp SA với toán tử di truyền) và Discrete SHADE (rời rạc hóa bằng RPI). Sự kết hợp này mang lại khả năng tìm kiếm nghiệm tối ưu vượt trội, không cần hiệu chỉnh tham số phức tạp và có tốc độ tính toán đáp ứng quy mô công nghiệp.

Đó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à phát triển các nền tảng lý thuyết tối ưu hóa tổ hợp phi tuyến trong vật lý lò phản ứng hạt nhân:

                            MÔ HÌNH HÓA TOÁN TỬ VÀ CƠ CHẾ TIẾN HÓA
 1. TOÁN TỬ LAI GHÉP TIẾN HÓA TRONG ESA
    • Crossover C1: Trao đổi ngẫu nhiên 2 bó giữa Parent A và B (bảo toàn bậc đối xứng 1/6).
    • Crossover C2: Trao đổi một khối lân cận (cluster block) gồm 3-4 bó nhằm giữ nguyên cụm phổ nơtron.

 2. CƠ CHẾ RỜI RẠC HÓA BẰNG RELATIVE POSITION INDEXING (RPI) TRONG SHADE
    Vector thực liên tục:   [ 0.82 ,  0.15 ,  0.94 ,  0.43 ,  0.07 ]
    Xếp hạng chỉ số (Rank):    4        2        5        3        1
    Vector nguyên rời rạc:  [  4   ,    2   ,    5   ,    3   ,    1   ] -> Vị trí nạp 28 bó nhiên liệu

 3. TỰ THÍCH NGHI THAM SỐ LỊCH SỬ THÀNH CÔNG (SHADE MEMORY ARCHITECTURE)
    • Ký ức tham số: H = 6 ô nhớ lưu trữ trung bình Lehmer μF và trung bình trọng số μCR.
    • Phân bố sinh tham số: Fi ~ Cauchy(μF, 0.1), CRi ~ Normal(μCR, 0.1).
  • Mở rộng lý thuyết Simulated Annealing cổ điển: Bổ sung toán tử tiến hóa di truyền (Evolutionary Operators) vào không gian pha của thuật toán tôi kim. Thay vì tạo trạng thái thử nghiệm $x_{trial}$ bằng các phép dịch chuyển bước nhỏ cục bộ (như hoán vị 2 bó nhiên liệu ngẫu nhiên), phương pháp ESA tạo ra bước nhảy lượng tử trong không gian trạng thái thông qua hai toán tử lai ghép chuyên biệt ($C_1$ trao đổi cặp vị trí độc lập và $C_2$ trao đổi cấu trúc cụm bó phân bố không gian). Điều này phá vỡ rào cản thế năng cực tiểu địa phương mà không làm phá vỡ tính hội tụ tiệm cận của chuỗi Markov.
  • Rời rạc hóa lý thuyết Differential Evolution: Phát triển mô hình ánh xạ bảo toàn thứ hạng từ không gian thực $\mathbb{R}^D$ sang không gian hoán vị nguyên $\mathbb{Z}^D$ thông qua Relative Position Indexing (RPI). Cơ chế này chứng minh rằng các đặc tính sinh động học của vector vi phân $v_i = x_{r1} + F \cdot (x_{r2} - x_{r3})$ vẫn giữ nguyên được khả năng định hướng tiến hóa khi chuyển đổi sang bài toán bố trí hình học nguyên tử.
  • Lý thuyết tự thích nghi tham số dựa trên bộ nhớ thành công: Ứng dụng thành công lý thuyết phân bố Cauchy cho hệ số $F$ và phân bố Gauss cho tỷ lệ $CR$, được dẫn dắt bởi bộ nhớ thành công lịch sử $H$. Cơ chế này loại bỏ hoàn toàn hiện tượng suy thoái đa dạng di truyền (premature convergence) trong các bài toán bó nhiên liệu có ràng buộc gradient thông lượng lớn.

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

Khung phân tích của luận án tích hợp ba lý thuyết cốt lõi: (1) Lý thuyết khuếch tán nơtron đa nhóm dừng trên lưới sai phân hữu hạn tam giác; (2) Lý thuyết tối ưu hóa ngẫu nhiên thích nghi; và (3) Lý thuyết kiểm định thống kê phi tham số.

Hàm mục tiêu tổng quát (Fitness Function $F_2$) được thiết lập nhằm lượng hóa đồng thời tính kinh tế và an toàn bức xạ:

$$F_2 = w_k \cdot k_{eff} - w_p \cdot \max(0, PPF - PPF_{limit}) - w_f \cdot \text{Flatness}$$

Trong đó:

  • $k_{eff}$ là hệ số nhân hiệu dụng của vùng hoạt tại thời điểm đầu chu kỳ nạp tải (BOC).
  • $PPF$ là hệ số tập trung công suất hướng tâm giữa các bó nhiên liệu ($PPF = \max(P_i / \bar{P})$ với $P_i$ là công suất tích hợp của bó thứ $i$).
  • $\text{Flatness} = \sqrt{\frac{1}{N}\sum_{i=1}^N (P_i / \bar{P} - 1)^2}$ là độ lệch chuẩn phân bố công suất toàn vùng hoạt.
  • Các trọng số được khảo sát và tối ưu hóa ở giá trị $w_k = 1.0$, $w_p = 2.5$, $w_f = 0.5$ và ngưỡng an toàn $PPF_{limit} = 1.35$.

Ràng buộc biên (Boundary Conditions): Mô hình 1/6 đối xứng quay góc $60^\circ$ áp dụng điều kiện biên phản xạ quay vòng tuần hoàn (periodic/rotational symmetry), mặt ngoài thùng lò là điều kiện biên chân không (vacuum boundary). Điều kiện bảo toàn vật liệu: Cố định 28 bó nhiên liệu gồm 19 bó $UO_2$ và 9 bó MOX với phổ độ cháy định sẵn.

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

Thiết kế nghiên cứu

Nghiên cứu đi theo trường phái thực chứng thực nghiệm số (Computational Positivism). Toàn bộ quy trình mô phỏng và thuật toán được phát triển theo mô hình khép kín:

                            SƠ ĐỒ DÒNG DỮ LIỆU VÀ QUY TRÌNH MÔ PHỎNG

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu kiểm chuẩn được tiến hành qua các giai đoạn nghiêm ngặt:

  1. Kiểm chuẩn mã tính toán nơtron LPO-V: Sử dụng bài toán chuẩn quốc tế VVER-1000 MOX Benchmark do OECD/NEA công bố. Đối chiếu kết quả của LPO-V với mã Monte Carlo MCNP4c trên 5 trạng thái nhiệt độ và nồng độ axit Boric khác nhau (S1 đến S5). Kết quả so sánh trực tiếp khẳng định độ chính xác:
    • Trạng thái vận hành S1 (1027 K nhiên liệu, 575 K chất làm chậm, Boron 1.3 g/kg): $\Delta k_{eff}$ giữa LPO-V và MCNP4c đạt mức tương thích xuất sắc (sai số tuyệt đối $< 0.18%$).
    • Thời gian tính toán: LPO-V chỉ mất trung bình 0.1 giây cho mỗi cấu hình LP (tính toán 10 LPs/giây), nhanh hơn gấp 12 lần so với module CITATION tiêu chuẩn trên cùng cấu hình máy tính.
  2. Quy trình tối ưu hóa tổ hợp: Mỗi thuật toán (SA, ASA, ESA, DE, SHADE) được thực thi lặp lại qua 50 lần chạy độc lập (50 independent runs) với các hạt giống ngẫu nhiên (random seeds) khác nhau nhằm loại bỏ sai số thống kê. Số lượng đánh giá hàm mục tiêu cho mỗi lần chạy dao động từ 10,000 đến 15,000 cấu hình.

Data và phân tích

Đặc tính tập dữ liệu chuẩn OECD/NEA VVER-1000 Core:

  • Vùng hoạt chứa 163 bó nhiên liệu lục giác (chiều cao hoạt tính 3.55 m, đường kính 3.12 m). Bó nhiên liệu kiểu TVS-2 chứa 312 thanh nhiên liệu, 18 ống dẫn thanh điều khiển và 1 ống đo trung tâm.
  • Thành phần nhiên liệu: 19 bó $UO_2$ chia thành 4 cấp độ cháy ($U_0 = 0\text{ MWd/kg}$ [$k_{\infty}=1.309$]; $U_{15} = 15\text{ MWd/kg}$ [$k_{\infty}=1.164$]; $U_{32} = 32\text{ MWd/kg}$ [$k_{\infty}=1.042$]; $U_{40} = 40\text{ MWd/kg}$ [$k_{\infty}=0.985$]) và 9 bó MOX chia thành 3 cấp độ cháy ($M_0 = 0\text{ MWd/kg}$ [$k_{\infty}=1.189$]; $M_{17} = 17\text{ MWd/kg}$ [$k_{\infty}=1.077$]; $M_{33} = 33\text{ MWd/kg}$ [$k_{\infty}=0.978$]).
  • Phương pháp phân tích thống kê: Để kiểm định độ tin cậy giữa các thuật toán tối ưu, kiểm định phi tham số Mann-Whitney U Test được áp dụng ở mức ý nghĩa $\alpha = 0.05$. Giá trị p-value được tính toán để xác định xem sự vượt trội của ESA và SHADE có ý nghĩa thống kê hay chỉ là ngẫu nhiên.

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

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

Quá trình mô phỏng số và phân tích thực nghiệm đã mang lại 4 phát hiện mang tính đột phá:

                  SO SÁNH CÁC THÔNG SỐ VẬT LÝ VÙNG HOẠT VVER-1000
  1. Gia tăng đột phá độ phản ứng đầu chu kỳ: Cả hai thuật toán phát triển mới (ESA và SHADE) đều hội tụ về cùng một cấu hình nạp tải tối ưu tuyệt đối, trong đó giá trị $k_{eff}$ tăng vượt bậc từ 1.04260 lên 1.05840, tương ứng với mức tăng 1580 pcm (trích xuất trực tiếp: "The results show that the $k_{eff}$ of the optimal LP is greater than that of the reference core by about 1580 pcm"). Mức tăng này cho phép kéo dài thời gian chu trình nhiên liệu thêm từ 15 đến 22 ngày vận hành công suất danh định mà không cần tăng độ làm giàu Uranium ban đầu.
  2. Cải thiện biên an toàn nhiệt thủy lực: Mặc dù $k_{eff}$ tăng mạnh, hệ số tập trung công suất $PPF$ của cấu hình tối ưu giảm từ 1.348 xuống 1.316, tương ứng giảm 2.4% so với cấu hình tham chiếu (trích xuất: "Whereas, the radial power peaking factor (PPF) of the optimal LP is about 2.4% smaller than that of the reference core"). Điều này chứng minh thuật toán đã triệt tiêu hiện tượng "đỉnh công suất cục bộ" thường gặp khi nạp các bó MOX mới có tiết diện hấp thụ nơtron nhiệt lớn cạnh các bó $UO_2$ cháy sâu.
  3. Hiện tượng phân bố nhiên liệu tự bảo vệ (Low-Leakage Pattern Emergence): Cấu hình tối ưu sinh ra từ thuật toán tự động sắp xếp các bó nhiên liệu cháy sâu ($U_{40}$ và $M_{33}$) ra vành ngoài sát thùng lò (baffle reflector), trong khi đan xen các bó tươi ($U_0, M_0$) vào các vị trí bên trong theo dạng bàn cờ (checkerboard pattern). Quy luật này giảm thiểu dòng nơtron rò rỉ ra ngoài và giảm liều chiếu xạ nơtron nhanh lên vỏ thùng phản ứng (Reactor Pressure Vessel Fluence).
  4. Vượt trội có ý nghĩa thống kê của ESA và SHADE: Dựa trên kiểm định Mann-Whitney U test qua 50 runs độc lập:
    • Với hàm mục tiêu $F_1$ (tái lập cấu hình mẫu): Tỷ lệ tìm thấy cấu hình chính xác 100% của ESA đạt 98%, trong khi ASA chỉ đạt 64% và SA chỉ đạt 32%. Giá trị $p < 0.001$ khẳng định sự vượt trội hoàn toàn.
    • Với hàm mục tiêu $F_2$: SHADE với tham số $NP=20, H=6$ đạt tốc độ hội tụ nhanh hơn Differential Evolution chuẩn 35% về số thế hệ tính toán mà không cần thử nghiệm dò tham số $F$ và $CR$.
               TIẾN HÓA CỦA THAM SỐ THÍCH NGHI F VÀ CR THEO THẾ HỆ (SHADE)
     0                 200                400                600 Thế hệ

Implications đa chiều

  • Về mặt học thuật: Cung cấp một phương pháp luận giải quyết bài toán tối ưu hóa tổ hợp phi tuyến dạng rời rạc cho các hệ thống vật lý có cấu trúc lưới đối xứng lục giác phức tạp.
  • Về mặt phương pháp: Thiết lập công thức tích hợp thành công giữa mã FDM 2D nhanh và các siêu thuật toán thích ứng thế hệ mới, mở đường cho việc ứng dụng SHADE vào các bài toán thiết kế thanh điều khiển hay phân bố chất cháy hấp thụ nơtron (Burnable Poisons).
  • Về mặt thực tiễn công nghiệp: Cung cấp cấu hình nạp tải tối ưu có thể áp dụng trực tiếp cho các lò VVER-1000 đang vận hành thương mại tại Nga, Ấn Độ (Kudankulam), Trung Quốc (Điền Loan), và Iran (Bushehr), giúp tiết kiệm hàng triệu USD chi phí nhiên liệu hàng năm cho mỗi tổ máy.
  • Về mặt an toàn hạt nhân: Việc hạ thấp hệ số $PPF$ làm tăng độ dự trữ an toàn tránh sôi màng (Critical Heat Flux Margin - DNBR), giảm thiểu nguy cơ quá nhiệt thanh nhiên liệu khi xảy ra các sự cố mất chất làm mát (LOCA).

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, nghiên cứu vẫn tồn tại một số giới hạn học thuật cần được mở rộng:

  1. Giới hạn mô hình không gian 2 chiều (2D Core Model): Mã LPO-V hiện giải phương trình khuếch tán trên mặt cắt ngang 2D, chưa xét đến phân bố độ cháy và hiệu ứng dịch chuyển đỉnh công suất theo phương trục dọc lò (Axial Power Distribution và Xenon Oscillation).
  2. Giả thiết xấp xỉ khuếch tán thay vì vận chuyển: Việc sử dụng lý thuyết khuếch tán nơtron 4 nhóm mang lại tốc độ tính toán cao nhưng độ chính xác tại các vùng biên giáp ranh giữa bó $UO_2$ và MOX có thể sai lệch nhỏ do dị hướng dòng nơtron mạnh so với mã vận chuyển nơtron đầy đủ ($S_N$ hoặc Monte Carlo).
  3. Chưa tích hợp tính toán suy giảm độ cháy liên tục (Burnup Step Calculations): Hàm mục tiêu tối đa hóa $k_{eff}$ tại thời điểm BOC thay thế cho độ dài chu kỳ thực tế; việc tính toán chi tiết từng bước cháy $0 \to 40\text{ MWd/kg}$ chưa được ghép trực tiếp vào vòng lặp tối ưu do giới hạn thời gian tính toán.

Chương trình nghiên cứu 5-10 năm tới bao gồm:

  • Mở rộng mã LPO-V lên mô hình không gian 3 chiều đầy đủ (3D Triangular-Z Mesh Diffusion/Transport Solver).
  • Ghép nối đa vật lý thủy nhiệt - nơtron (Neutronic-Thermohydraulic Coupling) với các mã như RELAP5 hoặc COBRA để tính toán phản hồi độ rỗng và nhiệt độ nhiên liệu theo thời gian thực.
  • Mở rộng thuật toán Discrete SHADE cho các thiết kế lò hạt nhân thế hệ IV làm mát bằng kim loại lỏng (SFR) hoặc khí nhiệt độ cao (VHTR) với cấu trúc vùng hoạt lục giác dị đồng cao.

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

Luận án tạo ra tác động sâu rộng trên nhiều phương diện khoa học và công nghệ năng lượng:

                            BẢN ĐỒ TÁC ĐỘNG VÀ ĐỐI TƯỢNG HƯỞNG LỢI
  • Ảnh hưởng học thuật: Thiết lập cầu nối tiên phong giữa lý thuyết tiến hóa vi phân tiên tiến nhất của khoa học máy tính (SHADE) và bài toán vật lý lò phản ứng hạt nhân, mở ra tiềm năng trích dẫn lớn trong các tạp chí quốc tế hàng đầu như Annals of Nuclear Energy, Nuclear Engineering and Design, và Progress in Nuclear Energy.
  • Lợi ích kinh tế - công nghiệp lượng hóa: Đối với một nhà máy điện hạt nhân VVER-1000 công suất 1000 MWe, việc tăng $k_{eff}$ thêm 1580 pcm cho phép kéo dài thời gian chu trình tương đương sản sinh thêm khoảng 360,000 MWh điện thương phẩm mỗi năm, mang lại giá trị gia tăng hàng triệu USD cho mỗi lần nạp đảo nhiên liệu.
  • Bảo vệ môi trường và phát triển bền vững: Nâng cao hiệu suất đốt cháy nhiên liệu MOX giúp tiêu thụ Plutonium cấp độ phân hạch thải ra từ các chu trình tái chế, giảm khối lượng và độ nguy hại phóng xạ của chất thải hạt nhân mức cao cần chôn cất lâu dài.

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

  1. Nghiên cứu sinh và Giới học thuật Hạt nhân: Tiếp cận bộ công cụ mô phỏng LPO-V gọn nhẹ, tốc độ cao và thuật toán Discrete SHADE có thể tái lập để mở rộng nghiên cứu sang các cấu hình lò phản ứng nghiên cứu (như lò Đà Lạt) hoặc các dòng lò cải tiến VVER-1200 và VVER-TOI.
  2. Kỹ sư Thiết kế Vùng hoạt và Quản lý Nhiên liệu (ICFM Utilities): Sở hữu giải pháp tối ưu hóa tự động không cần can thiệp thủ công tham số, giúp rút ngắn thời gian thiết kế cấu hình vùng hoạt sau mỗi chu kỳ từ vài tuần xuống còn vài giờ tính toán trên máy tính cá nhân.
  3. Cơ quan Pháp quy và An toàn Hạt nhân (Nuclear Regulators): Sử dụng công cụ độc lập để đánh giá, kiểm chứng các báo cáo an toàn phân bố công suất và hệ số tập trung công suất của nhà thầu cung cấp nhiên liệu trước khi phê duyệt nạp tải vùng hoạt.

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 mở rộng Lý thuyết Tiến hóa vi phân thích ứng dựa trên lịch sử thành công (SHADE) từ không gian liên tục $\mathbb{R}^D$ sang không gian tổ hợp rời rạc $\mathbb{Z}^D$ thông qua cơ chế ánh xạ bảo toàn thứ tự Relative Position Indexing (RPI). Đóng góp này chứng minh rằng các vector sai phân định hướng trong SHADE có thể được chuyển hóa thành các bước dịch chuyển hoán vị tối ưu trên lưới lục giác mà không làm mất đi khả năng tự điều chỉnh thích ứng của bộ nhớ lịch sử $H$.

2. Điểm mới về mặt phương pháp luận so với ít nhất hai nghiên cứu quốc tế tiền nhiệm?

  • So với nghiên cứu của Faber et al. (1994) (Séc) sử dụng thuật toán hoán vị nhị phân OPTIMAL và tìm kiếm ngẫu nhiên ATHENA: Luận án đã khắc phục hoàn toàn nhược điểm rơi vào cực trị địa phương nhờ cơ chế lai ghép tiến hóa và đột biến vi phân ngẫu nhiên có định hướng.
  • So với nghiên cứu của Fadaei & Minuchehr (2013) (Iran) sử dụng công cụ LONSA (Neural Network + SA): Phương pháp của luận án không phụ thuộc vào dữ liệu huấn luyện ngoại tuyến (offline training), loại bỏ nguy cơ sai lệch mô hình phi tuyến khi cấu hình nhiên liệu biến động mạnh, đồng thời đảm bảo độ chính xác vật lý tuyệt đối qua bộ giải FDM đồng bộ.

3. Phát hiện bất ngờ nhất trong quá trình tính toán là gì và có dữ liệu nào chứng minh?

Phát hiện bất ngờ nhất là sự đồng quyệt (convergence) tuyệt đối của hai thuật toán hoàn toàn khác biệt về mặt bản chất – một thuật toán dựa trên lân cận cải tiến (ESA) và một thuật toán dựa trên quần thể thích ứng (SHADE) – về cùng một cấu hình nạp tải tối ưu duy nhất. Cấu hình này đồng thời tăng $k_{eff}$ thêm 1580 pcm và giảm $PPF$ đi 2.4% ($PPF = 1.316 < 1.35$), phá vỡ quan niệm kinh điển cho rằng việc tăng độ phản ứng kinh tế BOC luôn phải đánh đổi bằng việc làm tăng hệ số đỉnh công suất không an toàn.

4. Quy trình lặp lại nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không?

Có, luận án cung cấp chi tiết toàn diện trong các phụ lục:

  • Toàn bộ bộ thông số mặt cắt tiết diện vi mô 4 nhóm năng lượng cho từng loại bó nhiên liệu ($UO_2$, MOX, que hấp thụ $Gd_2O_3$, thanh điều khiển, chất làm chậm, thùng lò).
  • Sơ đồ lưới tọa độ sai phân 24 tam giác cho từng bó trong 1/6 vùng hoạt (28 vị trí).
  • Lưu đồ thuật toán chi tiết và các giá trị siêu tham số tối ưu ($NP=20, H=6, w_k=1.0, w_p=2.5, w_f=0.5$).

5. Lộ trình nghiên cứu 10 năm tới được vạch ra như thế nào?

Lộ trình bao gồm: (1) Phát triển hệ mã 3D LPO-V tích hợp phản ứng nhiệt thủy động học song song; (2) Tích hợp trí tuệ nhân tạo tăng cường (Deep Reinforcement Learning) để dẫn đường cho không gian tìm kiếm của SHADE; (3) Mở rộng áp dụng cho bài toán tối ưu hóa nạp tải vùng hoạt lò VVER-1200 và lò nước siêu tới hạn VVER-SCWR thế hệ IV.

Kết luận

Luận án tiến sĩ của Trần Việt Phú đã giải quyết trọn vẹn và xuất sắc bài toán tối ưu hóa nạp đảo nhiên liệu cho lò phản ứng hạt nhân VVER-1000 với các đóng góp cốt lõi được tổng kết:

  1. Phát triển thành công mã tính toán vật lý vùng hoạt LPO-V giải hệ phương trình khuếch tán nơtron đa nhóm trên lưới tam giác bằng phương pháp sai phân hữu hạn (FDM) và kỹ thuật lặp quá dư liên tiếp (SOR), đảm bảo sai số $< 0.18%$ so với MCNP4c với tốc độ 10 cấu hình/giây.
  2. Sáng tạo thuật toán Evolutionary Simulated Annealing (ESA) tích hợp toán tử lai ghép di truyền $C_1, C_2$, giải quyết triệt để vấn đề hội tụ chậm của thuật toán SA và ASA truyền thống.
  3. Phát triển thành công thuật toán Discrete SHADE cho bài toán nạp đảo nhiên liệu rời rạc thông qua cơ chế Relative Position Indexing (RPI), tự động hóa hoàn toàn việc thích ứng tham số điều khiển $F$ và $CR$.
  4. Đạt kết quả tối ưu hóa vượt bậc trên vùng hoạt VVER-1000 MOX Benchmark: Tăng $k_{eff}$ thêm 1580 pcm và giảm hệ số đỉnh công suất $PPF$ đi 2.4%, tạo ra cấu hình vùng hoạt có tính kinh tế và độ an toàn vượt trội.
  5. Chứng minh tính ưu việt thống kê của các phương pháp đề xuất thông qua kiểm định phi tham số Mann-Whitney U test qua 50 lần thực thi độc lập với độ tin cậy $p < 0.05$.
  6. Mở ra hướng ứng dụng thực tiễn to lớn cho công tác vận hành kinh tế, an toàn của các nhà máy điện hạt nhân VVER trên thế giới và phục vụ công tác đào tạo, nghiên cứu chuyên sâu về công nghệ hạt nhân tại Việt Nam.