Tổng quan luận án

Biến đổi khí hậu toàn cầu kéo theo sự gia tăng của các hiện tượng thiên tai cực đoan như bão lũ, vỡ đập, sóng thần và sạt lở đất. Để mô phỏng, đánh giá và dự báo tác động của các hiện tượng dòng chảy này, việc nghiên cứu giải số các hệ phương trình đạo hàm riêng mô tả động lực học chất lỏng giữ vai trò nền tảng. Trong số đó, hệ phương trình nước nông (shallow water equations) và các dạng mở rộng có tính đến gradient nhiệt độ (hệ Ripa) hoặc địa hình đáy biến thiên là những mô hình toán học phản ánh trực tiếp các quy luật thủy động lực.

Về mặt giải tích, hệ hyperbolic các định luật cân bằng với nguồn phi bảo toàn là sự tổng quát hóa của hệ các định luật bảo toàn cổ điển (vốn chỉ bảo toàn khối lượng, động lượng hoặc năng lượng). Đối với các hệ hyperbolic tổng quát có vế phải chứa số hạng nguồn phi bảo toàn, việc tìm nghiệm giải tích chính xác cho bài toán Cauchy/Riemann gặp nhiều rào cản toán học phức tạp. Khi sử dụng các lược đồ sai phân chuẩn hoặc phương pháp số truyền thống, sự xuất hiện của số hạng nguồn phi bảo toàn thường gây ra dao động giả, sai số tăng cao hoặc làm mất tính ổn định khi làm mịn lưới tính toán. Do đó, yêu cầu xây dựng các phương pháp số xấp xỉ nghiệm yếu đảm bảo tính cân bằng (well-balanced), duy trì tính dương của các biến vật lý và đạt bậc hội tụ cao là khoảng trống nghiên cứu trọng tâm mà luận án hướng tới.

Luận án tiến sĩ toán học đề 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 tại Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia Thành phố Hồ Chí Minh (chuyên ngành Toán giải tích, mã số 9460102, khóa đào tạo 2020), dưới sự hướng dẫn khoa học của PGS.TS. Mai Đức Thành và TS. Nguyễn Thị Thu Vân.

Mục tiêu nghiên cứu của luận án được xác định cụ thể gồm 03 mục tiêu:

  1. Xây dựng các phương pháp số cho mô hình hệ các phương trình nước nông (hệ nước nông với nhiệt độ một chiều và hệ nước nông với đáy biến thiên trong không gian hai chiều).
  2. Chứng minh các lược đồ được xây dựng bắt được sóng tĩnh và thỏa mãn các tính chất định tính quan trọng (tính cân bằng, tính luôn dương của chiều cao nước $h$ và nhiệt độ $\theta$).
  3. Đưa ra các kiểm định số để chứng minh lược đồ bắt được sóng tĩnh, nghiệm số hội tụ về nghiệm chính xác, kể cả trong các trường hợp nghiệm có cấu trúc cộng hưởng và tương tác sóng phức tạp.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: Các hệ phương trình hyperbolic phi tuyến có số hạng nguồn phi bảo toàn, cụ thể là:
    • Hệ phương trình nước nông với gradient nhiệt độ ngang (hệ Ripa một chiều);
    • Hệ phương trình nước nông với địa hình đáy biến thiên trong không gian hai chiều.
  • Phạm vi nghiên cứu:
    • Không gian: Khảo sát bài toán Cauchy và bài toán Riemann trong không gian một chiều ($x \in \mathbb{R}$) và không gian hai chiều ($(x, y) \in \mathbb{R}^2$).
    • Miền trạng thái: Khảo sát đầy đủ các trạng thái dòng chảy bao gồm miền dưới tới hạn (subcritical: $Fr < 1$), miền trên tới hạn (supercritical: $Fr > 1$), đường cong tới hạn (critical: $Fr = 1$) và vùng xảy ra hiện tượng cộng hưởng (resonance).
    • Thời gian thực hiện và hoàn thành: Khóa đào tạo 2020–2024, bảo vệ năm 2024.

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

