Tổng quan về luận án

Quản lý nhiên liệu trong vùng hoạt (In-core Fuel Management - ICFM) là một trong những bài toán tối ưu hóa đa mục tiêu phức tạp và có tính quyết định nhất trong kỹ thuật lò phản ứng hạt nhân. Trong chu trình vận hành của lò phản ứng nước áp lực thế hệ mới, chi phí nhiên liệu chiếm khoảng 15–25% tổng chi phí sản xuất điện ròng, trong khi cấu hình nạp tải nhiên liệu (Loading Pattern - LP) trực tiếp chi phối cả tính kinh tế (độ dài chu trình cháy, mức độ tiêu thụ nhiên liệu) lẫn giới hạn an toàn nhiệt thủy lực (hệ số đỉnh công suất xuyên tâm - $PPF$). 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. TS. Trần Hoài Nam và GS. Akio Yamamoto tại Viện Năng lượng Nguyên tử Việt Nam (VAEI / VINATOM), đã giải quyết bài toán cốt lõi này thông qua công trình mang tính tiên phong: "Nghiên cứu tối ưu thay đảo nhiên liệu lò phản ứng hạt nhân VVER-1000" (Mã số chuyên ngành: 9.44.01.06 - Vật lý Nguyên tử và Hạt nhân).

Trong thực tế thiết kế vùng hoạt lò phản ứng VVER-1000 chứa 30% nhiên liệu oxit hỗn hợp (Mixed Oxide - MOX), không gian tìm kiếm cấu hình khả dĩ là một tổ hợp khổng lồ: đối với vùng hoạt 163 bó nhiên liệu không xét tính đối xứng, số lượng cấu hình lên tới $1.9 \times 10^{280}$; ngay cả khi áp dụng đối xứng quay $1/6$ và gom nhóm tương đồng về độ cháy, không gian tìm kiếm vẫn đạt mức $4.5 \times 10^{10}$ cấu hình. Khoảng trống nghiên cứu (research gap) cốt lõi nằm ở chỗ: các thuật toán tối ưu hóa truyền thống như tìm kiếm trực tiếp (Direct Search), trao đổi nhị phân (Binary Exchange) hay thậm chí thuật toán mô phỏng tôi kim cổ điển (Simulated Annealing - SA) và giải thuật di truyền (Genetic Algorithm - GA) thường bị mắc kẹt tại các cực trị địa phương (local optima), tốc độ hội tụ chậm, hoặc đòi hỏi chi phí tính toán vật lý nơ-tron quá lớn nếu không có một công cụ mô phỏng chuyên dụng cho hình học mạng lưới tam giác lục giác đặc thù của công nghệ VVER.

Luận án đã thiết lập và kiểm chứng 2 câu hỏi nghiên cứu cùng hệ giả thuyết khoa học:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để xây dựng một bộ giải nơ-tron 2 chiều nhanh, chính xác, tương thích hoàn toàn với hình học lưới tam giác của bó nhiên liệu lục giác VVER nhằm phục vụ hàng trăm nghìn lần đánh giá cấu hình?
    • Giả thuyết 1 (H1): Mô hình sai phân hữu hạn (FDM) giải phương trình khuếch tán 4 nhóm năng lượng kết hợp thuật toán lặp quá giãn liên tiếp (SOR) sẽ mang lại độ chính xác tương đương Monte Carlo MCNP4c nhưng rút ngắn thời gian tính toán xuống dưới 0.1 giây/cấu hình.
  • Câu hỏi nghiên cứu 2 (RQ2): Cơ chế toán tử lai ghép và tự thích ứng thông số điều khiển nào có khả năng vượt qua rào cản hội tụ cục bộ trong không gian rời rạc phi tuyến tính của bài toán ICFM?
    • Giả thuyết 2 (H2): Sự tích hợp toán tử trao đổi chéo (Crossover) vào quy trình tôi kim (ESA) và việc áp dụng cơ chế bộ nhớ lịch sử thích ứng tham số kết hợp kỹ thuật lập chỉ số vị trí tương đối (RPI) trong tiến hóa vi phân (SHADE) sẽ vượt trội có ý nghĩa thống kê ($p < 0.05$) so với SA, ASA và DE tiêu chuẩn.

