Tổng quan luận án

Luận án Phó Tiến sĩ Toán Lý của tác giả Nguyễn Cam có tiêu đề chính thức: "Bài toán ngược và lời giải xấp xỉ ổn định ứng dụng vào phương trình tích phân", được hoàn thành năm 1995 tại Khoa Toán, Trường Đại học Sư phạm Thành phố Hồ Chí Minh, dưới sự hướng dẫn khoa học của Tiến sĩ Quốc gia Trần Văn Tấn (Trường Đại học Sư phạm Thành phố Hồ Chí Minh) và Giáo sư Tiến sĩ Đặng Đình Áng (Trường Đại học Tổng hợp Thành phố Hồ Chí Minh). Công trình thuộc chuyên ngành Giải tích hàm, mã hiệu 1.01, được phân tích và đánh giá chuyên môn với sự tham gia của Phó Giáo sư Dương Minh Đức cùng các đồng nghiệp tại hai cơ sở đào tạo nêu trên.

Về tính cấp thiết và khoảng trống nghiên cứu, luận án đặt trọng tâm vào lý thuyết bài toán ngược không chỉnh (ill-posed inverse problems). Theo tác giả trình bày, từ năm 1960 trở đi, các bài toán ngược không chỉnh bắt đầu được khảo sát sâu rộng trên thế giới thông qua các công trình nền tảng của A. N. Tikhonov, M. M. Lavrentiev và J. L. Lions. Lĩnh vực này có nhiều ứng dụng quan trọng trong y học, kỹ nghệ, địa vật lý và kỹ thuật. Tuy nhiên, theo nghĩa Hadamard, các bài toán ngược thường không thỏa mãn điều kiện tồn tại nghiệm, nghiệm không duy nhất hoặc nghiệm không phụ thuộc liên tục vào dữ liệu ban đầu. Sự mất ổn định này dẫn đến việc các dao động hoặc sai số nhỏ trong dữ liệu đo đạc thực tế có thể gây ra sai lệch vô hạn trong nghiệm giải tích. Do đó, yêu cầu nghiên cứu đặt ra là phải xây dựng các phương pháp chỉnh hóa (regularization), tìm lời giải xấp xỉ ổn định và xác lập điều kiện duy nhất nghiệm.

Khoảng trống cụ thể mà luận án lựa chọn giải quyết bao gồm hai bài toán ngược phi tuyến:

  1. Bài toán xác định miền cho một phương trình vi phân elliptic cấp hai trong không gian hai chiều $\mathbb{R}^2$ với điều kiện biên Cauchy không đầy đủ.
  2. Bài toán Stefan ngược (inverse Stefan problem) một chiều và hai chiều trong lý thuyết truyền nhiệt chuyển pha.

Mục tiêu nghiên cứu của luận án được xác định rõ qua ba nhiệm vụ cốt lõi:

  1. Xây dựng phương pháp chỉnh hóa và thuật toán tìm lời giải xấp xỉ ổn định cho hai bài toán ngược nêu trên;
  2. Khảo sát và chứng minh tính duy nhất nghiệm của các bài toán nghiên cứu;
  3. Đánh giá sai số giải tích giữa lời giải xấp xỉ ổn định và lời giải chính xác khi dữ liệu đầu vào đo được bị sai lệch với sai số thực nghiệm $\varepsilon$.

Đố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 loại elliptic bậc hai dạng $\Delta u - cu = 0$ ($c \ge 0$), phương trình biến phân liên kết, bài toán Stefan truyền nhiệt với biên tự do, phương trình tích phân loại một và phương trình dạng tích chập.
  • Phạm vi không gian: Miền giới hạn $G \subset \mathbb{R}^2$ không đơn liên được bao bọc bởi hai đường cong: biên ngoài $\Gamma$ đã biết và biên trong $\gamma$ cần xác định, cả hai đều có dạng hình sao ("Star-shaped") đối với một điểm quy chiếu $O \in \mathbb{R}^2$. Bài toán Stefan ngược được khảo sát trong không gian một chiều ($n=1$) và hai chiều ($n=2$).
  • Phạm vi giải tích hàm: Không gian hàm Sobolev $H^1(G) = W^{1,2}(G)$, $H_0^1(G)$, không gian đối ngẫu $[H^1(G)]^* = \mathcal{L}(H^1(G), \mathbb{R})$, không gian hàm liên tục và khả vi $C(\bar{G})$, $C^1(\bar{G})$, $C^2(G)$, cùng không gian Banach các phép biến đổi $C_b^1(\mathcal{O})$ với các đạo hàm riêng bị chặn.

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

