Tổng quan luận án

Luận án Phó Tiến sĩ Toán Lý với đề tài "Đạo hàm trung bình và các phương pháp phối hợp để giải một số bài toán biên" được thực hiện bởi nghiên cứu sinh tại Khoa Toán - Cơ - Tin học, Trường Đại học Tổng hợp Hà Nội vào năm 1993, dưới sự hướng dẫn khoa học của thầy Phan Văn Hạp. Luận án thuộc chuyên ngành Toán học tính toán, mã số chuyên ngành 1-01-07.

Trong thực tế tính toán cơ học, vật lý và kỹ thuật môi trường, nhiều bài toán biên xuất hiện các yếu tố "không trơn" như: nghiệm có bước nhảy hoặc gián đoạn, các hệ số trong phương trình vi phân bị gián đoạn, miền tính toán có biên không trơn (biên kỳ dị, có điểm góc), hoặc các quá trình không dừng diễn ra trên miền phức tạp. Khi giải các bài toán này, các công cụ vi tích phân cổ điển và đạo hàm theo nghĩa Gâteaux thường gặp bế tắc do đòi hỏi tính khả vi và độ trơn cao của hàm số. Khoảng trống nghiên cứu được luận án tập trung giải quyết là việc xây dựng, phát triển công cụ giải tích phù hợp – cụ thể là khái niệm đạo hàm trung bình tích phân – kết hợp đồng bộ với các phương pháp giải số và giải tích (phương pháp sai phân, phép biến đổi Laplace, phương pháp phần tử biên, phương pháp đường thẳng) nhằm thiết lập nghiệm giải tích chính xác hoặc nghiệm xấp xỉ số có độ tin cậy cao cho các lớp bài toán biên không trơn.

Mục tiêu nghiên cứu của luận án được xác định cụ thể qua các nhiệm vụ sau:

  1. Nghiên cứu cơ sở lý thuyết của đạo hàm trung bình tích phân, xác lập mối quan hệ tương hỗ giữa đạo hàm trung bình với các toán tử biến đổi tích phân, đặc biệt là phép biến đổi Laplace một chiều và hai chiều.
  2. Ứng dụng đạo hàm trung bình tích phân để giải đúng các bài toán biên mô tả quá trình khuếch tán và chuyển dịch vật chất từ nguồn ô nhiễm điểm có nghiệm không trơn và nghiệm gián đoạn.
  3. Xây dựng các lược đồ sai phân mới dựa trên đạo hàm trung bình để giải lớp bài toán biên một chiều có hệ số gián đoạn (loại 1), đồng thời chứng minh chặt chẽ tính xấp xỉ, tính ổn định và sự hội tụ của nghiệm sai phân về nghiệm đúng.
  4. Kết hợp đạo hàm trung bình với phương pháp phần tử biên (Boundary Element Method - BEM) dựa trên nguyên lý cân bằng sai số để giải bài toán biên elliptic (phương trình Laplace, phương trình Poisson) với điều kiện biên hỗn hợp trên miền có biên kỳ dị (chứa điểm góc).
  5. Mở rộng phát triển các phương pháp phối hợp (biến đổi Laplace phối hợp phần tử biên, phương pháp đường thẳng) để giải bài toán biên phương trình truyền dẫn - parabolic không dừng trên miền có biên kỳ dị.

Đối tượng nghiên cứu của luận án là các phương trình vi phân đạo hàm riêng và phương trình vi phân thường kèm theo các điều kiện biên loại Dirichlet, Neumann hoặc hỗn hợp (Robin), đặc trưng bởi các hàm số gián đoạn, nghiệm không trơn và miền hình học có biên góc cạnh. Phạm vi nghiên cứu bao gồm các bài toán một chiều, hai chiều và ba chiều trong không gian Euclid, cùng các quá trình không dừng theo biến thời gian $t > 0$.

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