Khung lý thuyết của nghiên cứu tích hợp lý thuyết khuếch tán nơ-tron nhiều nhóm năng lượng, cơ học thống kê phi cân bằng Metropolis-Hastings và tính toán tiến hóa ngẫu nhiên thích ứng. Kết quả định lượng đột phá của luận án xác nhận 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}$ ở đầu chu trình (BOC) thêm 1580 pcm (tương đương kéo dài đáng kể thời gian chu trình vận hành) đồ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 tham chiếu tiêu chuẩn OECD/NEA Benchmark, tạo ra biên độ an toàn nhiệt vượt trội cho hệ thống thanh nhiên liệu.

Literature Review và Positioning

Lịch sử tối ưu hóa cấu hình nạp tải nhiên liệu hạt nhân đã trải qua nhiều giai đoạn phát triển từ các phương pháp tất định sơ khai đến các thuật toán heuristic và meta-heuristic hiện đại. Các dòng nghiên cứu chính trong y văn thế giới bao gồm:

  1. Dòng thuật toán dựa trên lân cận (Neighborhood-based Search): Tiêu biểu là thuật toán Mô phỏng tôi kim (Kirkpatrick et al., 1983; Parks, 1990; Kropaczek & Turinsky, 1991) và Tìm kiếm Tabu (Tabu Search - Glover, 1986). SA cho phép chấp nhận các nghiệm kém hơn theo xác suất phân bố Boltzmann $P = \exp(-\Delta E / T)$, giúp thuật toán có khả năng thoát khỏi các cực tiểu/cực đại địa phương. Tuy nhiên, nhược điểm cốt tử của SA truyền thống là tốc độ hội tụ rất chậm do cơ chế sinh nghiệm thử nghiệm chỉ dựa trên các phép tráo đổi nhị phân (Binary Exchange) hoặc tam phân (Ternary Exchange) ngẫu nhiên cục bộ. Để khắc phục, Ingber (1993) và sau đó là các nghiên cứu của Stevens et al. (1995) đã phát triển phương pháp Tôi kim Thích ứng (Adaptive Simulated Annealing - ASA) với lịch làm nguội động và danh sách cấu hình hạn chế.
  2. Dòng thuật toán dựa trên quần thể (Population-based Search): Bao gồm Giải thuật Di truyền (GA - DeChaine & Feltus, 1995; Chapot et al., 1999; Yamamoto, 2004), Tối ưu bầy đàn (PSO - Alim et al., 2008; Khoshahval et al., 2014), Tối ưu đàn kiến (ACO - Babazadeh et al., 2009) và Tiến hóa vi phân (Differential Evolution - DE - Storn & Price, 1997). DE nổi bật với khả năng khám phá không gian tìm kiếm mạnh mẽ thông qua các toán tử đột biến véc-tơ hiệu và lai ghép nhị phân. Tuy nhiên, DE vốn được thiết kế cho biến liên tục thực $\mathbb{R}^D$, do đó việc ánh xạ vào không gian tổ hợp rời rạc nguyên của bài toán nạp tải nhiên liệu luôn gặp thách thức lớn về bảo toàn tính hợp lệ của cấu hình (không trùng lặp vị trí bó nhiên liệu).
  3. Dòng thuật toán thích ứng nâng cao: Tanabe & Fukunaga (2013) đã tạo ra bước đột phá khi đề xuất thuật toán SHADE (Success-History based Adaptive Differential Evolution), sử dụng bộ nhớ lịch sử $H$ để tự động điều chỉnh hệ số khuếch đại đột biến $F$ và tỷ lệ lai ghép $CR$ dựa trên phân bố Cauchy và phân bố chuẩn của các cá thể thành công trong các thế hệ trước.
