Tổng quan về luận án

Bài toán xác định nguồn (inverse source problem) cho phương trình đạo hàm riêng loại parabolic là một trong những chủ đề trọng tâm của lý thuyết bài toán ngược hiện đại, đóng vai trò mô hình hóa cốt lõi cho các quá trình truyền nhiệt nghịch, lan truyền chất ô nhiễm trong khí quyển và thủy vực, khuếch tán sinh học và tài chính lượng tử. Về bản chất toán học, đây là bài toán đặt không chỉnh theo nghĩa Hadamard (Hadamard ill-posed problem). Như văn bản luận án đã định nghĩa: "Một bài toán được gọi là đặt chỉnh theo nghĩa Hadamard nếu nó thỏa mãn cả ba điều kiện: i) Nghiệm của bài toán luôn tồn tại; ii) Nghiệm của bài toán là duy nhất; iii) Nghiệm phụ thuộc liên tục vào dữ liệu của bài toán. Nếu ít nhất một trong ba điều kiện trên không được thỏa mãn, thì bài toán được gọi là đặt không chỉnh". Khi bài toán không ổn định, một nhiễu loạn ngẫu nhiên vô cùng nhỏ $\delta > 0$ trong dữ liệu đo đạc thực nghiệm $g^\delta(x)$ tại thời điểm cuối $t = T$ sẽ dẫn đến sai số bùng nổ vô hạn trong nghiệm nguồn cần khôi phục $f(x)$.

Dữ liệu cuối t=T: g^δ(x) = u(x,T) + Nhiễu δ

Khoảng trống học thuật (research gap) cốt lõi mà luận án của tác giả Lương Duy Nhật Minh (Trường Đại học Vinh, dưới sự hướng dẫn của PGS. Nguyễn Văn Đức và TS. Nguyễn Trung Thành) giải quyết gồm ba bài toán biên - giá trị ban đầu chưa từng có lời giải hoàn chỉnh trong y văn quốc tế:

  1. Bài toán 1 (ISP1): Xác định nguồn nhiệt dạng tách biến $f(x)h(t)$ cho phương trình parabolic với hệ số khuếch tán phụ thuộc thời gian $a(t)$ trong không gian Hilbert vô hạn chiều $L^2(\mathbb{R}^n)$: $$\frac{\partial u}{\partial t}(x, t) = a(t)\Delta u(x, t) + f(x)h(t), \quad (x,t) \in \mathbb{R}^n \times (0, T)$$
  2. Bài toán 2 (ISP2): Xác định thành phần nguồn $f(x)$ cho phương trình parabolic bậc phân theo cả biến thời gian $\alpha \in (0, 1)$ (đạo hàm Caputo) và biến không gian $\gamma > 0$ (toán tử Laplace bậc phân $(-\Delta)^{\gamma/2}$): $$\partial_t^\alpha u(x, t) + (-\Delta)^{\frac{\gamma}{2}} u(x, t) = f(x)h(t), \quad (x,t) \in \mathbb{R}^n \times (0, T)$$
  3. Bài toán 3 (ISP3): Xác định nguồn $F \in X$ cho phương trình parabolic trừu tượng trong không gian Banach $X$ với toán tử sinh $-A$ tạo ra một nửa nhóm giải tích ${S(t)}_{t \ge 0}$: $$u'(t) + Au(t) = F, \quad t \in (0, T); \quad u(0) = 0; \quad u(T) = g \in X$$

Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học được chuẩn hóa:

  • RQ1: Tồn tại hay không các đánh giá ổn định có điều kiện (conditional stability estimates) kiểu Hölder cho các bài toán xác định nguồn khi hệ số khuếch tán phụ thuộc thời gian $a(t)$ hoặc khi toán tử vi phân mang tính dị thường bậc phân?
  • RQ2: Phương pháp chỉnh hóa làm nhuyễn (mollification regularization method) với nhân Dirichlet nhiều chiều có thiết lập được tốc độ hội tụ kiểu Hölder tối ưu dưới các luật chọn tham số tiên nghiệm (a priori) và hậu nghiệm (a posteriori)?
  • RQ3: Làm thế nào để khắc phục sự sụp đổ của phép biến đổi Fourier trong không gian Banach trừu tượng nhằm chỉnh hóa bài toán ngược xác định nguồn thông qua lý thuyết nửa nhóm toán tử?

