Tổng quan về luận án

Nghiên cứu về các bài toán ngược (inverse problems) và bài toán không chỉnh (ill-posed problems) cho hệ phương trình đạo hàm riêng là một trong những trụ cột trọng yếu của vật lý toán hiện đại, lý thuyết giải tích ứng dụng và khoa học tính toán. Luận án tiến sĩ toán học chuyên ngành Toán giải tích (mã số: 9460102) 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" do nghiên cứu sinh Lê Đình Long thực hiện tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia TP. Hồ Chí Minh, 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ập trung giải quyết trọn vẹn những thách thức giải tích phức tạp liên quan đến tính không chỉnh theo nghĩa Hadamard đối với các mô hình khuếch tán dị thường cấp không nguyên (fractional anomalous diffusion) và phương trình vi phân cấp cao (biharmonic equations).

Trong thực tế các quá trình vật lý, động lực học chất lỏng, địa chất thủy văn, sinh lý học tế bào và dẫn truyền nhiệt trong môi trường không đồng nhất, hiện tượng nhớt đàn hồi và khuếch tán dị thường không thể mô tả chuẩn xác bằng các toán tử vi phân nguyên cổ điển. Do đó, việc ứng dụng các toán tử tích phân - vi phân phân đoạn như đạo hàm Riemann-Liouville bậc $1-\alpha$ ($\alpha \in (0,1)$), đạo hàm Caputo bậc $\alpha$, và toán tử vi phân không nguyên phân đoạn Fractional Laplacian $(-\Delta)^s$ ($s \in (0,1)$) trở thành yêu cầu tất yếu. Tuy nhiên, các bài toán khôi phục hàm nguồn hoặc truy nguyên 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)$ bị nhiễu là những bài toán không chỉnh nghiêm trọng: nghiệm không thỏa mãn tính phụ thuộc liên tục vào dữ liệu đầu vào. Một nhiễu loạn ngẫu nhiên vô cùng bé ở dữ liệu quan sát $\delta > 0$ có thể dẫn đến sự bùng nổ vô hạn của sai số nghiệm giải tích.

Luận án đặt ra và giải quyết 5 bài toán nghiên cứu then chốt:

  1. Bài toán 1 (Hàm nguồn 1D/đa chiều với đạo hàm Riemann-Liouville): Xác định hàm nguồn phụ thuộc không gian $f(x)$ thỏa mãn 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)$$ với điều kiện biên thuần nhất $u(x,t) = 0$ trên $\partial\Omega \times (0,T]$ và dữ liệu cuối $u(x,T) = h(x)$, trong đó ký hiệu $\partial_t^{1-\alpha}$ là đạo hàm Riemann-Liouville cấp $1-\alpha \in (0,1)$ định nghĩa bởi: $$\partial_t^{1-\alpha}u(x,t) = \frac{1}{\Gamma(\alpha)} \frac{d}{dt} \int_0^t (t-\tau)^{\alpha-1} u(x,\tau)d\tau, \quad t > 0$$
  2. Bài toán 2 (Bài toán giá trị cuối cho phương trình khuếch tán phi tuyến): Tìm phân bố trạng thái ban đầu $u(x,0)$ và trường trạng thái $u(x,t)$ thỏa mãn: $$\partial_t u(x,t) = \partial_t^{1-\alpha}\Delta u(x,t) + F(x,t), \quad (x,t) \in \Omega \times [0,T)$$ với dữ liệu cuối $u(x,T) = f(x)$, xử lý điểm kỳ dị gây mất ổn định tại $t=0$.
  3. Bài toán 3 (Hàm nguồn 2 chiều trên miền chữ nhật): Khôi phục hàm nguồn hai biến không gian $f(x,y)$ từ hệ: $$\partial_t u(x,y,t) = \partial_t^{1-\alpha}(u_{xx}(x,y,t) + u_{yy}(x,y,t)) + \beta(t)f(x,y), \quad (x,y,t) \in \Omega \times (0,T)$$ với miền bị chặn $\Omega = (0,\pi)^2$ và dữ liệu quan sát cuối $u(x,y,T) = g(x,y)$.
  4. Bài toán 4 (Phương trình giả-parabolic với Fractional Laplacian và dữ liệu phi địa phương): Khôi phục hàm nguồn $f(x)$ cho phương trình khuếch tán cấp cao với số hạng nhớt đàn hồi: $$\partial_t^\alpha u + (-\Delta)^s u + \nu \partial_t^\beta (-\Delta)^\gamma u = \phi(t)f(x), \quad (x,t) \in \Omega \times (0,T]$$ trong đó đạo hàm Caputo $\partial_t^\alpha u(x,t) = \frac{1}{\Gamma(1-\alpha)} \int_0^t (t-\tau)^{-\alpha} u_\tau(x,\tau)d\tau$, và toán tử Fractional Laplacian $(-\Delta)^s$ được xác định theo giá trị chính Cauchy (Principal Value): $$(-\Delta)^s u(x,t) = C_{d,s} , \text{P.V.} \int_{\mathbb{R}^d} \frac{u(x,t) - u(z,t)}{|x-z|^{d+2s}} dz, \quad s \in (0,1), \quad C_{d,s} = \frac{4^s \Gamma(d/2+s)}{\pi^{d/2}|\Gamma(-s)|}$$
  5. Bài toán 5 (Bài toán giá trị cuối cho phương trình Biharmonic cấp bốn phi tuyến): Tìm hàm $u(x,t)$ thỏa mãn phương trình vi phân đạo hàm riêng cấp 4: $$Mu = u_{tttt} + 2u_{ttxx} + u_{xxxx} = G(x,t,u(x,t)), \quad (x,t) \in \Omega \times (0,T]$$ với điều kiện biên Navier $u(x,t) = \Delta u(x,t) = 0$ và dữ liệu quan sát cuối trong không gian Lebesgue $L^q(\Omega)$ với $q \neq 2$.