Trong phần dẫn nhập và các chương chuyên đề, tác giả Nguyễn Cam đã hệ thống hóa các hướng nghiên cứu liên quan đến bài toán biên tự do, phương trình biến phân và bài toán ngược không chỉnh:

  • Cơ sở lý thuyết chỉnh hóa bài toán ngược: Luận án viện dẫn các công trình kinh điển từ thập niên 1960 của Tikhonov, Lavrentiev và Lions về các nguyên lý chỉnh hóa bài toán không chỉnh. Tác giả áp dụng trực tiếp tư tưởng hàm điều chuẩn (regularizing functional) và phương pháp Tikhonov để xử lý tính mất ổn định của bài toán Stefan ngược.
  • Lý thuyết phương trình biến phân và biên tự do: Tác giả điểm lại các công trình về bài toán biên tự do chính quy (regular free boundary problem). Khi được xử lý theo hướng thuận, bài toán biên tự do thường được chuyển về dạng bất đẳng thức biến phân hoặc phương trình biến phân trên không gian Sobolev. Luận án trích dẫn công trình của J. F. Rodrigues (1987) [32] làm chỗ dựa phương pháp luận để chuyển phương trình elliptic với điều kiện biên hỗn hợp sang phương trình biến phân dạng song tuyến tính và thiết lập tính bức (coercivity).
  • Các định lý giải tích hàm và giải tích hiện đại: Tác giả sử dụng các kết quả giải tích hàm thực nghiệm từ tài liệu của H. Brezis (1983) [5], bao gồm bất đẳng thức Poincaré, công thức đổi biến số cho các hàm thuộc không gian Sobolev dưới tác động của vi đồng phôi lớp $C^1$, và định lý về vết (trace theorem) của Sobolev.
  • Tính duy nhất nghiệm của bài toán Cauchy đối với phương trình elliptic: Tác giả điểm lại các kết quả nghiên cứu về tính duy nhất của bài toán Cauchy cho toán tử elliptic bậc hai của các nhà toán học quốc tế, tiêu biểu là định lý của M. H. Protter (được trích lược trong văn bản dưới ký hiệu E. Protter / E. Rodrigues [32]), khẳng định rằng nếu nghiệm của phương trình elliptic thuần nhất triệt tiêu cùng đạo hàm pháp tuyến trên một phần biên mở $S$ thì nghiệm phải đồng nhất bằng 0 trên toàn bộ miền thông suốt.
  • Các hướng nghiên cứu ứng dụng thực tế: Luận án liên hệ cấu trúc toán học của bài toán xác định miền với các bài toán vật lý kỹ thuật cụ thể đã được công bố trong y văn:
    1. Bài toán xác định hình dáng bình chứa khi biết hình dáng của plasma, vốn là bài toán ngược của bài toán xác định hình dạng plasma trong vật lý lò phản ứng nhiệt hạch (viện dẫn tài liệu [12], [14]);
    2. Bài toán xác định hình dạng và vị trí vết nứt trong môi trường đàn hồi tuyến tính bằng phương pháp điện thế bề mặt (electric potential method) (viện dẫn tài liệu [25]);
    3. Bài toán dẫn nhiệt ngược: tìm phân bố nhiệt độ và biên dạng trong $\gamma$ nơi nhiệt độ đạt trạng thái đẳng nhiệt $u|\gamma = 0$ khi chỉ đo được nhiệt độ $u|\Gamma$ và thông lượng nhiệt $\frac{\partial u}{\partial n}|_{\Gamma_0}$ trên biên ngoài.

Vị trí của luận án được xác lập thông qua việc giải quyết tính không chỉnh theo nghĩa Hadamard của bài toán xác định miền. Thay vì khảo sát bài toán biên tự do thuận thông thường, luận án nghiên cứu trường hợp dữ liệu biên bị thiếu trên một phần biên và dư thừa (điều kiện Cauchy) trên một phần biên khác. Luận án giải quyết khoảng trống này bằng cách đưa bài toán hình học biên tự do về bài toán tìm hệ số biến phân trên một miền chuẩn cố định thông qua ánh xạ tọa độ, biến đổi bài toán thành việc tối ưu hóa một phiếm hàm liên tục trên tập compact hữu hạn chiều, từ đó chứng minh tính ổn định và tính duy nhất của cặp nghiệm $(G, u)$.

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