Khung lý thuyết nền tảng tích hợp giải tích hàm hiện đại, lý thuyết nửa nhóm toán tử tuyến tính (analytic semigroups), giải tích Fourier, giải tích bậc phân (fractional calculus) với các hàm đặc biệt Mittag-Leffler ($E_{\alpha,\beta}$), định lý biểu diễn Paley-Wiener và nguyên lý bù sai số Morozov. Đóng góp đột phá của luận án được định lượng qua 03 công trình công bố trên các tạp chí quốc tế thuộc danh mục SCIE uy tín cao: Inverse Problems in Science and Engineering (SCIE, IF: 1.314), Applicable Analysis (SCIE, IF: 1.107) và Applied Numerical Mathematics (SCIE, IF: 1.290).

Literature Review và Positioning

Nghiên cứu về bài toán xác định nguồn của phương trình parabolic khởi nguồn từ các công trình đặt nền móng của Cannon (1968, [13]), chứng minh tính duy nhất nghiệm nguồn nhiệt $f(x)$ với quan sát trên một phần biên $\partial u / \partial n = g$. Tuy nhiên, Cannon (1968) đã chỉ ra rằng bằng cách chọn dãy dữ liệu quan sát biên hội tụ đều về 0 ($a_n = \sqrt{n}$), chuẩn của nghiệm nguồn $|f_n|_{L^2}$ bùng nổ tiến ra vô cùng, minh chứng trực quan cho tính đặt không chỉnh nghiêm trọng. Tiếp đó, Yamamoto (1993, [130]) khảo sát nguồn dạng tách biến $\sigma(t)f(x_1, x_2)$ trong miền bị chặn hai chiều $\Omega = (0,1)^2$ với $\sigma(0) \neq 0$, chứng minh tính duy nhất và ổn định nghiệm dựa trên nguyên lý cực đại và mở rộng giải tích.

       LỊCH SỬ PHÁT TRIỂN & ĐỊNH VỊ HỌC THUẬT CỦA LUẬN ÁN

  Cannon (1968), Yamamoto (1993)
  [Nguồn nhiệt cổ điển, miền bị chặn, hệ số hằng số]
  Yi & Murio (2004), Đặng Đức Trọng và cs (2006, 2009)
  [Làm nhuyễn 1D-2D, Fourier truncation, miền bị chặn]
  Dou, Fu, Yang (2009); Hasanov (2012, 2013); Zhang & Wei (2013)
  [Chỉnh hóa Tikhonov, Tựa đảo, Gradient liên hợp, Bậc phân thời gian 1D]

Trong dòng chảy phương pháp chỉnh hóa, Yi & Murio (2004, [137, 138]) tiên phong ứng dụng phương pháp làm nhuyễn Gauss/Dirichlet cho bài toán một và hai chiều. Đặng Đức Trọng cùng các cộng sự (2006, 2009, [120, 121]) phát triển phương pháp chặt cụt chuỗi và tích phân Fourier cho nguồn $\varphi(t)f(x,y)$ trong miền compact. Dou, Fu và Yang (2009, [27]) áp dụng phương pháp tựa đảo (quasi-reversibility) cho phương trình truyền nhiệt 1D nhưng chỉ thu được đánh giá tiên nghiệm. Hasanov và Pektas (2012, 2013, [55, 56]) ứng dụng phương pháp biến phân và gradient liên hợp nhưng gặp khó khăn khi chứng minh tốc độ hội tụ giải tích.

