Tổng quan về luận án

Luận án tiến sĩ chuyên ngành Vật lý lý thuyết và Vật lý toán (Mã số: 62440103) của tác giả Phan Huy Thiện, với đề tài "Nghiên cứu một số bài toán vật lý và cơ học bằng phương pháp đạo hàm trung bình tích phân và phương pháp hỗn hợp", được thực hiện dưới sự hướng dẫn khoa học của GS. Nguyễn Xuân Hãn và GS. Phan Văn Hạp tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội. Công trình đặt trọng tâm vào giải quyết các bài toán biên phức tạp, miền biên phi tuyến hoặc chứa điểm kỳ dị đứt gãy trong vật lý toán, động lực học lò phản ứng hạt nhân và cơ học đàn hồi.

Khoảng trống nghiên cứu (Research Gap) cốt lõi xuất phát từ sự bất cập của các công cụ vi phân cổ điển (Gateaux, Frechet) và phương pháp giải tích - số truyền thống (như sai phân hữu hạn FDM, phần tử hữu hạn FEM) khi đối diện với lớp phương trình vi tích phân có biên hoặc hệ số không trơn. Như văn bản luận án đã chỉ rõ: "Đối với các bài toán cơ học và vật lý có biên phức tạp chứa đựng các điểm đứt gãy thì những phương pháp thông thường khó cho lời giải đảm bảo độ chính xác thậm chí không thể áp dụng được." Tại các điểm gãy hoặc vùng kỳ dị, việc chia lưới giải tích số thường xuyên phải thay đổi, làm tăng đột biến chi phí thuật toán và gây mất ổn định nghiệm.

Để giải quyết triệt để vấn đề này, luận án tập trung vào 4 câu hỏi nghiên cứu (Research Questions) then chốt:

  1. RQ1: Làm thế nào để xây dựng một toán tử khả vi tổng quát cho phép nhất thể hóa các hàm không trơn, chứa điểm góc hoặc bước nhảy gián đoạn trên toàn miền khảo sát?
  2. RQ2: Khả năng ứng dụng toán tử đạo hàm trung bình tích phân (IAD) trong việc tìm nghiệm xấp xỉ của hệ phương trình Wigner – Poisson phi tuyến trong vật lý plasma lượng tử tĩnh điện đạt hiệu quả ra sao?
  3. RQ3: Phương pháp IAD kết hợp đại số ma trận toán tử giải quyết bài toán khuếch tán neutron đa nhóm và neutron trễ trong lõi lò phản ứng hạt nhân tĩnh/động với độ chính xác và tốc độ hội tụ như thế nào so với phương pháp Nodal và Monte Carlo?
  4. RQ4: Cơ chế phối hợp giữa phương pháp IAD, phép biến hình bảo giác và hệ phương trình tích phân kỳ dị để giải bài toán nhiễu xạ sóng đàn hồi trong môi trường khuyết góc nêm đạt độ chuẩn xác ra sao?

Hệ thống giả thuyết nghiên cứu (Hypotheses) được xác lập:

  • H1: Toán tử IAD cho phép bảo toàn tính khả vi tại các điểm góc và điểm gián đoạn hữu hạn hai phía, loại bỏ sự phụ thuộc vào tính trơn lý tưởng của miền biên.
  • H2: Thuật toán hàm mũ ma trận giải tích $\exp(hA)$ kết hợp IAD tối ưu hóa thời gian tính toán và tiết kiệm bộ nhớ máy tính trong mô phỏng thông lượng neutron so với phương pháp thử xác suất Monte Carlo.
  • H3: Việc quy dẫn bài toán biên kỳ dị nêm đàn hồi về phương trình vi tích phân giải bằng IAD cung cấp nghiệm giải tích số ổn định mà không cần tái phân chia lưới phức tạp.

