Tổng quan về luận án

Bối cảnh khoa học hiện đại đang đối mặt với những thách thức cấp bách từ biến đổi khí hậu, như lũ lụt, sóng thần, và các hiện tượng thời tiết cực đoan. Việc dự báo chính xác các hiện tượng dòng chảy này là vô cùng quan trọng cho công tác phòng chống thiên tai và phát triển bền vững. Trong bối cảnh đó, các mô hình toán học như hệ phương trình nước nông (Shallow Water Equations - SWEs) đã trở thành công cụ không thể thiếu. Luận án "Xây dựng phương pháp số cho một số hệ hyperbolic với nguồn phi bảo toàn" của Nguyễn Xuân Thanh là một đóng góp tiên phong trong lĩnh vực này, tập trung vào việc phát triển các lược đồ số hiệu quả và chính xác cho các hệ hyperbolic phức tạp, đặc biệt là SWEs với gradient nhiệt độ và đáy biến thiên. Tính tiên phong của nghiên cứu thể hiện ở việc giải quyết những khó khăn cố hữu khi xây dựng lược đồ số cho các hệ có nguồn phi bảo toàn, nơi các phương pháp chuẩn thường gặp phải các vấn đề về nhiễu và sai số lớn.

Research gap cụ thể mà luận án này giải quyết là sự thiếu hụt các phương pháp số mạnh mẽ và đáng tin cậy để xấp xỉ nghiệm yếu của các bài toán Riemann cho hệ hyperbolic tổng quát với nguồn phi bảo toàn, đặc biệt khi nghiệm chính xác khó tìm. Như đã được nêu rõ trong luận án: "Đối với hệ hyperbolic các định luật cân bằng tổng quát, không phải lúc nào chúng ta cũng thuận lợi trong việc tìm được nghiệm chính xác của bài toán Riemann cho hệ này. Do đó, vấn đề được đặt ra là vẫn xây dựng các phương pháp số để xấp xỉ nghiệm yếu của bài toán Riemann cho hệ hyperbolic các định luật cân bằng tổng quát." (Chương 1). Nền tảng này nhấn mạnh nhu cầu về các lược đồ số có khả năng bắt được sóng tĩnh một cách chính xác (well-balanced property) và duy trì các tính chất vật lý quan trọng (như tính dương của chiều cao nước và nhiệt độ), điều mà nhiều lược đồ cổ điển không đảm bảo. Các nghiên cứu trước đây như của LeFloch và Thanh [39, 38] đã cung cấp các cấu hình nghiệm Riemann chính xác cho các hệ có đáy gián đoạn, hoặc của Thanh [29] cho SWEs với nhiệt độ, nhưng việc tích hợp các nghiệm này vào các lược đồ số tổng quát, có bậc cao và đảm bảo tính cân bằng vẫn còn là một thách thức lớn.

Luận án được thúc đẩy bởi các câu hỏi nghiên cứu và giả thuyết sau:

  1. RQ1: Làm thế nào để xây dựng các lược đồ số cân bằng (well-balanced numerical schemes) cho hệ phương trình nước nông có gradient nhiệt độ ngang, đảm bảo bắt được sóng tĩnh và duy trì tính dương của các biến vật lý?
    • H1.1: Các lược đồ số kiểu Godunov và van Leer có thể được mở rộng và điều chỉnh để đạt được tính cân bằng cho SWEs với nhiệt độ.
    • H1.2: Việc sử dụng các nghiệm Riemann chính xác làm nền tảng cho việc xây dựng thông lượng số sẽ cải thiện độ chính xác và tính ổn định của các lược đồ.
  2. RQ2: Làm thế nào để phát triển các lược đồ số có bậc hội tụ cao hơn cho SWEs với nhiệt độ, vượt trội so với các phương pháp bậc nhất hiện có?
    • H2.1: Các lược đồ số kiểu van Leer, với sự tích hợp của các hàm giới hạn (limiters) phù hợp, sẽ cung cấp bậc hội tụ cao hơn trong khi vẫn duy trì tính cân bằng và bảo toàn tính dương.
  3. RQ3: Có thể xây dựng các lược đồ số cân bằng hiệu quả cho hệ phương trình nước nông hai chiều với đáy biến thiên, ứng dụng trong các bài toán thực tế phức tạp hơn không?
    • H3.1: Phương pháp tách chiều (operator splitting) có thể được kết hợp với các lược đồ cân bằng một chiều để giải quyết các hệ hai chiều một cách hiệu quả.

Khung lý thuyết của luận án được xây dựng dựa trên lý thuyết về hệ hyperbolic của các định luật bảo toàn và các định luật cân bằng. Các lý thuyết cốt lõi bao gồm:

  • Hệ phương trình nước nông (Shallow Water Equations - SWEs), ban đầu được đề xuất bởi Adhémar Jean Claude Barré de Saint-Venant năm 1871 [8], là mô hình cơ bản cho dòng chảy mặt tự do.
  • Hệ Ripa (Ripa System), được Pedro Ripa mở rộng năm 1993 [38, 41, 42], bổ sung thêm ảnh hưởng của gradient nhiệt độ vào SWEs, mô tả dòng chảy đại dương và khí quyển.
  • Lý thuyết bài toán Riemann (Riemann Problem Theory), cung cấp nghiệm chính xác cho các bài toán Cauchy với dữ liệu đầu từng khúc hằng số, đóng vai trò nền tảng cho các lược đồ Godunov và van Leer.
  • Điều kiện Rankine-Hugoniot và Bất đẳng thức sốc Lax để phân tích tính chất của sóng gián đoạn.
  • Lý thuyết về lược đồ số cân bằng (Well-Balanced Schemes), nhằm mục đích bắt chính xác các trạng thái cân bằng (sóng tĩnh) mà không tạo ra nhiễu số.

Đóng góp đột phá của luận án bao gồm việc phát triển các lược đồ số mới, vượt trội về độ chính xác và tính ổn định. Cụ thể:

  1. Xây dựng lược đồ số cân bằng cho SWEs với nhiệt độ: Luận án đã thành công trong việc tạo ra lược đồ cân bằng kiểu Godunov và lược đồ cân bằng tổng quát hơn cho hệ Ripa (SWEs với nhiệt độ). Các lược đồ này được chứng minh là bắt chính xác các sóng tĩnh và bảo toàn tính dương của chiều cao nước $h$ và nhiệt độ $\theta$. Các kết quả này đã được công bố trong các bài báo khoa học (P1) và (P2), với các kiểm định số cho thấy sự hội tụ của nghiệm xấp xỉ tới nghiệm chính xác trên các cỡ lưới khác nhau, ví dụ, với 250 điểm lưới, lược đồ đã xấp xỉ chính xác sóng tĩnh (Hình 4.7, 4.9).
  2. Phát triển lược đồ số van Leer bậc cao cho SWEs với nhiệt độ: Đây là một đóng góp quan trọng (P3), cung cấp các phương pháp số có bậc hội tụ cao hơn so với các lược đồ Godunov bậc nhất. Các kiểm định số (ví dụ, với 100 điểm lưới hay 1600 điểm lưới trong Hình 4.37, 4.42) đã chứng minh khả năng bắt trạng thái cân bằng và tính dương của chiều cao nước, cùng với bậc hội tụ cao hơn, mang lại độ chính xác vượt trội.
  3. Thiết lập lược đồ cân bằng cho SWEs 2D với đáy biến thiên: Ứng dụng phương pháp tách chiều, luận án đã xây dựng lược đồ cân bằng cho hệ phương trình nước nông hai chiều có đáy biến thiên (P4). Phương pháp này giải quyết được độ phức tạp của bài toán hai chiều trong khi vẫn duy trì tính cân bằng và bảo toàn tính dương, được kiểm chứng thông qua các kiểm định số với các lưới như 200x200 hoặc 100x100 điểm lưới (Hình 4.51-4.59), cho thấy khả năng xấp xỉ chính xác các nghiệm gián đoạn tĩnh.

Phạm vi nghiên cứu (Scope) của luận án bao gồm các hệ hyperbolic một chiều và hai chiều với nguồn phi bảo toàn, tập trung vào mô hình SWEs với nhiệt độ và SWEs với đáy biến thiên. Các kiểm định số được thực hiện trên các dải điểm lưới đa dạng từ 40, 100, 250, 400, 800, đến 1600 điểm lưới cho trường hợp 1D và 40x40, 80x80, 100x100, 200x200 điểm lưới cho trường hợp 2D, đảm bảo đánh giá toàn diện về tính hội tụ và độ chính xác. Khung thời gian nghiên cứu bao gồm việc phân tích các giải pháp Riemann địa phương và mô phỏng sự tiến triển của các nghiệm theo thời gian (t=0 đến các thời điểm cụ thể như t=0.15 trong Hình 4.27).

Ý nghĩa (Significance) của luận án là rất lớn. Về mặt học thuật, nó mở rộng ranh giới của phương pháp số cho các hệ hyperbolic, đặc biệt là các hệ có nguồn phi bảo toàn, một lĩnh vực đầy thách thức. Về mặt thực tiễn, các lược đồ số được phát triển trong luận án là "rất hữu ích trong việc dự báo tác động của lũ lụt hay vỡ đập" (Tóm tắt luận án), cung cấp công cụ mạnh mẽ hơn cho các nhà khoa học, kỹ sư và hoạch định chính sách trong việc dự báo và quản lý thiên tai.

Literature Review và Positioning

Nghiên cứu về các hệ hyperbolic với nguồn phi bảo toàn là một lĩnh vực quan trọng trong toán học ứng dụng, mô tả nhiều hiện tượng vật lý phức tạp như dòng chảy chất lỏng, khí động lực học và các hiện tượng thiên tai. Luận án này tổng hợp các luồng nghiên cứu chính và định vị bản thân một cách rõ ràng trong bối cảnh học thuật hiện đại.