Về phương trình vi phân bậc phân, Zhang và Wei (2013, [127]) và Yang cùng cộng sự (2014, 2015, [133, 134, 135]) khảo sát bài toán truyền nhiệt bậc phân theo thời gian $\partial_t^\alpha u - u_{xx} = f(t)$ trong không gian 1D bằng phương pháp làm nhuyễn và Fourier cutoff. Tuy nhiên, kết quả hậu nghiệm trong Yang et al. (2014) chỉ đạt đánh giá sai số tốc độ chậm kiểu Logarit ($(\ln(1/\delta))^{-\mu}$). Mặt khác, các công trình của Li et al. (2013, [83]) và Sakamoto & Yamamoto (2011, [110]) nghiên cứu phương trình bậc phân không gian - thời gian nhưng không đưa ra được thuật toán số và ví dụ kiểm chứng thực nghiệm.

Đối với không gian Banach trừu tượng, các công bố kinh điển của Prilepko, Piskarev, Shaw (2000, [101]) và Hasanov (2013, [57]) chỉ chứng minh sự tồn tại, duy nhất nghiệm mà hoàn toàn bỏ ngỏ đánh giá tốc độ hội tụ và mô phỏng số.

Tiêu chí so sánh Dou, Fu & Yang (2009) [27] Yang, Fu et al. (2014, 2015) [133, 135] Luận án của Lương Duy Nhật Minh (2020)
Không gian hình học Miền 1D hữu hạn $(0, 1)$ Nửa trục 1D $[0, \infty)$ Toàn không gian Hilbert $L^2(\mathbb{R}^n)$ & Banach $X$
Cấu trúc phương trình Parabolic cổ điển, hệ số hằng Bậc phân chỉ theo biến thời gian $\partial_t^\alpha$ Bậc phân toàn phần $\partial_t^\alpha + (-\Delta)^{\gamma/2}$ & Hệ số $a(t)$
Đánh giá sai số hậu nghiệm Không thiết lập được Tốc độ chậm kiểu Logarit Tốc độ nhanh kiểu Hölder tối ưu
Khung không gian Banach Không áp dụng được Không áp dụng được Giải quyết trọn vẹn qua nửa nhóm giải tích

Đó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 đặt không chỉnh của Tikhonov & Arsenin (1977) và Lavrentiev (1967), xác lập các bất đẳng thức ổn định có điều kiện (conditional stability estimates) kiểu Hölder trên các tập compact hàm trơn Sobolev $H^p(\mathbb{R}^n)$ và tập hàm nguyên dạng mũ $||\cdot||_q$.

Về mặt giải tích điều hòa và không gian Sobolev, nếu dữ liệu nguồn thực tế thỏa mãn điều kiện tiên nghiệm $|f|{H^p} \le E_p$, với dữ liệu nhiễu $|g - g^\delta|{L^2} \le \delta$, đánh giá ổn định đạt được có dạng: $$|f|_{L^2} \le C_1^{\frac{p}{2(p+2)}} E_p^{\frac{2}{p+2}} \delta^{\frac{p}{p+2}}$$ với hằng số $C_1 = C_1(a, h, T) > 0$. Khi nghiệm thuộc không gian các hàm giải tích dạng mũ $|f|q \le E_q$, tính ổn định chuyển sang dạng Hölder mũ: $$|f|{L^2} \le (H(0))^{-\frac{q}{q+b}} \delta^{\frac{q}{q+b}} E_q^{\frac{b}{q+b}}, \quad \text{trong đó } b = \int_0^T a(t)dt, \quad H(0) = \int_0^T h(t)dt$$

                   CẤU TRÚC ĐÁNH GIÁ ỔN ĐỊNH KIỂU HÖLDER
                                
   Dữ liệu nhiễu đo đạc: ‖g - g^δ‖_{L²} ≤ δ
  Giả thiết Tiên nghiệm 1:       Giả thiết Tiên nghiệm 2:
  f ∈ Hᵖ(Rⁿ), ‖f‖_{Hᵖ} ≤ E_p      |||f|||_q ≤ E_q (Hàm giải tích)
  Ổn định kiểu lũy phân:        Ổn định kiểu Hölder mũ:
  ‖f‖ ≤ C₁ · E_p^(2/(p+2)) · δ^(p/(p+2))   ‖f‖ ≤ H(0)^(-q/(q+b)) · δ^(q/(q+b)) · E_q^(b/(q+b))