Hệ phương trình nước nông một chiều lần đầu tiên được thiết lập năm 1871 bởi Adhémar Jean Claude Barré de Saint-Venant, đặt nền móng cho các nghiên cứu thủy động lực học dòng chảy bề mặt tự do. Đến năm 1993, Pedro Ripa đã mở rộng hệ Saint-Venant bằng cách tích hợp sự biến thiên của nhiệt độ nước, tạo nên mô hình nước nông với gradient nhiệt độ ngang (thường gọi là hệ Ripa), ứng dụng phổ biến trong mô hình hóa dòng chảy đại dương và khí quyển. Bên cạnh đó, từ công trình khởi đầu về khí động lực học của Theodor Meyer năm 1908, nhiều mô hình liên quan như dòng chảy hai pha nén (Baer-Nunziato), dòng chảy chất lỏng trong ống có tiết diện biến thiên cũng lần lượt được phát triển.

Về phương diện giải tích nghiệm yếu và phương pháp số cho các hệ hyperbolic phi bảo toàn, luận án điểm lại các nhóm công trình tiêu biểu:

  • Nghiệm Riemann giải tích: Nhóm tác giả LeFloch và Thanh đã công bố các cấu hình nghiệm Riemann cho hệ nước nông với đáy gián đoạn cũng như dòng chảy chất lỏng trong ống có tiết diện gián đoạn; tác giả Thanh (2009) đã phân tích và thiết lập cấu hình nghiệm Riemann giải tích cho hệ nước nông với nhiệt độ. Nhóm Andrianov và Warnecke thiết lập nghiệm Riemann cho mô hình dòng chảy nén trong ống và dòng chảy hai pha Baer-Nunziato.
  • Lược đồ cân bằng (Well-balanced schemes): Greenberg và LeRoux đã đề xuất nguyên lý lược đồ số cân bằng cho số hạng nguồn trong phương trình hyperbolic dựa trên kỹ thuật sai phân hữu hạn (xử lý nguồn bị hấp thụ). Chinnayya, LeRoux và Seguin xây dựng lược đồ cân bằng thể tích hữu hạn cho hệ nước nông có đáy. Nhóm Thanh, Fazlul và Ismail phát triển lược đồ sai phân cân bằng cho hệ nước nông đáy tùy ý. Kröner và Thanh thiết lập lược đồ cân bằng cho dòng lưu chất trong ống tiết diện biến thiên.
  • Lược đồ kiểu Godunov và vấn đề cộng hưởng: LeFloch và Thanh công bố lược đồ cân bằng kiểu Godunov cho hệ nước nông đáy biến thiên trong miền cộng hưởng. Issacson và Temple khảo sát sự hội tụ của phương pháp Godunov cỡ $2 \times 2$ cho hệ cân bằng phi tuyến cộng hưởng tổng quát. Schwendeman cùng các cộng sự, và nhóm Ambrose nghiên cứu phương pháp Godunov cho mô hình dòng chảy hai pha nén.
  • Lược đồ cấp cao và mở rộng hai chiều: Quang và Thanh xây dựng lược đồ cân bằng kiểu van Leer cho hệ nước nông đáy biến thiên. Gallardo, Parés và Castro công bố lược đồ thể tích hữu hạn bậc cao cân bằng cho hệ nước nông đáy biến thiên có xét vùng khô. Han và Li thiết lập lược đồ sai phân hữu hạn WENO cân bằng cho mô hình Ripa. Chertock, Kurganov và Liu công bố lược đồ thuận gió (upwind) cho hệ nước nông với nhiệt độ trong không gian hai chiều.

Khoảng trống nghiên cứu được luận án giải quyết: Kế thừa các kết quả giải tích về bài toán Riemann của hệ Ripa, luận án hoàn thiện giải thuật tính toán tường minh nghiệm Riemann cục bộ làm nền tảng xây dựng lược đồ Godunov; thiết lập lược đồ cân bằng hấp thụ nguồn phi bảo toàn kết hợp các thông lượng số chuẩn (Lax-Friedrichs, Lax-Wendroff, FORCE); phát triển lược đồ cấp cao kiểu van Leer sử dụng các hàm giới hạn độ dốc; và mở rộng thuật toán cân bằng cho hệ nước nông hai chiều đáy biến thiên thông qua kỹ thuật tách chiều (dimensional splitting), đi kèm chứng minh toán học chặt chẽ về tính bảo toàn trạng thái dừng và tính dương của chiều cao nước.

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 phương trình đạo hàm riêng hyperbolic phi tuyến kết hợp giải thuật số trị hiện đại:

Khung lý thuyết và mô hình toán học

  1. Hệ phương trình Ripa 1D (nước nông với nhiệt độ): $$\begin{cases} \partial_t h + \partial_x (hu) = 0 \ \partial_t (hu) + \partial_x \left( hu^2 + \frac{1}{2}g h^2 \theta \right) = -gh\theta \partial_x a \ \partial_t (h\theta) + \partial_x (uh\theta) = 0 \end{cases}$$ Bổ sung phương trình dừng $\partial_t a = 0$ để đưa về dạng hệ chuẩn tắc phi bảo toàn: $\partial_t U + A(U)\partial_x U = 0$ với vector trạng thái $U = (h, hu, h\theta, a)^T$.
  2. Cấu trúc trị riêng và vector riêng: Ma trận hệ số $A(U)$ có 04 trị riêng thực: $$\lambda_1(U) = u - c, \quad \lambda_2(U) = u, \quad \lambda_3(U) = u + c, \quad \lambda_4(U) = 0$$ với vận tốc sóng âm $c = \sqrt{gh\theta}$. Trường đặc trưng 1 và 3 là thuần phi tuyến ($\nabla \lambda_i \cdot r_i \neq 0$); trường đặc trưng 2 và 4 là suy biến tuyến tính ($\nabla \lambda_i \cdot r_i = 0$).
  3. Số Froude tổng quát và phân vùng không gian trạng thái: Số Froude được xác định bởi: $$Fr(U) = \frac{|u|}{\sqrt{gh\theta}} = \frac{|u|}{c}$$
    • Miền $G_1 = {U: \lambda_4 < \lambda_1 < \lambda_2 < \lambda_3}$ và $G_4 = {U: \lambda_1 < \lambda_2 < \lambda_3 < \lambda_4}$ ứng với dòng chảy trên tới hạn ($Fr > 1$).
    • Miền $G_2 = {U: \lambda_1 < \lambda_4 < \lambda_2 < \lambda_3}$ và $G_3 = {U: \lambda_1 < \lambda_2 < \lambda_4 < \lambda_3}$ ứng với dòng chảy dưới tới hạn ($Fr < 1$).
    • Mặt tới hạn $C^{\pm} = {U: \lambda_{1,3}(U) = \lambda_4(U) = 0}$ ứng với $Fr = 1$.
  4. Hệ thức bước nhảy và nghiệm sóng dừng phi bảo toàn: Qua sóng tiếp xúc tĩnh (tiếp xúc-4, vận tốc $\sigma = 0$), nghiệm thỏa mãn hệ thức bước nhảy phi bảo toàn: $$[hu] = 0, \quad \left[ \frac{u^2}{2} + g\theta(h+a) \right] = 0, \quad [\theta] = 0$$ Từ đó dẫn đến phương trình bậc ba xác định chiều cao nước $h$: $$\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$$ Nghiệm sóng dừng được chọn duy nhất nhờ tiêu chuẩn chấp nhận đơn điệu (Monotonicity Criterion - MC) dựa trên giới hạn đáy cực đại $a_{\max}(U_0) = a_0 + h_0 + \frac{u_0^2}{2g\theta_0} - \frac{3}{2}\left(\frac{h_0^2 u_0^2}{g\theta_0}\right)^{1/3}$.