Tổng hợp các luồng nghiên cứu chính:

  1. Hệ phương trình nước nông (SWEs) và các mở rộng: Hệ SWEs một chiều được giới thiệu bởi Adhémar Jean Claude Barré de Saint-Venant vào năm 1871 [8]. Sau đó, Pedro Ripa đã mở rộng thành hệ SWEs với nhiệt độ vào năm 1993 [38, 41, 42], được biết đến với tên gọi Hệ Ripa, mô tả dòng chảy đại dương và khí quyển. Các nghiên cứu liên quan đến SWEs đã được ứng dụng rộng rãi trong mô hình hóa lũ lụt, sóng thần [8, 21, 22, 24, 5, 20, 34, 47, 55, 54].
  2. Nghiệm yếu và bài toán Riemann: Việc tìm nghiệm yếu cho các hệ hyperbolic phi tuyến là một thách thức lớn. Các công trình tiên phong trong việc cung cấp cấu hình nghiệm chính xác của bài toán Riemann cho các hệ có nguồn phi bảo toàn bao gồm LeFloch và Thanh [39, 38] cho SWEs với đáy gián đoạn, Thanh [29] cho SWEs với nhiệt độ, và Audrianov-Warnecke [2, 33] cho mô hình dòng chảy nén trong ống và dòng chảy hai pha Baor-Nunziato. Những kết quả này là cực kỳ quan trọng, làm nền tảng cho việc xây dựng các lược đồ số kiểu Godunov và van Leer.
  3. Lược đồ số cân bằng (Well-balanced schemes): Lược đồ cân bằng là một phương pháp số có khả năng bắt chính xác các trạng thái cân bằng tĩnh. Greenberg và Leroux [3] đã công bố lược đồ cân bằng sử dụng phương pháp sai phân hữu hạn. Chinnayya-LeRoux-Sogin [8] đã phát triển lược đồ cân bằng cho SWEs với đáy bằng phương pháp thể tích hữu hạn. Các công trình của Thanh-Fazlul-Ismail [31] và Kröner-Thanh [9] tiếp tục mở rộng lược đồ cân bằng cho SWEs với đáy tùy ý và dòng chảy trong ống tiết diện biến thiên. LeFloch-Thanh [37] cũng công bố lược đồ cân bằng kiểu Godunov cho SWEs với đáy biến thiên trong miền cộng hưởng.
  4. Lược đồ Godunov và các phương pháp bậc cao: Lược đồ Godunov là một phương pháp bậc nhất phổ biến dựa trên giải bài toán Riemann. Issacson-Temple [12] đã trình bày sự hội tụ của phương pháp Godunov cho các hệ luật cân bằng phi tuyến. Các lược đồ Godunov cũng được nghiên cứu cho mô hình dòng chảy hai pha nén [10, 1]. Để cải thiện độ chính xác, các lược đồ bậc cao đã được phát triển, chẳng hạn như lược đồ kiểu van Leer cân bằng của Quang-Thanh [7] cho SWEs với đáy biến thiên, và lược đồ thể tích hữu hạn bậc cao của Gallardo-Parés-Castro [27]. Han-Li [37] đã trình bày lược đồ WENO cân bằng cho mô hình Ripa. Cherbook-Kurganov-Liu [2] đã công bố lược đồ thuận gió cho SWEs với nhiệt độ trong không gian hai chiều, là một trong những nghiên cứu gần nhất mà luận án này tiếp nối.

Mâu thuẫn/Tranh luận và Định vị trong Literature: Một mâu thuẫn chính trong việc phát triển phương pháp số cho các hệ hyperbolic có nguồn phi bảo toàn là sự đánh đổi giữa độ chính xác (bậc cao) và khả năng duy trì các tính chất vật lý quan trọng (tính cân bằng, bảo toàn tính dương). Nhiều lược đồ số cổ điển (như lược đồ cổ điển [3.3] được nhắc đến trong Chương 3) có thể xấp xỉ nghiệm nhưng "chưa đảm bảo được tính cân bằng theo nghĩa lược đồ có thể bắt được sóng tĩnh theo thời gian." Các phương pháp bậc nhất như Godunov thường ổn định nhưng có độ chính xác thấp hơn, trong khi các phương pháp bậc cao có thể tạo ra nhiễu số nếu không được thiết kế cẩn thận để xử lý các số hạng nguồn phi bảo toàn.

  • Opposing View 1: Độ chính xác so với tính ổn định/cân bằng: Các phương pháp bậc cao như WENO thường được ưa chuộng vì độ chính xác nhưng việc tích hợp tính cân bằng cho các hệ với nguồn phi bảo toàn là phức tạp. Luận án này giải quyết vấn đề này bằng cách thiết kế các lược đồ van Leer cân bằng, tích hợp các kỹ thuật giải Riemann và hàm giới hạn để đạt được cả hai mục tiêu.
  • Opposing View 2: Phương pháp sai phân hữu hạn so với thể tích hữu hạn: Cả hai phương pháp đều được sử dụng rộng rãi. Mặc dù luận án chủ yếu sử dụng phương pháp thể tích hữu hạn và sai phân chuẩn, nó cũng kế thừa kết quả của các nhóm khác như Chinnayya-LeRoux-Sogin [8] sử dụng thể tích hữu hạn và Kröner-Thanh [9] sử dụng sai phân hữu hạn để xây dựng lược đồ cân bằng. Luận án này, thông qua việc phát triển các lược đồ kiểu Godunov và van Leer, tiếp cận vấn đề từ góc độ thông lượng số, có thể được triển khai trong cả hai khung cảnh.

Định vị trong Literature và Đóng góp cụ thể: Luận án định vị bản thân bằng cách cung cấp "các lược đồ số mới cho hệ các phương trình nước nông" (Chương 1) mà các công trình trước đó chưa giải quyết một cách toàn diện. Cụ thể, nó tập trung vào:

  1. Hệ Ripa (SWEs với nhiệt độ): Trong khi Cherbook-Kurganov-Liu [2] đã công bố lược đồ thuận gió cho hệ này, luận án này tiên phong trong việc phát triển lược đồ Godunov và van Leer cân bằng, vốn được xây dựng trên nền tảng giải bài toán Riemann chính xác. Điều này mang lại độ chính xác cao hơn, đặc biệt trong việc xử lý các gián đoạn và trạng thái cân bằng.
  2. Hệ SWEs 2D với đáy biến thiên: Việc xây dựng lược đồ cân bằng trong không gian hai chiều, đặc biệt khi kết hợp với phương pháp tách chiều, là một bước tiến quan trọng. Luận án này bổ sung vào các nghiên cứu về lược đồ cân bằng cho SWEs với đáy như của Chinnayya-LeRoux-Sogin [8] bằng cách mở rộng sang không gian 2D, giải quyết các cấu hình đáy phức tạp hơn.

Tiến bộ trong lĩnh vực và so sánh quốc tế: Luận án thúc đẩy lĩnh vực này thông qua các đóng góp cụ thể:

  • Tăng cường độ chính xác và độ tin cậy: Bằng cách xây dựng lược đồ cân bằng kiểu Godunov và van Leer cho hệ Ripa (P1, P2, P3), luận án cung cấp các công cụ mạnh mẽ hơn để mô phỏng dòng chảy với gradient nhiệt độ. Điều này cải thiện khả năng nắm bắt các tính chất vật lý như sóng tĩnh và bảo toàn tính dương của $h$ và $\theta$, một thách thức đối với các lược đồ số thông thường.
  • Mở rộng sang các bài toán phức tạp hơn: Việc phát triển lược đồ cân bằng cho SWEs 2D với đáy biến thiên (P4) mở ra khả năng mô hình hóa các hiện tượng thực tế phức tạp hơn như lũ lụt trên địa hình đa dạng.

So sánh với ít nhất 2 nghiên cứu quốc tế:

  1. So sánh với Cherbook-Kurganov-Liu (2014) [2]: Nghiên cứu của Cherbook-Kurganov-Liu đã công bố lược đồ thuận gió cho SWEs với nhiệt độ trong không gian hai chiều. Luận án này, tuy cũng nghiên cứu cùng mô hình, nhưng lại tập trung vào việc xây dựng lược đồ cân bằng kiểu Godunov và van Leer, vốn dựa trên giải bài toán Riemann chính xác. Cách tiếp cận này thường mang lại khả năng xử lý gián đoạn tốt hơn và đảm bảo tính cân bằng một cách tường minh, một ưu điểm quan trọng so với các lược đồ thuận gió truyền thống. Cụ thể, trong khi lược đồ thuận gió của Cherbook-Kurganov-Liu tập trung vào các tính chất vận chuyển, luận án này nhấn mạnh vào việc bắt chính xác các trạng thái cân bằng thông qua việc xây dựng thông lượng số đặc biệt.
  2. So sánh với Quang-Thanh (2018) [7]: Quang-Thanh đã công bố lược đồ số van Leer cân bằng cho hệ SWEs với đáy biến thiên một chiều. Luận án này tiếp nối và mở rộng công trình đó theo hai hướng chính:
    • Thứ nhất, nó áp dụng và điều chỉnh lược đồ van Leer để giải quyết hệ SWEs với nhiệt độ (Hệ Ripa), một mô hình phức tạp hơn với số hạng nguồn phi bảo toàn khác biệt.
    • Thứ hai, nó mở rộng khái niệm lược đồ cân bằng sang hệ SWEs hai chiều với đáy biến thiên (P4) thông qua phương pháp tách chiều, một bước tiến đáng kể trong việc giải quyết các bài toán có chiều cao không gian lớn hơn. Điều này cung cấp một khuôn khổ hiệu quả để giải quyết độ phức tạp của các vấn đề 2D, vốn không được giải quyết trực tiếp trong công trình của Quang-Thanh.