Luận án vận dụng hệ thống công cụ giải tích hàm giải tích biến phân và lý thuyết phương trình đạo hàm riêng. Các cấu trúc lý thuyết, định lý và phương pháp thực tế được triển khai đồng bộ như sau:

Nhóm công cụ Khái niệm / Định lý sử dụng trong văn bản Nguồn / Căn cứ văn bản Vai trò trong luận án
Không gian hàm $H^1(\Omega) = W^{1,2}(\Omega)$, $H_0^1(\Omega)$, $L^2(\Gamma)$, $C_b^1(\mathcal{O})$ Chuẩn Sobolev và chuẩn Banach $C_b^1$ Xác lập không gian năng lượng cho nghiệm trạng thái và không gian tham số cho miền biến hình
Bất đẳng thức giải tích Bất đẳng thức Poincaré; Mở rộng bất đẳng thức Poincaré cho miền có biên mở $\Gamma_0$ Tài liệu [5], Rodrigues (1987) [32] Chứng minh tính tương đương giữa chuẩn Sobolev $|u|{H^1}$ và chuẩn gradient $|\nabla u|{L^2}$ trên không gian con đóng
Định lý vết Toán tử vết $R: H^1(\Omega) \to L^2(\partial\Omega)$ liên tục duy nhất Định lý giải tích Sobolev Xác định giá trị biên của hàm suy rộng trên $\Gamma$ và $\Gamma_0$
Phương trình biến phân Định lý Lax - Milgram Rodrigues [32] Chứng minh sự tồn tại và duy nhất nghiệm của phương trình biến phân elliptic $a^\Psi(v, w-v) = \langle L, w-v \rangle$
Giải tích phi tuyến Định lý ánh xạ ngược, Định lý điểm bất động Banach Lý thuyết giải tích Banach Chứng minh tính khả đảo địa phương của vi đồng phôi và tính co của ánh xạ trong các lân cận
Lý thuyết xấp xỉ Bổ đề Leray và hệ quả xấp xỉ tập compact Giải tích hàm Xây dựng dãy toán tử chiếu liên tục $P_m: K \to E_m$ xấp xỉ tập compact $K$ bằng các không gian hữu hạn chiều
Công thức vi tích phân Công thức Green và Công thức Green mở rộng cho hàm $H^1$ Rodrigues [32] Biến đổi tích phân từng phần và chứng minh tính duy nhất nghiệm

Phương pháp nghiên cứu thực tế của luận án:

  • Phương pháp biến đổi vi đồng phôi (Mapping method): Luận án cố định hóa miền biến đổi $G$ bằng cách coi $G$ là ảnh của một miền chuẩn $\Omega$ cố định thông qua phép song ánh vi phôi $\Psi: \Omega \to G$, với $\Psi \in C_b^1(\mathcal{O})$, thỏa mãn $\Psi(x, y) = (x, y)$ trên biên ngoài $\Gamma$. Kỹ thuật này chuyển bài toán biên tự do phi tuyến phức tạp thành bài toán tìm hàm biến đổi $\Psi$ xuất hiện trong các hệ số của toán tử vi phân.
  • Phương pháp biến phân và tối ưu hóa phiếm hàm (Variational & Optimization method): Dẫn xuất phương trình trạng thái về dạng biến phân trên $\Omega$. Tác giả chứng minh dạng song tuyến tính $a^\Psi(u, v)$ có tính chất bức (coercive) trên không gian con $W = {v \in H^1(\Omega): v = 0 \text{ trên } \partial\Omega \setminus \Gamma}$, từ đó suy ra sự tồn tại duy nhất của nghiệm trạng thái $v_\Psi$. Tiếp đó, luận án thiết lập phiếm hàm mục tiêu bình phương tối thiểu: $$J(\Psi) = \int_{\Gamma_0} (v_\Psi|_{\Gamma_0} - f)^2 d\Gamma$$ và đưa bài toán xác định miền về bài toán tìm cực trị của $J(\Psi)$ trên tập compact $K \subset C_b^1(\mathcal{O})$.
  • Phương pháp xấp xỉ Galerkin và rời rạc hóa hữu hạn chiều: Để xây dựng thuật toán ổn định, không gian hàm trạng thái $V$ được xấp xỉ bằng dãy không gian con affin hữu hạn chiều ${V_n}$, còn tập compact $K$ được xấp xỉ bởi các tập hữu hạn chiều $K_l = I_l(K) \subset S_l$ thông qua bổ đề Leray. Bài toán cực trị ban đầu được thay thế bằng bài toán cực trị hữu hạn chiều: tìm $\Psi_{nl} \in K_l$ làm cực tiểu $J_n(\Psi)$.
  • Phương pháp chỉnh hóa Tikhonov và tích phân tích chập: Đối với bài toán Stefan ngược, phương pháp tiếp cận bao gồm việc thiết lập các phương trình truyền nhiệt với điều kiện biên di động, chuyển đổi bài toán về hệ phương trình tích phân Fredholm loại một, đưa về dạng tích chập và áp dụng toán tử chỉnh hóa Tikhonov để khắc phục tính không chỉnh.

