Tổng quan luận án

Luận án tiến sĩ "Bài toán giá trị cuối cho một số phương trình đạo hàm riêng" của nghiên cứu sinh Danh Hứa Quốc Nam được 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, thuộc ngành Toán giải tích (mã số: 9460102) dưới sự hướng dẫn khoa học của PGS. Nguyễn Huy Tuấn.

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

Trong giải tích toán học, phương trình đạo hàm riêng đóng vai trò quan trọng trong việc mô hình hóa các hiện tượng tự nhiên, vật lý, sinh học, môi trường và công nghệ xử lý ảnh. Các bài toán đạo hàm riêng được chia thành hai nhóm chính: bài toán thuận và bài toán ngược. Khác với bài toán thuận đã có hệ thống lý thuyết hoàn chỉnh, các bài toán ngược (ngược thời gian và xác định hàm nguồn) chỉ mới được tập trung nghiên cứu từ thập niên 60–70 của thế kỷ XX. Đa số các bài toán ngược đều rơi vào trạng thái không chỉnh theo nghĩa Hadamard, tức là nghiệm không thỏa mãn một trong ba điều kiện: tồn tại, duy nhất hoặc ổn định liên tục theo dữ liệu quan sát. Khi dữ liệu đầu vào bị nhiễu, sai số nhỏ có thể dẫn đến sai lệch lớn ở nghiệm đầu ra.

Khoảng trống nghiên cứu được luận án chỉ ra bao gồm:

  1. Các kết quả giải tích cho bài toán ngược thời gian đối với hệ phương trình parabolic chứa số hạng phi địa phương Kirchhoff còn hạn chế, cấu trúc phi tuyến dạng tích phân theo thời gian gây nhiều trở ngại trong tính toán và ước lượng.
  2. Các công trình khảo sát bài toán Cauchy cho phương trình elliptic đa phần chỉ dừng lại ở sai số quan sát trong không gian $L^2$, thiếu vắng các nghiên cứu khi dữ liệu bị nhiễu trong không gian tổng quát $L^p$ ($p \neq 2$).
  3. Mô hình bi-parabolic mô tả sự dẫn nhiệt hoàn chỉnh hơn phương trình parabolic cổ điển, nhưng việc khôi phục hàm nguồn bằng phương pháp chỉnh hóa Tikhonov cho phương trình này vẫn cần các đánh giá sai số lý thuyết tường minh.
  4. Nghiên cứu bài toán ngược thời gian cho phương trình sóng dầm mạnh trong điều kiện dữ liệu quan sát bị nhiễu ngẫu nhiên rời rạc (statistical discrete noise) còn rất ít được công bố.

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

Mục tiêu nghiên cứu của luận án được xác định rõ qua các nội dung:

  1. Xây dựng và thiết lập nghiệm chỉnh hóa cho 04 bài toán giá trị cuối/bài toán ngược đối với các phương trình đạo hàm riêng cụ thể.
  2. Chứng minh tính không chỉnh theo nghĩa Hadamard của các mô hình toán học tương ứng.
  3. Đánh giá sai số hội tụ giữa nghiệm chỉnh hóa đã đề xuất và nghiệm chính xác của bài toán dưới các điều kiện trơn tiên nghiệm khác nhau.
  4. Xây dựng thuật toán và minh họa một số kết quả mô phỏng số nhằm kiểm chứng cơ sở lý thuyết.

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

  • Đối tượng nghiên cứu: Nghiệm của 04 mô hình toán học cụ thể:
    • Bài toán 1: Bài toán ngược thời gian cho hệ phương trình parabolic có chứa số hạng Kirchhoff.
    • Bài toán 2: Bài toán giá trị cuối cho phương trình elliptic trong không gian $L^p$ ($p \neq 2$).
    • Bài toán 3: Bài toán xác định hàm nguồn cho phương trình bi-parabolic.
    • Bài toán 4: Bài toán giá trị cuối cho phương trình sóng dầm mạnh với dữ liệu bị nhiễu ngẫu nhiên rời rạc.
  • Phạm vi nghiên cứu: Lĩnh vực phương trình vi phân đạo hàm riêng, tập trung vào kỹ thuật chỉnh hóa bài toán ngược thời gian và bài toán xác định hàm nguồn trên miền bị chặn $\Omega \subset \mathbb{R}^N$ ($N \ge 1$) hoặc miền hình hộp chữ nhật $(0, \pi)^2 \times (0, M)$, trong khoảng thời gian xác định $t \in [0, T]$.

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

