Tổng quan về luận án

Nghiên cứu bài toán ngược (Inverse Problems) và bài toán không chỉnh (Ill-posed Problems) cho các phương trình đạo hàm riêng phân số (Fractional Partial Differential Equations - FPDEs) là một trong những địa hạt tiên phong của toán học ứng dụng hiện đại. Luận án tiến sĩ chuyên ngành Toán giải tích (Mã số: 9460102) của nghiên cứu sinh Lê Đình Long với đề tài "Xác định điều kiện đầu và hàm nguồn cho một số phương trình khuếch tán" dưới sự hướng dẫn khoa học của TS. Bùi Lê Trọng Thanh và TS. Nguyễn Anh Triết tại Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia Thành phố Hồ Chí Minh, đã giải quyết căn cơ các bài toán nhận dạng hàm nguồn và khôi phục trạng thái ban đầu trong các quá trình truyền nhiệt và khuếch tán dị thường (anomalous diffusion).

Trong thực tiễn vật lý, cơ học môi trường không đồng nhất, thủy văn nước ngầm và mô hình sinh lý học tế bào, các quá trình khuếch tán thường không tuân theo quy luật cổ điển Fickian mà thể hiện tính chất khuếch tán chậm (sub-diffusion) với hiện tượng ghi nhớ lịch sử dài hạn (memory effect) và liên kết không gian phi địa phương (non-local spatial interactions). Các hiện tượng này được mô hình hóa chính xác thông qua đạo hàm phân số Riemann-Liouville và Caputo. Tuy nhiên, việc xác định hàm nguồn $f(x)$ hoặc khôi phục trạng thái ban đầu $u(x,0)$ từ dữ liệu đo đạc tại thời điểm cuối $u(x,T) = h(x)$ là bài toán không chỉnh nghiêm trọng theo nghĩa Hadamard: nghiệm không tồn tại, không duy nhất hoặc không phụ thuộc liên tục vào dữ liệu đo đạc. Một sai số nhiễu cực nhỏ $\delta > 0$ ở dữ liệu biên-cuối $|h^\delta - h|_{L^2(\Omega)} \le \delta$ sẽ bị khuếch đại bùng nổ theo cấp số nhân trong quá trình giải ngược, làm sai lệch hoàn toàn nghiệm số.

                  ┌────────────────────────────────────────────────────────┐
                  │ Dữ liệu quan sát tại t = T: u(x, T) = h(x)             │
                  │ Nhiễu thực nghiệm: ||h_delta - h||_{L^2} <= delta      │
                  └──────────────────────────┬─────────────────────────────┘
                                             │
                                             ▼
                  ┌────────────────────────────────────────────────────────┐
                  │ BÀI TOÁN NGƯỢC KHÔNG CHỈNH THEO NGHĨA HADAMARD         │
                  │   Toán tử tích phân ngược suy biến phổ:                │
                  │   E_{alpha,alpha}(-lambda_n t^alpha) -> 0 khi n -> oo  │
                  └──────────────────────────┬─────────────────────────────┘
                                             │
               ┌─────────────────────────────┼─────────────────────────────┐
               ▼                             ▼                             ▼
┌─────────────────────────────┐┌───────────────────────────┐┌─────────────────────────────┐
│  Phương pháp Tựa Giá trị    ││   Phương pháp Chặt cụt    ││     Phương pháp Lặp         │
│         Biên (QBV)          ││      Chuỗi Fourier        ││   Fractional Landweber      │
│  alpha_reg*u(x,0)+u(x,T)=h  ││    P_N u = sum_{n=1}^N    ││  f^{k+1} = f^k - tau K*(...)│
└──────────────┬──────────────┘└─────────────┬─────────────┘└──────────────┬──────────────┘
               │                             │                             │
               └─────────────────────────────┼─────────────────────────────┘
                                             ▼
                  ┌────────────────────────────────────────────────────────┐
                  │ ĐÁNH GIÁ SAI SỐ HỘI TỤ TỐI ƯU TRÊN L^2(Omega) / H^s   │
                  │ - Quy tắc chọn tham số Tiên nghiệm (A-priori Rule)     │
                  │ - Quy tắc chọn tham số Hậu nghiệm (Morozov Discrepancy)│
                  │   Tốc độ hội tụ: ||u_reg - u_exact|| <= C * delta^p    │
                  └────────────────────────────────────────────────────────┘

Khoảng trống nghiên cứu (Research Gap) then chốt mà luận án trực tiếp giải quyết bao gồm:

  1. Thiếu hụt các khung phân tích định lượng về sự không chỉnh của bài toán khôi phục hàm nguồn phụ thuộc không gian đối với phương trình khuếch tán chứa đạo hàm Riemann-Liouville bậc phân số kết hợp toán tử elliptic suy biến tổng quát.
  2. Các công trình trước đây của Wei và Wang (2014), Liu và Yamamoto (2011) chủ yếu tập trung vào đạo hàm Caputo hoặc phương trình cấp nguyên truyền thống, chưa mở rộng đầy đủ sang toán tử Fractional Laplacian $(-\Delta)^s$ với $s \in (0,1)$ kết hợp điều kiện tích phân phi địa phương theo thời gian.
  3. Hạn chế trong việc thiết lập đánh giá sai số chỉnh hóa sắc bén (sharp error estimates) khi dữ liệu quan sát nằm ngoài không gian Hilbert $L^2(\Omega)$, cụ thể là trường hợp không gian Banach $L^q(\Omega)$ ($q \neq 2$) đối với phương trình Biharmonic phi tuyến bậc bốn.

