Tổng quan luận án

Luận án Phó tiến sĩ Toán Lý chuyên ngành Toán học tính toán (Mã số: 1-01-07) của tác giả tại Khoa Toán - Cơ - Tin học, Trường Đại học Tổng hợp Hà Nội (năm 1993), dưới sự hướng dẫn của Thầy hướng dẫn Phan Văn Háp, mang tiêu đề: "Đạ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".

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

Trong thực tế mô hình hóa các quá trình vật lý, cơ học và môi trường (như hiện tượng khuếch tán chất ô nhiễm trong không khí từ nguồn điểm, hiện tượng truyền nhiệt trong môi trường không đồng chất), các phương trình vi phân thường chứa các hệ số gián đoạn, nghiệm có điểm góc hoặc bước nhảy (nghiệm không trơn), hoặc miền tính toán có biên không trơn chứa điểm góc (biên kỳ dị). Tại các điểm đặc biệt này, các khái niệm đạo hàm cổ điển cũng như đạo hàm Gâteaux không tồn tại. Việc áp dụng các phương pháp giải tích truyền thống hoặc phương pháp số thông thường gặp nhiều trở ngại về mặt cơ sở giải tích và tính hội tụ. Luận án đặt ra yêu cầu nghiên cứu việc sử dụng 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 toán tử (biến đổi Laplaxơ, phương pháp sai phân biến phân, 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 và nghiệm số cho các lớp bài toán biên đặc thù nêu trên.

Mục tiêu và câu hỏi nghiên cứu

Mục tiêu tổng quát của luận án là ứng dụng khái niệm đạo hàm trung bình tích phân và phát triển các phương pháp phối hợp để giải quyết ba lớp bài toán biên: bài toán có nghiệm không trơn, bài toán có hệ số không trơn và bài toán trên miền có biên kỳ dị (cho cả trường hợp dừng và không dừng).

Các vấn đề nghiên cứu cốt lõi được định hình gồm:

  1. Xác định mối liên hệ giải tích giữa đạo hàm trung bình tích phân với phép biến đổi Laplaxơ đối với các hàm có điểm góc hoặc gián đoạn loại 1; kiểm chứng tính đúng đắn của phương pháp phối hợp này trên mô hình khuếch tán ô nhiễm môi trường.
  2. Thiết lập lược đồ sai phân dựa trên đạo hàm trung bình cho phương trình vi phân có hệ số không trơn một chiều, chứng minh tính ổn định và xác định cấp hội tụ của nghiệm xấp xỉ.
  3. Ứng dụng phương pháp cân bằng sai số và đạo hàm trung bình để xây dựng phương pháp phần tử biên (Boundary Element Method - BEM) cho phương trình Laplaxơ và Poatxông với điều kiện biên hỗn hợp trên miền có biên kỳ dị.
  4. Mở rộng các phương pháp phối hợp (biến đổi Laplaxơ, phương pháp đường thẳng) để giải bài toán biên không dừng (dạng parabolic) trên miền có biên kỳ dị.

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

  • Đối tượng nghiên cứu: Phương trình vi phân đạo hàm riêng và phương trình vi phân thường cấp hai (phương trình khuếch tán - dịch chuyển, phương trình elliptic dạng Laplaxơ, Poatxông và phương trình không dừng dạng parabolic); các hệ số gián đoạn; nghiệm không trơn; miền có biên kỳ dị (chứa điểm góc).
  • Phạm vi không gian và thời gian: Không gian một chiều ($x \in \mathbb{R}$ hoặc $[0,1]$), không gian hai chiều và ba chiều ($\Omega \subset \mathbb{R}^2, \mathbb{R}^3$); miền thời gian $t > 0$.
  • Điều kiện biên: Điều kiện suy giảm ở vô hạn, điều kiện biên Dirichlet, Neumann và điều kiện biên tổ hợp/hỗn hợp (Robin).

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

Trong phần dẫn nhập của từng chương, tác giả đã tổng thuật các hướng nghiên cứu liên quan trực tiếp đến đề tài:

Các hướng nghiên cứu trước được điểm lại

  1. Hướng nghiên cứu bài toán biên có nghiệm không trơn trong ô nhiễm môi trường: Tác giả Mactruc [3] đã đưa ra các mô hình toán học mô tả quá trình khuếch tán và dịch chuyển của thực thể trong không khí với vận tốc gió khác không và bằng không, sinh ra bởi nguồn ô nhiễm có kích thước nhỏ so với phạm vi ảnh hưởng. Mô hình này dẫn đến các phương trình vi phân có vế phải chứa hàm Dirac $\delta(x-x_0)$ và nghiệm giải tích bị gián đoạn loại 1 hoặc có điểm góc tại vị trí nguồn.
  2. Hướng nghiên cứu bài toán biên có hệ số không trơn: Các tác giả Ladujenxkaia [5], Ianhenkô [8], Xamarxki [7] đã nghiên cứu giải các bài toán biên có hệ số không trơn bằng phương pháp xây dựng các lược đồ sai phân dựa trên các nguyên lý biến phân.
  3. Hướng nghiên cứu bài toán biên trên miền có biên kỳ dị: Các tác giả Ianhenkô [3], Xamarxki [7], Gôdunốp [2] đã khảo sát bài toán biên cho phương trình Laplaxơ và Poatxông với các điều kiện biên Dirichlet và Neumann bằng các phương pháp khai triển tổng, phương pháp phân rã.
  4. Cơ sở phương pháp luận bổ trợ: Khái niệm biến phân trung bình tích phân thứ nhất được kế thừa từ tài liệu [13]; phương pháp nghiệm xấp xỉ và phương pháp cân bằng sai số dựa trên tài liệu [1]; kỹ thuật lựa chọn tham số trong biến đổi Laplaxơ ngược tham khảo tài liệu [9].

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

  • Các phương pháp sai phân biến phân cổ điển đòi hỏi kỹ thuật giải tích phức tạp khi khử các bước nhảy của hệ số. Luận án sử dụng trực tiếp khái niệm đạo hàm trung bình kết hợp hệ hàm cơ sở để đơn giản hóa việc thiết lập lược đồ sai phân ba điểm.
  • Các công trình trước về miền biên kỳ dị chủ yếu tập trung vào điều kiện biên thuần nhất Dirichlet hoặc Neumann. Luận án giải quyết bài toán với điều kiện biên hỗn hợp dạng tổ hợp $\alpha(x)u(x) + \beta(x)\frac{\partial u}{\partial n} = \Phi(x)$ trên biên có điểm góc bằng phương pháp phần tử biên.
  • Luận án xây dựng quy trình phối hợp biến đổi Laplaxơ hoặc phương pháp đường thẳng với phương pháp phần tử biên để xử lý bài toán phương trình vi phân không dừng trên miền có biên kỳ dị.

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

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

  1. Đạo hàm trung bình tích phân:

    • Xét ánh xạ $F: X \to Y$ ($X, Y$ là các không gian Banach), $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ì $\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 $x_0$.
    • Nếu $\delta F(x_0, h) = Ah$ với $A$ là toán tử tuyến tính giới nội, thì $F$ được gọi là khả vi trung bình, có đạo hàm trung bình theo tích phân tại $x_0$, ký hiệu $F'(x_0, h) = Ah$.
    • Tính chất cơ bản: Tuyến tính $(\alpha F + \beta G)'(x_0) = \alpha F'(x_0) + \beta G'(x_0)$. Một hàm có thể không có đạo hàm Gâteaux tại $x_0$ nhưng vẫn có đạo hàm trung bình (ví dụ $F(x) = |x|$ tại $x_0=0$ có $F'(0)=0$). Nếu hàm có đạo hàm hữu hạn hai phía tại $x_0$ thì: $$F'(x_0) = \frac{F'(x_0+0) + F'(x_0-0)}{2}$$

  2. Phương pháp cân bằng sai số (Error Balancing Method):

    • Để tìm nghiệm xấp xỉ $u^* = \sum_{i=1}^n \alpha_i \varphi_i + \alpha_0$ của phương trình toán tử $L(u) = b$ trong miền $\Omega$ với điều kiện biên $S(u)=s$ trên $\Gamma_1$, $G(u)=g$ trên $\Gamma_2$, sai số $R = L(u^*) - b$ được cân bằng thông qua tích phân trọng lượng: $$\int_\Omega R \cdot W d\Omega = 0$$ trong đó $W$ là hàm trọng lượng.

  3. Phép biến đổi toán tử Laplaxơ:

    • Biến đổi Laplaxơ một chiều: $\mathcal{L}[f(t), s] = \int_0^{+\infty} f(t) e^{-st} dt$.
    • Biến đổi Laplaxơ hai chiều: $\mathcal{L}2[f(t), s] = \int{-\infty}^{+\infty} f(t) e^{-st} dt = \mathcal{L}[f(t), s] + \mathcal{L}[f(-t), -s]$.
    • Mối liên hệ với đạo hàm trung bình: Nếu $f(t)$ giới nội, có đạo hàm trung bình và có các điểm gián đoạn $t_1, t_2, \dots$ với giới hạn hai phía hữu hạn: $$\mathcal{L}[f'(t), s] = sF(s) - f(0+0) - \sum_i \exp(-t_i s) [f(t_i+0) - f(t_i-0)]$$

Phương pháp nghiên cứu thực tế

Luận án thuộc hướng nghiên cứu toán học tính toán, kết hợp giải tích lý thuyết và giải thuật số:

  • Phương pháp phối hợp giải tích - toán tử: Áp dụng biến đổi Laplaxơ hai chiều trực tiếp lên phương trình vi phân có nghiệm không trơn chứa đạo hàm trung bình để tìm cấu trúc nghiệm giải tích.
  • Phương pháp sai phân hữu hạn (FDM) biến phân: Sử dụng phép nhân vô hướng với hệ hàm cơ sở nón $\omega_k(x)$ và tích phân từng phần hai lần để chuyển phương trình vi phân có hệ số không trơn về hệ đại số tuyến tính ma trận 3 đường chéo.
  • Phương pháp phần tử biên (BEM): Sử dụng nghiệm cơ bản của toán tử Laplaxơ, kỹ thuật cách ly điểm kỳ dị bằng lùi biên (mặt cầu bán kính $\varepsilon \to 0$) để thiết lập phương trình tích phân biên; rời rạc hóa biên thành các phần tử hằng số và tuyến tính.
  • Phương pháp đường thẳng (Method of Lines): Rời rạc hóa biến thời gian bằng sai phân để đưa phương trình parabolic về chuỗi phương trình dạng elliptic tại mỗi bước thời gian, sau đó áp dụng BEM để giải.

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

LUẬN ÁN: ĐẠO HÀM TRUNG BÌNH VÀ CÁC PHƯƠNG PHÁP PHỐI HỢP

Chương 1: Mở đầu

Chương này trình bày cơ sở lý thuyết toán học phục vụ toàn bộ luận án:

  • Định nghĩa giải tích của đạo hàm trung bình tích phân thứ nhất $\delta F(x_0, h)$ và toán tử khả vi trung bình $F'(x_0, h) = Ah$.
  • Chứng minh các tính chất giải tích: tính tuyến tính, tính chất qua điểm góc, tính chất qua điểm gián đoạn và công thức liên hệ đạo hàm hai phía.
  • Trình bày khái niệm nghiệm xấp xỉ và phương pháp cân bằng sai số theo hệ hàm đầy đủ ${\varphi_i}$.
  • Trình bày định nghĩa, điều kiện tồn tại, công thức biến đổi ngược của phép biến đổi Laplaxơ một chiều và hai chiều; thiết lập định lý liên hệ giữa biến đổi Laplaxơ với đạo hàm trung bình tích phân khi hàm gốc có bước nhảy.

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

Chương 2 khảo sát các bài toán biên xuất phát từ mô hình toán học ô nhiễm môi trường do Mactruc [3] đề xuất:

  1. Bài toán biên dạng phức tạp có nghiệm gián đoạn (§1):

    • Phương trình khuếch tán và chuyển dịch của thực thể có vận tốc gió: $$\mu \frac{d^2\Phi(x)}{dx^2} - \alpha \frac{d\Phi(x)}{dx} + \beta [\Phi(x) - \lambda u^(x-x_0)] = \eta + 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) = \lambda$ ($\lambda > 0$, $\mu > 0$, $u^$ là hàm đơn vị Heaviside, $\delta$ là hàm Dirac).
    • Nghiệm giải tích của bài toán bị gián đoạn loại 1 tại $x = x_0$, do đó không có đạo hàm Gâteaux nhưng có đạo hàm trung bình tích phân.
    • Định lý 2: Chứng minh rằng nghiệm tìm được bằng phương pháp phối hợp đạo hàm trung bình tích phân và biến đổi Laplaxơ hai chiều hoàn toàn trùng khớp với nghiệm giải tích: $$\Phi(x) = \begin{cases} \frac{Q}{\sqrt{4B\mu + \alpha^2}} \exp\left[ -\left(\frac{\sqrt{4B\mu+\alpha^2}}{2\mu} - \frac{\alpha}{2\mu}\right)(x-x_0) \right] + \lambda, & x > x_0 \ \frac{Q}{\sqrt{4B\mu + \alpha^2}} \exp\left[ \left(\frac{\sqrt{4B\mu+\alpha^2}}{2\mu} + \frac{\alpha}{2\mu}\right)(x_0-x) \right], & x < x_0 \end{cases}$$

  2. Bài toán biên dạng phức tạp có nghiệm không trơn (§2):

    • Xét trường hợp $\lambda = 0$: $\mu \frac{d^2\Phi}{dx^2} - \alpha \frac{d\Phi}{dx} + \beta \Phi = \mu Q \delta(x-x_0)$, $\lim_{x \to \pm\infty} \Phi(x) = 0$.
    • Mệnh đề 1: Nghiệm tìm theo phương pháp phối hợp đạo hàm trung bình và biến đổi Laplaxơ trùng với nghiệm giải tích, đồ thị nghiệm liên tục nhưng có điểm góc tại $x = x_0$ khi $\alpha > 0$.

  3. Bài toán biên dạng đơn giản có nghiệm không trơn (§3):

    • Xét trường hợp vận tốc không khí bằng không ($\alpha = 0$): $-\mu \frac{d^2\Phi}{dx^2} + \beta \Phi(x) = Q \delta(x-x_0)$, $\lim_{x \to \pm\infty} \Phi(x) = 0$.
    • Mệnh đề 2: Nghiệm giải tích đối xứng qua $x_0$ và trùng với nghiệm tìm bằng phương pháp phối hợp: $\Phi(x) = \frac{Q}{2\sqrt{\beta\mu}} \exp\left(-\sqrt{\frac{\beta}{\mu}}|x-x_0|\right)$.

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 phương pháp số giải bài toán biên một chiều có hệ số gián đoạn:

  1. Bài toán một chiều dạng đơn giản (§1):

    • Xét phương trình: $-\frac{d}{dx}\left[p(x)\frac{d\Phi(x)}{dx}\right] + q(x)\Phi(x) = f(x)$, $0 < x < 1$, điều kiện biên $\Phi(0) = \Phi(1) = 0$.
    • Các hệ số $p(x), q(x)$ và nghiệm $\Phi(x)$ tồn tại đạo hàm trung bình và giới nội.
    • Luận án thiết lập lưới ${x_k}, {x_{k+1/2}}$ và hệ hàm cơ sở nón $\omega_k(x)$ xác định bởi: $$\omega_k(x) = \begin{cases} \frac{x - x_{k-1}}{\Delta_{k-1/2}}, & x \in [x_{k-1}, x_k] \ \frac{x_{k+1} - x}{\Delta_{k+1/2}}, & x \in [x_k, x_{k+1}] \ 0, & x \notin [x_{k-1}, x_{k+1}] \end{cases}$$
    • Thực hiện tích phân từng phần hai lần phương trình nhân với $\omega_k(x)$ để loại bỏ đạo hàm của hệ số không trơn $p(x)$, dẫn đến hệ phương trình đại số tuyến tính: $$A\Phi = F$$ trong đó ma trận $A$ là ma trận 3 đường chéo chính, đối xứng và xác định dương.
    • Định lý 3: Với bước lưới $h = \max_k |x_{k+1} - x_k|$ đủ bé, hệ phương trình sai phân có nghiệm duy nhất và hội tụ về nghiệm đúng của bài toán liên tục với đánh giá sai số: $$|\Phi_h - \Phi^*|_F \le c h^{1/2} \quad (c > 0)$$
    • Định lý 4: Nếu giả thiết thêm $\Phi(x)$ là hàm không tăng thì cấp hội tụ của nghiệm sai phân đạt cấp 1 theo chuẩn năng lượng mở rộng: $$|\Phi|F^2 = \sum{i=1}^{n-1} \Delta_i \left( \sum_{j=1}^i \Delta_j f_j \right)^2$$
    • Luận án cũng đưa ra lược đồ sai phân cải tiến tại công thức (1.20) với các hệ số chứa đạo hàm trung bình tại các nửa bước lưới.

  2. Bài toán một chiều dạng tổng quát (§2):

    • Xét phương trình: $-\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)$ với $p, r, q, \Phi$ là các hàm gián đoạn loại 1.
    • Sử dụng đa thức xấp xỉ bậc $m$ hoặc hệ hàm lượng giác $\Phi_{Ni}(x) = \sum_{k=0}^{2N_i} C_k u_k(x)$ trên từng đoạn lưới $[x_i, x_{i+1}]$.
    • Định lý 5: Thiết lập đánh giá sai số của nghiệm xấp xỉ $\Phi^(x)$ so với nghiệm đúng: $$|\Phi - \Phi^|_{L_2[0,1]} \le \frac{\text{const}}{N^\alpha}$$ trong đó $N = \min N_i$, và $\alpha$ là độ trơn của nghiệm.

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