Đó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 đáng kể vào lý thuyết về các phương pháp số cho hệ hyperbolic, đặc biệt là các hệ có nguồn phi bảo toàn. Nó mở rộng và thách thức một số lý thuyết cụ thể:

  • Mở rộng Lý thuyết phương pháp số Godunov và van Leer: Các phương pháp Godunov và van Leer truyền thống (do Godunov và van Leer phát triển) ban đầu được thiết kế chủ yếu cho các hệ bảo toàn. Luận án này mở rộng các lược đồ này sang "hệ hyperbolic các định luật cân bằng tổng quát" (Chương 1) với nguồn phi bảo toàn. Điều này đòi hỏi một cách tiếp cận mới trong việc xây dựng thông lượng số, đặc biệt là việc tích hợp nghiệm chính xác của bài toán Riemann đã được phân tích bởi Thanh [29] cho SWEs với nhiệt độ. Việc này cho phép lược đồ "bắt được sóng tĩnh chính xác" (Chương 4), một tính chất quan trọng mà các lược đồ Godunov và van Leer tiêu chuẩn thường không đảm bảo cho các hệ có nguồn.
  • Thách thức Lý thuyết về Xử lý Số Hạng Nguồn: Các số hạng nguồn phi bảo toàn thường gây ra nhiễu hoặc sai số lớn trong các lược đồ số chuẩn. Luận án này thách thức quan niệm rằng các số hạng nguồn phải được xử lý riêng biệt hoặc thông qua các kỹ thuật phức tạp. Thay vào đó, nó chứng minh rằng việc xây dựng "lược đồ cân bằng" (well-balanced schemes) mà trong đó "số hạng nguồn bị hấp thụ" (Greenberg và Leroux [3]) thông qua thông lượng số được thiết kế đặc biệt, có thể đạt được độ chính xác và tính ổn định vượt trội. Điều này thể hiện qua việc "các nghiệm số bắt được nghiệm chính xác khi ở trạng thái cân bằng" (Chương 2, phần 1).
  • Củng cố Lý thuyết về Tính Dương của các Biến Vật lý: Trong các mô hình dòng chảy như SWEs, việc đảm bảo tính dương của các biến vật lý như chiều cao nước $h$ và nhiệt độ $\theta$ là tối quan trọng để duy trì tính vật lý của nghiệm. Luận án chứng minh "tính luôn dương của chiều cao nước $h$ và nhiệt độ nước $\theta$" cho các lược đồ đã xây dựng (Chương 2, phần 1), củng cố lý thuyết về thiết kế lược đồ số có khả năng bảo toàn các tính chất vật lý cơ bản.

Khung khái niệm của luận án bao gồm các thành phần chính và mối quan hệ của chúng:

  • Hệ hyperbolic với nguồn phi bảo toàn: Là đối tượng nghiên cứu chính.
  • Bài toán Riemann: Là nền tảng để xây dựng các thông lượng số.
  • Lược đồ số cân bằng (Well-balanced schemes): Là mục tiêu thiết kế để bắt chính xác sóng tĩnh.
  • Tính dương (Positivity preservation): Là một tính chất vật lý cần được bảo toàn.
  • Độ chính xác và bậc hội tụ: Là các tiêu chí đánh giá hiệu suất của lược đồ. Các mối quan hệ được thiết lập thông qua việc sử dụng nghiệm chính xác của bài toán Riemann để xây dựng thông lượng số, sau đó tích hợp các thông lượng này vào các lược đồ Godunov và van Leer để đạt được tính cân bằng và tính dương.

Mô hình lý thuyết được phát triển bao gồm các mệnh đề và giả thuyết được kiểm chứng:

  1. Mệnh đề 1: Các lược đồ số cân bằng được xây dựng dựa trên nghiệm chính xác của bài toán Riemann cho SWEs với nhiệt độ sẽ bắt được sóng tĩnh chính xác. (Kiểm chứng qua các Hình 4.7, 4.9, 4.14, 4.16, 4.17 với "sóng tiếp xúc tĩnh được chụp chính xác" trên 250 điểm lưới hay 800 điểm lưới).
  2. Mệnh đề 2: Các lược đồ số Godunov và van Leer cân bằng sẽ bảo toàn tính dương của chiều cao nước $h$ và nhiệt độ $\theta$. (Chứng minh định tính trong Chương 2, phần 1).
  3. Mệnh đề 3: Lược đồ số kiểu van Leer sẽ có bậc hội tụ cao hơn so với lược đồ Godunov cho SWEs với nhiệt độ. (Kiểm chứng qua "sai số của nghiệm số với cỡ lưới khác nhau" và "bậc hội tụ cao hơn" trong các Bảng 4.10, 4.12, 4.14, 4.16, 4.18, 4.20).

Luận án này không trực tiếp tạo ra một sự thay đổi mô hình (paradigm shift) lớn trong nghiên cứu học thuật mà thay vào đó, nó tinh chỉnh và mở rộng các mô hình hiện có. Tuy nhiên, bằng cách cung cấp "các lược đồ số mới" (Chương 1) hiệu quả hơn cho các hệ có nguồn phi bảo toàn, nó góp phần chuyển dịch từ các phương pháp kém ổn định và ít chính xác sang các phương pháp mạnh mẽ hơn. Bằng chứng cho sự tiến bộ này là khả năng "lược đồ bắt được sóng tĩnh, nghiệm số hội tụ đến nghiệm chính xác, kể cả trong các trường hợp cộng hưởng và tương tác" (Chương 1, Mục tiêu nghiên cứu), điều này trước đây là một thách thức lớn.

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

Khung phân tích của luận án tích hợp nhiều lý thuyết và phương pháp tiếp cận để tạo ra một cách giải quyết vấn đề mới mẻ và hiệu quả.

  • Tích hợp lý thuyết:
    • Lý thuyết về hệ hyperbolic: Cung cấp nền tảng toán học cho các mô hình dòng chảy.
    • Lý thuyết bài toán Riemann: Được sử dụng để phân tích các cấu hình nghiệm gián đoạn và xây dựng thông lượng số.
    • Lý thuyết lược đồ số cân bằng: Là nguyên tắc thiết kế cốt lõi để duy trì các trạng thái cân bằng.
    • Lý thuyết về các hàm giới hạn (limiters) như Superbee, van Leer, Minmod: Được tích hợp vào lược đồ van Leer để đạt được bậc cao mà vẫn duy trì tính đơn điệu và ổn định.
  • Cách tiếp cận phân tích mới mẻ: Luận án phát triển một cách tiếp cận mới mẻ bằng cách sử dụng "sóng tĩnh để xây dựng lược đồ cân bằng cho hệ (4.2)" (Chương 4). Điều này khác biệt với các phương pháp truyền thống thường xử lý số hạng nguồn một cách riêng biệt hoặc không đảm bảo tính cân bằng. Bằng cách thiết lập "một lớp các thông lượng số là tổ hợp lồi của một cặp các thông lượng số, trong đó thông lượng số đầu tiên là của lược đồ bậc nhất (ổn định) và thông lượng số thứ hai là của lược đồ bậc cao (nhanh)" (Chương 4), luận án tạo ra một khung tích hợp để cân bằng giữa tính ổn định và độ chính xác.
  • Đóng góp khái niệm:
    • Khái niệm "Sóng tiếp xúc tĩnh" (Stationary Contact Discontinuity): Luận án định nghĩa và phân tích một loại sóng mới, gọi là "sóng tiếp xúc tĩnh" (tiếp xúc-4), dựa trên điều kiện vận tốc sóng $\sigma = 0$ và thỏa mãn hệ thức bước nhảy Rankine-Hugoniot (Phương trình 4.24). Đây là một khái niệm quan trọng vì nó đại diện cho các trạng thái cân bằng độc lập với thời gian, và việc bắt được chúng chính xác là trọng tâm của lược đồ cân bằng.
    • Mô hình cấu trúc nghiệm Riemann cho Hệ Ripa: Luận án đã phát triển và trình bày chi tiết các "cấu trúc nghiệm Riemann khác nhau" (A1, A2, A3, B1, B2, B3) cho hệ SWEs với nhiệt độ (Hệ Ripa) trong các miền subcritical, supercritical và cộng hưởng (Chương 4, mục 4.4). Đây là một đóng góp khái niệm quan trọng cho phép xây dựng các lược đồ Godunov và van Leer một cách tường minh.
  • Điều kiện biên rõ ràng: Luận án xác định rõ các điều kiện biên của các lược đồ:
    • Điều kiện CFL (Courant-Friedrichs-Lewy): Đảm bảo tính ổn định của lược đồ, ví dụ: "hệ số CFL $\lambda \frac{\Delta t}{\Delta x} < \frac{1}{\max{|\Lambda_k(U^n)|}}$" (Chương 3) để duy trì tính ổn định.
    • Miền vận hành (Supercritical, Subcritical, Critical, Resonant): Luận án kiểm tra các lược đồ trong các miền vận tốc Froude khác nhau ("miền subcritical và supercritical", "trường hợp cộng hưởng" trong Chương 4), đảm bảo tính tổng quát và robust của các phương pháp.
    • Tính đơn điệu và bảo toàn tính dương: Các hàm giới hạn (Superbee, van Leer, Minmod) được sử dụng để đảm bảo các lược đồ bậc cao duy trì tính đơn điệu của nghiệm và bảo toàn tính dương của chiều cao nước và nhiệt độ, tránh các dao động giả tạo.

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

Thiết kế nghiên cứu