Các hướng nghiên cứu trước

Luận án đã tổng thuật hệ thống các công trình nền tảng và phương pháp chỉnh hóa:

  • Phương pháp tựa đảo (Quasi-Reversibility): Được đề xuất bởi Robert Lattès và Jacques-Louis Lions (1967) qua việc xấp xỉ toán tử $A$ bởi $A_\varepsilon = A - \varepsilon A^2$ cho phương trình vi phân trừu tượng $u_t + Au = 0$. Đến năm 1974, Ralph Edwin Showalter phát triển phương pháp này với toán tử $A_\varepsilon = A(I + \varepsilon A)^{-1}$ ($\varepsilon > 0$).
  • Phương pháp tựa giá trị biên (Quasi Boundary Value): Đề xuất bởi Ralph Edwin Showalter (1983) thông qua kỹ thuật làm nhiễu điều kiện biên không địa phương $u(T) + \varepsilon u(0) = g$, được Gordon Wayne Clark và Seth Fredric Oppenheimer (1994) khảo sát chi tiết. Mohamed Denche và Khaled Bessila đưa ra biến thể thay thế điều kiện cuối bằng $u(T) - \varepsilon u_t(0) = g$.
  • Phương pháp chỉnh hóa Tikhonov: Tổng thuật dựa trên công trình của Andreas Kirsch về phương pháp biến phân tìm phần tử cực tiểu hàm mục tiêu $J_\alpha(x) = |Kx - y|^2 + \alpha |x|^2$. Các công trình của Yang và cộng sự, Wei và cộng sự, Xiong và Fu đã vận dụng Tikhonov giải quyết các bài toán khuếch tán phân số và phương trình Helmholtz.
  • Phương pháp chặt cụt chuỗi Fourier: Nguyễn Huy Tuấn và các cộng sự đã ứng dụng thành công kỹ thuật loại bỏ tần số cao cho các lớp bài toán parabolic và elliptic phi tuyến.
  • Nghiên cứu về mô hình Kirchhoff: Bắt nguồn từ mô hình dao động sợi dây của Gustav Kirchhoff (1883), mở rộng từ phương trình d'Alembert (1747, 1749) và các bài toán phương trình phi địa phương của Lê Minh Triết cùng các nhà toán học khác.
  • Nghiên cứu phương trình elliptic và sóng dầm: Các nghiên cứu của Hongwu Zhang, Xiaoju Zhang, Renhu Wang, Võ Anh Khoa (hệ elliptic sine-Gordon), Nguyễn Hữu Cần, Erkan Nane, Mokhtar Kirane, Nguyễn Đăng Minh về hệ vi phân tắt dần mạnh, phương trình điện báo và phương trình Chen-Russell.
Tác giả / Năm Phương pháp / Mô hình Đối tượng khảo sát
Robert Lattès, Jacques-Louis Lions (1967) Phương pháp tựa đảo ($A - \varepsilon A^2$) Bài toán ngược thời gian trừu tượng
Ralph Edwin Showalter (1974, 1983) Tựa đảo ($A(I+\varepsilon A)^{-1}$) và Tựa giá trị biên Bài toán ngược thời gian parabolic
Andreas Kirsch Chỉnh hóa Tikhonov Bài toán ngược toán tử compact
Gustav Kirchhoff (1883) Mô hình phi địa phương Phương trình sóng phi tuyến
Hongwu Zhang, Renhu Wang Chỉnh hóa Lavrentiev và Tikhonov Phương trình elliptic nửa tuyến tính
Nguyễn Huy Tuấn và cs Chặt cụt Fourier, phương trình tích phân Bài toán Cauchy phi tuyến, phân số

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