Khung lý thuyết (Theoretical Framework) tích hợp lý thuyết giải tích hàm phi cổ điển (đạo hàm trung bình tích phân Phan Văn Hạp, 1994), lý thuyết vận chuyển Boltzmann, lý thuyết khuếch tán neutron Fermi, cơ học lượng tử không gian pha Wigner (1932) và lý thuyết đàn hồi giải tích phức. Luận án công bố 08 công trình khoa học (gồm 5 bài báo trên tạp chí quốc tế), đóng góp trực tiếp vào thiết kế an toàn lò phản ứng hạt nhân và giải tích nứt gãy cơ học.

Literature Review và Positioning

Lịch sử toán tính toán và vật lý hạt nhân ghi nhận nhiều nhánh nghiên cứu chính trong việc giải các bài toán biên vật lý toán:

[Nhánh 1: Vi phân cổ điển & Lưới số]
[Nhánh 2: Tính toán Lò phản ứng]    [Nhánh 3: Động học Không gian pha]
       [Hệ hình mới của Luận án]

Nhánh thứ nhất tập trung vào phương pháp sai phân hữu hạn (FDM) và phần tử hữu hạn (FEM) cổ điển được phát triển bởi các nhà toán học tính toán như Samarskii (1971), Richardson (1911). Các phương pháp này đòi hỏi hàm nghiệm và biên hình học phải thỏa mãn điều kiện đủ trơn (Sobolev spaces $H^k$ với $k \ge 2$). Khi áp dụng vào môi trường không đồng nhất, biên đứt gãy, sai số tại nút lưới tăng vọt, đòi hỏi kỹ thuật adaptive mesh refinement vô cùng tốn kém tài nguyên máy tính.

Nhánh thứ hai phát triển các công cụ tính toán lò phản ứng hạt nhân dựa trên phương trình vận chuyển Boltzmann và lý thuyết khuếch tán neutron. Nổi bật là phương pháp Nodal đa nhóm của Lawrence (1986) và phương pháp Nodal biến phân (Stacey, 2007; Hennart, 1995), dựa trên phép đồng nhất hóa ô mạng (homogenization) hai giai đoạn. Song song đó, phương pháp Monte Carlo (Spanier & Gelbard, 1969; Brown, 2003) giải phương trình vận chuyển bằng mô phỏng ngẫu nhiên từng lịch sử chuyển động của neutron.

Tranh luận học thuật cốt lõi diễn ra giữa hai trường phái:

  • Trường phái Nodal: Ưu tiên tốc độ bằng cách làm phẳng (trung bình hóa) thông lượng trong từng ô mạng thô, nhưng đánh mất độ chính xác cục bộ tại các vị trí thanh điều khiển, thanh hấp thụ gadolinium hoặc biên phản xạ phức tạp.
  • Trường phái Monte Carlo: Đạt độ chính xác vật lý vi mô cao nhưng tiêu tốn thời gian tính toán khổng lồ, phụ thuộc nặng nề vào việc lấy mẫu thống kê ban đầu và không thể cung cấp bức tranh phân bố liên tục tức thời của toàn hệ thống trong các kịch bản dừng khẩn cấp.

Nhánh thứ ba nghiên cứu phương trình Wigner lượng tử (Wigner, 1932; Markowich et al., 1990) và hệ liên kết Wigner – Poisson trong động học plasma lượng tử, nơi tính phi tuyến và tích phân kì dị gây cản trở các thuật toán cầu phương thông thường. Nhánh thứ tư là lý thuyết bài toán biên đàn hồi nêm kỳ dị (Muskhelishvili, 1953; Bogy, 1972) sử dụng hệ phương trình tích phân kỳ dị.