Thuật toán / Công cụ Đối tượng lò / Cấu hình Phương pháp cốt lõi Hạn chế chính trong y văn Đóng góp vượt trội của Luận án
OPTIMAL / ATHENA (Czech Republic) Lò VVER-440 (Dukovany) Binary Exchange & Random Search thuần túy Dễ rơi vào cực trị địa phương, số lượng cấu hình đánh giá hạn chế Tích hợp toán tử Crossover toàn cục (ESA) giúp vượt rào cản thế năng địa phương
SCAM-W / CIGARO (Czech Republic) Lò VVER-1000 (Temelin) Standard Genetic Algorithms (GA) kết hợp bộ giải thương mại Hiệu suất hội tụ phụ thuộc vào tham số cố định, chi phí tính toán vật lý cao Phát triển discrete SHADE thích ứng động $F, CR$ kết hợp bộ giải nhanh LPO-V
LONSA Tool (Iran - Fadaei et al.) Lò VVER-1000 (Bushehr) Mạng nơ-ron nhân tạo kết hợp Simulated Annealing Mạng nơ-ron đòi hỏi tập dữ liệu huấn luyện lớn, sai số xấp xỉ công suất cục bộ Mô phỏng vật lý trực tiếp FDM 4 nhóm chuẩn xác, không bị phụ thuộc dữ liệu offline
LPO-V / ESA & SHADE (Luận án Trần Việt Phú) Lò VVER-1000 (OECD/NEA MOX Benchmark) Tri-mesh FDM + ESA (C1/C2) + Discrete SHADE (RPI) Cần tiếp tục mở rộng sang mô hình 3D toàn chu trình cháy Đạt mức tăng $k_{eff}$ +1580 pcm, giảm $PPF$ -2.4%, kiểm chứng vững chắc qua Mann-Whitney U-test

Các tranh luận khoa học lớn trong y văn tập trung vào hai quan điểm đối lập: một bên ủng hộ các phương pháp heuristic cục bộ tăng cường (như ASA cải tiến) nhằm tiết kiệm số lần gọi mô phỏng vật lý nơ-tron; bên kia ủng hộ các thuật toán tiến hóa quần thể có tính đa dạng di truyền cao (như GA, DE) nhưng phải đối mặt với chi phí tính toán bùng nổ. Luận án của Trần Việt Phú định vị chính xác vào giao điểm này: phát triển giải pháp lai ghép ESA (kết hợp khả năng tìm kiếm sâu của SA với sức mạnh trao đổi cấu trúc quy mô lớn của toán tử GA Crossover) và Discrete SHADE (sử dụng cơ chế Relative Position Indexing để chuyển đổi không gian liên tục sang hoán vị số nguyên), đồng thời tự xây dựng bộ giải nơ-tron LPO-V chuyên biệt chạy ở tốc độ cao để giải phóng hoàn toàn nghẽn cổ chai tính toá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 phong phú hệ thố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 thông qua các đóng góp nền tảng:

  1. Mở rộng lý thuyết tôi kim trong không gian hoán vị cấu hình nạp tải: Thách thức tiên đề truyền thống của Kirkpatrick rằng nghiệm thử nghiệm mới trong SA chỉ nên tạo ra từ các nhiễu loạn nhỏ cục bộ (Local Perturbation / Binary Exchange). Luận án chứng minh bằng toán học và thực nghiệm số rằng việc đưa toán tử trao đổi khối ngẫu nhiên (Block Crossover - C2) và tráo đổi hai bó ngẫu nhiên giữa hai cấu hình cha mẹ (Two-assembly Crossover - C1) vào lõi thuật toán tôi kim (ESA) giúp phá vỡ các rào cản thế năng liên kết giữa các bó nhiên liệu có độ cháy khác nhau, tạo bước nhảy lượng tử trong không gian trạng thái mà không làm mất tính ổn định nhiệt động của hệ thống.
  2. Hình thức hóa cơ chế thích ứng SHADE trong miền tổ hợp rời rạc: DE truyền thống yêu cầu 3 tham số cố định ($N_P, F, CR$). Luận án đã tích hợp cơ chế thích ứng dựa trên lịch sử thành công của Tanabe & Fukunaga vào bài toán vật lý hạt nhân, biến đổi $F$ và $CR$ thành các biến nội sinh tự điều chỉnh theo tiến trình hội tụ thế hệ.
  3. Mô hình hóa toán học bài toán tối ưu hóa đa mục tiêu: Hàm thích nghi (Fitness Function - FF) được chuẩn hóa chặt chẽ dưới dạng: $$\text{Maximize } \mathcal{F}1 = \frac{1}{\sum{i=1}^{N_{FA}} |P_i - P_{i}^{ref}| + \epsilon}$$ $$\text{Maximize } \mathcal{F}2 = k{eff} \cdot \left(1 - w_p \cdot \max(0, PPF - PPF_{limit})\right) \cdot \left(1 + w_f \cdot \text{Flatness}\right)$$ trong đó các trọng số phạt $w_p = 2.5$ và trọng số làm phẳng $w_f = 0.05$ được khảo sát tham số tối ưu hóa nhằm tạo ra sự cân bằng hoàn hảo giữa hiệu suất kinh tế nơ-tron và an toàn nhiệt.

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

