Tổng quan luận án
Hệ phương trình Navier-Stokes-Voigt (đôi khi được viết là Navier-Stokes-Voight) được A. P. Oskolkov giới thiệu lần đầu tiên vào năm 1973 như một mô hình toán học mô tả chuyển động của chất lưu nhớt đàn hồi tuyến tính không nén được. Về sau, hệ phương trình này được Y. Cao, E. M. Lunasin và E. S. Titi (2006) đề xuất như một mô hình chính quy hóa cho hệ phương trình Navier-Stokes ba chiều phục vụ mô phỏng số trực tiếp. Hệ Navier-Stokes-Voigt thuộc nhóm các mô hình $\alpha$ trong cơ học chất lưu với ưu điểm không đòi hỏi các điều kiện biên nhân tạo bổ sung ngoài điều kiện biên Dirichlet để đảm bảo tính đặt đúng toàn cục. Bên cạnh đó, phương trình còn được ứng dụng trong kỹ thuật xử lý và phục hồi hình ảnh (image inpainting).
Dù các kết quả về tính tồn tại duy nhất nghiệm, nghiệm tiệm cận và dáng điệu dài hạn (tập hút) của hệ Navier-Stokes-Voigt trong miền bị chặn, miền không bị chặn và không gian $\mathbb{R}^3$ đã được nhiều tác giả nghiên cứu, lý thuyết điều khiển tối ưu đối với hệ phương trình Navier-Stokes-Voigt ba chiều phi dừng vẫn chưa được khảo sát trước thời điểm luận án thực hiện. Đây là khoảng trống nghiên cứu chính mà tác giả Trần Minh Nguyệt lựa chọn giải quyết trong luận án tiến sĩ toán học chuyên ngành Phương trình vi phân và tích phân (mã số 9.03), dưới sự hướng dẫn khoa học của GS. Cung Thế Anh tại Trường Đại học Sư phạm Hà Nội (2019).
Mục tiêu nghiên cứu của luận án được xác định cụ thể thông qua việc chứng minh tính tồn tại nghiệm tối ưu và thiết lập hệ điều kiện tối ưu cần bậc một cùng điều kiện tối ưu đủ bậc hai cho ba bài toán:
- Bài toán (P1) - Điều khiển phân bố tối ưu (A distributed optimal control problem): Khảo sát hệ phương trình Navier-Stokes-Voigt phi dừng ba chiều với phiếm hàm mục tiêu dạng toàn phương, trong đó biến điều khiển đóng vai trò là lực tác động thuộc một tập lồi, đóng, khác rỗng bất kỳ.
- Bài toán (P2) - Điều khiển thời gian tối ưu (A time optimal control problem): Khảo sát bài toán điều khiển thời gian cho hệ Navier-Stokes-Voigt phi dừng ba chiều với tập điều khiển chấp nhận được là tập lồi, đóng, khác rỗng bất kỳ.
- Bài toán (P3) - Điều khiển biên tối ưu (An optimal boundary control problem): Khảo sát bài toán điều khiển biên phi tuyến cho hệ Navier-Stokes-Voigt ba chiều với phiếm hàm mục tiêu dạng toàn phương và biến điều khiển biên thỏa mãn các điều kiện tương thích.
Đối tượng nghiên cứu là hệ phương trình vi phân đạo hàm riêng Navier-Stokes-Voigt phi dừng ba chiều trong miền không gian mở, bị chặn $\Omega \subset \mathbb{R}^3$ với biên $\Gamma$ Lipschitz cục bộ hoặc trơn cấp $C^m$, trong khoảng thời gian hữu hạn $(0, T)$ với hình trụ không - thời gian $Q = \Omega \times (0, T)$. Luận án tập trung chủ yếu vào không gian ba chiều (trường hợp mang ý nghĩa vật lý thực tế nhất), đồng thời các kết quả và kỹ thuật chứng minh vẫn giữ nguyên tính đúng đắn trong không gian hai chiều.
Tổng quan tài liệu và vị trí của luận án
Lý thuyết điều khiển tối ưu cho các phương trình vi phân đạo hàm riêng trong cơ học chất lưu bắt đầu phát triển mạnh từ những năm 1980 với các công trình của A. V. Fursikov (1980, 1981) về sự tồn tại nghiệm điều khiển đối với hệ phương trình Navier-Stokes. Năm 1990, F. Abergel và R. Temam công bố công trình nghiên cứu về các điều kiện tối ưu đầu tiên cho dòng chảy chất lưu, mở ra hướng nghiên cứu phân tích tối ưu hóa cho hệ phương trình Navier-Stokes.
Các hướng nghiên cứu tiền đề được luận án điểm lại và kế thừa gồm:
- Nghiên cứu về phương trình Navier-Stokes-Voigt: A. P. Oskolkov (1973) thiết lập hệ phương trình mô tả chất lỏng Kelvin-Voigt; Y. Cao, E. M. Lunasin, E. S. Titi (2006) khảo sát tính chính quy hóa và giải thuật số; C. Bertozzi, A. L. Bertozzi (2001), B. Garcia-Archilla, J. Novo, E. S. Titi (1999) phân tích các mô hình $\alpha$; V. K. Kalantarov, B. Levine, E. S. Titi cùng nhiều tác giả khác phân tích tính giải được và tập hút; các nghiên cứu trên toàn không gian $\mathbb{R}^3$ về nghiệm và tốc độ suy giảm nghiệm của các tác giả [4, 56, 75].
- Điều khiển phân bố cho hệ Navier-Stokes: Các công trình của H. O. Fattorini, S. S. Sritharan (1994), M. Hinze (2000), D. Wachsmuth (2006, 2007) cho trường hợp không có ràng buộc trạng thái; G. Wang (1998) và C. Liu (2001) cho trường hợp có ràng buộc trạng thái.
- Điều khiển thời gian cho hệ Navier-Stokes: Các nghiên cứu của V. Barbu (1998) và E. Fernández-Cara (1999).
- Điều khiển biên cho hệ Navier-Stokes: Khảo sát trường hợp dừng bởi M. D. Gunzburger, L. S. Hou, T. P. Svobodny (1991, 1992), J. P. Raymond (2006); trường hợp phi dừng bởi E. Casas (1997), A. V. Fursikov, O. Y. Imanuvilov (1998), J. P. Raymond (2007), cùng nguyên lý cực đại Pontryagin của B. Rousselet và các tổng quan trong luận án tiến sĩ, chuyên khảo của D. Wachsmuth (2006), K. Kunisch và cộng sự.
Vị trí của luận án: Trong khi các bài toán điều khiển tối ưu cho phương trình Navier-Stokes kinh điển đã được khảo sát tương đối phong phú, bài toán điều khiển tối ưu cho hệ Navier-Stokes-Voigt ba chiều chưa từng được công bố trước đó. Do tính chất của hệ Navier-Stokes-Voigt có chứa số hạng đạo hàm bậc cao theo thời gian và không gian $-\alpha^2 \Delta y_t$, cấu trúc giải tích của bài toán đòi hỏi việc xây dựng không gian hàm phù hợp, phân tích các toán tử phi tuyến và phương trình đối ngẫu tương ứng. Luận án giải quyết trọn vẹn bài toán tồn tại nghiệm, thiết lập phương trình liên hợp, đưa ra bất đẳng thức biến phân bậc một và thiết lập điều kiện đủ bậc hai thông qua tính tăng bậc hai (quadratical growth).
Cơ sở lý thuyết và phương pháp nghiên cứu
Luận án sử dụng khung lý thuyết giải tích hàm hiện đại, lý thuyết phương trình đạo hàm riêng phi tuyến và giải tích lồi:
- Không gian hàm và không gian hàm trừu tượng: Không gian Lebesgue $L^p(\Omega)$, $L^2_0(\Omega)$, không gian Sobolev $W^{m,p}(\Omega)$, $H^m(\Omega)$, $H_0^1(\Omega)$, không gian vết trên biên $H^s(\Gamma)$, $H^{1/2}(\Gamma)$. Không gian hàm solenoidal $\mathcal{V} = {y \in (C_0^\infty(\Omega))^3 : \text{div } y = 0}$, cùng các bao đóng $H$ (trong $(L^2(\Omega))^3$) và $V$ (trong $(H_0^1(\Omega))^3$). Các không gian hàm nhận giá trị trong không gian Banach/Hilbert: $L^p(0, T; X)$, $C([0, T]; X)$ và $W^{1,p}(0, T; X) = {y \in L^p(0, T; X) : y_t \in L^p(0, T; X)}$.
- Toán tử và dạng song tuyến/tam tuyến: Dạng tam tuyến $b(u, v, w) = \sum_{i,j=1}^3 \int_\Omega u_i \frac{\partial v_j}{\partial x_i} w_j dx$, toán tử tuyến tính liên tục $A: V \to V'$ xác định bởi $\langle Au, v \rangle_{V',V} = (\nabla u, \nabla v)$, và toán tử phi tuyến $B: V \times V \to V'$ xác định bởi $\langle B(u, v), w \rangle_{V',V} = b(u, v, w)$.
- Bất đẳng thức giải tích: Bất đẳng thức Hölder tổng quát, bất đẳng thức Poincaré, bất đẳng thức Ladyzhenskaya ba chiều $|v|_{L^4(\Omega)} \le 2^{1/2} |v|^{1/4} |\nabla v|^{3/4}$, bổ đề Gronwall vi phân và tích phân, bất đẳng thức Young chứa tham số $\epsilon$.
- Giải tích lồi: Khái niệm nón chuẩn tắc $N_U(u)$, nón tiếp tuyến cực $T_U(u)$, nón các hướng chấp nhận được $F_U(u)$ và định lý đồng nhất $F_U(u) = T_U(u)$ trên tập lồi.
- Định lý nhúng compact: Định lý Rellich-Kondrachov, các định lý nhúng compact kiểu Aubin-Lions và Simon: $W^{1,2}(0, T; V) \hookrightarrow L^2(Q)$, $W^{1,2}(0, T; L^2(\Omega)) \hookrightarrow C([0, T]; V')$, $W^{1,2}(0, T; H^{1/2}(\Gamma)) \hookrightarrow C([0, T]; L^2(\Gamma))$.
Phương pháp nghiên cứu trong luận án là phương pháp toán học thuần túy (định tính và giải tích phương trình vi phân):
- Chứng minh tồn tại nghiệm tối ưu: Sử dụng phương pháp dãy cực tiểu hóa kết hợp tính đóng yếu, tính compact của các phép nhúng không gian hàm trừu tượng để chuyển qua giới hạn yếu trong các số hạng phi tuyến.
- Thiết lập điều kiện tối ưu cần bậc một: Thay vì dùng phương pháp nhân tử Lagrange hình thức, tác giả chọn hướng tiếp cận tính trực tiếp đạo hàm theo hướng của phiếm hàm mục tiêu qua ánh xạ điều khiển - trạng thái (control-to-state mapping $S: u \mapsto y$), kết hợp giải hệ phương trình tuyến tính hóa và hệ phương trình trạng thái liên hợp (trực đối ngẫu) để rút ra bất đẳng thức biến phân dạng tường minh.
- Thiết lập điều kiện tối ưu đủ bậc hai: Sử dụng phương pháp phản chứng (contradiction argument) mở rộng kỹ thuật của D. Wachsmuth để phân tích dãy hướng suy biến, chứng minh tính tăng bậc hai cục bộ của phiếm hàm mục tiêu trong lân cận chuẩn $L^2(Q)$.
| Yếu tố phương pháp | Nội dung triển khai trong luận án |
|---|---|
| Mô hình phương trình trạng thái | Hệ phi dừng Navier-Stokes-Voigt ba chiều: $y_t - \nu \Delta y - \alpha^2 \Delta y_t + (y \cdot \nabla)y + \nabla p = u$ |
| Không gian trạng thái | $W^{1,2}(0, T; V)$ liên tục nhúng vào $C([0, T]; V)$ |
| Không gian điều khiển | $L^2(Q)$ đối với điều khiển phân bố; tập lồi đóng $U_{ad}$ |
| Phương pháp giải tích tối ưu | Đạo hàm theo hướng trực tiếp, hệ phương trình liên hợp, phương pháp phản chứng |
Nội dung chính theo từng chương
Chương 1: Kiến thức chuẩn bị và các kết quả bổ trợ (Preliminaries and Auxiliary Results)
Chương 1 tổng hợp các định nghĩa, không gian hàm Sobolev, không gian véc-tơ solenoidal, các bất đẳng thức giải tích cơ bản và tính chất nhúng liên tục, nhúng compact. Trọng tâm của chương nằm ở các kết quả giải tích liên quan trực tiếp đến hệ Navier-Stokes-Voigt:
- Tính giải được của hệ Navier-Stokes-Voigt thuần nhất: Định nghĩa nghiệm yếu của hệ (1.19) trên $(0, T)$ trong không gian $W^{1,2}(0, T; V)$ và nhắc lại định lý tồn tại duy nhất nghiệm yếu $y \in W^{1,2}(0, T; V) \cap C([0, T]; V)$ với mọi $y_0 \in V$ và $f \in L^2(0, T; V')$.
- Tính giải được và tính chất nghiệm của hệ tuyến tính hóa: Khảo sát hệ phương trình tuyến tính hóa (1.20): $$z_t - \nu \Delta z - \alpha^2 \Delta z_t + (z \cdot \nabla)\bar{y} + (\bar{y} \cdot \nabla)z + \nabla q = m$$ Chứng minh hệ có duy nhất nghiệm yếu $z \in W^{1,2}(0, T; V)$, thỏa mãn đánh giá tiên nghiệm: $$|z|{W^{1,2}(0, T; V)} \le C |m|{L^2(0, T; V')}$$ Đồng thời chứng minh tính chính quy cao hơn: nếu $m \in W^{1,2}(0, T; V')$ thì $z \in C^1([0, T]; V)$; nếu dãy $m_k \to m$ trong $C([0, T]; V')$ thì $z_{kt} \to z_t$ trong $C([0, T]; V)$.
- Tính chất liên tục của toán tử phi tuyến: Chứng minh rằng nếu $y_n \rightharpoonup y$ yếu trong $W^{1,2}(0, T; V)$ thì $B(y_n, y_n) \to B(y, y)$ mạnh trong $L^2(0, T; V')$.
Chương 2: Bài toán điều khiển phân bố tối ưu (A Distributed Optimal Control Problem)
Chương 2 giải quyết bài toán (P1) với phiếm hàm mục tiêu: $$J(y, u) = \frac{\alpha_T}{2} \int_\Omega |y(x, T) - y_T(x)|^2 dx + \frac{\alpha_Q}{2} \iint_Q |y(x, t) - y_Q(x, t)|^2 dxdt + \frac{\gamma}{2} \iint_Q |u(x, t)|^2 dxdt$$ trong đó $y_0 \in V$, $y_T \in V$, $y_Q \in L^2(Q)$, $\alpha_T, \alpha_Q \ge 0$ (ít nhất một hệ số dương), $\gamma > 0$ và $U_{ad}$ là tập con lồi, đóng, khác rỗng trong $L^2(Q)$.
- Sự tồn tại nghiệm tối ưu: Định lý 2.1 chứng minh bài toán luôn có nghiệm tối ưu toàn cục $\bar{u} \in U_{ad}$ cùng trạng thái tương ứng $\bar{y} \in W^{1,2}(0, T; V)$ thông qua tính bị chặn của dãy cực tiểu $(y_n, u_n)$ trong $W^{1,2}(0, T; V) \times L^2(Q)$ và tính lồi nửa liên tục dưới của phiếm hàm $J$.
- Điều kiện cần bậc một: Xây dựng hệ phương trình trạng thái liên hợp (2.9): $$\begin{cases} -w_t - \nu \Delta w + \alpha^2 \Delta w_t - (\bar{y} \cdot \nabla)w + (\nabla \bar{y})^T w + \nabla p = \alpha_Q(\bar{y} - y_Q) \ \text{div } w = 0 \ w|{\Gamma} = 0 \ w(T) - \alpha^2 \Delta w(T) + \nabla q = \alpha_T(\bar{y}(T) - y_T) \end{cases}$$ Chứng minh hệ liên hợp có duy nhất nghiệm yếu $w \in W^{1,2}(0, T; V)$. Định lý 2.2 thiết lập điều kiện cần bậc một qua bất đẳng thức biến phân: $$\iint_Q (w + \gamma \bar{u}) \cdot h , dxdt \ge 0, \quad \forall h \in T{U_{ad}}(\bar{u})$$ tương đương $\iint_Q (w + \gamma \bar{u}) \cdot (v - \bar{u}) dxdt \ge 0, \forall v \in U_{ad}$.
- Các trường hợp cụ thể:
- Ràng buộc hộp (Box constraints): $u_{a,i}(x, t) \le u_i(x, t) \le u_{b,i}(x, t)$, nghiệm tối ưu có biểu diễn phép chiếu $\bar{u}i(x, t) = P{[u_{a,i}, u_{b,i}]} \left( -\frac{w_i(x, t)}{\gamma} \right)$.
- Không có ràng buộc ($U_{ad} = L^2(Q)$): $\bar{u}(x, t) = -\frac{1}{\gamma} w(x, t)$.
- Ràng buộc hình cầu $|u|_{L^2(Q)} \le R$: $\bar{u}(x, t) = -\frac{1}{\gamma + s} w(x, t)$ với $s \ge 0$.
- Điều kiện đủ bậc hai: Định lý 2.3 chỉ ra rằng nếu cặp chấp nhận được $\bar{v} = (\bar{y}, \bar{u})$ thỏa mãn điều kiện cần bậc một cùng điều kiện đủ bậc hai: $$q(z, h) = \frac{\alpha_T}{2} \int_\Omega |z(T)|^2 dx + \frac{\alpha_Q}{2} \iint_Q |z|^2 dxdt + \frac{\gamma}{2} \iint_Q |h|^2 dxdt - \int_0^T b(z(t), z(t), w(t)) dt > 0$$ với mọi $h \in (T_{U_{ad}}(\bar{u}) \cap C(\bar{u})) \setminus {0}$, thì tồn tại $\varepsilon > 0, \rho > 0$ sao cho: $$J(y, u) \ge J(\bar{y}, \bar{u}) + \varepsilon |u - \bar{u}|^2_{L^2(Q)}$$ với mọi $u \in U_{ad}$ thỏa mãn $|u - \bar{u}|_{L^2(Q)} \le \rho$.
Chương 3: Bài toán điều khiển thời gian tối ưu (A Time Optimal Control Problem)
Chương 3 trình bày các kết quả về bài toán điều khiển hệ Navier-Stokes-Voigt từ trạng thái ban đầu $y_0 \in V$ đạt tới trạng thái mục tiêu xác định trong thời gian $T$ ngắn nhất:
- Phát biểu bài toán: Thiết lập bài toán cực tiểu hóa thời gian $T > 0$ sao cho tồn tại hàm điều khiển $u \in U_{ad}$ đưa quỹ đạo nghiệm của hệ từ $y_0$ đến tập mục tiêu thỏa mãn các ràng buộc hình học và ràng buộc trạng thái.
- Sự tồn tại của nghiệm tối ưu thời gian: Chứng minh tính tồn tại của thời gian tối ưu $T^*$ và điều khiển tối ưu $\bar{u}$ tương ứng trên cơ sở tính liên tục của hàm quỹ đạo trạng thái theo biến thời gian và tính compact yếu của tập các điều khiển chấp nhận được.
- Điều kiện tối ưu cần bậc nhất và đủ bậc hai: Thiết lập hệ điều kiện cực trị cho bài toán điều khiển thời gian, đưa ra cấu trúc của phương trình liên hợp trên khoảng thời gian biến thiên $(0, T^*)$ và khảo sát các điều kiện đủ đảm bảo tính tối ưu cục bộ ngặt.
Chương 4: Bài toán điều khiển biên tối ưu (An Optimal Boundary Control Problem)
Chương 4 khảo sát bài toán điều khiển biên phi tuyến, trong đó biến điều khiển tác động trực tiếp lên biên $\Gamma$ của miền không gian $\Omega$:
- Tính giải được của hệ Navier-Stokes-Voigt với điều kiện biên không thuần nhất: Khảo sát hệ với điều kiện biên $y = g$ trên $\Gamma \times (0, T)$. Sử dụng toán tử nâng (lifting operator) từ không gian biên $H^{1/2}(\Gamma)$ vào $H^1(\Omega)$ để chuyển bài toán về dạng thuần nhất tương đương và chứng minh tính tồn tại duy nhất nghiệm yếu.
- Sự tồn tại nghiệm điều khiển biên tối ưu: Khẳng định sự tồn tại nghiệm tối ưu biên $\bar{u}_\Gamma$ cho phiếm hàm mục tiêu toàn phương với các điều kiện tương thích đối với dữ kiện biên.
- Hệ điều kiện tối ưu cần bậc một và bậc hai: Thiết lập phương trình liên hợp biên, tính toán đạo hàm Fréchet/đạo hàm theo hướng trên không gian vết và đưa ra hệ điều kiện cần bậc một (Mục 4.4.1) và bậc hai (Mục 4.4.2).
- Điều kiện đủ bậc hai: Thiết lập định lý về điều kiện đủ bậc hai (Mục 4.5) đảm bảo nghiệm tối ưu biên thỏa mãn tính tăng bậc hai cục bộ theo chuẩn $L^2(\Sigma)$ (với $\Sigma = \Gamma \times (0, T)$).
Kết quả và những đóng góp mới
Luận án tự nêu các đóng góp khoa học mới trong phần Mở đầu, Lời cam đoan và Kết luận:
- Tính mới về đối tượng và kết quả: Là công trình đầu tiên nghiên cứu một cách hệ thống lý thuyết điều khiển tối ưu đối với hệ phương trình Navier-Stokes-Voigt phi dừng ba chiều trong miền bị chặn.
- Giải quyết trọn vẹn ba bài toán điều khiển cơ bản:
- Chứng minh sự tồn tại nghiệm tối ưu toàn cục cho cả ba dạng bài toán: điều khiển phân bố (P1), điều khiển thời gian (P2) và điều khiển biên không thuần nhất (P3).
- Thiết lập hệ phương trình liên hợp thích hợp và rút ra hệ điều kiện tối ưu cần bậc một dưới dạng bất đẳng thức biến phân tường minh cho từng bài toán.
- Thiết lập các định lý điều kiện tối ưu đủ bậc hai và chứng minh đặc trưng tăng bậc hai cục bộ của phiếm hàm mục tiêu trong lân cận nghiệm tối ưu.
- Đóng góp về phương pháp luận tiếp cận: Sử dụng phương pháp tính đạo hàm theo hướng trực tiếp thay cho phương pháp nhân tử Lagrange truyền thống, giúp cấu trúc điều kiện bậc hai có dạng tường minh, không phụ thuộc vào đạo hàm cấp hai phức tạp của hàm Lagrange. Đồng thời, kỹ thuật phản chứng xử lý thành công tính chất phi tuyến trong không gian ba chiều mà không vướng phải các khó khăn về tính chính quy như hệ Navier-Stokes cổ điển.
- Công bố khoa học: Các kết quả nghiên cứu trong luận án được công bố trong 02 bài báo trên các tạp chí quốc tế chuyên ngành thuộc danh mục ISI/Scopus:
- Bài báo [CT1] công bố trên tạp chí Numerical Functional Analysis and Optimization (nội dung Chương 2).
- Bài báo [CT2] công bố trên tạp chí Applied Mathematics and Optimization (nội dung Chương 3).
- Bài báo [CT3] hoàn thành bản thảo gửi đăng (nội dung Chương 4).
Hạn chế và hướng nghiên cứu tiếp
Dựa trên phạm vi và nội dung được tác giả trình bày tại phần Kết luận và Hướng phát triển (Conclusion and Future Work), các giới hạn và hướng mở rộng bao gồm:
- Hạn chế về cấu trúc điều khiển và điều kiện cần bậc hai: Trong bài toán điều khiển phân bố (Chương 2), do tác giả xét tập điều khiển chấp nhận được $U_{ad}$ là tập lồi, đóng bất kỳ (không giới hạn ở dạng hộp), khi điểm tối ưu nằm trên biên của $U_{ad}$, nón tiếp tuyến không thể xấp xỉ được bằng các hướng chấp nhận được như trường hợp ràng buộc hộp, dẫn đến việc chưa thiết lập được điều kiện tối ưu cần bậc hai cho trường hợp tập lồi tổng quát này.
- Hướng nghiên cứu tiếp được đề xuất:
- Nghiên cứu bài toán điều khiển tối ưu đối với hệ phương trình Navier-Stokes-Voigt có chứa ràng buộc trạng thái (state constraints) trên miền không - thời gian.
- Khảo sát bài toán điều khiển thưa (sparse optimal control problems) đối với hệ Navier-Stokes-Voigt, trong đó phiếm hàm mục tiêu chứa số hạng chuẩn $L^1$ để tạo ra các điều khiển có giá mang nhỏ gọn theo không gian hoặc thời gian.
- Xây dựng và phân tích các giải thuật số, giải pháp xấp xỉ phần tử hữu hạn để tính toán nghiệm số cho các bài toán điều khiển tối ưu đã thiết lập lý thuyết trong luận án.
Giá trị tham khảo
Luận án là tài liệu chuyên khảo có giá trị chuyên môn đối với:
- 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: Đặc biệt là hướng Phương trình vi phân và tích phân, Lý thuyết điều khiển tối ưu và Cơ học chất lưu lý thuyết.
- Giảng viên và nhà nghiên cứu giải tích PDE: Phần kỹ thuật đánh giá tiên nghiệm, tính chất toán tử tam tuyến $b(u, v, w)$, phép nhúng compact trừu tượng và kỹ thuật thiết lập phương trình liên hợp ngược thời gian (với điều kiện cuối $t = T$) tại Chương 1 và Chương 2 là mẫu chuẩn mực có thể áp dụng cho các mô hình $\alpha$ khác trong thủy động lực học.
- Các nhóm nghiên cứu mô phỏng số và tối ưu hóa dòng chảy: Cung cấp cơ sở giải tích lý thuyết vững chắc (tính đặt đúng, điều kiện nghiệm tối ưu, dạng toán tử chiếu) làm tiền đề cho việc xây dựng thuật toán tối ưu hóa số và bài toán tối ưu hóa điều khiển dòng chảy nhớt đàn hồi.
Câu hỏi thường gặp
1. Điểm khác biệt cơ bản về mặt toán học giữa hệ Navier-Stokes-Voigt và hệ Navier-Stokes kinh điển là gì?
Hệ Navier-Stokes-Voigt chứa thêm số hạng $-\alpha^2 \Delta y_t$ (trong đó $\alpha \ne 0$ là tham số độ dài đặc trưng cho tính đàn hồi của chất lưu). Số hạng này đóng vai trò chính quy hóa phương trình, giúp hệ có tính đặt đúng toàn cục và nghiệm yếu duy nhất trong không gian ba chiều $W^{1,2}(0, T; V)$ mà không cần đặt thêm điều kiện biên nhân tạo nào ngoài điều kiện biên Dirichlet.
2. Tại sao tác giả không sử dụng phương pháp nhân tử Lagrange truyền thống để tìm điều kiện tối ưu?
Tác giả lựa chọn tính toán trực tiếp đạo hàm theo hướng của phiếm hàm mục tiêu thông qua ánh xạ điều khiển - trạng thái. Cách tiếp cận này tự nhiên, dựa trên ý tưởng cực trị giải tích cổ điển, giúp biểu diễn dạng tường minh của điều kiện tối ưu bậc hai mà không phải biểu diễn qua đạo hàm cấp hai của hàm Lagrange, dù giá trị định lượng cuối cùng của hai cách là tương đương.
3. Tại sao trong Chương 2 chưa thiết lập được điều kiện tối ưu cần bậc hai cho tập điều khiển tổng quát?
Vì tác giả chọn tập điều khiển chấp nhận được $U_{ad}$ là một tập lồi, đóng bất kỳ trong $L^2(Q)$. Trong trường hợp điểm cực trị $\bar{u}$ thuộc biên của $U_{ad}$, không thể xấp xỉ mọi hướng trong nón tiếp tuyến $T_{U_{ad}}(\bar{u}) \cap C(\bar{u})$ bằng dãy các hướng chấp nhận được thuộc $F_{U_{ad}}(\bar{u}) \cap C(\bar{u})$ như khi xét tập ràng buộc dạng hộp (box constraints).
4. Dạng phương trình liên hợp xuất hiện trong bài toán điều khiển phân bố (Chương 2) có cấu trúc như thế nào?
Phương trình liên hợp là hệ phương trình vi phân ngược thời gian đối với hàm véc-tơ $w$:
$$-w_t - \nu \Delta w + \alpha^2 \Delta w_t - (\bar{y} \cdot \nabla)w + (\nabla \bar{y})^T w + \nabla p = \alpha_Q(\bar{y} - y_Q)$$
với điều kiện solenoidal $\text{div } w = 0$, điều kiện biên thuần nhất $w|_\Gamma = 0$ và điều kiện tại thời điểm cuối $w(T) - \alpha^2 \Delta w(T) + \nabla q = \alpha_T(\bar{y}(T) - y_T)$.
5. Kết quả điều khiển tối ưu biên ở Chương 4 xử lý điều kiện biên không thuần nhất bằng kỹ thuật gì?
Tác giả sử dụng định lý vết và toán tử nâng liên tục (continuous lifting operator) từ không gian vết $H^{1/2}(\Gamma)$ vào không gian Sobolev $H^1(\Omega)$ để chuyển bài toán có điều khiển biên không thuần nhất về bài toán tương đương có điều kiện biên Dirichlet thuần nhất, từ đó áp dụng các định lý giải được trong không gian hàm solenoidal.
Kết luận
Luận án tiến sĩ của tác giả Trần Minh Nguyệt là công trình giải tích toán học nghiên cứu hệ thống về lý thuyết điều khiển tối ưu cho hệ phương trình Navier-Stokes-Voigt phi dừng ba chiều. Tác giả đã giải quyết trọn vẹn bài toán tồn tại nghiệm tối ưu, thiết lập hệ điều kiện tối ưu cần bậc một và điều kiện tối ưu đủ bậc hai cho ba lớp bài toán: điều khiển phân bố, điều khiển thời gian và điều khiển biên. Toàn bộ các kết quả lý thuyết trong luận án đều có chứng minh giải tích chính xác, nhất quán và đã được công bố trên các tạp chí quốc tế uy tín trong lĩnh vực tối ưu hóa và giải tích hàm ứng dụng.