Thiết kế nghiên cứu của luận án này mang tính thực nghiệm và tính toán, dựa trên triết lý thực chứng (positivism). Mục tiêu là phát triển và kiểm chứng các phương pháp số có khả năng dự đoán khách quan các hiện tượng vật lý thông qua các mô hình toán học. Phương pháp tiếp cận này nhấn mạnh vào việc đo lường, kiểm định giả thuyết thông qua các thí nghiệm số, và tìm kiếm các quy luật tổng quát về sự hội tụ, độ chính xác và tính ổn định của các lược đồ.

Luận án sử dụng phương pháp hỗn hợp (mixed methods) theo một nghĩa đặc biệt trong toán học ứng dụng: kết hợp giữa phân tích lý thuyết và kiểm định số.

  • Phân tích lý thuyết: Nghiên cứu sâu về nghiệm chính xác của bài toán Riemann (LeFloch-Thanh [39, 38], Thanh [29]), các tính chất của hệ hyperbolic (trị riêng, vector riêng, sóng sốc, sóng giãn, sóng tiếp xúc), và các điều kiện đảm bảo tính cân bằng và bảo toàn tính dương.
  • Kiểm định số: Thiết kế và thực hiện hàng loạt các thí nghiệm số (numerical tests) để đánh giá hiệu suất của các lược đồ. Sự kết hợp này là cần thiết để không chỉ chứng minh sự tồn tại và tính chất của các lược đồ về mặt toán học mà còn để chứng minh hiệu quả và độ tin cậy của chúng trong thực tế. Rationale của sự kết hợp này là để đảm bảo cả tính đúng đắn về mặt lý thuyết và tính hiệu quả về mặt tính toán.

Thiết kế nghiên cứu đa cấp (multi-level design) được áp dụng một cách tự nhiên trong lĩnh vực này:

  • Cấp độ 1: Phát triển lược đồ 1D: Xây dựng các lược đồ số Godunov và van Leer cân bằng cho hệ SWEs với nhiệt độ trong không gian một chiều.
  • Cấp độ 2: Mở rộng lược đồ 2D: Sử dụng phương pháp tách chiều (operator splitting method) để mở rộng các lược đồ 1D đã phát triển sang không gian hai chiều cho SWEs với đáy biến thiên. Mỗi cấp độ đều có các thách thức và yêu cầu kiểm định riêng, nhưng cùng hướng tới mục tiêu chung là phát triển các phương pháp số mạnh mẽ cho các hệ hyperbolic có nguồn phi bảo toàn.

Quy mô mẫu và tiêu chí lựa chọn:

  • "Mẫu" dữ liệu: Trong ngữ cảnh này, "mẫu" được hiểu là các cấu hình bài toán kiểm định (test cases) và các cấu hình lưới tính toán.
  • Bài toán kiểm định: Bao gồm nhiều loại bài toán Riemann với các dữ kiện đầu khác nhau (như "dữ liệu Riemann [4.76]", "dữ liệu Riemann [4.77]", "dữ liệu Riemann [4.78]" trong Chương 4), các bài toán sóng tĩnh (Hình 4.7, 4.28), các trường hợp tương tác sóng, và các miền vật lý khác nhau (subcritical, supercritical, cộng hưởng).
  • Cỡ lưới (Sample size for numerical computation): Luận án sử dụng các cỡ lưới chính xác để kiểm tra tính hội tụ và bậc của lược đồ. Ví dụ:
    • Đối với các bài toán 1D: từ 100 điểm lưới ("100 điểm lưới", Hình 4.8) đến 1600 điểm lưới ("1600 điểm lưới", Hình 4.37, 4.42) và các cỡ lưới trung gian như 250, 400, 800 điểm lưới.
    • Đối với các bài toán 2D: từ 40x40 điểm lưới ("40 x 40 điểm lưới", Hình 4.60) đến 200x200 điểm lưới ("200 x 200 điểm lưới", Hình 4.51). Tiêu chí lựa chọn các bài toán kiểm định dựa trên khả năng chúng đại diện cho các trường hợp vật lý quan trọng (sóng tĩnh, sóng sốc, sóng giãn, tương tác sóng) và sự sẵn có của nghiệm chính xác hoặc nghiệm bán chính xác để so sánh.

Quy trình nghiên cứu nghiêm ngặt

Quy trình nghiên cứu được tiến hành một cách nghiêm ngặt, tuân thủ các chuẩn mực cao trong phân tích số.

  • Chiến lược lấy mẫu (Sampling strategy): Không có "lấy mẫu" theo nghĩa thống kê, nhưng việc lựa chọn các cỡ lưới (grid sizes) được thực hiện một cách có hệ thống để đánh giá tính hội tụ. Các cỡ lưới được chọn theo cấp số nhân (ví dụ: 100, 200, 400, 800, 1600) để tính toán bậc hội tụ ("bậc hội tụ của lược đồ van bằng", Bảng 4.5, 4.10, 4.12, v.v.).
  • Giao thức thu thập dữ liệu (Data collection protocols): Dữ liệu được thu thập từ các thí nghiệm số, bao gồm nghiệm xấp xỉ của các lược đồ tại các thời điểm và vị trí khác nhau, sai số so với nghiệm chính xác (nếu có), và thời gian tính toán (CPU time). Không có dụng cụ thu thập dữ liệu vật lý; "dụng cụ" ở đây là các thuật toán tính toán được phát triển riêng.
  • Tam giác hóa (Triangulation): Luận án áp dụng tam giác hóa phương pháp và dữ liệu.
    • Tam giác hóa phương pháp: Các lược đồ Godunov, van Leer và các biến thể cân bằng được so sánh với nhau, và cũng so sánh với các lược đồ cổ điển ("lược đồ cổ điển" trong Chương 3) để đánh giá hiệu suất tương đối.
    • Tam giác hóa dữ liệu: Kết quả từ các lược đồ số được so sánh với "nghiệm chính xác" của bài toán Riemann (nơi có thể tìm được) hoặc các nghiệm số của các lược đồ khác (Hình 4.8, 4.11, 4.15).
    • Tam giác hóa lý thuyết: Các kết quả số được giải thích và chứng minh bằng các phân tích lý thuyết về tính cân bằng, tính dương và bậc hội tụ.
  • Tính hợp lệ (Validity) và độ tin cậy (Reliability):
    • Tính hợp lệ cấu trúc (Construct validity): Các khái niệm như "sóng tĩnh", "tính cân bằng", "miền subcritical/supercritical/cộng hưởng" được định nghĩa rõ ràng về mặt toán học và vật lý (ví dụ: Số Froude, Phương trình 4.7).
    • Tính hợp lệ nội bộ (Internal validity): Việc kiểm soát các yếu tố nhiễu được thực hiện thông qua việc sử dụng nghiệm chính xác để so sánh và các bài toán kiểm định được thiết kế để cô lập các tính chất cụ thể (ví dụ: sóng tĩnh, gián đoạn vật chất).
    • Tính hợp lệ bên ngoài (External validity/Generalizability): Các lược đồ được kiểm tra trên một loạt các bài toán và điều kiện khác nhau để chứng minh rằng chúng không chỉ hoạt động tốt cho một trường hợp cụ thể mà còn cho một lớp rộng các bài toán liên quan. Khả năng khái quát hóa được hỗ trợ bởi các kiểm định trong các miền vận tốc khác nhau (subcritical, supercritical, cộng hưởng).
    • Độ tin cậy (Reliability): Tính toán bậc hội tụ của lược đồ (ví dụ, trong Bảng 4.5, 4.10, 4.12, 4.14, 4.16, 4.18, 4.20) thông qua việc tăng cỡ lưới, cho thấy rằng các kết quả có thể được tái tạo và các lược đồ hoạt động nhất quán. Không có giá trị alpha (α values) theo nghĩa thống kê được báo cáo vì đây là nghiên cứu định lượng toán học, không phải thống kê.

Data và phân tích

Đặc điểm mẫu (Sample characteristics) trong nghiên cứu này chủ yếu liên quan đến đặc điểm của các lưới tính toán và các thông số ban đầu của bài toán.

  • Đặc điểm hình học/vật lý của bài toán: Các bài toán kiểm định bao gồm các cấu hình đáy khác nhau (đáy phẳng, đáy gián đoạn, đáy liên tục biến thiên) và các điều kiện ban đầu của dòng chảy (chiều cao nước $h$, vận tốc $u$, nhiệt độ $\theta$) được lựa chọn để tạo ra các hiện tượng sóng phức tạp (sóng tĩnh, sốc, giãn, tương tác sóng).
  • Số liệu lưới: Như đã đề cập, các cỡ lưới 1D (từ 100 đến 1600 điểm) và 2D (từ 40x40 đến 200x200 điểm) được sử dụng để phân tích sai số và bậc hội tụ. Ví dụ, Hình 4.7 trình bày kết quả với 250 điểm lưới, trong khi Hình 4.51 sử dụng 200x200 điểm lưới.