Nguồn dữ liệu của luận án là dữ liệu toán học giải tích: các hàm dữ liệu biên cho trước $f \in L^2(\Gamma)$, $h \in L^2(\Gamma_0)$ và sai số đo đạc mô phỏng $\varepsilon$ trên các đại lượng này.

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

Luận án có kết cấu gồm phần Mở đầu, phần Ký hiệu và Một số kết quả sử dụng, tiếp đến là hai phần nội dung chính:

Phần Một: Bài toán xác định miền cho một phương trình Elliptic

Phần Một là trọng tâm nghiên cứu của luận án, được chia thành hai mục lớn gồm Sự chỉnh hóa và Sự duy nhất nghiệm.

I. Sự chỉnh hóa

Tác giả phát biểu bài toán toán tử: Tìm miền giới hạn không đơn liên $G \subset \mathbb{R}^2$ và hàm trạng thái $u$ xác định trên $G$ thỏa mãn phương trình elliptic: $$\Delta u - cu = 0 \quad \text{trong } G \quad (c(x) \ge 0)$$ cùng các điều kiện biên:

  • Điều kiện Dirichlet trên biên ngoài: $u = f$ trên $\Gamma$;
  • Điều kiện thông lượng Neumann (Cauchy) trên phần biên mở $\Gamma_0 \subset \Gamma$ ($\Gamma_0 \ne \emptyset$): $\frac{\partial u}{\partial n} = h$ trên $\Gamma_0$;
  • Điều kiện Dirichlet đồng nhất trên biên trong cần tìm $\gamma = \partial G \setminus \Gamma$: $u = 0$ trên $\gamma$.

Để giải quyết tính không chỉnh, tác giả thiết lập mô hình hình học: Miền $G$ giới hạn bởi hai đường cong $\Gamma$ và $\gamma$ thuộc lớp $C^1$, đều có tính chất "Star-shaped" đối với điểm gốc $O$. Biểu diễn trong hệ tọa độ cực $(\theta, \rho)$:

  • $\Gamma = {(\theta, R(\theta)): \theta \in [0, 2\pi]}$
  • $\gamma = {(\theta, r(\theta)): \theta \in [0, 2\pi]}$
  • Miền chuẩn cố định $\Omega$ được giới hạn bởi $\Gamma$ và một đường tròn $\gamma_0 = {(\theta, r_0(\theta)): \theta \in [0, 2\pi]}$ nằm phía trong $\Gamma$.

Phép vi phôi $\Psi$ biến đổi $\Omega$ thành $G$ được xây dựng tường minh qua hệ thức bán kính cực: $$\Psi(\theta, \rho) = \left(\theta, r(\theta) + \frac{R(\theta) - r(\theta)}{R(\theta) - r_0(\theta)} [\rho - r_0(\theta)]\right)$$ và phép biến đổi ngược $\phi = \Psi^{-1}$ khả tích lớp $C^1$.

Bằng công thức đổi biến số cho hàm Sobolev, bài toán được quy về dạng biến phân trên miền cố định $\Omega$: $$a^\Psi(v, w - v) = \langle L, w - v \rangle \quad \forall w \in V, \ v \in V$$ trong đó dạng song tuyến tính $a^\Psi(v, w) = \int_\Omega \left( \sum_{i,j=1}^2 a_{ij}^\Psi \frac{\partial v}{\partial x_i}\frac{\partial w}{\partial x_j} + c^\Psi v w \right) dx$ phụ thuộc trực tiếp vào ma trận đạo hàm riêng và định thức Jacobi của $\Psi$. Nhờ bất đẳng thức Poincaré mở rộng (Rodrigues [32]), tác giả chứng minh $a^\Psi$ là bức trên không gian con $W = {v \in H^1(\Omega): v = 0 \text{ trên } \partial\Omega \setminus \Gamma}$, suy ra phương trình biến phân có nghiệm duy nhất $v_\Psi \in V$.

