Tổng quan luận án

Bản luận án tiến sĩ 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 bởi nghiên cứu sinh Phan Huy Thiệ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. Luận án được bảo vệ tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội vào năm 2016, thuộc chuyên ngành Vật lý lý thuyết và vật lý toán (mã số: 62440103).

Tính cấp thiết và khoảng trống nghiên cứu

Trong nghiên cứu vật lý toán và cơ học lý thuyết, nhiều hiện tượng tự nhiên vi mô và vĩ mô được mô hình hóa dưới dạng các phương trình vi phân thường, vi phân đạo hàm riêng hoặc phương trình vi tích phân với các điều kiện biên tĩnh, động hoặc hỗn hợp. Các phương pháp toán học giải gần đúng truyền thống thường phát huy hiệu quả khi hệ số, nghiệm hoặc miền biên có tính chất trơn và liên tục. Tuy nhiên, trong thực tiễn kỹ thuật và vật lý (chẳng hạn như quá trình khuếch tán neutron trong lò phản ứng hạt nhân, truyền nhiệt, thủy động học hay bài toán cơ học đàn hồi), các bài toán thường chứa đựng hệ số, nghiệm hoặc biên không trơn, xuất hiện các điểm đứt gãy hoặc miền kỳ dị.

Đối với các bài toán có biên phức tạp chứa điểm đứt gãy, các phương pháp cổ điển như sai phân hữu hạn (FDM) vấp phải trở ngại lớn do thuật toán đòi hỏi phải liên tục thay đổi phân vùng chia lưới tại các vùng kỳ dị, gây tốn kém tài nguyên tính toán và khó cài đặt trên máy tính. Trong lĩnh vực vật lý lò phản ứng hạt nhân, việc tính toán phân bố thông lượng neutron và hệ số nhân hiệu dụng ($k_{eff}$) đóng vai trò then chốt nhằm bảo đảm an toàn vận hành, đặc biệt trong các tình huống ứng phó sự cố (như bài học từ sự cố Chernobyl hay Fukushima). Trong khi đó, ở lĩnh vực cơ học đàn hồi, bài toán nhiễu xạ sóng trong môi trường đàn hồi có biên cứng bị hổng một góc $\alpha$ vẫn chưa có lời giải hoàn chỉnh bằng các phương pháp giải tích số tiêu chuẩn. Khoảng trống nghiên cứu được tác giả xác định là sự thiếu vắng một công cụ toán học có khả năng nhất thể hóa các khái niệm vi phân và đạo hàm trên toàn miền (kể cả biên kỳ dị yếu), giúp xử lý số liệu đầu vào không trơn và giải gần đúng các phương trình vật lý toán một cách ổn định, giảm khối lượng tính toán.

Mục tiêu nghiên cứu

Mục tiêu nghiên cứu tổng quát và các mục tiêu cụ thể được luận án xác định bao gồm:

  1. Đề xuất việc sử dụng phương pháp đạo hàm trung bình tích phân (Integrally Averaged Derivative - IAD) phối hợp với các phương pháp số khác (phương pháp lặp, phương pháp hàm riêng và trị riêng, phương pháp sai phân, phương pháp tham số liên tục, phương pháp biến đổi Fourier và Laplace) để tìm nghiệm gần đúng cho các bài toán vật lý và cơ học có miền biên chứa điểm đứt gãy hoặc kỳ dị.
  2. Ứng dụng phương pháp IAD để khảo sát nghiệm xấp xỉ của phương trình Wigner và hệ liên kết Wigner – Poisson trong vật lý plasma lượng tử tĩnh điện.
  3. Áp dụng phương pháp IAD kết hợp toán tử ma trận để giải phương trình khuếch tán neutron đa nhóm trong lò phản ứng hạt nhân ở trạng thái dừng, xác định thông lượng neutron và các điều kiện tới hạn; so sánh hiệu quả tính toán với phương pháp Nodal và phương pháp Monte Carlo.
  4. Áp dụng phương pháp IAD giải bài toán cơ học về nhiễu xạ sóng trong môi trường đàn hồi thuần nhất có biên cứng bị hổng một góc $\alpha$ thông qua việc đưa bài toán về hệ phương trình vi tích phân kỳ dị.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Các phương trình vật lý toán mô tả hiện tượng vật lý vi mô và cơ học vĩ mô, cụ thể gồm: phương trình Wigner dạng elliptic trong plasma lượng tử; phương trình vận chuyển Boltzmann và phương trình khuếch tán neutron (đơn nhóm, đa nhóm, có xét đến neutron trễ); phương trình sóng Helmholtz và phương trình vi tích phân mô tả nhiễu xạ sóng trong môi trường đàn hồi.
  • Phạm vi nghiên cứu:
    • Về mặt thời gian và trạng thái: Luận án tập trung chủ yếu vào bài toán dừng (trạng thái tĩnh, không phụ thuộc thời gian hoặc đạo hàm theo thời gian bằng 0), khảo sát sự phân bố neutron ở nhiệt độ danh định $20^\circ\text{C}$ đối với neutron nhiệt.
    • Về mặt không gian và hình học: Các mô hình hình học đối xứng tiêu chuẩn một chiều và hai chiều bao gồm tấm phẳng vô hạn độ dày $H$, hình trụ hữu hạn/vô hạn bán kính $R$, hình cầu bán kính $R$, và mặt phẳng đàn hồi thuần nhất bị cắt một góc hình nêm $\alpha$.

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu trước được tác giả điểm lại