Các câu hỏi và giả thuyết nghiên cứu được xác định rõ ràng:

  • RQ1: Cấu trúc toán tử phổ của bài toán khôi phục hàm nguồn bậc phân số suy biến như thế nào khi các giá trị riêng $\lambda_n \to \infty$?
  • RQ2: Liệu phương pháp tựa giá trị biên sửa đổi (Quasi-Boundary Value Method) và phương pháp Fractional Landweber có đảm bảo tốc độ hội tụ bậc tối ưu Hölder-type dưới cả quy tắc chọn tham số tiên nghiệm ($a$-priori) lẫn hậu nghiệm ($a$-posteriori)?
  • RQ3: Làm thế nào để thiết lập sự ổn định nghiệm trong không gian $L^q(\Omega)$ cho bài toán giá trị cuối của phương trình Biharmonic dạng $\mathcal{M}u = u_{tttt} + 2u_{ttxx} + u_{xxxx} = G(x,t,u(x,t))$ với hàm nguồn phi tuyến thỏa điều kiện Lipschitz địa phương?

Khung lý thuyết của nghiên cứu được xây dựng vững chắc trên nền tảng Giải tích hàm hiện đại, Lý thuyết toán tử tự liên hợp nửa dưới trong không gian Hilbert scale $H^s(\Omega)$, Giải tích phân số (Fractional Calculus) với các hàm đặc biệt Gamma $\Gamma(z)$, Mittag-Leffler hai tham số $E_{\alpha,\beta}(z)$, và Lý thuyết chỉnh hóa bài toán ngược Hadamard-Tikhonov. Nghiên cứu mang lại đóng góp đột phá với 5 công trình công bố trên các tạp chí quốc tế hàng đầu thuộc danh mục ISI/Web of Science (Mathematical Methods in the Applied Sciences, Q1, IF = 3.1; Computational and Applied Mathematics, Q2, IF = 3.0; Journal of Inequalities and Applications, Q2, IF = 4.0; Acta Mathematica Sinica, Q2, IF = 0.9; AIMS Mathematics, Q2, IF = 2.1).

Literature Review và Positioning

Lịch sử nghiên cứu các bài toán parabolic ngược thời gian khởi nguồn từ công trình kinh điển của John (1955) khi giải bài toán Cauchy cho phương trình truyền nhiệt ngược bằng phương pháp số và chứng minh tính ổn định có điều kiện trên tập các hàm dương bị chặn. Bước ngoặt phương pháp luận xuất hiện khi Tikhonov (1963) công bố nguyên lý chỉnh hóa biến phân (Tikhonov Regularization), biến các bài toán toán tử không chỉnh thành bài toán cực tiểu hóa phiếm hàm khả vi chặt chẽ. Tiếp đó, Lattes và Lions (1967) phát triển phương pháp tựa đảo (Quasi-reversibility Method), Ewing (1977) ứng dụng không gian Sobolev, và Showalter (1974) hoàn thiện phương pháp tựa giá trị biên (Quasi-Boundary Value - QBV) cho phương trình tiến hóa thuần nhất.

┌────────────────────────────────────────────────────────────────────────────────────────┐
│                   TIẾN TRÌNH PHÁT TRIỂN LÝ THUYẾT BÀI TOÁN NGƯỢC                       │
├────────────────────────────────────────────────────────────────────────────────────────┤
│ 1955: F. John         -> Bài toán truyền nhiệt ngược thời gian (Cauchy problem)        │
│ 1963: A. N. Tikhonov  -> Lý thuyết chỉnh hóa biến phân kinh điển                       │
│ 1967: Lattes & Lions  -> Phương pháp Tựa đảo (Quasi-reversibility method)              │
│ 1974: R. E. Showalter -> Phương pháp Tựa giá trị biên (Quasi-boundary value method)   │
│ 2011: Liu & Yamamoto  -> Bài toán ngược thời gian cho FPDEs cấp phân số 1D             │
│ 2014: Wei & Wang      -> Chỉnh hóa xác định hàm nguồn phân số qua tích phân loại 1     │
│ 2016: Tuấn, Long et al-> Khôi phục điều kiện đầu qua QBV cải tiến                      │
│ 2023: Lê Đình Long    -> Đột phá đồng thời: Riemann-Liouville, Fractional Laplacian    │
│      (Luận án này)       $(-\Delta)^s$, Fractional Landweber, QBV trên $L^q(\Omega)$   │
└────────────────────────────────────────────────────────────────────────────────────────┘

Trong dòng chảy nghiên cứu quốc tế về phương trình vi phân phân số, hai trường phái chính đã hình thành:

  1. Trường phái tiếp cận Biến phân và Tikhonov: Đại diện bởi Ting Wei và Jun Xian (2019), sử dụng phương pháp biến phân yếu kết hợp thuật toán Conjugate Gradient để tối thiểu hóa sai số số học.
  2. Trường phái tiếp cận Giải tích phổ và Toán tử chỉnh hóa trực tiếp: Đại diện bởi Ji-Jun Liu và Masahiro Yamamoto (2011), tập trung khai triển chuỗi nghiệm Fourier thông qua biểu diễn giải tích của toán tử giải (resolvent operator) và hạt nhân Mittag-Leffler.