Luận án giải quyết trực tiếp 4 khoảng trống:

  1. Xây dựng nghiệm chỉnh hóa cho hệ parabolic chứa số hạng Kirchhoff phi địa phương phụ thuộc $|\nabla u|{L^2}$ và $|\nabla v|{L^2}$ vốn dẫn đến nghiệm biểu diễn dưới dạng số mũ chứa tích phân phi tuyến.
  2. Vượt qua rào cản kỹ thuật của không gian Hilbert $L^2$ thông thường để thiết lập đánh giá sai số cho bài toán elliptic với sai số dữ liệu trong không gian $L^p$ ($p \neq 2$) bằng phép nhúng Hilbert scale.
  3. Giải quyết bài toán ngược xác định hàm nguồn $f(x)$ cho phương trình bi-parabolic cấp 4 theo không gian bằng phương pháp Tikhonov.
  4. Thiết lập thuật toán chỉnh hóa và đánh giá sai số ngẫu nhiên cho phương trình sóng dầm mạnh với dữ liệu rời rạc bị nhiễu thống kê trong không gian Bochner.

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

Khung lý thuyết và khái niệm

Luận án sử dụng định nghĩa bài toán không chỉnh theo Hadamard cho phương trình toán tử $Ku = g$: bài toán gọi là không chỉnh nếu vi phạm một trong ba điều kiện: (1) Nghiệm $u \in X$ tồn tại với mọi $g \in Y$; (2) Nghiệm $u$ là duy nhất; (3) Nghiệm $u$ ổn định liên tục theo $g$ (tức $|g - g_\epsilon|Y \to 0 \implies |u - u\epsilon|_X \to 0$).

Hệ thống không gian hàm được sử dụng:

  • Không gian Lebesgue $L^2(\Omega)$ với tích vô hướng $\langle\cdot,\cdot\rangle$ và chuẩn cảm sinh.
  • Không gian các hàm liên tục $C([0,T]; X)$, không gian $C^m[0,T]; X$, và không gian có trọng số $C^\alpha_\omega((0,T]; X)$ với chuẩn: $$|v|{C^\alpha\omega} := \sup_{0 < t \le T} t^\alpha |v(t)|_X < \infty$$
  • Không gian tích phân $L^p(0,T; X)$ với chuẩn $|v|_{L^p(0,T;X)} = \left(\int_0^T |v(t)|_X^p dt\right)^{1/p}$ ($1 \le p < \infty$).
  • Hệ cơ sở trực chuẩn ${\varphi_p}{p=1}^\infty$ của toán tử $-\Delta$ trên miền bị chặn $\Omega$ thỏa mãn $-\Delta \varphi_p = \lambda_p \varphi_p$ trong $\Omega$ và $\varphi_p = 0$ trên $\partial\Omega$, với $0 < \lambda_1 \le \lambda_2 \le \dots \le \lim{p\to\infty} \lambda_p = +\infty$.
  • Không gian Hilbert scale $H^s(\Omega)$ ($s \ge 0$): $$H^s(\Omega) = \left{f \in L^2(\Omega) : \sum_{p=1}^\infty \lambda_p^s \langle f, \varphi_p \rangle^2 < \infty\right}, \quad |f|{H^s(\Omega)}^2 = \sum{p=1}^\infty \lambda_p^s \langle f, \varphi_p \rangle^2$$ với tính chất nhúng $H^{s'}(\Omega) \hookrightarrow H^s(\Omega) \hookrightarrow L^2(\Omega)$ ($s' \ge s \ge 0$).
  • Lớp hàm Gevrey $G_{\sigma, \nu}$ ($\sigma, \nu > 0$): $$G_{\sigma, \nu} = \left{f \in L^2(\Omega) : |f|{G{\sigma,\nu}}^2 = \sum_{p=1}^\infty e^{\sigma \lambda_p} \lambda_p^\nu \langle f, \varphi_p \rangle^2 < +\infty\right}$$
  • Không gian Bochner $L^p_B \equiv L^p(\tilde{\Omega}, \mathcal{F}, \mu; B)$ trên không gian xác suất $(\tilde{\Omega}, \mathcal{F}, \mu)$ với chuẩn $|u|_{L^p_B} = \left(\mathbb{E}|u|_B^p\right)^{1/p}$.