Phương pháp số thực tế

  1. Lược đồ cân bằng cho định luật bảo toàn có nguồn (1D): Xấp xỉ phương trình $\partial_t U + \partial_x F(U) = S(U)\partial_x H(x)$ theo công thức: $$U_j^{n+1} = U_j^n - \frac{\Delta t}{\Delta x} \left[ f(U_j^n, U_{j+1,-}^n) - f(U_{j-1,+}^n, U_j^n) \right]$$ trong đó $U_{j+1,-}^n$ và $U_{j-1,+}^n$ là các trạng thái tĩnh trung gian bên trái và bên phải giải từ sóng dừng để nguồn phi bảo toàn bị triệt tiêu hoàn toàn tại trạng thái cân bằng.
  2. Các hàm thông lượng số chuẩn:
    • Thông lượng Lax-Friedrichs ($f^{LF}$): $$f^{LF}(U_1, U_2) = \frac{1}{2}\left[F(U_1) + F(U_2)\right] - \frac{\Delta x}{2\Delta t}(U_2 - U_1)$$
    • Thông lượng Lax-Wendroff ($f^{LW}$): $$f^{LW}(U_1, U_2) = \frac{1}{2}\left[F(U_1) + F(U_2)\right] - \frac{\Delta t}{2\Delta x} A(U_1, U_2)\left[F(U_2) - F(U_1)\right]$$
    • Thông lượng FORCE ($f^{FORCE}$): $$f^{FORCE}(U_1, U_2) = \frac{1}{2}f^{LF}(U_1, U_2) + \frac{1}{2}f^{LW}(U_1, U_2)$$
  3. Phương pháp Godunov: Sử dụng nghiệm chính xác của bài toán Riemann $U_{exact}(x/t; U_L, U_R)$ để xác định thông lượng số: $$U_j^{n+1} = U_j^n - \frac{\Delta t}{\Delta x} \left[ F(U_{exact}(0; U_j^n, U_{j+1}^n)) - F(U_{exact}(0; U_{j-1}^n, U_j^n)) \right]$$ thỏa mãn điều kiện CFL: $\frac{\Delta t}{\Delta x}\max_k |\lambda_k(U_j^n)| \le \frac{1}{2}$.
  4. Phương pháp cấp cao kiểu van Leer: Tái dựng tuyến tính các trạng thái biên ô lưới tại $t_{n+1/2}$ sử dụng các hàm giới hạn độ dốc $\phi(p)$ nhằm loại bỏ hiện tượng giả dao động quanh mặt gián đoạn:
    • Hàm superbee: $\phi(p) = \max(0, \min(1, 2p), \min(p, 2))$
    • Hàm van Leer: $\phi(p) = \frac{p + |p|}{1 + |p|}$
    • Hàm minmod: $\phi(p) = \max(0, \min(p, 1))$
  5. Phương pháp tách chiều (Dimensional Splitting) cho hệ 2D: Tách hệ nước nông 2 chiều: $$\partial_t U + \partial_x F(U) + \partial_y H(U) = S_x(U)\partial_x a + S_y(U)\partial_y a$$ thành hai bài toán 1D kế tiếp trên lưới đều $\Delta x \times \Delta y$:
    • Bước 1 (theo phương $x$): Giải $\partial_t U + \partial_x F(U) = S_x(U)\partial_x a$ trong khoảng $\Delta t$ để tìm trạng thái trung gian $U^{n+1/2}$.
    • Bước 2 (theo phương $y$): Giải $\partial_t U + \partial_y H(U) = S_y(U)\partial_y a$ với dữ liệu đầu $U^{n+1/2}$ để thu được $U^{n+1}$.

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

Chương 1: Mở đầu

Chương này giới thiệu tổng quan bối cảnh thực tiễn của biến đổi khí hậu và sự cần thiết của các mô hình toán học dự báo thiên tai (lũ lụt, vỡ đập, sóng thần). Tác giả trình bày khái quát về hệ hyperbolic các định luật cân bằng với nguồn phi bảo toàn, làm rõ khó khăn trong việc tìm nghiệm giải tích và xấp xỉ số. Từ đó, xác định rõ mục đích nghiên cứu, 03 mục tiêu cụ thể, đối tượng khảo sát (hệ Ripa 1D và hệ nước nông 2D đáy biến thiên) cùng bố cục 5 chương của luận án.

Chương 2: Tổng quan