Đối với phương trình bậc phân toàn phần không - thời gian, luận án đã vận dụng các tính chất giải tích sâu của hàm Mittag-Leffler hai tham số $E_{\alpha, \beta}(z) = \sum_{k=0}^\infty \frac{z^k}{\Gamma(\alpha k + \beta)}$ để kiểm soát toán tử chuyển tiếp giải tích, chứng minh rằng tốc độ hội tụ của nghiệm chỉnh hóa đạt bậc tối ưu (optimal order) trên tập nguồn bị chặn: $$\Delta(\delta, S_\delta) \le c \inf_S \Delta(\delta, S) \quad (c \ge 1)$$

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

Khung phân tích của luận án tích hợp chặt chẽ 4 trụ cột lý thuyết toán học:

  1. Lý thuyết nửa nhóm giải tích (Analytic Semigroup Theory - Pazy 1983, Engel & Nagel 2000): Xây dựng biểu diễn nghiệm qua toán tử tích phân nửa nhóm $u(t) = \int_0^t S(t-s)F ds$, chuyển hóa toán tử không bị chặn $A$ thành toán tử chỉnh hóa bị chặn dạng resolvent $(I + \alpha A^p)^{-1}$.
  2. Lý thuyết biến đổi Fourier và Định lý Paley-Wiener: Xây dựng không gian các hàm nguyên dạng mũ $M_{\nu, 2}(\mathbb{R}^n)$ có biến đổi Fourier với giá mang compact trong lăng trụ đa chiều $\Delta_\nu = [-\nu, \nu]^n$.
  3. Lý thuyết xấp xỉ nhân Dirichlet đa chiều: Định nghĩa toán tử làm nhuyễn đa biến dạng tích chập: $$S_\nu(g)(x) := \left(\frac{2}{\pi}\right)^{\frac{n}{2}} (D_\nu * g)(x) = \frac{1}{\pi^n} \int_{\mathbb{R}^n} \left(\prod_{j=1}^n \frac{\sin(\nu y_j)}{y_j}\right) g(x-y) dy$$ Toán tử này hoạt động như một bộ lọc thông thấp (low-pass filter) lý tưởng: $$\mathcal{F}S_\nu(g) = \widehat{g}(\xi) \cdot \chi_{M_\nu}(\xi)$$ triệt tiêu toàn bộ các thành phần tần số cao $|\xi| > \nu$ vốn là nơi sai số nhiễu bị khuếch đại vô hạn.
  4. Giải tích vi tích phân bậc phân (Fractional Calculus): Vận dụng đạo hàm Caputo bậc $\alpha \in (0, 1)$ và toán tử vi phân giả (pseudodifferential operators) Laplace phân số $(-\Delta)^{\gamma/2}$ có biểu diễn Fourier $\mathcal{F}((-\Delta)^{\gamma/2} v)(\xi) = |\xi|^\gamma \widehat{v}(\xi)$.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ lập trường triết học duy lý bản thể luận và chuẩn mực tiên nghiệm của Toán học thuần túy kết hợp Toán ứng dụng (Positivism / Mathematical Formalism). Toàn bộ các định lý về tính tồn tại, duy nhất, đánh giá ổn định và tốc độ hội tụ đều được suy diễn logic nghiêm ngặt (rigorous mathematical proofs).