Tổng quan các công trình nghiên cứu trước đây được tác giả tổng hợp và phân tích qua các nhánh chính sau:

  • Nghiên cứu về bài toán ô nhiễm môi trường và nghiệm không trơn: Tác giả Mactruc [3] đã mô hình hóa các bài toán lan truyền chất ô nhiễm môi trường phát ra từ nguồn ô nhiễm có kích thước nhỏ so với phạm vi ảnh hưởng, dẫn đến các phương trình vi phân chứa hàm xung Dirac $\delta(x - x_0)$ và hàm bước nhảy đơn vị $u^*(x - x_0)$. Các nghiệm của mô hình này có đặc tính gián đoạn loại 1 hoặc có điểm góc (không trơn), làm mất đi tính khả vi Gâteaux thông thường.
  • Nghiên cứu về phương pháp số cho bài toán có hệ số không trơn: Các tác giả Ladujenxkaia [5], Ianhenko [8], và Xamarxki [7] đã phát triển nhiều phương pháp xây dựng lược đồ sai phân dựa trên các nguyên lý biến phân và phương pháp thể tích hữu hạn/cân bằng thông lượng để xử lý bài toán có hệ số không trơn hoặc gián đoạn.
  • Nghiên cứu về bài toán trên miền biên kỳ dị: Các công trình của Godunov [2], Ianhenko [3], Xamarxki [7] đã giải quyết các bài toán biên với điều kiện Dirichlet hoặc Neumann trên miền không trơn bằng các kỹ thuật như khai triển tổng, phương pháp phân rã miền.
  • Cơ sở phương pháp cân bằng sai số và đạo hàm mở rộng: Các khái niệm nghiệm xấp xỉ và phương pháp cân bằng sai số (Weighted Residuals) được tiếp thu từ tài liệu [1], kết hợp với lý thuyết về đạo hàm trung bình tích phân trong không gian hàm từ tài liệu [13], cùng các kỹ thuật chọn tham số biến đổi ngược Laplace từ tài liệu [9].

Vị trí học thuật của luận án: Luận án không đi theo lối mòn của các phương pháp giải tích cổ điển vốn đòi hỏi hàm trơn, cũng không hoàn toàn phụ thuộc vào các phép xấp xỉ sai phân biến phân truyền thống. Bằng việc sử dụng đạo hàm trung bình tích phân làm công cụ trung tâm, tác giả đã tạo ra một cầu nối nhất quán giữa giải tích hàm phi cổ điển với các thuật toán số (lược đồ sai phân, phần tử biên, sai phân thời gian). Cách tiếp cận này cho phép khắc phục hiện tượng mất đạo hàm tại các điểm gián đoạn và điểm góc của biên, mang lại lời giải chính xác cho các bài toán truyền chất suy biến và thiết lập các hệ đại số tuyến tính đối xứng, xác định dương cho các bài toán biên kỳ dị phức tạp.

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

Luận án xây dựng trên hệ thống lý thuyết toán học giải tích và các phương pháp tính số chặt chẽ:

Khái niệm Đạo hàm trung bình tích phân

Xét ánh xạ $F: X \to Y$, trong đó $X, Y$ là các không gian Banach (B-không gian), $x_0 \in X$. Nếu với mọi $h \in X$, tồn tại giới hạn: $$\lim_{\delta \to 0} \frac{1}{2\delta} \int_{-\delta}^{\delta} \frac{F(x_0 + th) - F(x_0)}{t} dt = \delta F(x_0, h)$$ thì đại lượng $\delta F(x_0, h)$ được gọi là biến phân trung bình (theo tích phân) thứ nhất của hàm $F$ tại điểm $x_0$.

