Tổng quan về luận án

Luận án Tiến sĩ Toán học với đề tài "Some Optimal Control Problems for Navier-Stokes-Voigt Equations" (Một số bài toán điều khiển tối ưu đối với hệ phương trình Navier-Stokes-Voigt) do Nghiên cứu sinh Trần Minh Nguyệt thực hiện dưới sự hướng dẫn khoa học của GS. Cung Thế Anh tại Khoa Toán - Tin học, Trường Đại học Sư phạm Hà Nội (chuyên ngành Phương trình Vi phân và Tích phân, Mã số: 9.03, năm bảo vệ: 2019). Công trình nghiên cứu này đặt nền móng lý thuyết tiên phong trong lĩnh vực lý thuyết điều khiển tối ưu hiện đại dành cho hệ phương trình thủy động lực học phi tuyến mô tả chất lưu đàn nhớt tuyến tính không nén được (linear viscoelastic incompressible fluids).

                 +-------------------------------------------------------------+
                 |            HỆ PHƯƠNG TRÌNH NAVIER-STOKES-VOIGT             |
                 |      y_t - νΔy - α²Δy_t + (y·∇)y + ∇p = u,  ∇·y = 0         |
                 +-------------------------------------------------------------+
                                                |
            +-----------------------------------+-----------------------------------+
            |                                   |                                   |
            v                                   v                                   v
+-----------------------+           +-----------------------+           +-----------------------+
|    BÀI TOÁN (P1)      |           |    BÀI TOÁN (P2)      |           |    BÀI TOÁN (P3)      |
|  Điều khiển phân bố   |           |  Điều khiển thời gian |           |   Điều khiển biên     |
| (Distributed Control) |           | (Time-Optimal Control)|           |   (Boundary Control)  |
+-----------------------+           +-----------------------+           +-----------------------+
            |                                   |                                   |
            +-----------------------------------+-----------------------------------+
                                                |
                                                v
                 +-------------------------------------------------------------+
                 |                    HỆ THỐNG PHÂN TÍCH TOÁN HỌC              |
                 |  - Tính giải được yếu trong không gian W^{1,2}(0,T; V)      |
                 |  - Hệ liên hợp phân rã ngược thời gian                      |
                 |  - Điều kiện cần bậc một (Bất đẳng thức biến phân)          |
                 |  - Điều kiện đủ bậc hai (Tăng trưởng bậc hai q(z,h) > 0)    |
                 +-------------------------------------------------------------+

Về bối cảnh khoa học, hệ phương trình Navier-Stokes-Voigt (NSV) được giới thiệu lần đầu bởi Oskolkov (1973) nhằm mô tả chuyển động của chất lưu nhớt đàn hồi Kelvin-Voigt. Sau đó, Cao, Lunasin và Titi (2006) đã đề xuất mô hình này như một phương pháp chỉnh hóa hữu hiệu ($\alpha$-regularization) đối với hệ phương trình Navier-Stokes ba chiều phục vụ mô phỏng số trực tiếp (DNS). Điểm ưu việt vượt trội của hệ NSV so với các mô hình $\alpha$ khác (như Camassa-Holm hay Leray-$\alpha$) nằm ở chỗ: hệ bảo toàn tính đặt chỉnh toàn cục (global well-posedness) trong không gian ba chiều mà không cần áp đặt bất kỳ điều kiện biên nhân tạo phức tạp nào ngoài điều kiện biên Dirichlet thuần nhất hoặc không thuần nhất tiêu chuẩn.

Khoảng trống nghiên cứu (research gap): Trước thời điểm luận án được công bố, các công trình quốc tế chủ yếu tập trung vào sự tồn tại nghiệm, tính chính quy và tập hút (attractors) của hệ NSV trong miền giải tích bị chặn hoặc phân rã nghiệm trong toàn không gian $\mathbb{R}^3$ (tiêu biểu như các nghiên cứu của Oskolkov, Cao, Lunasin, Titi, Kalantarov, Zelik). Trong khi đó, lý thuyết điều khiển tối ưu (optimal control theory) cho hệ động lực chất lưu Navier-Stokes kinh điển đã phát triển mạnh mẽ qua các công trình của Fursikov, Abergel, Temam, Casas, Tröltzsch, Wachsmuth, Barbu, nhưng hoàn toàn chưa có bất kỳ công bố nào giải quyết bài toán điều khiển tối ưu cho hệ Navier-Stokes-Voigt ba chiều (3D). Luận án đã lấp đầy khoảng trống học thuật này bằng cách xây dựng hệ thống giải tích hoàn chỉnh cho ba bài toán then chốt:

  1. Bài toán (P1) - Điều khiển phân bố (Distributed Optimal Control): Cực tiểu hóa hàm mục tiêu dạng toàn phương với biến điều khiển $u$ đóng vai trò trường lực ngoại lai thuộc tập lồi, đóng, khác rỗng $U_{ad} \subset L^2(Q)$.
  2. Bài toán (P2) - Điều khiển thời gian tối ưu (Time-Optimal Control): Đưa trạng thái vận tốc dòng chất lưu từ cấu hình ban đầu về trạng thái mục tiêu trong thời gian ngắn nhất dưới ràng buộc điều khiển lồi bất kỳ.
  3. Bài toán (P3) - Điều khiển biên tối ưu (Boundary Optimal Control): Tối ưu hóa trường vận tốc thông qua tác động điều khiển đặt trên biên $\Gamma$ thỏa mãn các điều kiện tương thích thủy động lực học.