Các định lý phụ trợ nòng cốt gồm Định lý điểm bất động Banach trên không gian metric đầy đủ và Bất đẳng thức Grönwall dạng tích phân.

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

  • Phương pháp tra cứu và kế thừa: Hệ thống hóa các kết quả nghiên cứu trong và ngoài nước về phương trình đạo hàm riêng và lý thuyết bài toán ngược.
  • Phương pháp phân loại và tổng hợp lý thuyết: Xây dựng biểu diễn nghiệm tích phân (nghiệm nhẹ) qua chuỗi Fourier trực chuẩn, sau đó phân tích các thành phần gây mất ổn định tại tần số cao.
  • Kỹ thuật chỉnh hóa:
    • Chặt cụt chuỗi Fourier: Giới hạn số hạng điều hòa $n \le N(\delta)$ để loại bỏ các số mũ dương phân kỳ.
    • Chỉnh hóa Tikhonov: Tối thiểu hóa phiếm hàm mục tiêu nhằm xác định nghiệm chỉnh hóa ổn định $x_\alpha = (\alpha I + K^K)^{-1} K^ y$.
    • Ước lượng phi tham số: Xử lý dữ liệu nhiễu ngẫu nhiên thống kê rời rạc.
  • Kỹ thuật đánh giá sai số: Sử dụng phép nhúng giữa không gian $L^p$ và Hilbert scale, kết hợp nguyên lý ánh xạ co để kiểm soát các toán tử phi tuyến.

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

Chương 1: Mở đầu

Chương này giới thiệu tổng quan về đề tài, nêu bật tính cấp thiết của việc nghiên cứu các phương trình đạo hàm riêng xuất phát từ thực tiễn mô phỏng khoa học kỹ thuật. Tác giả trình bày lý do lựa chọn đề tài, xác định rõ 04 đối tượng nghiên cứu, phạm vi giải tích toán học và phương pháp tiếp cận. Chương 1 cũng liệt kê danh mục 04 công trình khoa học của nghiên cứu sinh đã công bố trên các tạp chí quốc tế:

  • [P1] Mathematical Methods in the Applied Sciences (ISI-Q1)
  • [P2] Applicable Analysis (ISI-Q2)
  • [P3] Discrete and Continuous Dynamical Systems - Series S (ISI-Q2)
  • [P4] Advances in Difference Equations (tên mới: Advances in Continuous and Discrete Models, ISI-Q2)

Đồng thời, chương này tổng kết các báo cáo chuyên môn tại Đại hội Toán học toàn quốc lần thứ IX (Nha Trang, 2018), Hội nghị khoa học lần thứ XII Trường ĐHKHTN (2020), Hội nghị Toán học Miền Trung - Tây Nguyên lần thứ IV (Huế, 2022) và Seminar Bộ môn Giải tích (2023).

Chương 2: Tổng quan

Chương 2 trình bày chi tiết lịch sử và các phương pháp chỉnh hóa bài toán ngược tuyến tính và phi tuyến (Tựa đảo, Tựa giá trị biên, Chặt cụt chuỗi Fourier, Tikhonov). Chương này dành một phần quan trọng để khảo sát ứng dụng của bài toán ngược thời gian trong xử lý và khôi phục ảnh số bị nhiễu Gaussian thông qua biến đổi Fourier rời rạc (DFT).

Tác giả thực nghiệm trên dữ liệu ảnh kích thước $400 \times 256$ pixel cho logo hai trường đại học:

  • Logo Trường Đại học Thủ Dầu Một: Dữ liệu gốc bị làm mờ bởi nhiễu Gaussian được khôi phục với các tham số chỉnh hóa $\alpha \in {1 \times 10^{-1}, 2 \times 10^{-1}, 4 \times 10^{-1}}$, thời gian CPU chạy tương ứng là 10.838 giây, 16.532 giây và 20.370 giây.
  • Logo Trường Đại học Khoa học Tự nhiên - ĐHQG TP.HCM: Kết quả khôi phục tương ứng với thời gian chạy CPU là 8.428 giây ($\alpha = 0.1$), 9.023 giây ($\alpha = 0.2$) và 13.435 giây ($\alpha = 0.4$).

Chương 2 khép lại bằng việc tổng quan tình hình nghiên cứu chi tiết cho 04 bài toán mục tiêu của luận án.

Chương 3: Phương pháp nghiên cứu