Các tranh luận học thuật sâu sắc xoay quanh tính tối ưu của tốc độ hội tụ khi chọn tham số chỉnh hóa. Trong khi quy tắc tiên nghiệm yêu cầu thông tin trơn $a$-priori $|u|_{H^p(\Omega)} \le E$ (vốn khó khả thi trong đo đạc thực tế), quy tắc hậu nghiệm dựa trên nguyên lý sai số Morozov (Morozov’s Discrepancy Principle) thường gặp thách thức lớn về sự hội tụ phi tuyến khi phương trình chứa toán tử phi địa phương.

Luận án của Lê Đình Long định vị chính xác tại giao điểm của các tranh luận này. Bằng cách kết hợp phép biến đổi tích phân Laplace-Fourier, định lý nhúng Sobolev, và tính chất tiệm cận của hàm Mittag-Leffler: $$E_{\alpha,\beta}(-z) = -\sum_{k=1}^m \frac{(-z)^{-k}}{\Gamma(\beta - \alpha k)} + \mathcal{O}(|z|^{-m-1}), \quad |z| \to \infty$$ tác giả đã nâng tầm các kết quả của Ting Wei và Jun-Gang Wang (2014) từ phương trình khuếch tán phân số một chiều với toán tử Laplace chuẩn sang phương trình khuếch tán đa chiều với toán tử elliptic suy biến tổng quát $A = -\sum_{i,j=1}^d \partial_{x_i}(a_{ij}(x)\partial_{x_j}) + c(x)$, đồng thời thiết lập các chặn trên sai số hội tụ tối ưu mà các nghiên cứu trước chưa đạt được.

Tiêu chí So sánh Nghiên cứu của Liu & Yamamoto (2011) Công trình của Wei & Wang (2014) Luận án Lê Đình Long (2023)
Loại đạo hàm thời gian Caputo fractional $\partial_t^\alpha$ ($0 < \alpha < 1$) Caputo fractional $\partial_t^\alpha$ Riemann-Liouville $\partial_t^{1-\alpha}$ và Caputo $\partial_t^\alpha$
Toán tử không gian 1D Classical Laplacian $-\frac{d^2}{dx^2}$ Bounded domain 2D Laplacian $-\Delta$ Elliptic tổng quát $\mathcal{A}$, Fractional Laplacian $(-\Delta)^s$
Không gian hàm quan sát Không gian Hilbert $L^2(0,1)$ Không gian Hilbert $L^2(\Omega)$ Mở rộng trên cả $L^2(\Omega)$ và không gian Banach $L^q(\Omega)$ ($q \neq 2$)
Phương pháp chỉnh hóa Quasi-reversibility Quasi-Boundary Value cổ điển QBV cải tiến, Fourier Truncation, Fractional Landweber
Quy tắc chọn tham số Tiên nghiệm ($a$-priori) Tiên nghiệm và Hậu nghiệm Chứng minh đầy đủ $a$-priori & $a$-posteriori với chặn sai số sắc bé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 sâu sắc lý thuyết bài toán ngược cho phương trình vi phân trừu tượng trong không gian Banach và Hilbert scale. Nghiên cứu đã chứng minh toán học tính không chỉnh nghiệm của bài toán xác định hàm nguồn $f(x)$ thông qua việc chỉ ra toán tử giải tích phân loại một có phổ kỳ dị tích lũy tại vô hạn.

Cụ thể, đối với bài toán khôi phục hàm nguồn cho hệ: $$\partial_t u(x,t) = \partial_t^{1-\alpha}\Delta u(x,t) + f(x)q(t), \quad (x,t) \in \Omega \times (0,T)$$ nghiệm dưới dạng chuỗi Fourier được biểu diễn qua tích phân hàm Mittag-Leffler: $$u_n(T) = u_n(0)E_{\alpha,1}(-\lambda_n T^\alpha) + f_n \int_0^T (T-\tau)^{\alpha-1} E_{\alpha,\alpha}(-\lambda_n(T-\tau)^\alpha) q(\tau) d\tau$$ Khi $n \to \infty$, giá trị riêng $\lambda_n \to \infty$, thừa số tích phân phân rã về 0 với tốc độ tiệm cận $\mathcal{O}(\lambda_n^{-1})$, khiến cho hệ số Fourier của hàm nguồn $f_n \approx \lambda_n h_n$ bị khuếch đại vô hạn khi dữ liệu $h(x)$ bị nhiễu.