Tiếp đó, tác giả đưa ra chuỗi 5 bổ đề kỹ thuật phục vụ việc chứng minh sự hội tụ:

  • Bổ đề 1: Cho dãy ${\Psi_k} \subset Q$, nếu $\Psi_k \to \Psi$ trong chuẩn $C_b^1(\mathcal{O})$ thì các hệ số ma trận $a_{ij}^{\Psi_k}$ hội tụ đều về $a_{ij}^\Psi$ trên $\bar{\Omega}$.
  • Bổ đề 2: Đạo hàm riêng của ánh xạ ngược $(D_i \phi_j^k) \circ \Psi_k$ hội tụ đều về $(D_i \phi_j) \circ \Psi$ trên $\bar{\Omega}$.
  • Bổ đề 3 & Bổ đề 4: Chứng minh tính ổn định hình học của lân cận: tồn tại tập mở $E$ sao cho $\Psi_k(\bar{\Omega}) \subset E \subset \Psi_k(\mathcal{O})$ với mọi chỉ số $k$ đủ lớn.
  • Bổ đề 5: Dựa trên định lý giá trị trung bình và định lý điểm bất động cho ánh xạ co (hệ số co $k = 1/2$), chứng minh các lân cận cầu mở được bảo toàn qua dãy ánh xạ vi phôi khả đảo.

Các kết quả then chốt của mục Sự chỉnh hóa được đúc kết qua ba định lý toán học:

  • Định lý 1: Khi $\Psi_k \to \Psi$ trong $C_b^1(\mathcal{O})$ và số chiều xấp xỉ Galerkin $n_k \to \infty$, nghiệm xấp xỉ hữu hạn chiều $v_k = v_{n_k}^{\Psi_k}$ hội tụ mạnh trong chuẩn Sobolev $W$ về nghiệm chính xác $v_\Psi$.
  • Định lý 2 và Định lý 2': Chứng minh bài toán cực trị rời rạc trên không gian hữu hạn chiều luôn giải được; dãy nghiệm xấp xỉ ${\Psi_{nl}}$ chứa dãy con hội tụ về nghiệm cực trị $\Psi_0$ của bài toán ban đầu; khoảng cách Hausdorff giữa tập nghiệm xấp xỉ $A_{nl}$ và tập nghiệm lý thuyết $A$ tiến dần về 0: $$\lim_{n, l \to \infty} \rho(A_{nl}, A) = 0$$
  • Định lý 3 (Tính ổn định của phương pháp chỉnh hóa): Khi dữ liệu đo biên bị sai số nhưng hội tụ về dữ liệu chính xác: $f \to f_0$ trong chuẩn $L^2(\Gamma)$ và $h \to h_0$ trong chuẩn $L^2(\Gamma_0)$, tập nghiệm cực tiểu $A^{fh}$ và tập nghiệm rời rạc $A_{nl}^{fh}$ phụ thuộc liên tục vào cặp dữ liệu $(f, h)$ và hội tụ về nghiệm chuẩn $(f_0, h_0)$ theo độ đo khoảng cách tập nghiệm.

II. Sự duy nhất nghiệm

Tác giả khảo sát tính duy nhất nghiệm của bài toán Cauchy cho phương trình vi phân elliptic bậc hai dạng: $$Lu \equiv \sum_{i,j=1}^n a_{jk}(x)\frac{\partial^2 u}{\partial x_j \partial x_k} + \sum_{i=1}^n b_i(x)\frac{\partial u}{\partial x_i} + c(x)u = 0$$ Dựa trên kết quả nền tảng của M. H. Protter, tác giả thiết lập mệnh đề: Nếu $u \in C^2(G) \cap C(\bar{G})$ thỏa mãn đồng thời $u = 0$ và $\frac{\partial u}{\partial n} = 0$ trên một phần biên mở $S \subset \partial G$ thì $u \equiv 0$ trên toàn bộ miền liên thông $G$. Từ đó suy ra tính duy nhất của nghiệm trạng thái $u$ khi miền $G$ đã biết trước.

Đối với bài toán ngược khi biên $\gamma$ chưa biết, luận án xây dựng chứng minh cho tính duy nhất của cả cặp $(G, u)$: Xét hai miền $G_1$ và $G_2$ cùng có biên ngoài $\Gamma$ và hai biên trong tương ứng $\gamma_1, \gamma_2$ đều thỏa mãn tính chất "Star-shaped" đối với cùng một điểm cực $O$. Bằng cách sử dụng công thức Green, nguyên lý cực đại và tính giải tích của nghiệm phương trình elliptic, tác giả chứng minh rằng nếu hai cặp nghiệm $(G_1, u_1)$ và $(G_2, u_2)$ cùng tương thích với một cặp dữ liệu Cauchy $(f, h)$ trên $\Gamma_0$ và điều kiện triệt tiêu trên biên trong thì bắt buộc $G_1 \equiv G_2$ (tức $\gamma_1 \equiv \gamma_2$) và $u_1 \equiv u_2$.