Thiết kế đa tầng (Multi-level theoretical design) bao phủ từ không gian hàm Hilbert cổ điển $L^2(\mathbb{R}^n)$, mở rộng lên không gian Sobolev tổng quát $H^p(\mathbb{R}^n)$, không gian hàm giải tích $M_{\eta, 2}(\mathbb{R}^n)$ và đạt mức tổng quát hóa cao nhất trên không gian Banach vô hạn chiều $X$.

     THIẾT KẾ ĐA TẦNG CỦA KHUNG PHƯƠNG PHÁP NGHIÊN CỨU

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu gồm 4 giai đoạn chuẩn hóa:

  1. Biểu diễn giải tích nghiệm: Chuyển đổi bài toán vi phân parabolic thuận sang phương trình đại số trong miền tần số Fourier hoặc biểu diễn dạng toán tử tích phân nửa nhóm.
  2. Thiết lập bất đẳng thức ổn định: Sử dụng bất đẳng thức Hölder giải tích hàm, bất đẳng thức Minkowski và tính chất đơn điệu của hàm Mittag-Leffler để chặn cận dưới của toán tử truyền tải $I(\theta, \tau) = \int_0^T e^{-\tau \int_s^T a(t)dt} h(s)ds$.
  3. Thiết kế toán tử chỉnh hóa: Xây dựng nghiệm xấp xỉ $f^\delta_\nu(x)$ từ dữ liệu đo chứa nhiễu $g^\delta(x)$.
  4. Xác lập luật chọn tham số chỉnh hóa $\nu(\delta)$:
    • Luật tiên nghiệm (A priori): Chọn tham số $\nu$ chỉ phụ thuộc mức độ nhiễu $\delta$, cụ thể: $$\nu = \left(\frac{1}{2nb} \ln\left(\frac{E_q}{\delta}\right)\right)^{\frac{1}{2}}$$
    • Luật hậu nghiệm (A posteriori - Nguyên lý bù sai số Morozov): Chọn tham số dừng $\nu(\delta, g^\delta)$ sao cho sai dư đạt ngưỡng sai số: $$|S_\nu(g^\delta) - g^\delta|_{L^2} = \tau \delta, \quad \tau > 1$$

Data và phân tích

Trong khoa học bài toán ngược tính toán, dữ liệu thực nghiệm được mô phỏng số (synthetic data generation) với các mức nhiễu kiểm soát nghiêm ngặt. Dữ liệu chính xác $g(x)$ được tạo từ nghiệm giải tích, sau đó pha tạp nhiễu Gauss ngẫu nhiên: $$g^\delta(x) = g(x) + \delta \cdot \text{rand}(x)$$ với các mức sai số tương đối $\epsilon = \frac{|g - g^\delta|{L^2}}{|g|{L^2}} \in {10^{-1}, 10^{-2}, 10^{-3}, 10^{-4}}$.

Thuật toán số được lập trình trên nền tảng MATLAB và Python (NumPy, SciPy), tận dụng thuật toán Biến đổi Fourier Nhanh (Fast Fourier Transform - FFT/IFFT) trên lưới không gian $N = 256 \times 256$ hoặc $512 \times 512$ điểm nút. Sai số tuyệt đối $L^2$ và sai số tương đối (Relative Error) giữa nghiệm chính xác $f_{\text{exact}}$ và nghiệm chỉnh hóa $f^\delta_{\text{reg}}$ được tính toán để kiểm chứng tốc độ hội tụ lý thuyết.

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

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

  1. Chứng minh tính duy nhất và ổn định Hölder cho hệ số $a(t)$ biến thiên: Luận án chứng minh rằng tính chất đặt không chỉnh của phương trình parabolic không bị phá vỡ bởi hệ số phụ thuộc thời gian $a(t)$ (miễn là thỏa mãn $0 < \underline{a} \le a(t) \le \overline{a}$). Nghiệm nguồn $f(x)$ phụ thuộc liên tục kiểu Hölder vào dữ liệu cuối $g(x)$.
  2. Khám phá tốc độ hội tụ kiểu Hölder thay vì Logarit cho bài toán bậc phân: Khác với các công bố trước đây của Yang et al. (2014) chỉ đạt tốc độ Logarit khi dùng luật chọn tham số hậu nghiệm, luận án chứng minh rằng phương pháp làm nhuyễn Dirichlet vẫn đảm bảo tốc độ hội tụ kiểu Hölder lũy thừa cho phương trình bậc phân không - thời gian.
  3. Đạt bậc hội tụ tối ưu (Optimal Order): Luận án chứng minh sai số lớn nhất $\Delta(\delta, S_\nu)$ giữa nghiệm chỉnh hóa và nghiệm thực tế trên tập compact nguồn $M_{\phi, E}$ đạt chặn dưới lý thuyết thông tin (information-theoretic lower bound), nghĩa là không một thuật toán chỉnh hóa nào khác có thể đạt bậc hội tụ cao hơn.
  4. Giải pháp nửa nhóm giải tích cho không gian Banach: Thiết lập thành công công thức chỉnh hóa trong không gian Banach trừu tượng $X$ mà không cần bất kỳ giả định nào về tính Hilbert hay sự tồn tại của tích vô hướng: $$F_\alpha^\delta = A(I - S(T))^{-1} (I + \alpha A^p)^{-1} g^\delta$$ đạt độ hội tụ $|F - F_\alpha^\delta|_X = \mathcal{O}(\delta^{\frac{p}{p+1}})$.
                      SO SÁNH HIỆU NĂNG HỘI TỤ
                                
  Sai số ‖f_exact - f_reg‖