Tác giả đã đề xuất và chứng minh thành công 4 định lý hội tụ nền tảng, thiết lập các bất đẳng thức năng lượng mới trong không gian Sobolev $H^\gamma(\Omega)$:

  • Mệnh đề 1 (Ước lượng chặn dưới toán tử giải): Tồn tại hằng số $M_1(\alpha, T) > 0$ sao cho: $$\int_0^T (T-\tau)^{\alpha-1} E_{\alpha,\alpha}(-\lambda_n(T-\tau)^\alpha) d\tau \ge \frac{M_1}{1 + \lambda_n T^\alpha}$$
  • Mệnh đề 2 (Đánh giá sai số tiên nghiệm QBV): Khi chọn tham số chỉnh hóa $\epsilon(\delta) = \delta^{\frac{2}{p+2}}$, nghiệm chỉnh hóa $f_\epsilon^\delta$ thỏa mãn: $$|f_\epsilon^\delta - f|_{L^2(\Omega)} \le C_1 \delta^{\frac{p}{p+2}}$$ với $p > 0$ là chỉ số độ trơn của hàm nguồn trong không gian $H^p(\Omega)$.
  • Mệnh đề 3 (Đánh giá sai số hậu nghiệm theo Morozov): Thiết lập sự hội tụ mạnh khi tham số $\epsilon = \epsilon(\delta, h^\delta)$ được xác định qua phương trình phi tuyến $|K f_\epsilon^\delta - h^\delta|_{L^2} = \tau \delta$ ($\tau > 1$).
  • Mệnh đề 4 (Ổn định nghiệm phương trình Biharmonic phi tuyến): Chứng minh tính duy nhất và liên tục của nghiệm chỉnh hóa trong không gian $C([0,T]; L^q(\Omega))$ dưới điều kiện hàm nguồn $G(x,t,u)$ thỏa mãn bất đẳng thức tăng trưởng cục bộ.
                    ┌────────────────────────────────────────────────────────┐
                    │               KHUNG PHÂN TÍCH TỔNG QUÁT                │
                    └──────────────────────────┬─────────────────────────────┘
                                               │
               ┌───────────────────────────────┴───────────────────────────────┐
               ▼                                                               ▼
┌──────────────────────────────┐                              ┌──────────────────────────────┐
│  Toán tử Phân số Thời gian   │                              │   Toán tử Không gian         │
│  - Riemann-Liouville         │                              │   - Elliptic suy biến A      │
│    \partial_t^{1-\alpha}     │                              │   - Fractional Laplacian     │
│  - Caputo \partial_t^\alpha  │                              │     (-\Delta)^s              │
└──────────────┬───────────────┘                              └──────────────┬───────────────┘
               │                                                             │
               └───────────────────────────────┬─────────────────────────────┘
                                               │
                                               ▼
                    ┌────────────────────────────────────────────────────────┐
                    │   Toán tử Giải tích phân Tích chập Mittag-Leffler       │
                    │   K f = \int_0^T (T-\tau)^{\alpha-1}                   │
                    │         * E_{\alpha,\alpha}(-A(T-\tau)^\alpha) f d\tau │
                    └──────────────────────────┬─────────────────────────────┘
                                               │
                                               ▼
                    ┌────────────────────────────────────────────────────────┐
                    │   Không gian Hàm và Phép Nhúng Sobolev:                │
                    │   L^2(\Omega) \subset H_0^1(\Omega) \subset H^p(\Omega)│
                    │   Mở rộng Banach: L^q(\Omega) (q \neq 2)               │
                    └────────────────────────────────────────────────────────┘

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

Khung phân tích của luận án tích hợp nhuần nhuyễn 3 trụ cột lý thuyết: Lý thuyết phổ toán tử vi phân elliptic đối xứng xác định dương, Lý thuyết giải tích phân số với các hàm chuyển tiếp Mittag-Leffler, và Lý thuyết xấp xỉ chỉnh hóa không gian Hilbert. Cách tiếp cận mới lạ nằm ở việc biến đổi bài toán biên phi địa phương phức tạp thành các phương trình toán tử tích phân loại một trong không gian mở rộng, sau đó sử dụng toán tử tựa biên sửa đổi chứa tham số suy giảm lũy thừa để dập tắt sự bất ổn định ở các tần số cao trong chuỗi Fourier.

Điều kiện biên và giới hạn ứng dụng (Boundary Conditions) được xác định tường minh: miền không gian $\Omega \subset \mathbb{R}^d$ ($d \in {1,2,3}$) là miền bị chặn với biên trơn $\partial\Omega \in C^2$, toán tử elliptic thỏa mãn điều kiện elliptic đều $\sum_{i,j=1}^d a_{ij}(x)\xi_i\xi_j \ge \mu |\xi|^2$ ($\mu > 0$), và hàm trọng số thời gian $q(t) \in C([0,T])$ thỏa mãn $q(t) \ge q_0 > 0$ hầu khắp nơi trên $[0,T]$.

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

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan khoa học thực chứng suy diễn (Positivist / Deductive Paradigm) của toán học giải tích chính xác. Thiết kế nghiên cứu là sự kết hợp hữu cơ giữa phân tích định lý giải tích thuần túy (rigorous mathematical proofs) và tính toán số mô phỏng (computational numerical experiments) nhằm kiểm chứng và xác nhận các chặn sai số lý thuyết.

Thiết kế đa cấp độ (Multi-level design) được cấu trúc chặt chẽ qua 4 bước:

  1. Khảo sát dạng nghiệm nhẹ (Mild solution) thông qua tích phân Fourier-Laplace.
  2. Thiết lập định lý không chỉnh Hadamard với bằng chứng giải tích về sự bùng nổ nghiệm.
  3. Xây dựng toán tử chỉnh hóa $\mathcal{R}_\epsilon$ tương ứng với từng phương pháp (QBV, Fourier Truncation, Fractional Landweber).
  4. Phân tích sự hội tụ và tốc độ hội tụ $\lim_{\delta \to 0} |\mathcal{R}{\epsilon(\delta)} h^\delta - f|{L^2} = 0$.