Khung phân tích của luận án là sự hợp nhất của 3 trụ cột lý thuyết chuyên sâu:

Phương pháp tiếp cận giải tích mới mẻ này xác định rõ ràng các điều kiện biên: áp dụng đối xứng quay góc $60^\circ$ (1/6 core symmetry), điều kiện biên phản xạ (Reflective Boundary Condition) tại các mặt phẳng đối xứng bên trong và điều kiện biên chân không / bề mặt tự do (Free Surface Boundary Condition) tại mặt ngoài của vỏ lò phản ứng hạt nhân.

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

Thiết kế nghiên cứu

Nghiên cứu theo đuổi trường phái nhận thức luận thực chứng (Positivism / Computational Physics Paradigm) với độ tin cậy được bảo chứng qua mô hình toán số tất định kết hợp tối ưu hóa ngẫu nhiên stochastics. Đối tượng thử nghiệm trực tiếp là bài toán chuẩn quốc tế OECD/NEA VVER-1000 MOX Benchmark Core.

Quy trình nghiên cứu và Phát triển mã LPO-V

Quy trình nghiên cứu được triển khai theo chuỗi giá trị tính toán khép kín và nghiêm ngặt:

  1. Chuẩn bị số liệu tiết diện vi mô và vĩ mô (Group Constants Generation): Sử dụng mô-đun PIJ trong hệ mã chuẩn SRAC (Japan Atomic Energy Agency) với thư viện hạt nhân liên tục ENDF/B-VII.0. Cấu trúc 107 nhóm năng lượng (62 nhóm nhanh, 45 nhóm nhiệt) được tính toán va chạm xác suất trong ô mạng mạng lục giác 400 vùng hình học, sau đó rút gọn (collapsed) thành hệ tiết diện vĩ mô 4 nhóm năng lượng (3 nhóm nhanh, 1 nhóm nhiệt) phục vụ tính toán vùng hoạt.

  2. Thiết lập bộ giải LPO-V (Loading Pattern Optimization of VVER): Tác giả tự lập trình mã nguồn giải hệ phương trình khuếch tán đa nhóm dừng: $$-\nabla \cdot (D_g \nabla \phi_g) + \Sigma_{r,g} \phi_g = S_g$$ $$S_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$$ Phương pháp sai phân hữu hạn (FDM) được thiết lập trên lưới tam giác 2D với 24 lưới tam giác trên mỗi bó nhiên liệu lục giác. Dòng nơ-tron thực giữa hai lưới tam giác lân cận $i$ và $k$ được tính toán qua hệ số khuếch tán liên kết: $$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 đồng thời bằng phương pháp Lặp quá giãn liên tiếp (Successive Over-Relaxation - SOR) với hệ số gia tốc tối ưu hóa, kết hợp kỹ thuật lặp nguồn nơ-tron (Power Iteration).

  3. Xác minh (Verification) mã LPO-V: Bộ giải LPO-V được đối chuẩn trực tiếp với mã Monte Carlo MCNP4c trên 5 trạng thái nhiệt động khác nhau của vùng hoạt benchmark (S1: Trạng thái làm việc đầy tải; S2: Nhiệt độ đồng nhất; S3: Trạng thái nguội nồng độ Boron cao; S4: Trạng thái làm việc không Boron; S5: Nhiệt độ đồng nhất không Boron).

  4. Kiểm chứng thống kê phi tham số: Để tránh hiện tượng ngẫu nhiên do khởi tạo hạt giống (seed), mọi thuật toán (SA, ASA, DE, ESA, SHADE) đều được thực thi lặp lại 50 lần chạy độc lập (50 independent runs). Sự sai khác về hiệu năng hội tụ và chất lượng nghiệm được thẩm định thông qua kiểm định Mann-Whitney U-test ở mức ý nghĩa thống kê chuẩn $\alpha = 0.05$.

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

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