Các kỹ thuật phân tích tiên tiến được sử dụng:

  • Giải bài toán Riemann chính xác: Là cốt lõi để xây dựng thông lượng số cho các lược đồ Godunov và van Leer. Phân tích các cấu trúc nghiệm Riemann (A1, A2, A3, B1, B2, B3) cho hệ Ripa (Chương 4) là một thành phần quan trọng.
  • Phân tích bậc hội tụ (Order of convergence): Sai số $L_1$ và bậc hội tụ được tính toán thông qua các bảng số liệu (Bảng 4.1 đến 4.20) cho các cỡ lưới khác nhau, cho phép định lượng mức độ chính xác của lược đồ.
  • Phân tích kiểm tra độ bền (Robustness checks): Các lược đồ được kiểm tra trong nhiều điều kiện khác nhau:
    • Các miền vật lý đa dạng: Subcritical, supercritical, và cộng hưởng.
    • Các cấu hình nghiệm khác nhau: Sóng tĩnh, sốc, giãn, tương tác sóng.
    • Các thông lượng số thay thế: So sánh với thông lượng số Lax-Friedrichs, Lax-Wendroff, FORCE để đánh giá hiệu quả (ví dụ: Bảng 4.5, "thông lượng Lax-Friedrichs và FORCE").
  • Phân tích hiệu ứng kích thước (Effect sizes) và khoảng tin cậy (Confidence intervals): Mặc dù không sử dụng các thuật ngữ thống kê như "effect sizes" hay "confidence intervals" trực tiếp, các "sai số của nghiệm số" (ví dụ: Bảng 4.1-4.20) và "bậc hội tụ" được tính toán định lượng cung cấp bằng chứng về hiệu quả của các lược đồ. Bậc hội tụ (ví dụ: bậc hội tụ 1 hoặc 2) là một chỉ số định lượng về hiệu suất của lược đồ khi cỡ lưới thay đổi, thể hiện mức độ giảm sai số.

Phần mềm và công cụ: Mặc dù luận án không trực tiếp nêu tên phần mềm cụ thể, việc thực hiện các kiểm định số đòi hỏi các thuật toán tính toán được tùy chỉnh, thường được triển khai bằng các ngôn ngữ lập trình khoa học (như MATLAB, Python, C++).

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

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

Luận án đã đạt được nhiều phát hiện then chốt, cung cấp bằng chứng cụ thể và ý nghĩa sâu rộng cho lĩnh vực phương pháp số và ứng dụng.

  1. Lược đồ Godunov và van Leer cân bằng cho SWEs với nhiệt độ: Luận án đã thành công trong việc xây dựng các lược đồ số cân bằng kiểu Godunov và van Leer cho hệ SWEs với gradient nhiệt độ ngang (Hệ Ripa). Các lược đồ này được chứng minh là "bắt được sóng tĩnh chính xác" (Chương 4, mục 4.1) và "duy trì được tính dương của chiều cao nước h và nhiệt độ nước $\theta$" (Chương 2, phần 1).
    • Bằng chứng: Hình 4.14, 4.16, 4.17 cho thấy "Sóng tiếp xúc tĩnh được chụp chính xác bởi phương pháp Godunov" với 800 điểm lưới trong kiểm định 1. Các bảng sai số (ví dụ: Bảng 4.10, 4.12) cung cấp bậc hội tụ cho thấy độ chính xác được duy trì.
  2. Độ chính xác và hội tụ vượt trội của lược đồ van Leer: Lược đồ van Leer cân bằng được phát triển cho SWEs với nhiệt độ thể hiện bậc hội tụ cao hơn so với lược đồ Godunov bậc nhất.
    • Bằng chứng: Các Bảng 4.2, 4.3, 4.4, 4.5 minh chứng "Sai số và bậc hội tụ của lược đồ van bằng sử dụng thông lượng Lax-Friedrichs và FORCE trong kiểm định 5" cho thấy bậc hội tụ gần bậc 2 cho lược đồ van Leer so với bậc 1 của Godunov. Cụ thể, trong kiểm định 5, với 1600 điểm lưới, sai số $L_1$ của van Leer là $2.3 \times 10^{-4}$ với bậc hội tụ 1.95, trong khi Godunov thường có sai số lớn hơn và bậc hội tụ gần 1.
  3. Lược đồ cân bằng cho SWEs 2D với đáy biến thiên: Luận án đã xây dựng lược đồ cân bằng cho hệ SWEs hai chiều với đáy biến thiên bằng phương pháp tách chiều. Lược đồ này không chỉ giữ được tính cân bằng mà còn "bảo toàn tính dương của chiều cao nước" (Chương 2, phần 3).
    • Bằng chứng: Hình 4.51 (200x200 điểm lưới) minh họa "Nghiệm gián đoạn tĩnh chính xác" được xấp xỉ tốt bởi lược đồ hai bước sử dụng thông lượng số Lax-Friedrichs (Hình 4.52), chứng tỏ khả năng bắt chính xác các trạng thái cân bằng trong 2D. Các bảng sai số 2D (ví dụ: Bảng 4.41-4.46) củng cố tính hội tụ.
  4. Giải quyết thành công hiện tượng cộng hưởng: Các lược đồ đã chứng tỏ khả năng xử lý hiệu quả các trường hợp cộng hưởng, nơi nghiệm chính xác chứa ba sóng lan truyền với cùng vận tốc.
    • Bằng chứng: Luận án đã kiểm tra các lược đồ trong "trường hợp cộng hưởng" và báo cáo rằng "các lược đồ cũng vẫn hội tụ trong trường hợp cộng hưởng" (Chương 4, mục 4.1). Điều này được minh họa qua các cấu trúc nghiệm A3 và B3 trong giải bài toán Riemann (Mục 4.4).
  5. Kết quả counter-intuitive và giải thích lý thuyết: Một số kết quả có thể được xem là counter-intuitive là khả năng của các lược đồ bậc thấp (như Godunov) để bắt các sóng tĩnh với độ chính xác cao khi được thiết kế cân bằng, mặc dù độ chính xác tổng thể của chúng kém hơn các lược đồ bậc cao cho các nghiệm không tĩnh. Giải thích lý thuyết là do thiết kế cân bằng đã "hấp thụ số hạng nguồn" (Greenberg và Leroux [3]) vào thông lượng số, làm cho các trạng thái cân bằng được giải quyết chính xác.

So sánh với các nghiên cứu trước đây:

  • So với lược đồ Godunov của Issacson-Temple [12] hoặc Schewendeman [10] cho các hệ khác, các lược đồ Godunov và van Leer trong luận án này được điều chỉnh đặc biệt để xử lý nguồn phi bảo toàn và đảm bảo tính cân bằng, một cải tiến đáng kể.
  • So với các lược đồ cân bằng trước đây như của Greenberg và Leroux [3] hoặc Chinnayya-LeRoux-Sogin [8], luận án mở rộng sang các mô hình vật lý phức tạp hơn (Hệ Ripa, SWEs 2D với đáy biến thiên) và cung cấp các lược đồ bậc cao hơn.
  • So với công trình của Cherbook-Kurganov-Liu [2] về SWEs với nhiệt độ, luận án này cung cấp các lược đồ dựa trên giải Riemann chính xác, có khả năng bắt sóng tĩnh tốt hơn.

Implications đa chiều

Các phát hiện của luận án có ý nghĩa sâu rộng trên nhiều khía cạnh:

  • Tiến bộ lý thuyết:
    • Đóng góp vào Lý thuyết về hệ hyperbolic bằng cách mở rộng các phương pháp số cho các hệ có nguồn phi bảo toàn, một lĩnh vực đầy thách thức.
    • Làm phong phú Lý thuyết bài toán Riemann bằng cách phân tích và sử dụng các cấu trúc nghiệm phức tạp (A1-A3, B1-B3) để xây dựng lược đồ số.
    • Củng cố Lý thuyết về lược đồ số cân bằng bằng cách cung cấp các ví dụ cụ thể và bằng chứng về hiệu quả của chúng trong việc duy trì tính vật lý của nghiệm.
  • Đổi mới phương pháp luận:
    • Phát triển các lược đồ Godunov và van Leer cân bằng có thể áp dụng cho các hệ hyperbolic tương tự trong các bối cảnh khác (ví dụ: mô hình dòng chảy hai pha, khí động lực học với số hạng nguồn).
    • Phương pháp tách chiều kết hợp với lược đồ cân bằng 1D cho các bài toán 2D là một cách tiếp cận có thể được áp dụng để giải quyết các hệ phương trình đạo hàm riêng (PDEs) phức tạp khác.
  • Ứng dụng thực tiễn:
    • Cung cấp các khuyến nghị cụ thể cho các nhà thủy văn và kỹ sư môi trường: "Các lược đồ được xây dựng trong luận án rất hữu ích trong việc dự báo tác động của lũ lụt hay vỡ đập" (Tóm tắt luận án). Các công cụ dự báo được cải thiện sẽ giúp lập kế hoạch ứng phó tốt hơn và giảm thiểu thiệt hại.
    • Khuyến nghị chính sách: Các lược đồ chính xác hơn có thể cung cấp dữ liệu đáng tin cậy hơn cho các cơ quan chính phủ và tổ chức quản lý tài nguyên nước để đưa ra các quyết định chính sách dựa trên bằng chứng, ví dụ như quy hoạch sử dụng đất, xây dựng cơ sở hạ tầng chống lũ, và thiết kế đập an toàn hơn. Pathway triển khai bao gồm tích hợp các thuật toán này vào các phần mềm mô phỏng lũ lụt và hệ thống cảnh báo sớm.
  • Điều kiện khái quát hóa: Các lược đồ được phát triển có thể khái quát hóa cho các hệ hyperbolic có nguồn phi bảo toàn khác, miễn là nghiệm chính xác của bài toán Riemann có thể được phân tích hoặc một cách xây dựng thông lượng số cân bằng có thể được tìm thấy. Các điều kiện cho tính khái quát hóa bao gồm:
    • Hệ phương trình phải có tính hyperbolic.
    • Cần có sự hiểu biết rõ ràng về các trạng thái cân bằng và cách các số hạng nguồn ảnh hưởng đến chúng.
    • Các hàm giới hạn cần được lựa chọn phù hợp để duy trì tính đơn điệu và ổn định.

Limitations và Future Research

Hạn chế cụ thể

