Tổng quan luận án
Luận án tiến sĩ toán học của tác giả Trần Thị Quỳnh Chi, với đề tài "Dáng điệu tiệm cận nghiệm của một số phương trình và hệ phương trình parabolic phi tuyến", được thực hiện tại Trường Đại học Sư phạm Hà Nội dưới sự hướng dẫn khoa học của GS. TS. Cung Thế Anh, thuộc chuyên ngành Phương trình vi phân và tích phân (mã số: 9 46 01 03), bảo vệ năm 2023.
Tính cấp thiết và khoảng trống nghiên cứu
Các phương trình đạo hàm riêng (PDEs) loại parabolic xuất hiện nhiều trong việc mô hình hóa các quá trình vật lý và sinh học, tiêu biểu như các mô hình biến động quần thể, các quá trình truyền nhiệt và hiện tượng khuếch tán. Sau khi xác lập tính đặt đúng (well-posedness) của bài toán, việc nghiên cứu dáng điệu tiệm cận của nghiệm khi thời gian tiến ra vô cùng ($t \to +\infty$) giữ vai trò quan trọng nhằm xác định trạng thái giới hạn và dự đoán xu thế vận động dài hạn của các hệ động lực. Do không gian trạng thái của các phương trình đạo hàm riêng là vô hạn chiều, lý thuyết tập hút (attractor theory) được sử dụng làm công cụ trung tâm để nắm bắt thông tin tiệm cận.
Từ việc điểm luận các kết quả đã có, tác giả xác định ba khoảng trống nghiên cứu chính:
- Đối với phương trình parabolic suy biến tựa tuyến tính chứa toán tử dạng $p$-Laplace có trọng, phần lớn các công trình trước đây chỉ khảo sát với số hạng phi tuyến $f(u)$ bị khống chế bởi điều kiện tăng trưởng đa thức kiểu Sobolev. Vấn đề nghiên cứu sự tồn tại của tập hút toàn cục khi phần phi tuyến có bậc tăng trưởng tùy ý vẫn là một bài toán mở.
- Đối với phương trình parabolic suy biến tựa tuyến tính không ôtônôm (chứa tham số thời gian trong số hạng ngoại lực $g(x,t)$) xét trên toàn không gian $\mathbb{R}^N$, các kết quả nghiên cứu trước đó đòi hỏi hàm trọng $a(x)$ phải thỏa mãn điều kiện $(H_{\alpha,\beta})$, tức là $a(x)$ phải tăng nhanh hơn $|x|^\beta$ khi $|x| \to \infty$ để bảo đảm tính compắc của phép nhúng Sobolev có trọng. Việc thay thế điều kiện hạn chế dáng điệu tại vô cùng bằng một điều kiện khả tích địa phương mới nhằm bảo đảm tính compắc tiệm cận lùi trên toàn miền $\mathbb{R}^N$ là bài toán chưa có lời giải.
- Đối với hệ phương trình phản ứng khuếch tán, các kết quả hiện có về tập hút chủ yếu dừng lại ở trường hợp các số hạng phi tuyến thỏa mãn tăng trưởng đa thức hoặc tiêu hao kiểu Sobolev. Khi số hạng phi tuyến mang đặc tính tăng trưởng kiểu mũ (chẳng hạn $f(u) \sim e^u$), các điều kiện đa thức kinh điển không còn thỏa mãn, đòi hỏi phải thiết lập các đánh giá giải tích mới.
Mục tiêu nghiên cứu
Luận án xác định các mục tiêu nghiên cứu cụ thể sau:
- Thiết lập tính đặt đúng và chứng minh sự tồn tại của tập hút toàn cục trong các không gian hàm $L^2(\Omega)$ và không gian Sobolev có trọng $M_0^{1,p}(\Omega, a)$ đối với một lớp phương trình parabolic suy biến tựa tuyến tính chứa toán tử $p$-Laplace với phần phi tuyến tăng trưởng tùy ý trên miền bị chặn.
- Thiết lập tính đặt đúng và chứng minh sự tồn tại của tập hút lùi (pullback attractor) trong các không gian $L^2(\mathbb{R}^N)$, $L^q(\mathbb{R}^N)$ và $M_0^{1,p}(\mathbb{R}^N, a)$ đối với lớp phương trình parabolic suy biến tựa tuyến tính không ôtônôm trên toàn không gian $\mathbb{R}^N$ với điều kiện mới về hệ số khuếch tán $a(x)$.
- Nghiên cứu tính đặt đúng, sự tồn tại tập hút toàn cục trong $(L^2(\Omega))^2$, cấu trúc tập hút toàn cục đối với hệ gradient, đồng thời chứng minh sự tồn tại và tính ổn định mũ của nghiệm dừng đối với hệ phương trình phản ứng khuếch tán có phần phi tuyến kiểu mũ.
Đối tượng và phạm vi nghiên cứu
- Đối tượng nghiên cứu: Nghiệm yếu và dáng điệu tiệm cận của nghiệm (tập hút toàn cục, tập hút lùi, nghiệm dừng và tính ổn định mũ) của các phương trình và hệ phương trình parabolic phi tuyến.
- Phạm vi nghiên cứu: Luận án tập trung vào ba bài toán giải tích cụ thể:
- Nội dung 1: Lớp phương trình parabolic suy biến tựa tuyến tính tự trị: $$u_t - \text{div}(a(x)|\nabla u|^{p-2}\nabla u) + f(u) = g(x), \quad x \in \Omega, ; t > 0$$ với điều kiện biên Dirichlet thuần nhất $u(x,t) = 0$ trên $\partial\Omega$ và điều kiện ban đầu $u(x,0) = u_0(x) \in L^2(\Omega)$, trong đó $\Omega \subset \mathbb{R}^N$ ($N \ge 2$) là miền bị chặn có biên trơn, $2 \le p \le N$, $f$ tăng trưởng tùy ý.
- Nội dung 2: Lớp phương trình parabolic suy biến tựa tuyến tính không ôtônôm: $$\frac{\partial u}{\partial t} - \text{div}(a(x)|\nabla u|^{p-2}\nabla u) + \lambda |u|^{p-2}u + f(u) = g(x,t), \quad x \in \mathbb{R}^N, ; a < t$$ với điều kiện ban đầu $u(x,a) = u_a(x) \in L^2(\mathbb{R}^N)$, trong đó $N \ge 3$, $2 \le p < N$, $\lambda > 0$, phần phi tuyến $f$ tăng trưởng cấp đa thức bậc $q \ge 2$.
- Nội dung 3: Lớp hệ phương trình phản ứng khuếch tán: $$\begin{cases} \frac{\partial u}{\partial t} - \Delta u + f(u,v) = g_1(x), & x \in \Omega, ; t > 0 \ \frac{\partial v}{\partial t} - \Delta v + h(u,v) = g_2(x), & x \in \Omega, ; t > 0 \end{cases}$$ với điều kiện biên Dirichlet thuần nhất trên $\partial\Omega$, điều kiện ban đầu $u(x,0)=u_0(x), v(x,0)=v_0(x)$, trong đó $\Omega \subset \mathbb{R}^N$ là miền bị chặn có biên trơn, ngoại lực $g_1, g_2 \in L^2(\Omega)$, các số hạng phi tuyến $f, h$ mang tính chất tăng trưởng kiểu mũ.
Tổng quan tài liệu và vị trí của luận án
Các hướng nghiên cứu tiền đề
Lý thuyết tập hút toàn cục đối với các hệ động lực vô hạn chiều sinh bởi các phương trình đạo hàm riêng parabolic đã được trình bày có hệ thống trong các chuyên khảo của Babin và Vishik (1992 [12]), Chepyzhov và Vishik (2002 [17]), Temam (1988 [49]), Robinson (2001 [44]). Ban đầu, phần lớn các kết quả đạt được tập trung vào phương trình parabolic vô hướng nửa tuyến tính hoặc tựa tuyến tính không suy biến.
Đối với trường hợp suy biến, các công trình nghiên cứu phát triển qua các mốc:
- P. Caldiroli và R. Musina (2000 [10]) nghiên cứu bài toán elliptic chứa toán tử suy biến $-\text{div}(a(x)\nabla u)$ trên miền $\Omega \subset \mathbb{R}^N$ với hàm trọng $a(x) \ge 0$ đo được, có thể triệt tiêu tại hữu hạn điểm. Nhóm tác giả xây dựng không gian Sobolev có trọng $D_0^1(\Omega; a)$, thiết lập các bất đẳng thức Sobolev có trọng và các định lý nhúng kiểu Sobolev và Rellich làm tiền đề giải tích.
- Kê (2013 [4], 2015 [5]) mở rộng toán tử của Caldiroli và Musina sang trường hợp tựa tuyến tính chứa toán tử $p$-Laplace suy biến $-\text{div}(a(x)|\nabla u|^{p-2}\nabla u)$ trong phương trình $u_t - \text{div}(a|\nabla u|^{p-2}\nabla u) + f(x,u) = g$, xây dựng các không gian Sobolev có trọng tương thích cho trường hợp $p \ne 2$.
- H. Li và cộng sự (2016 [27]) nghiên cứu phương trình parabolic suy biến chứa toán tử $-\text{div}(a(x)\nabla u)$ với điều kiện mở rộng: tập suy biến $\Sigma = {x \in \Omega: a(x) = 0}$ là một tập con đóng có độ đo Lebesgue bằng 0 ($m(\Sigma) = 0$), cho phép tập suy biến chứa vô số điểm. Các tác giả xét số hạng phi tuyến $f(u)$ thỏa mãn đánh giá đa thức hai phía: $$C_1|u|^p - C_0 \le f(u)u \le C_1|u|^p + C_0, \quad 2 \le p$$ và thiết lập sự tồn tại của tập hút toàn cục trong $L^p(\Omega)$ và $H_0^{1,a}(\Omega)$.
- Chai (2019 [41]) tiếp tục mở rộng bài toán trên với điều kiện tích phân yếu hơn $\int_\Omega [a(x)]^{-1/\gamma} dx < \infty$ ($\gamma \in (0, p-1)$), chỉ ra sự tồn tại của nghiệm entropy và tập hút toàn cục trong không gian $L^1(\Omega)$.
- Về phương trình không ôtônôm và tập hút lùi: Khi số hạng ngoại lực phụ thuộc thời gian $g(x,t)$, hệ chuyển thành quá trình hai tham số $U(t,a)$. Khái niệm tập hút lùi được phát triển trong các nghiên cứu [28, 31, 34, 35, 46, 47, 54] và các bài báo liên quan [1, 5, 23, 24]. Trên miền không bị chặn, các công trình trước đây đều phải áp đặt giả thiết $(H_{\alpha,\beta})$: $$\liminf_{|x|\to\infty} \frac{a(x)}{|x|^\beta} > 0 \quad \text{với } \beta > p + \frac{N}{2}(p-2)$$
- Về hệ phương trình phản ứng khuếch tán: Các công trình [13, 16, 25, 37, 38, 52, 58] khảo sát sự tồn tại và tính chất của tập hút khi các hàm phi tuyến thỏa mãn điều kiện tăng trưởng đa thức và tiêu hao kiểu Sobolev. Khi $f(u) = e^u$, hàm tăng trưởng nhanh hơn bất kỳ đa thức nào nên các kết quả này không thể áp dụng trực tiếp.
Khoảng trống luận án lựa chọn giải quyết
Luận án định vị nghiên cứu vào ba điểm then chốt:
- Giải quyết bài toán phương trình parabolic suy biến tựa tuyến tính chứa toán tử $p$-Laplace có trọng trong trường hợp hàm phi tuyến $f(u)$ có tốc độ tăng trưởng tùy ý (chỉ yêu cầu điều kiện tiêu hao dưới dạng vi phân $f'(s) \ge -A$ và đại số $f(s)s \ge \mu |s|^p - c_1$), khắc phục giới hạn tăng trưởng đa thức của các công trình [27, 41].
- Thay thế điều kiện hạn chế hành vi tại vô cùng $(H_{\alpha,\beta})$ bằng một giả thiết mới $(H_a)$ đòi hỏi tính khả tích địa phương bậc cao của hàm nghịch đảo hệ số khuếch tán trên các hình vành khăn, thiết lập tập hút lùi trên toàn miền $\mathbb{R}^N$ khi không có phép nhúng Sobolev compắc thông thường.
- Mở rộng lý thuyết tập hút toàn cục và nghiên cứu tính ổn định của nghiệm dừng đối với hệ phương trình phản ứng khuếch tán có số hạng phi tuyến kiểu mũ (exponential growth).
Cơ sở lý thuyết và phương pháp nghiên cứu
Cơ sở lý thuyết và các không gian hàm
Luận án vận dụng lý thuyết giải tích hàm phi tuyến, lý thuyết không gian Sobolev có trọng và lý thuyết hệ động lực vô hạn chiều:
-
Không gian Sobolev có trọng trên miền bị chặn $M_0^{1,p}(\Omega, a)$: Cho $\Omega \subset \mathbb{R}^N$ ($N \ge 2$) là miền bị chặn, $a \in L_{loc}^1(\Omega)$, $a(x) > 0$ hầu khắp nơi trừ tập đóng $\Sigma$ có độ đo $m(\Sigma) = 0$. Giả thiết tồn tại $\alpha \in (0, p)$ sao cho $\int_\Omega [a(x)]^{-N/\alpha} dx < \infty$. Không gian $M_0^{1,p}(\Omega, a)$ được định nghĩa là bao đóng của $C_0^\infty(\Omega)$ theo chuẩn: $$|u|{M_0^{1,p}(\Omega, a)} = \left(\int\Omega a(x)|\nabla u|^p dx\right)^{1/p}$$ Không gian đối ngẫu tương ứng được ký hiệu là $M^{-1,p'}(\Omega, a)$ với $\frac{1}{p} + \frac{1}{p'} = 1$.
-
Không gian Sobolev có trọng trên toàn không gian $M_0^{1,p}(\mathbb{R}^N, a)$: Là phần bổ sung đủ của $C_0^\infty(\mathbb{R}^N)$ theo chuẩn: $$|u|{M_0^{1,p}(\mathbb{R}^N, a)} = \left(\int{\mathbb{R}^N} \left(a(x)|\nabla u|^p + |u|^p\right) dx\right)^{1/p}$$
-
Lý thuyết tập hút toàn cục: Xét nửa nhóm liên tục ${S(t)}_{t \ge 0}$ trên không gian pha Banach $X$. Tập con $\mathcal{A} \subset X$ là tập hút toàn cục nếu thỏa mãn:
- $\mathcal{A}$ đóng và bị chặn trong $X$;
- $\mathcal{A}$ bất biến hoàn toàn: $S(t)\mathcal{A} = \mathcal{A}, ; \forall t \ge 0$;
- $\mathcal{A}$ hút mọi tập bị chặn $B \subset X$: $\lim_{t \to +\infty} \text{dist}(S(t)B, \mathcal{A}) = 0$. Luận án áp dụng định lý: Một hệ động lực tiêu hao và compắc tiệm cận thì tồn tại duy nhất một tập hút toàn cục kết nối $\mathcal{A} = \omega(B_0)$, với $B_0$ là tập hấp thụ bị chặn.
-
Lý thuyết tập hút lùi (Pullback Attractor): Xét quá trình hai tham số ${U(t, a)}{t \ge a, a \in \mathbb{R}}$ trên không gian Banach $X$ thỏa mãn $U(t,s)U(s,a) = U(t,a)$ và $U(a,a) = \text{Id}$. Họ tập bị chặn $\hat{\mathcal{A}} = {A(t)}{t \in \mathbb{R}}$ là tập hút lùi duy nhất nếu $A(t)$ compắc, bất biến ($U(t,a)A(a) = A(t)$) và hút lùi mọi tập bị chặn $D \subset X$ theo nghĩa: $$\lim_{a \to -\infty} \text{dist}(U(t,a)D, A(t)) = 0$$ Tính compắc tiệm cận lùi được kiểm tra thông qua điều kiện $\omega$-compắc giới hạn lùi sử dụng độ đo không compắc Kuratowski: $$\alpha(B) = \inf{\delta > 0 : B \text{ có một phủ mở hữu hạn gồm các tập có đường kính } \le \delta}$$
Phương pháp nghiên cứu
- Phương pháp xấp xỉ Galerkin: Xây dựng hệ nghiệm xấp xỉ hữu hạn chiều dựa trên hệ hàm riêng trực chuẩn của toán tử trong $L^2$. Sử dụng định lý Peano để giải hệ phương trình vi phân thường xấp xỉ.
- Phương pháp compắc và đánh giá tiên nghiệm: Thiết lập các chặn giải tích đều theo tham số thời gian và thứ tự xấp xỉ trong các không gian $L^2$, $L^p$ và không gian Sobolev có trọng. Sử dụng bổ đề compắc Aubin-Lions-Simon (Định lý 1.5 của luận án) để trích dãy con hội tụ mạnh.
- Kỹ thuật cắt Nemytskii: Do $f(u)$ tăng trưởng tùy ý, nghiệm $u$ không thể dùng trực tiếp làm hàm thử trong tích phân. Tác giả sử dụng hàm cắt $B_m(s) = \max(-m, \min(s, m))$ và toán tử Nemytskii tương ứng $\hat{B}_m(w)(x) = B_m(w(x))$ kết hợp với không gian Orlicz để thực hiện phép qua giới hạn suy rộng.
- Phương pháp đơn điệu: Áp dụng tính chất đơn điệu của toán tử tựa tuyến tính để xác định giới hạn phi tuyến.
- Hệ thống các bất đẳng thức giải tích: Áp dụng bất đẳng thức Cauchy với số thực dương $\varepsilon$, bất đẳng thức Young, bất đẳng thức Hölder, bất đẳng thức Poincaré, bất đẳng thức Gronwall dạng vi phân, tích phân và dạng đều.
Nguồn tư liệu và dữ liệu
Luận án thuộc hình thức toán học lý thuyết; toàn bộ dữ liệu là các cấu trúc giải tích, định lý, bổ đề và các chuỗi ước lượng tiên nghiệm được chứng minh trực tiếp từ các tiên đề toán học và các công trình khoa học trong danh mục tài liệu tham khảo.
Nội dung chính theo từng chương
Chương 1: Kiến thức chuẩn bị
Chương 1 hệ thống hóa các định nghĩa, ký hiệu và kết quả giải tích cơ sở:
- Trình bày định nghĩa các lớp không gian Sobolev có trọng $M_0^{1,p}(\Omega, a)$ trên miền bị chặn và $M_0^{1,p}(\mathbb{R}^N, a)$ trên toàn không gian.
- Định nghĩa các toán tử suy biến: $$L_{p,a}u := -\text{div}(a(x)|\nabla u|^{p-2}\nabla u)$$ và toán tử trên toàn không gian $\mathbb{R}^N$: $$\mathcal{L}_{p,a}u := -\text{div}(a(x)|\nabla u|^{p-2}\nabla u) + \lambda |u|^{p-2}u$$
- Trình bày lý thuyết nửa nhóm liên tục, định nghĩa tập hút toàn cục, khái niệm tính compắc tiệm cận, cấu trúc tập $\omega$-giới hạn và các điều kiện đủ để tồn tại tập hút toàn cục.
- Thiết lập khung lý thuyết tập hút lùi cho quá trình hai tham số $U(t,a)$, khái niệm tính liên tục mạnh - yếu (norm-to-weak continuity), họ tập hấp thụ lùi và mối liên hệ giữa tính compắc tiệm cận lùi với độ đo Kuratowski $\alpha$.
- Tập hợp các bất đẳng thức giải tích: Cauchy, Young, Hölder, Gronwall vi phân, Gronwall tích phân, Gronwall đều và Poincaré.
- Nêu các bổ đề giải tích hàm quan trọng: Bổ đề compắc Aubin-Lions-Simon (Định lý 1.5) về phép nhúng compact của không gian $E_{p,r} = {u \in L^p(0,T; B_0), u_t \in L^r(0,T; B_2)}$ vào $L^p(0,T; B_1)$; Bổ đề 1.3 về tính hội tụ yếu của dãy phi tuyến trong $L^p(M)$; Bổ đề 1.4 về tính chất hội tụ tích phân dưới tác động của hàm thử $\xi \in C_0([0,T]; M_0^{1,p} \cap L^\infty)$.
Chương 2: Tập hút toàn cục đối với một lớp phương trình parabolic suy biến tựa tuyến tính với số hạng phi tuyến bậc tùy ý
Chương này dẫn giải dựa theo công trình khoa học số [1] của tác giả, nghiên cứu bài toán: $$u_t - \text{div}(a(x)|\nabla u|^{p-2}\nabla u) + f(u) = g(x), \quad x \in \Omega, ; t > 0$$ với $u(x,t) = 0$ trên $\partial\Omega$ và $u(x,0) = u_0(x) \in L^2(\Omega)$, $2 \le p \le N$.
Các giả thiết gồm:
- (H1): $a \in L_{loc}^1(\Omega)$, $a(x) = 0$ trên $\Sigma$ ($m(\Sigma) = 0$), $a(x) > 0$ trên $\Omega \setminus \Sigma$ và $\int_\Omega [a(x)]^{-N/\alpha} dx < \infty$ với $\alpha \in (0, p)$.
- (H2): $f \in C^1(\mathbb{R})$, $f'(s) \ge -A$ và $f(s)s \ge \mu |s|^p - c_1$ với $c_1, A, \mu > 0$ (khi $p=2$ đòi hỏi $0 < \mu < c_0$).
- (H3): $g \in L^2(\Omega)$.
Các phát hiện và kết quả giải tích của Chương 2:
- Mệnh đề nhúng 2.1:
- Phép nhúng $M_0^{1,p}(\Omega, a) \hookrightarrow W_0^{1,\beta}(\Omega)$ liên tục nếu $1 \le \beta \le \frac{pN}{N+\alpha}$.
- Phép nhúng $M_0^{1,p}(\Omega, a) \hookrightarrow L^r(\Omega)$ liên tục nếu $1 \le r \le p^*_\alpha := \frac{pN}{N-p+\alpha}$.
- Phép nhúng $M_0^{1,p}(\Omega, a) \hookrightarrow L^r(\Omega)$ compắc nếu $1 \le r < p^*\alpha$. Từ đó suy ra tồn tại hằng số $c_0 > 0$ thỏa mãn $c_0|u|{L^2(\Omega)}^2 \le |u|_{M_0^{1,p}(\Omega, a)}^p$.
- Mệnh đề về tính compắc của toán tử: Phép nhúng $D(L_{p,a}) \hookrightarrow M_0^{1,p}(\Omega, a)$ là compắc.
- Sự tồn tại và duy nhất nghiệm yếu toàn cục: Sử dụng phương pháp Galerkin và hàm cắt Nemytskii $\hat{B}m(w)(x) = B_m(w(x))$ để chứng minh sự tồn tại nghiệm yếu $u \in C([0,T]; L^2(\Omega)) \cap L^p(0,T; M_0^{1,p}(\Omega, a))$ với $u_t \in L^{p'}(0,T; M^{-1,p'}(\Omega, a)) + L^1(\Omega_T)$. Đánh giá chênh lệch nghiệm $w = u_1 - u_2$ dẫn tới bất đẳng thức: $$|w(t)|{L^2(\Omega)}^2 \le |w(0)|{L^2(\Omega)}^2 + 2A \int_0^t |w(s)|{L^2(\Omega)}^2 ds$$ khẳng định tính duy nhất và sự phụ thuộc liên tục vào điều kiện ban đầu.
- Tập hút toàn cục:
- Thiết lập sự tồn tại tập hấp thụ trong $L^2(\Omega)$: Tồn tại bán kính $\rho_0 > 0$ sao cho với mọi tập bị chặn $B \subset L^2(\Omega)$, tồn tại $T_1(B)$ để $|S(t)u_0|_{L^2(\Omega)}^2 \le \rho_0$ với mọi $t \ge T_1$.
- Định lý 2.1: Nửa nhóm ${S(t)}{t \ge 0}$ sở hữu một tập hút toàn cục $\mathcal{A}{L^2}$ trong $L^2(\Omega)$.
- Bằng cách nhân phương trình với $u_t$ và tích phân, tác giả thu được: $$\frac{d}{dt}\left[\frac{1}{p}\int_\Omega a(x)|\nabla u|^p dx + \int_\Omega F(u)dx - \int_\Omega gu dx\right] = -|u_t|{L^2(\Omega)}^2 \le 0$$ Từ đó xây dựng tập hấp thụ bị chặn trong miền xác định $D(L{p,a})$.
- Định lý 2.2: Nửa nhóm ${S(t)}{t \ge 0}$ có một tập hút toàn cục $\mathcal{A}{M_0^{1,p}}$ trong không gian Sobolev có trọng $M_0^{1,p}(\Omega, a)$.
Chương 3: Tập hút lùi đối với lớp phương trình parabolic suy biến tựa tuyến tính không ôtônôm
Chương này dẫn giải theo công trình khoa học số [2] của tác giả, khảo sát bài toán trên toàn miền $\mathbb{R}^N$ ($N \ge 3$): $$\frac{\partial u}{\partial t} - \text{div}(a(x)|\nabla u|^{p-2}\nabla u) + \lambda |u|^{p-2}u + f(u) = g(x,t), \quad x \in \mathbb{R}^N, ; a < t$$ với $u(x,a) = u_a(x) \in L^2(\mathbb{R}^N)$, $2 \le p < N$, $\lambda > 0$.
Các điều kiện giải tích bao gồm:
- Điều kiện phi tuyến (F): Hàm $f \in C^1(\mathbb{R})$ thỏa mãn: $$\alpha_1 |u|^q - \beta_1 |u|^p + \gamma_1 u^2 \le f(u)u \le \alpha_2 |u|^q + \beta_2 |u|^p + \gamma_2 u^2, \quad q \ge 2, ; \beta_1 < \lambda$$
- Điều kiện mới $(H_a)$ về hệ số khuếch tán: Hàm đo được không âm $a \in L_{loc}^1(\mathbb{R}^N)$, tồn tại $\alpha \in (0, p)$ sao cho $\liminf_{x \to z} \frac{a(x)}{|x-z|^\alpha} > 0$, đồng thời thỏa mãn một trong hai điều kiện:
- Tồn tại $M_0 > 0$ sao cho $\sup_{m \ge M_0} \sup_{m \le |x| \le 2m} a(x) < \infty$; hoặc
- Tồn tại $M_0 > 0$ sao cho $\sup_{m \ge M_0} \int_{m \le |x| \le 2m} [a(x)]^{-\frac{2p-2}{p-2}} dx < \infty$.
- Điều kiện ngoại lực $(l)$: Hàm $g(x,t)$ thỏa mãn điều kiện khả tích lùi thích hợp.
Các phát hiện và kết quả giải tích của Chương 3:
- Định lý 3.1: Bài toán có nghiệm yếu duy nhất $u(t) = U(t,a)u_a$, thỏa mãn ước lượng tiêu hao: $$|u(t)|{L^2(\mathbb{R}^N)}^2 \le C e^{-\gamma_1(t-a)}|u_a|{L^2(\mathbb{R}^N)}^2 + e^{-\gamma_1 t}\int_a^t e^{\gamma_1 s}|g(s)|_{L^2(\mathbb{R}^N)}^2 ds$$
- Bổ đề 3.1: Quá trình $U(t,a)$ có một họ tập hấp thụ lùi bị chặn $\hat{\mathcal{B}} = {B(t)}{t \in \mathbb{R}}$ trong không gian giao $M_0^{1,p}(\mathbb{R}^N, a) \cap L^q(\mathbb{R}^N) \cap L^2(\mathbb{R}^N)$, nghĩa là tồn tại $r_0(t) > 0$ và $a_0(t, B) < t$ sao cho: $$|u(t)|{M_0^{1,p}(\mathbb{R}^N, a)}^p + |u(t)|{L^q(\mathbb{R}^N)}^q + |u(t)|{L^2(\mathbb{R}^N)}^2 \le r_0(t), \quad \forall a \le a_0$$
- Thiết lập tập hút lùi: Do xét trên toàn không gian $\mathbb{R}^N$, phép nhúng Sobolev không còn tính compắc. Tác giả kết hợp kỹ thuật phân rã miền không gian, đánh giá đuôi tích phân tại vô cùng và sử dụng độ đo Kuratowski để chứng minh tính compắc tiệm cận lùi.
- Luận án chứng minh sự tồn tại của tập hút lùi duy nhất $\hat{\mathcal{A}}$ lần lượt trong các không gian:
- Trong không gian $L^2(\mathbb{R}^N)$;
- Trong không gian $L^q(\mathbb{R}^N)$;
- Trong không gian Sobolev có trọng $M_0^{1,p}(\mathbb{R}^N, a)$.
Chương 4: Sự tồn tại và dáng điệu tiệm cận nghiệm đối với một lớp hệ phương trình phản ứng khuếch tán với số hạng phi tuyến kiểu mũ
Chương này nghiên cứu hệ phản ứng khuếch tán gồm hai phương trình liên kết: $$\begin{cases} \frac{\partial u}{\partial t} - \Delta u + f(u,v) = g_1(x), & x \in \Omega, ; t > 0 \ \frac{\partial v}{\partial t} - \Delta v + h(u,v) = g_2(x), & x \in \Omega, ; t > 0 \end{cases}$$ với $u=v=0$ trên $\partial\Omega$, trong đó số hạng phi tuyến $f, h$ mang tính chất tăng trưởng kiểu mũ (tiêu biểu như trường hợp chứa thành phần $e^u$).
Các phát hiện và kết quả giải tích của Chương 4:
- Chứng minh tính đặt đúng toàn cục: Thiết lập sự tồn tại và tính duy nhất của nghiệm yếu $(u, v)$ trong không gian thích hợp thông qua phương pháp xấp xỉ và kỹ thuật đánh giá năng lượng.
- Chứng minh sự tồn tại của tập hút toàn cục đối với hệ động lực tương ứng trong không gian tích $(L^2(\Omega))^2$.
- Phân tích cấu trúc tập hút toàn cục trong trường hợp hệ gradient: chỉ ra tập hút toàn cục trùng với bao đóng của tập các đa tạp không bền xuất phát từ các điểm dừng của hệ.
- Nghiệm dừng và tính ổn định: Chứng minh định lý về sự tồn tại của nghiệm dừng (stationary solutions) đối với hệ phản ứng khuếch tán, đồng thời thiết lập các điều kiện đủ để nghiệm dừng đạt tính ổn định mũ (exponential stability).
Bảng tổng hợp cấu trúc các bài toán và kết quả giải tích chính trong luận án
| Đặc điểm so sánh | Chương 2 | Chương 3 | Chương 4 |
|---|---|---|---|
| Dạng bài toán | Phương trình parabolic suy biến tự trị | Phương trình parabolic suy biến không ôtônôm | Hệ phương trình phản ứng khuếch tán tự trị |
| Miền không gian | Miền bị chặn $\Omega \subset \mathbb{R}^N$ ($N \ge 2$) | Toàn không gian $\mathbb{R}^N$ ($N \ge 3$) | Miền bị chặn $\Omega \subset \mathbb{R}^N$ |
| Toán tử không gian | $L_{p,a}u = -\text{div}(a(x)|\nabla u|^{p-2}\nabla u)$ | $\mathcal{L}_{p,a}u = -\text{div}(a|\nabla u|^{p-2}\nabla u) + \lambda |u|^{p-2}u$ | Hệ Laplace liên kết $(-\Delta u, -\Delta v)$ |
| Hệ số khuếch tán | $a(x) \ge 0$, triệt tiêu trên $\Sigma$, $\int_\Omega a^{-N/\alpha} dx < \infty$ | $a(x) \ge 0$, thỏa mãn điều kiện tích phân địa phương $(H_a)$ | $a(x) \equiv 1$ (Laplace thông thường) |
| Dạng số hạng phi tuyến | $f(u)$ tăng trưởng bậc tùy ý | $f(u)$ tăng trưởng cấp đa thức bậc $q \ge 2$ | $f(u,v), h(u,v)$ tăng trưởng kiểu mũ |
| Khung động lực học | Nửa nhóm liên tục ${S(t)}_{t \ge 0}$ | Quá trình hai tham số ${U(t,a)}_{t \ge a}$ | Nửa nhóm liên tục trên không gian tích |
| Không gian tồn tại tập hút | $L^2(\Omega)$ và $M_0^{1,p}(\Omega, a)$ | $L^2(\mathbb{R}^N)$, $L^q(\mathbb{R}^N)$, $M_0^{1,p}(\mathbb{R}^N, a)$ | Không gian tích $(L^2(\Omega))^2$ |
| Kết quả bổ sung | Xử lý triệt để hàm thử bằng cắt Nemytskii | Khắc phục thiếu compact trên $\mathbb{R}^N$ không dùng $(H_{\alpha,\beta})$ | Cấu trúc hệ gradient; sự tồn tại và ổn định mũ nghiệm dừng |
Kết quả và những đóng góp mới
Luận án tự nêu các đóng góp mới về mặt lý luận giải tích toán học như sau:
-
Về phương trình parabolic suy biến với phi tuyến tăng trưởng tùy ý: Thiết lập sự tồn tại và duy nhất nghiệm yếu toàn cục, đồng thời chứng minh sự tồn tại của tập hút toàn cục trong các không gian $L^2(\Omega)$ và không gian Sobolev có trọng $M_0^{1,p}(\Omega, a)$ cho lớp phương trình parabolic suy biến tựa tuyến tính chứa toán tử $p$-Laplace có trọng. Điểm mới căn bản là loại bỏ hoàn toàn các giả thiết hạn chế cận trên của bậc tăng trưởng đa thức đối với hàm phi tuyến $f(u)$ vốn tồn tại trong các công trình tiền đề [27, 41].
-
Về phương trình parabolic suy biến không ôtônôm trên toàn không gian $\mathbb{R}^N$: Thiết lập giả thiết mới $(H_a)$ cho hệ số khuếch tán $a(x)$ dựa trên tính khả tích địa phương bậc cao của hàm $a(x)^{-1}$ trên các hình vành khăn thay thế cho giả thiết $(H_{\alpha,\beta})$ về tốc độ tăng trưởng ra vô cùng. Trên cơ sở đó, luận án chứng minh thành công sự tồn tại của tập hút lùi duy nhất trong $L^2(\mathbb{R}^N)$, $L^q(\mathbb{R}^N)$ và $M_0^{1,p}(\mathbb{R}^N, a)$, khắc phục trở ngại mất tính compắc của phép nhúng Sobolev trên miền không bị chặn.
-
Về hệ phương trình phản ứng khuếch tán với phi tuyến kiểu mũ: Mở rộng các kết quả nghiên cứu tập hút toàn cục đối với hệ phản ứng khuếch tán từ lớp phi tuyến đa thức/Sobolev sang lớp hàm phi tuyến có tốc độ tăng trưởng kiểu mũ. Xác định cấu trúc tập hút toàn cục trong trường hợp hệ gradient, đồng thời chứng minh sự tồn tại của nghiệm dừng và thiết lập điều kiện đủ để nghiệm dừng đạt tính ổn định mũ.
Kiến nghị và định hướng ứng dụng
Về mặt lý thuyết, các kết quả của luận án cung cấp các công cụ phân tích giải tích để tiếp tục xử lý các bài toán biên đối với phương trình đạo hàm riêng phi tuyến chứa các toán tử suy biến hoặc các số hạng kích động phức tạp. Về mặt thực tiễn, khung giải tích tiệm cận này đóng góp cơ sở toán học để tính toán và dự báo các trạng thái cân bằng trong các mô hình vật lý môi trường không đồng nhất, mô hình khuếch tán nhiệt suy biến và các quá trình sinh thái học mô tả tương tác hai chủng loài.
Hạn chế và hướng nghiên cứu tiếp
Dựa trên các nhận xét nội tại và bình luận cuối các chương của tác giả:
-
Hạn chế: Trong Chương 2, kết quả về tập hút toàn cục đối với phương trình suy biến tựa tuyến tính chứa toán tử $p$-Laplace với phần phi tuyến tăng trưởng tùy ý mới chỉ đạt được trên miền bị chặn $\Omega$. Khi miền $\Omega$ không bị chặn, việc mất tính compắc của phép nhúng Sobolev có trọng kết hợp với sự hiện diện của phần phi tuyến tăng trưởng tùy ý tạo ra các trở ngại kỹ thuật mà phương pháp hiện tại của chương chưa bao quát được.
-
Hướng nghiên cứu tiếp:
- Nghiên cứu mở rộng bài toán của Chương 2 cho trường hợp miền $\Omega$ không bị chặn trong $\mathbb{R}^N$.
- Mở rộng các kết quả tiệm cận của hệ phương trình phản ứng khuếch tán với phần phi tuyến kiểu mũ (Chương 4) sang trường hợp hệ chứa toán tử $p$-Laplace suy biến hoặc xét trên các miền không bị chặ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: Thuộc các chuyên ngành Phương trình vi phân và tích phân, Giải tích toán học, Toán ứng dụng và Lý thuyết hệ động lực vô hạn chiều.
- Giảng viên và chuyên gia nghiên cứu PDEs: Sử dụng làm tài liệu tham khảo phục vụ giảng dạy các học phần chuyên đề giải tích hàm phi tuyến, lý thuyết phương trình đạo hàm riêng và các phương pháp định tính trong hệ động lực.
Các phần có giá trị tham khảo kỹ thuật cao:
- Kỹ thuật cắt Nemytskii trong Chương 2: Cung cấp phương pháp chi tiết xử lý phép qua giới hạn trong không gian Orlicz khi phương trình phi tuyến không thỏa mãn bất đẳng thức tăng trưởng đa thức.
- Kỹ thuật đánh giá địa phương và độ đo Kuratowski trong Chương 3: Cung cấp quy trình toán học để khắc phục sự thiếu vắng phép nhúng compắc Sobolev trên miền không bị chặn $\mathbb{R}^N$ thông qua giả thiết tích phân địa phương $(H_a)$.
Câu hỏi thường gặp
1. Điểm khác biệt mấu chốt giữa điều kiện $(H_a)$ trong Chương 3 và điều kiện $(H_{\alpha,\beta})$ của các nghiên cứu trước đây là gì?
Trong các nghiên cứu trước (như [5, 10]), điều kiện $(H_{\alpha,\beta})$ bắt buộc hệ số khuếch tán $a(x)$ phải tăng nhanh ra vô cùng với tốc độ lớn hơn $|x|^\beta$ ($\beta > p + \frac{N}{2}(p-2)$) khi $|x| \to \infty$ nhằm ép phép nhúng Sobolev có trọng thành compắc. Điều kiện $(H_a)$ do tác giả đề xuất thay thế toàn bộ yêu cầu về dáng điệu tại vô cùng bằng điều kiện khả tích địa phương bậc cao $\sup_{m \ge M_0} \int_{m \le |x| \le 2m} [a(x)]^{-\frac{2p-2}{p-2}} dx < \infty$ trên các miền vành khăn, cho phép áp dụng bài toán cho các lớp hàm trọng không tăng ra vô cùng.
2. Tại sao tác giả phải dùng ánh xạ cắt Nemytskii $\hat{B}_m$ trong bước chứng minh sự tồn tại nghiệm ở Chương 2?
Do phần phi tuyến $f(u)$ có tốc độ tăng trưởng tùy ý (chỉ giả thiết $f'(s) \ge -A$ và $f(s)s \ge \mu |s|^p - c_1$), biểu thức $f(u)$ không nhất thiết thuộc không gian đối ngẫu $M^{-1,p'}(\Omega, a)$ mà chỉ thuộc $L^1(\Omega_T)$. Do đó, hàm nghiệm $u$ không thể dùng trực tiếp làm hàm thử trong phương trình tích phân. Ánh xạ cắt Nemytskii $\hat{B}_m(w)(x) = B_m(w(x))$ với $B_m(s) = \max(-m, \min(s, m))$ làm trơn và chặn hàm thử, cho phép thực hiện việc qua giới hạn khi $m \to \infty$ và $s \to 0$ mà không vi phạm tính khả tích.
3. Trong trường hợp nào luận án sử dụng khái niệm tập hút lùi thay vì tập hút toàn cục?
Khái niệm tập hút lùi được sử dụng trong Chương 3 khi bài toán xét trong trường hợp không ôtônôm, nghĩa là ngoại lực $g(x,t)$ phụ thuộc trực tiếp vào thời gian $t$. Khi đó toán tử nghiệm không còn lập thành nửa nhóm một tham số ${S(t)}{t \ge 0}$ có tính bất biến theo phép tịnh tiến thời gian, mà tạo thành một quá trình hai tham số $U(t,a)$. Do đó, tập hút toàn cục cổ điển không còn áp dụng được, buộc phải chuyển sang nghiên cứu tập hút lùi $\hat{\mathcal{A}} = {A(t)}{t \in \mathbb{R}}$ khi thời điểm bắt đầu lùi xa về âm vô cùng ($a \to -\infty$).
4. Luận án giải quyết trở ngại toán học nào khi mở rộng sang số hạng phi tuyến kiểu mũ $f(u) \sim e^u$ ở Chương 4?
Các nghiên cứu trước đây về hệ phản ứng khuếch tán thường yêu cầu các số hạng phi tuyến phải bị chặn trên bởi đa thức kiểu Sobolev $|f(u)| \le C(1 + |u|^{p-1})$ để đảm bảo toán tử phi tuyến ánh xạ liên tục từ không gian năng lượng vào không gian đối ngẫu. Hàm kiểu mũ $e^u$ tăng nhanh hơn mọi đa thức, phá vỡ cấu trúc đánh giá Sobolev tiêu chuẩn. Luận án đã xây dựng lại các đánh giá tiên nghiệm thích hợp trong không gian tích $(L^2(\Omega))^2$ để đảm bảo sự tồn tại của nghiệm toàn cục và tập hấp thụ.
5. Kết quả về nghiệm dừng trong Chương 4 bao gồm những nội dung cụ thể nào?
Trong Chương 4, tác giả không chỉ chứng minh định lý tồn tại nghiệm dừng đối với hệ phương trình phản ứng khuếch tán phi tuyến kiểu mũ mà còn thiết lập các điều kiện giải tích đủ để nghiệm dừng này có tính ổn định mũ (exponential stability), nghĩa là mọi nghiệm xuất phát gần nghiệm dừng sẽ hội tụ về nghiệm dừng với tốc độ phân rã theo hàm mũ khi thời gian tiến ra vô cùng.
Kết luận
Luận án tiến sĩ của tác giả Trần Thị Quỳnh Chi đã giải quyết một cách có hệ thống các bài toán về tính đặt đúng và dáng điệu tiệm cận của nghiệm đối với ba lớp phương trình và hệ phương trình parabolic phi tuyến. Bằng việc kết hợp các phương pháp giải tích hàm hiện đại như xấp xỉ Galerkin, ánh xạ cắt Nemytskii, bổ đề compắc Aubin-Lions-Simon và độ đo không compắc Kuratowski, công trình đã thiết lập các kết quả giải tích mới về sự tồn tại của tập hút toàn cục với phi tuyến tăng trưởng tùy ý, tập hút lùi dưới điều kiện trọng mới $(H_a)$ trên toàn không gian $\mathbb{R}^N$, và tập hút cùng tính ổn định của hệ phản ứng khuếch tán phi tuyến kiểu mũ. Đây là một công trình nghiên cứu lý thuyết hoàn chỉnh, đóng góp kết quả trực tiếp cho chuyên ngành Phương trình vi phân và tích phân.