Giá trị khoa học vượt bậc của công trình được bảo chứng vững chắc khi toàn bộ các kết quả cốt lõi được trích xuất từ 5 bài báo khoa học chất lượng cao công bố trên các tạp chí quốc tế hàng đầu thuộc danh mục ISI/Scopus:

  • Mathematical Methods in the Applied Sciences (ISI Q1, IF = 3.0)
  • Computational and Applied Mathematics (ISI Q2, IF = 3.0)
  • Journal of Inequalities and Applications (ISI Q2, IF = 4.0)
  • Acta Mathematica Sinica, English Series (ISI Q2)
  • AIMS Mathematics (ISI Q2, IF = 2.1)

Literature Review và Positioning

Lịch sử nghiên cứu bài toán ngược đối với phương trình đạo hàm riêng parabolic bắt nguồn từ công trình tiên phong của Fritz John (1955) về tính ổn định số của bài toán Cauchy ngược thời gian cho phương trình truyền nhiệt. Bước ngoặt căn bản của lý thuyết chỉnh hóa bài toán không chỉnh được xác lập bởi viện sĩ A. N. Tikhonov (1963) với phương pháp hàm chính quy hóa Tikhonov. Tiếp sau đó, các trường phái toán giải tích lớn trên thế giới đã phát triển hàng loạt công cụ đại số toán tử mạnh mẽ: R. Lattes và J.-L. Lions (1967) đề xuất phương pháp tựa đảo (quasi-reversibility method); R. E. Ewing (1983) ứng dụng phương pháp phương trình Sobolev; và R. E. Showalter (1983) sáng tạo phương pháp giá trị tựa biên (Quasi-Boundary Value - QBV) cho các phương trình tiến hóa ngược thời gian.

                    TIẾN TRÌNH LỊCH SỬ VÀ VỊ TRÍ CỦA LUẬN ÁN
                    
 1955       1963          1967             1983            2010s           2023
 Fritz      A.N.          Lattes &         Showalter &     Ting Wei,       LÊ ĐÌNH LONG
 John       Tikhonov      Lions            Ewing           Yamamoto        (LUẬN ÁN TIẾN SĨ)
 Cauchy    Chỉnh hóa     Tựa đảo          Tựa biên        Mô hình         Mở rộng hệ
 ngược    toán tử       (Quasi-          (QBV) &         khuếch tán      Riemann-Liouville,
 thời     không         reversibility)   phương trình    phân đoạn       Caputo, Fractional
 gian     chỉnh                          Sobolev         cơ bản          Laplacian & L^q

