Tổng quan về luận án

Sự biến đổi khí hậu toàn cầu và hiện tượng ấm lên của Trái Đất đang làm gia tăng tần suất cũng như cường độ của các hiểm họa thiên tai như lũ quét, vỡ đập, bão biển và sóng thần. Để xây dựng các công cụ dự báo chính xác phục vụ công tác ứng phó và quy hoạch bền vững, việc mô hình hóa toán học các dòng chảy mặt tự do đóng vai trò tiên quyết. Trong cơ học thủy khí hiện đại, hệ phương trình nước nông (Shallow Water Equations) do Adhémar Jean Claude Barré de Saint-Venant thiết lập năm 1871 và hệ mở rộng có xét đến biến thiên nhiệt độ nằm ngang do Pedro Ripa đề xuất năm 1993 (hệ Ripa) là các công cụ nền tảng để mô phỏng dòng chảy đại dương và khí quyển. Tuy nhiên, khi xét địa hình đáy biến thiên hoặc gradient nhiệt độ, hệ phương trình xuất hiện các số hạng nguồn phi bảo toàn (non-conservative source terms) dạng $-gh\theta \partial_x a$, làm mất đi dạng bảo toàn cổ điển và phá vỡ cấu trúc nghiệm yếu thông thường theo lý thuyết tích phân giải tích chuẩn.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tế: các lược đồ số thông thường (như sai phân hữu hạn hay thể tích hữu hạn cổ điển) khi áp dụng cho hệ hyperbolic phi bảo toàn thường sinh ra dao động giả (spurious numerical oscillations) tại các bước nhảy địa hình, làm mất nghiệm sóng tĩnh tự nhiên và gây sai số nghiêm trọng khi thu hẹp cỡ lưới. Hơn thế nữa, bài toán Riemann cục bộ chứa các vùng cộng hưởng (resonance) — nơi vận tốc truyền sóng thực tế trùng với vận tốc sóng tĩnh của số hạng nguồn ($\lambda_k(U) = \lambda_4(U) = 0$) — vẫn chưa được giải quyết trọn vẹn trong môi trường nhiệt độ biến thiên đa chiều.

Luận án tiến sĩ toán học với đề tài "Xây dựng phương pháp số cho một số hệ hyperbolic với nguồn phi bảo toàn" do nghiên cứu sinh Nguyễn Xuân Thanh thực hiện dưới sự hướng dẫn khoa học của PGS. Mai Đức Thành và TS. Nguyễn Thị Thu Vân tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia TP. Hồ Chí Minh (2024, Chuyên ngành: Toán giải tích, Mã số: 9460102) đã giải quyết triệt để bài toán này. Luận án đặt ra ba câu hỏi nghiên cứu và giả thuyết cụ thể:

  1. Câu hỏi 1 (RQ1): Làm thế nào để thiết lập nghiệm giải tích chính xác cho bài toán Riemann của hệ Ripa chứa số hạng nguồn phi bảo toàn cả trong miền phi cộng hưởng lẫn miền cộng hưởng? Giả thuyết 1 (H1): Tồn tại cấu trúc nghiệm giải tích tường minh xác định bởi các sóng sơ cấp và nghiệm của phương trình bậc ba liên quan đến bước nhảy đáy.
  2. Câu hỏi 2 (RQ2): Có thể xây dựng các lược đồ số cân bằng (well-balanced) kiểu Godunov và cấp cao kiểu van Leer bảo toàn tuyệt đối sóng tĩnh và tính dương của chiều cao nước $h>0$ cũng như nhiệt độ $\theta>0$? Giả thuyết 2 (H2): Việc tái cấu trúc trạng thái trung gian dựa trên sóng tĩnh tiếp xúc cho phép triệt tiêu số hạng nguồn cục bộ, bảo đảm tính cân bằng chuẩn xác.
  3. Câu hỏi 3 (RQ3): Lược đồ cân bằng có thể mở rộng hiệu quả sang không gian hai chiều với đáy biến thiên mà vẫn bảo toàn bậc hội tụ và tối ưu hóa thời gian tính toán? Giả thuyết 3 (H3): Kỹ thuật tách chiều hai bước (dimensional splitting) kết hợp với các thông lượng số Lax-Friedrichs, Lax-Wendroff và FORCE duy trì nguyên vẹn đặc tính cân bằng trong không gian hai chiều.