Chương 2 phân tích, tổng hợp và đánh giá chi tiết các công trình nghiên cứu trong và ngoài nước liên quan trực tiếp đến đề tài. Nội dung tập trung vào:

  • Lịch sử hình thành hệ phương trình nước nông Saint-Venant (1871) và hệ Ripa (1993).
  • Các kết quả về cấu hình nghiệm giải tích Riemann cho hệ phi bảo toàn của LeFloch, Thanh, Andrianov, Warnecke.
  • Các phương pháp số cân bằng cho định luật bảo toàn có số hạng nguồn của Greenberg–LeRoux, Chinnayya–LeRoux–Seguin, Thanh–Fazlul–Ismail, Kröner–Thanh.
  • Lược đồ kiểu Godunov trong miền cộng hưởng của LeFloch–Thanh, Issacson–Temple.
  • Các lược đồ cấp cao và mở rộng 2D của Quang–Thanh, Gallardo–Parés–Castro, Han–Li, Chertock–Kurganov–Liu. Chương này định vị rõ đóng góp của luận án là bước phát triển tiếp nối các công trình trên nhằm đưa ra các lược đồ cân bằng, Godunov, van Leer và tách chiều 2D có khả năng bắt sóng dừng chính xác.

Chương 3: Phương pháp nghiên cứu

Chương 3 trình bày các cơ sở phương pháp luận toán học và công thức tính toán làm nền tảng cho việc xây dựng các lược đồ số ở Chương 4:

  • Mục 3.1: Thiết lập phương pháp số cân bằng cho một định luật bảo toàn vô hướng có nguồn, cơ chế đưa các trạng thái tĩnh trung gian vào thông lượng số để triệt tiêu nguồn vi phân tại biên ô.
  • Mục 3.2: Trình bày thuật toán Godunov cho định luật bảo toàn dựa trên nghiệm giải tích bài toán Riemann và điều kiện ổn định CFL.
  • Mục 3.3: Xây dựng phương pháp cấp cao sử dụng các hàm giới hạn góc dốc (van Leer, superbee, minmod).
  • Mục 3.4: Xây dựng phương pháp tách chiều hai bước cho định luật bảo toàn trong không gian hai chiều, thiết lập công thức cập nhật theo phương $x$ và phương $y$.

Chương 4: Kết quả nghiên cứu và phân tích, đánh giá, thảo luận

Đây là chương trọng tâm của luận án, được chia thành 3 phần nghiên cứu chính:

1. Xây dựng một số lược đồ số cho hệ phương trình nước nông với nhiệt độ (Mục 4.1)

  • Tính chất giải tích: Khảo sát ma trận Jacobi, 04 trường đặc trưng, tính chất sóng sốc Lax ($S_1, S_3$), sóng giãn ($R_1, R_3$), sóng tiếp xúc vật chất ($W_2$), sóng tiếp xúc tĩnh ($W_4$) và thiết lập phương trình bậc ba $\psi(h) = 0$.
  • Cấu trúc nghiệm Riemann giải tích: Xây dựng đầy đủ các cấu trúc nghiệm Riemann chính xác:
    • Trường hợp A ($U_L \in G_1 \cup C^+$): Gồm Cấu trúc A1 (trạng thái trái/phải đều supercritical), Cấu trúc A2 (trái supercritical, phải subcritical), Cấu trúc A3 (hiện tượng cộng hưởng chứa 3 sóng lan truyền cùng vận tốc 0).
    • Trường hợp B ($U_L \in G_2$): Gồm Cấu trúc B1, B2, B3 ứng với các chuyển tiếp trạng thái khác nhau.
  • Xây dựng lược đồ và chứng minh toán học:
    • Thiết lập lược đồ cân bằng và lược đồ Godunov.
    • Chứng minh định lý: Lược đồ bắt chính xác các sóng tĩnh và duy trì tính luôn dương của chiều cao nước $h$ và nhiệt độ $\theta$.
  • Kiểm định số: Thực hiện 8 bài kiểm định số (Kiểm định 1 đến 8) đánh giá sai số chuẩn $L^1$ và bậc hội tụ trên các cỡ lưới khác nhau (từ 100 đến 1600 điểm lưới), bao gồm cả bài toán tương tác sóng và cộng hưởng.