Chương 4 nghiên cứu giải pháp số cho các bài toán biên elliptic trên miền đóng giới nội $\Omega \subset \mathbb{R}^n$ có biên $\Gamma$ chứa điểm góc (biên kỳ dị), chịu điều kiện biên tổ hợp: $$\alpha(x)u(x) + \beta(x)\frac{\partial u(x)}{\partial n} = \Phi(x), \quad x \in \Gamma$$

  1. Phương trình Laplaxơ $\nabla^2 u(x) = 0$ (§1):

    • Áp dụng phương pháp cân bằng sai số với hàm trọng lượng $W(\xi, x) = u^(\xi, x)$ là nghiệm cơ bản của phương trình $\nabla^2 u^ = -\delta(\xi, x)$: $$u^*(\xi, x) = \begin{cases} \frac{1}{2\pi} \ln \frac{1}{r(\xi, x)}, & (n = 2) \ \frac{1}{4\pi r(\xi, x)}, & (n = 3) \end{cases}$$
    • Tích phân từng phần hai lần đưa về phương trình tích phân: $$u(\xi) + \int_\Gamma u(x) \frac{\partial u^(\xi, x)}{\partial n(x)} d\Gamma(x) = \int_\Gamma \frac{\partial u(x)}{\partial n(x)} u^(\xi, x) d\Gamma(x)$$
    • Khi cho điểm quan sát $\xi \to \Gamma$, tích phân chứa kỳ dị yếu được xử lý bằng cách lùi biên qua hình cầu tâm $\xi$ bán kính $\varepsilon \to 0$, thu được 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 = \frac{\partial u}{\partial n}$, $q^* = \frac{\partial u^*}{\partial n}$.
    • Chia biên $\Gamma$ thành $N$ phần tử biên $\Gamma_1, \dots, \Gamma_N$. Luận án xây dựng ma trận cho 3 mô hình rời rạc hóa:
      • $u, q$ là hằng số trên mỗi phần tử $\Gamma_j$ (chọn 1 nút tại trung điểm).
      • $u, q$ biến thiên tuyến tính trên mỗi phần tử $\Gamma_j$ (chọn 2 nút tại hai đầu mút).
      • $u$ biến thiên tuyến tính, $q$ là hằng số trên mỗi phần tử.
    • Hợp nhất với điều kiện biên dẫn đến hệ phương trình đại số ma trận: $$HU = GQ \quad \text{kết hợp} \quad EU + FQ = D \implies 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: Chứng minh sự tồn tại duy nhất nghiệm của hệ phương trình đạ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 liên tục khi $u, q$ biến thiên ít trên từng phần tử.
    • Định lý 8: Khi $\beta(x) = 0$ (điều kiện Dirichlet), ma trận $G$ là ma trận đối xứng.

  2. Phương trình Poatxông $\nabla^2 u(x) = b(x)$ (§2):

    • Xuất hiện thêm số hạng tích phân thể tích $\int_\Omega b(x) u^*(\xi, x) d\Omega(x)$ được tính bằng phương pháp số.
    • Đưa bài toán về hệ đại số dạng: $$B + HU = GQ \quad \text{và} \quad EU + FQ = D$$
    • Định lý 9, 10, 11: Khẳng định sự hội tụ của nghiệm, tính duy nhất nghiệm và tính đối xứng của hệ phương trình 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 các kỹ thuật trên cho phương trình parabolic không dừng: $$\nabla^2 u - a \frac{\partial u}{\partial t} - bu = f(x,t), \quad t > 0, x \in \Omega$$ 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.

  1. Phương pháp biến đổi Laplaxơ (§1):

    • Tác động biến đổi Laplaxơ theo biến thời gian: $\bar{u}(x, \lambda) = \int_0^{+\infty} u(x, t) e^{-\lambda t} dt$, chuyển phương trình về dạng elliptic tham số: $$\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 2 chiều biểu diễn qua hàm Bessel biến dạng loại 2 bậc 0 ($K_0$): $$u^*(\xi, x, \lambda) = \frac{1}{2\pi} K_0\left(\sqrt{a\lambda + b} \cdot r(\xi, x)\right)$$
    • Thiết lập phương trình tích phân biên và giải số cho $N$ giá trị tham số $\lambda$ chọn trước.
    • Nghiệm trong miền thời gian được xấp xỉ bằng chuỗi hàm mũ qua phép biến đổi ngược: $$u(\xi, t) = 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ừ các giá trị $\lambda$ đã chọn.

  2. Phương pháp đường thẳng (§2):

    • Rời rạc hóa đạo hàm thời gian bằng sai phân: $\frac{\partial u}{\partial t} \approx \frac{u(x, t+\Delta t) - u(x, t)}{\Delta t}$.
    • Phương trình tại mỗi bước thời gian trở thành: $$\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 tương ứng có dạng hàm Bessel $K_0\left(\sqrt{\frac{a}{\Delta t} + b} \cdot r\right)$ trong $\mathbb{R}^2$ hoặc $K_{1/2}$ trong $\mathbb{R}^3$.
    • Tại mỗi bước thời gian $t_k = t_{k-1} + \Delta t$, giải phương trình tích phân biên để tìm $u(x, t_k)$, lấy kết quả làm điều kiện ban đầu cho bước tiếp theo.

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

