Tổng quan luận án
Luận án tiến sĩ toán học của nghiên cứu sinh Đặng Quang Long, chuyên ngành Toán ứng dụng (mã số: 9 46 01 12), với đề tài "Sự tồn tại, duy nhất nghiệm và phương pháp lặp giải một số bài toán biên cho phương trình vi phân phi tuyến cấp ba", được thực hiện dưới sự hướng dẫn khoa học của GS. Nguyễn Đông Anh tại Học viện Khoa học và Công nghệ thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam vào năm 2023.
Nhiều bài toán trong các lĩnh vực cơ học, vật lý, sinh học, môi trường, kỹ thuật như uốn dầm, dẫn nhiệt, dòng nước ngầm, đàn hồi nhiệt và vật lý plasma dẫn đến việc giải các bài toán biên (Boundary Value Problems - BVPs) cho phương trình vi phân thường phi tuyến (ODE) cấp cao, phương trình vi tích phân (IDE) và phương trình vi phân hàm (FDE). Việc tìm nghiệm chính xác giải tích cho các bài toán này chỉ thực hiện được trong một số rất ít trường hợp đặc biệt; phần lớn các trường hợp đòi hỏi phải tìm nghiệm xấp xỉ bằng các phương pháp số.
Khoảng trống nghiên cứu được tác giả xác định xuất phát từ những hạn chế chính của các công bố trước đó:
- Đối với bài toán biên phương trình vi phân phi tuyến cấp ba và cấp bốn, các công trình trước thường áp dụng các điều kiện tồn tại nghiệm rất ngặt và khó kiểm chứng trên thực tế (như điều kiện kiểu Nagumo trên toàn không gian, điều kiện tăng tuyến tính ở vô cùng, điều kiện đơn điệu tăng hoặc điều kiện Lipschitz một phía trên toàn bộ $\mathbb{R}$ hoặc $\mathbb{R}^2$). Đa phần các nghiên cứu chỉ dừng lại ở việc chứng minh sự tồn tại thuần túy lý thuyết mà không đưa ra nghiệm giải tích cụ thể hoặc phương pháp tính toán tương ứng.
- Đối với các phương pháp số giải bài toán biên phi tuyến cấp ba, các tác giả trước đây (như các phương pháp dùng spline hoặc sai phân) chủ yếu tập trung thiết lập sơ đồ rời rạc và giả định rằng hệ đại số phi tuyến nhận được sẽ được giải bằng một phương pháp lặp đã biết, bỏ qua việc phân tích sai số phát sinh từ quá trình lặp và không đưa ra được đánh giá sai số tổng hợp (total error) kết hợp giữa sai số lặp liên tục và sai số rời rạc hóa.
- Đối với bài toán biên với điều kiện tích phân, phương trình vi tích phân và phương trình vi phân hàm, các kết quả công bố trước đây còn rời rạc, một số công trình có sai sót trong việc chứng minh cấp hội tụ hoặc đặt ra các giả thiết Lipschitz không thỏa mãn trong miền không bị chặn.
Mục tiêu nghiên cứu của luận án là nghiên cứu sự tồn tại, tính duy nhất, tính dương, tính đơn điệu của nghiệm và xây dựng phương pháp lặp giải tích - số hữu hiệu để tìm nghiệm cho một số bài toán biên đối với phương trình vi phân phi tuyến cấp cao, phương trình vi tích phân và phương trình vi phân hàm.
Các nội dung nghiên cứu chính được đánh số theo văn bản luận án gồm:
- Nội dung 1 (Content 1): Sự tồn tại, tính duy nhất của nghiệm và các phương pháp lặp cho một số bài toán biên đối với phương trình vi phân phi tuyến cấp ba.
- Nội dung 2 (Content 2): Sự tồn tại, tính duy nhất của nghiệm và các phương pháp lặp cho một số bài toán đối với phương trình vi phân phi tuyến cấp ba và cấp bốn với điều kiện biên tích phân.
- Nội dung 3 (Content 3): Sự tồn tại, tính duy nhất của nghiệm và các phương pháp lặp cho một số bài toán biên đối với phương trình vi tích phân và phương trình vi phân hàm.
Đối tượng và phạm vi nghiên cứu:
- Đối tượng nghiên cứu: Lớp bài toán biên hai điểm cho phương trình vi phân phi tuyến đầy đủ cấp ba dạng $u'''(t) = f(t, u(t), u'(t), u''(t))$; bài toán biên tích phân cho phương trình vi phân phi tuyến cấp ba và cấp bốn; bài toán biên cho phương trình vi tích phân phi tuyến cấp bốn và phương trình vi phân hàm pantograph có trễ tỉ lệ.
- Phạm vi nghiên cứu: Xét trên các đoạn hữu hạn (chủ yếu là đoạn $[0, 1]$ và $[a, b]$); không gian hàm liên tục $C[0, 1]$, không gian metric đầy đủ và không gian lưới đều $\bar{\omega}_h = {t_i = ih, h = 1/N, i = 0, 1, \dots, N}$; việc kiểm chứng tính toán thực nghiệm số được lập trình thực hiện trên môi trường phần mềm MATLAB.
Tổng quan tài liệu và vị trí của luận án
Luận án đã tổng quan hệ thống các công trình nghiên cứu trong và ngoài nước liên quan đến các hướng đề tài:
Về bài toán biên hai điểm cho phương trình vi phân phi tuyến cấp ba:
- Các phương pháp nghiên cứu tính định tính (sự tồn tại, duy nhất, tính dương) trước đây sử dụng phương pháp nghiệm trên - nghiệm dưới và kỹ thuật đơn điệu (Yao và Feng [19], Feng và Liu [17], [7, 15, 16, 18]), nguyên lý thác triển Leray-Schauder (Hopkins và Kosmatov [20]), lý thuyết điểm bất động trên nón (Guo và Yang [21, 46], Sun và cộng sự [24], Li và Li [22], Bai [23]), hoặc điều kiện dải chắn (Kelevedjiev và Todorov [25]).
- Các phương pháp giải tích và số: nhúng hàm Green vào phép lặp Picard hoặc Krasnoselskii-Mann [26], phương pháp phân rã Adomian [27, 28], khai triển chuỗi Taylor (He 2020 [29]), phương pháp spline bậc ba (Al Said và cộng sự [30]), spline bậc bốn (Noor và cộng sự [31], El-Danaf [34]), B-spline bậc bốn và spline bậc năm [32, 33], spline đa thức bậc bốn của Pandey (2016 [35]), sơ đồ sai phân của Pandey (2017 [36]), và phương pháp spline kết hợp của Chaurasia và cộng sự (2019 [37]).
Về bài toán biên với điều kiện biên tích phân:
- Phương trình vi phân cấp ba có điều kiện tích phân được khảo sát bởi Boucherif và cộng sự [44] (dùng toán tử thu gọn và toán tử co), Sun và Li (2010 [24], dùng định lý Krasnoselskii), Guo, Liu và Liang (2012 [45], dùng chỉ số điểm bất động trên nón), Guo và Yang (2013 [46]), Guendouz và cộng sự [47].
- Phương trình cấp bốn có điều kiện tích phân được nghiên cứu bởi Benaicha và Haddouchi [48], Zhang và Ge [58], Li và cộng sự (2013 [54]), Lv và cộng sự (2015 [55]).
Về bài toán biên cho phương trình vi tích phân (IDE) và vi phân hàm (FDE):
- Đối với IDE cấp hai và cấp bốn: Tahernezhad và Jalilian (2020 [65], dùng exponential spline đạt bậc sai số $O(h^2)$), Chen và cộng sự [60, 69] (dùng Galerkin đa tỷ lệ), Singh và Wazwaz [63] (Adomian kết hợp hàm Green cho phương trình Volterra), nghiên cứu cho phương trình Fredholm [59], phương pháp lặp và phổ cho IDE phi tuyến [12, 68], công trình của Đặng và Nguyễn [11] cho IDE cấp bốn với số hạng phi địa phương $M(\int_0^L |u'(t)|^2 dt)$, Wang [66] cho bài toán cấp bốn phát sinh từ mô hình cầu treo [70, 71] và lý thuyết lượng tử [72].
- Đối với FDE: Reutskiy (2015 [74], phương pháp collocation cho phương trình pantograph tuyến tính), Bica và cộng sự (2016 [75], dùng xấp xỉ liên tiếp với cubic spline), Khuri và Sayfy (2018 [76], dùng hàm Green lặp), và các thuật toán trí tuệ tính toán/mạng nơ-ron [77, 78].
Khoảng trống và vị trí học thuật của luận án: Luận án định vị giải quyết các hạn chế tồn đọng thông qua một hướng tiếp cận thống nhất:
- Chuyển đổi các bài toán biên về phương trình toán tử đối với số hạng phi tuyến $f$ ở vế phải phương trình vi phân thay vì thiết lập toán tử đối với hàm nghiệm cần tìm $u$.
- Thiết lập các điều kiện đủ đơn giản, dễ kiểm chứng trong các miền bị chặn $D_M$ thay vì đòi hỏi các điều kiện ngặt trên toàn không gian $\mathbb{R}^n$.
- Khắc phục các sai sót học thuật trong các công bố quốc tế trước đây, cụ thể: chỉ ra việc Bica và cộng sự [75] ngầm giả thiết sai về tính liên tục của các đạo hàm riêng cấp 3 và cấp 4 của hàm Green $\partial^3 G/\partial s^3, \partial^4 G/\partial s^4$ trên toàn miền $[a, b] \times [a, b]$ (dẫn đến ước lượng sai cấp hội tụ, đã phải đính chính trong [79]), và chỉ ra mâu thuẫn trong kết luận về sự tồn tại nghiệm dương của Li và Li [22].
- Xây dựng các thuật toán lặp rời rạc kết hợp các công thức cầu phương thích hợp (hình thang và Simpson cải biên) và thiết lập đánh giá sai số toàn phần kết hợp giữa quá trình lặp và kích thước lưới $h$.
Cơ sở lý thuyết và phương pháp nghiên cứu
Luận án sử dụng các công cụ lý thuyết toán học giải tích, giải tích hàm phi tuyến và phương pháp số:
Các định lý điểm bất động:
- Định lý điểm bất động Schauder: Giả sử $U$ là tập con lồi, đóng, bị chặn trong không gian Banach $X$ và $T: U \to U$ là toán tử compact. Khi đó $T$ có ít nhất một điểm bất động.
- Định lý điểm bất động Krasnoselskii (trên nón): Áp dụng cho toán tử hoàn toàn liên tục trên nón $P$ trong không gian Banach để chứng minh sự tồn tại nghiệm dương thông qua điều kiện dãn nón và nén nón dạng chuẩn.
- Định lý điểm bất động Banach (Nguyên lý ánh xạ co): Áp dụng cho toán tử co $T: M \subset X \to M$ với hệ số co $q \in (0, 1)$ trên tập đóng $M$ của không gian metric đầy đủ $(X, d)$, đi kèm các công thức đánh giá sai số tiên nghiệm $d(x_n, x) \le \frac{q^n}{1-q} d(x_0, x_1)$ và sai số hậu nghiệm $d(x_{n+1}, x) \le \frac{q}{1-q} d(x_n, x_{n+1})$.
Biểu diễn hàm Green và các chuẩn tích phân: Nghiệm của bài toán tuyến tính tương ứng được biểu diễn qua hàm Green $G(t, s)$ và các đạo hàm riêng $G_1(t, s) = G't(t, s)$, $G_2(t, s) = G''{tt}(t, s)$: $$u(t) = \int_0^1 G(t, s)\phi(s)ds, \quad u'(t) = \int_0^1 G_1(t, s)\phi(s)ds, \quad u''(t) = \int_0^1 G_2(t, s)\phi(s)ds$$ Các hằng số bị chặn được xác định bởi: $$M_0 = \max_{0 \le t \le 1} \int_0^1 |G(t, s)|ds, \quad M_1 = \max_{0 \le t \le 1} \int_0^1 |G_1(t, s)|ds, \quad M_2 = \max_{0 \le t \le 1} \int_0^1 |G_2(t, s)|ds$$ Miền compact tương ứng với số $M > 0$ được định nghĩa: $$D_M = {(t, x, y, z) \mid 0 \le t \le 1, |x| \le M_0 M, |y| \le M_1 M, |z| \le M_2 M}$$
Phương pháp số và kỹ thuật rời rạc:
- Sử dụng lưới đều $\bar{\omega}_h = {t_i = ih, h = 1/N, i = 0, 1, \dots, N}$.
- Áp dụng công thức cầu phương hình thang (Trapezoidal rule) có bậc chính xác $O(h^2)$ và công thức Simpson tổng quát có bậc chính xác $O(h^4)$.
- Xây dựng công thức Simpson cải biên (Modified Simpson rule) xử lý riêng tại các điểm nút chẵn và lẻ để khắc phục tính không trơn của đạo hàm hàm Green trên đường chéo $s = t$, đem lại bậc hội tụ rời rạc $O(h^3)$.
- Tiêu chuẩn dừng của quá trình lặp: $|\Phi_{k+1} - \Phi_k| \le TOL$ với các mức dung sai $TOL \in {10^{-4}, 10^{-6}, 10^{-10}}$ hoặc $|U_m - u| \le h^2$.
Nội dung chính theo từng chương
Chương 1: Kiến thức chuẩn bị (Preliminaries)
Chương này trình bày các kiến thức nền tảng phục vụ cho các kết quả ở các chương tiếp theo:
- Trình bày chi tiết định lý điểm bất động Brouwer, Schauder, Krasnoselskii trên nón và Banach, cùng các công thức đánh giá sai số.
- Thiết lập định nghĩa, điều kiện tồn tại duy nhất và dạng giải tích tường minh của hàm Green $G(x, t)$ cho các bài toán biên mẫu:
- Bài toán cấp hai: $u''(x) = \phi(x), u(0)=0, u(1)=0$, hàm Green $G(x, t) = x(1-t)$ khi $0 \le x \le t \le 1$ và $G(x, t) = t(1-x)$ khi $0 \le t \le x \le 1$.
- Bài toán cấp ba: $u'''(x) = \phi(x), u(0)=0, u'(0)=0, u'(1)=0$, hàm Green $G(x, t) = \frac{1}{2}x^2(t-1)$ khi $0 \le x \le t \le 1$ và $G(x, t) = \frac{1}{2}(2xt - x^2 - t)$ khi $0 \le t \le x \le 1$.
- Bài toán cấp bốn: $u^{(4)}(x) = \phi(x), u(0)=u(1)=u''(0)=u''(1)=0$, hàm Green $G(x, t) = \frac{1}{6}t(x-1)(t^2 - 2x + x^2)$ khi $0 \le t \le x \le 1$.
- Nhắc lại các công thức cầu phương đại số: công thức hình thang đơn và dạng lặp có số hạng dư bậc hai $-\frac{1}{12}(b-a)h^2 f''(\xi)$; công thức Simpson đơn và dạng lặp có số hạng dư bậc bốn $-\frac{1}{180}(b-a)h^4 f^{(4)}(\xi)$.
Chương 2: Kết quả tồn tại và phương pháp lặp cho bài toán biên hai điểm phi tuyến cấp ba (Existence results and iterative method for two-point third order nonlinear BVPs)
Chương 2 trình bày các kết quả về lý thuyết định tính và thuật toán số cho phương trình vi phân phi tuyến đầy đủ cấp ba: $$u'''(t) = f(t, u(t), u'(t), u''(t)), \quad 0 < t < 1$$ với các điều kiện biên tuyến tính hai điểm tổng quát $B_1[u] = B_2[u] = B_3[u] = 0$.
Phương pháp toán tử và các định lý tồn tại, duy nhất:
- Luận án đưa bài toán biên về phương trình điểm bất động $A\phi = \phi$ đối với toán tử $A\phi(t) = f(t, u(t), u'(t), u''(t))$ với $u(t), u'(t), u''(t)$ biểu diễn qua hàm Green.
- Định lý 2.2 (Tồn tại nghiệm): Nếu tồn tại $M > 0$ sao cho $|f(t, x, y, z)| \le M$ trên $D_M$, thì bài toán có ít nhất một nghiệm $u(t)$ thỏa mãn $|u(t)| \le M_0 M, |u'(t)| \le M_1 M, |u''(t)| \le M_2 M$.
- Định lý 2.3 (Nghiệm không âm, đơn điệu): Xác định điều kiện đủ về dấu của hàm Green $\sigma(G)$ và đạo hàm $\sigma(G_1)$ để nghiệm $u(t)$ đơn điệu tăng hoặc đơn điệu giảm và không âm.
- Định lý 2.4 (Duy nhất nghiệm): Nếu $f$ thỏa mãn điều kiện Lipschitz với các hằng số $L_0, L_1, L_2$ trên $D_M$ và thỏa mãn điều kiện co: $$q := L_0 M_0 + L_1 M_1 + L_2 M_2 < 1$$ thì bài toán có duy nhất một nghiệm trong $D_M$.
- Thuật toán lặp liên tục: Bắt đầu từ $\phi_0(t) = 0$, tính $u_k, y_k, z_k$ theo các tích phân chứa hàm Green, sau đó cập nhật $\phi_{k+1}(t) = f(t, u_k(t), y_k(t), z_k(t))$. Tốc độ hội tụ thỏa mãn đánh giá: $$|u_k - u| \le M_0 p_k, \quad |u'_k - u'| \le M_1 p_k, \quad |u''_k - u''| \le M_2 p_k, \quad \text{với } p_k = \frac{q^k}{1-q}|\phi_1 - \phi_0|$$
Xây dựng và đánh giá hai phương pháp lặp rời rạc:
- Phương pháp lặp rời rạc 1 (Method 1): Rời rạc hóa tích phân bằng công thức hình thang, chứng minh sai số toàn phần đạt cấp $O(h^2)$: $$|U_k - u| \le M_0 p_k d + O(h^2), \quad |Y_k - u'| \le M_1 p_k d + O(h^2), \quad |Z_k - u''| \le M_2 p_k d + O(h^2)$$
- Phương pháp lặp rời rạc 2 (Method 2): Rời rạc hóa tích phân bằng công thức Simpson cải biên cho điểm nút chẵn và lẻ, chứng minh sai số toàn phần đạt cấp $O(h^3)$: $$|U_k - u| \le M_0 p_k d + O(h^3), \quad |Y_k - u'| \le M_1 p_k d + O(h^3), \quad |Z_k - u''| \le M_2 p_k d + O(h^3)$$
Các ví dụ minh họa và kiểm chứng số liệu:
- Ví dụ 2.1: $u'''(t) = -e^{u(t)}$ với $u(0)=u'(0)=u'(1)=0$. Các hằng số $M_0 = 1/12, M_1 = 1/8, M_2 = 1/2$. Chọn $M = 1.5, L_0 = e^{1.5/12} \approx 1.1331 \implies q \approx 0.0944 < 1$. Sau 5 bước lặp trên lưới $h=0.01$, phương pháp lặp cho nghiệm dương đơn điệu bị chặn bởi $0.125$.
- Ví dụ 2.2: $u'''(t) = -\frac{5u^3+4u+3}{u^2+1}$, chọn $M = 4.1$, phương pháp lặp hội tụ sau 8 bước lặp.
- Ví dụ 2.3: $u'''(t) = -e^{u(t)} - e^{u'(t)}$, chọn $M = 2.7$, hội tụ sau 9 bước lặp.
- Ví dụ 2.4: $u'''(t) = -\frac{1}{36}(u')^2 + \frac{1}{24}uu'' + \frac{1}{4}t^2 - 6$, với $u(0)=u'(0)=u''(1)=0$. Nghiệm chính xác là $u(t) = -t^3 + 3t^2$. Với $h=0.01$, độ lệch cực đại đạt $3.3429 \times 10^{-6}$ sau 5 bước lặp.
- Ví dụ 2.5: $u'''(t) = \frac{1}{18}(u')^2 - \frac{1}{12}uu'' + \frac{1}{2}t + \frac{11}{2}$, nghiệm chính xác $u(t) = t^3 - 3t^2 + 3t$. Với $h=0.01$, độ lệch cực đại đạt $3.9318 \times 10^{-6}$ sau 6 bước lặp.
- Ví dụ 2.7 (So sánh với Pandey [35]): Bài toán $u'''(x) = x^4 u(x) - u^2(x) + g(x)$ có nghiệm chính xác $u^*(x) = (x-1)\sin(x)$.
| Số điểm lưới ($N$) | Số bước lặp ($K$) | Sai số Method 1 ($Errortrap$) | Bậc hội tụ Method 1 | Sai số Method 2 ($ErrorSimp$) | Bậc hội tụ Method 2 | Kết quả của Pandey [35] |
|---|---|---|---|---|---|---|
| 8 | 7 | $9.3248 \times 10^{-4}$ | - | $7.6710 \times 10^{-5}$ | - | $1.2373 \times 10^{-3}$ |
| 16 | 7 | $2.3168 \times 10^{-4}$ | 2.0089 | $1.0345 \times 10^{-5}$ | 2.8906 | $1.4312 \times 10^{-4}$ |
| 32 | 7 | $5.7828 \times 10^{-5}$ | 2.0023 | $1.3411 \times 10^{-6}$ | 2.9475 | $1.2941 \times 10^{-5}$ |
| 64 | 7 | $1.4451 \times 10^{-5}$ | 2.0006 | $1.6983 \times 10^{-7}$ | 2.9813 | $2.2619 \times 10^{-6}$ |
| 128 | 7 | $3.6122 \times 10^{-6}$ | 2.0002 | $2.1352 \times 10^{-8}$ | 2.9917 | - |
- Ví dụ 2.8 (So sánh với Pandey [36]): Bài toán $u'''(x) = -x u''(x) - 6x^2 + 3x - 6$ có nghiệm chính xác $u(x) = x^2(\frac{3}{2} - x)$.
| Số điểm lưới ($N$) | Số bước lặp ($K$) | Sai số Method 1 ($Errortrap$) | Sai số Method 2 ($ErrorSimp$) | Sai số Pandey [36] | Số lần lặp Pandey [36] |
|---|---|---|---|---|---|
| 8 | 11 | $1.1578 \times 10^{-4}$ | $1.2023 \times 10^{-5}$ | $1.3417 \times 10^{-3}$ | 53 |
| 16 | 11 | $2.8687 \times 10^{-5}$ | $1.5645 \times 10^{-6}$ | $3.4172 \times 10^{-4}$ | 5 |
| 32 | 11 | $7.1557 \times 10^{-6}$ | $1.9904 \times 10^{-7}$ | $8.6315 \times 10^{-5}$ | 3 |
| 64 | 11 | $1.7879 \times 10^{-6}$ | $2.5083 \times 10^{-8}$ | $2.1691 \times 10^{-5}$ | 4 |
- Ví dụ 2.9: Phương trình $u'''(x) = -e^{u(x)} - e^{u'(x)} - \frac{1}{10}(u''(x))^2$ với $u(0)=u'(0)=u'(1)=0$, nghiệm giải tích chưa biết. Các điều kiện co thỏa mãn với $M=3, L_0=1.2840, L_1=1.4550, L_2=0.3000 \implies q = 0.4389 < 1$, phương pháp lặp xác lập nghiệm số ổn định.
Chương 3: Kết quả tồn tại và phương pháp lặp cho một số phương trình vi phân thường phi tuyến với điều kiện biên tích phân (Existence results and iterative method for some nonlinear ODEs with integral boundary conditions)
Chương 3 mở rộng lý thuyết và thuật toán cho các bài toán biên phi tuyến chứa điều kiện biên tích phân:
- Khảo sát bài toán biên tích phân cấp ba đầy đủ: $$u'''(t) = f(t, u(t), u'(t), u''(t)), \quad 0 < t < 1$$ với các điều kiện biên tích phân dạng $u(0)=0, u'(0) - a u''(0) = \int_0^1 h_1(u(s), u'(s))ds, u'(1) + b u''(1) = \int_0^1 h_2(u(s), u'(s))ds$ hoặc $u(0)=u'(0)=0, u(1)=\int_0^1 g(t)u(t)dt$. Thiết lập hệ toán tử phi tuyến tương ứng và chứng minh sự tồn tại, duy nhất nghiệm thông qua định lý điểm bất động Schauder và Banach.
- Khảo sát bài toán biên tích phân cấp bốn đầy đủ: $$u^{(4)}(t) = f(t, u(t), u'(t), u''(t), u'''(t)), \quad t \in [0, 1]$$ với điều kiện biên chứa số hạng tích phân tại biên. Xây dựng phương pháp lặp trên mức liên tục và thiết lập sơ đồ lặp rời rạc, chứng minh sự hội tụ của thuật toán số và tính toán các ví dụ thực nghiệm.
Chương 4: Kết quả tồn tại và phương pháp lặp cho phương trình vi tích phân và phương trình vi phân hàm (Existence results and iterative method for integro-differential and functional differential equations)
Chương 4 tập trung vào hai lớp phương trình mở rộng:
- Phương trình vi tích phân phi tuyến (IDE): Nghiên cứu bài toán biên cho phương trình vi tích phân cấp bốn phi tuyến dạng Fredholm và Volterra. Xây dựng toán tử tích phân tương đương, thiết lập điều kiện tồn tại duy nhất nghiệm và thuật toán số xấp xỉ.
- Phương trình vi phân hàm (FDE): Xét phương trình vi phân hàm có trễ dạng pantograph. Luận án phân tích chi tiết và chỉ ra thiếu sót của các nghiên cứu trước (như Bica và cộng sự [75]), thiết lập các điều kiện Lipschitz chặt chẽ trên miền compact $D_M$, chứng minh định lý tồn tại và duy nhất nghiệm, đồng thời thiết lập phương pháp lặp giải tích - số với đánh giá sai số hội tụ chính xác.
Kết quả và những đóng góp mới
Luận án tự tổng kết các đóng góp học thuật và thực tiễn chính:
- Đóng góp mới về phương pháp tiếp cận: Đề xuất hướng tiếp cận chuyển đổi các bài toán biên đối với phương trình vi phân phi tuyến cấp ba, cấp bốn, phương trình có điều kiện biên tích phân, IDE và FDE về các phương trình toán tử đối với số hạng phi tuyến $f$ ở vế phải (thay vì đối với hàm cần tìm $u$). Cách tiếp cận này làm giảm nhẹ đáng kể độ phức tạp của các điều kiện giải được.
- Thiết lập các điều kiện định tính khả thi: Đưa ra các điều kiện đủ đơn giản, dễ kiểm chứng trong miền bị chặn $D_M$ để đảm bảo sự tồn tại, duy nhất, tính dương và tính đơn điệu của nghiệm. Các điều kiện này yếu hơn và thuận tiện hơn khi kiểm tra so với các điều kiện Lipschitz hoặc điều kiện tăng toàn cục trên toàn bộ không gian của các tác giả trước.
- Phát hiện và đính chính các hạn chế học thuật trong tài liệu quốc tế: Chỉ ra giả định không hợp lệ về tính liên tục của đạo hàm hàm Green cấp cao trong công trình của Bica và cộng sự [75], đồng thời chỉ ra kết luận chưa chính xác về sự tồn tại nghiệm dương trong công trình của Li và Li [22].
- Xây dựng các phương pháp lặp số và đánh giá sai số tổng hợp: Thiết lập hai phương pháp lặp rời rạc (Method 1 dùng quy tắc hình thang đạt bậc hội tụ 2, Method 2 dùng quy tắc Simpson cải biên đạt bậc hội tụ 3). Thiết lập công thức tường minh cho sai số toàn phần kết hợp giữa sai số lặp lý thuyết $M_i p_k d$ và sai số bước lưới rời rạc $O(h^2)$ hoặc $O(h^3)$.
- Công bố khoa học: Các kết quả chính của luận án đã được công bố trong 6 bài báo khoa học (từ [AL1] đến [AL6]).
Hạn chế và hướng nghiên cứu tiếp
Hạn chế của luận án:
- Luận án giới hạn phạm vi khảo sát trên các miền hữu hạn (đoạn $[0, 1]$ hoặc $[a, b]$) và các bài toán biên hai điểm với các điều kiện biên tuyến tính hoặc tích phân xác định.
- Các điều kiện co và điều kiện tồn tại nghiệm dựa trên các hằng số Lipschitz $L_i$ và các giá trị cực đại tích phân của hàm Green $M_i$, đòi hỏi hàm phi tuyến $f$ phải thỏa mãn tính liên tục hoặc khả vi đến cấp nhất định trong miền bị chặn $D_M$.
Hướng nghiên cứu tiếp:
- Mở rộng phương pháp tiếp cận toán tử cho các bài toán biên đa điểm (multi-point BVPs) hoặc các bài toán biên trên khoảng không bị chặn.
- Áp dụng phương pháp lặp đề xuất cho các phương trình vi phân có cấp cao hơn (cấp 5, cấp 6) và các hệ phương trình vi phân phi tuyến phát sinh từ các mô hình vật lý phức tạp.
- Nghiên cứu mở rộng phương pháp lặp rời rạc với các công thức cầu phương bậc cao hơn hoặc kết hợp với các kỹ thuật thích nghi lưới.
Giá trị tham khảo
Luận án có giá trị tham khảo chuyên môn cho các đối tượng và nội dung cụ thể:
- Nghiên cứu sinh và học viên cao học ngành Toán ứng dụng, Toán giải tích và Giải tích số: Tham khảo kỹ thuật chuyển đổi bài toán biên về phương trình toán tử cho số hạng phi tuyến, cách sử dụng các định lý điểm bất động (Schauder, Banach, Krasnoselskii) trên miền bị chặn và kỹ thuật đánh giá sai số toàn phần.
- Giảng viên và nhà nghiên cứu trong lĩnh vực phương pháp số và mô hình hóa: Tham khảo thuật toán lặp rời rạc 1 và rời rạc 2 (đặc biệt là công thức Simpson cải biên xử lý điểm kỳ dị/không trơn của hàm Green trên đường chéo) để lập trình giải số các phương trình vi phân phi tuyến cấp cao, phương trình có điều kiện tích phân và phương trình vi tích phân.
Câu hỏi thường gặp
1. Hướng tiếp cận toán tử của luận án có điểm gì khác biệt cốt lõi so với các nghiên cứu trước đây?
Thay vì đặt phương trình toán tử trực tiếp cho hàm cần tìm $u(t)$ (đòi hỏi các điều kiện ngặt trên toàn không gian $\mathbb{R}^n$), luận án chuyển bài toán biên về phương trình toán tử điểm bất động $A\phi = \phi$ đối với chính hàm phi tuyến $\phi(t) = f(t, u(t), u'(t), u''(t))$. Nhờ đó, việc khảo sát tính compact và tính co của toán tử được thực hiện trong không gian $C[0, 1]$ trên miền đóng bị chặn $D_M$, giúp các điều kiện tồn tại và duy nhất nghiệm trở nên đơn giản và dễ kiểm chứng hơn.
2. Phương pháp lặp rời rạc 1 và Phương pháp lặp rời rạc 2 khác nhau như thế nào về kỹ thuật tính và cấp chính xác?
Phương pháp 1 sử dụng công thức cầu phương hình thang thông thường để tính gần đúng các tích phân chứa hàm Green tại mỗi bước lặp, đạt cấp chính xác rời rạc là $O(h^2)$. Phương pháp 2 sử dụng công thức Simpson cải biên phân chia xử lý riêng cho các nút lưới chẵn ($i = 2m$) và nút lưới lẻ ($i = 2m+1$) để bù đắp điểm gián đoạn của đạo hàm hàm Green trên đường chéo $s = t$, nâng cấp chính xác rời rạc lên $O(h^3)$.
3. Đánh giá sai số toàn phần (total error) của phương pháp lặp rời rạc được cấu thành từ những thành phần nào?
Sai số toàn phần của nghiệm rời rạc $U_k$ so với nghiệm chính xác $u$ gồm hai thành phần độc lập: sai số của quá trình lặp liên tục $M_0 p_k d$ (với $p_k = \frac{q^k}{1-q}$, suy giảm theo cấp số nhân với công bội co $q < 1$) và sai số rời rạc hóa của công thức tích phân số ($O(h^2)$ đối với Method 1 hoặc $O(h^3)$ đối với Method 2). Khi hai thành phần này đạt trạng thái cân bằng, việc tiếp tục tăng số điểm lưới $N$ (giảm $h$) sẽ không làm tăng thêm độ chính xác của nghiệm số.
4. Luận án đã chỉ ra sai sót học thuật nào trong công trình của Bica và cộng sự (2016) và kết luận của Li và Li (2020)?
Luận án chỉ ra rằng trong [75], Bica và cộng sự đã mặc định sai khi coi các đạo hàm riêng $\partial^3 G/\partial s^3, \partial^4 G/\partial s^4$ liên tục trên $[a, b] \times [a, b]$, trong khi thực tế $\partial^3 G/\partial s^3$ bị gián đoạn tại đường chéo $s = t$, dẫn đến ước lượng sai bậc hội tụ $O(h^4)$ (sai sót này đã phải đính chính trong [79]). Đồng thời, tại Ví dụ 2.6 cho phương trình $u''' = u^3 + u(u')^2 + u(u'')^2$, luận án chứng minh bài toán chỉ có nghiệm tầm thường $u(t) \equiv 0$ là nghiệm không âm duy nhất, bác bỏ kết luận về sự tồn tại nghiệm dương của Li và Li trong [22].
5. Kết quả thực nghiệm số của luận án thể hiện ưu điểm gì so với phương pháp của Pandey (2016, 2017)?
So với phương pháp spline đa thức bậc bốn của Pandey [35] (sử dụng phép lặp Newton-Raphson) và sơ đồ sai phân của Pandey [36] (sử dụng phép lặp Gauss-Seidel đòi hỏi tới 53 lần lặp ở lưới thô), phương pháp lặp của luận án đạt sai số nhỏ hơn từ 1 đến 3 bậc độ lớn với cùng số điểm chia $N$, đồng thời số bước lặp cố định và khối lượng tính toán thấp hơn.
Kết luận
Luận án của nghiên cứu sinh Đặng Quang Long đã giải quyết một cách có hệ thống các bài toán biên đối với phương trình vi phân phi tuyến cấp ba, cấp bốn, phương trình có điều kiện biên tích phân, phương trình vi tích phân và phương trình vi phân hàm. Thông qua việc đề xuất phương pháp toán tử cho số hạng phi tuyến kết hợp các định lý điểm bất động trên miền bị chặn, tác giả đã thiết lập các điều kiện tồn tại, duy nhất, tính dương và đơn điệu của nghiệm một cách tường minh, dễ kiểm chứng. Đồng thời, luận án đã xây dựng hai sơ đồ lặp số đạt cấp chính xác $O(h^2)$ và $O(h^3)$, cung cấp các đánh giá sai số toàn phần chặt chẽ và kiểm chứng độ chính xác thông qua các thử nghiệm tính toán trên MATLAB.