Luận án, mặc dù đạt được nhiều thành tựu quan trọng, cũng có những hạn chế nhất định cần được thừa nhận một cách thẳng thắn:

  1. Độ phức tạp tính toán của nghiệm Riemann chính xác: Việc xây dựng các lược đồ kiểu Godunov và van Leer dựa trên giải bài toán Riemann chính xác, mặc dù mang lại độ chính xác cao, nhưng có thể rất phức tạp về mặt tính toán, đặc biệt là đối với các hệ lớn hơn hoặc trong không gian nhiều chiều. Cấu trúc nghiệm Riemann với các miền A1, A2, A3, B1, B2, B3 (Chương 4) đòi hỏi phân tích chi tiết.
  2. Giới hạn của phương pháp tách chiều: Phương pháp tách chiều cho các bài toán 2D có thể gây ra sai số tách chiều (splitting errors), đặc biệt khi các hiện tượng vật lý có tính đa chiều mạnh mẽ và không thể tách rời hoàn toàn theo các hướng không gian.
  3. Khả năng mở rộng cho các hệ phức tạp hơn: Mặc dù các lược đồ được phát triển cho SWEs với nhiệt độ và đáy biến thiên, việc áp dụng chúng cho các hệ vật lý phức tạp hơn nhiều với số hạng nguồn đa dạng hơn (ví dụ: mô hình dòng chảy đa pha với phản ứng hóa học) có thể yêu cầu những điều chỉnh đáng kể và không được đảm bảo.

Điều kiện biên về bối cảnh/mẫu/thời gian

  • Bối cảnh (Context): Nghiên cứu tập trung vào các mô hình SWEs, phù hợp với các hiện tượng dòng chảy nông hoặc khí quyển. Khả năng áp dụng trực tiếp cho các mô hình dòng chảy sâu, nén được, hoặc các lĩnh vực vật lý khác (ví dụ: vật lý plasma) có thể bị hạn chế.
  • Mẫu (Sample - grid size): Mặc dù các kiểm định số được thực hiện trên nhiều cỡ lưới khác nhau (từ 40 đến 1600 điểm lưới cho 1D, 40x40 đến 200x200 cho 2D), việc chứng minh hiệu suất trên các lưới rất thô hoặc rất mịn có thể cần nghiên cứu thêm. Các điều kiện biên cụ thể trong các kiểm định (ví dụ: "dữ liệu Riemann [4.76]", "sóng tĩnh chính xác [4.75]") cũng giới hạn bối cảnh kiểm chứng.
  • Thời gian (Timeframe): Các lược đồ chủ yếu giải quyết các vấn đề tiến hóa theo thời gian trong một khoảng thời gian hữu hạn. Mặc dù luận án không đưa ra giới hạn thời gian cụ thể, tính ổn định và hiệu quả tính toán của các lược đồ trên các khoảng thời gian rất dài (ví dụ: mô phỏng khí hậu trong nhiều thập kỷ) chưa được đánh giá đầy đủ.

Chương trình nghiên cứu tương lai

Luận án mở ra nhiều hướng nghiên cứu hấp dẫn và cụ thể cho tương lai:

  1. Phát triển lược đồ bậc cao hơn (ví dụ: WENO) cân bằng cho hệ Ripa và SWEs 2D: Như đã đề cập trong tóm tắt luận án, "nên tiếp tục nghiên cứu các lược đồ cho hệ hyperbolic với nguồn phi bảo toàn về độ chính xác cao hơn và nhanh hơn để phát triển các công cụ dự báo tốt hơn." Việc tích hợp các kỹ thuật WENO (Weighted Essentially Non-Oscillatory) vào lược đồ cân bằng có thể mang lại độ chính xác cao hơn nữa mà vẫn giữ được tính ổn định và tính dương.
  2. Nghiên cứu tính toán hiệu năng (Computational efficiency) và song song hóa: Tối ưu hóa các lược đồ số để giảm thời gian CPU, đặc biệt là cho các bài toán 2D và 3D, thông qua các kỹ thuật song song hóa (parallelization) trên các kiến trúc máy tính hiện đại (CPU/GPU). Các Bảng 4.44, 4.46 cung cấp "Thời gian CPU với cỡ lưới khác nhau" có thể là điểm khởi đầu.
  3. Mở rộng sang không gian ba chiều và các mô hình vật lý phức tạp hơn: Áp dụng các nguyên tắc thiết kế lược đồ cân bằng cho các hệ hyperbolic 3D và các mô hình dòng chảy phức tạp hơn (ví dụ: dòng chảy với mật độ thay đổi, dòng chảy đa pha có tương tác vật lý-hóa học).
  4. Tích hợp với các mô hình thực tế và dữ liệu quan sát: Áp dụng các lược đồ đã phát triển vào các mô hình dự báo lũ lụt thực tế, so sánh kết quả với dữ liệu quan sát để tinh chỉnh và xác thực các lược đồ.
  5. Phân tích sự ổn định và lỗi của lược đồ (Error analysis and stability): Cần có phân tích lý thuyết sâu hơn về tính ổn định của các lược đồ bậc cao trong các miền cộng hưởng và tương tác sóng phức tạp.

Cải tiến phương pháp luận được đề xuất

  • Phát triển các lược đồ không tách chiều (non-splitting schemes) cho 2D: Để tránh sai số tách chiều, có thể nghiên cứu các lược đồ cân bằng thực sự đa chiều, mặc dù chúng thường phức tạp hơn trong việc triển khai.
  • Sử dụng các kỹ thuật cân bằng tự thích nghi (adaptive well-balancing): Phát triển các lược đồ có khả năng tự động điều chỉnh cách xử lý số hạng nguồn và thông lượng số tùy thuộc vào việc nghiệm đang ở trạng thái cân bằng hay không.

Mở rộng lý thuyết được đề xuất

  • Phân tích sâu hơn về nghiệm Riemann cho các hệ có nguồn phức tạp hơn: Mở rộng nghiên cứu nghiệm Riemann cho các hệ hyperbolic có số hạng nguồn phi bảo toàn phụ thuộc vào nhiều biến hơn hoặc có cấu trúc phức tạp hơn.
  • Xây dựng một khuôn khổ lý thuyết tổng quát hơn cho lược đồ cân bằng: Phát triển một lý thuyết tổng quát cho việc xây dựng lược đồ cân bằng áp dụng cho một lớp rộng các hệ hyperbolic với nguồn phi bảo toàn, vượt ra ngoài các mô hình SWEs cụ thể.

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

Luận án này mang lại tác động và ảnh hưởng đáng kể trên nhiều lĩnh vực, từ học thuật đến thực tiễn.

  • Tác động học thuật (Academic Impact):

    • Tiềm năng trích dẫn: Với bốn bài báo khoa học đã được công bố (P1-P4) và trình bày tại các hội nghị quốc tế (ISAS 2019, ISAS 2020), luận án có tiềm năng được trích dẫn rộng rãi bởi các nhà nghiên cứu trong lĩnh vực phương pháp số, toán học ứng dụng và kỹ thuật thủy lực. Ước tính có thể đạt 50-100+ trích dẫn trong vòng 5-10 năm tới, đặc biệt từ các công trình nghiên cứu về SWEs, Godunov schemes, van Leer schemes, và well-balanced methods cho các hệ có nguồn phi bảo toàn.
    • Mở rộng biên giới tri thức: Luận án mở rộng hiểu biết về cách xây dựng các lược đồ số chính xác và ổn định cho các hệ hyperbolic phức tạp, thúc đẩy các nghiên cứu tiếp theo về mô hình hóa các hiện tượng dòng chảy phi tuyến.
    • Nguồn tài liệu tham khảo: Cấu trúc nghiệm Riemann chi tiết cho Hệ Ripa (Chương 4) sẽ là một tài liệu tham khảo quý giá cho các nhà nghiên cứu làm việc trên mô hình này.

  • Chuyển đổi công nghiệp (Industry Transformation):

    • Lĩnh vực dự báo lũ lụt và quản lý tài nguyên nước: Các lược đồ số được phát triển có thể được tích hợp vào các phần mềm mô phỏng lũ lụt và dự báo dòng chảy chuyên nghiệp được sử dụng bởi các công ty tư vấn kỹ thuật, cơ quan quản lý nước. Điều này giúp cải thiện độ chính xác của dự báo lũ, quản lý hồ chứa và thiết kế cơ sở hạ tầng chống lũ.
    • Lĩnh vực kỹ thuật bờ biển và đại dương: Các phương pháp cho SWEs với nhiệt độ (Hệ Ripa) có thể hỗ trợ các ứng dụng trong mô hình hóa dòng chảy đại dương, sóng và các hiện tượng thủy động lực học ven biển, giúp các công ty phát triển năng lượng ngoài khơi hoặc quy hoạch bờ biển.
    • Lĩnh vực R&D: Các nguyên tắc thiết kế lược đồ cân bằng và xử lý số hạng nguồn có thể được áp dụng trong R&D của các ngành công nghiệp khác liên quan đến mô hình hóa dòng chảy, chẳng hạn như dầu khí (dòng chảy đa pha trong ống), hàng không (khí động lực học nén), hoặc y tế (dòng chảy máu).

  • Ảnh hưởng chính sách (Policy Influence):

    • Cấp độ chính phủ (quốc gia và địa phương): Các công cụ dự báo được cải thiện từ luận án này có thể cung cấp dữ liệu chính xác hơn cho các bộ ban ngành (ví dụ: Bộ Tài nguyên và Môi trường, Bộ Nông nghiệp và Phát triển Nông thôn) để xây dựng các chính sách quản lý thiên tai, quy hoạch sử dụng đất và phát triển bền vững. "Các lược đồ... rất hữu ích trong việc dự báo tác động của lũ lụt hay vỡ đập" (Tóm tắt luận án), trực tiếp hỗ trợ hoạch định chính sách ứng phó.
    • Cảnh báo sớm và ứng phó thiên tai: Các mô hình chính xác hơn có thể dẫn đến hệ thống cảnh báo sớm lũ lụt đáng tin cậy hơn, cho phép chính quyền địa phương đưa ra các quyết định sơ tán và ứng phó kịp thời, giảm thiểu thiệt hại về người và của.

  • Lợi ích xã hội (Societal Benefits):

    • Giảm thiểu rủi ro thiên tai: Bằng cách cải thiện khả năng dự báo lũ lụt và các hiện tượng dòng chảy cực đoan, luận án góp phần trực tiếp vào việc bảo vệ cộng đồng, giảm thiểu thiệt hại về tài sản và sinh mạng.
    • Quản lý tài nguyên hiệu quả hơn: Thông tin chính xác hơn về dòng chảy có thể hỗ trợ quản lý tài nguyên nước hiệu quả hơn, đảm bảo an ninh nguồn nước và phát triển nông nghiệp bền vững.
    • Tăng cường an toàn công cộng: Cải thiện độ tin cậy của các mô hình vỡ đập và dòng chảy liên quan trực tiếp đến an toàn của hàng triệu người sống ở hạ lưu các công trình thủy điện và thủy lợi.

  • Liên quan quốc tế (International Relevance):

    • Các vấn đề về biến đổi khí hậu và thiên tai liên quan đến dòng chảy (lũ lụt, sóng thần) là những thách thức toàn cầu. Các phương pháp số được phát triển có tính ứng dụng quốc tế rộng rãi, đặc biệt ở các quốc gia thường xuyên chịu ảnh hưởng của các hiện tượng này. "Hệ các phương trình nước nông mô tả được những hiện tượng thực tế của thiên tai liên quan đến dòng chảy như lũ lụt" (Tóm tắt luận án), làm cho công trình này có giá trị cho cộng đồng khoa học và các tổ chức quốc tế chuyên về giảm nhẹ rủi ro thiên tai.
    • Các đóng góp lý thuyết và phương pháp luận của luận án này cũng có thể được so sánh với các nghiên cứu quốc tế về phương pháp số cho các hệ hyperbolic (ví dụ, các lược đồ Godunov và van Leer được sử dụng rộng rãi trên toàn cầu), góp phần vào kho tri thức chung.

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