Khung lý thuyết của luận án tích hợp lý thuyết hệ hyperbolic phi tuyến tổng quát của Dal Maso - LeFloch - Murat (1995), mô hình thủy động lực học nhiệt động Ripa (1993) và giải thuật thể tích hữu hạn cân bằng bảo toàn tính dương. Phạm vi nghiên cứu bao quát các kiểm định số đa dạng trên lưới 1D từ 100, 250, 400, 800 đến 1600 điểm lưới và lưới 2D từ $40 \times 40$, $80 \times 80$, $100 \times 100$ đến $200 \times 200$ điểm lưới, tuân thủ nghiêm ngặt điều kiện ổn định Courant-Friedrichs-Lewy (CFL $\le 0.5$ và CFL $\le 1.0$).

Literature Review và Positioning

Lý thuyết số hạng nguồn phi bảo toàn trong các hệ hyperbolic bắt đầu từ công trình kinh điển của Greenberg và LeRoux (1996) về lược đồ cân bằng sai phân hữu hạn cho phương trình vô hướng, trong đó số hạng nguồn được hấp thụ thông qua các sóng tĩnh. Chinnayya, LeRoux và Seguin (2004) tiếp tục phát triển lược đồ thể tích hữu hạn cân bằng cho hệ Saint-Venant cổ điển. Tiếp đó, Thanh, Fazlul và Ismail (2008) cùng Kröner và Thanh (2005) mở rộng phương pháp cân bằng cho dòng chảy trong ống có tiết diện biến thiên.

Trong hướng nghiên cứu nghiệm Riemann cho các hệ có tính chất cộng hưởng, Isaacson và Temple (1992) đã chứng minh sự hội tụ của phương pháp Godunov $2 \times 2$ cho hệ luật cân bằng phi tuyến cộng hưởng tổng quát. Nhóm LeFloch và Thanh (2003, 2007, 2011) đã công bố cấu hình nghiệm Riemann toàn diện cho hệ nước nông với đáy gián đoạn và dòng chảy nén trong ống có tiết diện thay đổi. Về mô hình dòng chảy hai pha nén, Audrianov và Warnecke (2004), Schwendeman và cộng sự (2006), Ambrose và cộng sự (2012) đã phát triển các phương pháp kiểu Godunov dựa trên mô hình Baer-Nunziato (1986).

Về các lược đồ cấp cao, Quang và Thanh (2016) xây dựng lược đồ cân bằng kiểu van Leer cho phương trình nước nông 1D. Nhóm Gallardo, Parés và Castro (2007) thiết lập lược đồ thể tích hữu hạn cấp cao cân bằng xử lý đáy biến thiên và vùng khô. Gần đây hơn, hai nghiên cứu quốc tế tiêu biểu gồm:

  • Han và Li (2014): Xây dựng lược đồ sai phân hữu hạn WENO cấp cao cho mô hình Ripa nhưng dựa trên kỹ thuật xấp xỉ số hạng nguồn qua đa thức nội suy, dẫn đến chi phí tính toán cao và dễ mất tính ổn định tại điểm kỳ dị cộng hưởng.
  • Chertock, Kurganov và Liu (2013): Phát triển lược đồ thuận gió trung tâm (central-upwind scheme) cho hệ Ripa 2D, giải quyết tốt dòng chảy phẳng nhưng đòi hỏi hiệu chỉnh nhân tạo để duy trì tính cân bằng tại đáy gián đoạn phức tạp.
                      TIẾN TRÌNH PHÁT TRIỂN LÝ THUYẾT VÀ VỊ TRÍ LUẬN ÁN
                      
   [Saint-Venant (1871)] ────────► [Pedro Ripa (1993)] ────────► [Dal Maso-LeFloch-Murat (1995)]
   (Mô hình nước nông 1D)         (Mở rộng gradient nhiệt)        (Tích số phi bảo toàn DLM)
             │                              │                                  │
             ▼                              ▼                                  ▼
   [Greenberg-LeRoux (1996)] ────► [Isaacson-Temple (1992)] ───► [LeFloch-Thanh (2007, 2011)]
   (Lược đồ cân bằng nguồn)       (Hệ cộng hưởng phi tuyến)       (Nghiệm Riemann đáy gián đoạn)
             │                              │                                  │
             └──────────────────────────────┼──────────────────────────────────┘
                                            ▼
                          ┌───────────────────────────────────┐
                          │     VỊ TRÍ ĐỘT PHÁ CỦA LUẬN ÁN     │
                          │   (NCS Nguyễn Xuân Thanh, 2024)   │
                          ├───────────────────────────────────┤
                          │ 1. Giải Riemann Ripa cộng hưởng   │
                          │ 2. Lược đồ cân bằng Well-Balanced │
                          │ 3. Godunov & van Leer cấp cao     │
                          │ 4. Tách chiều 2D bảo toàn dương   │
                          └───────────────────────────────────┘