Chương Phương trình toán tử Đặc trưng kỳ dị / phi trơn Phương pháp phối hợp sử dụng Kết quả toán học chính
Chương 2 $\mu \Phi'' - \alpha \Phi' + \beta \Phi = \eta + Q\delta$ Nghiệm gián đoạn loại 1 hoặc có điểm góc tại $x_0$ Đạo hàm trung bình tích phân + Biến đổi Laplaxơ 2 chiều Định lý 2; Mệnh đề 1, 2: Nghiệm phối hợp trùng nghiệm giải tích
Chương 3 $-(p\Phi')' + q\Phi = f$ Hệ số $p(x), q(x)$ và nghiệm gián đoạn loại 1 Đạo hàm trung bình + Sai phân biến phân với hàm cơ sở nón Định lý 3, 4: Ma trận 3 đường chéo đối xứng; hội tụ cấp $h^{1/2}$ và cấp 1
Chương 4 $\nabla^2 u = 0$ và $\nabla^2 u = b(x)$ Miền $\Omega$ có biên $\Gamma$ chứa điểm góc; biên hỗn hợp Phương pháp cân bằng sai số + BEM (nghiệm cơ bản $\ln(1/r)$, $1/r$) Định lý 6, 7, 8, 9, 10, 11: Hội tụ theo chuẩn năng lượng; hệ $HU=GQ$
Chương 5 $\nabla^2 u - a u_t - bu = f(x,t)$ Bài toán không dừng (parabolic); biên $\Gamma$ có điểm góc Biến đổi Laplaxơ / Phương pháp đường thẳng + BEM (hàm Bessel $K_\nu$) Giải tích - số: Phương trình tích phân biên tham số và bước thời gian

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

Đóng góp về mặt lý luận

  1. Phát triển ứng dụng giải tích của đạo hàm trung bình tích phân: Thiết lập mối liên hệ toán học chặt chẽ giữa đạo hàm trung bình tích phân và phép biến đổi Laplaxơ đối với các hàm có bước nhảy gián đoạn hữu hạn. Chứng minh tính tương đương giữa nghiệm tìm bằng phương pháp toán tử này với nghiệm giải tích của phương trình truyền chất ô nhiễm từ nguồn điểm.
  2. Lý thuyết sai phân cho bài toán có hệ số không trơn: Xây dựng thành công kỹ thuật tích phân từng phần dựa trên đạo hàm trung bình để khử đạo hàm của hệ số gián đoạn. Chứng minh các định lý về tính xác định dương, tính đối xứng của ma trận hệ số, tính duy nhất nghiệm, tính ổn định và đánh giá sai số hội tụ ($\mathcal{O}(h^{1/2})$ tổng quát và cấp 1 khi nghiệm đơn điệu).
  3. Mở rộng phương pháp phần tử biên: Thiết lập cơ sở giải tích cho phương pháp phần tử biên áp dụng trên miền có biên chứa điểm góc với điều kiện biên tổ hợp (Robin) cho phương trình Laplaxơ và Poatxông. Chứng minh sự hội tụ của nghiệm xấp xỉ phần tử biên theo chuẩn năng lượng.

Đóng góp về mặt thực tiễn và giải thuật tính toán

  1. Quy trình số hóa phần tử biên linh hoạt: Xây dựng chi tiết công thức tính toán các hệ số ma trận $H$ và $G$ cho các phần tử biên với các phép gán nghiệm hằng số và tuyến tính, cho phép xử lý hiệu quả các miền hình học phức tạp có góc nhọn hoặc đường gấp khúc khép kín.
  2. Thuật toán giải phương trình không dừng trên miền biên kỳ dị: Đề xuất hai thuật toán số phối hợp hoàn chỉnh: (1) Thuật toán kết hợp biến đổi Laplaxơ - BEM đưa về việc giải hệ phương trình biên tại một số điểm tham số và nội suy hàm mũ ngược; (2) Thuật toán đường thẳng - BEM giải tuần tự các bài toán biên elliptic tựa dừng tại từng bước thời gian sử dụng hàm Bessel biến dạng.

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

Dựa trên phạm vi và các giả thiết được nêu trong văn bản luận án:

  • Hạn chế về mô hình: Khái niệm đạo hàm trung bình tích phân được luận án sử dụng cố định theo định nghĩa tích phân đối xứng xác định trên khoảng $[-\delta, \delta]$.
  • Hạn chế về không gian và hệ số: Các bài toán biên có hệ số không trơn trong Chương 3 chủ yếu được khảo sát chi tiết trên miền một chiều $[0, 1]$. Các bài toán không dừng trong Chương 5 xét với các hệ số hằng số $a, b$.
  • Hướng mở rộng: Mở rộng phương pháp sai phân dựa trên đạo hàm trung bình cho các bài toán biên hai và ba chiều có hệ số gián đoạn trên các mặt phân chia phức tạp; nghiên cứu mở rộng phương pháp cho các phương trình đạo hàm riêng phi tuyến.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên khảo dành cho:

  • 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, Giải tích số: Tham khảo kỹ thuật chứng minh tính ổn định và hội tụ của lược đồ sai phân cho bài toán có hệ số gián đoạn (Chương 3); cách thiết lập phương trình tích phân biên và xử lý kỳ dị tích phân bằng kỹ thuật lùi biên (Chương 4).
  • Giảng viên và nhà nghiên cứu Cơ học tính toán, Vật lý toán: Tham khảo phương pháp giải bài toán truyền nhiệt, khuếch tán không dừng trên miền có hình học phức tạp chứa điểm góc mà không cần làm trơn biên (Chương 5).
  • Kỹ sư và nhà mô hình hóa môi trường: Tham khảo cơ sở toán học và giải pháp giải tích cho bài toán lan truyền chất ô nhiễm không khí từ nguồn điểm gây ra nghiệm gián đoạn (Chương 2).

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

1. Tại sao hàm $F(x) = |x|$ không có đạo hàm Gâteaux tại $x=0$ nhưng lại có đạo hàm trung bình tích phân?

Theo định nghĩa đạo hàm Gâteaux, giới hạn tỉ sai phân khi $t \to 0$ từ hai phía trái và phải có giá trị khác nhau (bằng $-|h|$ và $+|h|$), nên không tồn tại giới hạn duy nhất. Tuy nhiên, theo định nghĩa đạo hàm trung bình tích phân: $$F'(0) = \lim_{\delta \to 0} \frac{1}{2\delta} \int_{-\delta}^{\delta} |h| \text{sign}(t) dt$$ Do hàm dưới dấu tích phân $\text{sign}(t)$ là hàm lẻ trên đoạn đối xứng $[-\delta, \delta]$, tích phân luôn triệt tiêu bằng 0 với mọi $\delta > 0$, dẫn đến $F'(0) = 0$.

2. Nghiệm của bài toán Mactruc trong Chương 2 có đặc tính gì tại vị trí nguồn ô nhiễm $x = x_0$?

Tại điểm $x = x_0$, hàm nghiệm $\Phi(x)$ có bước nhảy gián đoạn loại 1 với độ lệch $\lambda$ giữa hai phía $x > x_0$ và $x < x_0$ (khi $\lambda > 0$), hoặc đồ thị có điểm góc nhọn liên tục nhưng đạo hàm cấp một gián đoạn (khi $\lambda = 0, \alpha > 0$). Do đó, nghiệm không khả vi cổ điển nhưng thỏa mãn phương trình vi phân khi đạo hàm được hiểu theo nghĩa đạo hàm trung bình tích phân.

3. Ma trận $A$ thu được từ lược đồ sai phân trong Chương 3 có những tính chất đại số nào?

Ma trận $A$ thu được từ lược đồ sai phân ba điểm (khi tích phân từng phần với hệ hàm cơ sở $\omega_k(x)$) là ma trận vuông cấp $(n-1)$, có cấu trúc ba đường chéo chính, mang tính đối xứng và xác định dương. Tính xác định dương này bảo đảm hệ phương trình đại số luôn có nghiệm duy nhất.

4. Tích phân biên chứa kỳ dị yếu tại $\xi = x$ trong Chương 4 được xử lý bằng kỹ thuật nào?

Tác giả sử dụng kỹ thuật lùi biên: cô lập điểm kỳ dị $\xi \in \Gamma$ bằng một hình cầu tâm $\xi$ bán kính $\varepsilon$, ký hiệu phần biên của hình cầu nằm trong miền là $\Gamma_\varepsilon^$. Sau đó, tích phân được tách thành tích phân trên phần biên ngoài $\Gamma - \Gamma_\varepsilon$ và trên $\Gamma_\varepsilon^$. Khi cho $\varepsilon \to 0$, tích phân trên $\Gamma_\varepsilon^*$ sinh ra hệ số tự do $C(\xi)$ (với $C(\xi) = 1/2$ trên biên trơn và phụ thuộc vào góc trong tại điểm góc), đưa phương trình về dạng phương trình tích phân Fredholm.

5. Phương pháp đường thẳng trong Chương 5 biến đổi bài toán parabolic thành bài toán nào?

Phương pháp đường thẳng rời rạc hóa đạo hàm thời gian $\frac{\partial u}{\partial t}$ bằng sai phân tiến bước $\Delta t$. Phép biến đổi này khử biến thời gian và chuyển phương trình parabolic ban đầu thành một chuỗi phương trình vi phân elliptic dừng dạng Helmholtz/Poatxông đối với ẩn số $u(x, t+\Delta t)$ tại từng bước thời gian, với nghiệm cơ bản được biểu diễn qua hàm Bessel biến dạng loại hai $K_0$ hoặc $K_{1/2}$.

Kết luận

Luận án của tác giả đã hệ thống hóa và ứng dụng công cụ đạo hàm trung bình tích phân để giải quyết các bài toán biên có tính chất phi trơn và kỳ dị trong toán lý tính toán. Bằng việc kết hợp đạo hàm trung bình với biến đổi Laplaxơ, phương pháp sai phân biến phân và phương pháp phần tử biên, công trình đã thiết lập các nghiệm giải tích chính xác cho mô hình ô nhiễm môi trường và đề xuất các lược đồ số ổn định, hội tụ cho phương trình Laplaxơ, Poatxông cũng như phương trình không dừng trên miền có biên kỳ dị. Các kết quả của luận án đóng góp cơ sở giải tích và thuật toán số hoàn chỉnh cho lớp bài toán biên chứa hệ số hoặc biên không trơn.