Phần Hai: Bài toán Stefan ngược

Nội dung Phần Hai được tác giả tóm lược trong phần Mở đầu với các bước giải tích cụ thể:

  • Mô hình hóa bài toán: Xét bài toán Stefan ngược một chiều ($n=1$) và hai chiều ($n=2$). Trong bài toán này, biên tự do phân chia pha di động theo thời gian, đại diện cho quá trình đông kết hoặc nóng chảy trong vật lý nhiệt.
  • Chuyển đổi tích phân: Đưa hệ phương trình đạo hàm riêng parabol với điều kiện biên tự do về phương trình tích phân loại một. Sau đó, bằng các phép biến đổi hàm và biến số tích phân, phương trình tích phân loại một được quy đổi thành phương trình tích chập.
  • Chỉnh hóa và đánh giá sai số: Áp dụng nguyên lý điều chuẩn Tikhonov đối với phương trình tích chập nhằm khắc phục sự mất ổn định của toán tử tích phân loại một. Luận án thiết lập công thức nghiệm xấp xỉ ổn định và dẫn xuất giải tích đánh giá sai số giữa nghiệm xấp xỉ chỉnh hóa và nghiệm chính xác khi dữ liệu đo đạc thực nghiệm bị nhiễu với độ lệch sai số $\varepsilon$.

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

Luận án tự đúc kết các đóng góp mới về mặt lý luận giải tích hàm và toán ứng dụng:

  1. Xây dựng quy trình chỉnh hóa biến phân cho bài toán xác định miền elliptic: Thiết lập thành công phương pháp vi phôi để chuyển một bài toán biên tự do không chỉnh thành bài toán tìm hệ số trong phương trình biến phân elliptic trên một miền chuẩn cố định;
  2. Thiết lập tính giải được và tính hội tụ hữu hạn chiều: Chứng minh phiếm hàm bình phương tối thiểu $J(\Psi)$ đạt cực trị trên tập compact $K \subset C_b^1(\mathcal{O})$; chứng minh tính hội tụ của sơ đồ xấp xỉ Galerkin kết hợp bổ đề Leray từ các không gian hữu hạn chiều về nghiệm của bài toán liên tục;
  3. Chứng minh tính ổn định của tập nghiệm: Xác lập Định lý 3 chứng minh sự phụ thuộc liên tục của tập lời giải xấp xỉ theo dữ liệu biên bị nhiễu $f$ và $h$ trong chuẩn không gian Hilbert $L^2$, giải quyết tính không chỉnh theo nghĩa Hadamard;
  4. Khẳng định tính duy nhất của cặp nghiệm $(G, u)$: Chứng minh toán học chặt chẽ tính duy nhất đồng thời của hình dạng miền $G$ và hàm thế $u$ dưới giả thiết biên hình sao (Star-shaped), mở rộng phạm vi áp dụng so với các bài toán Cauchy truyền thống vốn chỉ xét tính duy nhất của $u$ trên miền cố định;
  5. Giải thuật giải tích cho bài toán Stefan ngược: Đưa thành công bài toán Stefan ngược một chiều và hai chiều về phương trình tích chập, áp dụng phương pháp Tikhonov xây dựng nghiệm ổn định và thu được đánh giá sai số tường minh theo cấp của sai số đo $\varepsilon$.