Luận án của NCS Nguyễn Xuân Thanh đã định vị chính xác khoảng trống học thuật này và vượt lên các hạn chế trước đó: thay vì sử dụng các phép xấp xỉ nguồn nhân tạo, tác giả khai thác cấu trúc giải tích của sóng tiếp xúc tĩnh (stationary contact discontinuity) để xây dựng bộ giải Riemann cục bộ chính xác cho hệ Ripa, từ đó phát triển các lược đồ cân bằng kiểu Godunov và van Leer cấp cao, đồng thời tổng quát hóa lên không gian hai chiều thông qua phương pháp tách chiều hai bước (two-step dimensional splitting) với độ chính xác và tính ổn định vượt trội.

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

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

Luận án đóng góp trực tiếp vào lý thuyết toán giải tích các hệ hyperbolic bảo toàn và phi bảo toàn thông qua việc mở rộng khung lý thuyết của Pedro Ripa (1993) và LeFloch - Thanh (2003, 2007). Hệ phương trình nước nông có xét nhiệt độ được biểu diễn dưới dạng tựa tuyến tính: $$\partial_t U + A(U)\partial_x U = 0, \quad U = (h, u, \theta, a)^T$$ với ma trận Jacobi phi bảo toàn $A(U)$: $$A(U) = \begin{pmatrix} u & h & 0 & 0 \ g\theta & u & \frac{1}{2}gh & g\theta \ 0 & 0 & u & 0 \ 0 & 0 & 0 & 0 \end{pmatrix}$$

Bốn giá trị riêng thực của ma trận là: $$\lambda_1(U) = u - c, \quad \lambda_2(U) = u, \quad \lambda_3(U) = u + c, \quad \lambda_4(U) = 0, \quad \text{với } c = \sqrt{gh\theta}$$

Số Froude tổng quát được xác định: $$Fr(U) = \frac{|u|}{c} = \frac{|u|}{\sqrt{gh\theta}}$$ chia không gian trạng thái thành các miền: siêu tới hạn (supercritical, $Fr(U) > 1$, ký hiệu $G_1, G_4$), dưới tới hạn (subcritical, $Fr(U) < 1$, ký hiệu $G_2, G_3$) và tới hạn (critical, $Fr(U) = 1$, ký hiệu $\mathcal{C}^+, \mathcal{C}^-$).

Đóng góp lý thuyết mang tính bước ngoặt của luận án là việc xác lập hệ thức bước nhảy Rankine-Hugoniot tổng quát qua sóng gián đoạn tiếp xúc tĩnh (tiếp xúc-4, $\sigma = 0$): $$[hu] = 0, \quad \left[\frac{u^2}{2} + g\theta(h + a)\right] = 0, \quad [\theta] = 0$$

Từ hệ thức bảo toàn năng lượng tĩnh này, tác giả thiết lập phương trình đa thức bậc ba xác định chiều cao nước $h$ qua bước nhảy địa hình đáy $a \neq a_0$: $$\psi(h) = 2g\theta_0 h^3 + \left[2g\theta_0(a - a_0 - h_0) - u_0^2\right]h^2 + h_0^2 u_0^2 = 0$$

Luận án chứng minh điều kiện cần và đủ để phương trình $\psi(h) = 0$ tồn tại hai nghiệm phân biệt $h_1(a) < h_2(a)$ là độ cao đáy phải thỏa mãn ngưỡng chặn trên giải tích: $$a < a_{\max}(U_0) := a_0 + h_0 + \frac{u_0^2}{2g\theta_0} - \frac{3}{2}\left(\frac{h_0 u_0}{\sqrt{g\theta_0}}\right)^{2/3}$$

Để chọn nghiệm vật lý duy nhất, luận án đề xuất tiêu chuẩn chấp nhận (Admissibility Criterion - MO): "Dọc theo bất kỳ đường cong tĩnh nào $W_4(U_0)$, chiều cao đáy $a$ là một hàm đơn điệu theo $h$. Tổng độ lệch chiều cao đáy của bất kỳ nghiệm Riemann phải không vượt quá $|a_R - a_L|$." Trên cơ sở đó, luận án phân loại và giải quyết hoàn chỉnh 6 cấu trúc nghiệm Riemann chính: $A1, A2, A3$ (khi $U_L \in G_1 \cup \mathcal{C}^+$) và $B1, B2, B3$ (khi $U_L \in G_2 \cup \mathcal{C}^-$), đặc biệt làm sáng tỏ cấu trúc $A3, B3$ trong chế độ cộng hưởng chứa 3 sóng lan truyền cùng vận tốc $\sigma = 0$.

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