Nếu tại $x_0$ tồn tại toán tử tuyến tính giới nội $A$ sao cho $\delta F(x_0, h) = Ah$, thì hàm $F$ được gọi là khả vi trung bình và có đạo hàm trung bình theo tích phân tại $x_0$, ký hiệu là $F'(x_0, h) = Ah$. Khái niệm này sở hữu các tính chất nền tảng:

  • Tính tuyến tính: $(F + G)'(x_0) = F'(x_0) + G'(x_0)$ và $(\alpha F)'(x_0) = \alpha F'(x_0)$.
  • Tính bao quát hơn đạo hàm Gâteaux: Một hàm số có thể không tồn tại đạo hàm Gâteaux tại một điểm nhưng vẫn tồn tại đạo hàm trung bình tích phân (ví dụ hàm $F(x) = |x|$ tại $x_0 = 0$ không có đạo hàm Gâteaux nhưng có đạo hàm trung bình $F'(0) = 0$).
  • Xử lý điểm góc và bước nhảy: Hàm có điểm góc hoặc không liên tục tại $x_0$ vẫn có thể có đạo hàm trung bình. Cụ thể, nếu $F$ có đạo hàm hữu hạn hai phía tại $x_0$, thì đạo hàm trung bình được xác định bằng trung bình cộng của đạo hàm trái và đạo hàm phải: $$F'(x_0) = \frac{F'(x_0 + 0) + F'(x_0 - 0)}{2}$$

Phương pháp cân bằng sai số (Weighted Residual Method)

Để tìm nghiệm xấp xỉ của phương trình toán tử vi phân $L(U) = b$ trong miền $\Omega$ với các điều kiện biên $S(U) = s$ trên $\Gamma_1$ và $G(U) = g$ trên $\Gamma_2$, nghiệm xấp xỉ được biểu diễn dưới dạng tổ hợp $U^* = \sum_{i=1}^n \alpha_i \phi_i + \alpha_0$ trên hệ hàm cơ sở đầy đủ ${\phi_i}$. Bằng cách sử dụng các hàm trọng lượng $W$, tích phân sai số có trọng số được triệt tiêu trên toàn miền: $$\int_\Omega [L(U^*) - b] W d\Omega = 0$$

Các phép biến đổi toán tử tích phân

  • Biến đổi Laplace một chiều: $F(s) = \mathcal{L}[f(t), s] = \int_0^{+\infty} f(t) e^{-st} dt$. Đối với hàm bị giới nội có đạo hàm trung bình và chứa các điểm gián đoạn hữu hạn $t_1, t_2, \dots$, công thức biến đổi của đạo hàm trung bình được xác lập: $$\mathcal{L}[f'(t), s] = sF(s) - f(0+0) - \sum_{i} \exp(-t_i s) [f(t_i + 0) - f(t_i - 0)]$$
  • Biến đổi Laplace hai chiều: $\mathcal{L}2[f(t), s] = \int{-\infty}^{+\infty} f(t) e^{-st} dt$, liên hệ với biến đổi một chiều qua biểu thức $\mathcal{L}_2[f(t), s] = \mathcal{L}[f(t), s] + \mathcal{L}[f(-t), -s]$.

Phương pháp phần tử biên và phương pháp đường thẳng

  • Phương pháp phần tử biên (BEM): Sử nghiệm cơ bản của phương trình vi phân kết hợp tích phân từng phần hai lần theo công thức Green để đưa phương trình vi phân trong miền $\Omega$ về phương trình tích phân thuần túy trên biên $\Gamma$, cho phép xử lý biên có điểm góc bằng cách chia biên thành $N$ phần tử biên $\Gamma_j$.
  • Phương pháp đường thẳng (Method of Lines): Rời rạc hóa đạo hàm thời gian bằng sai phân hữu hạn để biến đổi phương trình parabolic không dừng thành một chuỗi các bài toán elliptic dừng tại từng bước thời gian, giải tiếp bằng phương pháp phần tử biên.

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

Chương 1: Mở đầu

Chương này dành cho việc trình bày hệ thống các định nghĩa, tính chất của đạo hàm trung bình tích phân và các khái niệm bổ trợ toán học phục vụ cho toàn bộ các chương tiếp theo của luận án.

Tác giả đưa ra định nghĩa giải tích chuẩn xác của đạo hàm trung bình tích phân, chứng minh chi tiết định lý về đạo hàm trung bình của hàm có đạo hàm hai phía tại điểm kỳ dị $x_0$. Chương 1 trình bày phương pháp nghiệm xấp xỉ, nguyên lý cân bằng sai số, và khảo sát tường minh các phép biến đổi toán tử gồm biến đổi Laplace một chiều, biến đổi Laplace hai chiều, phép biến đổi Laplace ngược, cùng điều kiện hội tụ tuyệt đối của các tích phân. Định lý thiết lập mối quan hệ hữu cơ giữa phép biến đổi Laplace và đạo hàm trung bình được thiết lập, chứng minh tính đúng đắn khi áp dụng cho các hàm gián đoạn có bước nhảy hữu hạn.

Chương 2: Các bài toán biên có nghiệm không trơn

Chương 2 tập trung giải quyết các bài toán biên xuất phát từ mô hình toán học của Mactruc [3] mô tả quá trình khuếch tán và di chuyển của thực thể trong không khí với vận tốc gió và nguồn thải điểm.

  • Bài toán biên phức tạp có nghiệm gián đoạn: Phương trình vi phân có dạng: $$\nu \frac{d^2 \Phi(x)}{dx^2} - \alpha \frac{d\Phi(x)}{dx} + \beta [\Phi(x) - \chi u^(x - x_0)] = \chi_1 + Q \delta(x - x_0)$$ với điều kiện biên $\lim_{x \to -\infty} \Phi(x) = 0$, $\lim_{x \to +\infty} \Phi(x) = \chi$ ($\chi > 0$), trong đó $\nu > 0, \alpha, \beta, Q$ là các hằng số, $u^(x-x_0)$ là hàm bước nhảy Heaviside và $\delta(x-x_0)$ là hàm Dirac. Nghiệm của bài toán bị gián đoạn loại 1 tại điểm nguồn $x = x_0$. Bằng cách coi các đạo hàm trong phương trình là đạo hàm trung bình tích phân và áp dụng biến đổi Laplace hai chiều, tác giả thiết lập Định lý 2 khẳng định nghiệm tìm được bằng phương pháp phối hợp hoàn toàn trùng khớp với nghiệm giải tích chính xác: $$\Phi(x) = \begin{cases} \frac{Q}{\nu \sqrt{\frac{\alpha^2}{\nu^2} + \frac{4\beta}{\nu}}} \exp\left[\left(\frac{\alpha}{2\nu} - \sqrt{\frac{\alpha^2}{4\nu^2} + \frac{\beta}{\nu}}\right)(x - x_0)\right] + \chi, & x > x_0 \ \frac{Q}{\nu \sqrt{\frac{\alpha^2}{\nu^2} + \frac{4\beta}{\nu}}} \exp\left[\left(\frac{\alpha}{2\nu} + \sqrt{\frac{\alpha^2}{4\nu^2} + \frac{\beta}{\nu}}\right)(x - x_0)\right], & x < x_0 \end{cases}$$

  • Bài toán biên dạng phức tạp với nghiệm không trơn: Phương trình $\nu \frac{d^2 \Phi(x)}{dx^2} - \alpha \frac{d\Phi(x)}{dx} + \beta \Phi(x) = Q \delta(x - x_0)$ với điều kiện $\lim_{x \to \pm\infty} \Phi(x) = 0$. Mệnh đề 1 khẳng định nghiệm thu được từ phương pháp phối hợp đạo hàm trung bình và biến đổi Laplace trùng với nghiệm giải tích (nghiệm liên tục nhưng có điểm góc tại $x_0$).

  • Bài toán biên dạng đơn giản có nghiệm không trơn: Khi vận tốc khối không khí triệt tiêu ($\alpha = 0$), phương trình có dạng $\nu \frac{d^2 \Phi(x)}{dx^2} + \beta \Phi(x) = Q \delta(x - x_0)$. Mệnh đề 2 chứng minh phương pháp phối hợp cho nghiệm chính xác mô tả sự khuếch tán đối xứng hai phía quanh điểm $x_0$.

Chương 3: Các bài toán biên có hệ số không trơn

Chương 3 xây dựng các lược đồ sai phân hữu hạn mới cho các bài toán biên một chiều có các hệ số $p(x), q(x), r(x)$ gián đoạn loại 1, ứng dụng khái niệm đạo hàm trung bình.

  • Bài toán một chiều dạng đơn giản: $$-\frac{d}{dx}\left[p(x)\frac{d\Phi(x)}{dx}\right] + q(x)\Phi(x) = f(x), \quad 0 < x < 1$$ với $\Phi(0) = \Phi(1) = 0$. Chia đoạn $[0, 1]$ thành lưới không đều với các nút ${x_k}$ và các điểm giữa ${x_{k+1/2}}$. Sử dụng hệ hàm cơ sở ${\omega_k(x)}$ dạng chóp tuyến tính từng đoạn, tích phân từng phần hai lần phương trình có trọng số theo đạo hàm trung bình dẫn đến hệ đại số sai phân dạng ma trận $A\Phi = F$, trong đó $A$ là ma trận ba đường chéo, đối xứng và xác định dương.

    • Định lý 3: Chứng minh rằng với bước lưới $h = \max_k |x_{k+1} - x_k|$ đủ bé, hệ sai phân có nghiệm duy nhất và hội tụ về nghiệm đúng của bài toán xuất phát. Tác giả đánh giá chi tiết tính xấp xỉ $|(A\Phi)_h - f_h| \le Mh^{1/2} + Nh^{3/2} + Kh^{1/2} \le C h^{1/2}$ và chứng minh tính ổn định $|\Phi| \le \frac{1}{\rho_0} |F|$.
    • Định lý 4: Nếu bổ sung giả thiết hàm $\Phi(x)$ không tăng, cấp hội tụ của nghiệm sai phân về nghiệm đúng đạt cấp 1 đối với chuẩn năng lượng mở rộng.
    • Chương này cũng đề xuất và kiểm chứng lược đồ sai phân cải tiến thứ hai (công thức 1.20) cho các hàm hệ số biến thiên.

  • Bài toán một chiều dạng tổng quát: $$-\frac{d}{dx}\left[p(x)\frac{d\Phi(x)}{dx}\right] + r(x)\frac{d\Phi(x)}{dx} + q(x)\Phi(x) = f(x)$$ Trên từng đoạn $[x_i, x_{i+1}]$, nghiệm được xấp xỉ bằng hệ đa thức $g_k(x)$ bậc $m$ hoặc hệ hàm cơ sở lượng giác $\omega_{k,i}(x)$.

    • Định lý 5: Với các điều kiện ghép nối tại nút $\Phi_{N_0}(0) = \Phi_{N_n}(1) = 0$ và $\Phi_{N_{i-1}}(x_i) = \Phi_{N_i}(x_i)$, độ lệch của nghiệm gần đúng so với nghiệm đúng được ước lượng qua công thức: $$|\Phi - \Phi^*|_{L_2[0,1]} \le \text{const} \cdot N^{-\alpha}$$ trong đó $N = \min N_i$ và $\alpha$ là độ trơn từng khúc của nghiệm.

Chương 4: Các bài toán biên có biên kỳ dị

Chương 4 kết hợp đạo hàm trung bình với phương pháp phần tử biên (BEM) để giải các bài toán biên elliptic trên miền đóng $\Omega \subset \mathbb{R}^n$ có biên $\Gamma$ chứa các điểm góc (biên kỳ dị).

  • Bài toán với phương trình Laplace: $\nabla^2 u(x) = 0$ trong $\Omega$, với điều kiện biên tổ hợp tuyến tính $\alpha(x)u(x) + \beta(x)\frac{\partial u(x)}{\partial n} = \gamma(x)$ trên $\Gamma$.

    • Tác giả sử dụng hàm trọng lượng là nghiệm cơ bản $u^(\xi, x) = -\frac{1}{2\pi}\ln r(\xi, x)$ (bài toán 2D) hoặc $u^(\xi, x) = \frac{1}{4\pi r(\xi, x)}$ (bài toán 3D).
    • Tích phân từng phần hai lần và xử lý kỳ dị yếu tại biên bằng kỹ thuật bao quanh điểm kỳ dị $\xi$ bởi hình cầu bán kính $\epsilon$ rồi cho $\epsilon \to 0$, thiết lập phương trình tích phân biên: $$C(\xi) u(\xi) + \int_\Gamma u(x) q^(\xi, x) d\Gamma(x) = \int_\Gamma q(x) u^(\xi, x) d\Gamma(x)$$ với $q(x) = \frac{\partial u(x)}{\partial n}$, $q^(\xi, x) = \frac{\partial u^(\xi, x)}{\partial n(x)}$.
    • Rời rạc hóa biên $\Gamma$ thành $N$ phần tử biên $\Gamma_j$, đưa phương trình về hệ phương trình ma trận đại số $HU = GQ$, kết hợp điều kiện biên $EU + FQ = D$ thành hệ thống nhất $AY = B$.
    • Định lý 6: Chứng minh nghiệm của hệ đại số $AY = B$ hội tụ về nghiệm đúng của bài toán biên theo chuẩn năng lượng.
    • Định lý 7: Khi các hàm $u, q$ thay đổi không nhiều trên mỗi phần tử biên, sự tồn tại duy nhất nghiệm của hệ đại số tương đương với sự tồn tại duy nhất nghiệm của bài toán biên ban đầu.
    • Định lý 8: Với điều kiện biên Dirichlet thuần túy ($\beta(x) = 0$), hệ phương trình đại số thu được là hệ đối xứng.
    • Tác giả khảo sát chi tiết ba trường hợp rời rạc hóa: $u, q$ là hằng số trên mỗi phần tử biên; $u, q$ có dạng tuyến tính; $u$ có dạng tuyến tính và $q$ là hằng số.

  • Bài toán với phương trình Poisson: $\nabla^2 u(x) = b(x)$ trên miền có biên kỳ dị. Phương trình tích phân biên xuất hiện thêm tích phân khối $\int_\Omega b(x) u^*(\xi, x) d\Omega$.

    • Hệ đại số rời rạc có dạng $B + HU = GQ$ kết hợp $EU + FQ = D$.
    • Định lý 9, Định lý 10, Định lý 11: Thiết lập các kết quả tương ứng về tính hội tụ theo chuẩn năng lượng, tính giải được duy nhất và tính đối xứng của hệ phương trình ma trận khi $\beta(x) = 0$.

Chương 5: Một số phương pháp phối hợp giải bài toán biên không dừng trên miền có biên kỳ dị

Chương 5 mở rộng bài toán sang miền thời gian $t > 0$, xét phương trình truyền dẫn/truyền nhiệt không dừng: $$\nabla^2 u - a \frac{\partial u}{\partial t} - b u = f(x, t), \quad x \in \Omega, \ t > 0$$ với điều kiện ban đầu $u(x, 0) = u_0(x)$ và điều kiện biên tổ hợp $\alpha(x)u(x, t) + \beta(x)\frac{\partial u(x, t)}{\partial n} = 0$ trên biên $\Gamma$ có điểm góc kỳ dị.

  • Phương pháp biến đổi Laplace kết hợp phần tử biên:

    • Áp dụng biến đổi Laplace theo thời gian $t$, đưa phương trình về dạng elliptic phụ thuộc tham số $\lambda$: $$\nabla^2 \bar{u}(x, \lambda) - (a\lambda + b)\bar{u}(x, \lambda) = \bar{f}(x, \lambda) - a u_0(x)$$
    • Nghiệm cơ bản tương ứng trong không gian hai chiều chứa hàm Bessel biến dạng loại 2 bậc 0: $u^*(\xi, x, \lambda) = \frac{1}{2\pi} K_0\left(\sqrt{a\lambda + b} , r(\xi, x)\right)$.
    • Thiết lập phương trình tích phân biên cho $N$ giá trị tham số $\lambda_n$ chọn trước, giải bằng phương pháp phần tử biên để tìm $\bar{u}(\xi, \lambda_n)$, sau đó thực hiện phép biến đổi Laplace ngược xấp xỉ dạng: $$u(\xi, t) \approx u(\xi, \infty) + \sum_{n=1}^N a_n(\xi) \exp[-b_n(\xi) t]$$ trong đó các hệ số $a_n(\xi), b_n(\xi)$ được xác định từ hệ đại số xây dựng theo hướng dẫn của tài liệu [9].

  • Phương pháp đường thẳng (Method of Lines):

    • Xấp xỉ đạo hàm thời gian bằng sai phân hữu hạn lùi $\frac{\partial u}{\partial t} \approx \frac{u(x, t+\Delta t) - u(x, t)}{\Delta t}$, chuyển phương trình về bài toán Helmholtz tại mỗi bước thời gian: $$\nabla^2 u(x, t+\Delta t) - \left(\frac{a}{\Delta t} + b\right) u(x, t+\Delta t) = f(x, t+\Delta t) - \frac{a}{\Delta t} u(x, t)$$
    • Nghiệm cơ bản có dạng liên quan đến hàm Bessel $K_{1/2}$. Phương trình tích phân biên được giải tuần tự tại mỗi bước thời gian $t_k = t_0 + k\Delta t$, sử dụng nghiệm của bước trước làm điều kiện nguồn cho bước tiếp theo.

  • Phương pháp lặp tích phân: Trình bày thuật toán lặp giải phương trình tích phân biên không dừng qua các phép xấp xỉ liên tiếp.

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

Luận án đã đạt được các kết quả học thuật và đóng góp mới sau:

  1. Về mặt lý thuyết toán học:

    • Hệ thống hóa và mở rộng việc áp dụng khái niệm đạo hàm trung bình tích phân cho các hàm không trơn, hàm có điểm góc và hàm gián đoạn loại 1.
    • Thiết lập các định lý toán học xác lập mối quan hệ giữa đạo hàm trung bình tích phân và các phép biến đổi toán tử tích phân (biến đổi Laplace một chiều và hai chiều).

  2. Về mặt phương pháp tính và giải thuật số:

    • Đề xuất phương pháp phối hợp đạo hàm trung bình với biến đổi Laplace để giải chính xác các bài toán biên lan truyền ô nhiễm môi trường suy biến có nghiệm không trơn và nghiệm gián đoạn, chứng minh nghiệm trùng khớp hoàn toàn với nghiệm giải tích.
    • Xây dựng thành công các họ lược đồ sai phân mới cho bài toán biên có hệ số không trơn dựa trên đạo hàm trung bình. Chứng minh chặt chẽ tính duy nhất, tính ổn định và sự hội tụ của nghiệm sai phân với cấp hội tụ cụ thể theo bước lưới $h$ và bậc đa thức xấp xỉ.
    • Phát triển hoàn chỉnh phương pháp phối hợp đạo hàm trung bình và phương pháp phần tử biên (BEM) để giải lớp bài toán biên hỗn hợp cho phương trình Laplace và phương trình Poisson trên các miền có biên kỳ dị (chứa điểm góc), chứng minh sự hội tụ theo chuẩn năng lượng và tính đối xứng của ma trận hệ số.
    • Đề xuất quy trình tính toán phối hợp giải bài toán biên không dừng trên miền biên kỳ dị thông qua việc kết hợp phần tử biên với biến đổi Laplace ngược và phương pháp đường thẳng (sai phân thời gian).

Bảng tổng hợp các lớp bài toán và phương pháp giải trong luận án

Lớp bài toán Đặc trưng toán học Phương pháp phối hợp đề xuất Kết quả giải thuật chính
Bài toán khuếch tán ô nhiễm (Chương 2) Nghiệm gián đoạn loại 1 hoặc có điểm góc; chứa hàm $\delta(x)$ và $u^*(x)$. Đạo hàm trung bình + Biến đổi Laplace hai chiều Nghiệm giải tích chính xác; trùng khớp với cấu trúc nghiệm Mactruc (Định lý 2, Mệnh đề 1 & 2).
Bài toán hệ số gián đoạn (Chương 3) Hệ số $p(x), q(x), r(x)$ gián đoạn loại 1 trên $[0, 1]$. Đạo hàm trung bình + Lược đồ sai phân chóp/đa thức Ma trận 3 đường chéo đối xứng; nghiệm hội tụ bậc $h^{1/2}$ và bậc 1 (Định lý 3, 4, 5).
Bài toán biên Laplace / Poisson (Chương 4) Miền $\Omega$ có biên $\Gamma$ kỳ dị (điểm góc), điều kiện biên hỗn hợp. Đạo hàm trung bình + Phương pháp phần tử biên (BEM) Đưa về hệ đại số $AY=B$; hội tụ theo chuẩn năng lượng; đối xứng khi $\beta(x)=0$ (Định lý 6-11).
Bài toán không dừng (Chương 5) Phương trình parabolic trên miền biên kỳ dị, $t > 0$. BEM + Biến đổi Laplace / BEM + Phương pháp đường thẳng Giải chuỗi bài toán elliptic với hàm cơ bản Bessel $K_0, K_{1/2}$; nội suy ngược thời gian.

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

Căn cứ trên các nội dung và phạm vi nghiên cứu được trình bày trong văn bản luận án, các hạn chế và hướng nghiên cứu mở rộng bao gồm:

  • Phạm vi mô hình: Các bài toán khảo sát chủ yếu tập trung vào các phương trình vi phân tuyến tính (phương trình khuếch tán tuyến tính, phương trình Laplace, Poisson, truyền dẫn tuyến tính). Luận án chưa xét đến các phương trình phi tuyến mạnh hoặc hệ số phụ thuộc phi tuyến vào nghiệm $\Phi(x)$ hoặc $u(x)$.
  • Cấu trúc hình học của biên: Việc khảo sát số cho phần tử biên ở Chương 4 và Chương 5 chủ yếu thực hiện cho các biên kỳ dị dạng đường gấp khúc khép kín (2D) hoặc mặt đa diện (3D). Các dạng kỳ dị phức tạp hơn như kỳ dị đỉnh nón vô hạn hoặc biên fractal chưa nằm trong phạm vi giải thuật trực tiếp.
  • Hướng phát triển tiếp theo:
    • Mở rộng ứng dụng đạo hàm trung bình tích phân cho các hệ phương trình đạo hàm riêng phi tuyến nhiều chiều.
    • Nghiên cứu sâu hơn thuật toán lặp tích phân trong Chương 5 để nâng cao độ chính xác và tốc độ hội tụ khi thực hiện phép biến đổi Laplace ngược số.
    • Tối ưu hóa việc chia lưới thích nghi tự động tại lân cận các điểm góc kỳ dị và điểm gián đoạn của nguồn thải.

Giá trị tham khảo

Luận án mang lại giá trị tham khảo chuyên môn phục vụ các nhóm đối tượng:

  • Nghiên cứu sinh và học viên cao học chuyên ngành Toán học tính toán, Toán ứng dụng: Luận án cung cấp một tài liệu tham khảo về kỹ thuật áp dụng giải tích phi cổ điển (đạo hàm trung bình) vào việc thiết lập và chứng minh các định lý hội tụ, tính ổn định của các phương pháp sai phân và phương pháp phần tử biên.
  • Giảng viên và nhà nghiên cứu giải tích số: Các kỹ thuật xử lý tích phân chứa kỳ dị yếu tại biên điểm góc (Chương 4) và phương pháp giải phương trình không dừng bằng biến đổi tích phân kết hợp phần tử biên (Chương 5) là tài liệu mẫu mực về phương pháp luận toán tính toán.
  • Kỹ sư cơ học, vật lý kỹ thuật và mô phỏng môi trường: Các công thức nghiệm giải tích giải quyết bài toán khuếch tán từ nguồn ô nhiễm điểm (Chương 2) và các lược đồ sai phân cho môi trường có hệ số biến tính không liên tục (Chương 3) có thể trực tiếp ứng dụng vào các phần mềm tính toán lan truyền chất ô nhiễm và truyền nhiệt trong vật liệu dị tính.

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

1. Khái niệm đạo hàm trung bình tích phân khác biệt như thế nào so với đạo hàm Gâteaux thông thường?
Đạo hàm trung bình tích phân được định nghĩa thông qua giới hạn của tích phân trung bình biến phân $\lim_{\delta \to 0} \frac{1}{2\delta} \int_{-\delta}^{\delta} \frac{F(x_0+th) - F(x_0)}{t} dt$. Nhờ phép lấy tích phân làm trơn, đạo hàm trung bình tích phân có thể tồn tại tại những điểm mà hàm số có điểm góc hoặc thậm chí gián đoạn (nơi đạo hàm Gâteaux hoàn toàn không tồn tại). Ví dụ, hàm $F(x) = |x|$ không có đạo hàm Gâteaux tại $x=0$, nhưng có đạo hàm trung bình bằng 0 tại điểm này.

2. Làm thế nào luận án chứng minh được sự hội tụ của lược đồ sai phân cho bài toán có hệ số không trơn trong Chương 3?
Trong Chương 3, tác giả chứng minh sự hội tụ thông qua hai bước độc lập theo nguyên lý giải tích số: (i) Chứng minh tính xấp xỉ của lược đồ bằng cách đánh giá sai số $|(A\Phi)_h - f_h| \le Mh^{1/2} + Nh^{3/2} + Kh^{1/2} \le Ch^{1/2}$; và (ii) Chứng minh tính ổn định của toán tử sai phân dựa trên tính xác định dương của ma trận $A$, từ đó suy ra nghiệm sai phân hội tụ duy nhất về nghiệm đúng của bài toán xuất phát theo chuẩn năng lượng (Định lý 3).

3. Tại sao ma trận hệ phương trình đại số của phương pháp phần tử biên trong Chương 4 lại có tính chất đối xứng khi $\beta(x) = 0$?
Theo Định lý 8 của luận án, khi $\beta(x) = 0$, bài toán biên chuyển thành bài toán Dirichlet thuần túy. Lúc này, việc rời rạc hóa phương trình tích phân biên dẫn đến ma trận ảnh hưởng $G$ xây dựng từ nghiệm cơ bản $u^*(\xi, x) = -\frac{1}{2\pi}\ln r(\xi, x)$ mang tính đối xứng hình học tự nhiên qua khoảng cách $r(\xi, x) = r(x, \xi)$, làm cho toàn bộ hệ đại số trở thành hệ phương trình đối xứng.

4. Luận án giải quyết tính kỳ dị yếu của tích phân biên tại các điểm nút $\xi = x$ bằng cách nào?
Khi điểm quan sát $\xi$ nằm trên biên $\Gamma$ và tiến tới điểm lấy tích phân $x$, tích phân chứa nhân $\frac{\partial u^}{\partial n}$ và $u^$ xuất hiện kỳ dị yếu. Tác giả giải quyết bằng phương pháp hình học: cô lập điểm kỳ dị $\xi$ bằng một hình cầu nhỏ tâm $\xi$ bán kính $\epsilon$, chia biên thành phần trơn $\Gamma - \Gamma_\epsilon$ và phần biên cầu $\Gamma_\epsilon$, sau đó lấy giới hạn giải tích khi $\epsilon \to 0$ để triệt tiêu phần kỳ dị và thu được hệ số góc $C(\xi)$ cùng tích phân theo nghĩa giá trị chính Cauchy.

5. Trong Chương 5, nghiệm cơ bản nào được sử dụng khi áp dụng phương pháp phần tử biên cho phương trình sau biến đổi Laplace?
Sau khi áp dụng biến đổi Laplace theo thời gian, phương trình vi phân không dừng chuyển về phương trình dạng Helmholtz biến dạng $\nabla^2 \bar{u} - (a\lambda + b)\bar{u} = \bar{F}$. Trong không gian hai chiều, nghiệm cơ bản tương ứng của phương trình này là hàm Bessel biến dạng loại hai bậc không: $u^*(\xi, x, \lambda) = \frac{1}{2\pi} K_0\left(\sqrt{a\lambda + b} , r(\xi, x)\right)$.

Kết luận

Luận án "Đạo hàm trung bình và các phương pháp phối hợp để giải một số bài toán biên" của nghiên cứu sinh dưới sự hướng dẫn của thầy Phan Văn Hạp (1993) là một công trình khoa học chuyên sâu thuộc chuyên ngành Toán học tính toán. Bằng việc xây dựng và khai thác công cụ đạo hàm trung bình tích phân kết hợp với các phương pháp giải tích và giải số (biến đổi Laplace, sai phân hữu hạn, phần tử biên, phương pháp đường thẳng), luận án đã giải quyết thành công các bài toán biên phức tạp có nghiệm gián đoạn, hệ số không trơn và miền tính toán có biên kỳ dị. Các định lý về sự tồn tại, tính duy nhất, tính ổn định và sự hội tụ của các nghiệm giải tích cũng như nghiệm xấp xỉ số được chứng minh tường minh, tạo nên hệ thống phương pháp luận hoàn chỉnh và có giá trị khoa học cao trong lĩnh vực giải tích số và toán ứng dụng.