Chương này định nghĩa toán học chính xác về tính không chỉnh theo nghĩa Hadamard, thiết lập cấu trúc các không gian hàm nền tảng ($L^2$, $L^p(0,T;X)$, $C^m[0,T;X]$, $H^s(\Omega)$, lớp Gevrey $G_{\sigma, \nu}$, Bochner $L^p_B$). Luận án trình bày hai định lý bổ trợ nòng cốt là Định lý điểm bất động Banach và Bất đẳng thức Grönwall. Quy trình nghiên cứu được vạch ra theo 2 bước: thiết lập nghiệm nhẹ và chứng minh tính không chỉnh; tiếp đó xây dựng nghiệm chỉnh hóa và đánh giá sai số hội tụ dựa trên các giả thiết trơn tiên nghiệm.

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à phần nội dung cốt lõi của luận án, được chia thành 4 phân mục tương ứng với 4 bài báo công bố:

4.1 Chỉnh hóa nghiệm bài toán ngược thời gian cho hệ phương trình parabolic có chứa số hạng Kirchhoff

Khảo sát hệ phương trình trên $\Omega \times (0,T)$: $$\begin{cases} \frac{\partial u}{\partial t} = L\left(|\nabla u|{L^2(\Omega)}, |\nabla v|{L^2(\Omega)}\right) \Delta u \ \frac{\partial v}{\partial t} = L\left(|\nabla u|{L^2(\Omega)}, |\nabla v|{L^2(\Omega)}\right) \Delta v \end{cases}$$ với điều kiện Dirichlet thuần nhất và dữ liệu cuối $u(x,T) = f(x)$, $v(x,T) = g(x)$. Hàm $L$ thỏa mãn $0 < M_0 \le L(z_1, z_2) \le M_1$ và tính chất liên tục Lipschitz với hằng số $K_l > 0$.

  • Nghiệm nhẹ: Được thiết lập dưới dạng chuỗi: $$u(x,t) = \sum_{n=1}^\infty \exp\left(\lambda_n \int_t^T L(|\nabla u(\cdot,s)|{L^2}, |\nabla v(\cdot,s)|{L^2}) ds\right) \langle f, e_n \rangle e_n(x)$$
  • Tính không chỉnh: Tác giả chọn dãy dữ liệu thử nghiệm $f_m(x) = g_m(x) = \frac{e_m(x)}{\lambda_m}$. Khi $m \to +\infty$, sai số dữ liệu $|(f_m, g_m)|{L^2 \times L^2} \to 0$, nhưng chuẩn nghiệm lại phân kỳ: $$|(u_m, v_m)|{L^\infty(0,T; H^1(\Omega))}^2 \ge \frac{2 \exp(2\lambda_m M_0 T)}{\lambda_m} \to +\infty$$
  • Chỉnh hóa và đánh giá sai số: Sử dụng phương pháp chặt cụt Fourier tại $N(\delta)$, nghiệm chỉnh hóa $(u_{N,\delta}, v_{N,\delta})$ được chứng minh tồn tại duy nhất qua nguyên lý ánh xạ co. Khi chọn $N(\delta)$ thỏa mãn $\frac{1-\nu}{2TM_1} \ln(1/\delta) \le \lambda_N \le \frac{1-\nu}{TM_1} \ln(1/\delta)$ với $0 < \nu < 1$, sai số hội tụ đạt bậc logarit: $$|u_{N,\delta} - u|{H^1}^2 + |v{N,\delta} - v|_{H^1}^2 \le \mathcal{O}\left(\max\left{\left(\ln\frac{1}{\delta}\right)^{-2\gamma}, \delta^{2\nu} \ln\frac{1}{\delta}\right}\right)$$

4.2 Chỉnh hóa nghiệm bài toán giá trị cuối cho phương trình elliptic trong không gian $L^p$

Khảo sát phương trình elliptic $\Delta u = u_{zz} + u_{xx} + u_{yy} = G(x,y,z,u)$ trên $(0,\pi)^2 \times (0,M)$ với điều kiện biên thuần nhất, $u(x,y,M) = f(x,y)$ và $u_z(x,y,M) = 0$.

  • Luận án xây dựng không gian Hilbert scale trên miền hai chiều và kỹ thuật nhúng Sobolev vào không gian $L^p$ ($p \neq 2$).
  • Khảo sát hoàn chỉnh cả trường hợp nguồn tuyến tính $G(x,y,z)$ (phương trình Poisson, Helmholtz) và phi tuyến $G(x,y,z,u)$ (phương trình sine-Gordon).
  • Đánh giá sai số hội tụ được xác lập trong chuẩn $L^p$ nhờ phương pháp chặt cụt chuỗi Fourier.