Hệ thống giả thuyết và câu hỏi nghiên cứu định lượng của luận án:

  • Câu hỏi nghiên cứu 1 (RQ1): Dưới cấu trúc metric của không gian hàm Sobolev-Bochner, nghiệm yếu của phương trình trạng thái NSV 3D có tính giải được duy nhất và phụ thuộc liên tục vào hàm điều khiển hay không?
  • Câu hỏi nghiên cứu 2 (RQ2): Phương trình liên hợp (adjoint equation) tương ứng với số hạng điều chỉnh thời gian $-\alpha^2 \Delta y_t$ có dạng toán học như thế nào để thiết lập điều kiện cần bậc một thông qua đạo hàm theo hướng?
  • Câu hỏi nghiên cứu 3 (RQ3): Làm thế nào để thiết lập điều kiện đủ bậc hai tường minh (explicit second-order sufficient conditions) đảm bảo tính cực tiểu địa phương chặt với tốc độ tăng trưởng bậc hai mà không phụ thuộc vào phương pháp nhân tử Lagrange trừu tượng?

Khung lý thuyết của luận án tích hợp chặt chẽ: Lý thuyết phương trình vi phân đạo hàm riêng phi tuyến hiện đại, Không gian hàm Solenoidal ($H, V$), Lý thuyết nhúng compact Sobolev-Aubin-Lions-Simon, và Giải tích lồi vô hạn chiều.

Literature Review và Positioning

Lịch sử nghiên cứu điều khiển tối ưu các phương trình trong cơ học chất lưu bắt đầu từ những năm 1980 với các công trình nền tảng của Fursikov (1980, 1981), khi ông chứng minh các định lý tồn tại nghiệm tối ưu cho phương trình Navier-Stokes. Bước ngoặt lớn xuất hiện trong công trình kinh điển của Abergel và Temam (1990), nơi hệ điều kiện tối ưu bậc một lần đầu tiên được thiết lập cho dòng chảy chất lưu không nén được. Kể từ đó, nhiều nhánh nghiên cứu chuyên sâu đã hình thành:

                            TIẾN TRÌNH PHÁT TRIỂN HỌC THUẬT
                            
  Oskolkov (1973)               Fursikov (1980s)              Abergel & Temam (1990)
  Khởi xướng hệ NSV     --->   Tồn tại nghiệm điều khiển  --->   Hệ điều kiện tối ưu
  cho chất lưu đàn nhớt         Navier-Stokes                  cho dòng chất lưu
                                      |
                                      v
  Cao, Lunasin, Titi (2006)     Barbu (1998), Casas & Tröltzsch (2000s)
  Mô hình hóa α-regularization --->   Điều khiển tối ưu thời gian,
  cho 3D Navier-Stokes DNS           biên & phân bố cho 2D/3D NS
                                      |
                                      v
                    LUẬN ÁN TIẾN SĨ TRẦN MINH NGUYỆT (2019)
                    - Tồn tại toàn cục nghiệm điều khiển tối ưu 3D NSV
                    - Hệ liên hợp ngược chiều thời gian
                    - Điều kiện đủ bậc hai tường minh qua đạo hàm trực tiếp
  1. Nhóm nghiên cứu điều khiển phân bố (Distributed Control): Được phát triển bởi Casas (1997), Hinze (2000), Tröltzsch (2010), và đặc biệt là D. Wachsmuth (2006). Hầu hết các công trình này nghiên cứu phương trình Navier-Stokes 2D hoặc các trường hợp đặc biệt của 3D khi không có ràng buộc trạng thái. Trong trường hợp có ràng buộc trạng thái, Wang (2004) và Liu (2008) đã khảo sát các cấu trúc phức tạp của không gian nhân tử Lagrange đo được.
  2. Nhóm nghiên cứu điều khiển tối ưu thời gian (Time-Optimal Control): Barbu (1998) và Fernández-Cara (1999) đã mở đường cho bài toán điều khiển thời gian tối ưu đối với hệ Navier-Stokes kinh điển bằng cách sử dụng các nguyên lý cực đại Pontryagin dạng vô hạn chiều.
  3. Nhóm nghiên cứu điều khiển biên (Boundary Optimal Control): Nghiên cứu trên hệ dừng bởi Gunzburger, Hou, Svobodny (1991, 1992) và trên hệ không dừng bởi Raymond (2006), Casas (2006), Fursikov (2004).

Tranh luận khoa học then chốt: Đối với hệ Navier-Stokes 3D tiêu chuẩn, bài toán tồn tại nghiệm toàn cục và duy nhất của nghiệm trơn (smooth/weak solutions) vẫn là một bài toán thiên niên kỷ mở. Do đó, khi nghiên cứu điều khiển tối ưu cho hệ Navier-Stokes 3D, các nhà toán học buộc phải:

  • Giả định nghiệm tồn tại cục bộ theo thời gian, hoặc
  • Giới hạn dữ liệu ban đầu phải "đủ nhỏ", hoặc
  • Nghiên cứu trong không gian nghiệm yếu không duy nhất (non-unique weak solutions) dẫn đến việc sử dụng tập hút quỹ đạo vô cùng phức tạp.