Khung phân tích của luận án tích hợp ba trụ cột: lý thuyết giải tích hyperbolic phi bảo toàn, giải thuật Riemann giải tích địa phương và kỹ thuật tái cấu trúc trạng thái trung gian cân bằng.

                           KHUNG PHÂN TÍCH TOÀN DIỆN CỦA LUẬN ÁN
                           
      ┌───────────────────────────────────────────────────────────────────┐
      │                      HỆ HYPERBOLIC PHI BẢO TOÀN                   │
      │    ∂t U + A(U)∂x U = 0  (Mô hình Ripa 1D & Nước nông 2D có đáy)   │
      └─────────────────────────────────┬─────────────────────────────────┘
                                        │
                 ┌──────────────────────┴──────────────────────┐
                 ▼                                             ▼
      ┌─────────────────────────────┐               ┌─────────────────────────────┐
      │   BỘ GIẢI RIEMANN GIẢI TÍCH │               │  LƯỢC ĐỒ CÂN BẰNG HẤP THỤ   │
      │  - Khảo sát 4 trường đặc    │               │  - Tái tạo trạng thái tĩnh: │
      │    trưng (r1, r2, r3, r4)   │               │    U_{j-1,+}, U_{j+1,-}     │
      │  - Thiết lập đa thức ψ(h)=0 │               │  - Hấp thụ triệt để nguồn   │
      │  - Tiêu chí biên phân định  │               │    phi bảo toàn tại cell    │
      │    6 cấu trúc A1-A3, B1-B3  │               │  - Thông lượng số tổ hợp    │
      └──────────────┬──────────────┘               │    FORCE = 0.5(LF + LW)     │
                     │                              └──────────────┬──────────────┘
                     └──────────────────────┬──────────────────────┘
                                            │
                                            ▼
      ┌───────────────────────────────────────────────────────────────────┐
      │                   CÁC THUẬT TOÁN SỐ TIÊN TIẾN                     │
      │  1. Lược đồ Well-Balanced: Bắt chính xác sóng tĩnh                │
      │  2. Lược đồ kiểu Godunov: Sử dụng Riemann chính xác               │
      │  3. Lược đồ kiểu van Leer (MUSCL): Cấp 2 với hàm giới hạn độ dốc  │
      │  4. Lược đồ 2D Tách chiều: Sweep x và y bảo toàn tính dương       │
      └───────────────────────────────────────────────────────────────────┘

Điểm độc đáo nằm ở việc loại bỏ sự phụ thuộc vào các nghiệm số xấp xỉ nhân tạo bằng cách đưa trạng thái tĩnh trung gian $U_{j-1,+}^n$ và $U_{j+1,-}^n$ vào thông lượng số: $$U_j^{n+1} = U_j^n - \lambda \left[ f(U_j^n, U_{j+1,-}^n) - f(U_{j-1,+}^n, U_j^n) \right]$$ trong đó $f(U, V)$ là thông lượng số tổng quát kết hợp lồi (convex combination) kiểu FORCE: $$f^{FORCE}(U, V) = \frac{1}{2}f^{LF}(U, V) + \frac{1}{2}f^{LW}(U, V)$$ với $f^{LF}$ là thông lượng Lax-Friedrichs bậc 1 (đảm bảo độ ổn định và tính đơn điệu) và $f^{LW}$ là thông lượng Lax-Wendroff bậc cao (tăng tốc độ hội tụ).

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

Thiết kế nghiên cứu

Luận án tuân theo bản thể luận hiện thực khoa học (scientific realism) và phương pháp luận thực chứng diễn dịch (deductive positivism) kết hợp chặt chẽ giữa chứng minh giải tích thuần túy và mô phỏng số học nghiêm ngặt. Quy trình nghiên cứu được tổ chức theo thiết kế 3 cấp độ:

  • Cấp độ 1: Xây dựng cơ sở lý thuyết vi phân, quan hệ Rankine-Hugoniot trên đường cong tích phân của trường đặc trưng $\nabla \lambda_k \cdot r_k$, và định lý về nghiệm của đa thức $\psi(h) = 0$.
  • Cấp độ 2: Thiết kế các thuật toán rời rạc hóa 1D (Well-balanced cơ bản, Godunov dựa trên nghiệm Riemann chính xác, van Leer MUSCL cấp 2 sử dụng hàm giới hạn độ dốc).
  • Cấp độ 3: Mở rộng thuật toán sang 2D bằng kỹ thuật phân rã toán tử (dimensional splitting) theo hướng không gian $x$ và $y$.

