Tổng quan luận án
Bài toán ngược (inverse problems) và bài toán không chỉnh (ill-posed problems) đối với phương trình đạo hàm riêng là hướng nghiên cứu có vai trò quan trọng trong toán giải tích và toán ứng dụng. Trong các mô hình vật lý, hóa học, cơ học, động lực học chất lỏng và sinh học tế bào, hiện tượng khuếch tán dị thường (anomalous diffusion) và các quá trình truyền nhiệt với hiệu ứng nhớt hoặc môi trường không đồng nhất thường được mô tả bởi các phương trình vi phân cấp không nguyên (fractional differential equations). Khi mô hình hóa các quá trình này, hai bài toán nghịch đảo cơ bản nảy sinh là: xác định hàm nguồn (source identification problem) và xác định điều kiện ban đầu hay bài toán giá trị cuối (backward/final value problem).
Trong thực tế thực nghiệm, các đại lượng trạng thái tại thời điểm cuối hoặc điều kiện biên được đo đạc bằng các thiết bị vật lý chuyên dụng và luôn chứa đựng sai số đo đạc với mức nhiễu $\delta > 0$. Theo định nghĩa của Hadamard, các bài toán ngược này nhìn chung đều không chỉnh: nghiệm có thể 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 ban đầu. Một sai số rất nhỏ trong dữ liệu quan sát có thể dẫn đến sai số rất lớn ở nghiệm tìm kiếm, gây mất ổn định nghiêm trọng trong tính toán số. Do đó, việc xây dựng các phương pháp chỉnh hóa (regularization methods) nhằm khôi phục nghiệm ổn định và thiết lập tốc độ hội tụ của nghiệm chỉnh hóa về nghiệm chính xác là yêu cầu cấp thiết.
Luận án tiến sĩ toán học 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", thực hiện 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 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 các bài toán cụ thể sau:
- Bài toán 1: Tìm hàm nguồn phụ thuộc không gian $f(x)$ cho phương trình khuếch tán với đạo hàm Riemann-Liouville bậc $1 - \alpha \in (0,1)$:
$$\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 quan sát tại thời điểm cuối $u(x,T) = h(x)$ trong $\Omega$.
- Bài toán 2: Xác định trạng thái ban đầu $u(x,0)$ và phân tích bài toán giá trị cuối phi địa phương cho phương trình khuếch tán phi tuyến với đạo hàm Riemann-Liouville bậc $1 - \alpha \in (0,1)$:
$$\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 $u(x,t) = 0$ trên $\partial\Omega \times [0,T)$ và $u(x,T) = f(x)$ trong $\Omega$.
- Bài toán 3: Khôi phục hàm nguồn $f(x,y)$ cho phương trình khuếch tán với đạo hàm Riemann-Liouville bậc $1 - \alpha \in (0,1)$ trong miền hai chiều $\Omega = (0,\pi)^2$:
$$\partial_t u(x,y,t) = \partial_t^{1-\alpha} (u_{xx} + u_{yy}) + \beta(t)f(x,y), \quad (x,y,t) \in \Omega \times (0,T)$$
với $u(x,y,t) = 0$ trên $\partial\Omega \times (0,T]$ và dữ liệu cuối $u(x,y,T) = g(x,y)$.
- Bài toán 4: Xác định hàm nguồn $f(x)$ cho phương trình giả - parabol (pseudo-parabolic) với đạo hàm Caputo bậc $\alpha \in (0,1)$ và toán tử Laplacian phân số $(-\Delta)^s$ ($0 < s < 1$):
$$\partial_t^\alpha u + (-\Delta)^s u + \nu \partial_t^\alpha (-\Delta)^s u = \varphi(t)f(x), \quad (x,t) \in \Omega \times (0,T]$$
kết hợp dữ liệu tích phân thời gian phi địa phương.
- Bài toán 5: Chỉnh hóa bài toán giá trị cuối cho phương trình biharmonic với hàm nguồn tuyến tính và phi tuyến:
$$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 $u(x,t) = \Delta u(x,t) = 0$ trên $\partial\Omega \times (0,T)$ và dữ liệu quan sát đo đạc thuộc không gian $L^q(\Omega)$ ($q \neq 2$).
Đối tượng nghiên cứu của luận án là các bài toán xác định hàm nguồn và bài toán giá trị cuối đối với phương trình đạo hàm riêng cấp phân số và phương trình vi phân cấp cao. Phạm vi nghiên cứu thuộc lĩnh vực phương trình vi phân - đạo hàm riêng, lý thuyết toán tử elliptic, giải tích hàm trên các không gian Hilbert $L^2(\Omega)$, $H^\mu(\Omega)$, thang không gian Sobolev $H_0^1(\Omega) \cap H^2(\Omega)$ và không gian Lebesgue $L^q(\Omega)$ trên miền bị chặn $\Omega \subset \mathbb{R}^d$ ($d \ge 1$) với biên trơn $\partial\Omega$ trong khoảng thời gian hữu hạn $t \in [0,T]$.
Tổng quan tài liệu và vị trí của luận án
Nghiên cứu về bài toán ngược và bài toán không chỉnh đã hình thành và phát triển qua nhiều giai đoạn với các công trình nền tảng:
- John (1955) nghiên cứu bài toán Cauchy cho phương trình truyền nhiệt ngược thời gian và chứng minh tính ổn định trên tập các hàm số dương bị chặn.
- Tikhonov (1963) đề xuất phương pháp chỉnh hóa Tikhonov đặt nền móng cho khoa học tính toán bài toán ngược.
- Lattes và Lions (1967) đưa ra phương pháp tựa đảo (quasi-reversibility).
- Ewing (1979) sử dụng phương pháp phương trình Sobolev để chỉnh hóa bài toán thuần nhất.
- Showalter đề xuất phương pháp tựa giá trị biên (quasi-boundary value - QBV).
Đối với trường hợp phương trình truyền nhiệt và parabolic cổ điển ($\alpha = 1$):
- Trong, Dinh và Long (2005) khảo sát bài toán nhiệt ngược không thuần nhất $u_t - u_{xx} = -\varphi(t)f(x)$ dựa trên biến đổi Fourier.
- Trọng, Định và Nam (2009) nghiên cứu xác định cặp $(u,f)$ cho phương trình nhiệt hai chiều trên $(0,1)^2$ bằng phương pháp nội suy và chặt cụt chuỗi Fourier.
- Nam, Trọng và Tuấn (2010) giải bài toán nhiệt ngược hai chiều bằng phương pháp chặt cụt chuỗi Fourier, đạt ước lượng sai số cấp $(\ln(e/\delta))^{-p/2}$ trong $H^p(\Omega)$.
- Trọng và Tuấn (2011, 2014) khảo sát bài toán parabolic ngược phi tuyến với toán tử tự liên hợp dương và hàm nguồn Lipschitz địa phương bằng phương pháp chặt cụt Fourier và tựa phương trình chính.
- Tuấn, Trọng và Quân (2012) khảo sát bài toán nhiệt ngược phi tuyến hai chiều bằng phương pháp tựa giá trị biên phi địa phương.
Đối với trường hợp phương trình khuếch tán với đạo hàm cấp phân số ($\alpha \neq 1$):
- Ji Jun Liu và Yamamoto (2011) khảo sát bài toán ngược thời gian một chiều cho phương trình đạo hàm phân số bằng phương pháp tựa phương trình chính.
- Jun-Gang Wang, Yu-Bin Zhou và Ting Wei (2013), Ting Wei và Jun-Gang Wang (2014) đề xuất phương pháp giá trị tựa biên với quy tắc lựa chọn tham số tiên nghiệm và hậu nghiệm để xác định dữ liệu ban đầu và hàm nguồn phụ thuộc không gian.
- Ting Wei và cộng sự (2016) nghiên cứu hàm nguồn phụ thuộc thời gian trong phương trình khuếch tán phân đoạn từ dữ liệu Cauchy.
- Ting Wei và Jun Xian (2019) giải bài toán khuếch tán phân số bằng phương pháp biến phân kết hợp chỉnh hóa Tikhonov và thuật toán gradient liên hợp.
- Tuấn, Long, Thịnh và Thạnh Tran (2016) khôi phục điều kiện đầu cho phương trình khuếch tán cấp phân số bằng phương pháp tựa giá trị biên.
- Tuấn và cộng sự (2017) khảo sát bài toán khôi phục hàm nguồn với dữ liệu tại một thời điểm cố định bằng phương pháp chặt cụt tích phân và biến đổi Fourier.
- Trong, Hai và Dien (2019) chứng minh tính chỉnh khi $s > 0$ và tính không chỉnh khi $s = 0$ cho bài toán khuếch tán phân đoạn không-thời gian $\partial_t^\alpha u + \kappa_s L_s u = f(x,t)$.
- Điện, Hải, Việt và Trọng (2020) nghiên cứu bài toán ngược nguồn nhiệt phụ thuộc thời gian trong thang không gian Hilbert $H^\mu(\mathbb{R}^n)$ bằng phương pháp Tikhonov.
- Lực, Long, Dumitru Baleanu và Ravi P. Agarwal (2021) khảo sát bài toán xác định hàm nguồn với điều kiện tích phân phi địa phương $\int_0^T u(x,t)dt = \varphi(x)$, sử dụng phương pháp Fractional Landweber, đạt tốc độ hội tụ $\delta^{\frac{r}{r+1}}$ theo quy tắc tiên nghiệm và hậu nghiệm.
Khoảng trống nghiên cứu được luận án lựa chọn giải quyết là: hệ thống hóa và giải quyết trọn vẹn bài toán xác định hàm nguồn và bài toán giá trị cuối đối với phương trình khuếch tán cấp phân số sử dụng đạo hàm Riemann-Liouville, Caputo và toán tử Laplacian phân số $(-\Delta)^s$; thiết lập biểu thức nghiệm tích phân (mild solution); chứng minh tường minh tính không chỉnh theo nghĩa Hadamard; đề xuất các sơ đồ chỉnh hóa gồm phương pháp tựa giá trị biên (QBV), phương pháp chặt cụt Fourier, và phương pháp Fractional Landweber; chứng minh các đánh giá sai số hội tụ dưới cả hai quy tắc chọn tham số tiên nghiệm và hậu nghiệm; đồng thời mở rộng việc phân tích bài toán giá trị cuối cho phương trình biharmonic với hàm nguồn tuyến tính và phi tuyến trong không gian dữ liệu $L^q(\Omega)$ tổng quát ($q \neq 2$).
Cơ sở lý thuyết và phương pháp nghiên cứu
Luận án sử dụng khung phân tích giải tích hàm, lý thuyết phương trình đạo hàm riêng và lý thuyết bài toán không chỉnh theo nghĩa Hadamard. Một bài toán được định nghĩa là chỉnh theo nghĩa Hadamard nếu thỏa mãn đồng thời ba điều kiện:
- Nghiệm tồn tại;
- Nghiệm là duy nhất;
- Nghiệm phụ thuộc liên tục vào dữ liệu ban đầu.
Nếu vi phạm ít nhất một trong ba điều kiện trên, bài toán là không chỉnh.
Các công cụ toán học và khái niệm cơ sở được sử dụng gồm:
- Hàm Gamma:
$$\Gamma(z) = \int_0^\infty t^{z-1} e^{-t} dt, \quad z > 0$$
- Hàm Mittag-Leffler hai tham số: Với $\alpha > 0, \beta \in \mathbb{R}$,
$$E_{\alpha,\beta}(z) = \sum_{k=0}^\infty \frac{z^k}{\Gamma(k\alpha + \beta)}, \quad z \in \mathbb{C}$$
Các trường hợp đặc biệt: $E_{1,1}(z) = e^z$, $E_{2,1}(z) = \cosh\sqrt{z}$, $E_{2,2}(z) = \frac{\sinh\sqrt{z}}{\sqrt{z}}$.
- Đạo hàm cấp phân số Riemann-Liouville: Với $\alpha \in (0,1)$,
$$\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$$
- Đạo hàm cấp phân số Caputo: Với $\alpha \in (0,1)$,
$$\partial_t^\alpha u(x,t) = \frac{1}{\Gamma(1-\alpha)} \int_0^t (t-\tau)^{-\alpha} u_\tau(x,\tau) d\tau, \quad t > 0$$
- Toán tử Laplacian phân số $(-\Delta)^s$ ($0 < s < 1$):
$$(-\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 C_{d,s} = \frac{4^s \Gamma(d/2+s)}{\pi^{d/2} |\Gamma(-s)|}$$
- Toán tử Elliptic đều, đối xứng và xác định dương $L$: Xác định trên miền $D(L) = H_0^1(\Omega) \cap H^2(\Omega) \subset L^2(\Omega)$. Phổ của $L$ là dãy các giá trị riêng thực dương $0 < \lambda_1 \le \lambda_2 \le \dots \le \lambda_n \to \infty$, tương ứng với hệ hàm riêng trực chuẩn ${e_n}_{n=1}^\infty$ tạo thành cơ sở trực giao trong $L^2(\Omega)$.
- Thang không gian Hilbert $H^\mu(\Omega)$:
$$H^\mu(\Omega) = \left{ \varphi \in L^2(\Omega) : \sum_{n=1}^\infty \lambda_n^{2\mu} |(\varphi, e_n)|^2 < +\infty \right}, \quad |\varphi|{H^\mu(\Omega)} = \left( \sum{n=1}^\infty \lambda_n^{2\mu} |(\varphi, e_n)|^2 \right)^{1/2}$$
Phương pháp nghiên cứu thực tế bao gồm:
- Sử dụng phương pháp tách biến, biến đổi Laplace và biến đổi Fourier để tìm dạng nghiệm nhẹ (mild solution) của bài toán dưới dạng chuỗi trực giao.
- Phân tích sự suy giảm số mũ hoặc sự triệt tiêu của mẫu số tại giá trị riêng lớn $\lambda_n \to \infty$ để chỉ ra tính không chỉnh.
- Áp dụng 3 kỹ thuật chỉnh hóa:
- Phương pháp tựa giá trị biên (Quasi-Boundary Value Method);
- Phương pháp chặt cụt chuỗi Fourier (Fourier Truncation Method);
- Phương pháp Fractional Landweber (Landweber cấp phân số).
- Thiết lập quy tắc lựa chọn tham số chỉnh hóa tiên nghiệm ($a$-priori) và hậu nghiệm ($a$-posteriori) theo nguyên lý sai lệch (Morozov discrepancy principle).
- Nguồn dữ liệu kiểm chứng: Dữ liệu quan sát lý thuyết kết hợp với mô phỏng sai số nhân tạo mức $\delta > 0$ trong không gian $L^2(\Omega)$ và $L^q(\Omega)$ ($q \neq 2$).
Nội dung chính theo từng chương
Chương 1: Mở đầu
Chương 1 trình bày tính cấp thiết, bối cảnh khoa học, lý do chọn đề tài, đối tượng, phạm vi nghiên cứu cùng ý nghĩa khoa học và thực tiễn của luận án. Nội dung luận án được tổng hợp từ 5 công trình khoa học đã công bố trên các tạp chí quốc tế thuộc danh mục ISI (Web of Science):
- Mathematical Methods in the Applied Sciences (ISI, Q1, IF = 3);
- Computational and Applied Mathematics (ISI, Q2, IF = 3);
- Journal of Inequalities and Applications (ISI, Q2, IF = 4);
- Acta Mathematica Sinica, English Series (ISI, Q2, IF = 0);
- AIMS Mathematics (ISI, Q2, IF = 2).
Chương 2: Tổng quan
Chương 2 phân tích tình hình nghiên cứu trong nước và quốc tế đối với hai nhóm bài toán chính: bài toán xác định hàm nguồn (Loại I) và bài toán giá trị cuối (Loại II) cho phương trình khuếch tán cổ điển ($\alpha = 1$) và phân số ($\alpha \neq 1$). Chương này xác định rõ mục tiêu nghiên cứu và cấu trúc 5 bài toán cụ thể cần giải quyết trong luận án.
Chương 3: Phương pháp nghiên cứu
Chương 3 trình bày các phương pháp tiếp cận: phương pháp đọc tài liệu, phân loại và hệ thống hóa lý thuyết; thiết lập các định nghĩa về tính không chỉnh Hadamard, các tính chất vi tích phân của hàm Gamma, hàm Mittag-Leffler $E_{\alpha,\beta}(z)$, đạo hàm Riemann-Liouville, đạo hàm Caputo, hệ toán tử elliptic đối xứng xác định dương, hệ cơ sở trực chuẩn và định lý phép nhúng Sobolev.
Chương 4: Kết quả nghiên cứu và phân tích, đánh giá, thảo luận
Chương 4 là chương nội dung trọng tâm của luận án, được chia thành 5 phần tương ứng với các kết quả nghiên cứu:
4.1 Khôi phục hàm nguồn ngược cho phương trình khuếch tán với đạo hàm Riemann-Liouville
Xem xét phương trình $\partial_t u(x,t) = \partial_t^{1-\alpha} \Delta u(x,t) + f(x)q(t)$ với điều kiện biên Dirichlet thuần nhất và dữ liệu đo cuối $u(x,T) = h(x)$.
- Khai triển nghiệm chuỗi Fourier $u(x,t) = \sum_{n=1}^\infty u_n(t) e_n(x)$, áp dụng biến đổi Laplace chuyển phương trình về miền tần số:
$$\tau \widehat{u}_n(\tau) - u_n(0) = -\lambda_n \tau^{1-\alpha} \widehat{u}_n(\tau) + \widehat{F}_n(\tau)$$
- Xây dựng nghiệm nhẹ của hàm nguồn $f(x)$ và chỉ ra nguyên nhân gây không chỉnh: mẫu số chứa tích phân hàm Mittag-Leffler $\int_0^T E_{\alpha,1}(-\lambda_n (T-\tau)^\alpha) d\tau$ tiến về 0 khi $\lambda_n \to \infty$, dẫn đến việc sai số đo đạc $|h^\delta - h|_{L^2(\Omega)} \le \delta$ bị khuếch đại vô hạn.
- Đề xuất phương pháp chỉnh hóa tựa giá trị biên (QBV), xây dựng nghiệm xấp xỉ $f_\epsilon(x)$, chứng minh tính ổn định và thiết lập tốc độ hội tụ $|f_\epsilon - f|_{L^2(\Omega)} \to 0$ khi $\delta \to 0$ theo quy tắc chọn tham số tiên nghiệm.
4.2 Chỉnh hóa bài toán giá trị cuối cho phương trình khuếch tán với đạo hàm khác đạo hàm cổ điển
Xem xét phương trình $\partial_t u(x,t) = \partial_t^{1-\alpha} \Delta u(x,t) + F(x,t)$ trên $T := \Omega \times [0,T)$ với $u(x,T) = f(x)$.
- Chứng minh tính không chỉnh tại thời điểm ban đầu $t = 0$.
- Với hàm nguồn $F \in L^\infty(0,T; L^2(\Omega)) \cap X(0,T; L^2(\Omega))$, áp dụng phương pháp tựa giá trị biên để xây dựng nghiệm chỉnh hóa $u_{\epsilon(\delta)}(x,0)$.
- Chứng minh sự hội tụ của nghiệm chỉnh hóa về nghiệm chính xác dưới cả hai quy tắc chọn tham số: tiên nghiệm ($a$-priori) và hậu nghiệm ($a$-posteriori).
4.3 Khôi phục hàm nguồn cho phương trình khuếch tán trong miền hai chiều
Khảo sát hệ phương trình trên miền $\Omega = (0,\pi)^2$:
$$\partial_t u(x,y,t) = \partial_t^{1-\alpha} (u_{xx} + u_{yy}) + \beta(t)f(x,y), \quad u(x,y,T) = g(x,y)$$
- Áp dụng phương pháp chặt cụt chuỗi Fourier (Fourier Truncation Method) với ngưỡng cắt $N(\delta) \in \mathbb{N}^*$.
- Thiết lập công thức nghiệm chỉnh hóa $f_N(x,y) = \sum_{n,m=1}^N f_{n,m} \sin(nx)\sin(my)$.
- Chứng minh các đánh giá sai số hội tụ $|f_N - f|_{L^2(\Omega)}$ theo quy tắc chọn tham số tiên nghiệm và hậu nghiệm.
4.4 Chỉnh hóa bài toán xác định hàm nguồn cho phương trình giả - parabol với dữ liệu phi địa phương
Khảo sát phương trình giả - parabol chứa đạo hàm Caputo $\partial_t^\alpha$ và toán tử Laplacian phân số $(-\Delta)^s$ kết hợp số hạng nhớt $\nu \partial_t^\alpha (-\Delta)^s u$.
- Chứng minh tính không chỉnh của bài toán xác định hàm nguồn $f(x)$ với điều kiện đo phi địa phương theo thời gian.
- Đề xuất phương pháp Fractional Landweber (phương pháp lặp Landweber bậc không nguyên) để xây dựng nghiệm chỉnh hóa.
- Thiết lập đánh giá sai số và tốc độ hội tụ theo nguyên lý sai lệch tiên nghiệm và hậu nghiệm.
4.5 Chỉnh hóa bài toán giá trị cuối cho phương trình biharmonic với dữ liệu trong $L^q(\Omega)$
Khảo sát phương trình biharmonic:
$$Mu = u_{tttt} + 2u_{ttxx} + u_{xxxx} = G(x,t,u(x,t))$$
với điều kiện biên $u = \Delta u = 0$ trên $\partial\Omega$, $u(x,T) = f(x)$, $u_t(x,T) = 0$, $\Delta u(x,T) = g(x)$, $\Delta u_t(x,T) = 0$.
- Mở rộng không gian dữ liệu quan sát sang không gian Lebesgue $L^q(\Omega)$ tổng quát với $q \neq 2$, thỏa mãn mức nhiễu:
$$|f^\delta - f|{L^q(\Omega)} + |g^\delta - g|{L^q(\Omega)} + |G^\delta - G|_{L^\infty(0,T; L^q(\Omega))} < \delta$$
- Sử dụng phương pháp chặt cụt chuỗi Fourier để thiết lập bài toán chỉnh hóa cho cả hai trường hợp: hàm nguồn tuyến tính $G(x,t)$ và hàm nguồn phi tuyến $G(x,t,u)$.
- Đưa ra các đánh giá sai số hội tụ trong không gian nghiệm phù hợp.
Tổng hợp các bài toán nghiên cứu và phương pháp chỉnh hóa trong luận án
| STT |
Bài toán nghiên cứu |
Phương trình vi phân & Toán tử |
Phương pháp chỉnh hóa |
Không gian & Quy tắc tham số |
| 1 |
Xác định hàm nguồn $f(x)$ |
$\partial_t u = \partial_t^{1-\alpha}\Delta u + f(x)q(t)$ (Riemann-Liouville) |
Tựa giá trị biên (QBV) |
$L^2(\Omega)$, Tiên nghiệm |
| 2 |
Bài toán giá trị cuối $u(x,0)$ |
$\partial_t u = \partial_t^{1-\alpha}\Delta u + F(x,t)$ (Riemann-Liouville) |
Tựa giá trị biên (QBV) |
$L^2(\Omega)$, Tiên nghiệm & Hậu nghiệm |
| 3 |
Xác định hàm nguồn 2D $f(x,y)$ |
$\partial_t u = \partial_t^{1-\alpha}\Delta u + \beta(t)f(x,y)$ trên $(0,\pi)^2$ |
Chặt cụt Fourier |
$L^2(\Omega)$, Tiên nghiệm & Hậu nghiệm |
| 4 |
Xác định nguồn giả - parabol |
$\partial_t^\alpha u + (-\Delta)^s u + \nu \partial_t^\alpha (-\Delta)^s u = \varphi(t)f(x)$ |
Fractional Landweber |
$L^2(\Omega)$, Tiên nghiệm & Hậu nghiệm |
| 5 |
Giá trị cuối biharmonic |
$u_{tttt} + 2u_{ttxx} + u_{xxxx} = G(x,t,u)$ (Tuyến tính & Phi tuyến) |
Chặt cụt Fourier |
$L^q(\Omega)$ ($q \neq 2$), Tiên nghiệm & Hậu nghiệm |
Kết quả thử nghiệm số minh họa
Luận án trình bày hệ thống các bảng số liệu và đồ thị mô phỏng số (từ Bảng 4.1 đến Bảng 4.4 và từ Hình 4.1 đến Hình 4.17):
- Bảng 4.1: Sai số giữa nghiệm chính xác và nghiệm xấp xỉ tại các điểm $x \in {0.1, \dots}$.
- Bảng 4.2 và Bảng 4.3: Sai số giữa nghiệm chính xác và nghiệm chỉnh hóa tại $y = 0.3$, $y = 0.4$ với số bước lặp $N = 50, 100$.
- Bảng 4.4: Sai số tại mức nhiễu $\delta = 0.015$ dưới các giá trị $\alpha \in {0.2, 0.4, 0.5, 0.6, 0.8}$.
- Các hình vẽ (Hình 4.1 – 4.17) thể hiện bề mặt đồ thị nghiệm chính xác, nghiệm chỉnh hóa và sai số tuyệt đối dưới quy tắc chọn tham số tiên nghiệm và hậu nghiệm trong miền $D_\Omega \times D_T \in (0,1] \times (0,1)$.
Kết quả và những đóng góp mới
Các kết quả và đóng góp mới của luận án được công bố trên các tạp chí khoa học quốc tế bao gồm:
- Làm rõ bản chất toán học của tính không chỉnh: Luận án chỉ ra nguyên nhân gây ra sự không chỉnh của bài toán khôi phục hàm nguồn và bài toán giá trị cuối đối với phương trình khuếch tán cấp phân số thông qua việc phân tích phổ của toán tử và tính triệt tiêu của tích phân hàm Mittag-Leffler.
- Đề xuất và phát triển các phương pháp chỉnh hóa hiệu quả: Đề xuất phương pháp tựa giá trị biên (QBV), phương pháp chặt cụt Fourier, và phương pháp Fractional Landweber để giải quyết các bài toán ngược phân số 1D, 2D và bài toán chứa toán tử Laplacian phân số $(-\Delta)^s$.
- Thiết lập tốc độ hội tụ tối ưu: Thiết lập tường minh các đánh giá sai số giữa nghiệm chỉnh hóa và nghiệm chính xác dưới cả hai chiến lược: quy tắc chọn tham số tiên nghiệm ($a$-priori) và quy tắc chọn tham số hậu nghiệm ($a$-posteriori).
- Mở rộng bài toán giá trị cuối cho phương trình biharmonic: Khảo sát thành công bài toán giá trị cuối cho phương trình biharmonic trong cả hai trường hợp hàm nguồn tuyến tính và phi tuyến với dữ liệu quan sát đo đạc trong không gian Lebesgue $L^q(\Omega)$ ($q \neq 2$).
- Kiểm chứng bằng mô phỏng số: Xây dựng các thuật toán tính toán và cung cấp các ví dụ số minh họa, kiểm chứng tính đúng đắn của các định lý và ước lượng lý thuyết.
- Ứng dụng thực tiễn và đào tạo: Các kết quả nghiên cứu phục vụ trực tiếp công tác nghiên cứu khoa học, thảo luận seminar chuyên môn và ứng dụng trong đào tạo học viên cao học, nghiên cứu sinh tại Khoa Toán - Tin học, Trường Đại học Khoa học Tự nhiên, ĐHQG-HCM.
Hạn chế và hướng nghiên cứu tiếp
Dựa trên các nội dung chưa khảo sát trong phạm vi luận án, tác giả đề xuất 4 hướng nghiên cứu tiếp theo trong tương lai:
- Hướng 1: Khảo sát bài toán khôi phục hàm nguồn phụ thuộc theo biến thời gian $f(t)$ thay vì biến không gian.
- Hướng 2: Khảo sát bài toán giá trị cuối cho các mô hình khuếch tán sử dụng các dạng đạo hàm cấp phân số khác như đạo hàm Atangana-Baleanu-Caputo (ABC), đạo hàm Caputo-Fabrizio (CF) hoặc đạo hàm conformable.
- Hướng 3: Mở rộng khảo sát các mô hình khuếch tán cụ thể với hàm nguồn phi tuyến và dữ liệu quan sát đo đạc trong thang không gian $L^q$.
- Hướng 4: Khảo sát các bài toán phương trình đạo hàm riêng cấp phân số với dữ liệu chứa nhiễu ngẫu nhiên (random noise / white noise).
Giá trị tham khảo
Luận án là tài liệu tham khảo cho:
- Nghiên cứu sinh, học viên cao học và giảng viên ngành Toán giải tích, Toán ứng dụng, chuyên ngành Phương trình vi phân và tích phân, Bài toán ngược và không chỉnh.
- Các nhà nghiên cứu trong lĩnh vực vật lý toán và khoa học tính toán quan tâm đến kỹ thuật ước lượng hàm Mittag-Leffler, kỹ thuật chỉnh hóa Fourier, phương pháp tựa giá trị biên và sơ đồ lặp Landweber phân số.
- Phần phương pháp nghiên cứu (Chương 3) và các kỹ thuật đánh giá sai số tiên nghiệm/hậu nghiệm (Chương 4) có giá trị phương pháp luận cho việc xử lý các phương trình đạo hàm riêng cấp không nguyên.
Câu hỏi thường gặp
1. Vì sao các bài toán xác định hàm nguồn và bài toán giá trị cuối trong luận án là không chỉnh theo nghĩa Hadamard?
Trả lời: Trong biểu thức nghiệm dạng chuỗi Fourier, các hệ số nghiệm chứa mẫu số là các tích phân liên quan đến hàm Mittag-Leffler $E_{\alpha,1}(-\lambda_n t^\alpha)$. Khi chỉ số $n \to \infty$, giá trị riêng $\lambda_n \to \infty$ làm cho mẫu số tiến nhanh về 0. Khi dữ liệu đo đạc bị nhiễu bởi mức sai số $\delta > 0$, thành phần nhiễu ở tần số cao bị khuếch đại vô hạn, làm nghiệm không phụ thuộc liên tục vào dữ liệu ban đầu.
2. Luận án sử dụng những phương pháp chỉnh hóa chính nào?
Trả lời: Luận án sử dụng ba phương pháp chỉnh hóa chính:
- Phương pháp tựa giá trị biên (Quasi-Boundary Value Method);
- Phương pháp chặt cụt chuỗi Fourier (Fourier Truncation Method);
- Phương pháp Fractional Landweber (Landweber bậc không nguyên).
3. Sự khác biệt giữa quy tắc chọn tham số tiên nghiệm và hậu nghiệm trong luận án là gì?
Trả lời: Quy tắc chọn tham số tiên nghiệm ($a$-priori) xác định tham số chỉnh hóa $\epsilon$ hoặc ngưỡng cắt $N$ chỉ dựa trên mức độ nhiễu $\delta$ và độ trơn tiên nghiệm của nghiệm chính xác. Quy tắc chọn tham số hậu nghiệm ($a$-posteriori) xác định tham số chỉnh hóa dựa trực tiếp trên dữ liệu đo đạc thực tế thông qua nguyên lý sai lệch Morozov, giúp thuật toán thích ứng khi không biết trước độ trơn của nghiệm.
4. Điểm mở rộng trong nghiên cứu phương trình biharmonic ở phần 4.5 của luận án là gì?
Trả lời: Luận án đã mở rộng khảo sát bài toán giá trị cuối cho phương trình biharmonic $Mu = u_{tttt} + 2u_{ttxx} + u_{xxxx} = G(x,t,u)$ từ không gian chuẩn Hilbert $L^2(\Omega)$ thông thường sang không gian tổng quát $L^q(\Omega)$ ($q \neq 2$) cho cả hai trường hợp hàm nguồn tuyến tính và hàm nguồn phi tuyến.
5. Kết quả của luận án đã được công bố trên những ấn phẩm khoa học nào?
Trả lời: Luận án được tổng hợp từ 5 công trình công bố trên các tạp chí quốc tế thuộc danh mục ISI: Mathematical Methods in the Applied Sciences (Q1, IF = 3), Computational and Applied Mathematics (Q2, IF = 3), Journal of Inequalities and Applications (Q2, IF = 4), Acta Mathematica Sinica, English Series (Q2, IF = 0), và AIMS Mathematics (Q2, IF = 2).
Kết luận
Luận án tiến sĩ của tác giả Lê Đình Long đã nghiên cứu và giải quyết trọn vẹn các bài toán ngược và không chỉnh về xác định điều kiện ban đầu và khôi phục hàm nguồn cho các phương trình khuếch tán với đạo hàm cấp phân số Riemann-Liouville, Caputo, toán tử Laplacian phân số $(-\Delta)^s$ và phương trình biharmonic. Bằng việc phát triển các phương pháp tựa giá trị biên, chặt cụt chuỗi Fourier và Fractional Landweber, công trình đã thiết lập các ước lượng sai số hội tụ tối ưu dưới cả hai quy tắc chọn tham số tiên nghiệm và hậu nghiệm, đồng thời kiểm chứng tính chính xác qua các thử nghiệm số. Các kết quả của luận án đã được công bố trên 5 tạp chí quốc tế uy tín, đóng góp thiết thực cho sự phát triển của lý thuyết bài toán ngược và phương trình đạo hàm riêng cấp không nguyên.