Cấu trúc nghiệm Riemann Miền trạng thái ban đầu Thành phần cấu trúc sóng Đặc điểm tính toán
Cấu trúc A1 $U_L \in G_1, U_R \in G_1$ $W_4(U_L, U_1) + W_1(U_1, U_2) + W_2(U_2, U_3) + W_3(U_3, U_R)$ Cả hai trạng thái thuộc miền supercritical
Cấu trúc A2 $U_L \in G_1, U_R \in G_2$ $W_1(U_L, U_1) + W_4(U_1, U_2) + W_2(U_2, U_3) + W_3(U_3, U_R)$ Chuyển tiếp từ supercritical sang subcritical
Cấu trúc A3 $U_L \in G_1, U_R \approx C^0$ $W_4(U_L, U_1) + S_1^0(U_1, U_2) + W_4(U_2, U_3) + W_2(U_3, U_4) + W_3(U_4, U_R)$ Hiện tượng cộng hưởng (3 sóng vận tốc 0)
Cấu trúc B1 $U_L \in G_2, U_R \in G_1$ $R_1(U_L, U_1) + W_4(U_1, U_2) + W_1(U_2, U_3) + W_2(U_3, U_4) + W_3(U_4, U_R)$ Chuyển tiếp từ subcritical sang supercritical

2. Lược đồ số cấp cao cho hệ phương trình nước nông với nhiệt độ (Mục 4.2)

  • Thiết lập lược đồ cân bằng cấp cao kiểu van Leer kết hợp các hàm giới hạn slope limiter (van Leer, superbee, minmod).
  • Chứng minh tính bảo toàn trạng thái dừng, tính dương của $h$ và $\theta$.
  • Tiến hành các kiểm định số (Kiểm định 1 đến 8) với dữ liệu đầu trơn và dữ liệu Riemann gián đoạn; phân tích sai số $|E|_{L^1}$ và kiểm chứng bậc hội tụ cao vượt trội so với lược đồ cấp một và lược đồ cổ điển.

3. Lược đồ số cho hệ phương trình nước nông trong không gian hai chiều (Mục 4.3)

  • Thiết lập hệ 2D nước nông đáy biến thiên: $$\begin{cases} \partial_t h + \partial_x (hu) + \partial_y (hv) = 0 \ \partial_t (hu) + \partial_x \left( hu^2 + \frac{1}{2}gh^2 \right) + \partial_y (huv) = -gh\partial_x a \ \partial_t (hv) + \partial_x (huv) + \partial_y \left( hv^2 + \frac{1}{2}gh^2 \right) = -gh\partial_y a \end{cases}$$
  • Xây dựng lược đồ cân bằng hai bước sử dụng phương pháp tách chiều kết hợp các thông lượng số Lax-Friedrichs, Lax-Wendroff và FORCE.
  • Thực hiện các kiểm định số 2D (Kiểm định 1, 2, 3) trên các lưới $40 \times 40, 80 \times 80, 100 \times 100, 200 \times 200$:
    • Kiểm tra khả năng bắt nghiệm gián đoạn tĩnh chính xác;
    • Kiểm tra bài toán planar Riemann;
    • Đo đạc thời gian CPU và đánh giá bậc hội tụ số trị giữa các dạng thông lượng.

Chương 5: Kết luận và kiến nghị

Chương tổng kết toàn bộ các kết quả đạt được của luận án và đề xuất các định hướng nghiên cứu tiếp theo về việc nâng cao bậc chính xác và tối ưu hóa thời gian tính toán cho các hệ hyperbolic phi bảo toàn phức tạp.

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

Các đóng góp mới về mặt khoa học của luận án

Theo công bố chính thức trong văn bản luận án, các kết quả mới bao gồm:

  1. Xây dựng thành công lược đồ cân bằng cho hệ các phương trình nước nông có gradient nhiệt độ ngang (hệ Ripa).
  2. Xây dựng lược đồ số kiểu Godunov cho hệ các phương trình nước nông với nhiệt độ dựa trên giải thuật tính toán tường minh nghiệm giải tích Riemann.
  3. Xây dựng lược đồ số cấp cao kiểu van Leer cho hệ các phương trình nước nông với nhiệt độ, đạt bậc hội tụ cao hơn so với lược đồ Godunov và lược đồ cổ điển.
  4. Xây dựng lược đồ cân bằng cho hệ các phương trình nước nông với đáy biến thiên trong không gian hai chiều bằng phương pháp tách chiều.