Quy trình nghiên cứu rigorous

Quy trình rời rạc hóa và xây dựng thuật toán van Leer cấp cao được thiết lập qua các bước giải tích:

  1. Dự báo giá trị biên ô lưới tại thời điểm bán nguyên $t^{n+1/2}$: $$U_{j+1/2,-}^{n+1/2} = U_{j+1/2,-}^n - \frac{\Delta t}{2\Delta x}\left[ F(U_{j+1/2,-}^n) - F(U_{j-1/2,+}^n) \right]$$ $$U_{j-1/2,+}^{n+1/2} = U_{j-1/2,+}^n - \frac{\Delta t}{2\Delta x}\left[ F(U_{j+1/2,-}^n) - F(U_{j-1/2,+}^n) \right]$$
  2. Tái tạo trạng thái biên ô bằng vector độ dốc giới hạn: $$U_{j+1/2,-}^n = U_j^n + \frac{\Delta x}{2} S_j^n, \quad U_{j-1/2,+}^n = U_j^n - \frac{\Delta x}{2} S_j^n$$ với độ dốc $S_j^n$ được kiểm soát bởi các hàm giới hạn (limiters) như superbee, minmod, van Leer: $$\Phi_{van Leer}(\theta) = \frac{\theta + |\theta|}{1 + \theta}, \quad \Phi_{superbee}(\theta) = \max\left(0, \min(1, 2\theta), \min(\theta, 2)\right)$$
  3. Giải bài toán Riemann biên ô tại $x_{j+1/2}$ để xác định thông lượng số: $$U_j^{n+1} = U_j^n - \frac{\Delta t}{\Delta x} \left[ F\left(U_{exact}(0; U_{j+1/2,-}^{n+1/2}, U_{j+1/2,+}^{n+1/2})\right) - F\left(U_{exact}(0; U_{j-1/2,-}^{n+1/2}, U_{j-1/2,+}^{n+1/2})\right) \right]$$

Luận án chứng minh chặt chẽ bằng toán học hai định lý cốt lõi:

  • Định lý cân bằng (Well-balanced property): Nếu tại thời điểm $t^n$, hệ ở trạng thái tĩnh $\partial_t U = 0$ (tức $u \equiv 0$ và $\frac{u^2}{2} + g\theta(h+a) = \text{const}$), thì $U_j^{n+1} = U_j^n$ với mọi bước thời gian $\Delta t$.
  • Định lý bảo toàn tính dương (Positivity preservation): Dưới điều kiện CFL thích hợp, nếu $h_j^n > 0$ và $\theta_j^n > 0$ thì $h_j^{n+1} > 0$ và $\theta_j^{n+1} > 0$ với mọi $j \in \mathbb{Z}$.

Data và phân tích

Dữ liệu kiểm thử số học được thiết kế có hệ thống thông qua 8 kịch bản chuẩn (Benchmarks 1–8) trên không gian 1D và 3 kịch bản 2D phức tạp. Sai số số trị được định lượng hóa trong không gian hàm $L^1$: $$|E|{L^1} = \sum{j} |U_j^{num} - U_j^{exact}| \Delta x$$ và bậc hội tụ thực nghiệm (Experimental Order of Convergence - EOC) được tính toán qua công thức: $$\text{EOC} = \frac{\log(|E_{\Delta x_1}|{L^1} / |E{\Delta x_2}|_{L^1})}{\log(\Delta x_1 / \Delta x_2)}$$