Trong phần tổng quan, tác giả hệ thống hóa các hướng tiếp cận tính toán lý thuyết và phương pháp số từ các công trình đã công bố trong và ngoài nước:

  • Phương pháp đạo hàm trung bình tích phân: Xuất phát từ công trình của GS. Phan Văn Hạp công bố từ năm 1994 về việc định nghĩa đạo hàm trung bình tích phân và áp dụng giải các bài toán biên có nghiệm không trơn, bài toán tích phân kỳ dị tuyến tính và phi tuyến.
  • Lý thuyết plasma lượng tử và phương trình Wigner: Phương trình Liouville lượng tử do Wigner đưa ra năm 1932 từ phương trình Schrödinger, cùng các nghiên cứu về hệ Wigner – Poisson trong môi trường plasma tĩnh điện dựa trên lý thuyết trường trung bình.
  • Lý thuyết vận chuyển và khuếch tán neutron: Các mô hình lý thuyết kinh điển về lò phản ứng hạt nhân của Fermi (lý thuyết làm chậm liên tục theo độ tuổi), lý thuyết vận chuyển hạt của Boltzmann, và các phép gần đúng khuếch tán đơn nhóm, đa nhóm.
  • Các phương pháp số giải phương trình lò phản ứng:
    • Phương pháp sai phân hữu hạn (FDM) và phần tử hữu hạn (FEM): Thường áp dụng cho các hệ tọa độ $X-Y$ hoặc $R-Z$, gặp khó khăn khi hình học phức tạp hoặc khi xuất hiện điểm kỳ dị của thông lượng neutron.
    • Phương pháp rời rạc hóa phương trình vận chuyển: Phương pháp $S_N$ (rời rạc theo hướng) và phương pháp $P_N$ (rời rạc theo tọa độ không gian).
    • Phương pháp Nodal và Nodal biến phân: Kỹ thuật tính toán hai giai đoạn (đồng nhất hóa ô mạng thanh nhiên liệu ở cấp vi mô bằng lý thuyết vận chuyển, sau đó giải phương trình khuếch tán trên lưới thưa toàn vùng lõi). Kỹ thuật này đã được phát triển để xử lý các lõi lò không đồng nhất thế hệ mới nhưng vẫn tồn tại sai số tại các ô lưới có thanh hấp thụ mạnh (như thanh gadolinium).
    • Phương pháp ngẫu nhiên Monte Carlo: Mô phỏng hành vi chuyển động, tán xạ, hấp thụ của từng neutron theo phép thử xác suất. Phương pháp này đòi hỏi thời gian tính toán rất lớn, bộ nhớ cao và kết quả phụ thuộc vào mẫu ngẫu nhiên ban đầu, khó quan sát phân bố thông lượng toàn hệ thống theo thời gian cố định.
  • Bài toán nhiễu xạ sóng cơ học: Các công trình nghiên cứu cơ học biến dạng giải bài toán nhiễu xạ sóng đàn hồi có khuyết tật hình học, trong đó tác giả trích dẫn trực tiếp để đối chiếu số liệu là công trình của tác giả [95] và tài liệu [130].

Khoảng trống luận án lựa chọn giải quyết