Các phân tích số liệu thực nghiệm trong luận án đã mang lại 5 phát hiện mang tính đột phá cho ngành vật lý lò phản ứng:

           Hệ số nhân hiệu dụng keff (BOC)             Hệ số đỉnh công suất xuyên tâm (PPF)
  1. Hiệu năng vượt bậc của cấu hình nạp tải tối ưu: Cấu hình tối ưu được tìm thấy bởi cả hai phương pháp ESA và discrete SHADE đạt hệ số nhân nơ-tron hiệu dụng $k_{eff} = 1.12495$, tăng vượt trội 1580 pcm so với cấu hình tham chiếu của OECD/NEA Benchmark ($k_{eff}^{ref} = 1.10915$). Đồng thời, hệ số đỉnh công suất bó $PPF$ giảm từ $1.354$ xuống còn $1.321$ (giảm 2.4%), đáp ứng hoàn toàn giới hạn an toàn nhiệt $PPF < 1.35$.
  2. Sự vượt trội của toán tử Crossover trong thuật toán ESA: Trong phép thử nghiệm tái tạo cấu hình chuẩn ($\mathcal{F}_1$), phương pháp ESA (đặc biệt là biến thể ESA-C1A4 sử dụng lai ghép 2 bó kết hợp 4 phép tráo đổi ngẫu nhiên) đạt tỷ lệ tái tạo thành công 100% trong 50 lần chạy độc lập, trong khi SA nguyên bản chỉ đạt 18% và ASA chỉ đạt 68%.
  3. Tính thích ứng linh hoạt của discrete SHADE: Khảo sát tham số điều khiển cho thấy kích thước quần thể $N_P = 20$ và kích thước bộ nhớ lịch sử $H = 6$ là cấu hình tối ưu nhất cho bài toán nạp tải VVER-1000. Các tham số đột biến $F$ và lai ghép $CR$ tự động điều chỉnh mượt mà từ giai đoạn phân tán ban đầu hướng đến các giá trị hội tụ ổn định quanh $F \approx 0.55$ và $CR \approx 0.85$ ở các thế hệ sau.
  4. Bằng chứng thống kê từ Mann-Whitney U Test: Kết quả kiểm định xác nhận sự khác biệt giữa ESA/SHADE so với SA/ASA có giá trị xác suất $p$-value $< 0.001$, khẳng định tính ưu việt của các thuật toán phát triển mới không phải do yếu tố ngẫu nhiên mà xuất phát từ cơ chế toán tử cải tiến.
  5. Hiệu suất tính toán của mã LPO-V: Thời gian tính toán một cấu hình nạp tải trên máy tính cá nhân tiêu chuẩn chỉ mất xấp xỉ vài chục mili-giây, nhanh hơn đáng kể so với mô-đun CITATION tiêu chuẩn trong khi vẫn duy trì sai số $k_{eff}$ dưới 200 pcm so với MCNP4c.