Tại Việt Nam, các chuyên gia hàng đầu về bài toán ngược như Đặng Đức Trọng, Đinh Nho Hào, Nguyễn Huy Tuấn, Đặng Văn Vinh, Phạm Hoàng Quân đã có nhiều đóng góp quan trọng:

  • Đặng Đức Trọng, Đinh Nho Hào và Lê Đình Long (2005, 2009) đã khảo sát bài toán nhiệt ngược không thuần nhất $u_t - u_{xx} = -\beta(t)f(x)$ và $u_t - \Delta u = g(t)f(x,y)$ trên miền hai chiều $\Omega = (0,1)^2$, áp dụng biến đổi Fourier và kỹ thuật chặt cụt chuỗi để khôi phục nguồn không trơn.
  • Nguyễn Huy Tuấn và các cộng sự (2010, 2011, 2012, 2016, 2017) đã mở rộng thành công các phương pháp QBV phi địa phương, phương pháp chặt cụt Fourier cho phương trình parabolic ngược thời gian phi tuyến $u_t + Au(t) = f(t,u(t))$ và phương trình khuếch tán phân số đa chiều $\partial_t^\alpha u(x,t) = \Delta u(x,t) + F(x,t)$.
  • Đặng Đức Trọng, Nguyễn Thành Điện, Nguyễn Hữu Cần và Trần Bảo Ngọc (2019, 2020) nghiên cứu bài toán khuếch tán không thuần nhất trong thang đo Hilbert $H^\mu(\mathbb{R}^n)$ với đạo hàm Caputo và toán tử phi địa phương.

Trên bình diện quốc tế, các nghiên cứu tiêu biểu của nhóm tác giả Ting Wei & Jun-Gang Wang (2014) tại Applied Mathematical ModellingInverse Problems đã giải quyết bài toán hàm nguồn ngược cho phương trình khuếch tán cấp phân số với hệ số biến thiên thông qua phương pháp tích phân loại một; Ji-Jun Liu & Masahiro Yamamoto (2011) công bố công trình mẫu mực về bài toán ngược thời gian cho phương trình vi phân cấp phân số một chiều; Ting Wei & Jun Xian (2019) sử dụng phương pháp biến phân Tikhonov kết hợp thuật toán gradient liên hợp.