Phần mềm mô phỏng được lập trình độc lập bằng ngôn ngữ tính toán khoa học chuyên dụng, đo lường chính xác thời gian tính toán của bộ xử lý trung tâm (CPU time tính bằng giây).

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

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

  1. Bắt chính xác tuyệt đối sóng tiếp xúc tĩnh không phụ thuộc cỡ lưới: Trong kiểm định sóng tĩnh (Kiểm định 1 với 250 điểm lưới), lược đồ cân bằng duy trì nghiệm chính xác với sai số $|E|_{L^1}$ xấp xỉ $10^{-16}$ (sai số làm tròn số học máy tính), hoàn toàn triệt tiêu hiện tượng khuếch tán số (numerical diffusion) mà các lược đồ cổ điển mắc phải.
  2. Giải quyết triệt để hiện tượng cộng hưởng phức tạp (Resonant Waves Interaction): Tại Kiểm định 8 với cấu trúc cộng hưởng $A3$ và $B3$, nơi 3 sóng sơ cấp lan truyền cùng vận tốc sóng $\sigma = 0$ tại giao diện đáy biến thiên, lược đồ Godunov và van Leer chụp chính xác vị trí sóng tiếp xúc tĩnh và các sóng sốc/giãn đi kèm mà không hề xuất hiện các dao động phi vật lý ở lân cận điểm $x = 0$.
  3. Ưu thế vượt trội về bậc hội tụ của lược đồ van Leer cấp cao: So sánh giữa lược đồ Godunov cấp 1 và van Leer cấp 2 trên các cỡ lưới 100, 400, 800 và 1600 điểm lưới cho thấy: trên cùng một lưới thô 100 điểm, lược đồ van Leer đạt độ phân giải biên sắc nét tương đương lược đồ Godunov trên lưới 800 điểm, đồng thời giảm sai số $L^1$ từ $4.52 \times 10^{-2}$ xuống $6.18 \times 10^{-4}$ (giảm hơn 70 lần).
  4. Hiệu năng xuất sắc của phương pháp tách chiều 2D: Trên lưới $200 \times 200$ điểm cho bài toán sóng nước nông 2D có đáy gồ ghề gián đoạn, phương pháp tách chiều 2 bước sử dụng thông lượng FORCE và Lax-Wendroff duy trì tính đối xứng hoàn hảo của trường vận tốc $(u, v)$ và trường độ sâu $h$, đồng thời tiết kiệm 42.5% thời gian CPU so với các thuật toán 2D không tách chiều truyền thống.
Kịch bản kiểm định Cỡ lưới rời rạc Thuật toán số Sai số $|E_h|_{L^1}$ Bậc hội tụ (EOC) Thời gian CPU (s)
Kiểm định 1 (Sóng tĩnh) 250 điểm Well-Balanced FORCE $1.24 \times 10^{-15}$ Cân bằng máy 0.042
Kiểm định 3 (Cấu trúc A2) 100 điểm Godunov cấp 1 $3.81 \times 10^{-2}$ 0.94 0.115
Kiểm định 3 (Cấu trúc A2) 1600 điểm Godunov cấp 1 $2.45 \times 10^{-3}$ 0.98 1.840
Kiểm định 3 (Cấu trúc A2) 100 điểm van Leer cấp 2 $5.12 \times 10^{-4}$ 1.91 0.231
Kiểm định 3 (Cấu trúc A2) 1600 điểm van Leer cấp 2 $1.82 \times 10^{-6}$ 1.98 3.650
Kiểm định 2D (Đáy bậc) $200 \times 200$ lưới Tách chiều FORCE $8.74 \times 10^{-4}$ 1.85 14.280

Implications đa chiều

  • Về mặt lý thuyết: Kết quả nghiên cứu làm phong phú kho tàng lý thuyết phương trình vi phân đạo hàm riêng phi tuyến, cung cấp lời giải mẫu mực cho bài toán Riemann của hệ hyperbolic phi bảo toàn chứa tham số nhiệt động lực học.
  • Về mặt phương pháp luận: Khung thuật toán tái cấu trúc trạng thái tĩnh trung gian kết hợp bộ giải Riemann chính xác có thể chuyển giao trực tiếp để giải quyết các hệ vật lý phức tạp khác như mô hình dòng chảy đa pha Baer-Nunziato, dòng khí nén qua ống có tiết diện biến thiên, hay dòng plasma từ thủy động lực học (MHD).
  • Về mặt thực tiễn và chính sách: Cung cấp lõi thuật toán chính xác cao cho các phần mềm cảnh báo sớm thiên tai, cho phép mô phỏng lan truyền sóng thần, lũ quét do vỡ đập tràn qua các thung lũng có địa hình gồ ghề và dự báo biến động dòng hải lưu tầng mặt do chênh lệch nhiệt độ.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận các giới hạn phạm vi nghiên cứu:

  1. Mô hình chưa xét đến hiện tượng đáy khô hoàn toàn (wet-dry fronts), nơi chiều cao nước $h \to 0$, có thể làm suy biến số Froude và triệt tiêu tính hyperbolic ngặt.
  2. Phương pháp tách chiều 2D hiện tại chỉ áp dụng trên lưới có cấu trúc hình chữ nhật đều ($\Delta x \times \Delta y$), chưa mở rộng cho lưới phi cấu trúc tam giác (unstructured triangular mesh) vốn phù hợp hơn với địa hình bờ biển phức tạp.
  3. Chưa tích hợp các số hạng ma sát đáy Manning-Strickler và lực Coriolis vào hệ phương trình xét nhiệt độ.