┌─────────────────────────────────────────────────────────────────────────────────────────┐
│                           QUY TRÌNH THIẾT KẾ NGHIÊN CỨU                                 │
├─────────────────────────────────────────────────────────────────────────────────────────┤
│ [GIAI ĐOẠN 1] THIẾT LẬP DẠNG NGHIỆM MILD                                                │
│  - Áp dụng biến đổi Laplace theo thời gian t và khai triển Fourier theo cơ sở {e_n(x)}   │
│  - Biểu diễn nghiệm qua hàm Mittag-Leffler hai tham số E_{\alpha, \beta}                │
├─────────────────────────────────────────────────────────────────────────────────────────┤
│ [GIAI ĐOẠN 2] CHỨNG MINH TÍNH KHÔNG CHỈNH THEO HADAMARD                                 │
│  - Đánh giá phổ toán tử giải: \lambda_n \to \infty => mẫu số phân rã về 0                │
│  - Trích xuất hằng số bùng nổ sai số: ||f^\delta - f|| -> \infty                        │
├─────────────────────────────────────────────────────────────────────────────────────────┤
│ [GIAI ĐOẠN 3] XÂY DỰNG NGHIỆM CHỈNH HÓA VÀ THUẬT TOÁN                                   │
│  - Thiết lập phương trình tựa biên / toán tử chiếu P_N / bước lặp Fractional Landweber  │
│  - Xác lập quy tắc chọn tham số tiên nghiệm \epsilon(\delta) và hậu nghiệm Morozov      │
├─────────────────────────────────────────────────────────────────────────────────────────┤
│ [GIAI ĐOẠN 4] MÔ PHỎNG SỐ HỌC VÀ KIỂM CHỨNG TRÊN MATLAB                                 │
│  - Rời rạc hóa lưới không gian - thời gian: N_x = 100, N_t = 100                        │
│  - Bơm nhiễu ngẫu nhiên Gauss: \delta \in {0.001, 0.015, 0.05}                          │
│  - Đo sai số tương đối và vẽ đồ thị hội tụ xác thực định lý lý thuyết                   │
└─────────────────────────────────────────────────────────────────────────────────────────┘

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu số mô phỏng được chuẩn hóa nghiêm ngặt:

  • Dữ liệu chính xác được sinh từ các nghiệm giải tích mẫu $u_{exact}(x,t)$ và hàm nguồn chính xác $f_{exact}(x)$.
  • Dữ liệu quan sát bị nhiễu được tạo lập bằng cách cộng nhiễu ngẫu nhiên phân phối đều/Gauss: $$h^\delta(x) = h(x) + \delta \cdot (2\cdot \text{rand}(x) - 1)$$ thỏa mãn điều kiện chuẩn $|h^\delta - h|_{L^2(\Omega)} \le \delta$.
  • Giao thức tam giác hóa (Triangulation) được thực hiện thông qua việc so sánh đối chéo kết quả giữa 3 phương pháp chỉnh hóa độc lập: Phương pháp Chặt cụt Fourier, Phương pháp Tựa giá trị biên, và Phương pháp Lặp Fractional Landweber.
  • Độ tin cậy và tính hợp lệ (Validity & Reliability) được khẳng định qua việc các cận sai số giải tích chặn trên luôn bảo toàn tính đúng đắn trên toàn bộ miền biến thiên của bậc đạo hàm $\alpha \in (0,1)$.

Data và phân tích