Ngược lại, mô hình Navier-Stokes-Voigt với tham số đàn hồi $\alpha > 0$ đã khắc phục triệt để bế tắc này nhờ số hạng tán xạ bậc cao $-\alpha^2 \Delta y_t$. Nghiên cứu của Cao, Lunasin, và Titi (2006) chứng minh rằng hệ NSV 3D luôn có nghiệm yếu duy nhất trên mọi khoảng thời gian $[0, T]$.

Định vị của luận án: So với công trình của D. Wachsmuth (2006) trên hệ 2D Navier-Stokes với ràng buộc hộp (box constraints), luận án này mở rộng bài toán sang không gian 3D thực tế hơn đối với hệ NSV, đồng thời mở rộng tập điều khiển chấp nhận được $U_{ad}$ từ dạng ràng buộc hộp sang tập lồi, đóng, khác rỗng bất kỳ trong không gian Hilbert $L^2(Q)$. So với phương pháp giải tích đa trị của Barbu (1998) trên hệ Navier-Stokes, luận án đưa ra một phương pháp tiếp cận trực tiếp và tự nhiên hơn dựa trên việc tính đạo hàm theo hướng của phiếm hàm mục tiêu kết hợp cấu trúc giải tích lồi của nón tiếp tuyến (tangent cone) và nón pháp tuyến (normal cone).

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đã mở rộng và hoàn thiện lý thuyết điều khiển tối ưu phi tuyến cho lớp phương trình đạo hàm riêng giả hyperbolic-parabolic thông qua các đóng góp then chốt:

  1. Xây dựng Lý thuyết Giải tích Khả vi của Ánh xạ Điều khiển - Trạng thái ($S: u \mapsto y$): Luận án chứng minh rằng ánh xạ nghiệm $S$ khả vi theo hướng từ $L^2(Q)$ vào không gian năng lượng $W^{1,2}(0, T; V)$, với đạo hàm theo hướng $S'(\bar{u}; h) = z$ chính là nghiệm yếu duy nhất của hệ phương trình tuyến tính hóa: $$z_t - \nu\Delta z - \alpha^2 \Delta z_t + (z \cdot \nabla)\bar{y} + (\bar{y} \cdot \nabla)z + \nabla q = h$$
  2. Thiết lập Cấu trúc Toán học của Phương trình Liên hợp (Adjoint System): Do sự xuất hiện của số hạng $\alpha^2 \Delta y_t$, phương trình liên hợp phải mang dấu ngược lại đối với các đạo hàm thời gian để bảo toàn tính ổn định ngược: $$-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 cuối kỳ chặt chẽ trong không gian đối ngẫu $V'$: $$(I + \alpha^2 A)w(T) = \alpha_T (\bar{y}(T) - y_T)$$ Toán tử $I + \alpha^2 A: V \to V'$ là một đẳng cấu Banach-Hilbert, đảm bảo việc xác định duy nhất trạng thái cuối $w(T) \in V$.
  3. Bất đẳng thức Biến phân Tối ưu Bậc một: Thiết lập điều kiện cần thông qua nón cực của các tiếp tuyến $T_{U_{ad}}(\bar{u})$: $$\iint_Q (w + \gamma \bar{u}) \cdot h , dxdt \ge 0, \quad \forall h \in T_{U_{ad}}(\bar{u})$$ tương đương với biểu diễn hình chiếu tổng quát: $$\iint_Q (w + \gamma \bar{u}) \cdot (v - \bar{u}) , dxdt \ge 0, \quad \forall v \in U_{ad}$$
                CƠ CHẾ PHÂN TÍCH TĂNG TRƯỞNG BẬC HAI q(z,h)
                
+-------------------------------------------------------------------------+
|     Dạng toàn phương:                                                   |
|     q(z,h) = α_T/2 ||z(T)||^2 + α_Q/2 ||z||^2 + γ/2 ||h||^2            |
|              - ∫_0^T b(z(t), z(t), w(t)) dt                             |
+-------------------------------------------------------------------------+
                                    |
                                    v
+-------------------------------------------------------------------------+
|     Giả thiết điều kiện đủ bậc hai: q(z,h) > 0, ∀h ∈ (T_U ∩ C)\{0}      |
+-------------------------------------------------------------------------+
                                    |
                                    v
+-------------------------------------------------------------------------+
|     KẾT QUẢ ĐỘT PHÁ: TÍNH TĂNG TRƯỞNG TOÀN PHƯƠNG ĐỊA PHƯƠNG            |
|     J(v) ≥ J(v̄) + ε ||u - ū||_{L^2(Q)}^2                                 |
+-------------------------------------------------------------------------+

Khung phân tích độc đáo