4.3 Chỉnh hóa nghiệm bài toán xác định hàm nguồn cho phương trình bi-parabolic

Khảo sát phương trình dẫn nhiệt cấp 4: $$u_{tt}(x,t) + 2\Delta u_t(x,t) + \Delta^2 u(x,t) = \psi(t)f(x), \quad (x,t) \in \Omega \times (0,T)$$ với điều kiện ban đầu thuần nhất $u(x,0) = u_t(x,0) = 0$ và dữ liệu đo cuối $u(x,T) = h(x)$.

  • Dữ liệu quan sát thực tế $(h_\varepsilon, \psi_\varepsilon)$ có sai số $|\psi_\varepsilon - \psi|{L^\infty(0,T)} \le \varepsilon$ và $|h\varepsilon - h|_{L^2(\Omega)} \le \varepsilon$.
  • Áp dụng phương pháp Tikhonov để tìm nghiệm chỉnh hóa ổn định $f^{\varepsilon,\alpha}$.
  • Tác giả lập bảng số liệu kiểm chứng sai số tương đối và tuyệt đối tại các tham số $M = 40, T = 1, \tau = 1.2$ (Bảng 4.1 trong luận án) và đồ thị so sánh hàm nguồn trước và sau chỉnh hóa (Hình 4.4).

4.4 Chỉnh hóa nghiệm bài toán giá trị cuối cho phương trình sóng dầm mạnh với dữ liệu bị nhiễu ngẫu nhiên rời rạc

Khảo sát phương trình: $$u_{tt} - \Delta u + b\Delta^2 u - a\Delta u_t = F(x,t,u(x,t)), \quad (x,t) \in \Omega \times [0,T]$$ với $u(x,T) = g(x)$, $u_t(x,T) = h(x)$.

  • Thiết lập khái niệm nghiệm nhẹ và xấp xỉ các hệ số Fourier từ tập dữ liệu rời rạc bị nhiễu ngẫu nhiên.
  • Sử dụng phương pháp chặt cụt chuỗi Fourier kết hợp ước lượng phi tham số thống kê trong không gian Bochner.
  • Minh họa đồ thị nghiệm chính xác $u_{ex}$ và nghiệm chỉnh hóa $u_{reg}$ tại thời điểm $t = 0.8$ với số điểm đo $n \in {50, 100, 150}$ (Hình 4.5) và bảng sai số định lượng (Bảng 4.2).

Chương 5: Kết luận và kiến nghị

Chương cuối tóm tắt toàn bộ các kết quả lý thuyết về tính không chỉnh, xây dựng nghiệm chỉnh hóa và đánh giá sai số cho 04 bài toán. Đồng thời đề xuất các hướng mở rộng nghiên cứu trong tương lai.

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

Đóng góp mới về mặt khoa học và lý luận

  1. Chứng minh tính không chỉnh: Chỉ rõ sự vi phạm điều kiện phụ thuộc liên tục của nghiệm vào dữ liệu ban đầu theo nghĩa Hadamard cho hệ parabolic Kirchhoff, phương trình elliptic, bi-parabolic và sóng dầm mạnh.
  2. Đề xuất phương pháp chỉnh hóa phù hợp:
    • Ứng dụng thành công phương pháp chặt cụt chuỗi Fourier cho Bài toán 1, Bài toán 2 và Bài toán 4.
    • Ứng dụng phương pháp chỉnh hóa Tikhonov để giải quyết bài toán xác định hàm nguồn cho phương trình bi-parabolic ở Bài toán 3.
  3. Thiết lập đánh giá sai số tường minh:
    • Chứng minh tốc độ hội tụ đạt bậc logarit đối với hệ parabolic Kirchhoff.
    • Thiết lập đánh giá sai số hội tụ trong chuẩn $L^p$ ($p \neq 2$) cho phương trình elliptic thông qua phép nhúng Hilbert scale.
    • Xác lập tính ổn định và hội tụ của nghiệm hàm nguồn theo chuẩn $L^2$ cho phương trình bi-parabolic.
    • Ước lượng sai số kỳ vọng toán học trong không gian Bochner cho phương trình sóng dầm mạnh với dữ liệu quan sát ngẫu nhiên rời rạc.

