Tổng quan luận án
Luận án tiến sĩ toán học của nghiên cứu sinh Nguyễn Thanh Hường, với đề tài "Giải gần đúng một số bài toán biên phi tuyến cho phương trình vi phân cấp bốn", thuộc chuyên ngành Toán ứng dụng (mã số 9 46 01 12), được hoàn thành năm 2019 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, dưới sự hướng dẫn khoa học của GS. Đặng Quang Á và TS. Vũ Vinh Quang.
Nhiều bài toán trong vật lý và cơ học, tiêu biểu là lý thuyết uốn cong của dầm và của bản trên nền đàn hồi, được mô hình hóa bởi các bài toán biên cho phương trình vi phân thường hoặc phương trình đạo hàm riêng phi tuyến cấp bốn. Các điều kiện biên tại các biên giới hạn rất đa dạng, bao gồm dạng gối - tựa đơn giản, ngàm - tự do, điều kiện Dirichlet, điều kiện biên tổ hợp, tuần hoàn và các điều kiện biên phi tuyến. Về mặt phân loại, phương trình vi phân cấp bốn được chia thành hai nhóm: phương trình vi phân cấp bốn địa phương và phương trình vi phân cấp bốn không địa phương (phương trình có chứa thành phần tích phân, thường được gọi là phương trình loại Kirchhoff).
Trước khi luận án được thực hiện, các phương pháp tiếp cận chủ yếu trong các tài liệu chuyên ngành gồm có phương pháp biến phân, phương pháp nghiệm trên và nghiệm dưới, lý thuyết bậc tô pô (bậc Brouwer, Leray-Schauder), phương pháp chuỗi và phương pháp sai phân hữu hạn. Tuy nhiên, các cách tiếp cận này tồn tại các khoảng trống kỹ thuật:
- Phương pháp biến phân thường chỉ thiết lập sự tồn tại của nghiệm hoặc đa nghiệm, đòi hỏi các giả thiết khắt khe về điều kiện tăng trưởng tại vô cùng của hàm vế phải, không khẳng định được tính duy nhất và không cung cấp phương pháp giải hay ví dụ nghiệm cụ thể.
- Phương pháp nghiệm trên và nghiệm dưới đòi hỏi sự tồn tại tiên nghiệm của cặp nghiệm trên - nghiệm dưới (vốn rất khó tìm trong thực tế) và buộc hàm phi tuyến phải thỏa mãn các điều kiện phức tạp như điều kiện Nagumo hoặc tính đơn điệu chặt chẽ.
- Phương pháp chuỗi thường công nhận tiên nghiệm tính giải được duy nhất để tìm nghiệm dưới dạng tổng chuỗi.
- Phương pháp sai phân hữu hạn truyền thống thường tiến hành rời rạc hóa ngay từ đầu mà bỏ qua phân tích định tính liên tục, dẫn đến khó khăn trong việc đánh giá sự ổn định, sự hội tụ và sai số giữa nghiệm đúng và nghiệm xấp xỉ.
Nhằm khắc phục các hạn chế trên, mục tiêu nghiên cứu của luận án được xác định gồm ba nội dung:
- Nghiên cứu định tính (sự tồn tại nghiệm, tính duy nhất nghiệm, tính dương của nghiệm) bằng cách sử dụng các định lý điểm bất động và nguyên lý cực đại, không cần đến điều kiện tăng trưởng tại vô cùng hay điều kiện Nagumo của hàm vế phải.
- Xây dựng các phương pháp lặp giải bài toán và chứng minh sự hội tụ của phương pháp.
- Đưa ra các ví dụ minh họa cho các kết quả lý thuyết, trong đó có những ví dụ thể hiện ưu thế của phương pháp đề xuất so với phương pháp của một số tác giả khác.
Đối tượng và phạm vi nghiên cứu của luận án tập trung vào:
- Năm bài toán biên cho phương trình vi phân thường phi tuyến cấp bốn địa phương và không địa phương trên đoạn hữu hạn $[0, 1]$ (hoặc $[0, L]$, $[a, b]$) với các điều kiện biên tổ hợp, Dirichlet, gối - tựa đơn giản và phi tuyến.
- Hai bài toán biên cho phương trình đạo hàm riêng phi tuyến cấp bốn (phương trình song điều hòa và phương trình song điều hòa loại Kirchhoff) trên miền giới nội $\Omega \subset \mathbb{R}^K$ ($K \ge 2$, thực nghiệm số xét trên miền phẳng hai chiều).
Tổng quan tài liệu và vị trí của luận án
Tổng quan tài liệu trong luận án phân tích chi tiết các hướng tiếp cận chính trong và ngoài nước đối với bài toán biên cấp bốn:
-
Hướng nghiên cứu sử dụng phương pháp biến phân:
- T. Ma (2000) nghiên cứu bài toán biên phương trình vi phân cấp bốn loại Kirchhoff $u^{(4)}(x) - M(c \int_0^1 |u'(s)|^2 ds)u''(x) + f(x, u(x)) = 0$ trên $(0, 1)$ với điều kiện biên $u'(0) = u'(1) = 0$, $u'''(0) = -g(u(0))$, $u'''(1) = g(u(1))$, chứng minh sự tồn tại nghiệm với giả thiết $F(x, t) \to +\infty$ khi $|t| \to \infty$. Đến năm 2003, T. Ma tiếp tục công bố kết quả cho bài toán dạng này trên $(0, L)$ với điều kiện biên $u(0)=u'(0)=u''(L)=0$ kèm các tham số chặn tăng trưởng $m_0, \alpha_0, \beta_0$.
- Heidarkhani và các cộng sự khảo sát bài toán cấp bốn địa phương $u^{(4)}(x) = \lambda f(x, u) + \mu g(x, u) + p(u)$ với điều kiện biên $u(0)=u'(0)=u''(1)=0$, $u'''(1)=h(u(1))$, sử dụng cấu trúc hàm Carathéodory và hàm Lipschitz nhưng đi kèm các giả thiết phức tạp về sự tăng trưởng ở vô cùng.
- Pei nghiên cứu bài toán biên Navier cho phương trình song điều hòa $\Delta^2 u = f(x, u)$ trên miền $\Omega \subset \mathbb{R}^K$ ($K > 4$), chứng minh sự tồn tại ít nhất ba nghiệm không tầm thường dựa trên các điều kiện về giá trị riêng $\lambda_1, \lambda_k$ của toán tử $(\Delta^2, H^2(\Omega) \cap H_0^1(\Omega))$.
- An xét phương trình song điều hòa loại Kirchhoff $\Delta^2 u = M(c \int_\Omega |\nabla u|^2 dx)\Delta u + f(x, u)$ với điều kiện biên $u = \Delta u = 0$ trên $\Gamma$, đòi hỏi các điều kiện tăng trưởng (B1')–(B4') đối với hàm $f(x, t)$.
-
Hướng nghiên cứu sử dụng phương pháp nghiệm trên và nghiệm dưới:
- Z. Bai (2007) thiết lập sự tồn tại nghiệm cho phương trình $u^{(4)} = f(x, u, u', u'', u''')$ trên $(0, 1)$ với điều kiện biên $u(0)=u'(1)=u''(0)=u'''(1)=0$, đòi hỏi hàm $f$ thỏa mãn điều kiện Nagumo theo cặp nghiệm trên - nghiệm dưới $\alpha, \beta$ có $\alpha'' \le \beta''$.
- H. Feng và cộng sự (2009) mở rộng cho bài toán có điều kiện biên tổ hợp $au''(0) - bu'''(0) = 0$, $cu''(1) + du'''(1) = 0$, với giả thiết hàm $f$ giảm theo $u, y$, tăng chặt theo $z$ và thỏa mãn điều kiện Nagumo.
- Y. Wang (2006) khảo sát bài toán $\Delta(k(x)\Delta u) = f(x, u, \Delta u)$ trên miền bị chặn $\Omega \subset \mathbb{R}^K$, đòi hỏi tính đơn điệu của hàm vế phải theo biến $u$.
-
Hướng nghiên cứu sử dụng các định lý điểm bất động và lý thuyết bậc:
- R. Chow (1984) đưa bài toán Dirichlet đối với $u^{(4)} = f(x, u, u', u'', u''')$ về toán tử tích phân chứa hàm Green và đa thức bậc ba $P_3(x)$, áp dụng nguyên lý Schauder để chỉ ra sự tồn tại và định lý Banach để chứng minh duy nhất nghiệm cùng dãy lặp Picard.
- T. Ma và cộng sự (2006) áp dụng định lý điểm bất động Krasnosel'skii trên nón để chứng minh sự tồn tại nghiệm dương cho phương trình Kirchhoff cấp bốn.
- P. Amster (2008) sử dụng định lý Leray-Schauder kết hợp bậc Brouwer cho bài toán $u^{(4)} - Au'' + g(x, u) = 0$.
- M. Chang (2014) dùng bậc Leray-Schauder chứng minh nghiệm cho bài toán biên tổ hợp dưới điều kiện Nagumo.
-
Hướng nghiên cứu của các tác giả trong nước:
- Đặng Quang Á và cộng sự (trong các công trình từ [17] đến [24]) nghiên cứu định tính và phương pháp lặp giải các bài toán biên cho phương trình vi phân thường và đạo hàm riêng cấp bốn.
- Phạm Kỳ Anh (trong [10], [11]) khảo sát tính giải được, cấu trúc tập nghiệm và xấp xỉ nghiệm của bài toán biên tuần hoàn.
- Ngô Thị Kim Quy (2017) nghiên cứu thành công bài toán biên cho phương trình vi phân thường cấp bốn địa phương với điều kiện biên tuyến tính dạng gối - tựa đơn giản và ngàm - tự do bằng cách chuyển về phương trình toán tử đối với thành phần phi tuyến trên miền giới nội.
Khoảng trống luận án giải quyết: Tiếp tục phát triển hướng tiếp cận chuyển bài toán về phương trình toán tử đối với hàm vế phải hoặc hàm trung gian trên một miền giới nội, luận án giải quyết bài toán với các điều kiện biên phức tạp hơn (điều kiện biên tổ hợp, điều kiện biên Dirichlet, điều kiện biên phi tuyến, và phương trình đạo hàm riêng song điều hòa Kirchhoff). Bằng cách này, luận án không đòi hỏi hàm phi tuyến phải thỏa mãn điều kiện tăng trưởng tại vô cùng hay điều kiện Nagumo, không cần tìm trước nghiệm trên và nghiệm dưới, đồng thời bảo đảm tính duy nhất của nghiệm và thiết lập thuật toán lặp hội tụ với tốc độ cấp số nhân.
Cơ sở lý thuyết và phương pháp nghiên cứu
Luận án kết hợp phân tích giải tích hàm định tính và tính toán số định lượng dựa trên các công cụ toán học:
-
Cơ sở lý thuyết giải tích hàm và điểm bất động:
- Định lý điểm bất động Brouwer (1912): Áp dụng cho ánh xạ liên tục trên tập con khác rỗng, lồi, compact trong không gian Euclid hữu hạn chiều $\mathbb{R}^K$.
- Định lý điểm bất động Schauder (1930): Áp dụng cho toán tử compact trên tập con khác rỗng, lồi, đóng, bị chặn trong không gian Banach $X$, dùng để khẳng định sự tồn tại nghiệm.
- Định lý điểm bất động Banach (1922): Áp dụng cho ánh xạ co $T$ trên không gian metric đầy đủ với hệ số co $q < 1$. Định lý bảo đảm sự tồn tại duy nhất của điểm bất động, tính hội tụ của dãy lặp xấp xỉ liên tiếp $x_{k+1} = T x_k$, cùng các đánh giá sai số tiên nghiệm: $$d(x_k, x) \le \frac{q^k}{1-q} d(x_0, x_1)$$ và đánh giá sai số hậu nghiệm: $$d(x_{k+1}, x) \le \frac{q}{1-q} d(x_k, x_{k+1}).$$
- Lý thuyết hàm Green: Thiết lập biểu diễn giải tích tường minh của nghiệm qua tích phân toán tử với hàm Green $G(x, t)$ đối với phương trình vi phân tuyến tính cấp hai và cấp bốn có điều kiện biên thuần nhất tương ứng chỉ có nghiệm tầm thường.
-
Phương pháp định tính:
- Chuyển bài toán vi phân ban đầu về phương trình toán tử $\phi = A\phi$ đối với hàm vế phải $\phi(x) = f(x, u, u', u'', u''')$ hoặc hàm trung gian trên không gian Banach các hàm liên tục $C[0, 1]$ trang bị chuẩn cực đại $|.|$.
- Sử dụng các hằng số tích phân cực đại của hàm Green và đạo hàm của nó để giới hạn miền biến thiên $D_M$ của các biến số $(x, u, u', u'', u''')$.
- Đặt điều kiện Lipschitz cục bộ trên miền đóng $B[O, M]$ để chứng minh toán tử $A$ là ánh xạ co từ $B[O, M]$ vào chính nó, suy ra sự tồn tại duy nhất của nghiệm.
- Áp dụng nguyên lý cực đại và tính chất không âm của hàm Green để chứng minh tính dương / không âm của nghiệm bài toán.
-
Phương pháp số và kỹ thuật rời rạc:
- Đạo hàm số: Dùng đa thức nội suy Lagrange tại ba điểm (sai số $O(h^2)$) và năm điểm (sai số đạo hàm cấp một là $O(h^4)$, đạo hàm cấp bốn là $O(h)$).
- Tích phân số: Công thức hình thang với sai số toàn phần $\frac{M_2 (b-a)}{12} h^2$; công thức Simpson một chiều với sai số $\frac{M_4 (b-a)}{2880} h^4$; công thức Simpson hai lớp trên lưới chữ nhật cho phương trình đạo hàm riêng với sai số $O(h^4) + O(k^4)$.
- Lược đồ sai phân cấp bốn: Xấp xỉ phương trình Poisson hai chiều $\Delta u = -f$ bằng toán tử sai phân cấp bốn $\Lambda u = ux_1 x_1 + ux_2 x_2$ thông qua vế phải điều chỉnh $\phi = f + \frac{h_1^2}{12} L_1 f + \frac{h_2^2}{12} L_2 f$.
- Phương pháp giải hệ phương trình lưới:
- Phương pháp truy đuổi (thuật toán Thomas) giải hệ phương trình ba điểm với ma trận ba đường chéo trội, gồm biến thể truy đuổi phải và truy đuổi từ hai phía, với khối lượng tính toán xấp xỉ $8N$ phép toán.
- Phương pháp rút gọn hoàn toàn (cyclic reduction) giải hệ phương trình véc tơ ba điểm $-Y_{j-1} + C Y_j - Y_{j+1} = F_j$ trên lưới chữ nhật $N = 2^n$ bước thông qua $n-1$ bước khử liên tiếp.
-
Nguồn tư liệu và công cụ thực nghiệm:
- Dữ liệu thực nghiệm được tạo ra từ việc mô phỏng tính toán các bài toán biên đã biết nghiệm giải tích chính xác và các bài toán chưa biết nghiệm giải tích.
- Môi trường lập trình và thử nghiệm số: Phần mềm MATLAB phiên bản 7.0 chạy trên máy tính cá nhân (cấu hình phần cứng: CPU Intel Core i3, bộ nhớ RAM 4GB).
Nội dung chính theo từng chương
Chương 1: Kiến thức bổ trợ
Chương 1 hệ thống hóa toàn bộ cơ sở giải tích và giải tích số phục vụ cho việc nghiên cứu định tính và định lượng ở hai chương tiếp theo:
- Trình bày chi tiết lý thuyết điểm bất động: Định nghĩa không gian có tính chất điểm bất động, Định lý điểm bất động Brouwer trên $\mathbb{R}^K$, Định lý điểm bất động Schauder cho toán tử compact, và Định lý điểm bất động Banach cho ánh xạ co trên không gian metric đầy đủ kèm các đánh giá sai số tiên nghiệm và hậu nghiệm.
- Xây dựng hàm Green cho các bài toán biên vi phân cấp hai và cấp bốn:
- Bài toán cấp hai Dirichlet $u''(x) = f(x)$, $u(0) = u(L) = 0$ có hàm Green: $$G(x, t) = \begin{cases} \frac{1}{L} x (t - L), & 0 \le x \le t \le L \ \frac{1}{L} t (x - L), & 0 \le t \le x \le L \end{cases}$$
- Bài toán cấp hai hỗn hợp $u''(x) = f(x)$, $u(0) = 0$, $u'(L) = 0$ có hàm Green $G(x, t) = -x$ khi $x \le t$ và $t - x$ khi $t \le x$.
- Bài toán cấp hai với điều kiện biên loại ba $au(0) - bu'(0) = 0$, $cu(1) + du'(1) = 0$ (với $\rho = ad + bc + ac \ne 0$) có hàm Green: $$G(x, t) = \begin{cases} \frac{1}{\rho} (ct - c - d)(b + ax), & 0 \le x \le t \le 1 \ \frac{1}{\rho} (b + at)(cx - c - d), & 0 \le t \le x \le 1 \end{cases}$$
- Bài toán cấp bốn $u^{(4)}(x) = f(x)$ với điều kiện biên tựa đơn giản $u(0) = u''(0) = u(L) = u''(L) = 0$, điều kiện biên Dirichlet $u(a) = u'(a) = u(b) = u'(b) = 0$, và điều kiện biên tổ hợp $u(0) = 0, u'(1) = 0, au''(0) - bu'''(0) = 0, cu''(1) + du'''(1) = 0$.
- Tổng hợp các công thức đạo hàm số và tích phân số:
- Công thức sai phân xấp xỉ đạo hàm cấp một với sai số $O(h^2)$ và $O(h^4)$ dựa trên nội suy Lagrange 3 điểm và 5 điểm.
- Công thức tích phân hình thang (sai số $O(h^2)$), công thức Simpson 1 chiều (sai số $O(h^4)$) và công thức Simpson 2 chiều trên miền chữ nhật $\Omega$ (sai số $O(h^4) + O(k^4)$).
- Thiết lập lược đồ sai phân độ chính xác cấp bốn cho phương trình Poisson hai chiều $\Delta u = -f$.
- Trình bày thuật toán giải hệ phương trình đại số tuyến tính dạng ba điểm:
- Phương pháp truy đuổi phải và truy đuổi từ hai phía cho hệ phương trình vô hướng ma trận ba đường chéo trội (chứng minh tính khả thi và ổn định theo Bổ đề 1.1).
- Thuật toán rút gọn hoàn toàn (cyclic reduction) giải hệ phương trình véc tơ ba điểm với $N = 2^n$ bước.
Chương 2: Sự tồn tại duy nhất nghiệm và phương pháp lặp giải bài toán biên cho phương trình vi phân thường phi tuyến cấp bốn
Chương 2 trình bày kết quả nghiên cứu định tính và thuật toán số cho năm bài toán biên đối với phương trình vi phân thường cấp bốn địa phương và không địa phương. Nội dung được chia thành hai phần chính: bài toán địa phương (gồm ba bài toán với điều kiện biên tổ hợp, điều kiện biên Dirichlet, và điều kiện biên phi tuyến) và bài toán không địa phương (gồm hai bài toán với điều kiện biên gối - tựa đơn giản và điều kiện biên phi tuyến).
Trọng tâm phân tích được thể hiện qua bài toán biên địa phương với điều kiện biên tổ hợp mô tả độ võng của dầm: $$u^{(4)}(x) = f(x, u(x), u'(x), u''(x), u'''(x)), \quad 0 < x < 1$$ thỏa mãn điều kiện biên: $$u(0) = 0, \quad u'(1) = 0, \quad au''(0) - bu'''(0) = 0, \quad cu''(1) + du'''(1) = 0$$ với các giả thiết $a, b, c, d \ge 0$, $\rho := ad + bc + ac > 0$, hàm $f: [0, 1] \times \mathbb{R}^4 \to \mathbb{R}$ liên tục.
Phương pháp nghiên cứu đưa bài toán về tìm điểm bất động của toán tử phi tuyến: $$(A\phi)(x) = f(x, u(x), u'(x), u''(x), u'''(x))$$ trong đó $u(x)$ là nghiệm duy nhất của phương trình vi phân tuyến tính $u^{(4)}(x) = \phi(x)$ với cùng điều kiện biên. Bằng cách hạ bậc thông qua biến phụ $v(x) = u''(x)$, bài toán được phân rã thành hai bài toán biên cấp hai liên tiếp:
- $v''(x) = \phi(x)$ với $av(0) - bv'(0) = 0$, $cv(1) + dv'(1) = 0$.
- $u''(x) = v(x)$ với $u(0) = 0$, $u'(1) = 0$.
Biểu diễn nghiệm được xác định qua các hàm Green tương ứng, từ đó tính toán các hằng số chặn tích phân: $$\rho_1 = \max_{0 \le x \le 1} \int_0^1 G(x, t) dt = \frac{1}{24} + \frac{2ad + bc + 6bd}{12\rho}$$ $$\rho_2 = \max_{0 \le x \le 1} \int_0^1 |G'x(x, t)| dt = \frac{1}{12} + \frac{ad + bc + 4bd}{4\rho}$$ $$\rho_3 = \max{0 \le x \le 1} \int_0^1 |G_1(x, t)| dt = \frac{ac + 2bc + 2bd}{2\rho}$$ $$\rho_4 = \max_{0 \le x \le 1} \int_0^1 |(G_1)'_x(x, t)| dt = \frac{ac + bc + bd}{\rho}.$$
Luận án định nghĩa miền đóng: $$D_M = \left{ (x, u, y, v, z) \in [0, 1] \times \mathbb{R}^4 ;\middle|; |u| \le \rho_1 M, |y| \le \rho_2 M, |v| \le \rho_3 M, |z| \le \rho_4 M \right}.$$
Các định lý định tính chính của chương gồm:
- Bổ đề 2.1 và Định lý 2.1: Giả sử tồn tại các hằng số $M > 0$ và $K_1, K_2, K_3, K_4 \ge 0$ sao cho $|f(x, u, y, v, z)| \le M$ trên miền $D_M$; đồng thời hàm $f$ thỏa mãn điều kiện Lipschitz theo bốn biến hàm trên $D_M$ với các hệ số $K_i$. Nếu: $$q = K_1 \rho_1 + K_2 \rho_2 + K_3 \rho_3 + K_4 \rho_4 < 1$$ thì toán tử $A$ là ánh xạ co trên quả cầu đóng $B[O, M] \subset C[0, 1]$. Bài toán biên có nghiệm duy nhất $u(x)$ thỏa mãn các đánh giá $|u| \le \rho_1 M$, $|u'| \le \rho_2 M$, $|u''| \le \rho_3 M$, $|u'''| \le \rho_4 M$.
- Định lý 2.2 (Tính dương của nghiệm): Nếu trên miền $D_M$, hàm $f$ thỏa mãn thêm điều kiện không âm $0 \le f(x, u, y, v, z) \le M$ và các điều kiện co của Bổ đề 2.1 được bảo đảm, thì nghiệm duy nhất của bài toán là nghiệm không âm ($u(x) \ge 0, \forall x \in [0, 1]$).
- Các kết quả tương tự được thiết lập cho trường hợp điều kiện biên Dirichlet, điều kiện biên phi tuyến, và hai bài toán phương trình cấp bốn không địa phương loại Kirchhoff.
Thực nghiệm số của Chương 2 được minh họa qua hệ thống từ Ví dụ 2.1 đến Ví dụ 2.25. Quá trình lặp kiểm tra sai số qua đại lượng $e(K)$ và tỷ số hội tụ thực nghiệm $r(K)$ cho cả hai trường hợp: đã biết nghiệm đúng giải tích và chưa biết trước nghiệm đúng (sử dụng xấp xỉ đầu $v_0 = 0$ hoặc $u_0 = 0, v_0 = 0$).
| Ví dụ số trong Chương 2 | Điều kiện xấp xỉ đầu | Nội dung kiểm tra | Dữ liệu bảng/hình tương ứng |
|---|---|---|---|
| Ví dụ 2.1 - 2.8 | $v_0$ xác định | Đánh giá sai số $e(K)$, tỷ số $r(K)$, nghiệm xấp xỉ | Bảng 2.1, 2.2; Hình 2.1 - 2.8 |
| Ví dụ 2.13 | $v_0 = 0$ | Kiểm tra sự hội tụ | Bảng 2.3, 2.4 |
| Ví dụ 2.14 - 2.16 | $v_0 = 0$ | Đánh giá sự hội tụ và đồ thị nghiệm xấp xỉ | Bảng 2.5 - 2.7; Hình 2.11 - 2.17 |
| Ví dụ 2.21 - 2.25 | $u_0 = 0, v_0 = 0$ | Đánh giá sự hội tụ phương trình không địa phương | Bảng 2.9 - 2.13; Hình 2.18 - 2.22 |
Chương 3: Sự tồn tại duy nhất nghiệm và phương pháp lặp giải bài toán biên phi tuyến cho phương trình đạo hàm riêng cấp bốn
Chương 3 phát triển kỹ thuật toán tử trung gian cho hai bài toán biên phi tuyến đối với phương trình đạo hàm riêng cấp bốn:
-
Bài toán biên phi tuyến cho phương trình song điều hòa: Khảo sát phương trình dạng: $$\Delta^2 u = f(x, u, \Delta u), \quad x \in \Omega$$ với các điều kiện biên tương ứng trên $\partial\Omega$. Luận án thiết lập định lý về sự tồn tại và duy nhất nghiệm, đồng thời phân tích tính dương của nghiệm khi hàm $f$ nhận giá trị không âm. Về mặt tính toán số, luận án đề xuất phương pháp lặp phân rã bậc hai ký hiệu là DIM2. Phương pháp này đưa phương trình song điều hòa về hệ hai bài toán Poisson cấp hai liên tiếp tại mỗi bước lặp, áp dụng lược đồ sai phân cấp bốn kết hợp phương pháp rút gọn hoàn toàn giải hệ đại số.
-
Bài toán biên phi tuyến cho phương trình song điều hòa loại Kirchhoff: Khảo sát bài toán có chứa thành phần tích phân phi tuyến không địa phương: $$\Delta^2 u = M\left( \int_\Omega |\nabla u|^2 dx \right) \Delta u + f(x, u), \quad x \in \Omega$$ dưới các điều kiện biên xác định. Luận án chứng minh sự tồn tại duy nhất nghiệm bằng định lý ánh xạ co Banach trên miền giới nội thích hợp và xây dựng thuật toán lặp số kết hợp công thức tích phân Simpson hai lớp.
-
Thực nghiệm tính toán và so sánh số liệu: Các thử nghiệm số trong Chương 3 được thực hiện trên lưới đều kích thước $65 \times 65$ nút (Bảng 3.1, Bảng 3.2). Kết quả tính toán đánh giá sai số $e(m)$ và tỷ số $r(m)$ (tương ứng từ Ví dụ 3.1 đến 3.6). Luận án đối chiếu trực tiếp hiệu năng của phương pháp DIM2 với phương pháp số nhanh nhất của tác giả Wang:
- Hình 3.1, 3.2, 3.3, 3.4 thể hiện đồ thị so sánh sai số $e(m)$ giữa phương pháp DIM2 (bên trái) và phương pháp nhanh nhất của Wang (bên phải) trong các ví dụ từ 3.1 đến 3.4, chứng minh tốc độ giảm sai số của phương pháp đề xuất.
- Hình 3.5 đến 3.15 biểu diễn đồ thị sai số $e(m)$, tỷ số $r(m)$ và hình dáng bề mặt của nghiệm xấp xỉ hai chiều trên miền tính toán.
Kết quả và những đóng góp mới
Các đóng góp mới của luận án được tác giả tự xác định bao gồm hai phương diện lý luận và thực tiễn tính toán:
-
Đóng góp về mặt lý luận:
- Luận án đề xuất một phương pháp tiếp cận chuyển các bài toán biên phi tuyến cấp bốn (cả phương trình vi phân thường và phương trình đạo hàm riêng, dạng địa phương và không địa phương) về phương trình toán tử đối với hàm vế phải hoặc hàm trung gian.
- Thiết lập được các điều kiện đủ đơn giản, dễ kiểm tra (dựa trên tính Lipschitz cục bộ và hằng số tích phân của hàm Green trên miền giới nội $D_M$) để bảo đảm sự tồn tại và duy nhất nghiệm cho:
- Năm bài toán biên cho phương trình vi phân thường phi tuyến cấp bốn với các điều kiện biên tổ hợp, Dirichlet, gối - tựa đơn giản và phi tuyến.
- Hai bài toán biên phi tuyến cho phương trình song điều hòa và phương trình song điều hòa loại Kirchhoff.
- Chứng minh được tính dương của nghiệm đối với ba bài toán: bài toán vi phân thường cấp bốn với điều kiện biên Dirichlet, bài toán vi phân thường cấp bốn với điều kiện biên tổ hợp, và bài toán biên cho phương trình song điều hòa.
- Khắc phục được các đòi hỏi hạn chế của các phương pháp truyền thống: không cần áp đặt điều kiện tăng trưởng tại vô cùng, không cần điều kiện Nagumo phức tạp, và không cần giả thiết về sự tồn tại tiên nghiệm của cặp nghiệm trên - nghiệm dưới.
- Đề xuất các phương pháp lặp tìm nghiệm ở mức độ bài toán liên tục và chứng minh sự hội tụ của phương pháp lặp với tốc độ cấp số nhân.
-
Đóng góp về mặt thực tiễn và tính toán số:
- Xây dựng thuật toán tính toán cụ thể kết hợp các lược đồ sai phân có độ chính xác cao: công thức đạo hàm số nội suy Lagrange bậc cao, công thức tích phân số Simpson một chiều và hai chiều (độ chính xác cấp bốn), và lược đồ sai phân cấp bốn cho phương trình Poisson.
- Ứng dụng các thuật toán giải hệ phương trình lưới hiệu quả cao: phương pháp truy đuổi phải và truy đuổi hai phía cho hệ ba điểm (chi phí khoảng $8N$ phép toán), phương pháp rút gọn hoàn toàn (cyclic reduction) cho hệ véc tơ ba điểm.
- Cung cấp hệ thống ví dụ tính toán số phong phú trên phần mềm MATLAB 7.0 cho cả bài toán đã biết nghiệm đúng giải tích và chưa biết nghiệm giải tích; đưa ra các so sánh định lượng trực tiếp minh họa ưu thế về tốc độ hội tụ của phương pháp đề xuất (phương pháp DIM2) so với phương pháp nhanh nhất của Wang.
Hạn chế và hướng nghiên cứu tiếp
Hạn chế của luận án:
- Phạm vi nghiên cứu giới hạn trong năm mô hình bài toán vi phân thường cấp bốn một chiều trên đoạn $[0, 1]$ (hoặc $[0, L]$, $[a, b]$) và hai mô hình bài toán đạo hàm riêng cấp bốn trên miền giới nội $\Omega$.
- Các điều kiện bảo đảm tính co của toán tử đòi hỏi hệ số Lipschitz $K_i$ của hàm vế phải và bán kính $M$ của miền giới nội phải thỏa mãn bất đẳng thức chặt $q = \sum K_i \rho_i < 1$. Nếu bài toán có tính phi tuyến mạnh vượt khỏi ngưỡng co này, phương pháp chưa áp dụng được.
- Trong tính toán số cho phương trình đạo hàm riêng cấp bốn (Chương 3), thực nghiệm số mới được tiến hành trên miền hình chữ nhật phẳng hai chiều với lưới đều sai phân kích thước $65 \times 65$ nút; chưa khảo sát trên các miền có biên cong tùy ý hoặc không gian ba chiều.
Hướng nghiên cứu tiếp được xác định từ phạm vi luận án:
- Mở rộng phương pháp toán tử trung gian cho các bài toán biên phi tuyến cấp bốn trên các miền hình học phức tạp, phi chữ nhật thông qua việc kết hợp với các kỹ thuật biến đổi tọa độ hoặc phương pháp phần tử hữu hạn.
- Nghiên cứu áp dụng cho các hệ phương trình vi phân phi tuyến bậc cao mô tả các mô hình vật lý phức tạp hơn trong cơ học kết cấu, chẳng hạn như hệ phương trình bản Von Kármán hoặc các bài toán dầm phi tuyến chịu tải trọng động phụ thuộc thời gian.
- Tối ưu hóa thuật toán tính toán song song khi giải hệ phương trình lưới trên các cấu trúc lưới có số nút lớn hơn.
Giá trị tham khảo
Luận án là tài liệu tham khảo hữu ích cho các nhóm đối tượng nghiên cứu và đào tạo sau:
- Nghiên cứu sinh, học viên cao học và giảng viên thuộc chuyên ngành Toán ứng dụng, Toán giải tích, Giải tích số và Phương pháp tính. Luận án cung cấp một quy trình mẫu mực về việc kết hợp giữa phân tích định tính (chứng minh tồn tại, duy nhất, nghiệm dương qua lý thuyết điểm bất động) và xây dựng thuật toán số giải gần đúng kèm đánh giá sai số.
- Các nhà nghiên cứu và kỹ sư trong lĩnh vực Cơ học tính toán, Cơ học kết cấu công trình, đặc biệt là những người làm việc với các mô hình uốn dầm, tấm và bản phẳng trên nền đàn hồi.
- Phần đáng tham khảo nhất trong luận án gồm:
- Chương 1: Quy trình tính toán tường minh hàm Green cho các bài toán biên vi phân cấp hai và cấp bốn với các điều kiện biên hỗn hợp và tổ hợp; thuật toán truy đuổi hai phía và thuật toán rút gọn hoàn toàn giải hệ phương trình sai phân.
- Mục 2.1.1 của Chương 2: Kỹ thuật phân rã phương trình cấp bốn thành hai bài toán cấp hai liên tiếp thông qua biến phụ $v = u''$, cách xây dựng miền giới nội $D_M$ và phương pháp tính toán chính xác các hằng số cực đại tích phân $\rho_1, \rho_2, \rho_3, \rho_4$.
- Chương 3: Thuật toán lặp DIM2 giải phương trình song điều hòa và song điều hòa Kirchhoff trên lưới sai phân hai chiều kết hợp công thức Simpson hai lớp.
Câu hỏi thường gặp
-
Phương trình vi phân cấp bốn không địa phương trong luận án được hiểu như thế nào và khác gì so với phương trình địa phương? Theo định nghĩa trong luận án, phương trình vi phân cấp bốn địa phương là phương trình mà giá trị của toán tử vi phân tại một điểm chỉ phụ thuộc vào giá trị của hàm và các đạo hàm tại chính điểm đó. Ngược lại, phương trình vi phân cấp bốn không địa phương (hay phương trình loại Kirchhoff) là phương trình có chứa thành phần tích phân trên toàn miền, ví dụ số hạng chứa tích phân $M\left( c \int_0^1 |u'(s)|^2 ds \right) u''(x)$ hoặc $M\left( \int_\Omega |\nabla u|^2 dx \right) \Delta u$.
-
Vì sao phương pháp toán tử trung gian của luận án khắc phục được hạn chế của phương pháp nghiệm trên - nghiệm dưới và điều kiện Nagumo? Phương pháp nghiệm trên - nghiệm dưới đòi hỏi phải tìm được trước một cặp hàm $\alpha(x) \le \beta(x)$ thỏa mãn các bất đẳng thức vi phân biên và đòi hỏi hàm vế phải phải thỏa mãn điều kiện Nagumo (chặn tốc độ tăng trưởng của đạo hàm). Luận án không tìm nghiệm trên/dưới mà đưa trực tiếp bài toán về phương trình toán tử đối với hàm vế phải $\phi$ trên một miền đóng giới nội $D_M$. Bằng cách thu hẹp miền khảo sát vào $D_M$, hàm phi tuyến chỉ cần thỏa mãn tính liên tục Lipschitz cục bộ, do đó loại bỏ hoàn toàn việc kiểm tra điều kiện Nagumo và điều kiện tăng trưởng tại vô cùng.
-
Các hằng số $\rho_1, \rho_2, \rho_3, \rho_4$ trong Bổ đề 2.1 được xác định từ đâu và có vai trò gì? Các hằng số $\rho_1, \rho_2, \rho_3, \rho_4$ là các giá trị cực đại trên $[0, 1]$ của tích phân các chuẩn của hàm Green $G(x, t)$, $G'_x(x, t)$, $G_1(x, t)$ và $(G_1)'_x(x, t)$ thu được khi phân rã bài toán cấp bốn thành hai bài toán cấp hai liên tiếp. Chúng đóng vai trò xác định kích thước biên của miền đóng $D_M$ ($|u| \le \rho_1 M$, $|u'| \le \rho_2 M$, $|u''| \le \rho_3 M$, $|u'''| \le \rho_4 M$) và cấu thành trực tiếp hệ số co $q = K_1 \rho_1 + K_2 \rho_2 + K_3 \rho_3 + K_4 \rho_4 < 1$ trong Định lý điểm bất động Banach.
-
Luận án đã sử dụng những thuật toán giải hệ phương trình đại số nào sau khi rời rạc hóa bài toán? Luận án sử dụng hai thuật toán đại số chuyên biệt:
- Phương pháp truy đuổi phải và truy đuổi từ hai phía để giải hệ phương trình vô hướng ba điểm có ma trận ba đường chéo trội (thu được khi rời rạc hóa phương trình vi phân thường cấp hai), với độ phức tạp tính toán thấp, xấp xỉ $8N$ phép toán.
- Phương pháp rút gọn hoàn toàn (cyclic reduction) với $n-1$ bước khử liên tiếp để giải hệ phương trình véc tơ ba điểm $-Y_{j-1} + C Y_j - Y_{j+1} = F_j$ thu được khi rời rạc hóa phương trình elliptic cấp hai trên lưới chữ nhật $N = 2^n$ bước.
- Thực nghiệm tính toán số trong luận án được kiểm chứng trên nền tảng nào và phương pháp đề xuất có ưu thế gì so với nghiên cứu trước? Các thực nghiệm tính toán số được lập trình trên môi trường MATLAB 7.0, vận hành trên máy tính cá nhân sử dụng CPU Intel Core i3, RAM 4GB. Trong Chương 3, phương pháp lặp DIM2 của luận án khi thử nghiệm trên lưới đều $65 \times 65$ nút đã được so sánh trực tiếp với phương pháp số nhanh nhất của Wang (2006). Kết quả biểu diễn qua đồ thị sai số $e(m)$ tại các hình 3.1 đến 3.4 cho thấy phương pháp DIM2 có tốc độ giảm sai số nhanh hơn rõ rệt qua các bước lặp.
Kết luận
Luận án của nghiên cứu sinh Nguyễn Thanh Hường đã giải quyết bài toán biên phi tuyến cho phương trình vi phân thường và phương trình đạo hàm riêng cấp bốn bằng cách tiếp cận đưa về phương trình toán tử đối với hàm vế phải hoặc hàm trung gian trên miền giới nội. Công trình đã chứng minh sự tồn tại, tính duy nhất và tính dương của nghiệm dưới các điều kiện Lipschitz cục bộ dễ kiểm tra, đồng thời xây dựng các phương pháp lặp giải gần đúng hội tụ cấp số nhân. Hệ thống tính toán thực nghiệm trên MATLAB kết hợp các lược đồ sai phân cấp cao đã kiểm chứng tính đúng đắn của lý thuyết và thể hiện ưu thế tính toán so với một số phương pháp đã công bố trước đó.