Khung phân tích của luận án kết hợp hài hòa ba trụ cột lý thuyết toán học cao cấp:

  1. Lý thuyết Không gian Solenoidal và Toán tử Navier-Stokes: Sử dụng không gian $V = {u \in H_0^1(\Omega)^3 : \text{div } u = 0}$ và $H = {u \in L^2(\Omega)^3 : \text{div } u = 0}$, cùng dạng tam tuyến tính liên tục $b(u,v,w) = \sum_{i,j=1}^3 \int_\Omega u_i \frac{\partial v_j}{\partial x_i} w_j dx$.
  2. Phương pháp Tiếp cận Trực tiếp qua Đạo hàm Phân tích (Direct Directional Derivative Approach): Thay vì sử dụng hàm Lagrange vô hạn chiều đòi hỏi tính chính quy cao của nhân tử, luận án khai triển trực tiếp hiệu $J(y, u) - J(\bar{y}, \bar{u})$ theo chuỗi Taylor dạng biến phân: $$J(y_k, u_k) = J(\bar{y}, \bar{u}) + t_k \iint_Q (w + \gamma \bar{u}) \cdot h_k , dxdt + t_k^2 q(z_k, h_k) + t_k^2 Z_k$$ trong đó $Z_k \to 0$ khi $t_k \to 0$, từ đó làm xuất hiện trực tiếp dạng song tuyến tính $q(z, h)$.
  3. Giải tích Lồi Hiện đại trên Tập Ràng buộc Tổng quát: Sử dụng tính chất của nón pháp tuyến $N_U(u)$ và nón tiếp tuyến $T_U(u)$, cho phép chứng minh định lý bao hàm cho cả trường hợp điều khiển bị chặn trong hình cầu đóng $\bar{B}_R \subset L^2(Q)$ và siêu hộp (multidimensional box constraints).

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Luận án tuân thủ nghiêm ngặt quan điểm nhận thức luận Duy lý Toán học Suy diễn (Mathematical Deductive Rationalism). Thiết kế nghiên cứu dựa trên việc giải tích hóa các phương trình đạo hàm riêng phi tuyến vô hạn chiều trên miền không gian compact $\Omega \subset \mathbb{R}^3$ có biên Lipschitz $\Gamma$, trong khoảng thời gian hữu hạn $(0, T)$.

+-------------------------------------------------------------------------+
|                  KHÔNG GIAN VÀ CẤU TRÚC GIẢI TÍCH                       |
+-------------------------------------------------------------------------+
|  Không gian trạng thái:      y ∈ W^{1,2}(0,T; V) ⊂ C([0,T]; V)          |
|  Không gian đạo hàm thời gian: y_t ∈ L^2(0,T; V)                        |
|  Không gian đối ngẫu:        V' = (H_0^1,div=0)'                        |
|  Không gian điều khiển:      u ∈ U_{ad} ⊆ L^2(Q) = L^2(0,T; L^2(Ω)³)    |
|  Hình lăng trụ thời-không:    Q = Ω × (0, T)                             |
+-------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

Quy trình chứng minh được chuẩn hóa qua các công cụ giải tích hàm phức tạp:

Công cụ / Bất đẳng thức Biểu thức toán học / Định mức Vai trò trong chứng minh
Bất đẳng thức Poincaré $|u|{L^2(\Omega)} \le C\Omega |\nabla u|_{L^2(\Omega)}$ Đồng nhất chuẩn $|u|_1$ và chuẩn $|u|$ trong không gian $V \subset H_0^1(\Omega)$
Bất đẳng thức Ladyzhenskaya $|v|{L^4(\Omega)} \le 2^{1/4} |v|{L^2(\Omega)}^{1/4} |\nabla v|_{L^2(\Omega)}^{3/4}$ Đánh giá chặn số hạng đối lưu phi tuyến $b(u,v,w)$ trong không gian 3 chiều
Định lý Nhúng Aubin-Lions-Simon $W^{1,2}(0,T; V) \hookrightarrow L^2(Q)$ (compact) Đảm bảo tính hội tụ mạnh của dãy nghiệm xấp xỉ trong phép lấy giới hạn phi tuyến
Bổ đề Gronwall tích phân $\eta(t) \le e^{\int_0^t \phi(s)ds} [\eta(0) + \int_0^t \psi(s)ds]$ Đánh giá tiên nghiệm tính bị chặn đều của năng lượng $|y_n(t)|^2 + \alpha^2 |\nabla y_n(t)|^2$
Tính liên tục yếu của $B(y,y)$ $B(y_n, y_n) \rightharpoonup B(y,y)$ trong $L^2(0,T; V')$ Chuyển qua giới hạn đối với số hạng phi tuyến $(y \cdot \nabla)y$ khi $y_n \rightharpoonup y$

Data và phân tích toán học

Trong nghiên cứu lý thuyết toán học thuần túy, "dữ liệu" là các cấu trúc biến phân và các dãy cực tiểu hóa (minimizing sequences).

  1. Thiết lập dãy cực tiểu hóa: Xét dãy ${u_n} \subset U_{ad}$ sao cho $\lim_{n \to \infty} J(y_n, u_n) = \inf_{(y,u)} J(y,u) = \bar{J}$.
  2. Đánh giá tiên nghiệm (A priori estimates): Nhân phương trình trạng thái với $y_n(t) \in V$, sử dụng tính chất triệt tiêu $b(y_n, y_n, y_n) = 0$, thu được: $$\frac{1}{2}\frac{d}{dt}\left(|y_n(t)|^2 + \alpha^2 |y_n(t)|^2\right) + \nu |y_n(t)|^2 = (u_n(t), y_n(t)) \le \frac{\nu}{2}|y_n(t)|^2 + \frac{C^2}{2\nu}|u_n(t)|^2$$ Từ đó suy ra tính bị chặn đều của ${y_n}$ trong $L^\infty(0,T; H) \cap L^2(0,T; V)$. Nhân tiếp phương trình với $y_{nt}(t) \in V$, sử dụng tính chặn của dạng tam tuyến tính: $$|b(y_n, y_n, y_{nt})| \le C |y_n|^{1/4} |y_n|^{3/4} |y_n|^{1/4} |y_n|1^{3/4} |y{nt}| \le \frac{\alpha^2}{2} |y_{nt}|^2 + C |y_n|^4$$ suy ra ${y_{nt}}$ bị chặn đều trong $L^2(0,T; V)$.
  3. Trích dãy hội tụ yếu và tính đóng yếu: Trích dãy con $u_{n'} \rightharpoonup \bar{u}$ trong $L^2(Q)$ và $y_{n'} \rightharpoonup \bar{y}$ trong $W^{1,2}(0,T;V)$. Do $U_{ad}$ lồi và đóng trong $L^2(Q)$ nên $U_{ad}$ đóng yếu, suy ra $\bar{u} \in U_{ad}$.
  4. Tính nửa liên tục dưới yếu (Weak Lower Semicontinuity): Do phiếm hàm $J(y,u)$ là phiếm hàm lồi và liên tục theo chuẩn nên nó nửa liên tục dưới yếu, dẫn đến $J(\bar{y}, \bar{u}) \le \liminf_{n' \to \infty} J(y_{n'}, u_{n'}) = \bar{J}$, chứng minh $\bar{u}$ là nghiệm tối ưu toàn cục.

Phát hiện đột phá và implications

Những phát hiện then chốt

Luận án đã đạt được 5 kết quả đột phá có giá trị học thuật xuất sắc:

+----------------------------------------------------------------------------------------------------+
|                                 5 PHÁT HIỆN ĐỘT PHÁ CỦA LUẬN ÁN                                    |
+----------------------------------------------------------------------------------------------------+
| [1] Tồn tại duy nhất nghiệm tối ưu toàn cục ū ∈ U_{ad} cho hệ NSV 3D không dừng                    |
| [2] Đặc trưng hóa giải tích hệ liên hợp ngược chiều thời gian với số hạng điều chỉnh α²Δw_t        |
| [3] Thiết lập công thức biểu diễn giải tích dạng hình chiếu cho điều khiển tối ưu                  |
|     - Ràng buộc hộp: ū_i(x,t) = P_{[u_{a,i}, u_{b,i}]}(-(1/γ) w_i(x,t))                            |
|     - Ràng buộc cầu: ū_i(x,t) = -(1/(γ + s)) w_i(x,t),  s ≥ 0                                      |
| [4] Chứng minh định lý điều khiển tối ưu thời gian và sự tồn tại thời gian tối ưu T*               |
| [5] Thiết lập định lý điều kiện đủ bậc hai với tốc độ tăng trưởng bậc hai nghiêm ngặt:             |
|     J(v) ≥ J(v̄) + ε ||u - ū||_{L^2(Q)}^2                                                           |
+----------------------------------------------------------------------------------------------------+
  1. Chứng minh Tồn tại Nghiệm Toàn cục trong 3D Không gian Thực: Khác với phương trình Navier-Stokes tiêu chuẩn (nơi bài toán 3D thường gặp bế tắc về tính duy nhất), nghiệm tối ưu của hệ Navier-Stokes-Voigt 3D tồn tại toàn cục và duy nhất trên bất kỳ khoảng thời gian hữu hạn nào mà không đòi hỏi giả thiết độ nhỏ của dữ liệu ban đầu.
  2. Khám phá Cơ chế Tán xạ Ngược của Phương trình Tuyến tính hóa và Liên hợp: Phát hiện ra rằng sự kết hợp giữa toán tử quán tính thời gian $y_t$ và tán xạ đàn nhớt $-\alpha^2 \Delta y_t$ tạo ra một cấu trúc đối ngẫu hoàn hảo, cho phép nghiệm của phương trình liên hợp đạt độ chính quy cao $w \in W^{1,2}(0,T;V) \cap C([0,T];V)$. Khi trường lực ngoại lai thuộc $W^{1,2}(0,T;V')$, nghiệm đạt tính chính quy $C^1([0,T];V)$.
  3. Biểu thức Giải tích Tường minh cho Nghiệm Tối ưu:
    • Đối với ràng buộc hộp: $\bar{u}i(x,t) = \min{u{b,i}(x,t), \max{u_{a,i}(x,t), -\frac{1}{\gamma}w_i(x,t)}}$.
    • Đối với không gian toàn phần $U_{ad} = L^2(Q)$: $\bar{u}(x,t) = -\frac{1}{\gamma}w(x,t)$ hầu khắp nơi.
  4. Giải quyết Trọn vẹn Bài toán Điều khiển Tối ưu Biên Không thuần nhất: Xây dựng thành công toán tử nâng (lifting operator) từ không gian biên Sobolev phân số $H^{1/2}(\Gamma)$ vào miền trong $H^1(\Omega)$, xử lý triệt để điều kiện tương thích div bằng 0 trên biên.
  5. Chứng minh Bằng Phản chứng Định lý Điều kiện Đủ Bậc hai: Thiết lập một chuỗi đánh giá vi phân cấp 2 chính xác thông qua bổ đề khai triển $y_k = \bar{y} + t_k z_k + \frac{t_k^2}{2} \delta_k + t_k^2 \eta_k$, chỉ ra rằng dạng toàn phương $q(z,h) > 0$ trên nón tới hạn $T_{U_{ad}}(\bar{u}) \cap C(\bar{u}) \setminus {0}$ là điều kiện đủ để nghiệm là cực tiểu địa phương chặt.

Implications đa chiều

  • Về mặt Lý thuyết Toán học: Đặt nền móng vững chắc cho giải tích biến phân trên các hệ phương trình cơ học chất lưu có tính đàn hồi nhớt. Kết quả mở rộng các định lý của Lions, Temam, Tröltzsch sang lớp toán tử Kelvin-Voigt.
  • Về mặt Phương pháp luận: Chứng minh rằng phương pháp đạo hàm trực tiếp hoàn toàn có thể thay thế phương pháp nhân tử Lagrange cổ điển trong việc rút ra các điều kiện tối ưu bậc hai tường minh, giảm bớt sự phụ thuộc vào các định lý hàm ẩn phức tạp trong không gian vô hạn chiều.
  • Về mặt Ứng dụng Thực tiễn và Công nghệ:
    • Kỹ thuật Giảm lực cản Thủy động học (Hydrodynamic Drag Reduction): Ứng dụng trong việc bơm phụ gia polymer mạch dài vào dòng chảy để kiểm soát nhiễu loạn và giảm ma sát nhớt trên vỏ tàu ngầm, máy bay.
    • Xử lý Ảnh số (Image Inpainting): Hệ NSV được sử dụng để nội suy, phục hồi các vùng ảnh bị khuyết tật hoặc xóa vật thể trong ảnh dựa trên cơ chế lan truyền vận tốc Navier-Stokes-Voigt.
    • Mô phỏng Dòng dung dịch Sinh học: Dòng máu chảy qua van tim nhân tạo hoặc mạch máu có tính đàn nhớt cao.

Limitations và Future Research

Hạn chế học thuật thừa nhận

  1. Giới hạn về Điều kiện Cần Bậc hai đối với Tập Ràng buộc Tổng quát: Luận án chỉ ra rằng đối với tập lồi $U_{ad}$ tổng quát, không thể xấp xỉ một hướng bất kỳ trong nón tiếp tuyến $T_{U_{ad}}(\bar{u}) \cap C(\bar{u})$ bằng một dãy các hướng khả thi trong $F_{U_{ad}}(\bar{u}) \cap C(\bar{u})$ khi nghiệm $\bar{u}$ nằm trên biên của $U_{ad}$. Do đó, điều kiện cần bậc hai chỉ thiết lập được đối với ràng buộc dạng hộp mà chưa tổng quát hóa cho mọi tập lồi đóng bất kỳ.
  2. Thiếu vắng Ràng buộc Trạng thái Điểm (Pointwise State Constraints): Nghiên cứu tập trung vào các ràng buộc trên biến điều khiển ($u \in U_{ad}$). Các bài toán có ràng buộc trực tiếp trên biến trạng thái (ví dụ: áp suất giới hạn hoặc gradient vận tốc không vượt quá ngưỡng phá hủy vật liệu) chưa được tích hợp vào hàm mục tiêu.
  3. Mô hình Xác định thuần túy (Deterministic Setting): Luận án chưa xét đến tác động của các nhiễu ngẫu nhiên trắng (Gaussian noise/Brownian motion) tác động lên hệ động lực chất lưu.

Chương trình nghiên cứu tương lai (Future Research Agenda)

                            ĐỊNH HƯỚNG PHÁT TRIỂN 5 - 10 NĂM
                            
  +-----------------------------------+     +-----------------------------------+
  |   Sparse Optimal Control (L1)     |     |  Feedback Stabilization & Riccati |
  |   Tối ưu hóa vị trí cảm biến      |     |  Điều khiển phản hồi thời gian    |
  |   thông qua chuẩn L^1 phi trơn    |     |  thực cho hệ NSV hỗn loạn         |
  +-----------------------------------+     +-----------------------------------+
                    |                                         |
                    +--------------------+--------------------+
                                         |
                                         v
  +-----------------------------------+     +-----------------------------------+
  |  Stochastic NSV Control           |     |  Scientific Machine Learning      |
  |  Điều khiển ngẫu nhiên với vi phân|     |  PINNs & Deep Operator Networks   |
  |  ngẫu nhiên Itô / McKean-Vlasov   |     |  giải bài toán ngược NSV          |
  +-----------------------------------+     +-----------------------------------+
  1. Bài toán Điều khiển Tối ưu Thưa (Sparse Optimal Control): Nghiên cứu hàm mục tiêu chứa số hạng chuẩn $|u|_{L^1(Q)}$ để tối ưu hóa việc bố trí hữu hạn các bộ truyền động (actuators) trên thực tế.
  2. Điều khiển Ổn định Phản hồi (Feedback Stabilization): Thiết lập phương trình Riccati phi tuyến vô hạn chiều để tìm hàm điều khiển phản hồi $u(t) = -K y(t)$ đưa vận tốc dòng chảy về trạng thái cân bằng ổn định tiệm cận.
  3. Phân tích Số và Ước lượng Sai số (Numerical Analysis & Error Estimates): Xây dựng sơ đồ phần tử hữu hạn không - thời gian (Space-Time Finite Element Method) rời rạc hóa bài toán tối ưu và chứng minh tốc độ hội tụ của nghiệm rời rạc về nghiệm liên tục.
  4. Mở rộng sang Hệ Điều khiển Tối ưu Ngẫu nhiên (Stochastic Optimal Control): Nghiên cứu hệ phương trình vi phân ngẫu nhiên NSV chịu kích động của tích phân Itô trong không gian Hilbert.

Tác động và ảnh hưởng

Tác động học thuật

Các kết quả nghiên cứu trong luận án đã được bình duyệt nghiêm ngặt và công bố trên các tạp chí quốc tế hàng đầu thuộc danh mục ISI/Scopus (Q1/Q2) trong lĩnh vực Tối ưu hóa và Toán ứng dụng:

  • Numerical Functional Analysis and Optimization (Taylor & Francis): Công bố bài báo [CT1] về bài toán điều khiển phân bố cho hệ NSV.
  • Applied Mathematics and Optimization (Springer): Công bố bài báo [CT2] về bài toán điều khiển thời gian tối ưu cho hệ NSV.
  • Công trình về điều khiển biên [CT3] được hoàn thiện và nộp đăng trên các tạp chí chuyên ngành phương trình đạo hàm riêng uy tín.

Các công trình này đã nhận được sự quan tâm và trích dẫn từ cộng đồng toán học giải tích quốc tế, khẳng định vị thế của trường phái giải tích phương trình đạo hàm riêng và điều khiển tối ưu Việt Nam tại Đại học Sư phạm Hà Nội.

Tác động thực tiễn và liên ngành

  • Kỹ nghệ Hàng không & Hàng hải: Cung cấp mô hình toán học giải thích hiện tượng giảm lực cản nhớt khi pha dung dịch polymer nồng độ thấp (hiệu ứng Toms), giúp thiết kế tối ưu hệ thống phun chất đàn nhớt chủ động trên bề mặt chuyển động.
  • Công nghiệp Polymer và Vật liệu Mềm: Hỗ trợ việc kiểm soát trường nhiệt - vận tốc trong quá trình đùn ép nhựa và sản xuất sợi tổng hợp, nơi chất lỏng phi Newton có tính nhớt đàn hồi chi phối toàn bộ động học.

Đối tượng hưởng lợi

+----------------------------------------------------------------------------------------------------+
|                                 ĐỐI TƯỢNG HƯỞNG LỢI VÀ GIÁ TRỊ NHẬN ĐƯỢC                           |
+----------------------------------------------------------------------------------------------------+
| 1. Nghiên cứu sinh & Nhà Toán học Lý thuyết:                                                       |
|    - Tiếp cận khung giải tích hoàn chỉnh về tính khả vi của ánh xạ điều khiển - trạng thái 3D NSV. |
|    - Phương pháp luận chứng minh điều kiện đủ bậc hai thông qua bất đẳng thức năng lượng.          |
|                                                                                                    |
| 2. Kỹ sư Tính toán Động lực học Chất lưu (CFD Engineers):                                          |
|    - Thuật toán giải quyết bài toán biên không thuần nhất bảo toàn tính solenoidal (div = 0).     |
|    - Công thức chiếu ū_i giúp lập trình nhanh các module điều khiển biên trong OpenFOAM/ANSYS.     |
|                                                                                                    |
| 3. Nhà nghiên cứu Thị giác Máy tính & Xử lý Ảnh (Computer Vision):                                 |
|    - Ứng dụng phương trình Navier-Stokes-Voigt trong thuật toán inpainting nội suy phục hồi ảnh số.|
|                                                                                                    |
| 4. Cơ quan Quản lý Khoa học & Viện Hàn lâm:                                                       |
|    - Cung cấp minh chứng cho sự phát triển của nghiên cứu cơ bản đạt chuẩn quốc tế tại Việt Nam.  |
+----------------------------------------------------------------------------------------------------+

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Trả lời: Đóng góp lý thuyết độc đáo nhất là việc thiết lập hệ điều kiện tối ưu bậc một và bậc hai tường minh cho hệ Navier-Stokes-Voigt 3D mà không cần sử dụng không gian nhân tử Lagrange trừu tượng. Luận án đã mở rộng lý thuyết điều khiển tối ưu hệ phương trình Navier-Stokes của Abergel & Temam (1990) và lý thuyết điều khiển tối ưu phương trình phi tuyến của Tröltzsch (2010) sang lớp phương trình chứa toán tử vi phân hỗn hợp bậc cao $-\alpha^2 \Delta y_t$, giải quyết triệt để rào cản về tính duy nhất nghiệm trong không gian 3 chiều.

2. Sự đổi mới về mặt phương pháp luận so với các nghiên cứu tiền nhiệm quốc tế là gì?

Trả lời: So với nghiên cứu của D. Wachsmuth (2006) trên hệ Navier-Stokes 2D vốn dựa nhiều vào cấu trúc hình học của ràng buộc hộp để lấy xấp xỉ nón tiếp tuyến, luận án đã:

  • Mở rộng không gian bài toán sang 3D nhờ tận dụng tính chất điều chỉnh đàn hồi $\alpha$.
  • Sử dụng phương pháp tính trực tiếp đạo hàm theo hướng của phiếm hàm mục tiêu kết hợp với đánh giá vi sai cấp 2 thông qua bổ đề Gronwall và bất đẳng thức Ladyzhenskaya. Phương pháp này cho phép nhận được trực tiếp dạng biến phân toàn phương $q(z,h)$ mà không cần chuyển qua dạng đạo hàm Fréchet cấp hai của hàm Lagrange.

3. Phát hiện nào trong luận án mang tính bất ngờ nhất về mặt toán học?

Trả lời: Phát hiện bất ngờ nhất là cơ chế tự điều chỉnh của phương trình liên hợp đối với số hạng vi phân thời gian. Trong khi phương trình trạng thái thuận chiều thời gian có số hạng $-\alpha^2 \Delta y_t$, thì trong phương trình liên hợp (vốn giải ngược thời gian từ $t = T$ về $t = 0$), số hạng này chuyển thành $+\alpha^2 \Delta w_t$. Khi đổi biến thời gian $\tau = T - t$, cấu trúc toán tử được bảo toàn hoàn hảo, cho phép nghiệm liên hợp $w$ thừa hưởng toàn bộ tính chính quy $C([0,T]; V) \cap W^{1,2}(0,T;V)$, đảm bảo tính đặt chỉnh của bài toán đối ngẫu vô hạn chiều.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) chặt chẽ không?

Trả lời: Có. Luận án xây dựng quy trình chứng minh toán học logic với cấu trúc tiên đề hóa hoàn chỉnh:

  • Định nghĩa rõ ràng các không gian hàm Sobolev và không gian phân số $H^s(\Gamma)$.
  • Từng bước đánh giá tiên nghiệm qua các tích phân năng lượng chi tiết (xem các phương trình từ (1.21) đến (1.24) và (2.6) đến (2.8)).
  • Quy trình lấy giới hạn yếu và chứng minh tính compact được dẫn dắt minh bạch thông qua định lý Rellich-Kondrachov và bổ đề Simon. Bất kỳ nhà nghiên cứu nào trong lĩnh vực PDE đều có thể kiểm tra và tái lập từng bước logic của các định lý.

5. Lộ trình nghiên cứu 10 năm được phác thảo từ luận án gồm những gì?

Trả lời: Luận án mở ra lộ trình phát triển ba giai đoạn:

  • Giai đoạn 1 (1-3 năm): Phát triển thuật toán số (Discontinuous Galerkin / Space-Time FEM) và xấp xỉ gradient liên hợp giải bài toán phân bố và điều khiển biên.
  • Giai đoạn 2 (3-6 năm): Nghiên cứu bài toán ổn định hóa phản hồi (Feedback Stabilization) thông qua lý thuyết phương trình Riccati đại số/vi phân cho hệ NSV.
  • Giai đoạn 3 (6-10 năm): Tích hợp mạng nơ-ron vật lý (Physics-Informed Neural Networks - PINNs) và Deep Operator Networks (DeepONet) để học toán tử điều khiển tối ưu thời gian thực cho hệ chất lưu đàn nhớt phức tạp trong công nghiệp.

Kết luận

Luận án Tiến sĩ của tác giả Trần Minh Nguyet đã hoàn thành xuất sắc các mục tiêu nghiên cứu đặt ra, tạo nên bước đột phá học thuật quan trọng trong giải tích điều khiển phương trình đạo hàm riêng. Các đóng góp cốt lõi bao gồm:

  1. Chứng minh Tồn tại Nghiệm Tối ưu Toàn cục: Khẳng định sự tồn tại của cặp trạng thái - điều khiển tối ưu $(\bar{y}, \bar{u}) \in W^{1,2}(0,T; V) \times U_{ad}$ cho bài toán điều khiển phân bố, bài toán điều khiển thời gian và bài toán điều khiển biên của hệ Navier-Stokes-Voigt 3D.
  2. Thiết lập Hệ Điều kiện Tối ưu Bậc một: Xây dựng tường minh phương trình liên hợp và chứng minh bất đẳng thức biến phân tối ưu $\iint_Q (w + \gamma \bar{u}) \cdot (v - \bar{u}) dxdt \ge 0$ trên tập điều khiển lồi đóng bất kỳ.
  3. Thiết lập Biểu thức Chiếu Giải tích: Đưa ra công thức hiện dạng cho điều khiển tối ưu $\bar{u}$ dưới các dạng ràng buộc hộp và ràng buộc hình cầu trong không gian $L^2(Q)$.
  4. Chứng minh Điều kiện Đủ Bậc hai Tường minh: Khẳng định tính tăng trưởng bậc hai địa phương chặt chẽ $J(v) \ge J(\bar{v}) + \varepsilon |u - \bar{u}|_{L^2(Q)}^2$ khi dạng toàn phương $q(z,h) > 0$ trên nón tới hạn.
  5. Giải quyết Tính Đặt chỉnh của Bài toán Biên Không thuần nhất: Xây dựng không gian nghiệm và toán tử nâng cho phép giải bài toán điều khiển biên với điều kiện tương thích bảo toàn khối lượng.
  6. Mở ra Ba Hướng Nghiên cứu Mới: Đặt nền tảng cho điều khiển thưa $L^1$, điều khiển phản hồi Riccati và giải tích số định lượng cho hệ chất lưu đàn nhớt Kelvin-Voigt trên trường quốc tế.