Positioning của luận án: Mặc dù các công trình quốc tế đã xử lý một số mô hình khuếch tán phân đoạn cơ bản với đạo hàm Caputo, nhưng việc phân tích toàn diện toán tử tích phân - vi phân Riemann-Liouville bậc $1-\alpha$ kết hợp số hạng nguồn phi địa phương, toán tử vi phân cấp không nguyên Fractional Laplacian $(-\Delta)^s$, và đặc biệt là bài toán giá trị cuối trong không gian phi Hilbert $L^q(\Omega)$ ($q \neq 2$) cho phương trình biharmonic vẫn còn nhiều khoảng trống học thuật lớn. Luận án của Lê Đình Long định vị chính xác tại giao điểm của các lý thuyết giải tích hàm hiện đại, thiết lập cầu nối giữa lý thuyết nửa nhóm toán tử, giải tích điều hòa và lý thuyết xấp xỉ số.


Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mang lại những đóng góp nền tảng cho lý thuyết toán giải tích phương trình đạo hàm riêng và bài toán ngược:

  1. Chứng minh tường minh tính không chỉnh theo nghĩa Hadamard: Bằng cách biểu diễn nghiệm giải tích yếu (mild solution) dưới dạng chuỗi Fourier thông qua các hàm đặc biệt Mittag-Leffler hai tham số $E_{\alpha,\beta}(z) = \sum_{k=0}^\infty \frac{z^k}{\Gamma(k\alpha+\beta)}$, luận án đã chỉ ra cơ chế suy biến phổ khi các giá trị riêng $\lambda_n \to +\infty$. Tốc độ triệt tiêu hàm mũ của toán tử giải tích thuận dẫn đến việc toán tử nghịch đảo không bị chặn, khiến bài toán khôi phục hàm nguồn và điều kiện đầu cực kỳ nhạy cảm với nhiễu dữ liệu đo.
  2. Thiết lập toán tử chỉnh hóa tối ưu cho lớp đạo hàm Riemann-Liouville: Khác với đạo hàm Caputo có tính chất bảo toàn điều kiện ban đầu cổ điển, đạo hàm Riemann-Liouville $\partial_t^{1-\alpha}$ tạo ra tính kỳ dị mạnh tại lân cận $t \to 0^+$. Luận án đã phát triển cấu trúc toán tử bù điều kiện biên cải biên, bảo đảm tính giải được duy nhất của bài toán xấp xỉ chính quy.
  3. Phát triển lý thuyết đánh giá sai số hội tụ trong không gian phi Hilbert $L^q(\Omega)$ ($q \neq 2$): Trong khi đại đa số các nghiên cứu quốc tế dừng lại ở không gian Hilbert $L^2(\Omega)$, luận án đã giải quyết thành công bài toán biên cấp bốn (biharmonic) $Mu = G(x,t,u)$ trên không gian Banach $L^q(\Omega)$ nhờ kết hợp định lý nhúng Sobolev $H_0^2(\Omega) \hookrightarrow L^q(\Omega)$ và các bất đẳng thức giải tích thực nâng cao.

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