Implications đa chiều

  • Về mặt lý thuyết: Thống nhất khung lý thuyết chỉnh hóa làm nhuyễn cho các lớp toán tử vi phân parabolic phi cổ điển; mở ra hướng tiếp cận xử lý toán tử sinh nửa nhóm trong không gian Banach.
  • Về mặt phương pháp luận: Cung cấp quy trình chọn tham số chỉnh hóa tự động hóa hoàn toàn bằng thuật toán giải số, loại bỏ tính chủ quan khi ước lượng tham số làm nhuyễn $\nu$.
  • Về ứng dụng thực tiễn:
    • Môi trường: Định vị chính xác tọa độ và cường độ xả thải của các nguồn gây ô nhiễm hóa chất/phóng xạ trong nước ngầm và không khí dựa trên dữ liệu trạm quan trắc nồng độ tại thời điểm $T$.
    • Kỹ thuật nhiệt & Vật liệu: Khôi phục sự phân bố nguồn nhiệt nội tại trong vật liệu composite có độ dẫn nhiệt biến đổi theo thời gian $a(t)$.
    • Y sinh học: Mô hình hóa sự khuếch tán bất thường (anomalous diffusion) của protein và tế bào ung thư tuân theo phương trình đạo hàm bậc phân.

Limitations và Future Research

  1. Giới hạn hình học toàn không gian $\mathbb{R}^n$: Phương pháp làm nhuyễn Dirichlet dựa trên biến đổi Fourier liên tục đòi hỏi miền xác định là toàn không gian $\mathbb{R}^n$. Việc áp dụng trực tiếp cho các miền hình học phức tạp hữu hạn $\Omega \subset \mathbb{R}^n$ với điều kiện biên hỗn hợp (Robin, Wentzell) đòi hỏi phải kết hợp phương pháp phần tử hữu hạn (FEM) hoặc mở rộng biên giải tích.
  2. Giả định nguồn tách biến: Hàm nguồn trong các mô hình được giả thiết có dạng tách biến $f(x)h(t)$ với $h(t)$ đã biết trước. Trường hợp hàm nguồn tổng quát phi tuyến $F(x, t, u(x,t))$ chưa được khảo sát trọn vẹn trong luận án.
  3. Chi phí tính toán trong không gian số chiều lớn: Khi số chiều $n \ge 3$, việc tính tích chập $n$-chiều của nhân Dirichlet $\prod_{j=1}^n \frac{\sin(\nu y_j)}{y_j}$ gặp hiện tượng bùng nổ số chiều (curse of dimensionality), đòi hỏi các kỹ thuật tính toán lưới thích nghi song song.

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

  • Mở rộng phương pháp sang phương trình parabolic phi tuyến và nửa tuyến tính (semilinear parabolic equations).
  • Phát triển thuật toán chỉnh hóa mạng nơ-ron sâu kết hợp vật lý (Physics-Informed Neural Networks - PINNs) để giải bài toán ngược xác định nguồn trong miền dị vật đục lỗ.
  • Khảo sát bài toán xác định đồng thời cả hệ số khuếch tán $a(t)$ và nguồn nhiệt $f(x)$ từ dữ liệu đo đạc biên bổ sung.

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

  • Tác động học thuật (Academic Impact): Các bài báo công bố từ luận án trên Inverse Problems in Science and Engineering, Applicable Analysis, và Applied Numerical Mathematics tạo dựng nền tảng trích dẫn vững chắc cho cộng đồng nghiên cứu bài toán ngược và giải tích bậc phân quốc tế. Dự báo công trình sẽ nhận được sự quan tâm và trích dẫn liên tục từ các nhóm nghiên cứu tại Hoa Kỳ, Đức, Nhật Bản, Trung Quốc và Việt Nam.
  • Chuyển giao công nghệ & Công nghiệp (Industry R&D): Cung cấp thuật toán cốt lõi cho các phần mềm mô phỏng địa kỹ thuật, mô hình hóa dòng chảy môi trường (như MODFLOW, COMSOL Multiphysics) trong module phân tích ngược nguồn ô nhiễm.
  • Ảnh hưởng chính sách & Xã hội: Cung cấp cơ sở khoa học định lượng chuẩn xác giúp các cơ quan quản lý tài nguyên môi trường truy cứu nguồn gốc sự cố xả thải công nghiệp trái phép dựa trên dấu vết quan trắc nồng độ ô nhiễm trễ.

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

                     MA TRẬN ĐỐI TƯỢNG HƯỞNG LỢI

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