Luận án chỉ ra rằng các phương pháp hiện hành chia thành hai thái cực: phương pháp giải tích thuần túy không xử lý được biên phức tạp hoặc số liệu đầu vào đứt gãy, trong khi các phương pháp số thuần túy (FDM, Nodal biến phân, Monte Carlo) đòi hỏi tài nguyên máy tính lớn hoặc gặp bất ổn định nghiệm tại điểm kỳ dị. Luận án lựa chọn giải quyết khoảng trống này bằng cách xây dựng phương pháp IAD kết hợp: đưa các bài toán vi phân trên miền có kỳ dị yếu về một mô hình toán học thống nhất, làm trơn số liệu đầu vào phi tuyến bằng trung bình hóa tích phân, và chuyển phương trình vi phân nhiều nhóm về hệ ma trận giải tích giúp tối ưu hóa thời gian tính toán và bộ nhớ máy tính.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết và khái niệm sử dụng

  1. Khái niệm Đạo hàm trung bình tích phân (IAD): Cho $X$ là không gian hữu hạn chiều $\mathbb{R}^n$, hàm $F: X \to X$, biến $x_0 \in X$ (thực hoặc phức) và hướng $h \in X$. Biến phân trung bình theo tích phân thứ nhất của hàm $F$ tại $x_0$ được định nghĩa: $$\delta F(x_0, h) = \lim_{t \to 0} \frac{1}{2t} \int_{-t}^t [F(x_0 + \tau h) - F(x_0)] d\tau$$ Nếu tồn tại ma trận vuông cấp $n \times n$ ký hiệu là $A$ sao cho $\delta F(x_0, h) = Ah$, thì hàm $F$ được gọi là khả vi trung bình theo tích phân tại $x_0$, ký hiệu $F'(x_0, h) = Ah$. Biểu thức mở rộng cho phép lấy tích phân lệch: $$\delta F(x_0, h) = \lim_{\epsilon \to 0, \tau \to 0} \frac{1}{2(\epsilon + \tau)} \int_{-\epsilon}^\tau [F(x_0 + t h) - F(x_0)] dt$$
  2. Lý thuyết vận chuyển Boltzmann và phép gần đúng khuếch tán: Sự cân bằng neutron trong lõi lò phản ứng được thiết lập dựa trên phương trình cân bằng số neutron có năng lượng từ $E$ đến $E + dE$: $$\frac{1}{v}\frac{\partial \Phi(r, E, t)}{\partial t} = D(E)\Delta \Phi(r, E, t) - \Sigma_a(E)\Phi(r, E, t) + Q_s + Q_f$$ Trong đó, tốc độ gia tăng do tán xạ $Q_s = \int_E \Sigma_s(E' \to E)\Phi(r, E', t)dE'$, và tốc độ sinh do phân hạch $Q_f = \chi(E) \int_0^\infty \nu \Sigma_f(E')\Phi(r, E', t)dE'$.
  3. Lý thuyết làm chậm Fermi (Fermi Age Theory): Mô tả quá trình neutron mất năng lượng liên tục qua đại lượng lethargy $u$. Mật độ dòng làm chậm $j(r, u, t) \approx \xi \Sigma_s(u) \Phi(r, u, t)$ thỏa mãn phương trình $\Delta j - \frac{\partial j}{\partial \tau} + Q = 0$, với $\tau$ là độ tuổi neutron ($\text{cm}^2$). Đại lượng $\tau$ liên hệ với bình phương trung bình khoảng cách dịch chuyển $r_0^2 = 6\tau$.
  4. Phương pháp giải bài toán nhiễu xạ bằng biến hình bảo giác: Sử dụng lý thuyết hàm biến phức và các phép biến đổi bảo giác để chuyển đổi miền biên kỳ dị hình nêm về biên chuẩn, đưa phương trình sóng về hệ phương trình vi tích phân kỳ dị.

Phương pháp nghiên cứu thực tế

  • Phương pháp định tính: Khảo sát tính chất giải tích của toán tử vi phân và toán tử tích phân; chứng minh định lý về tính bị chặn của nghiệm trong hệ phương trình Wigner – Poisson; phân tích các tính chất của đạo hàm IAD đối với hàm có điểm góc, hàm gián đoạn nhảy bậc và hàm có đạo hàm Gateaux.
  • Phương pháp định lượng và kỹ thuật tính toán:
    • Rời rạc hóa sai phân kết hợp IAD: Rời rạc hóa biến không gian trên lưới sai phân hữu hạn, kết hợp trung bình hóa tích phân tại các mắt lưới để làm trơn số liệu đầu vào.
    • Kỹ thuật toán tử ma trận giải tích: Đưa hệ phương trình vi phân khuếch tán đa nhóm theo thời gian về phương trình ma trận $\frac{d\Phi(t)}{dt} = A\Phi(t) + B$. Nghiệm tổng quát có dạng ma trận hàm mũ $\Phi(t_{m+1}) = \exp(hA)\Phi(t_m) + (\exp(hA) - I)A^{-1}B$.
    • Khai triển hàm mũ ma trận: Biểu diễn $\exp(hA) = \sum_{k=0}^K \alpha_k (hA)^k$ với $K = G + I - 1$ ($G$ là số nhóm năng lượng, $I$ là số nhóm neutron trễ). Dùng phép khử Gauss giải hệ phương trình xác định các hệ số $\alpha_k$ từ các trị riêng $\lambda_k$ của ma trận $A$.
    • Phương pháp cân bằng sai số (Method of Error Balancing): Tìm nghiệm xấp xỉ $U^* = \sum_{i=1}^n c_i \phi_i$ thỏa mãn hệ hàm cơ sở đầy đủ, cực tiểu hóa phiếm hàm sai số trên miền và trên biên.

Nguồn tư liệu và dữ liệu nghiên cứu

Dữ liệu đầu vào phục vụ mô phỏng số bao gồm: các hằng số tiết diện vi mô và vĩ mô của neutron ($\Sigma_t, \Sigma_s, \Sigma_a, \nu\Sigma_f$), hệ số khuếch tán $D$, phổ phân hạch $\chi$, các giá trị trị riêng của hàm Bessel bậc không $J_0(\xi)$, và các thông số vật liệu của mô hình benchmark lò phản ứng hạt nhân ở $20^\circ\text{C}$.

Bảng dưới đây tổng hợp các nghiệm hàm riêng và thông số hình học cho các cấu hình lò phản ứng được dẫn ra trong luận án:

Dạng hình học lò phản ứng Dạng toán tử Laplace $\Delta \Phi$ Nghiệm giải tích của hàm riêng $\Phi$ Thông số hình học $B^2$ Kích thước tới hạn ($B = \omega$)
Tấm phẳng vô hạn (độ dày $H$) $\frac{d^2\Phi}{dx^2}$ $\cos(Bx)$ $\left(\frac{\pi}{H}\right)^2$ $H_{th} = \frac{\pi}{\omega} = \frac{\pi L}{\sqrt{k_\infty - 1}}$
Hình trụ vô hạn (bán kính $R$) $\frac{1}{r}\frac{d}{dr}\left(r\frac{d\Phi}{dr}\right)$ $J_0(Br)$ $\left(\frac{2,405}{R}\right)^2$ $R_{th} = \frac{2,405}{\omega} = \frac{2,405 L}{\sqrt{k_\infty - 1}}$
Hình trụ hữu hạn (bán kính $R$, chiều cao $H$) $\frac{1}{r}\frac{\partial}{\partial r}\left(r\frac{\partial\Phi}{\partial r}\right) + \frac{\partial^2\Phi}{\partial z^2}$ $J_0(B_r r)\cos(B_z z)$ $\left(\frac{2,405}{R}\right)^2 + \left(\frac{\pi}{H}\right)^2$ Xác định qua hệ thức $\left(\frac{2,405}{R}\right)^2 + \left(\frac{\pi}{H}\right)^2 = \omega^2$
Hình cầu (bán kính $R$) $\frac{1}{r^2}\frac{d}{dr}\left(r^2\frac{d\Phi}{dr}\right)$ $\frac{\sin(Br)}{r}$ $\left(\frac{\pi}{R}\right)^2$ $R_{th} = \frac{\pi}{\omega} = \pi \sqrt{\frac{L^2}{k_\infty - 1}}$

Nội dung chính theo từng chương

CHƯƠNG 1 ĐẠO HÀM TRUNG BÌNH TÍCH PHÂN

Chương 1 thiết lập nền tảng toán học cho toàn bộ luận án, tập trung vào định nghĩa, tính chất của đạo hàm IAD, các khái niệm nghiệm xấp xỉ, phương pháp cân bằng sai số và ứng dụng giải phương trình Wigner.

Các tính chất cốt lõi của đạo hàm IAD được chứng minh:

  • Tính tuyến tính: $(F + G)'(x_0) = F'(x_0) + G'(x_0)$ và $(\alpha F)'(x_0) = \alpha F'(x_0)$.
  • Mối liên hệ với đạo hàm Gateaux: Nếu hàm $F$ có đạo hàm Gateaux tại $x_0$ thì cũng có đạo hàm trung bình tích phân tại đó và hai đạo hàm trùng nhau ($F'(x_0) = F'_G(x_0)$).
  • Khả năng xử lý điểm không khả vi thông thường: Tồn tại những hàm không có đạo hàm Gateaux nhưng tồn tại đạo hàm IAD. Ví dụ tiêu biểu là hàm trị tuyệt đối $F(x) = |x|$ tại điểm $x_0 = 0$: $$F'(0) = \lim_{t \to 0} \frac{1}{2t} \int_{-t}^t \text{sign}(\tau) d\tau = 0$$
  • Hàm có điểm góc hoặc gián đoạn loại 1: Nếu $F$ có đạo hàm hữu hạn hai phía tại $x_0$, thì đạo hàm trung bình tích phân tồn tại và bằng trung bình cộng của đạo hàm hai phía: $$F'(x_0) = \frac{F'(x_0 + 0) + F'(x_0 - 0)}{2}$$ Tính chất này cho phép đưa các phương trình có hàm số nhảy bậc hoặc hệ số đứt gãy về dạng khả vi trên toàn miền mà không làm mất tính chất vật lý của bài toán.

Ứng dụng vào phương trình Wigner: Phương trình Wigner trong không gian pha có dạng: $$\frac{\partial w}{\partial t} + v \nabla_x w - \nabla_x E w_v = \theta_v w$$ Áp dụng phép biến đổi Fourier $u = \mathcal{F}w$, tác giả đưa phương trình về dạng: $$\frac{\partial u}{\partial t} - i \text{div}_\xi \nabla_x u - E v u = \delta_v(x, \xi, t)u$$ Sử dụng phương pháp IAD kết hợp công thức cầu phương cho tập hợp các hạt rời rạc tại tọa độ $(x_i(t), v_i(t))$ với trọng số $w_i(t)$, tác giả thu được hệ phương trình vi phân đại số xác định trọng số: $$\frac{dw_i}{dt} = \sum_j w_j(x_i - x_j) \theta(x_j, v_i - v_j)w_j$$ Đồng thời, luận án chứng minh định lý về tính bị chặn của nghiệm đối với toán tử vận chuyển kết hợp phương trình Poisson $V = n - V$, tạo cơ sở toán học để áp dụng IAD cho các phương trình vận chuyển ở các chương sau.

CHƯƠNG II TÍNH TOÁN THÔNG LƯỢNG NEUTRON TRONG LÒ PHẢN ỨNG HẠT NHÂN

Chương 2 là nội dung ứng dụng trọng tâm của luận án, giải quyết bài toán cân bằng neutron trong lõi lò phản ứng hạt nhân.

Phương trình khuếch tán neutron nhiệt một nhóm: Trong trạng thái tĩnh ($\partial \Phi / \partial t = 0$), phương trình khuếch tán neutron nhiệt rút gọn về dạng: $$D\Delta \Phi - \Sigma_a \Phi + \nu \Sigma_f \Phi = 0 \iff \Delta \Phi + \omega^2 \Phi = 0$$ Với $\omega^2 = \frac{k_\infty - 1}{L^2}$ là thông số vật liệu lò phản ứng ($L = \sqrt{D/\Sigma_a}$ là độ dài khuếch tán, $k_\infty = \nu \Sigma_f / \Sigma_a$ là hệ số nhân vô hạn). Đẳng thức $\omega^2 = B^2$ là điều kiện tới hạn để duy trì phản ứng dây chuyền tự bão hòa.

Phép gần đúng theo độ tuổi Fermi: Đối với dải năng lượng neutron liên tục từ phân hạch đến nhiệt hóa, lý thuyết Fermi đưa ra phương trình trạng thái tới hạn: $$k_{eff} = k_\infty \frac{\exp(-B^2 \tau)}{1 + B^2 L^2} = 1$$ Trong đó, $\exp(-B^2 \tau)$ biểu diễn xác suất tránh rò rỉ của neutron nhanh trong quá trình làm chậm có độ tuổi $\tau$, và $1/(1 + B^2 L^2)$ là xác suất tránh rò rỉ của neutron nhiệt. Đối với các lò phản ứng có kích thước lớn, biểu thức được xấp xỉ thông qua diện tích vận chuyển $M^2 = L^2 + \tau$: $$k_{eff} = \frac{k_\infty}{1 + B^2 M^2}$$

Phương pháp toán tử ma trận IAD cho mô hình đa nhóm và neutron trễ: Tác giả thiết lập hệ phương trình động học đa nhóm có xét đến $I$ nhóm neutron trễ: $$\frac{1}{v_g}\frac{\partial \Phi_g(r, t)}{\partial t} = \nabla D_g \nabla \Phi_g - \Sigma_{ag}\Phi_g - \sum_{g' \neq g} \Sigma_{s}^{g \to g'}\Phi_g + \sum_{g' \neq g} \Sigma_{s}^{g' \to g}\Phi_{g'} + (1 - \beta)\chi_g \sum_{g'} \nu \Sigma_{fg'}\Phi_{g'} + \sum_{i=1}^I \lambda_i \chi_{gi} C_i$$ Bằng việc sử dụng sai phân hữu hạn và trung bình hóa tích phân theo IAD, hệ phương trình được chuyển về phương trình ma trận $\frac{d\Phi(t)}{dt} = A\Phi(t) + B$. Nghiệm ma trận hàm mũ được tính chính xác bằng phép khai triển đa thức bậc $K = G + I - 1$: $$\exp(hA) = \sum_{k=0}^K \alpha_k (hA)^k$$ Các hệ số $\alpha_k$ được xác định qua việc giải hệ phương trình đại số với các trị riêng $\lambda_k$ bằng thuật toán khử Gauss.

Khảo sát số học mô hình 1 chiều có nguồn: Tác giả khảo sát lò phản ứng 1 chiều độ dày $a$, có nguồn phân bố dạng hàm bước nhảy ở tâm: $$q(x) = \begin{cases} 1 & \text{khi } -a/2 \le x \le a/2 \ 0 & \text{khi } |x| > a/2 \end{cases}$$ Rời rạc hóa miền không gian thành $K$ khoảng với bước lưới $h$, hệ đại số ma trận $A\Phi = q$ được giải để thu được phân bố thông lượng với các giá trị $k_{eff}$ thay đổi từ $0,9$ đến $0,9999$. Kết quả chỉ ra rằng khi hệ ở xa trạng thái tới hạn ($k_{eff} = 0,9$), thông lượng chịu ảnh hưởng áp đảo bởi hình dạng nguồn $q(x)$. Khi hệ tiến gần đến ngưỡng tới hạn ($k_{eff} = 0,9999$), dạng hình học của thông lượng neutron trở nên ổn định, có dạng hình cosin trơn tru và hoàn toàn tương thích với kết quả từ các phương pháp lặp $S_N$ và phương pháp Monte Carlo, nhưng với thời gian tính toán nhanh hơn và không đòi hỏi xử lý ma trận thưa phức tạp.

CHƯƠNG III PHƯƠNG PHÁP GIẢI BÀI TOÁN NHIỄU XẠ SÓNG TRONG MÔI TRƯỜNG ĐÀN HỒI

Chương 3 áp dụng phương pháp IAD vào cơ học vật rắn biến dạng, cụ thể là bài toán nhiễu xạ sóng đàn hồi trong mặt phẳng $(r, \theta)$ bị hổng một góc hình nêm $\alpha$.

  • Mô hình toán học: Mặt phẳng đàn hồi chứa góc khuyết tạo ra điều kiện biên không trơn, tại đó các đạo hàm không gian thông thường của trường chuyển vị bị gián đoạn hoặc tiến tới vô hạn (kỳ dị biên).
  • Kỹ thuật biến hình bảo giác: Tác giả sử dụng phép biến hình bảo giác trong giải tích phức để ánh xạ miền góc nêm kỳ dị về miền nửa mặt phẳng chuẩn, từ đó thiết lập phương trình sóng dưới dạng hệ phương trình vi tích phân kỳ dị đối với hàm chuyển vị $u(x)$.
  • Xử lý số bằng IAD: Nhờ định nghĩa đạo hàm trung bình tích phân, các tích phân kỳ dị yếu trên biên được tích phân trực tiếp mà không cần loại bỏ điểm kỳ dị nhân tạo. Hệ phương trình được giải xấp xỉ bằng phương pháp số.
  • Đối chiếu kết quả: Đồ thị phân bố biên độ sóng nhiễu xạ $u(x)$ thu được bằng phương pháp IAD được so sánh với nghiệm trong tài liệu tham khảo [95]. Sự trùng khớp giữa hai đường cong khẳng định tính khả thi của phương pháp IAD trong việc giải quyết các bài toán cơ học có biên đứt gãy mà không cần dùng đến các lưới phần tử hữu hạn quá mịn tại đỉnh nêm.

Kết quả và những đóng góp mới

Đóng góp mới về mặt lý luận và phương pháp

  1. Hoàn thiện công cụ toán học IAD: Luận án đã hệ thống hóa và cụ thể hóa định nghĩa đạo hàm trung bình tích phân (IAD) cho các bài toán biên vật lý toán, chứng minh mối quan hệ giữa đạo hàm IAD với đạo hàm Gateaux và đạo hàm hai phía. Phương pháp này cho phép nhất thể hóa việc tính đạo hàm trên toàn miền chứa điểm góc, điểm đứt gãy hoặc kỳ dị yếu.
  2. Kỹ thuật trung bình hóa số liệu đầu vào: Đề xuất quy trình sử dụng IAD để làm trơn các bộ số liệu ban đầu bị rời rạc, phi tuyến hoặc có bước nhảy gián đoạn (như sự phân bố tiết diện tương tác trong các môi trường vật liệu không đồng nhất), đưa về dạng dữ liệu liên tục phục vụ tính toán số.
  3. Phát triển thuật toán toán tử ma trận cho động học lò phản ứng: Xây dựng thành công kỹ thuật giải tích ma trận giải phương trình khuếch tán đa nhóm có neutron trễ bằng cách khai triển ma trận hàm mũ $\exp(hA)$ thông qua các trị riêng và phép khử Gauss, khắc phục hiện tượng mất ổn định nghiệm số trong quá trình lặp.

Đóng góp thực tiễn và giải pháp đề xuất

  1. Tính toán tối ưu hóa thông lượng neutron: Ứng dụng thành công phương pháp IAD vào tính toán thông lượng neutron và hệ số nhân hiệu dụng cho các cấu hình lò phản ứng hình học cơ bản (tấm phẳng, hình trụ, hình cầu). Thuật toán IAD giúp giảm đáng kể thời gian tính toán và dung lượng bộ nhớ máy tính so với phương pháp chia lưới Nodal biến phân và phương pháp ngẫu nhiên Monte Carlo.
  2. Cung cấp công cụ dự báo an toàn lò phản ứng: Phương pháp giải nhanh thông lượng cho phép người vận hành lò phản ứng hạt nhân dự đoán sớm sự thay đổi trường bức xạ và điều chỉnh hệ số nhân hiệu dụng $k_{eff}$ về ngưỡng an toàn khi xảy ra biến động công suất hoặc sự cố cục bộ.
  3. Giải pháp cho bài toán nhiễu xạ cơ học: Đưa ra phương pháp giải số đơn giản cho bài toán nhiễu xạ sóng trong môi trường đàn hồi có góc khuyết $\alpha$ thông qua việc kết hợp biến hình bảo giác và IAD, loại bỏ sự phụ thuộc vào các phần mềm chia lưới phần tử hữu hạn cồng kềnh.

Danh mục công bố: Các kết quả của luận án đã được công bố trong 08 công trình khoa học gồm các bài báo trên tạp chí và kỷ yếu hội nghị chuyên ngành trong nước và quốc tế, trong đó có 05 bài báo đăng trên các tạp chí khoa học quốc tế.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của luận án

Căn cứ vào phạm vi và nội dung được trình bày trong văn bản, luận án tồn tại một số hạn chế sau:

  • Các bài toán lò phản ứng hạt nhân chủ yếu mới dừng lại ở việc khảo sát trạng thái dừng (bài toán tĩnh, $\partial \Phi / \partial t = 0$), chưa xét đến các quá trình chuyển tiếp động học phi tuyến phức tạp trong không gian ba chiều khi các thanh điều khiển dịch chuyển liên tục theo thời gian thực.
  • Mô hình hình học được giải số chi tiết trong luận án vẫn là các mô hình đơn giản hóa (mô hình một chiều, hình cầu, tấm phẳng hoặc hình trụ đồng nhất); chưa áp dụng trực tiếp cho toàn bộ cấu trúc không gian ba chiều thực tế của các bó nhiên liệu với hàng nghìn thanh nhiên liệu không đồng nhất.
  • Phép gần đúng khuếch tán vẫn dựa trên giả thiết tán xạ đẳng hướng; do đó tại các vùng biên giáp chân không hoặc vùng vật liệu hấp thụ rất mạnh (nơi quãng chạy tự do của neutron bị suy giảm đột ngột), sai số của lý thuyết khuếch tán vẫn lớn hơn so với phương trình vận chuyển Boltzmann chính xác.
  • Bài toán nhiễu xạ sóng cơ học mới khảo sát trên môi trường đàn hồi thuần nhất, đẳng hướng với một góc khuyết hình nêm duy nhất, chưa mở rộng cho môi trường dị hướng hoặc có hệ thống nhiều vết nứt tương tác.

Hướng nghiên cứu tiếp

  • Mở rộng phương pháp IAD để giải trực tiếp phương trình vận chuyển Boltzmann đầy đủ (phương pháp $S_N$ hoặc $P_N$ kết hợp IAD) trong không gian ba chiều nhằm xử lý chính xác các điều kiện biên chân không và vùng hấp thụ mạnh.
  • Phát triển các chương trình tính toán động học lò phản ứng phụ thuộc thời gian khi có phản hồi nhiệt độ và phản hồi thủy nhiệt.
  • Ứng dụng phương pháp IAD giải các bài toán cơ học rạn nứt, cơ học phá hủy với hệ thống vết nứt ngẫu nhiên trong môi trường đàn hồi phi đồng nhất.

Giá trị tham khảo

Bản luận án có giá trị tham khảo chuyên môn đối với các nhóm đối tượng cụ thể sau:

  • Nghiên cứu sinh và học viên cao học ngành Vật lý lý thuyết, Vật lý toán và Toán ứng dụng: Có thể tham khảo cơ sở lý thuyết toán học của đạo hàm trung bình tích phân (Chương 1), kỹ thuật xử lý các toán tử vi phân có hệ số không trơn, và phương pháp giải phương trình vi tích phân kỳ dị bằng cách phối hợp giải tích phức và giải tích số.
  • Giảng viên và nhà nghiên cứu trong lĩnh vực Vật lý hạt nhân và Kỹ thuật lò phản ứng: Sử dụng tài liệu này (đặc biệt là Chương 2) như tài liệu chuyên khảo về việc áp dụng lý thuyết độ tuổi Fermi, phương pháp đa nhóm và phương pháp toán tử ma trận để mô phỏng sự phân bố thông lượng neutron và tính toán kích thước tới hạn cho các dạng hình học lò khác nhau.
  • Kỹ sư cơ học và người làm chính sách an toàn hạt nhân: Tham khảo quy trình phân tích và đánh giá hệ số nhân hiệu dụng $k_{eff}$ để thiết kế hệ thống điều khiển lò phản ứng; tham khảo thuật toán xử lý sóng đàn hồi trong môi trường có khuyết tật (Chương 3) phục vụ kiểm tra không phá hủy vật liệu.

Câu hỏi thường gặp

1. Phương pháp IAD định nghĩa và xử lý đạo hàm tại các điểm không trơn hoặc gián đoạn như thế nào?

Theo định nghĩa trong Chương 1, phương pháp IAD sử dụng biến phân trung bình tích phân: $$\delta F(x_0, h) = \lim_{t \to 0} \frac{1}{2t}\int_{-t}^t [F(x_0 + \tau h) - F(x_0)]d\tau$$ Nếu hàm có đạo hàm hai phía tại điểm $x_0$ (như điểm góc hoặc điểm nhảy bậc gián đoạn loại 1), đạo hàm trung bình tích phân tồn tại và bằng trung bình cộng đạo hàm hai phía: $$F'(x_0) = \frac{F'(x_0 + 0) + F'(x_0 - 0)}{2}$$ Ví dụ, hàm $F(x) = |x|$ không có đạo hàm thông thường tại $x_0 = 0$, nhưng có đạo hàm IAD bằng $0$. Điều này giúp thuật toán xử lý trơn tru các điểm đứt gãy mà không làm gián đoạn tính toán số.

2. Công thức xác định bán kính tới hạn ($R_{th}$) cho lò phản ứng hạt nhân hình cầu được thiết lập ra sao?

Trong Chương 2, đối với lò phản ứng hình cầu không có chất phản xạ, phương trình sóng cho neutron nhiệt có dạng $\Delta \Phi + \omega^2 \Phi = 0$ với nghiệm riêng thỏa mãn điều kiện không âm ở tâm là $\Phi(r) = C_1 \frac{\sin(Br)}{r}$. Từ điều kiện biên trên biên ngoại suy $\Phi(R) = 0$, ta thu được giá trị riêng nhỏ nhất là $B = \pi / R$. Tại trạng thái tới hạn, thông số hình học bằng thông số vật liệu ($B = \omega$), suy ra bán kính tới hạn: $$R_{th} = \frac{\pi}{\omega} = \pi \sqrt{\frac{L^2}{k_\infty - 1}}$$

3. Trong lý thuyết độ tuổi Fermi, hệ số nhân hiệu dụng ($k_{eff}$) được biểu diễn qua các thành phần xác suất nào?

Theo biểu thức (2.52) và (2.53) trong Chương 2, hệ số nhân hiệu dụng được viết dưới dạng: $$k_{eff} = k_\infty \cdot P = k_\infty \cdot \frac{\exp(-B^2 \tau)}{1 + B^2 L^2}$$ Trong đó, xác suất tránh rò rỉ toàn phần $P$ gồm hai thành phần tích hợp: tử số $\exp(-B^2 \tau)$ là xác suất để neutron nhanh không bị rò rỉ ra ngoài trong quá trình làm chậm đến năng lượng nhiệt (với độ tuổi neutron $\tau$), và mẫu số $1 + B^2 L^2$ là xác suất để neutron nhiệt không bị rò rỉ trong quá trình khuếch tán trước khi bị hấp thụ (với $L^2$ là bình phương độ dài khuếch tán).

4. Kỹ thuật toán tử ma trận trong luận án giải quyết hệ phương trình động học đa nhóm như thế nào?

Sau khi rời rạc hóa biến không gian bằng sai phân và trung bình hóa theo IAD, hệ phương trình được đưa về dạng ma trận cấp $K = G + I - 1$: $$\frac{d\Phi(t)}{dt} = A\Phi(t) + B$$ Thay vì giải lặp số từng bước thời gian dễ gây tích tụ sai số, tác giả sử dụng nghiệm giải tích ma trận $\Phi(t_{m+1}) = \exp(hA)\Phi(t_m) + (\exp(hA) - I)A^{-1}B$. Hàm mũ ma trận $\exp(hA)$ được khai triển thành đa thức $\sum_{k=0}^K \alpha_k (hA)^k$, trong đó các trọng số $\alpha_k$ được tìm chính xác qua phép khử Gauss dựa trên các trị riêng $\lambda_k$ của ma trận $A$.

5. Bài toán nhiễu xạ sóng đàn hồi có góc hổng $\alpha$ ở Chương 3 được giải quyết bằng công cụ toán học nào?

Để giải quyết sự kỳ dị tại đỉnh góc hình nêm $\alpha$, tác giả không dùng phương pháp sai phân trực tiếp mà áp dụng phép biến hình bảo giác của giải tích hàm biến phức. Phép biến đổi này đưa miền biên gấp khúc về miền biên chuẩn, chuyển bài toán vi phân sóng đàn hồi thành hệ phương trình vi tích phân kỳ dị. Sau đó, phương pháp IAD được áp dụng để khử kỳ dị yếu và tìm nghiệm xấp xỉ bằng phương pháp số, cho kết quả hàm chuyển vị $u(x)$ trùng khớp với nghiệm kiểm chứng trong tài liệu [95].

Kết luận

Luận án tiến sĩ của tác giả Phan Huy Thiện đã phát triển và ứng dụng thành công phương pháp đạo hàm trung bình tích phân (IAD) kết hợp với các kỹ thuật ma trận, sai phân hữu hạn và biến đổi bảo giác để giải quyết lớp bài toán vật lý và cơ học có miền biên không trơn hoặc chứa điểm kỳ dị. Trong lĩnh vực vật lý hạt nhân, công trình đã xây dựng thuật toán ma trận giải phương trình khuếch tán neutron đa nhóm có neutron trễ, cho phép tính toán chính xác thông lượng neutron và hệ số tới hạn với chi phí thời gian và bộ nhớ tối ưu hơn so với các phương pháp Monte Carlo hay Nodal truyền thống. Trong cơ học đàn hồi, phương pháp IAD phối hợp với biến hình bảo giác đã giải quyết hiệu quả bài toán nhiễu xạ sóng trên miền biên có góc khuyết hình nêm. Toàn bộ các kết quả lý thuyết và mô phỏng số được bảo chứng thông qua 08 công trình khoa học đã công bố, đóng góp một hướng tiếp cận toán học tính toán hữu ích cho các bài toán biên phức tạp trong vật lý lý thuyết và cơ học ứng dụng.