Đóng góp về mặt thực tiễn và tính toán số

  • Xây dựng thuật toán và thực hiện mô phỏng số kiểm chứng cho các bài toán lý thuyết: mô hình khôi phục ảnh số bị nhiễu Gaussian qua biến đổi Fourier rời rạc, bài toán xác định hàm nguồn bi-parabolic, bài toán sóng dầm mạnh ngẫu nhiên với cỡ mẫu $n \in {50, 100, 150}$.
  • Các kết quả nghiên cứu đã được công bố trên 04 bài báo khoa học thuộc danh mục ISI (01 bài Q1 và 03 bài Q2).

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

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

  • Các đánh giá sai số chủ yếu dựa trên các giả định trơn tiên nghiệm của nghiệm chính xác (chẳng hạn nghiệm thuộc lớp Gevrey hoặc không gian Sobolev chỉ số cao).
  • Mô hình hệ Kirchhoff mới dừng lại ở trường hợp hệ thuần nhất hoặc nguồn đơn giản, chưa xử lý triệt để trường hợp hàm nguồn phi tuyến tổng quát.
  • Các mô phỏng số cho bài toán elliptic và sóng dầm chủ yếu thực hiện trên miền hình học có dạng chữ nhật hoặc đoạn một chiều $(0, \pi)$.

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

Luận án chỉ rõ 03 hướng nghiên cứu dự kiến tiếp tục phát triển:

  1. Nghiên cứu bài toán ngược cho mô hình Kirchhoff với hàm nguồn phi tuyến.
  2. Khảo sát bài toán ngược cho một số phương trình sóng dầm trong không gian $L^p$.
  3. Mở rộng các bài toán ngược có liên quan đến quá trình Wiener hoặc chuyển động Brown (Brownian motion).

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên khảo hữu ích 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: Cung cấp kỹ thuật giải tích hàm hiện đại trong việc xử lý các bài toán không chỉnh, kỹ thuật nhúng Hilbert scale vào không gian $L^p$, và phương pháp ánh xạ co trên không gian hàm có trọng số.
  • Các nhà nghiên cứu trong lĩnh vực Vật lý toán và Xử lý tín hiệu - hình ảnh: Cung cấp cơ sở toán học tường minh cho bài toán khôi phục tín hiệu, xử lý ảnh mờ/nhiễu và mô hình hóa quá trình dẫn nhiệt, truyền sóng tắt dần.
  • Nội dung đáng tham khảo nhất:
    • Kỹ thuật ước lượng sai số logarit cho bài toán parabolic phi địa phương Kirchhoff tại Mục 4.1.
    • Khung lý thuyết nhúng $L^p$ cho phương trình elliptic tại Mục 4.2.
    • Thuật toán Tikhonov xác định hàm nguồn bi-parabolic tại Mục 4.3.
    • Kỹ thuật kết hợp ước lượng thống kê phi tham số với khai triển Fourier cho dữ liệu nhiễu ngẫu nhiên rời rạc tại Mục 4.4.

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

1. Tại sao bài toán ngược thời gian cho hệ phương trình Kirchhoff lại phức tạp hơn bài toán ngược cho phương trình nhiệt thông thường?

Trong mô hình Kirchhoff, hệ số khuếch tán phụ thuộc vào chuẩn gradient $|\nabla u|{L^2(\Omega)}$ và $|\nabla v|{L^2(\Omega)}$. Mặc dù số hạng này chỉ phụ thuộc vào thời gian $t$, nhưng khi giải ngược thời gian, biểu thức nghiệm chứa số mũ tích phân $\exp\left(\lambda_n \int_t^T L(|\nabla u(\cdot,s)|, |\nabla v(\cdot,s)|) ds\right)$. Cấu trúc phi tuyến phi địa phương này khiến toán tử nghiệm phụ thuộc phi tuyến phức tạp vào chính hàm cần tìm ở cả tử số và mẫu số, đòi hỏi phải sử dụng định lý điểm bất động Banach trên không gian $L^\infty(0,T; H^1(\Omega))$ với các đánh giá tích phân quy nạp phức tạp.

2. Phương pháp chặt cụt chuỗi Fourier giải quyết tính không chỉnh của Bài toán 1 như thế nào?