1. Đóng góp lý thuyết độc đáo và có tính đột phá nhất của luận án là gì? Đó là việc thiết lập thành công các đánh giá ổn định có điều kiện kiểu Hölder và chứng minh tốc độ hội tụ đạt bậc tối ưu (optimal convergence order) cho bài toán parabolic bậc phân toàn phần không - thời gian $\partial_t^\alpha u + (-\Delta)^{\gamma/2} u = f(x)h(t)$ trong toàn không gian $L^2(\mathbb{R}^n)$. Luận án đã vượt qua rào cản kỹ thuật của tính phi địa phương (non-local property) của đạo hàm Caputo và toán tử Laplace phân số bằng cách khai thác triệt để các tính chất đơn điệu hoàn toàn của hàm giải tích Mittag-Leffler $E_{\alpha, 1}(-\lambda \xi^\gamma)$.

2. Đổi mới phương pháp luận của luận án khi so sánh với ít nhất 2 công trình quốc tế tiêu biểu?

  • So với công trình của Dou, Fu và Yang (2009, [27]): Công trình [27] chỉ áp dụng phương pháp tựa đảo cho phương trình truyền nhiệt 1D với luật chọn tham số tiên nghiệm, không có đánh giá hậu nghiệm. Luận án giải quyết trọn vẹn bài toán $n$-chiều với hệ số biến thiên $a(t)$, đồng thời thiết lập hoàn chỉnh cả luật tiên nghiệm và hậu nghiệm đạt tốc độ Hölder.
  • So với công trình của Yang và cộng sự (2014, [133]): Công trình [133] khảo sát bài toán bậc phân 1D nhưng với luật hậu nghiệm chỉ đạt tốc độ hội tụ Logarit cực kỳ chậm $\mathcal{O}((\ln(1/\delta))^{-\mu})$. Luận án đã cải tiến kỹ thuật chặn tích phân Fourier, nâng tốc độ hội tụ hậu nghiệm lên bậc Hölder lũy thừa $\mathcal{O}(\delta^{\frac{p}{p+2}})$, tạo bước nhảy vọt về hiệu năng số.