Về mặt ứng dụng và giải pháp kỹ thuật, luận án cung cấp cơ sở toán học để giải quyết ba bài toán thực tế:

  • Bài toán chuẩn đoán hình học vết nứt trong cấu kiện đàn hồi mà không cần phá hủy mẫu thông qua phương pháp điện thế bề mặt;
  • Bài toán thiết kế vỏ bình chứa bao quanh plasma trong vật lý lò phản ứng hạt nhân;
  • Bài toán xác định bề mặt đẳng nhiệt bên trong các vật thể truyền nhiệt khi chỉ có thể tiếp cận đo đạc thông lượng nhiệt và nhiệt độ ở vỏ biên ngoài.

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

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

  • Giả thiết hình học khắt khe: Toàn bộ các chứng minh về phép biến đổi vi phôi $\Psi$ và tính duy nhất nghiệm của miền $G$ đều phụ thuộc vào điều kiện biên $\Gamma$ và $\gamma$ là các đường cong "Star-shaped" đối với một điểm gốc quy chiếu $O$. Khi biên có hình học phức tạp, gấp nếp hoặc đa liên với nhiều đảo biên bên trong, phép đặt ẩn tọa độ cực không còn áp dụng được trực tiếp.
  • Giả thiết về tính compact: Việc chứng minh tồn tại nghiệm của bài toán cực trị đòi hỏi tập các biến đổi $K$ phải là tập compact nằm trong không gian Banach $C_b^1(\mathcal{O})$, điều này giới hạn lớp nghiệm miền trong một tập compact định trước.
  • Dạng phương trình elliptic: Phương trình trạng thái trong Phần Một được khảo sát với hệ số $c(x) \ge 0$ và toán tử Laplace chính quy, chưa mở rộng sang các hệ phương trình hoặc toán tử elliptic suy biến.
  • Phạm vi tính toán số: Luận án tập trung toàn diện vào các định lý giải tích hàm, chứng minh sự tồn tại, duy nhất, tính hội tụ và chặn sai số lý thuyết, văn bản chưa bao gồm các bảng tính số mô phỏng trên máy tính điện tử đối với các dữ liệu thực nghiệm cụ thể.

Hướng nghiên cứu tiếp suy biến từ phạm vi luận án:

  • Mở rộng lý thuyết xác định miền cho các hình học biên tổng quát không có tính chất Star-shaped;
  • Khảo sát bài toán xác định miền cho phương trình elliptic trong không gian ba chiều $\mathbb{R}^3$;
  • Xây dựng thuật toán số cụ thể và thử nghiệm tính toán trên máy tính cho phương pháp xấp xỉ Galerkin rời rạc kết hợp chỉnh hóa Tikhonov.

Giá trị tham khảo

Luận án là tài liệu chuyên khảo có giá trị học thuật chuyên sâu đối với các nhóm đối tượng nghiên cứu:

  • Nghiên cứu sinh và học viên cao học ngành Toán giải tích, Toán ứng dụng: Tham khảo kỹ thuật xử lý bài toán biên tự do bằng phương pháp biến đổi miền về hệ số phương trình biến phân; vận dụng bất đẳng thức Poincaré mở rộng của Rodrigues và bổ đề Leray để xấp xỉ tập compact trên không gian Banach.
  • Giảng viên và nhà nghiên cứu phương trình vi phân đạo hàm riêng: Tham khảo chuỗi bổ đề kỹ thuật (Bổ đề 1 đến Bổ đề 5) về sự hội tụ của vi đồng phôi và ma trận hệ số nghịch đảo, cũng như phương pháp luận chứng minh tính duy nhất của cặp nghiệm miền - hàm thế $(G, u)$.
  • Các nhà nghiên cứu trong lĩnh vực cơ học môi trường liên tục, địa vật lý và vật lý kỹ thuật: Phần Một cung cấp khung phân tích toán học cho bài toán xác định khuyết tật, vết nứt ngầm bằng phương pháp điện thế bề mặt và bài toán cân bằng biên plasma. Phần Hai cung cấp sơ đồ chuyển đổi bài toán truyền nhiệt chuyển pha Stefan về phương trình tích phân tích chập để xử lý bằng phương pháp Tikhonov kèm đánh giá sai số.

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

1. Tại sao bài toán xác định miền trong luận án được xếp vào loại bài toán không chỉnh theo nghĩa Hadamard?
Bài toán xác định miền cho phương trình elliptic trong luận án có điều kiện biên dạng Cauchy trên một phần biên mở $\Gamma_0$ (biết cả giá trị hàm $u=f$ và đạo hàm pháp tuyến $\frac{\partial u}{\partial n}=h$) nhưng lại khuyết hoàn toàn thông tin biên trên phần biên trong $\gamma$. Đây là bài toán Cauchy đối với phương trình elliptic, vốn không đảm bảo sự tồn tại nghiệm với dữ liệu đo bất kỳ, và khi có nghiệm thì nghiệm không phụ thuộc liên tục vào dữ liệu $f, h$. Một sai số nhỏ trong phép đo trên $\Gamma_0$ có thể gây ra biến dạng lớn hoặc mất nghiệm của biên $\gamma$.