Các thử nghiệm số được thực thi trên môi trường tính toán khoa học hiệu năng cao sử dụng phần mềm MATLAB R2022b kết hợp gói Symbolic Math Toolbox và Maple.

  • Không gian tính toán 1D: $\Omega = (0, \pi)$, chia lưới không gian $N_x = 100$ điểm nút, bước lưới $\Delta x = \pi / 100$.
  • Không gian tính toán 2D: $\Omega = (0,\pi) \times (0,\pi)$, kích thước lưới $N_x \times N_y = 100 \times 100$, tổng số điểm tính toán $10,000$ nút.
  • Bước thời gian: $T = 1.0$, phân đoạn thời gian $N_t = 100$, $\Delta t = 0.01$.
  • Các mức nhiễu thực nghiệm được khảo sát chi tiết: $\delta \in {0.001, 0.005, 0.015, 0.03, 0.05}$.
  • Phân tích độ nhạy (Robustness checks) được thực hiện đa dạng qua các giá trị bậc phân số: $\alpha = 0.1, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9$.
  • Thuật toán số giải hàm Mittag-Leffler sử dụng thuật toán tính toán tích phân đường viền ma trận cải tiến của Podlubny và Garrappa, đảm bảo sai số máy dưới $10^{-14}$.

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

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

  1. Minh chứng toán học về cơ chế bùng nổ sai số: Luận án đã trích xuất tường minh công thức kỳ dị phổ, chứng minh rằng khi không có chỉnh hóa, tỉ số giữa sai số nghiệm và sai số dữ liệu tăng trưởng theo bậc hàm mũ/đa thức bậc cao của giá trị riêng: $$\frac{|f^\delta - f|_{L^2}}{\delta} \ge C \cdot \lambda_N \to \infty \quad \text{khi } N \to \infty$$
  2. Khám phá tính phụ thuộc phi tuyến của tốc độ hội tụ vào bậc phân số $\alpha$: Kết quả thực nghiệm số chỉ ra một phát hiện trực quan sâu sắc: khi bậc đạo hàm $\alpha \to 1^-$ (tiệm cận khuếch tán cổ điển), sai số chỉnh hóa giảm dần và nghiệm có độ trơn cao hơn; ngược lại, khi $\alpha \to 0^+$ (khuếch tán cực chậm với hiệu ứng nhớ vĩ mô), tính không chỉnh gia tăng mạnh, đòi hỏi tham số chỉnh hóa $\epsilon$ phải được làm mịn với trọng số lớn hơn.
  3. Sự vượt trội của Phương pháp Tựa giá trị biên cải tiến: Tại mức nhiễu $\delta = 0.015$, bảng sai số số học trong luận án chứng minh nghiệm chỉnh hóa $f_\epsilon^\delta(x)$ bám sát nghiệm chính xác $f(x)$ với sai số tương đối trong không gian $L^2$ đạt dưới $2.35%$ trên toàn bộ miền $\Omega = (0,\pi)$, vượt trội hơn so với phương pháp Tikhonov bậc 0 truyền thống (sai số $> 6.8%$).
  4. Xác lập sự hội tụ đồng đều trên không gian $L^q(\Omega)$ cho phương trình Biharmonic: Đối với bài toán Biharmonic bậc bốn phi tuyến $\mathcal{M}u = G(x,t,u)$, luận án chứng minh thành công tính ổn định nghiệm trên không gian Banach $L^q(\Omega)$ ($1 < q < \infty$), phá vỡ giới hạn truyền thống vốn chỉ áp dụng được trên không gian Hilbert $L^2(\Omega)$.
Sai số L^2
   │
0.08│                                  * Không chỉnh hóa (Nhiễu bùng nổ)
0.07│                                 /
0.06│                                /
0.05│       * Tikhonov cổ điển      /
0.04│      /
0.03│     /
0.02│    /  * Phương pháp QBV của luận án (Sai số <= 2.35% tại delta = 0.015)
0.01│   /  * Phương pháp Fractional Landweber
   0└──┴─────────┴─────────┴─────────┴─────────┴─────────┴─────────►
      0.1       0.3       0.5       0.7       0.9       Bậc đạo hàm alpha

Implications đa chiều

  • Về mặt lý thuyết: Mở rộng biên giới của giải tích vô hạn chiều, bổ sung các công cụ sắc bén để nghiên cứu các lớp phương trình đạo hàm riêng phi tuyến chứa toán tử không địa phương.
  • Về phương pháp luận: Cung cấp bộ khung thuật toán chỉnh hóa hoàn chỉnh từ xấp xỉ liên tục (QBV) đến xấp xỉ rời rạc (Fourier Truncation) và lặp đại số (Fractional Landweber), có thể chuyển giao trực tiếp cho các lớp bài toán ngược trong cơ học chất lưu và truyền sóng địa chấn.
  • Về ứng dụng thực tiễn: Tạo cơ sở thuật toán vững chắc cho việc xử lý ảnh y sinh (chụp cắt lớp điện toán ngực, tái tạo hình ảnh MRI/PET dựa trên mô hình khuếch tán phân số mô sinh học), định lượng nguồn ô nhiễm ngầm trong thủy văn học, và chuẩn đoán rạn nứt kết cấu công trình chịu tải trọng Biharmonic.
  • Khuyến nghị triển khai: Các viện nghiên cứu địa kỹ thuật và xử lý tín hiệu y tế có thể tích hợp trực tiếp thuật toán QBV và nguyên lý sai số Morozov vào các phần mềm mô phỏng hiện hành nhằm loại bỏ nhiễu đo đạc cảm biến mà không làm mất các đặc trưng kỳ dị biên.

Limitations và Future Research

Luận án thừa nhận một cách khách quan và khoa học các giới hạn nghiên cứu:

  1. Giả định nguồn độc lập thời gian: Các mô hình xác định hàm nguồn trong Chương 4 chủ yếu xét hàm nguồn có dạng tách biến $F(x,t) = f(x)q(t)$ với $q(t)$ đã biết, chưa giải quyết trọn vẹn bài toán hàm nguồn tổng quát phụ thuộc đồng thời cả không gian và thời gian $f(x,t)$.
  2. Bản chất của dữ liệu nhiễu: Luận án khảo sát mô hình nhiễu tất định (Deterministic noise) với chặn chuẩn $L^2$. Trong các hệ thống vật lý phức tạp, nhiễu đo đạc thường có bản chất ngẫu nhiên trắng (White noise / Gaussian random field).
  3. Độ trơn hình học của biên: Miền không gian $\Omega$ được giả định có biên trơn $C^2$ hoặc dạng hình chữ nhật $(0,\pi)^d$, chưa mở rộng cho các miền biên Lipschitz fractal phức tạp.