Implications đa chiều

  • Về mặt học thuật lý thuyết: Mở ra mô hình toán học giải quyết bài toán tối ưu tổ hợp rời rạc nhiều mục tiêu có hàm mục tiêu phi tuyến phức tạp bằng cách kết hợp cơ chế thích ứng lịch sử tham số với ánh xạ vị trí tương đối (RPI).
  • Về mặt kỹ thuật công nghệ hạt nhân: Cung cấp công cụ tính toán nạp tải độc lập, tự chủ cho các kỹ sư vận hành nhà máy điện hạt nhân sử dụng công nghệ VVER (như VVER-1000, VVER-1200, VVER-TOI), cho phép tùy biến nạp tải cả nhiên liệu UO2 thông thường lẫn nhiên liệu tái chế MOX.
  • Về mặt kinh tế - xã hội: Việc tăng $1580\text{ pcm}$ độ phản ứng dư ở đầu chu trình tương đương với việc kéo dài chu kỳ nạp tải thêm hàng chục ngày hiệu dụng công suất định mức, giúp tiết kiệm hàng triệu USD chi phí nhiên liệu hạt nhân cho mỗi chu kỳ vận hành lò phản ứng.

Limitations và Future Research

Nhằm duy trì tính khách quan và chuẩn mực học thuật cao nhất, luận án thẳng thắn thừa nhận các giới hạn kỹ thuật nội tại:

  1. Mô hình tính toán 2 chiều (2D Core Model): Bộ giải LPO-V hiện tại giải phương trình khuếch tán trên mặt phẳng 2D $(x,y)$, chưa mô phỏng chi tiết sự phân bố thông lượng nơ-tron và độ cháy theo phương trục đứng 3D $(z)$.
  2. Giả định xấp xỉ BOC (Beginning of Cycle Proxy): Tối ưu hóa được thực hiện chủ yếu tại thời điểm đầu chu trình nạp tải (BOC) với giả định rằng hệ số nhân vô hạn $k_{\infty}$ suy giảm đơn điệu theo độ cháy, chưa tích hợp trực tiếp chu trình tính toán cháy đầy đủ (Full Core Burnup Depletion Calculation) qua từng bước thời gian.
  3. Bỏ qua hiệu ứng phản hồi nhiệt thủy lực ghép kênh (Thermal-Hydraulic Feedback): Nhiệt độ nhiên liệu và tỷ trọng chất làm mát được cố định theo các trạng thái benchmark, chưa ghép nối song song với bộ giải thủy nhiệt (Thermal-hydraulic subchannel code).

Chương trình nghiên cứu 5–10 năm tới

  • Mở rộng 3D toàn diện: Nâng cấp bộ giải LPO-V lên mô hình khuếch tán 3D lưới lăng trụ tam giác (3D Triangular-Z Mesh).
  • Tối ưu hóa đa chu trình nạp tải (Multi-cycle Optimization): Ghép nối thuật toán SHADE với mô-đun cháy nhiên liệu để theo dõi diễn biến kinh tế nhiên liệu qua 3–5 chu trình liên tiếp.
  • Ghép nối đa vật lý (Multi-physics Coupling): Tích hợp phản hồi nhiệt thủy lực thời gian thực để đánh giá chính xác giới hạn sôi màng nhiệt (DNBR) và phân bố nhiệt độ Doppler.

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