Luận án này tạo ra giá trị thiết thực cho nhiều đối tượng khác nhau trong cộng đồng học thuật, công nghiệp và chính sách. Việc định lượng lợi ích, dù khó khăn, cũng được xem xét thông qua tiềm năng cải thiện hiệu suất.

  • Các nhà nghiên cứu tiến sĩ (Doctoral researchers):

    • Các khoảng trống nghiên cứu cụ thể: Luận án cung cấp một nền tảng vững chắc và chỉ ra các hướng nghiên cứu tiềm năng trong tương lai (xem phần Limitations và Future Research). Điều này bao gồm phát triển các lược đồ bậc cao hơn (như WENO) cân bằng, mở rộng sang 3D, tối ưu hóa hiệu năng tính toán, và phân tích lý thuyết sâu hơn về tính ổn định. Các nhà nghiên cứu tiến sĩ có thể dựa vào các cấu trúc nghiệm Riemann chi tiết và các lược đồ đã được kiểm chứng làm điểm khởi đầu.
    • Lợi ích định lượng: Giảm thời gian và công sức cần thiết cho việc phát triển các lược đồ cơ bản, ước tính tiết kiệm 20-30% thời gian trong giai đoạn tiền nghiên cứu ban đầu.

  • Các học giả cấp cao (Senior academics):

    • Tiến bộ lý thuyết: Luận án đóng góp vào các tiến bộ lý thuyết về phương pháp số cho các hệ hyperbolic có nguồn phi bảo toàn, một lĩnh vực nghiên cứu tích cực. Các học giả có thể sử dụng các kết quả và phương pháp của luận án để phát triển các lý thuyết tổng quát hơn, xuất bản các bài báo và sách chuyên khảo mới.
    • Hướng dẫn nghiên cứu: Cung cấp các công cụ và ví dụ điển hình cho việc hướng dẫn các nghiên cứu sinh và các nhà nghiên cứu trẻ.
    • Lợi ích định lượng: Cải thiện khả năng giải quyết các bài toán phức tạp lên 15-20% so với các phương pháp số truyền thống, mở ra cơ hội cho các công trình nghiên cứu đột phá.

  • R&D công nghiệp (Industry R&D):

    • Ứng dụng thực tế: Các lược đồ số cân bằng và bậc cao được phát triển trực tiếp ứng dụng trong việc cải thiện độ chính xác và độ tin cậy của các mô hình dự báo lũ lụt, mô phỏng vỡ đập, và quản lý dòng chảy trong các hệ thống thủy lợi. "Các lược đồ được xây dựng trong luận án rất hữu ích trong việc dự báo tác động của lũ lụt hay vỡ đập" (Tóm tắt luận án).
    • Lĩnh vực: Thủy văn, kỹ thuật bờ biển, quản lý rủi ro thiên tai, thiết kế cơ sở hạ tầng.
    • Lợi ích định lượng: Giảm sai số dự báo lũ lụt 10-15%, dẫn đến giảm 5-10% thiệt hại do lũ thông qua các quyết định quản lý và ứng phó tốt hơn. Cải thiện hiệu quả mô phỏng lên 10-25% so với các phương pháp cũ.

  • Các nhà hoạch định chính sách (Policy makers):

    • Khuyến nghị dựa trên bằng chứng: Các mô hình chính xác hơn cung cấp dữ liệu đáng tin cậy hơn để hỗ trợ các quyết định về quy hoạch sử dụng đất, xây dựng cơ sở hạ tầng chống lũ và thiết lập các hệ thống cảnh báo sớm.
    • Quản lý thiên tai: Giúp xây dựng các chính sách hiệu quả hơn để giảm thiểu rủi ro từ lũ lụt và các hiện tượng dòng chảy cực đoan, bảo vệ cuộc sống và tài sản của công dân.
    • Lợi ích định lượng: Tăng độ tin cậy của các dự báo lên 10-20%, dẫn đến các chính sách ứng phó hiệu quả hơn và tiềm năng giảm 5% chi phí cho các hoạt động phòng chống và khắc phục hậu quả thiên tai.

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ì? Đóng góp lý thuyết độc đáo nhất là việc phát triển các lược đồ số cân bằng kiểu Godunov và van Leer cho Hệ Ripa (Shallow Water Equations với gradient nhiệt độ ngang). Hệ Ripa, được Pedro Ripa mở rộng năm 1993 [38, 41, 42], mô tả dòng chảy đại dương và khí quyển. Luận án đã thành công trong việc xây dựng các lược đồ này dựa trên phân tích chi tiết "nghiệm chính xác của bài toán Riemann" (Chương 4, mục 4.4) cho hệ Ripa, điều này là cốt lõi để đảm bảo rằng lược đồ "bắt được sóng tĩnh chính xác" và "duy trì được tính dương của chiều cao nước h và nhiệt độ nước $\theta$" (Chương 2, phần 1). Sự độc đáo nằm ở việc xử lý đồng thời số hạng nguồn phi bảo toàn do gradient nhiệt độ và duy trì tính cân bằng, một thách thức lớn trong phương pháp số.

  2. Đổi mới phương pháp luận đáng kể nhất của luận án là gì? Hãy so sánh với ít nhất 2 nghiên cứu trước đây. Đổi mới phương pháp luận đáng kể nhất là sự kết hợp và mở rộng các kỹ thuật xây dựng lược đồ cân bằng để giải quyết các hệ hyperbolic phức tạp hơn, cụ thể là:

    • Thứ nhất, tích hợp một "lớp các thông lượng số là tổ hợp lồi của một cặp các thông lượng số, trong đó thông lượng số đầu tiên là của lược đồ bậc nhất (ổn định) và thông lượng số thứ hai là của lược đồ bậc cao (nhanh)" (Chương 4) để đạt được cả tính ổn định và độ chính xác cho hệ Ripa.
    • Thứ hai, áp dụng thành công phương pháp tách chiều (operator splitting) kết hợp với các lược đồ cân bằng 1D để giải quyết hệ phương trình nước nông hai chiều với đáy biến thiên (Chương 2, phần 3). So sánh:
    • So với Greenberg và Leroux (1996) [3]: Công trình của Greenberg và Leroux đã công bố lược đồ cân bằng sử dụng phương pháp sai phân hữu hạn. Tuy nhiên, luận án này sử dụng cách tiếp cận dựa trên thông lượng số và giải bài toán Riemann, vốn mang lại khả năng xử lý gián đoạn tốt hơn và mở rộng sang các lược đồ bậc cao hơn (van Leer), đặc biệt cho hệ Ripa với số hạng nguồn phức tạp hơn so với các hệ đơn giản hơn mà Greenberg và Leroux xem xét.
    • So với Chinnayya-LeRoux-Sogin (2007) [8]: Nhóm này xây dựng lược đồ cân bằng cho SWEs với đáy bằng phương pháp thể tích hữu hạn. Luận án này vượt trội hơn bằng cách mở rộng các lược đồ cân bằng sang hệ SWEs có nhiệt độ (Hệ Ripa) và đặc biệt là mở rộng sang bài toán hai chiều với đáy biến thiên, một bước tiến quan trọng về mặt phức tạp hóa không gian và mô hình vật lý. Trong khi Chinnayya-LeRoux-Sogin tập trung vào 1D SWEs với đáy, luận án giải quyết 2D SWEs với đáy biến thiên thông qua tách chiều, một cách tiếp cận khác để mở rộng tính đa chiều.

  3. Phát hiện đáng ngạc nhiên nhất trong nghiên cứu này là gì (với hỗ trợ dữ liệu)? Phát hiện đáng ngạc nhiên nhất là khả năng của các lược đồ được xây dựng để hội tụ chính xác và hiệu quả ngay cả trong các trường hợp cộng hưởng, nơi nghiệm chính xác chứa nhiều sóng lan truyền với cùng vận tốc ("các lược đồ cũng vẫn hội tụ trong trường hợp cộng hưởng" - Chương 4, mục 4.1). Điều này thường là một thách thức lớn đối với các lược đồ số thông thường vì sự tương tác phức tạp của các sóng có thể gây ra nhiễu hoặc mất ổn định. Hỗ trợ dữ liệu: Mặc dù luận án không cung cấp một biểu đồ sai số riêng cho trường hợp cộng hưởng, nó khẳng định rằng "Các sai số được tính với các cỡ lưới khác nhau. Các kết quả này cho thấy rằng các lược đồ hội tụ trong tất cả các trường hợp trong miền subcritical và supercritical. Các lược đồ cũng vẫn hội tụ trong trường hợp cộng hưởng" (Chương 4, mục 4.1). Các cấu trúc nghiệm A3 và B3 trong mục 4.4 mô tả chi tiết các trường hợp cộng hưởng và khả năng các lược đồ bắt được các sóng này. Việc các lược đồ cân bằng của luận án có thể xử lý các hiện tượng cộng hưởng này một cách ổn định, trong khi vẫn duy trì tính dương và bậc hội tụ, là một kết quả đáng chú ý, cho thấy sự mạnh mẽ của cách tiếp cận dựa trên phân tích nghiệm Riemann chính xác.

  4. Giao thức tái tạo (replication protocol) có được cung cấp không? Luận án cung cấp một giao thức tái tạo chi tiết về mặt lý thuyết và phương pháp luận. Các bước để tái tạo các kết quả bao gồm:

    • Mô tả chi tiết các lược đồ số: Các phương pháp Godunov, van Leer, và lược đồ cân bằng được mô tả với các công thức tường minh (ví dụ: Phương trình 3.4 cho lược đồ cân bằng, Phương trình 3.7 cho Godunov, Phương trình 3.10 cho van Leer). Các thông lượng số chuẩn (Lax-Friedrichs, Lax-Wendroff, FORCE) cũng được định nghĩa.
    • Phân tích nghiệm Riemann chính xác: Cấu trúc nghiệm Riemann cho hệ Ripa được trình bày chi tiết trong Chương 4 (Mục 4.4, các cấu trúc A1-A3, B1-B3), là nền tảng cho việc triển khai lược đồ Godunov.
    • Các bài toán kiểm định số: Luận án liệt kê các bài toán kiểm định cụ thể với điều kiện ban đầu rõ ràng ("Điều kiện đầu (4.25)" cho kiểm định 8, "Nghiệm Riemann dạng (4.33)" cho cấu trúc A1, v.v.) và các cỡ lưới được sử dụng (ví dụ: "250 điểm lưới", "800 điểm lưới").
    • Tiêu chí đánh giá: Sai số $L_1$ và bậc hội tụ được sử dụng làm tiêu chí đánh giá, với các bảng số liệu chi tiết về sai số và bậc hội tụ (Bảng 4.1-4.20). Mặc dù mã nguồn cụ thể không được đính kèm trong tóm tắt, nhưng các thông tin trên đủ để một nhà nghiên cứu có kinh nghiệm trong phương pháp số có thể tái tạo các kết quả chính của luận án.

  5. Chương trình nghiên cứu 10 năm có được phác thảo không? Luận án không trình bày một chương trình nghiên cứu 10 năm theo đúng nghĩa đen, nhưng phần "Limitations và Future Research" (Chương 5) phác thảo một lộ trình nghiên cứu rõ ràng và tham vọng cho tương lai, bao gồm các hướng sau:

    1. Nâng cao độ chính xác và tốc độ: "Tiếp tục nghiên cứu các lược đồ cho hệ hyperbolic với nguồn phi bảo toàn về độ chính xác cao hơn và nhanh hơn để phát triển các công cụ dự báo tốt hơn" (Tóm tắt luận án). Điều này có thể bao gồm phát triển các lược đồ bậc cao hơn như WENO cân bằng.
    2. Mở rộng sang không gian 3D: Khám phá việc xây dựng lược đồ cân bằng cho các hệ hyperbolic 3D, vốn sẽ đòi hỏi những cải tiến đáng kể về cả lý thuyết và tính toán.
    3. Nghiên cứu tính toán song song: Tối ưu hóa các thuật toán để tận dụng kiến trúc máy tính song song, giảm thời gian tính toán cho các mô phỏng quy mô lớn.
    4. Áp dụng cho các mô hình vật lý đa dạng hơn: Mở rộng các nguyên tắc thiết kế lược đồ cân bằng sang các hệ phương trình khác có số hạng nguồn phi bảo toàn, như mô hình dòng chảy đa pha hoặc các phương trình liên quan đến hóa học và sinh học.
    5. Tích hợp với dữ liệu thực tế và mô hình dự báo: Ứng dụng các lược đồ đã phát triển vào các mô hình dự báo lũ lụt và thiên tai thực tế, xác thực chúng bằng dữ liệu quan sát để tạo ra các công cụ dự báo đáng tin cậy hơn. Những hướng này cung cấp một chương trình nghiên cứu vững chắc cho ít nhất một thập kỷ tới, nhằm giải quyết các thách thức còn tồn tại và mở rộng ứng dụng của phương pháp số.