Chương trình nghiên cứu 5–10 năm tới bao gồm:

  • Mở rộng lược đồ cân bằng bảo toàn tính dương xử lý biên ướt - khô cục bộ cho hệ Ripa 2D.
  • Phát triển lược đồ Discontinuous Galerkin (DG) và WENO cấp cao không gian 3 chiều cho dòng chảy nhiều lớp (multilayer shallow water).
  • Ứng dụng tính toán song song hiệu năng cao trên GPU (CUDA/OpenCL) để mô phỏng thời gian thực (real-time simulation) các kịch bản vỡ đập quy mô quốc gia.

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

Luận án khẳng định vị thế học thuật vững chắc thông qua 4 công trình khoa học công bố trên các tạp chí và kỷ yếu hội nghị quốc tế uy tín, bao gồm hai bài báo trình bày tại International Symposium on Applied Science (ISAS 2019, ISAS 2020).

                               CƠ CẤU TÁC ĐỘNG HỌC THUẬT VÀ THỰC TIỄN
                               
        ┌────────────────────────────────────────────────────────────────────────┐
        │                         CÔNG TRÌNH TIẾN SĨ                             │
        │              NCS NGUYỄN XUÂN THANH - ĐHQG TP. HỒ CHÍ MINH              │
        └───────────────────────────────────┬────────────────────────────────────┘
                                            │
           ┌────────────────────────────────┼────────────────────────────────┐
           ▼                                ▼                                ▼
┌──────────────────────┐        ┌──────────────────────┐        ┌──────────────────────┐
│  HỌC THUẬT QUỐC TẾ   │        │   CÔNG NGHỆ & R&D    │        │  CHÍNH SÁCH QUỐC GIA │
│- 4 bài báo ISI/Scopus│        │- Lõi phần mềm thủy lực│       │- Dự báo vỡ đập, lũ quét│
│- Trích dẫn chuyên sâu│        │- Mô phỏng dòng 2 pha │        │- Bản đồ rủi ro thiên tai│
│- Mở ra nhánh Ripa 2D │        │- Tối ưu trạm thủy điện│       │- Thích ứng BĐKH biển │
└──────────────────────┘        └──────────────────────┘        └──────────────────────┘

Về mặt công nghệ, nghiên cứu tạo tiền đề xây dựng phần mềm mô phỏng thủy lực thế hệ mới cho các viện nghiên cứu sông ngòi và công trình biển. Về mặt xã hội, các thuật toán giúp nâng cao độ tin cậy của bản đồ rủi ro ngập lụt, hỗ trợ cơ quan quản lý thiên tai ban hành lệnh sơ tán chính xác, giảm thiểu thiệt hại về người và tài sản trước các biến cố thời tiết cực đoan.

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

  • Nghiên cứu sinh và học giả Toán Ứng dụng: Thừa hưởng cấu hình nghiệm Riemann giải tích chi tiết và tiêu chuẩn chấp nhận (MO) để phát triển các đề tài giải tích số chuyên sâu.
  • Kỹ sư R&D Thủy khí & Phần mềm: Tiếp cận các công thức số tường minh của lược đồ van Leer và FORCE để nhúng trực tiếp vào các mã nguồn mô phỏng CFD thương mại hoặc mã nguồn mở (như OpenFOAM, Telemac).
  • Cơ quan Khí tượng Thủy văn và Quản lý Đê điều: Tiếp nhận cơ sở khoa học để chuẩn hóa mô hình dự báo sóng thần và tính toán an toàn hồ đập thủy điện.

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?

Đóng góp lý thuyết độc đáo nhất là việc phân tích hoàn chỉnh cấu trúc nghiệm Riemann cho hệ Ripa phi bảo toàn trong chế độ cộng hưởng (resonant regime, cấu trúc $A3, B3$), nơi xuất hiện sự tương tác đồng vận tốc của 3 sóng sơ cấp tại $\sigma = 0$. Luận án mở rộng trực tiếp lý thuyết hệ hyperbolic phi bảo toàn của Pedro Ripa (1993) và khung phân tích Riemann của LeFloch - Thanh (2007) bằng cách chứng minh giải tích sự tồn tại duy nhất của nghiệm qua nghiệm của đa thức bậc ba $\psi(h) = 0$ và ngưỡng hình học $a_{\max}(U_0)$.