Chương trình nghiên cứu tương lai (Future Research Agenda) được vạch ra với 4 hướng phát triển cụ thể:

  • Hướng 1: Khảo sát bài toán khôi phục hàm nguồn phụ thuộc hoàn toàn vào biến thời gian $f(t)$ hoặc hàm nguồn phi tuyến $f(x, u(x,t))$.
  • Hướng 2: Mở rộng nghiên cứu bài toán giá trị cuối cho các mô hình khuếch tán phân số với các toán tử vi phân nhân không kỳ dị mới như đạo hàm Atangana-Baleanu-Caputo (ABC), đạo hàm Caputo-Fabrizio (CF), và đạo hàm Conformable.
  • Hướng 3: Khảo sát các mô hình phi tuyến mạnh với dữ liệu quan sát đo đạc trong không gian Orlicz và Sobolev trọng số.
  • Hướng 4: Phát triển lý thuyết chỉnh hóa thống kê (Statistical Inverse Problems) cho phương trình đạo hàm riêng phân số với dữ liệu nhiễu ngẫu nhiên Ito-Wiener.

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

Công trình nghiên cứu của NCS Lê Đình Long tạo ra những tác động học thuật và thực tiễn sâu rộng:

  • Tác động học thuật quốc tế: 5 công trình công bố trên các tạp chí chuyên ngành quốc tế uy tín (ISI-Q1/Q2) với hệ số ảnh hưởng cao (tổng IF tích lũy > 13.0) khẳng định năng lực hội nhập đỉnh cao của toán học giải tích Việt Nam trên trường quốc tế. Các bài báo của tác giả mở ra hướng trích dẫn tiềm năng lớn cho cộng đồng nghiên cứu bài toán ngược FPDEs toàn cầu.
  • Đóng góp cho đào tạo và nghiên cứu trong nước: Luận án đóng vai trò là tài liệu chuyên khảo mẫu mực, phục vụ trực tiếp cho công tác đào tạo sau đại học, các chuyên đề tiến sĩ, và định hướng nghiên cứu cho các nhóm Seminar chuyên sâu về Phương trình đạo hàm riêng và Giải tích ứng dụng tại Trường Đại học Khoa học Tự nhiên - ĐHQG TP.HCM.
  • Chuyển giao công nghệ tính toán: Thuật toán chỉnh hóa đề xuất cung cấp giải pháp lõi cho các kỹ sư R&D trong việc xây dựng phần mềm mô phỏng quá trình thẩm thấu chất ô nhiễm trong lòng đất, kiểm tra không phá hủy vật liệu composite (Non-destructive Testing), và tối ưu hóa phân bố nhiệt trong thiết bị bán dẫn siêu nhỏ.

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

┌────────────────────────────────────────────────────────────────────────────────────────┐
│                                 ĐỐI TƯỢNG HƯỞNG LỢI                                    │
├────────────────────────────────┬───────────────────────────────────────────────────────┤
│ Nghiên cứu sinh & Nhà toán học │ Tiếp cận các bổ đề giải tích sắc bén, kỹ thuật ước    │
│ (Doctoral Researchers)         │ lượng hàm Mittag-Leffler, và phương pháp mở rộng L^q.  │
├────────────────────────────────┼───────────────────────────────────────────────────────┤
│ Chuyên gia Nghiên cứu Cao cấp  │ Khung phân tích tổng quát cho toán tử Fractional      │
│ (Senior Academics)             │ Laplacian kết hợp tích phân phi địa phương.           │
├────────────────────────────────┼───────────────────────────────────────────────────────┤
│ Kỹ sư R&D & Khoa học Dữ liệu   │ Thuật toán số MATLAB ổn định cao, triệt tiêu nhiễu đo │
│ (Industry R&D)                 │ đạc trong xử lý tín hiệu y tế và địa vật lý.          │
├────────────────────────────────┼───────────────────────────────────────────────────────┤
│ Nhà hoạch định Chính sách MT   │ Cơ sở khoa học chính xác để mô phỏng và ban hành quy │
│ (Environmental Policymakers)   │ chuẩn kiểm soát ô nhiễm nguồn nước ngầm.               │
└────────────────────────────────┴───────────────────────────────────────────────────────┘

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 lý thuyết độc đáo nhất là việc thiết lập thành công khung chỉnh hóa tựa giá trị biên sửa đổi (Modified Quasi-Boundary Value Method) cho phương trình khuếch tán chứa đạo hàm phân số Riemann-Liouville bậc $1-\alpha$ kết hợp toán tử elliptic suy biến tổng quát. Công trình đã mở rộng trực tiếp lý thuyết chỉnh hóa kinh điển của Showalter (1974) và Tikhonov (1963) từ phương trình đạo hàm riêng cấp nguyên sang phương trình cấp phân số vô hạn chiều, chứng minh tốc độ hội tụ đạt bậc tối ưu kiểu Hölder $|f_\epsilon^\delta - f|_{L^2} \le C \delta^{\frac{p}{p+2}}$ dưới các điều kiện nguồn trơn tổng quát.

2. Đột phá phương pháp luận so với các công trình quốc tế tiền nhiệm được thể hiện ra sao?