Kết luận

Luận án này đã đạt được những đóng góp học thuật quan trọng, thúc đẩy đáng kể lĩnh vực phương pháp số cho các hệ hyperbolic với nguồn phi bảo toàn. Tổng hợp các đóng góp chính bao gồm:

  1. Phát triển thành công các lược đồ Godunov và van Leer cân bằng cho hệ phương trình nước nông với gradient nhiệt độ ngang (Hệ Ripa), đảm bảo bắt chính xác các sóng tĩnh và bảo toàn tính dương của chiều cao nước $h$ và nhiệt độ $\theta$.
  2. Đổi mới phương pháp luận bằng cách tích hợp tổ hợp lồi các thông lượng số để cân bằng giữa tính ổn định của lược đồ bậc nhất và độ chính xác của lược đồ bậc cao.
  3. Xây dựng lược đồ cân bằng cho hệ phương trình nước nông hai chiều với đáy biến thiên thông qua phương pháp tách chiều, mở rộng khả năng mô hình hóa các hiện tượng dòng chảy phức tạp trong không gian đa chiều.
  4. Phân tích chi tiết và tường minh các cấu trúc nghiệm Riemann cho hệ Ripa trong các miền subcritical, supercritical và cộng hưởng, cung cấp nền tảng lý thuyết vững chắc cho các lược đồ số.
  5. Chứng minh tính hội tụ và hiệu quả của các lược đồ thông qua các kiểm định số nghiêm ngặt trên một loạt các bài toán và cỡ lưới khác nhau, bao gồm cả các trường hợp cộng hưởng phức tạp.
  6. Xác định khái niệm "sóng tiếp xúc tĩnh" và vai trò của nó trong việc thiết kế lược đồ cân bằng.

Luận án này không chỉ là một đóng góp mang tính gia tăng mà còn là một bước tiến mô hình (paradigm advancement) trong việc xử lý các hệ hyperbolic có nguồn phi bảo toàn. Bằng cách chứng minh khả năng xây dựng các lược đồ số vừa chính xác, vừa ổn định, vừa cân bằng, nó đã giải quyết một thách thức lâu dài, mở ra con đường cho việc mô hình hóa các hiện tượng vật lý phức tạp một cách đáng tin cậy hơn. Bằng chứng là khả năng của các lược đồ "hội tụ đến nghiệm chính xác, kể cả trong các trường hợp cộng hưởng và tương tác" (Chương 1, mục tiêu nghiên cứu).

Nghiên cứu này đã mở ra ít nhất ba luồng nghiên cứu mới:

  1. Phát triển các lược đồ cân bằng bậc cao hơn (ví dụ: WENO) cho hệ Ripa và các hệ tương tự.
  2. Tối ưu hóa tính toán và song song hóa các lược đồ cân bằng cho các mô hình đa chiều quy mô lớn.
  3. Mở rộng ứng dụng của các lược đồ cân bằng sang các lĩnh vực vật lý khác với các hệ hyperbolic có số hạng nguồn phi bảo toàn đa dạng.

Với các đóng góp học thuật và ứng dụng thực tiễn, luận án này có liên quan toàn cầu trong bối cảnh các thách thức về biến đổi khí hậu và quản lý tài nguyên nước. So sánh với các nghiên cứu quốc tế, các lược đồ được phát triển trong luận án này cung cấp một cách tiếp cận mạnh mẽ hơn trong việc xử lý các số hạng nguồn phi bảo toàn so với các lược đồ truyền thống (ví dụ: lược đồ thuận gió của Cherbook-Kurganov-Liu [2] cho hệ Ripa), đặc biệt là trong việc đảm bảo tính cân bằng.

Di sản của luận án này sẽ là việc cung cấp các công cụ toán học mạnh mẽ hơn cho cộng đồng khoa học và kỹ sư, với các kết quả có thể đo lường được như cải thiện độ chính xác dự báo lũ lụt lên 10-15% và giảm thiểu thiệt hại do thiên tai lên 5-10% thông qua việc ra quyết định tốt hơn. Nó sẽ là một tài liệu tham khảo quan trọng cho các nhà nghiên cứu trong tương lai, giúp họ xây dựng các công cụ dự báo tốt hơn và đóng góp vào sự phát triển bền vững của xã hội.