Nghiên cứu tạo ra dấu ấn khoa học và công nghệ rõ nét trong cộng đồng hạt nhân:

  • Tác động học thuật: Kết quả của luận án đã được công bố trên các tạp chí quốc tế chuyên ngành uy tín thuộc danh mục ISI/Scopus (như Annals of Nuclear Energy, Nuclear Engineering and Design), đóng góp giải pháp thuật toán chuẩn mực cho cộng đồng nghiên cứu ICFM thế giới.
  • Chuyển dịch công nghệ công nghiệp: Bộ công cụ LPO-V kết hợp thuật toán SHADE có khả năng ứng dụng trực tiếp vào các dự án vận hành an toàn lò phản ứng hạt nhân VVER-1000 và VVER-1200 trên toàn cầu (Nga, Trung Quốc, Ấn Độ, Đông Âu).
  • Lợi ích kinh tế lượng hóa: Tối ưu hóa cấu hình nạp tải giúp nâng cao độ sâu cháy trung bình của nhiên liệu phóng xạ xả thải, giảm thiểu khối lượng chất thải hạt nhân mức độ cao cần xử lý và lưu giữ địa chất dài hạn.

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

  • Nghiên cứu sinh và Giảng viên ngành Kỹ thuật Hạt nhân: Kế thừa bộ mã nguồn LPO-V và các thuật toán ESA, SHADE để phát triển các đề tài mở rộng về tối ưu hóa vùng hoạt và thiết kế bó nhiên liệu mới.
  • Kỹ sư Quản lý Vùng hoạt (Core Physicists): Sử dụng các cấu hình phân bố nạp tải tối ưu để áp dụng thực tế vào quá trình lập kế hoạch tiếp liệu hàng năm cho các tổ máy điện hạt nhân.
  • Cơ quan Quản lý An toàn Bức xạ và Hạt nhân (VARANS): Sở hữu công cụ tính toán độc lập nhằm thẩm định các báo cáo phân tích an toàn (SAR) và kiểm tra các giới hạn công suất $PPF$ do chủ đầu tư nhà máy đệ trình.

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

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và nó đã mở rộng lý thuyết nền tảng 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 (SHADE) sang Không gian Tổ hợp Rời rạc Đa mục tiêu trong Vật lý Vùng hoạt Lò phản ứng thông qua cơ chế lập chỉ số vị trí tương đối (Relative Position Indexing - RPI). Nghiên cứu đã giải quyết triệt để bài toán biến đổi không gian liên tục $\mathbb{R}^D$ sang hoán vị nguyên của 28 vị trí bó nhiên liệu đối xứng mà vẫn bảo toàn tính thích ứng tự động của bộ nhớ lịch sử $H$ đối với hai tham số $F$ và $CR$.

2. Sự đổi mới về phương pháp luận khi so sánh với ít nhất 2 nghiên cứu quốc tế tiền nhiệm?

So với nghiên cứu của Fadaei et al. (Iran) sử dụng công cụ LONSA (kết hợp mạng nơ-ron và SA truyền thống cho lò Bushehr VVER-1000) vốn dễ gặp sai số xấp xỉ thông lượng cục bộ, luận án tính toán trực tiếp thông lượng nơ-tron bằng phương pháp sai phân hữu hạn (FDM) giải phương trình khuếch tán 4 nhóm năng lượng trên lưới tam giác, đảm bảo tính tất định vật lý chuẩn xác. So với nghiên cứu của nhóm tác giả tại Séc (lò Dukovany và Temelin) sử dụng phần mềm OPTIMAL/SCAM-W với các toán tử di truyền cổ điển, phương pháp ESA của luận án đã kết hợp đột phá toán tử Block Crossover vào quy trình tôi kim, giúp tốc độ hội tụ nhanh hơn gấp nhiều lần và tránh hoàn toàn việc rơi vào cực trị địa phương.

3. Phát hiện gây ngạc nhiên nhất về mặt dữ liệu thực nghiệm là gì?

Phát hiện ấn tượng nhất là khả năng đồng thời tăng $k_{eff}$ lên tới 1580 pcm trong khi vẫn giảm được hệ số đỉnh công suất $PPF$ đi 2.4%. Trong vật lý lò phản ứng, hai mục tiêu này vốn có tính triệt tiêu lẫn nhau (để tăng $k_{eff}$ cần giảm rò rỉ nơ-tron bằng cách dồn nhiên liệu giàu vào tâm, nhưng điều này sẽ làm tăng vọt $PPF$). Thuật toán của luận án đã tìm ra cấu hình sắp xếp xen kẽ thông minh giữa 19 bó UO2 và 9 bó MOX có độ cháy khác nhau, tạo ra phân bố công suất phẳng đều tuyệt đối trên toàn bộ bán kính vùng hoạt.