Toàn bộ các kết quả chính của luận án đã được công bố trong 04 bài báo khoa học chuyên ngành (ký hiệu từ (P1) đến (P4)) và được báo cáo tại 02 hội nghị quốc tế: Hội nghị ISAS 2019 (International Symposium on Applied Science 2019) và Hội nghị ISAS 2020.

Khả năng ứng dụng và kiến nghị

  • Khả năng ứng dụng thực tiễn: Mô hình hệ phương trình nước nông mô tả các hiện tượng thủy văn thực tế liên quan đến dòng chảy thiên tai như lũ lụt, sóng thần, vỡ đập. Các lược đồ số cân bằng được xây dựng trong luận án cung cấp công cụ toán học hỗ trợ việc tính toán, mô phỏng và dự báo tác động của ngập lụt hay sự cố vỡ đập.
  • Kiến nghị: Tiếp tục nghiên cứu, mở rộng các lược đồ số cho hệ hyperbolic với nguồn phi bảo toàn nhằm đạt độ chính xác cao hơn và tốc độ tính toán nhanh hơn, hướng đến việc tích hợp vào các công cụ dự báo chuyên dụng.

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

  • Hạn chế: Luận án tập trung khảo sát trong khuôn khổ mô hình nước nông 1D có nhiệt độ và mô hình nước nông 2D đáy biến thiên; các kiểm định số chủ yếu thực hiện trên các miền hình học tiêu chuẩn và lưới tính toán có cấu trúc đều; chưa tích hợp điều kiện biên phức tạp cho địa hình thực tế quy mô lớn.
  • Hướng nghiên cứu tiếp theo:
    1. Phát triển các lược đồ số bậc cao hơn (như WENO, Discontinuous Galerkin) cho hệ hyperbolic với nguồn phi bảo toàn.
    2. Tối ưu hóa thuật toán song song nhằm giảm thiểu thời gian CPU khi tính toán trên lưới kích thước lớn trong không gian hai và ba chiều.
    3. Mở rộng phương pháp cho các mô hình đa pha, dòng chảy nén trong ống hoặc dòng chất lỏng có xét ma sát đáy phức tạp.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên khảo giá trị cho:

  • 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: Cung cấp phương pháp chi tiết về giải tích cấu trúc sóng Riemann, kỹ thuật xử lý số hạng nguồn phi bảo toàn và chứng minh tính ổn định định tính của lược đồ số.
  • Giảng viên và nhà nghiên cứu trong lĩnh vực Giải tích số, Cơ học thủy khí tính toán (CFD): Tham khảo giải thuật xây dựng lược đồ cân bằng (well-balanced), lược đồ cấp cao van Leer và kỹ thuật tách chiều 2D kết hợp thông lượng FORCE/Lax-Wendroff.
  • Kỹ sư và cán bộ nghiên cứu mô hình hóa thủy văn: Tham khảo nền tảng thuật toán mô phỏng các hiện tượng dòng chảy mặt tự do, dòng triều có biến thiên nhiệt độ và sóng lũ qua đáy biến thiên.

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

1. Nguồn phi bảo toàn trong hệ phương trình nước nông với nhiệt độ xuất phát từ đâu và gây khó khăn gì cho tính toán số?

Nguồn phi bảo toàn trong hệ Ripa xuất hiện ở vế phải của phương trình động lượng dưới dạng số hạng $-gh\theta \partial_x a$ (tích giữa gradient độ cao đáy $\partial_x a$ và các biến phụ thuộc $h, \theta$). Về mặt toán học, số hạng này không thể đưa về dạng đạo hàm toàn phần theo tọa độ không gian. Trong tính toán số, việc không bảo toàn này dẫn đến sự mất cân bằng giữa thông lượng và nguồn, gây ra các dao động giả phi vật lý hoặc làm sai lệch nghiệm khi kích thước lưới tính toán tiến dần về 0.

2. Tiêu chuẩn nào được sử dụng để xác định nghiệm duy nhất cho bài toán Riemann của hệ Ripa?