2. Phép biến đổi vi phôi $\Psi$ được xây dựng bằng cách nào để đưa bài toán biên tự do về miền cố định?
Tác giả tận dụng giả thiết biên ngoài $\Gamma$ và biên trong $\gamma$ là các đường cong "Star-shaped" đối với điểm gốc $O$. Trong tọa độ cực $(\theta, \rho)$, các biên được tham số hóa thành các hàm bán kính $R(\theta)$ và $r(\theta)$ thuộc lớp $C^1$. Tác giả chọn một miền chuẩn cố định $\Omega$ có biên trong là đường tròn bán kính $r_0(\theta)$ và thiết lập ánh xạ tuyến tính theo phương xuyên tâm: $$\Psi(\theta, \rho) = \left(\theta, r(\theta) + \frac{R(\theta) - r(\theta)}{R(\theta) - r_0(\theta)} [\rho - r_0(\theta)]\right)$$ Ánh xạ này là một song ánh lớp $C^1$, giữ nguyên các điểm trên biên $\Gamma$ và biến miền cố định $\Omega$ thành miền biến thiên $G$.

3. Phiếm hàm mục tiêu $J(\Psi)$ được thiết lập như thế nào và việc tìm nghiệm xấp xỉ được thực hiện ra sao?
Phiếm hàm mục tiêu được định nghĩa là độ lệch bình phương giữa vết nghiệm biến phân $v_\Psi$ trên phần biên $\Gamma_0$ và dữ liệu đo thực tế $f$: $$J(\Psi) = \int_{\Gamma_0} (v_\Psi|{\Gamma_0} - f)^2 d\Gamma$$ Nghiệm hình học $\Psi_0$ được tìm thông qua bài toán cực tiểu hóa $J(\Psi)$ trên tập compact $K \subset C_b^1(\mathcal{O})$. Để giải bài toán này, tác giả rời rạc hóa không gian Sobolev bằng phương pháp Galerkin ${V_n}$ và xấp xỉ tập compact $K$ bằng các tập compact hữu hạn chiều $K_l = I_l(K)$ dựa trên bổ đề Leray, tạo thành bài toán cực trị hữu hạn chiều $J_n(\Psi{nl}) \to \min$.

4. Nội dung bài toán Stefan ngược trong Phần Hai được giải quyết bằng các kỹ thuật toán học nào?
Trong Phần Hai, bài toán Stefan ngược một chiều và hai chiều được chuyển đổi từ hệ phương trình đạo hàm riêng mô tả quá trình truyền nhiệt với biên pha di động thành phương trình tích phân loại một. Sau đó, tác giả thực hiện phép đổi biến để đưa phương trình tích phân về phương trình dạng tích chập. Cuối cùng, phương pháp chỉnh hóa Tikhonov được sử dụng để xây dựng nghiệm xấp xỉ ổn định và xác lập công thức chặn sai số khi dữ liệu đo có sai số $\varepsilon$.

5. Điều kiện hình học "Star-shaped" đóng vai trò then chốt nào trong toàn bộ luận án?
Điều kiện "Star-shaped" (đường cong hình sao đối với điểm gốc $O$) đảm bảo rằng mỗi tia kẻ từ gốc $O$ chỉ cắt biên ngoài $\Gamma$ và biên trong $\gamma$ tại đúng một điểm duy nhất. Tính chất này cho phép biểu diễn các biên dưới dạng hàm số đơn trị tuần hoàn của góc cực $\theta \in [0, 2\pi]$ trong tọa độ cực, từ đó cho phép xây dựng tường minh phép vi phôi $\Psi$ và nghịch đảo $\phi = \Psi^{-1}$, đồng thời là nền tảng hình học để áp dụng công thức Green chứng minh tính duy nhất của cặp nghiệm $(G, u)$.

Kết luận

Luận án Phó Tiến sĩ Toán Lý của tác giả Nguyễn Cam đã giải quyết trọn vẹn về mặt giải tích bài toán xác định miền cho phương trình vi phân elliptic và bài toán Stefan ngược thông qua các phương pháp chỉnh hóa hiện đại. Bằng việc kết hợp phương pháp vi đồng phôi hình học, lý thuyết phương trình biến phân trên không gian Sobolev và phương pháp xấp xỉ compact hữu hạn chiều, công trình đã chứng minh thành công tính giải được, tính ổn định và tính duy nhất của nghiệm bài toán xác định miền. Các kết quả giải tích của luận án đóng góp cơ sở lý thuyết chuẩn xác cho việc mô hình hóa và xử lý các bài toán ngược biên tự do trong vật lý kỹ thuật, truyền nhiệt và cơ học phá hủy đàn hồi.