Vị trí của luận án (Positioning): Luận án định vị tại điểm giao thoa giữa giải tích toán học phi trơn và vật lý kỹ thuật tính toán. Luận án khắc phục hạn chế của FDM/FEM và phương pháp Nodal bằng cách áp dụng lý thuyết Đạo hàm trung bình tích phân (Integrally Averaged Derivative - IAD) do GS. Phan Văn Hạp đề xuất năm 1994, kết hợp với phương pháp toán tử ma trận và phép biến đổi giải tích phức. So với các nghiên cứu quốc tế của Stacey (2007) về Variational Nodal và Brown (2003) về Monte Carlo, phương pháp IAD trong luận án giúp đồng nhất hóa cấu trúc dữ liệu đầu vào không trơn, giảm bậc phương trình vi tích phân thành hệ đại số ma trận giải tích, duy trì tốc độ tính toán siêu nhanh mà vẫn bảo toàn độ chính xác tại các biên gián đoạ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 nền tảng giải tích toán học và vật lý lý thuyết thông qua việc hoàn thiện hóa lý thuyết đạo hàm trung bình tích phân:

  1. Mở rộng lý thuyết vi phân phi Gateaux: Trong không gian hữu hạn chiều $\mathbb{R}^n$, với $F: X \to X$, luận án chuẩn hóa định nghĩa biến phân trung bình tích phân thứ nhất: $$\lim_{t \to 0} \frac{1}{2t} \int_{-t}^{t} [F(x_0 + \tau h) - F(x_0)] d\tau = \delta F(x_0, h)$$ Khi $\delta F(x_0, h) = Ah$ với $A$ là ma trận vuông cấp $n \times n$, hàm $F$ khả vi trung bình tích phân tại $x_0$ ($F'{IAD}(x_0) = A$). Luận án chứng minh chặt chẽ: nếu hàm số có điểm góc hoặc gián đoạn hữu hạn hai phía tại $x_0$, đạo hàm Gateaux không tồn tại nhưng đạo hàm IAD luôn tồn tại và nhận giá trị trung bình đại số của đạo hàm hai phía: $$F'{IAD}(x_0) = \frac{F'(x_0 + 0) + F'(x_0 - 0)}{2}$$ Minh chứng thực nghiệm toán học chỉ ra hàm $F(x) = |x|$ không khả vi Gateaux tại $x_0 = 0$ nhưng tồn tại $F'_{IAD}(0) = 0$.

  2. Chuyển đổi hệ hình toán tử ma trận động học lò phản ứng: Thay vì xấp xỉ sai phân thời gian thông thường dễ dẫn đến mất ổn định số, luận án thiết lập nghiệm giải tích tổng quát dạng đóng cho hệ phương trình vi phân ma trận khuếch tán neutron: $$\frac{d\Phi(t)}{dt} = A\Phi(t) + B \implies \Phi(t_{m+1}) = \exp(hA)\Phi(t_m) + (\exp(hA) - I)A^{-1}B$$ Toán tử hàm mũ ma trận $\exp(hA)$ được khai triển thông qua hệ trị riêng $\lambda_k$ ($k = 0, 1, \dots, K$) và phép khử Gauss, chuyển đổi bài toán vi-tích phân 7 biến phức tạp của phương trình vận chuyển Boltzmann thành bài toán đại số tuyến tính khả thi cao.

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

Khung phân tích của luận án kết hợp đa tầng 3 lý thuyết nền tảng:

  • Lý thuyết IAD: Đóng vai trò bộ lọc trung bình hóa làm trơn dữ liệu đầu vào (tiết diện tán xạ vi phân $\Sigma_s$, tiết diện hấp thụ $\Sigma_a$, tiết diện phân hạch $\nu\Sigma_f$) và xử lý các điểm kỳ dị biên.
  • Lý thuyết toán tử sóng và hàm riêng Laplace: Khai triển thông lượng neutron và trường sóng đàn hồi theo hệ hàm điều hòa trực giao thỏa mãn điều kiện biên ngoại suy $R_{tn} = R_0 + 0.71\lambda_{tr}$.
  • Phép biến hình bảo giác (Conformal Mapping): Ánh xạ miền nêm đàn hồi bị hổng một góc $\alpha$ về nửa mặt phẳng phức tiêu chuẩn, đưa phương trình sóng Helmholtz về hệ phương trình tích phân kỳ dị tuyến tính.
[Dữ liệu & Biên kỳ dị: Góc nêm alpha, Tiết diện hạt nhân gián đoạn]
     [Toán tử IAD: Làm trơn tích phân & Đồng nhất hóa vi phân]
[Chương 2: Lò phản ứng hạt nhân]     [Chương 3: Cơ học đàn hồi]
         [Nghiệm xấp xỉ giải tích số độ chính xác cao]

Điều kiện biên giới hạn (Boundary Conditions) được xác lập nhất quán: trên biên ngoài vùng hoạt, thông lượng neutron $\Phi(R) = 0$; tại tâm đối xứng $r = 0$, nghiệm bị chặn; trên bề mặt tiếp xúc môi trường đàn hồi, các điều kiện ứng suất và chuyển vị được thỏa mãn qua phương trình cân bằng sai số với hàm trọng Sobolev.

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 hình dựa trên hệ hình thực chứng logic - diễn dịch (Positivism / Logic-Deductive Epistemology) đặc trưng của vật lý toán lý thuyết. Thiết kế nghiên cứu tuân thủ mô hình tính toán đa mức (multi-level computational modeling):

  • Mức vi mô: Tương tác tán xạ, hấp thụ, làm chậm neutron thông qua các đại lượng vi phân độ tuổi Fermi $\tau(E)$ và hàm phân bố không gian pha $w(x, v, t)$ của phương trình Liouville lượng tử Wigner.
  • Mức vĩ mô: Trường thông lượng neutron vô hướng $\Phi_g(r, t)$, mật độ công suất vùng lõi, hệ số nhân hiệu dụng $k_{eff}$, và trường chuyển vị sóng đàn hồi $u(x)$.

Quy trình nghiên cứu rigorous

Quy trình giải tích số được thiết kế gồm 4 giai đoạn nghiêm ngặt:

[Giai đoạn 1] Mô hình hóa vật lý & Xử lý số liệu đầu vào qua toán tử IAD
[Giai đoạn 2] Rời rạc hóa không gian - năng lượng (Biến đổi Fourier/Laplace/Sai phân IAD)
[Giai đoạn 3] Rút gọn đại số ma trận trị riêng/vector riêng & Tính exp(hA)
[Giai đoạn 4] Kiểm chứng đối chuẩn (Benchmark) với Monte Carlo và Nodal
  1. Giai đoạn 1 - Mô hình hóa vật lý: Thiết lập phương trình vận chuyển Boltzmann 7 chiều, sau đó rút gọn về hệ phương trình khuếch tán đa nhóm tĩnh và động trong không gian 3 chiều đối xứng.
  2. Giai đoạn 2 - Xử lý số liệu và rời rạc hóa: Sử dụng IAD để trung bình hóa các tham số vật liệu nhảy bậc. Rời rạc hóa biến không gian bằng sai phân kết hợp IAD, chia lưới $K$ khoảng với bước nhảy $h$.
  3. Giai đoạn 3 - Giải thuật ma trận toán tử: Xây dựng ma trận hệ số $A$ kích thước $(G + I) \times (G + I)$ (với $G$ nhóm năng lượng, $I$ nhóm neutron trễ). Tính toán ma trận cơ bản $\Phi(t)$ và ma trận nghịch đảo thông qua hệ nghiệm riêng.
  4. Giai đoạn 4 - Đối chuẩn (Triangulation & Validation): Đối chiếu kết quả thông lượng và giá trị tới hạn với các phương pháp chuẩn quốc tế: phương pháp Monte Carlo, phương pháp Nodal và phương pháp tọa độ rời rạc $S_N$.

Data và phân tích

  • Cấu hình hình học khảo sát: Mô hình lò phản ứng 1 chiều bản phẳng bề dày $H$; hình trụ 2 chiều bán kính $R$, chiều cao $H$; và hình cầu 3 chiều bán kính tới hạn $R_{th} = \pi / \omega$, với bán kính ngoại suy $R_{tn} = R_0 + 0.71\lambda_{tr}$.
  • Thông số năng lượng và vật liệu: Mô hình khuếch tán 2 nhóm năng lượng (nhóm nhanh $\Phi_1$ và nhóm nhiệt $\Phi_2$) kết hợp 1 nhóm neutron trễ tương đương. Khảo sát hệ số nhân hiệu dụng $k_{eff}$ trong dải biến thiên từ trạng thái dưới tới hạn đến sát tới hạn: $k_{eff} \in [0.9000; 0.9999]$.
  • Công cụ giải thuật: Sử dụng thuật toán khử Gauss, phép biến đổi tích phân Fourier ngược và phần mềm tính toán ma trận giải tích để giải hệ phương trình đại số tuyến tính $A\Phi = q$.

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

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

  1. Khả năng ổn định hóa dạng điệu thông lượng neutron sát ngưỡng tới hạn: Khi $k_{eff}$ tiến từ 0.9000 đến 0.9999, phương pháp IAD chỉ ra sự chuyển dịch hình thái thông lượng từ dạng phụ thuộc nguồn cưỡng bức sang dạng hàm điều hòa riêng hình sin ổn định. Tại $k_{eff} = 0.9999$, đồ thị phân bố thông lượng neutron đạt độ trùng khớp tuyệt đối với các phương pháp lặp chuẩn của lý thuyết vận chuyển.
  2. Nghiệm giải tích xấp xỉ cho hệ Wigner – Poisson: Luận án chứng minh định lý về chặn chuẩn không gian pha $\sup_{x, v} |u(x, v, t)| \le C(t)$ cho phương trình Wigner, cho phép áp dụng phương pháp IAD kết hợp xấp xỉ hàm delta qua hệ hạt rời rạc ${x_i(t), v_i(t), w_i(t)}$, biến đổi phương trình tích phân phi tuyến phức tạp thành hệ phương trình vi phân thường đối với trọng số $w_i$.
  3. Loại bỏ kỳ dị trong bài toán nhiễu xạ sóng đàn hồi nêm: Bằng cách kết hợp ánh xạ bảo giác với IAD, luận án đã khử hoàn toàn tính kỳ dị bậc cao tại đỉnh nêm góc $\alpha$, đưa ra đồ thị nhiễu xạ hàm chuyển vị $u(x)$ mượt mà, giải quyết bế tắc mà phương pháp phần tử hữu hạn thông thường gặp phải.
  4. Hiệu năng thuật toán vượt trội: Phương pháp IAD toán tử ma trận cho tốc độ tính toán phân bố thông lượng nhanh hơn đáng kể so với phương pháp Monte Carlo, đồng thời không đòi hỏi bộ nhớ dung lượng lớn như các phương pháp lưới tinh $S_N$/$P_N$.
Tiêu chí so sánh Phương pháp IAD (Luận án) Phương pháp Nodal cổ điển Phương pháp Monte Carlo (MC)
Xử lý biên kỳ dị / Không trơn Tự động làm trơn qua tích phân IAD Cần kỹ thuật homogenization bổ trợ Không áp dụng trực tiếp (lấy mẫu vi mô)
Thời gian tính toán Cực nhanh (nghiệm ma trận đóng) Nhanh (trên lưới thô) Rất chậm (phụ thuộc số hạt mẫu $N$)
Tài nguyên bộ nhớ Rất thấp Thấp đến trung bình Lớn (lưu trữ lịch sử tương tác)
Độ chính xác vùng dị thường Cao (không phát sinh sai số góc) Sai số cao tại thanh hấp thụ Cao nhưng dao động thống kê
Bản chất nghiệm Giải tích số tất định liên tục Rời rạc hóa nút mạng Xác suất ngẫu nhiên

Implications đa chiều

  • Về mặt lý thuyết: Mở ra hướng tiếp cận hoàn chỉnh cho giải tích phi trơn trong vật lý toán, chứng minh rằng toán tử trung bình tích phân là công cụ hữu hiệu thay thế đạo hàm Gateaux trong các môi trường vật chất phân lớp hoặc chứa kỳ dị cấu trúc.
  • Ứng dụng vận hành an toàn lò hạt nhân: Cung cấp thuật toán tính toán tức thời (real-time calculation) phân bố thông lượng và $k_{eff}$. Như luận án nhấn mạnh bối cảnh từ các thảm họa Chernobyl và Fukushima: "nhiệm vụ của người tính toán lò cần có một phương pháp tính nhanh nhất cho kết quả chính xác để đoán biết sự cố, và điều khiển lò ngay lập tức khắc phục sự cố ở miền có thông lượng cao để hệ số nhân hiệu dụng trở về ngưỡng an toàn cho phép".
  • Kỹ thuật cơ học công trình: Cung cấp công cụ giải chính xác ứng suất động học và trường tán xạ sóng xung quanh các vết nứt, khe hổng góc nêm trong môi trường đất đá và vật liệu composite.

Limitations và Future Research

  1. Giới hạn hình học: Luận án mới tập trung giải quyết chi tiết các mô hình đối xứng hình học cơ bản (tấm phẳng 1D, hình trụ 2D, hình cầu 3D). Các cấu hình bất đối xứng phức tạp 3 chiều toàn vùng lõi cần tiếp tục mở rộng.
  2. Giả định xấp xỉ năng lượng: Mô hình động học dừng lại ở phép gần đúng 2 nhóm neutron kết hợp 1 nhóm neutron trễ trung bình, chưa mở rộng trực tiếp sang phổ liên tục hàng trăm nhóm năng lượng.
  3. Hiệu ứng phản hồi nhiệt - thủy lực (Thermal-Hydraulic Feedback): Nghiên cứu tập trung vào trạng thái nơtron tĩnh và động học tuyến tính, chưa ghép nối hoàn chỉnh với phương trình truyền nhiệt dòng chảy chất làm mát phi tuyến.

Chương trình nghiên cứu tương lai (Future Agenda):

  • Phát triển phần mềm tính toán tự động 3D tích hợp IAD cho lõi lò phản ứng thế hệ mới (Gen-IV VHTR, SMR).
  • Ghép nối đa vật lý (Multiphysics Coupling): nơtron học - nhiệt thủy lực - độ bền cơ học nứt vỡ thanh nhiên liệu bằng phương pháp IAD hỗn hợp.
  • Mở rộng lý thuyết IAD cho các bài toán truyền hạt trong môi trường tán xạ bất đẳng hướng cao ($P_N$ bậc cao).

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

Công trình tạo ra những ảnh hưởng sâu sắc trên cả bình diện học thuật lẫn thực tiễn:

  • Tác động học thuật: Định hình một công cụ tính toán mới trong trường phái vật lý toán Việt Nam do GS. Phan Văn Hạp và GS. Nguyễn Xuân Hãn khởi xướng, được minh chứng qua 08 công trình công bố uy tín trong và ngoài nước.
  • Tác động công nghiệp năng lượng: Đóng góp trực tiếp vào cơ sở dữ liệu tính toán, thiết kế và thẩm định an toàn lò phản ứng hạt nhân nghiên cứu và nhà máy điện hạt nhân.
  • Đóng góp quốc tế: Cung cấp giải pháp thay thế hiệu quả cho các bài toán benchmark quốc tế về phương trình Wigner-Poisson và nhiễu xạ sóng đàn hồi trên miền nêm kỳ dị.

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

  • Nghiên cứu sinh & Nhà vật lý lý thuyết: Tiếp cận phương pháp toán tử IAD để giải các hệ phương trình đạo hàm riêng phi tuyến chứa điểm kỳ dị mà các giải tích cổ điển không xử lý được.
  • Kỹ sư thiết kế và vận hành lò phản ứng: Sở hữu thuật toán tính nhanh phân bố thông lượng neutron và hệ số nhân $k_{eff}$ phục vụ kiểm soát an toàn vùng hoạt thời gian thực.
  • Chuyên gia cơ học rạn nứt & Địa kỹ thuật: Ứng dụng mô hình nhiễu xạ sóng IAD để đánh giá trạng thái ứng suất tập trung tại đáy vết nứt công trình.
  • Cơ quan quản lý an toàn bức xạ và hạt nhân: Sử dụng các mô hình đánh giá độc lập để thẩm định an toàn thanh nhiên liệu và vùng hoạt lò phản ứng.

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ì?
    Trả lời: Đó là việc hoàn thiện và ứng dụng thành công lý thuyết Đạo hàm trung bình tích phân (IAD) của GS. Phan Văn Hạp vào vật lý toán. Luận án đã mở rộng khái niệm vi phân cho lớp hàm không trơn, chứng minh sự tồn tại đạo hàm IAD tại các điểm gãy, điểm gián đoạn và thiết lập công thức nghiệm ma trận mũ giải tích $\exp(hA)$ đóng cho hệ phương trình động học lò phản ứng.

  2. Đột phá phương pháp luận của luận án so với các nghiên cứu quốc tế trước đây?
    Trả lời: So với phương pháp Nodal (Lawrence, 1986) vốn làm mịn nhân tạo gây mất mát thông tin tại biên cục bộ và phương pháp Monte Carlo (Brown, 2003) tiêu tốn thời gian tính toán khổng lồ, phương pháp IAD của luận án thực hiện trung bình hóa tích phân ngay từ bước xử lý số liệu đầu vào, cho phép thu được nghiệm liên tục giải tích số chính xác cao mà không cần chia lưới tinh hay lấy mẫu ngẫu nhiên lặp lại.

  3. Phát hiện bất ngờ nhất về mặt dữ liệu tính toán là gì?
    Trả lời: Đó là sự ổn định hình thái thông lượng neutron khi hệ số nhân hiệu dụng tiến sát ngưỡng tới hạn $k_{eff} = 0.9999$. Tại trạng thái này, đồ thị thông lượng giải bằng IAD triệt tiêu hoàn toàn sự phụ thuộc vào dạng nhảy bậc của nguồn ngoài $q(x)$, tự động hội tụ về dạng phân bố chuẩn của hàm riêng Laplace.

  4. Luận án có cung cấp quy trình tái lập (Replication Protocol) không?
    Trả lời: Hoàn toàn có. Luận án cung cấp chi tiết hệ phương trình đại số ma trận, công thức tính toán hệ số suy biến ma trận qua phép khử Gauss, biểu thức giải tích thông số hình học $B^2$ cho 3 dạng hình học (cầu, trụ, phẳng), và các bước ánh xạ bảo giác cho bài toán nêm đàn hồi.

  5. Chương trình nghiên cứu 10 năm tiếp theo được định hình ra sao?
    Trả lời: Tập trung vào 3 trục: (1) Xây dựng mã nguồn tính toán song song 3D toàn vùng lõi cho lò phản ứng module nhỏ (SMR); (2) Phát triển gói phần mềm giải tích nứt gãy động lực học đàn - dẻo sử dụng IAD; (3) Mở rộng giải hệ phương trình Wigner lượng tử nhiều chiều trong công nghệ bán dẫn nano.

Kết luận

Luận án tiến sĩ của Phan Huy Thiện là một công trình nghiên cứu công phu, nghiêm túc và có hàm lượng khoa học cao, thể hiện qua 5 đóng góp cốt lõi:

  1. Xác lập hoàn chỉnh cơ sở giải tích toán học của phương pháp Đạo hàm trung bình tích phân (IAD) cho lớp bài toán biên không trơn và kỳ dị.
  2. Xây dựng thành công thuật toán xấp xỉ số cho hệ phương trình Wigner – Poisson phi tuyến trong plasma lượng tử.
  3. Thiết lập mô hình toán tử ma trận giải tích $\exp(hA)$ giải phương trình khuếch tán neutron đa nhóm và neutron trễ trong lõi lò phản ứng hạt nhân.
  4. Giải quyết trọn vẹn bài toán nhiễu xạ sóng đàn hồi trong môi trường khuyết góc nêm bằng cách kết hợp IAD và biến hình bảo giác.
  5. Chứng minh tính ưu việt vượt trội của IAD về tốc độ và độ chính xác so với phương pháp Nodal và Monte Carlo thông qua 08 công trình công bố uy tín.

Công trình không chỉ mở ra hệ hình tính toán mới cho chuyên ngành Vật lý lý thuyết và Vật lý toán, mà còn mang giá trị thực tiễn to lớn, đặt nền móng giải thuật vững chắc cho công tác tính toán và thiết kế an toàn lò phản ứng hạt nhân trong tương lai.