So với công trình của Ji-Jun Liu và Masahiro Yamamoto (2011) (chỉ giới hạn ở phương trình 1D với đạo hàm Caputo) và công trình của Ting Wei và Jun-Gang Wang (2014) (áp dụng trên toán tử Laplace chuẩn), luận án của Lê Đình Long đã tạo ra bước nhảy vọt phương pháp luận:

  • Xử lý đồng thời cả đạo hàm Riemann-Liouville và Caputo trên miền đa chiều $\mathbb{R}^d$ ($d \ge 1$).
  • Tích hợp thành công toán tử Fractional Laplacian $(-\Delta)^s$ ($s \in (0,1)$) với điều kiện tích phân phi địa phương theo thời gian $\int_0^T u(x,t)dt = \phi(x)$.
  • Phát triển phương pháp Fractional Landweber với số bước lặp được kiểm soát chặt chẽ bởi nguyên lý sai số Morozov cải tiến.

3. Phát hiện thực nghiệm số nào gây bất ngờ nhất và bằng chứng dữ liệu đi kèm?

Phát hiện số học ấn tượng nhất là sự suy giảm phi tuyến tính của sai số chỉnh hóa theo bậc phân số $\alpha$: khi $\alpha$ tăng từ $0.1$ đến $0.9$, độ nhạy của nghiệm đối với nhiễu $\delta = 0.015$ giảm hơn $65%$, và sai số giữa nghiệm chỉnh hóa với nghiệm chính xác giảm từ $0.0542$ xuống $0.0189$ tại $N = 50$. Điều này chứng minh bằng chứng thực nghiệm rằng tính chất nhớ dài hạn của môi trường (ứng với $\alpha$ nhỏ) là nguyên nhân vật lý cốt lõi gia tăng mức độ không chỉnh của bài toán ngược.

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 tường minh toàn bộ giao thức tái lập:

  • Công thức giải tích của các hàm thử nghiệm mẫu và cấu trúc hàm riêng $e_n(x) = \sqrt{\frac{2}{\pi}}\sin(nx)$.
  • Thuật toán xấp xỉ chuỗi Fourier với ngưỡng chặt cụt $N \in {50, 100}$.
  • Quy tắc điều chỉnh tham số tiên nghiệm $\epsilon = \delta^\theta$ và thuật toán giải phương trình phi tuyến xác định tham số hậu nghiệm Morozov. Toàn bộ mã nguồn thử nghiệm được lập trình chi tiết trên MATLAB R2022b.

5. Lộ trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Lộ trình 10 năm được xác lập qua 3 giai đoạn chiến lược:

  • Giai đoạn 1 (2024-2026): Hoàn thiện lý thuyết cho hàm nguồn biến thiên toàn phần $f(x,t)$ và các đạo hàm phi địa phương nhân suy biến (ABC, CF).
  • Giai đoạn 2 (2027-2030): Phát triển hệ thống bài toán ngược ngẫu nhiên (Stochastic Inverse Problems) kết hợp giải tích đa thang (Multiscale Analysis) và học máy khoa học (Physics-Informed Neural Networks - PINNs) cho FPDEs.
  • Giai đoạn 3 (2031-2034): Ứng dụng công nghiệp hóa thuật toán vào các bộ công cụ phần mềm chuẩn đoán y sinh và mô phỏng địa kỹ thuật quy mô quốc gia.

Kết luận

Luận án tiến sĩ của NCS Lê Đình Long là một công trình khoa học xuất sắc, mẫu mực và toàn diện trong chuyên ngành Toán giải tích. Những đóng góp cốt lõi được đúc kết qua 6 điểm nhấn:

  1. Thiết lập biểu diễn nghiệm tích phân giải tích chuẩn xác thông qua hệ hàm Mittag-Leffler hai tham số cho phương trình khuếch tán bậc phân số với đạo hàm Riemann-Liouville và Caputo.
  2. Chứng minh toán học chặt chẽ tính không chỉnh theo nghĩa Hadamard của bài toán xác định hàm nguồn và bài toán giá trị cuối, chỉ ra cơ chế suy biến phổ gây bùng nổ sai số.
  3. Đề xuất và phát triển thành công 3 phương pháp chỉnh hóa hiệu quả cao: Phương pháp Tựa giá trị biên cải tiến (QBV), Phương pháp Chặt cụt Fourier, và Phương pháp Fractional Landweber.
  4. Chứng minh các định lý hội tụ sắc bén, xác lập tốc độ hội tụ tối ưu dưới cả hai chiến lược chọn tham số tiên nghiệm ($a$-priori) và hậu nghiệm ($a$-posteriori) theo nguyên lý sai số Morozov.
  5. Tiên phong mở rộng khảo sát bài toán giá trị cuối cho phương trình Biharmonic phi tuyến trong không gian Banach $L^q(\Omega)$ ($q \neq 2$).
  6. Cung cấp hệ thống thử nghiệm số học phong phú trên MATLAB, minh chứng tính ổn định, độ tin cậy và khả năng ứng dụng thực tiễn vượt trội của các phương pháp đề xuất.

Công trình không chỉ làm phong phú thêm kho tàng lý thuyết bài toán ngược hiện đại mà còn đặt nền móng vững chắc cho các ứng dụng tính toán khoa học liên ngành trong tương lai.