3. Phát hiện thực nghiệm/toán học nào gây bất ngờ nhất trong quá trình nghiên cứu? Phát hiện bất ngờ nhất là phương pháp làm nhuyễn Dirichlet – vốn được xem là công cụ thuần túy của không gian Fourier Hilbert – có thể kết nối hoàn hảo với lý thuyết xấp xỉ hàm nguyên dạng mũ thông qua Định lý Paley-Wiener để xử lý bài toán có hệ số khuếch tán suy biến theo thời gian. Sự kết hợp này cho phép bảo toàn cấu trúc giải tích của nghiệm mà không làm phát sinh hiện tượng giả dao động Runge ở biên tần số.

4. Luận án có cung cấp đầy đủ giao thức tái lập nghiên cứu (Replication Protocol) không? Có. Toàn bộ cấu trúc thuật toán giải số, quy trình rời rạc hóa lưới không gian, công thức biến đổi Fourier rời rạc (FFT), thuật toán lặp tìm tham số dừng hậu nghiệm $\nu(\delta)$ theo nguyên lý Morozov, và các bộ tham số mô phỏng thử nghiệm số đều được mô tả chi tiết, rõ ràng trong các mục thuật toán của từng chương, cho phép tái lập 100% kết quả tính toán trên MATLAB/Python.

5. Bản đồ lộ trình nghiên cứu 10 năm (10-year research agenda) được vạch ra như thế nào? Lộ trình 10 năm gồm 3 giai đoạn:

  • Giai đoạn 1 (Năm 1-3): Mở rộng lý thuyết sang phương trình parabolic với toán tử elliptic cấp hai tổng quát $L(x,t)u = \sum a_{ij}(x,t) u_{x_i x_j}$ trong miền bị chặn $\Omega$ với biên Lipschitz.
  • Giai đoạn 2 (Năm 4-6): Nghiên cứu bài toán ngược nguồn ngẫu nhiên (stochastic inverse source problems) khi dữ liệu biên bị nhiễu bởi quá trình Wiener / chuyển động Brown phân số.
  • Giai đoạn 3 (Năm 7-10): Tích hợp phương pháp làm nhuyễn và lý thuyết nửa nhóm vào kiến trúc Deep Learning (Operator Learning / DeepONet / PINNs) để xây dựng hệ thống phần mềm giải bài toán ngược thời gian thực cho mạng lưới cảm biến thông minh.

Kết luận

  1. Xác lập hệ thống đánh giá ổn định có điều kiện kiểu Hölder cho ba lớp bài toán xác định nguồn parabolic phức tạp: phương trình có hệ số biến thiên theo thời gian $a(t)$ trong $L^2(\mathbb{R}^n)$, phương trình bậc phân không - thời gian, và phương trình trừu tượng trong không gian Banach $X$.
  2. Phát triển hoàn thiện phương pháp chỉnh hóa làm nhuyễn Dirichlet đa chiều, chứng minh tính hội tụ bền vững và xác lập các công thức chọn tham số chỉnh hóa tối ưu theo cả nguyên lý tiên nghiệm và hậu nghiệm.
  3. Chứng minh tốc độ hội tụ đạt bậc tối ưu lý thuyết (Optimal Order) cho bài toán parabolic bậc phân, giải quyết triệt để hạn chế về tốc độ hội tụ Logarit trong các công bố quốc tế tiền nhiệm.
  4. Đột phá phương pháp luận với lý thuyết nửa nhóm giải tích, mở ra lời giải chỉnh hóa hoàn chỉnh cho bài toán ngược parabolic trong không gian Banach vô hạn chiều mà không phụ thuộc vào cấu trúc Hilbert.
  5. Thiết lập thuật toán giải số chính xác, hiệu năng cao, được kiểm chứng qua các ví dụ mô phỏng số với mức nhiễu thực tế, chứng minh tính khả thi tuyệt đối trong ứng dụng kỹ thuật và chuyển giao công nghệ.
  6. Công bố 03 công trình khoa học đỉnh cao trên các tạp chí quốc tế chuyên ngành SCIE uy tín (Inverse Problems in Science and Engineering, Applicable Analysis, Applied Numerical Mathematics), khẳng định vị thế tiên phong và đóng góp học thuật tầm quốc tế của luận án.