4. Giao thức tái lập nghiên cứu (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ố chuẩn hóa: bộ thư viện tiết diện 4 nhóm năng lượng, thành phần đồng vị của 7 loại nhiên liệu và vật liệu cấu trúc trong phụ lục, cấu trúc lưới sai phân FDM 24 tam giác/bó, các phương trình vi phân sai phân hóa, cùng toàn bộ lưu đồ thuật toán và siêu tham số của ESA và discrete SHADE ($N_P = 20, H = 6, w_p = 2.5, w_f = 0.05$). Bất kỳ nhà nghiên cứu độc lập nào cũng có thể tái lập chính xác 100% kết quả mô phỏng.

5. Lộ trình nghiên cứu 10 năm tới được phác thảo như thế nào?

Lộ trình hướng tới việc hoàn thiện bộ công cụ tính toán đa vật lý tích hợp: chuyển đổi từ mô hình khuếch tán 2D sang mô hình vận chuyển nơ-tron 3D bậc cao (Transport Theory $S_N$ / Phương pháp đặc trưng MOC), ghép nối song song với mã nhiệt thủy lực đa kênh (Subchannel Thermal-Hydraulics), và mở rộng thuật toán sang các thế hệ lò tương lai như VVER-1200, VVER-TOI và lò phản ứng làm mát bằng nước siêu tới hạn VVER-SCWR (Gen-IV).

Kết luận

Luận án tiến sĩ của tác giả Trần Việt Phú đã đóng góp những giá trị khoa học xuất sắc và toàn diện cho chuyên ngành Vật lý Nguyên tử và Hạt nhân:

  1. Phát triển thành công mã LPO-V: Xây dựng hoàn chỉnh chương trình mô phỏng vật lý vùng hoạt chuyên dụng giải hệ phương trình khuếch tán nơ-tron đa nhóm dừng trên lưới tam giác 2D cho lò phản ứng VVER, đạt tốc độ tính toán vượt trội và sai số đối chuẩn chuẩn xác so với mã Monte Carlo quốc tế MCNP4c.
  2. Đột phá thuật toán ESA: Phát triển thành công phương pháp Mô phỏng tôi kim tiến hóa (ESA) tích hợp toán tử Crossover, nâng tỷ lệ tái tạo cấu hình tối ưu từ 18% (SA) lên 100% với độ tin cậy tuyệt đối.
  3. Tiên phong phát triển Discrete SHADE: Ứng dụng thành công thuật toán tiến hóa vi phân thích ứng dựa trên lịch sử tham số kết hợp kỹ thuật lập chỉ số vị trí tương đối (RPI) cho bài toán hoán vị bó nhiên liệu rời rạc.
  4. Hiệu quả tối ưu hóa vượt bậc: Tìm ra cấu hình nạp tải tối ưu cho lò VVER-1000 nạp 30% MOX giúp tăng $k_{eff}$ đầu chu trình thêm 1580 pcm và giảm hệ số đỉnh công suất $PPF$ 2.4% so với cấu hình chuẩn OECD/NEA.
  5. Chứng thực thống kê vững chắc: Thẩm định toàn diện hiệu năng của các thuật toán thông qua 50 lần chạy độc lập và kiểm định phi tham số Mann-Whitney U-test ($p < 0.05$).
  6. Mở ra các hướng nghiên cứu mới: Đặt nền móng vững chắc cho các nghiên cứu tiếp theo về tối ưu hóa đa chu trình nạp tải 3D, ghép nối đa vật lý nơ-tron - nhiệt thủy lực và làm chủ công nghệ tính toán an toàn hạt nhân tiên tiến.