Khung phân tích của luận án tích hợp một cách chặt chẽ ba trụ cột phương pháp luận toán học:

                            KHUNG PHÂN TÍCH TOÁN HỌC ĐỘC ĐÁO
  • Lý thuyết phổ toán tử elliptic tự liên hợp: Sử dụng toán tử đối xứng xác định dương $L: H_0^1(\Omega) \cap H^2(\Omega) \to L^2(\Omega)$ với dãy giá trị riêng thực dương tăng ngặt $0 < \lambda_1 \le \lambda_2 \le \dots \le \lambda_n \to +\infty$ và hệ hàm riêng trực chuẩn ${e_n}_{n=1}^\infty$.
  • Lý thuyết hàm Mittag-Leffler và phép biến đổi Laplace: Vận dụng hệ thức vi phân: $$\partial_t \left(t^{\alpha-1} E_{\alpha,\alpha}(-\lambda t^\alpha)\right) = t^{\alpha-2} E_{\alpha,\alpha-1}(-\lambda t^\alpha)$$ để tìm nghiệm mild chính xác dưới dạng chuỗi tích phân chập.
  • Nguyên lý bù nhiễu đa cấp (Multi-level Regularization Strategy): Xây dựng bộ ba phương pháp chính quy hóa: Phương pháp giá trị tựa biên cải biên (Modified QBV), Phương pháp chặt cụt phổ Fourier (Fourier Truncation), và Phương pháp lặp Fractional Landweber với tham số bậc không nguyên.

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ủ chặt chẽ lập trường nhận thức luận duy lý thực chứng toán học (Mathematical Rationalism / Positivism), kết hợp phép chứng minh giải tích suy diễn hình thức (Deductive Analytical Proofs) với các thuật toán kiểm chứng số học (Numerical Simulations). Thiết kế nghiên cứu được phân cấp đa tầng:

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu đảm bảo chuẩn mực giải tích hàm khắt khe bậc nhất:

  1. Thiết lập mô hình nhiễu quan sát tổng quát: Dữ liệu thực nghiệm thu nhận được từ các cảm biến vật lý luôn chứa sai số đo lường với mức độ nhiễu xác định $\delta > 0$: $$|\bar{h}\delta - h|{L^2(\Omega)} + |\bar{q}\delta - q|{L^2(0,T)} < \delta$$ hoặc đối với bài toán nguồn và thế phi tuyến: $$|\bar{f}\delta - f|{L^2(\Omega)} + |\bar{g}\delta - g|{L^2(\Omega)} + |G_\delta - G|_{L^\infty(0,T; L^2(\Omega))} < \delta$$
  2. Thiết lập điều kiện trơn nguồn (Source Condition): Giả định hàm nghiệm chính xác thuộc không gian Sobolev hoặc thang đo Hilbert $H^p(\Omega)$ ($p > 0$): $$\sum_{n=1}^\infty \lambda_n^{2p} |\langle f, e_n \rangle|^2 \le E^2 < +\infty$$ với $E > 0$ là biên tiên nghiệm của hàm nguồn hoặc điều kiện đầu.
  3. Quy tắc chọn tham số chỉnh hóa kép:
    • Quy tắc tiên nghiệm (A priori choice rule): Tham số chỉnh hóa $\epsilon = \epsilon(\delta)$ hoặc ngưỡng chặt cụt $N = N(\delta)$ được thiết lập hoàn toàn dựa trên mức nhiễu $\delta$ và tham số độ trơn $p$ để đạt tốc độ hội tụ tối ưu dạng Hölder: $$|f_\epsilon^\delta - f|_{L^2(\Omega)} \le C \delta^{\frac{p}{p+1}} E^{\frac{1}{p+1}}$$
    • Quy tắc hậu nghiệm (A posteriori choice rule - Morozov's Discrepancy Principle): Xác định tham số chỉnh hóa $\epsilon$ trực tiếp từ dữ liệu đo bị nhiễu thông qua việc giải phương trình sai dư phi tuyến mà không đòi hỏi biết trước tham số độ trơn $p$ của nghiệm chính xác.

Data và phân tích

Phần thực nghiệm số học trong luận án được lập trình và mô phỏng trên nền tảng tính toán khoa học chuyên dụng (MATLAB / Mathematica / Python NumPy):

  • Phân bố không gian - thời gian: Miền tính toán chuẩn $\Omega = (0,\pi)$ hoặc $\Omega = (0,\pi) \times (0,\pi)$, khoảng thời gian quan sát $T=1.0$.
  • Kích thước mẫu lưới rời rạc: Lưới không gian và thời gian với số điểm nút $N_x = 100$, $N_T = 100$ và $N = 50, 100, 200$.
  • Mức độ nhiễu nhân tạo: Mô phỏng các mức nhiễu ngẫu nhiên chuẩn $\delta \in {0.001, 0.005, 0.01, 0.015}$.
  • Đánh giá định lượng sai số: Sử dụng chuẩn $L^2$ rời rạc và chuẩn cực đại $L^\infty$ để đo khoảng cách sai số giữa nghiệm giải tích chính xác $f_{\text{exact}}$ và nghiệm xấp xỉ chỉnh hóa $f_{\text{reg}}^\delta$.

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

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

  1. Phát hiện 1: Thiết lập chặn dưới của sự suy biến toán tử nghịch đảo Luận án đã chứng minh định lượng rằng mẫu số trong hệ số phổ Fourier của nghiệm giảm theo hàm lũy thừa bậc phân đoạn kết hợp hàm Mittag-Leffler, thỏa mãn bất đẳng thức sắc bén: $$M_1(\lambda_1, \lambda_n) \le \int_0^T E_{\alpha,\alpha}(-\lambda_n(T-\tau)^\alpha) d\tau \le M_2(\lambda_n)$$ với $M_2(\lambda_n) \sim \mathcal{O}(\lambda_n^{-1})$ khi $n \to \infty$. Điều này khẳng định tính không chỉnh của bài toán khôi phục hàm nguồn có bậc phân đoạn $\alpha \in (0,1)$ dịu hơn bài toán truyền nhiệt cổ điển ($\alpha = 1$, suy biến hàm mũ $e^{-\lambda_n T}$), nhưng cấu trúc kỳ dị giải tích lại phức tạp hơn gấp nhiều lần.
  2. Phát hiện 2: Sự vượt trội của phương pháp Quasi-Boundary Value (QBV) Nghiệm chỉnh hóa $u_\epsilon(x,0)$ đạt được tốc độ hội tụ Logarithm bậc cao hoặc Hölder phụ thuộc trực tiếp vào tham số $\alpha$, khắc phục triệt để hiện tượng mất ổn định số tại thời điểm ban đầu $t=0$.
  3. Phát hiện 3: Tối ưu hóa bậc hội tụ của Fractional Landweber Đối với phương trình vi phân chứa toán tử Fractional Laplacian $(-\Delta)^s$ với $s \in (0,1)$, phương pháp Fractional Landweber với số bước lặp $k = k(\delta)$ lựa chọn theo nguyên tắc sai dư Morozov mang lại sai số xấp xỉ tối ưu cỡ $\mathcal{O}(\delta^{2s/(2s+d)})$, vượt trội so với các thuật toán gradient descent tiêu chuẩn.
  4. Phát hiện 4: Tính vững của nghiệm phương trình Biharmonic trong $L^q(\Omega)$ ($q \neq 2$) Mở rộng thành công kỹ thuật xấp xỉ Fourier cho bài toán cấp bốn phi tuyến với nguồn $G(x,t,u)$ thỏa mãn điều kiện Lipschitz địa phương, chứng minh sự hội tụ mạnh trong không gian $L^q(\Omega)$ ngay cả khi dữ liệu cuối chỉ thuộc $L^2(\Omega)$.

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Cung cấp khung giải tích hoàn chỉnh để xử lý các toán tử vi phân không nguyên bậc tùy ý, mở ra các công cụ phân tích hàm mới cho lớp phương trình tiến hóa dị thường.
  • Về mặt phương pháp tính: Đưa ra các thuật toán chỉnh hóa có tính ổn định cao, tốc độ tính toán nhanh, dễ dàng cài đặt song song trên các hệ thống tính toán hiệu năng cao (HPC).
  • Về mặt ứng dụng thực tiễn:
    • Địa kỹ thuật và môi trường: Khôi phục vị trí và công suất của các nguồn phát tán chất ô nhiễm ngầm trong tầng ngậm nước không đồng nhất.
    • Y sinh học và chẩn đoán hình ảnh: Mô hình hóa chính xác quá trình khuếch tán oxy và dược chất trong mô sinh học có cấu trúc fractal phức tạp.
    • Kỹ thuật nhiệt và vật liệu: Xác định phân bố nhiệt độ ban đầu trong các cấu kiện hợp kim composite chịu tải nhiệt cực hạn.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án cũng thẳng thắn chỉ ra các giới hạn nghiên cứu nội tại (Boundary conditions):

                       CÁC HƯỚNG MỞ RỘNG TRONG TƯƠNG LAI
  1. Giới hạn biến số của hàm nguồn: Luận án chủ yếu tập trung vào hàm nguồn tách biến phụ thuộc không gian $f(x)$ hoặc $f(x,y)$. Bài toán hàm nguồn phụ thuộc biến thời gian $f(t)$ hoặc hàm nguồn tổng quát toàn phần $f(x,t)$ đòi hỏi kỹ thuật đo đạc tích phân phi địa phương phức tạp hơn.
  2. Cấu trúc miền tính toán: Các thử nghiệm số chủ yếu thực hiện trên miền hình học bị chặn dạng khối chữ nhật hoặc hình hộp $\Omega \subset \mathbb{R}^d$ ($d=1,2,3$) với điều kiện biên Dirichlet thuần nhất. Miền hình học phi cấu trúc phức tạp cần đến sự hỗ trợ của phương pháp phần tử hữu hạn (FEM) phân đoạn.
  3. Mô hình nhiễu tất định: Dữ liệu quan sát giả định chịu mức nhiễu bị chặn trong chuẩn $L^2(\Omega)$. Trong thực tế, nhiễu đo lường thường mang tính chất ngẫu nhiên thống kê (Gaussian white noise).

Chương trình nghiên cứu 4 hướng tương lai do tác giả vạch ra:

  • Hướng 1: Khảo sát bài toán khôi phục hàm nguồn theo biến thời gian $f(t)$ từ dữ liệu tích phân năng lượng.
  • Hướng 2: Mở rộng nghiên cứu cho các loại đạo hàm phân đoạn không kỳ dị mới xuất hiện như đạo hàm Atangana-Baleanu Caputo (ABC), đạo hàm Caputo-Fabrizio (CF), và đạo hàm khả vi conformable.
  • Hướng 3: Khảo sát mô hình khuếch tán với hàm nguồn phi tuyến mạnh không thỏa mãn điều kiện Lipschitz toàn cục trong không gian $L^q$.
  • Hướng 4: Phát triển lý thuyết chỉnh hóa thống kê cho phương trình đạo hàm riêng phân đoạn ngẫu nhiên (Stochastic Fractional PDEs).

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

Luận án của NCS. Lê Đình Long tạo nên những xung lực phát triển học thuật mạnh mẽ cả trong nước và quốc tế:

  • 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í uy tín hàng đầu thuộc danh mục ISI (Mathematical Methods in the Applied Sciences, Computational and Applied Mathematics, Journal of Inequalities and Applications, Acta Mathematica Sinica, AIMS Mathematics) với chỉ số trích dẫn ước tính tăng trưởng liên tục trong cộng đồng nghiên cứu Inverse Problems toàn cầu.
  • Nâng cao vị thế cơ sở đào tạo: Đóng góp trực tiếp vào thành tích nghiên cứu khoa học xuất sắc của Khoa Toán - Tin học, Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia TP. Hồ Chí Minh.
  • Giá trị chuyển giao học thuật: Các thuật toán chỉnh hóa số đề xuất trong luận án có thể tích hợp trực tiếp vào các phần mềm mô phỏng vật lý địa kỹ thuật, công nghệ vật liệu mới và các gói thư viện giải thuật toán ứng dụng mã nguồn mở.

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

  • Nghiên cứu sinh và học viên cao học ngành Toán Giải tích, Toán Ứng dụng: Tiếp cận hệ thống phương pháp luận hoàn chỉnh, tài liệu tham khảo chuẩn mực về giải tích phân đoạn và kỹ thuật chỉnh hóa bài toán không chỉnh.
  • Các nhà nghiên cứu giải tích phương trình đạo hàm riêng: Thừa hưởng các đánh giá bất đẳng thức sắc bén liên quan đến hàm Mittag-Leffler, toán tử Fractional Laplacian và phép nhúng Sobolev.
  • Kỹ sư R&D trong các ngành công nghệ tính toán, địa vật lý, khí tượng: Sử dụng các thuật toán khôi phục dữ liệu thực nghiệm bị nhiễu phục vụ dự báo lan truyền ô nhiễm và truyền nhiệt phi tuyến.
  • Các giảng viên đại học: Giáo trình tham khảo chuyên đề cao học về "Bài toán ngược và bài toán không chỉnh trong vật lý toán".

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

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc thiết lập chặn dưới giải tích sắc bén của toán tử giải tích thuận cho phương trình chứa đạo hàm Riemann-Liouville bậc $1-\alpha$, đồng thời phát triển phương pháp Quasi-Boundary Value mở rộng trong không gian $L^q(\Omega)$ ($q \neq 2$). Kết quả này đã mở rộng trực tiếp lý thuyết chỉnh hóa tựa biên kinh điển của Showalter (1983) và phương pháp phương trình Sobolev của Ewing (1983) từ bài toán nguyên cổ điển sang bài toán phân đoạn không nguyên phi tuyến đa chiều.

2. Sự đổi mới về mặt phương pháp luận so với các nghiên cứu quốc tế trước đây?

So với các nghiên cứu của Ting Wei & Jun-Gang Wang (2014) hay Ji-Jun Liu & Masahiro Yamamoto (2011) vốn chỉ tập trung vào đạo hàm Caputo trên miền 1D hoặc không gian $L^2$, luận án đã:

  • Xử lý thành công đạo hàm Riemann-Liouville có độ kỳ dị cao tại $t=0$.
  • Tích hợp thành công thuật toán Fractional Landweber với quy tắc chọn tham số hậu nghiệm Morozov cho toán tử Fractional Laplacian không bị chặn $(-\Delta)^s$.
  • Đưa ra chứng minh chặt chẽ cho sự hội tụ của nghiệm trong thang đo không gian Banach $L^q(\Omega)$.

3. Phát hiện bất ngờ nhất về mặt hành vi số học của mô hình là gì?

Phát hiện bất ngờ nhất là khi tham số bậc phân đoạn $\alpha \to 1^-$, nghiệm chỉnh hóa của phương trình khuếch tán phân đoạn hội tụ mượt mà và bảo toàn liên tục về nghiệm của bài toán truyền nhiệt parabolic cổ điển, nhưng tốc độ hội tụ số thực nghiệm khi $\alpha$ nhỏ ($\alpha = 0.3$) lại ổn định hơn và ít nhạy cảm với mức nhiễu $\delta = 0.015$ hơn so với khi $\alpha$ tiến gần 1.

4. Luận án có cung cấp đầy đủ giao thức tái lập kết quả (Replication Protocol) không?

Có. Luận án cung cấp chi tiết toàn bộ các công thức biểu diễn nghiệm giải tích, hệ số phân rã phổ Fourier, cấu trúc toán tử xấp xỉ, các bước lặp thuật toán, quy tắc xác định tham số tiên nghiệm/hậu nghiệm và các thông số cụ thể của lưới sai phân ($N=50, 100$, mức nhiễu $\delta=0.015$), cho phép tái lập 100% các kết quả mô phỏng số.

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 nghiên cứu 10 năm tập trung vào việc giải quyết trọn vẹn 4 hướng mở: chuyển từ bài toán xác định nguồn không gian sang nguồn thời gian $f(t)$; mở rộng sang các toán tử vi phân không nguyên nhân phi kỳ dị (ABC, CF); xử lý bài toán trong miền hình học thực tế 3D với dữ liệu nhiễu ngẫu nhiên Gauss; và phát triển các thuật toán học máy giải tích kết hợp mạng nơ-ron vật lý (Physics-Informed Neural Networks - PINNs) để giải bài toán ngược quy mô lớn.


Kết luận

Luận án tiến sĩ "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" của tác giả Lê Đình Long là một công trình khoa học xuất sắc, mẫu mực và đồ sộ về mặt giải tích toán học. Những đóng góp vượt bậc của luận án được đúc kết qua các điểm then chốt:

  1. Chứng minh trọn vẹn bản chất không chỉnh theo nghĩa Hadamard cho 5 lớp phương trình vi phân phân đoạn và phương trình biharmonic cấp cao với đạo hàm Riemann-Liouville, Caputo và Fractional Laplacian.
  2. Đề xuất và làm chủ ba phương pháp chỉnh hóa hiện đại: Phương pháp giá trị tựa biên (QBV), phương pháp chặt cụt chuỗi Fourier, và phương pháp Fractional Landweber.
  3. Thiết lập tốc độ hội tụ tối ưu của nghiệm chỉnh hóa về nghiệm chính xác dưới cả hai chiến lược lựa chọn tham số tiên nghiệm (a priori) và hậu nghiệm (a posteriori) trong không gian $L^2(\Omega)$ và $L^q(\Omega)$ ($q \neq 2$).
  4. Kiểm chứng tính đúng đắn của lý thuyết thông qua hệ thống thử nghiệm số học nghiêm ngặt, minh chứng tính khả thi và độ chính xác cao của các thuật toán đề xuất dưới các mức nhiễu thực tế $\delta = 0.015$.
  5. Công bố 5 bài báo quốc tế chất lượng cao thuộc danh mục ISI (Q1, Q2), khẳng định năng lực nghiên cứu độc lập, chuyên sâu và vị thế khoa học xuất sắc của tập thể tác giả và cơ sở đào tạo trên trường quốc tế.
  6. Mở ra 4 hướng nghiên cứu liên ngành đột phá, đặt nền móng lý thuyết vững chắc cho các ứng dụng thực tiễn trong công nghệ mô phỏng truyền chất, truyền nhiệt và chẩn đoán y sinh học hiện đại.