2. Sự đổi mới về phương pháp luận số học so với các nghiên cứu quốc tế trước đây?

So với lược đồ WENO của Han và Li (2014) và lược đồ thuận gió của Chertock và cộng sự (2013), luận án không sử dụng các toán tử xấp xỉ nhân tạo cho số hạng nguồn. Thay vào đó, tác giả đưa trực tiếp các trạng thái tĩnh trung gian giải tích $U_{j-1,+}^n, U_{j+1,-}^n$ vào các thông lượng số chuẩn (Lax-Friedrichs, Lax-Wendroff, FORCE), tạo nên lược đồ cân bằng tự nhiên (natively well-balanced) với chi phí tính toán thấp hơn và triệt tiêu hoàn toàn sai số trôi tại các điểm gián đoạn địa hình.

3. Phát hiện số trị gây bất ngờ nhất có dữ liệu thực nghiệm chứng minh là gì?

Phát hiện bất ngờ nhất là tại Kiểm định 8 (cộng hưởng sóng), lược đồ số kiểu Godunov và van Leer bắt gọn sóng tiếp xúc tĩnh $C_4$ ngay trên lưới thô (100 điểm) mà không sinh ra bất kỳ đỉnh dao động giả nào tại điểm gián đoạn $x = 0$, trong khi các lược đồ thể tích hữu hạn truyền thống đều bị mất hội tụ hoặc tạo sóng phản xạ giả với biên độ sai số lên tới $35%$.

4. Quy trình tái lập (Replication Protocol) có được cung cấp minh bạch không?

Luận án cung cấp tường minh 100% các hệ số toán học, công thức hàm giới hạn (minmod, superbee, van Leer), công thức thông lượng số FORCE, các điều kiện ban đầu $(h_L, u_L, \theta_L, a_L)$ và $(h_R, u_R, \theta_R, a_R)$, cùng giới hạn CFL cho toàn bộ 8 kiểm định 1D và 3 kiểm định 2D, cho phép cộng đồng khoa học lập trình tái lập chính xác mọi kết quả mô phỏng.

5. Lộ trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Lộ trình hướng tới 3 trụ cột: (1) Tích hợp bài toán biên khô - ướt ($h \to 0$) vào mô hình Ripa 2D; (2) Phát triển lược đồ Discontinuous Galerkin cấp cao thích nghi lưới cục bộ (AMR) trên miền hình học thực tế; (3) Thương mại hóa lõi phần mềm mô phỏng dòng chảy ngập lụt hỗ trợ ra quyết định khẩn cấp trong phòng chống thiên tai.

Kết luận

Luận án tiến sĩ của NCS Nguyễn Xuân Thanh đã đạt được 5 đóng góp học thuật cốt lõi:

  1. Thiết lập thành công và giải tích hóa tường minh toàn bộ cấu hình nghiệm Riemann cho hệ phương trình nước nông có xét biến thiên nhiệt độ (hệ Ripa) chứa số hạng nguồn phi bảo toàn.
  2. Xây dựng thành công lớp lược đồ cân bằng (Well-balanced schemes) dựa trên thông lượng tổ hợp lồi FORCE, triệt tiêu hoàn toàn dao động số tại biên đáy gián đoạn.
  3. Thiết kế thành công lược đồ kiểu Godunov chính xác và lược đồ cấp cao kiểu van Leer (MUSCL), chứng minh chặt chẽ tính bảo toàn cân bằng và tính luôn dương của chiều cao nước $h > 0$ và nhiệt độ $\theta > 0$.
  4. Phát triển hoàn chỉnh lược đồ cân bằng 2D cho hệ nước nông có đáy biến thiên bằng kỹ thuật tách chiều hai bước tối ưu hóa thời gian tính toán CPU.
  5. Kiểm chứng toàn diện lý thuyết qua hệ thống kiểm định số phong phú từ 1D đến 2D, chứng minh tốc độ hội tụ thực nghiệm tiệm cận bậc 2 lý thuyết và độ ổn định tuyệt đối trong miền cộng hưởng.

Công trình không chỉ đánh dấu bước tiến quan trọng trong lĩnh vực toán giải tích và giải tích số tại Việt Nam mà còn cung cấp công cụ tính toán thủy động lực học có giá trị ứng dụng thực tiễn cao trên quy mô quốc tế.