Luận án sử dụng Tiêu chuẩn chấp nhận đơn điệu (Monotonicity Criterion - MC). Theo đó, dọc theo bất kỳ đường cong sóng tĩnh $W_4(U_0)$ nào, độ cao đáy $a$ phải là một hàm đơn điệu theo chiều cao nước $h$, và tổng độ lệch chiều cao đáy của nghiệm Riemann không được vượt quá bước nhảy đáy ban đầu $|a_L - a_R|$. Giới hạn này gắn liền với điều kiện $a \le a_{\max}(U_0)$ để phương trình bậc ba $\psi(h) = 0$ luôn đảm bảo tồn tại nghiệm thực phân biệt $h_1 < h_2$.

3. Cơ chế triệt tiêu nguồn vi phân trong lược đồ cân bằng của luận án hoạt động như thế nào?

Lược đồ cân bằng thay thế các trạng thái biên ô chuẩn $U_{j-1}, U_{j+1}$ bằng các trạng thái tĩnh trung gian $U_{j-1,+}^n, U_{j+1,-}^n$. Các trạng thái này được tính toán trực tiếp từ nghiệm của sóng dừng phi bảo toàn nối qua các điểm biên ô lưới $x_{j-1/2}$ và $x_{j+1/2}$. Khi dòng chảy ở trạng thái dừng (cân bằng tĩnh), các bước nhảy thông lượng sẽ triệt tiêu chính xác với nguồn vi phân, giúp lược đồ duy trì trạng thái dừng theo thời gian mà không sinh ra sai số số trị.

4. Hiện tượng cộng hưởng (resonance) trong hệ Ripa là gì và cấu trúc nghiệm Riemann tương ứng được mô tả ra sao?

Hiện tượng cộng hưởng xuất hiện khi trạng thái biên nằm sát đường cong tới hạn ($Fr \approx 1$), dẫn đến việc nghiệm giải tích Riemann xuất hiện 3 sóng thành phần khác nhau cùng lan truyền với vận tốc bằng 0 (Cấu trúc A3). Trong trường hợp này, cấu trúc nghiệm bao gồm: sóng tĩnh $W_4$ nối từ $U_L$ đến $U_1$, tiếp nối bởi sóng sốc $S_1^0$ vận tốc 0 đến $U_2$, sóng tĩnh $W_4$ thứ hai đến $U_3$, sóng tiếp xúc vật chất $W_2$ đến $U_4$, và kết thúc bằng sóng-3 ($W_3$) nối đến trạng thái phải $U_R$.

5. Phương pháp tách chiều (dimensional splitting) giải quyết hệ nước nông hai chiều như thế nào?

Phương pháp tách chiều chia việc cập nhật trạng thái $U_{j,k}^{n+1}$ trên lưới $2D$ thành hai bước 1D liên tiếp: Bước 1 cập nhật theo phương $x$ sử dụng lược đồ cân bằng 1D với tỷ số bước $\lambda_x = \Delta t/\Delta x$ để thu được $U_{j,k}^{n+1/2}$; Bước 2 sử dụng $U_{j,k}^{n+1/2}$ làm điều kiện đầu để cập nhật theo phương $y$ với tỷ số bước $\lambda_y = \Delta t/\Delta y$, qua đó thu được $U_{j,k}^{n+1}$. Mỗi bước đều tích hợp các thông lượng số chuẩn (Lax-Friedrichs, Lax-Wendroff, FORCE) giúp duy trì tính cân bằng trên toàn miền 2 chiều.

Kết luận

Luận án tiến sĩ của nghiên cứu sinh Nguyễn Xuân Thanh đã giải quyết trọn vẹn bài toán xây dựng và kiểm định các phương pháp số cho hệ phương trình hyperbolic có nguồn phi bảo toàn, tập trung vào hệ nước nông với gradient nhiệt độ và hệ nước nông hai chiều đáy biến thiên. Tác giả đã thiết lập thành công hệ thống lược đồ cân bằng, lược đồ Godunov dựa trên nghiệm Riemann giải tích, lược đồ cấp cao van Leer và lược đồ tách chiều 2D, đồng thời chứng minh đầy đủ tính cân bằng sóng tĩnh cũng như tính luôn dương của các biến vật lý. Kết quả nghiên cứu đóng góp cơ sở lý thuyết giải tích số quan trọng và cung cấp giải thuật tin cậy phục vụ mô phỏng, tính toán các bài toán thủy động lực học dòng chảy trong thực tiễn.