Tính không chỉnh của bài toán xuất phát từ việc thừa số $\exp(\lambda_n M_0 (T-t))$ tăng nhanh ra vô cùng khi số sóng $n \to \infty$. Phương pháp chặt cụt chuỗi Fourier giải quyết vấn đề này bằng cách chỉ giữ lại các tần số thấp $n \le N(\delta)$ và loại bỏ toàn bộ các thành phần dao động cao tần. Bằng cách chọn ngưỡng chặt cụt $N(\delta)$ phụ thuộc vào mức nhiễu $\delta$ theo công thức logarit $\lambda_N \sim \frac{1-\nu}{TM_1} \ln(1/\delta)$, nghiệm xấp xỉ đạt được sự cân bằng giữa sai số dữ liệu đo và sai số chặt cụt, đảm bảo nghiệm chỉnh hóa hội tụ về nghiệm chính xác với bậc logarit khi $\delta \to 0$.

3. Đóng góp đặc thù của luận án khi khảo sát phương trình elliptic trong không gian $L^p$ ($p \neq 2$) là gì?

Hầu hết các nghiên cứu trước đây về phương trình elliptic chỉ xét dữ liệu nhiễu trong không gian Hilbert $L^2$. Khi chuyển sang không gian Banach $L^p$ với $p \neq 2$, cấu trúc tích vô hướng trực chuẩn không còn áp dụng trực tiếp được. Luận án đã giải quyết khó khăn này bằng cách xây dựng không gian Hilbert scale trên miền chữ nhật $(0,\pi)^2$ và vận dụng các bất đẳng thức nhúng giữa $L^p$ và Hilbert scale, từ đó thiết lập được các đánh giá sai số hội tụ cho cả trường hợp tuyến tính (Poisson, Helmholtz) và phi tuyến (sine-Gordon).

4. Phương pháp chỉnh hóa Tikhonov trong Bài toán 3 được áp dụng cho mục tiêu gì và xử lý dữ liệu ra sao?

Phương pháp Tikhonov được áp dụng để xác định hàm nguồn không gian $f(x)$ trong phương trình bi-parabolic $u_{tt} + 2\Delta u_t + \Delta^2 u = \psi(t)f(x)$ khi biết dữ liệu cuối $u(x,T) = h(x)$. Luận án thiết lập toán tử chuyển tiếp tuyến tính compact từ $f$ sang $h$ và xây dựng nghiệm xấp xỉ qua việc cực tiểu hóa phiếm hàm Tikhonov $J_\alpha(f) = |Kf - h_\varepsilon|^2 + \alpha |f|^2$, giúp nghiệm $f^{\varepsilon,\alpha}$ ổn định và hội tụ về $f$ khi sai số dữ liệu $(\varepsilon)$ tiến về 0.

5. Xử lý nhiễu trong Bài toán 4 có điểm gì khác biệt so với các bài toán ngược truyền thống?

Trong Bài toán 4, dữ liệu đầu vào không phải là hàm liên tục với sai số xác định kiểu $|g - g_\delta| \le \delta$, mà là các giá trị đo rời rạc bị tác động bởi nhiễu ngẫu nhiên thống kê. Luận án đã kết hợp phương pháp chặt cụt chuỗi Fourier với kỹ thuật ước lượng phi tham số thống kê, phân tích sai số trên không gian Bochner $L^p_B$, và đánh giá sự hội tụ của nghiệm chỉnh hóa thông qua kỳ vọng toán học của bình phương sai số chuẩn.

Kết luận

Luận án tiến sĩ của tác giả Danh Hứa Quốc Nam là một công trình giải tích toán học nghiên cứu chuyên sâu về bài toán giá trị cuối và bài toán ngược cho 04 mô hình phương trình đạo hàm riêng (hệ parabolic Kirchhoff, elliptic trong $L^p$, bi-parabolic và sóng dầm mạnh). Bằng việc sử dụng linh hoạt các phương pháp chỉnh hóa chặt cụt Fourier, phương pháp Tikhonov kết hợp với các công cụ giải tích hàm và ước lượng thống kê, tác giả đã chứng minh tính không chỉnh, xây dựng nghiệm chỉnh hóa ổn định và thiết lập các đánh giá sai số hội tụ tường minh. Toàn bộ các kết quả lý thuyết chính đều được minh chứng qua các ví dụ mô phỏng số và được công bố trên 04 tạp chí quốc tế uy tín thuộc danh mục ISI.