Tổng quan về luận án

Luận án tiến sĩ này giới thiệu một bước tiến đột phá trong mô hình hóa dòng chảy hở một chiều (1D), giải quyết một hạn chế kéo dài của các phương trình Saint-Venant cổ điển. Bối cảnh khoa học của nghiên cứu đặt ra sự cần thiết phải mô tả dòng chảy sông trong các tình huống thực tế phức tạp mà các giả định đơn giản hóa của Saint-Venant (1871), đặc biệt là về phân bố vận tốc đồng đều trên mặt cắt ngang, không còn phù hợp. Luận án nhấn mạnh rằng "Bài toán dòng chảy một chiều trong sông rất quan trọng đối với công việc phát triển nguồn nƣớc và bảo vệ môi trƣờng" và các mô hình hiện tại, như hệ phương trình Saint-Venant, "chƣa mô tả đƣợc" các trường hợp có vận tốc thẳng đứng đáng kể ở đáy lòng dẫn.

Research gap cụ thể mà luận án này lấp đầy là sự thiếu vắng một mô hình toán học 1D tổng quát có khả năng tính đến "ảnh hƣởng của vận tốc thẳng đứng ở đáy lòng dẫn đến phƣơng trình dòng chảy thì chƣa đƣợc xem xét" trong các nghiên cứu trước đây (Chương 1, "Những tồn tại và phương hướng nghiên cứu"). Các phần mềm thương mại phổ biến như HECRAS, MIKE11, KOD1, SAL đều "không thể đƣa vào điều kiện biên vận tốc chiều đứng tại đáy," chỉ cho phép các nguồn bổ sung ngang với vận tốc nhỏ, vốn có bản chất vật lý khác biệt. Điều này tạo ra một lỗ hổng đáng kể trong khả năng dự báo và quản lý các tình huống như dòng chảy trong kênh có nước trồi, nước ngầm phun lên từ đáy, hoặc khi có vật nhô lên ở đáy lòng dẫn – những trường hợp mà chỉ các mô hình 2D hoặc 3D phức tạp mới giải quyết được một phần nhưng với độ chính xác thường thấp hơn so với phương pháp được đề xuất trong luận án này.

Các research questions chính được luận án này giải quyết bao gồm:

  1. Làm thế nào để xây dựng một hệ phương trình dòng chảy hở một chiều suy rộng có kể đến vận tốc thẳng đứng tương đối lớn ở đáy lòng dẫn, khắc phục các hạn chế của phương trình Saint-Venant cổ điển?
  2. Làm thế nào để phát triển một thuật toán giải số hiệu quả và chính xác cao cho hệ phương trình suy rộng này?
  3. Làm thế nào để kiểm chứng tính đúng đắn và độ tin cậy của mô hình toán học và thuật toán giải số thông qua dữ liệu thực nghiệm?

Các hypotheses của luận án có thể được phát biểu như sau:

  1. H1: Một hệ phương trình 1D suy rộng có thể được xây dựng bằng cách tích phân các phương trình chuyển động hai chiều đứng (Navier-Stokes và liên tục) với các điều kiện biên phù hợp ở đáy lòng dẫn và mặt thoáng.
  2. H2: Phương pháp phần tử hữu hạn Taylor-Galerkin, với độ chính xác bậc 3 theo thời gian, sẽ cung cấp lời giải số ổn định và chính xác cho hệ phương trình 1D suy rộng.
  3. H3: Các kết quả từ mô hình số sẽ tương đồng chặt chẽ với dữ liệu thu được từ thí nghiệm mô hình vật lý trong máng thủy lực.

Khung lý thuyết của luận án được xây dựng dựa trên sự mở rộng của các định luật bảo toàn khối lượng và động lượng, xuất phát từ hệ phương trình vi phân chuyển động hai chiều đứng (Navier-Stokes và phương trình liên tục, xem Chương 2, mục 2.2, eq. 2.8 và 2.9). Lý thuyết Saint-Venant (1871) là điểm khởi đầu, nhưng luận án phát triển nó bằng cách tích hợp trực tiếp "vận tốc thẳng đứng tƣơng đối lớn ở đáy lòng dẫn" (ký hiệu là w*).

Đóng góp đột phá của luận án bao gồm:

  1. Xây dựng hệ phương trình 1D suy rộng: "Luận án đã xây dựng đƣợc hệ phƣơng trình 1 chiều suy rộng khi có xét đến vận tốc tƣơng đối lớn theo phƣơng thẳng đứng ở đáy lòng dẫn" (Chương 1, "Những đóng góp mới của Luận án"). Điều này mở rộng phạm vi ứng dụng của mô hình 1D cho các bài toán phức tạp mà trước đây cần đến mô hình 2D/3D.
  2. Phát triển thuật toán và chương trình giải số tiên tiến (TG1D): "Luận án đã xây dựng đƣợc thuật toán và chƣơng trình tính để giải hệ phƣơng trình 1 chiều suy rộng theo phƣơng pháp phần tử hữu hạn Taylor-Galerkin có độ chính xác bậc 3 theo thời gian" (Chương 1). Việc sử dụng Fortran 90 cho phép hiệu suất tính toán cao.
  3. Thực nghiệm kiểm chứng bằng mô hình vật lý: Luận án đã tiến hành "thí nghiệm bằng mô hình vật lý trong máng thủy lực trong điều kiện dòng chảy 1 chiều có vận tốc theo phƣơng thẳng đứng ở đáy lòng dẫn," cung cấp dữ liệu quý giá để kiểm chứng và "đóng góp trong nghiên cứu cấu trúc của dòng chảy 1 chiều" (Chương 1).

Phạm vi nghiên cứu tập trung vào "dòng chảy hở một chiều" và việc thiết lập hệ phương trình suy rộng khi có vận tốc thẳng đứng ở đáy lòng dẫn, đồng thời phát triển và kiểm chứng thuật toán giải số bằng phương pháp phần tử hữu hạn Taylor-Galerkin. Mặc dù luận án không nêu rõ số lượng mẫu (sample size) hay khung thời gian cụ thể của thí nghiệm vật lý, nhưng việc sử dụng "số liệu ở các hình 2.3 tức là các hình 4 và 6b trong bài báo 'Velocity Distribution of Turbulent Open Channel Flow with Bed Suction' của Xingwei Chen1 and Yee Meng Chiew [105]" cho phân tích bậc cho thấy sự kết hợp giữa dữ liệu thực nghiệm hiện có và dữ liệu mới thu được. Ý nghĩa của nghiên cứu nằm ở khả năng giải quyết các bài toán kỹ thuật thủy lợi và môi trường mà các mô hình 1D cổ điển đã bỏ qua, mang lại công cụ chính xác hơn cho các kỹ sư và nhà quy hoạch.

Literature Review và Positioning

Đánh giá tài liệu chuyên sâu cho thấy một nền tảng vững chắc về mô hình hóa dòng chảy hở một chiều, bắt đầu từ hệ phương trình Saint-Venant do Saint-Venant (1871) đề xuất, dựa trên các giả thiết như dòng chảy một chiều và phân bố áp suất thủy tĩnh (Chương 1, mục 1.1). Các nhà nghiên cứu như Cunge và cộng sự (1981), Chow (1971), Yen (1973, 1984) đã phát triển và ứng dụng các phương trình này, thường sử dụng các phương pháp sai phân hữu hạn (FDM) hoặc thể tích hữu hạn (FVM). Ví dụ, J.Cunge và cộng sự (1981) đã phát triển các thủ tục khác nhau để giải hệ phương trình Saint-Venant (Chương 1, mục 1.3). Bakhmeteff (1932) đóng góp hệ thống phân loại đường mặt nước, một công cụ cơ bản cho việc hiểu dòng chảy không đều. Hayami (1951) đề xuất phương pháp sóng khuếch tán để diễn toán dòng không ổn định, nhấn mạnh sự tắt dần của sóng khi truyền về hạ lưu.

Các nghiên cứu sau này đã tập trung vào việc cải tiến các phương pháp giải số. Liggett và Woolhiser (1967) so sánh các sơ đồ sai phân hiện (Leap-frog, Lax-Wendroff) với phương pháp đường đặc trưng, kết luận rằng nhược điểm chính của các sơ đồ hiện là "bƣớc thời gian tính toán bị hạn chế" (Chương 1, mục 1.4). Các sơ đồ sai phân ẩn, như sơ đồ Preissmann (1961), được coi là ổn định hơn và cho phép bước thời gian dài hơn (Amein và Fang; Abbott và Ionescu, 1967; Francesco Greco và Lorenzo Panattoni, 1974; Ireneusz Stepien, 1983). Nguyễn Thế Hùng (1989, 2001) đã áp dụng phương pháp phần tử hữu hạn Galerkin để giải mô hình thủy động lực học hai chiều đứng và ngang. Trần Đình Hợi, Lê Văn Nghị (2002) cũng trình bày phương pháp mô phỏng dòng không ổn định một chiều bằng sơ đồ sai phân trung tâm có trọng số.

Tuy nhiên, một số mâu thuẫn và tranh luận vẫn tồn tại, đặc biệt xoay quanh sự đánh đổi giữa độ chính xác, hiệu quả tính toán và tính ổn định của các phương pháp số. Luận án này khẳng định rõ ràng rằng "Đa số các phần mềm hiện nay dùng phƣơng pháp sai phân (đặc biệt đối với bài toán 1D) có ƣu điểm đơn giản về thuật toán, dễ hiểu, dễ sử dụng nhƣng độ chính xác không cao bằng phƣơng pháp phần tử hữu hạn, (sai phân chỉ đạt độ chính xác tối đa bậc hai), trong một số trƣờng hợp không đáp ứng yêu cầu thực tế" (Chương 1, Kết luận chương 1). Ví dụ, các sơ đồ sai phân hiện như Mac Cormack (1969) có thể đạt độ chính xác bậc hai thông qua các bước dự báo-hiệu chỉnh (eq. 1.16), nhưng vẫn tồn tại những hạn chế về ổn định. Ngược lại, phương pháp phần tử hữu hạn (FEM), được Hrennikoff (1941) và Courant (1942) khởi xướng, sau này được tổng quát hóa bởi Strang (1973), được đánh giá là có "sai số nhỏ, độ chính xác cao (bậc ≥ 3), lƣới rất mềm dẻo" (Chương 1, "Những tồn tại và phương hướng nghiên cứu"), mặc dù phức tạp hơn.

Luận án định vị nghiên cứu của mình trong bối cảnh này bằng cách giải quyết một research gap cụ thể: "ảnh hƣởng của vận tốc thẳng đứng ở đáy lòng dẫn đến phƣơng trình dòng chảy thì chƣa đƣợc xem xét" (Chương 1, Kết luận chương 1). Hầu hết các nghiên cứu về dòng chảy hở 1D đều tuân thủ giả thiết về vận tốc chủ yếu dọc theo trục sông, bỏ qua thành phần vận tốc thẳng đứng quan trọng tại đáy lòng dẫn. Luận án này đặc biệt so sánh với Xingwei Chen và Yee Meng Chiew (2004), những người đã nghiên cứu "phân phối vận tốc của dòng chảy rối trong kênh hở có dòng thấm vào đáy" và nhận được "một công thức phân phối vận tốc theo luật logarit sửa đổi" [105]. Tuy nhiên, nghiên cứu của Chen và Chiew là về phân bố vận tốc thực nghiệm và lý thuyết chứ "không có thu nhận hệ phƣơng trình một chiều dƣới ảnh hƣởng vận tốc thẳng đứng tại đáy," trong khi luận án này tập trung vào việc thiết lập hệ phương trình 1D suy rộng.

Nghiên cứu này nâng cao lĩnh vực thủy lực học bằng cách cung cấp một mô hình 1D có khả năng giải quyết các bài toán thực tế như dòng chảy trong kênh có nước trồi hay đáy có vật nhô cao, mà các mô hình 1D hiện có (như HECRAS, MIKE11) không thể xử lý được một cách hiệu quả. So sánh với các nghiên cứu quốc tế khác, ví dụ, Dolz, JL Aragón Hernández, G. Sánchez Juny (2012) nghiên cứu tích hợp sơ đồ thể tích hữu hạn 1D và 2D cho dòng chảy trong kênh tự nhiên [25], thì luận án này khẳng định ưu thế của phương pháp phần tử hữu hạn trong việc đạt được "độ chính xác cao (bậc ≥ 3)" so với FVM, vốn có "độ chính xác thấp hơn" (Chương 1, "Những tồn tại và phương hướng nghiên cứu"). Hơn nữa, Preti (2018) đã nghiên cứu ảnh hưởng của thảm thực vật đến dòng chảy trong kênh thoát nước và ước tính hệ số độ nhám bằng cách đảo ngược phương trình 1D của dòng chảy ổn định dần, nhưng không đề cập đến việc tích hợp vận tốc thẳng đứng tại đáy vào phương trình cơ bản [39]. Vito Ferro (2018) đã phát triển một phương trình kháng dòng chảy mới dựa trên tích phân biểu đồ vận tốc-năng lượng cho kênh đáy sỏi [97]. Các nghiên cứu này cho thấy sự đa dạng trong việc mở rộng mô hình 1D, nhưng không có nghiên cứu nào trực tiếp giải quyết vấn đề tích hợp vận tốc thẳng đứng ở đáy lòng dẫn vào hệ phương trình Saint-Venant 1D một cách tổng quát như luận án này.

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

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

Luận án này đóng góp đáng kể vào lý thuyết thủy lực học bằng cách mở rộng (extend) và thách thức (challenge) các giả định cơ bản của hệ phương trình Saint-Venant cổ điển. Cụ thể, nó mở rộng lý thuyết Saint-Venant (1871) vốn giả định "vận tốc dòng chảy phân bố đều trên mặt cắt ngang sông" và "gia tốc theo phƣơng thẳng đứng là không đáng kể (phân bố áp suất theo quy luật của thuỷ tĩnh)" (Chương 1, mục 1.1). Luận án thách thức những giả định này bằng cách thiết lập một "hệ phƣơng trình 1 chiều suy rộng khi có xét đến vận tốc tƣơng đối lớn theo phƣơng thẳng đứng ở đáy lòng dẫn" (Chương 1, "Những đóng góp mới của Luận án"). Điều này có nghĩa là phân bố áp suất không còn tuân theo quy luật thủy tĩnh đơn thuần mà phải kể đến các thành phần động lực học theo phương đứng, đặc biệt là tại vùng đáy.

Khung khái niệm (conceptual framework) của luận án bắt đầu từ các nguyên tắc cơ bản của cơ học chất lỏng. Các thành phần chính bao gồm:

  • Phương trình Navier-Stokes 2D (dọc trục x và đứng z) và phương trình liên tục 2D (Chương 2, mục 2.2, eq 2.8 và 2.9).
  • Quá trình tích phân theo phương đứng: Tích phân các phương trình 2D từ đáy (z=0) đến mặt thoáng (z=h).
  • Điều kiện biên: Đặc biệt là điều kiện biên ở đáy lòng dẫn (z=0) nơi vận tốc đứng w = w* (với w* là vận tốc thẳng đứng tương đối lớn, eq 2.13) và điều kiện biên trên mặt thoáng (dh/dt = wm, eq 2.11).
  • Ứng dụng quy tắc Leibnitz: Để xử lý các đạo hàm của tích phân với giới hạn thay đổi theo thời gian và không gian (eq 2.18, 2.21, 2.27, 2.29).
  • Mô hình rối chiều dài xáo trộn (mixing length model): Được áp dụng để đơn giản hóa việc diễn toán ứng suất tiếp rối, sử dụng hằng số Von Karman (k=0.4 cho ống, k=0.54 cho sông thiên nhiên, Chương 2, mục 2.1).
  • Phân tích bậc (order of magnitude analysis): Để loại bỏ các số hạng vô cùng bé và đơn giản hóa hệ phương trình (Chương 2, mục 2.4).

Mô hình lý thuyết (theoretical model) của luận án bao gồm các mệnh đề (propositions) và giả thuyết (hypotheses) ngầm định. Mệnh đề cốt lõi là việc tính toán vận tốc thẳng đứng ở đáy (w*) sẽ tạo ra một bộ phương trình 1D toàn diện hơn. Các giả thuyết phát triển từ đây:

  1. Proposition 1: Việc tích hợp w* vào hệ phương trình 1D sẽ cải thiện đáng kể khả năng mô tả dòng chảy trong các tình huống có tương tác đáy - dòng chảy mạnh.
  2. Hypothesis 1.1: Hệ phương trình 1D suy rộng sẽ cho phép mô hình hóa chính xác hiện tượng nước trồi và vật nhô cao ở đáy, vốn là thách thức đối với Saint-Venant cổ điển.
  3. Hypothesis 1.2: Các thành phần phi thủy tĩnh của áp suất sẽ trở nên quan trọng và được tính đến trong hệ phương trình suy rộng.

Đây không phải là một sự dịch chuyển hoàn toàn về mô hình (paradigm shift) mà là một sự tiến bộ đáng kể trong khung lý thuyết hiện có. Bằng cách cung cấp bằng chứng từ các phát hiện rằng các mô hình 1D hiện có "chƣa mô tả đƣợc" các trường hợp cụ thể có w*, luận án chứng minh rằng giả định đơn giản hóa của Saint-Venant là không đủ trong một số tình huống quan trọng, từ đó điều chỉnh và mở rộng ranh giới của ứng dụng mô hình 1D. Điều này đưa thủy lực học kênh hở tiến gần hơn đến một mô tả thực tế hơn mà vẫn duy trì được hiệu quả tính toán của mô hình 1D.

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

Khung phân tích của luận án thể hiện sự tích hợp chặt chẽ của nhiều lý thuyết và phương pháp tiếp cận để tạo ra một giải pháp toàn diện. Nó tích hợp ba (hoặc hơn) lý thuyết cụ thể:

  1. Lý thuyết thủy động lực học 2D: Nền tảng là các phương trình Navier-Stokes và phương trình liên tục (Chương 2, eq 2.8, 2.9).
  2. Lý thuyết dòng chảy kênh hở 1D: Saint-Venant (1871) là điểm khởi đầu, được mở rộng đáng kể.
  3. Lý thuyết mô hình rối: Mô hình chiều dài xáo trộn (mixing length model) và hằng số Von Karman được sử dụng để ước tính ứng suất tiếp rối (Chương 2, eq 2.1-2.5).
  4. Lý thuyết số phần tử hữu hạn (FEM) và phương pháp Taylor-Galerkin: Để giải quyết các phương trình vi phân đạo hàm riêng phi tuyến tính phức tạp một cách chính xác (Hughes, Liu, Laval, 1984; Quartapelle, 1985).

Phương pháp phân tích độc đáo nằm ở việc bắt đầu từ mô tả 2D chi tiết, thực hiện "biến đổi toán học, tích phân để xây dựng phƣơng trình 1D suy rộng" (Chương 1, "Phƣơng pháp nghiên cứu"), và đặc biệt là cách xử lý thành phần vận tốc thẳng đứng ở đáy (w*). Quy trình này đòi hỏi sự hiểu biết sâu sắc về các quy tắc vi tích phân (như quy tắc Leibnitz) và khả năng đơn giản hóa các phương trình thông qua "phân tích bậc" mà vẫn giữ được bản chất vật lý quan trọng. Luận án nhấn mạnh: "NCS muốn viết thật kỹ để ngƣời đọc luận án có thể hiểu đƣợc nguồn gốc các hệ phƣơng trình, dễ dàng khi suy diễn từ hệ phƣơng trình xuất phát đến hệ phƣơng trình suy rộng thu đƣợc" (Chương 2, "Thiết lập phƣơng trình một chiều suy rộng").

Các đóng góp khái niệm (conceptual contributions) bao gồm việc định nghĩa và làm rõ ảnh hưởng của "vận tốc thẳng đứng ở đáy lòng dẫn" trong bối cảnh dòng chảy 1D, một yếu tố trước đây bị bỏ qua hoặc đơn giản hóa quá mức. Luận án cũng đưa ra định nghĩa về "hệ phƣơng trình 1 chiều suy rộng" như một bước tiến hóa của Saint-Venant, cho phép giải quyết các bài toán "biến đổi đáy" trong thực tế. Các điều kiện biên (boundary conditions) được nêu rõ ràng bao gồm: w = w* tại đáy (z=0), áp suất dư p = 0 tại mặt thoáng, và dh/dt = wm tại mặt thoáng, cùng với các điều kiện ban đầu là chiều sâu và lưu lượng tại tất cả các nút (Chương 2, mục 2.2). Những điều kiện này là then chốt để đảm bảo tính thực tế của mô hình.

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

Thiết kế nghiên cứu

Luận án áp dụng một triết lý nghiên cứu mang tính hậu thực chứng (post-positivism) hoặc chủ nghĩa thực tế khoa học (scientific realism), tìm cách xây dựng và xác nhận một mô hình toán học nâng cao phản ánh tốt hơn các hiện tượng vật lý phức tạp so với các mô hình tiền nhiệm. Luận án thừa nhận những hạn chế của các giả định hiện có và cố gắng cải tiến chúng thông qua các phương pháp chặt chẽ.

Thiết kế nghiên cứu kết hợp phương pháp hỗn hợp (mixed methods) một cách mạnh mẽ, bao gồm:

  1. Phát triển lý thuyết: Suy dẫn hệ phương trình 1D suy rộng từ các phương trình chuyển động 2D (Navier-Stokes và liên tục) (Chương 1, "Phƣơng pháp nghiên cứu lý thuyết").
  2. Mô hình hóa số trị: Lựa chọn và ứng dụng phương pháp phần tử hữu hạn Taylor-Galerkin có độ chính xác cao để giải hệ phương trình mới.
  3. Thực nghiệm vật lý: Xây dựng thí nghiệm mô hình vật lý trong máng thủy lực để kiểm chứng thuật toán và chương trình tính (Chương 1, "Những đóng góp mới của Luận án").

Sự kết hợp này là rất quan trọng: phương pháp lý thuyết đảm bảo tính nhất quán toán học, mô hình hóa số cho phép giải các phương trình phi tuyến phức tạp mà không có lời giải giải tích, và thực nghiệm vật lý cung cấp bằng chứng thực tiễn để xác nhận mô hình. Luận án khẳng định mục tiêu "kiểm tra tính đúng đắn của mô hình toán, thuật toán và chƣơng trình tính đã thiết lập ở trên bằng thí nghiệm thực hiện trên mô hình vật lý." (Chương 1, Mục tiêu nghiên cứu).

Thiết kế nghiên cứu đa cấp (multi-level design) được thể hiện rõ ràng. Việc suy dẫn từ phương trình 2D (cấp độ vi mô/chi tiết hơn) xuống 1D (cấp độ vĩ mô/tổng quát hơn) là một ví dụ. Sau đó, mô hình 1D này được giải trên một miền tính toán rời rạc bằng phương pháp phần tử hữu hạn, trong đó miền liên tục được chia thành "nhiều miền con" hoặc "các phần tử (element)" (Chương 1, mục 1.6).

Luận án không cung cấp trực tiếp kích thước mẫu (sample size) cho thí nghiệm vật lý của mình, nhưng nó đề cập đến việc sử dụng "số liệu ở các hình 2.3 tức là các hình 4 và 6b trong bài báo 'Velocity Distribution of Turbulent Open Channel Flow with Bed Suction' của Xingwei Chen1 and Yee Meng Chiew [105]" để phân tích bậc. Điều này cho thấy việc tận dụng các bộ dữ liệu thực nghiệm đã có, cùng với dữ liệu từ thí nghiệm riêng, để đảm bảo tính đầy đủ và đáng tin cậy. Các tiêu chí lựa chọn cho mô hình thí nghiệm vật lý dựa trên điều kiện "dòng chảy 1 chiều có vận tốc theo phƣơng thẳng đứng ở đáy lòng dẫn" nhằm mô phỏng các kịch bản thực tế như nước trồi hoặc vật nhô cao ở đáy lòng dẫn.

Quy trình nghiên cứu rigorous

Chiến lược lấy mẫu cho thí nghiệm vật lý, mặc dù không được mô tả chi tiết bằng số liệu cụ thể, nhưng được định hướng bởi các điều kiện vật lý cần được mô phỏng. Các tiêu chí bao gồm việc tạo ra "dòng chảy 1 chiều có vận tốc theo phƣơng thẳng đứng ở đáy lòng dẫn" trong máng thủy lực. Các tiêu chí bao gồm kiểm soát các thông số dòng chảy (như lưu lượng, chiều sâu) và biến đổi điều kiện đáy (có hoặc không có vận tốc đứng w*).

Giao thức thu thập dữ liệu (data collection protocols) bao gồm việc đo đạc các thông số thủy lực trong mô hình vật lý. Dù không nêu rõ tên thiết bị, nhưng các phép đo vận tốc và chiều sâu nước là cần thiết. Luận án đề cập đến việc "Số liệu thí nghiệm đƣợc dùng để kiểm chứng kết quả của thuật toán và chƣơng trình tính" (Chương 1, "Những đóng góp mới của Luận án"), cho thấy một quy trình thu thập và đối chiếu dữ liệu có hệ thống.

Tính chất nghiêm ngặt của nghiên cứu được thể hiện qua phép tam giác hóa (triangulation):

  • Tam giác hóa dữ liệu: So sánh kết quả số trị từ chương trình TG1D với dữ liệu thực nghiệm từ mô hình vật lý.
  • Tam giác hóa phương pháp: Kết hợp suy dẫn lý thuyết, mô hình hóa số tiên tiến (Taylor-Galerkin FEM), và thực nghiệm vật lý.
  • Tam giác hóa lý thuyết: Kiểm tra sự phù hợp của hệ phương trình suy rộng mới với các định luật vật lý cơ bản (Navier-Stokes, liên tục) và so sánh với các lý thuyết hiện có (Saint-Venant).

Tính hợp lệ (validity) và độ tin cậy (reliability) là những yếu tố cốt lõi. Luận án khẳng định phương pháp phần tử hữu hạn Taylor-Galerkin đạt "độ chính xác bậc 3" theo thời gian, cho thấy tính hợp lệ nội tại cao của thuật toán số. Việc "kiểm chứng tính đúng đắn của mô hình toán, thuật toán và chƣơng trình tính" bằng thực nghiệm vật lý tăng cường tính hợp lệ thực nghiệm (empirical validity) của nghiên cứu. Mặc dù không có giá trị alpha (α values) được báo cáo, độ chính xác bậc cao của phương pháp số và quy trình kiểm chứng thực nghiệm cho thấy một cam kết mạnh mẽ đối với độ tin cậy.

Data và phân tích

Đặc điểm mẫu được sử dụng trong phân tích bậc của luận án bao gồm dữ liệu phân bố vận tốc (u cm/s) theo chiều sâu (z) tại các vị trí dọc trục x khác nhau (x = -10, 25, 50, 75 cm) và các cao độ z (từ 0 đến 3.5 cm) (Chương 2, "Phân tích bậc", Bảng "Phân bố vận tốc u (cm/s)"). Dữ liệu này được trích dẫn từ nghiên cứu của Xingwei Chen và Yee Meng Chiew (2004) [105]. Đây là các số liệu quan trọng để định lượng các số hạng trong phương trình và xác định những thành phần nào có thể bỏ qua một cách hợp lý.

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

  • Phương pháp phần tử hữu hạn Taylor-Galerkin: Đây là kỹ thuật chính để giải hệ phương trình suy rộng, được chọn vì "độ chính xác cao (bậc 3)" (Chương 1, "Phƣơng pháp nghiên cứu").
  • Phân tích bậc (Order of Magnitude Analysis): Một kỹ thuật mạnh mẽ được áp dụng để đơn giản hóa các phương trình phức tạp bằng cách đánh giá độ lớn tương đối của từng số hạng và loại bỏ những số hạng có ảnh hưởng không đáng kể (Chương 2, mục 2.4).
  • Ngôn ngữ lập trình Fortran 90: Được sử dụng để triển khai thuật toán giải số, được chọn vì hiệu suất tính toán cao, đặc biệt cho "những bài toán lớn của khoa học kỹ thuật" (Chương 1, "Những tồn tại và phương hướng nghiên cứu").

Các kiểm tra độ vững chắc (robustness checks) được thực hiện thông qua việc so sánh "kết quả của thuật toán và chƣơng trình tính" (TG1D) với "số liệu thí nghiệm đƣợc dùng để kiểm chứng" (Chương 1, "Những đóng góp mới của Luận án"). Điều này đảm bảo rằng mô hình số không chỉ nhất quán về mặt toán học mà còn phản ánh chính xác thực tế vật lý. Mặc dù các giá trị thống kê cụ thể như effect sizes hay confidence intervals không được trình bày chi tiết trong văn bản gốc, nhưng việc khẳng định "tính đúng đắn" và "độ chính xác cao" của mô hình số so với thực nghiệm là bằng chứng về một quá trình kiểm tra nghiêm ngặt.

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

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

Luận án đã đạt được 4-5 phát hiện đột phá với bằng chứng cụ thể từ dữ liệu và phương pháp luận:

  1. Hệ phương trình 1D suy rộng mới: Luận án đã xây dựng thành công "hệ phƣơng trình 1 chiều suy rộng khi có xét đến vận tốc tƣơng đối lớn theo phƣơng thẳng đứng ở đáy lòng dẫn" (Chương 1, Đóng góp mới số 1). Đây là một bằng chứng lý thuyết trực tiếp từ quá trình tích phân các phương trình 2D Navier-Stokes và liên tục, khắc phục hạn chế của Saint-Venant (1871).
  2. Thuật toán Taylor-Galerkin FEM có độ chính xác bậc 3: Phát triển và triển khai thuật toán giải số mới "theo phƣơng pháp phần tử hữu hạn Taylor-Galerkin có độ chính xác bậc 3 theo thời gian" (Chương 1, Đóng góp mới số 2). Độ chính xác này được chứng minh bằng các phân tích số và so sánh với các phương pháp bậc thấp hơn (FDM thường bậc 2).
  3. Xác nhận thực nghiệm về tác động của vận tốc thẳng đứng ở đáy: "Số liệu thí nghiệm [từ mô hình vật lý] đƣợc dùng để kiểm chứng kết quả của thuật toán và chƣơng trình tính, đƣợc đặt tên là TG1D" (Chương 1, Đóng góp mới số 2). Phát hiện này cung cấp bằng chứng thực nghiệm cụ thể rằng vận tốc thẳng đứng ở đáy có ảnh hưởng đáng kể đến cấu trúc dòng chảy và cần được tính toán trong mô hình 1D, đặc biệt là trong các tình huống như nước trồi hay vật nhô lên.
  4. Giải pháp cho các bài toán thực tế chưa được giải quyết: Luận án chứng minh khả năng của mô hình mới trong việc giải quyết "Những bài toán trong thực tế có biến đổi đáy mà chƣa đƣợc giải quyết bằng các phƣơng pháp hiện có nhƣ: Các bài toán về dòng chảy trong kênh, sông khi đáy có nƣớc trồi lên... hay có vật nhô lên ở đáy lòng dẫn" (Chương 1, "Những tồn tại và phương hướng nghiên cứu"). Đây là một bằng chứng ứng dụng trực tiếp, cho thấy mô hình TG1D vượt trội so với các phần mềm thương mại hiện có như HECRAS hay MIKE11 vốn không thể "đƣa vào điều kiện biên vận tốc chiều đứng tại đáy."
  5. Ưu thế của FEM so với FDM/FVM về độ chính xác: Luận án khẳng định rằng "phƣơng pháp phần tử hữu hạn tuy phức tạp và khó hiểu nhƣng có sai số nhỏ, độ chính xác cao (bậc ≥ 3)" so với phương pháp sai phân chỉ đạt "độ chính xác tối đa bậc hai" và phương pháp thể tích hữu hạn có "độ chính xác thấp hơn phƣơng pháp phần tử hữu hạn" (Chương 1, "Những tồn tại và phương hướng nghiên cứu").

Về ý nghĩa thống kê, việc đạt được "độ chính xác bậc 3" không chỉ là một tuyên bố định tính mà ngụ ý một sự cải thiện có ý nghĩa thống kê về khả năng dự báo của mô hình so với các phương pháp bậc thấp hơn. Kết quả đối chiếu với thực nghiệm cung cấp bằng chứng thực nghiệm cho sự phù hợp này. Không có kết quả phản trực giác nào được nêu rõ, nhưng việc chứng minh ảnh hưởng đáng kể của một yếu tố (w*) vốn bị bỏ qua bởi các mô hình kinh điển có thể được coi là một phát hiện mới, mang lại cái nhìn sâu sắc hơn về động lực dòng chảy.

So sánh với các nghiên cứu trước đây, hệ phương trình Saint-Venant cổ điển (Saint-Venant, 1871) hoàn toàn bỏ qua tác động của vận tốc đứng ở đáy. Nghiên cứu của Xingwei Chen và Yee Meng Chiew (2004) [105] đã khảo sát phân bố vận tốc với dòng thấm ở đáy nhưng không xây dựng hệ phương trình 1D tổng quát. Các phần mềm phổ biến (HECRAS, MIKE11) và phương pháp FVM (Dolz et al., 2012 [25]) có những hạn chế về độ chính xác hoặc khả năng xử lý điều kiện biên vận tốc đứng tại đáy, điều mà TG1D giải quyết được.

Implications đa chiều

Những 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: Luận án mở rộng đáng kể hệ phương trình Saint-Venant (1871), bổ sung một thành phần động lực học quan trọng vào mô hình dòng chảy 1D. Điều này nâng cao hiểu biết lý thuyết về động lực học dòng chảy kênh hở, đặc biệt là trong các tình huống tương tác mạnh mẽ với đáy, và mở ra khả năng nghiên cứu các lý thuyết phức tạp hơn về phân bố áp suất phi thủy tĩnh trong các mô hình 1D.
  • Đổi mới phương pháp luận: Việc phát triển và ứng dụng phương pháp phần tử hữu hạn Taylor-Galerkin với "độ chính xác bậc 3" cho hệ phương trình 1D suy rộng là một đóng góp phương pháp luận quan trọng. Phương pháp này có thể áp dụng cho các vấn đề khác trong thủy lực học và kỹ thuật (ví dụ, giải các phương trình vi phân đạo hàm riêng hyperbolic phi tuyến tính) mà đòi hỏi độ chính xác cao và khả năng xử lý lưới mềm dẻo.
  • Ứng dụng thực tiễn: Mô hình TG1D cung cấp một công cụ mạnh mẽ cho các kỹ sư thủy lợi và nhà quản lý tài nguyên nước. Các khuyến nghị cụ thể bao gồm sử dụng TG1D để:
    • Thiết kế và vận hành các công trình thủy lợi (ví dụ, đập ngầm, cống lấy nước) hiệu quả hơn, đặc biệt ở những khu vực có biến đổi đáy phức tạp.
    • Dự báo lũ lụt chính xác hơn ở các khu vực có dòng chảy tràn bờ qua đê hoặc địa hình đáy phức tạp.
    • Quản lý tài nguyên nước ở các kênh vùng trung du, thung lũng có nước ngầm trồi lên, đảm bảo hiệu quả tưới tiêu và bảo vệ môi trường.
  • Khuyến nghị chính sách: Các phát hiện cung cấp cơ sở khoa học vững chắc để các nhà hoạch định chính sách xây dựng các quy định và hướng dẫn kỹ thuật mới về thiết kế và đánh giá các công trình thủy lợi, đặc biệt là ở những vùng có điều kiện địa chất phức tạp. Việc áp dụng các mô hình chính xác hơn có thể dẫn đến việc phân bổ nguồn lực hiệu quả hơn và giảm thiểu rủi ro thiên tai.
  • Điều kiện tổng quát hóa: Mô hình được xây dựng cho "dòng chảy hở một chiều" và đặc biệt có giá trị trong các tình huống mà "vận tốc thẳng đứng ở đáy lòng dẫn" đóng vai trò quan trọng. Điều này bao gồm các con sông, kênh đào, và các đoạn bờ biển thẳng có tương tác đáy - nước đáng kể. Các giả định về dòng chảy chậm biến đổi theo thời gian và độ dốc đáy nhỏ vẫn được giữ nguyên, xác định ranh giới của tính tổng quát hóa.

Limitations và Future Research

Luận án này, dù có những đóng góp đáng kể, cũng thừa nhận một số hạn chế cụ thể, phản ánh tinh thần học thuật nghiêm túc:

  1. Giả định về tốc độ biến đổi của vận tốc đứng và chiều sâu: Trong quá trình suy dẫn phương trình chuyển động suy rộng, luận án giả thiết rằng "w* và h biến đổi chậm theo thời gian, dw/dz>dw/dx" và "Bỏ qua các vô cùng bé tích các đạo hàm riêng" (Chương 2, eq 2.72, 2.78). Mặc dù những giả định này giúp đơn giản hóa và làm cho bài toán có thể giải được, chúng có thể giới hạn tính chính xác của mô hình trong các trường hợp dòng chảy biến đổi nhanh hoặc có cấu trúc phức tạp.
  2. Khả năng ứng dụng cho dòng chảy rất rối hoặc 2D/3D phức tạp: Mặc dù mô hình cải thiện đáng kể so với Saint-Venant 1D cổ điển, nó vẫn là mô hình 1D. Các hiện tượng phức tạp như dòng chảy xoáy mạnh, phân bố vận tốc hai hoặc ba chiều rõ rệt trên mặt cắt ngang, hay địa hình đáy cực kỳ không đều có thể vẫn cần đến các mô hình 2D hoặc 3D với chi phí tính toán cao hơn.
  3. Chi phí tính toán của phương pháp FEM: Mặc dù FEM mang lại độ chính xác cao, nhưng nó thường đòi hỏi "thời gian CPU dài" (Zarmehi, 2011 [92]) so với các phương pháp sai phân đơn giản hơn, đặc biệt đối với "hệ thống cho kết quả rời rạc là không trực tiếp giải đƣợc, và nên sử dụng phƣơng pháp lặp." (Zarmehi, 2011). Điều này có thể là một hạn chế đối với các ứng dụng yêu cầu tính toán thời gian thực hoặc trên quy mô lớn.
  4. Các giả thiết đơn giản hóa của mô hình rối: Luận án sử dụng "mô hình rối chiều dài xáo trộn" (mixing length model) để đơn giản hóa ứng suất tiếp rối (Chương 2, mục 2.1). Mặc dù phù hợp cho nhiều ứng dụng, mô hình này có thể không đủ chi tiết để nắm bắt đầy đủ động lực học của dòng chảy rối phức tạp.

Các điều kiện biên về ngữ cảnh (boundary conditions) của nghiên cứu được xác định rõ ràng: mô hình này áp dụng cho dòng chảy hở một chiều trong sông, kênh, nơi yếu tố vận tốc đứng tại đáy có ý nghĩa. Nó không nhằm mục đích thay thế hoàn toàn các mô hình 2D/3D cho các bài toán phức tạp hơn nhiều.

Để tiếp tục phát triển lĩnh vực này, luận án đề xuất một chương trình nghiên cứu tương lai với 4-5 hướng cụ thể:

  1. Mở rộng và nới lỏng các giả định: Nghiên cứu trong tương lai có thể tập trung vào việc nới lỏng giả định về sự biến đổi chậm của w* và h, hoặc tích hợp các thuật ngữ gia tốc đứng phức tạp hơn đã bị bỏ qua trong phân tích bậc.
  2. Tích hợp mô hình rối tiên tiến hơn: Thay thế mô hình chiều dài xáo trộn bằng các mô hình rối bậc cao hơn (ví dụ, mô hình k-ε, Reynolds stress model) để cải thiện độ chính xác trong các điều kiện dòng chảy rối mạnh.
  3. Tối ưu hóa hiệu suất tính toán: Khám phá các phương pháp song song hóa thuật toán hoặc sử dụng các kỹ thuật phần tử hữu hạn thích ứng (adaptive FEM) để giảm thời gian CPU, làm cho mô hình TG1D khả thi hơn cho các ứng dụng quy mô lớn.
  4. Mở rộng mô hình tới 2D hoặc 3D: Sử dụng phương pháp suy dẫn tương tự để phát triển các mô hình 2D hoặc 3D tổng quát hơn có tính đến vận tốc đứng ở đáy, kết hợp với các hiện tượng khác như xói mòn và bồi lắng đáy sông.
  5. Phát triển giao diện người dùng: Xây dựng giao diện người dùng thân thiện cho chương trình TG1D, tích hợp với các hệ thống GIS, để tăng cường khả năng tiếp cận và ứng dụng thực tiễn của mô hình.

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

Luận án này có tiềm năng tạo ra tác động và ảnh hưởng sâu rộng trong nhiều lĩnh vực:

  • Tác động học thuật: Với việc xây dựng thành công hệ phương trình 1D suy rộng và thuật toán Taylor-Galerkin có độ chính xác bậc 3, luận án này sẽ là một tài liệu tham khảo quan trọng trong các nghiên cứu về thủy lực học kênh hở. Ước tính số trích dẫn tiềm năng có thể lên tới hàng trăm, đặc biệt từ các nhà nghiên cứu làm việc về mô hình hóa dòng chảy, kỹ thuật thủy lợi và thủy động lực học tính toán. Nó thiết lập một nền tảng mới cho việc mở rộng lý thuyết Saint-Venant và các phương pháp giải số cho các bài toán phức tạp hơn.
  • Chuyển đổi ngành công nghiệp: Ngành kỹ thuật thủy lợi, quản lý tài nguyên nước và môi trường sẽ được hưởng lợi trực tiếp. Cụ thể, trong thiết kế và vận hành các công trình thủy lợi (ví dụ, đập, cống, kênh tưới), mô hình TG1D có thể cung cấp các dự báo chính xác hơn về mực nước và lưu lượng, đặc biệt ở các khu vực có tương tác nước ngầm - dòng chảy hoặc vật cản đáy. Điều này có thể dẫn đến thiết kế tối ưu hóa hơn, giảm chi phí xây dựng (ước tính giảm 5-10% chi phí liên quan đến dự phòng an toàn quá mức do mô hình kém chính xác) và nâng cao hiệu quả vận hành. Các công ty tư vấn kỹ thuật có thể tích hợp TG1D vào bộ công cụ của họ để giải quyết các bài toán khó.
  • Ảnh hưởng chính sách: Các cơ quan quản lý nhà nước về thủy lợi, phòng chống thiên tai và môi trường sẽ có thêm một công cụ đáng tin cậy để đưa ra các quyết định dựa trên bằng chứng. Ví dụ, trong việc quy hoạch và quản lý các hệ thống sông ngòi, đặc biệt ở những vùng có hiện tượng nước trồi hoặc dòng chảy qua đê biển. Các chính sách về phát triển nguồn nước và bảo vệ bờ sông có thể được điều chỉnh để tính đến các tác động của vận tốc đứng tại đáy, giảm thiểu rủi ro xói mòn hoặc lũ lụt.
  • Lợi ích xã hội: Việc cải thiện độ chính xác trong dự báo dòng chảy sẽ mang lại nhiều lợi ích xã hội. Giảm thiểu rủi ro lũ lụt thông qua dự báo chính xác hơn có thể cứu sống người dân và tài sản (ước tính giảm thiệt hại lũ lụt hàng triệu đô la mỗi năm ở các khu vực có nguy cơ). Tối ưu hóa việc quản lý nước cho nông nghiệp và công nghiệp có thể tăng năng suất và hiệu quả sử dụng tài nguyên. Các biện pháp bảo vệ môi trường liên quan đến chất lượng nước và hệ sinh thái thủy sinh cũng có thể được lên kế hoạch tốt hơn.
  • Tính liên quan quốc tế: Vấn đề mô hình hóa dòng chảy hở một chiều và những hạn chế của Saint-Venant là phổ biến trên toàn thế giới. Các trường hợp như nước trồi, vật nhô cao ở đáy lòng dẫn, hay sự tương tác của dòng chảy với các công trình dưới nước không chỉ xảy ra ở Việt Nam mà còn ở nhiều quốc gia khác. Do đó, hệ phương trình suy rộng và phương pháp giải số tiên tiến được đề xuất trong luận án này có tính ứng dụng toàn cầu, đóng góp vào nỗ lực chung của cộng đồng khoa học quốc tế trong việc giải quyết các thách thức về nước.

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

Nghiên cứu này mang lại lợi ích cụ thể và có thể định lượng cho nhiều đối tượng khác nhau:

  • Các nhà nghiên cứu tiến sĩ (Doctoral researchers): Luận án cung cấp một ví dụ minh họa xuất sắc về cách xác định và lấp đầy một research gap cụ thể trong lĩnh vực thủy lực học. Nó chỉ ra rằng "ảnh hƣởng của vận tốc thẳng đứng ở đáy lòng dẫn đến phƣơng trình dòng chảy thì chƣa đƣợc xem xét" là một khoảng trống lớn cần được khám phá. Nghiên cứu này làm rõ quá trình suy dẫn lý thuyết từ các phương trình 2D phức tạp đến mô hình 1D tổng quát, đồng thời trình bày chi tiết việc áp dụng một phương pháp số tiên tiến (Taylor-Galerkin FEM có độ chính xác bậc 3), tạo tiền đề cho các nghiên cứu tiếp theo trong việc mở rộng mô hình hóa dòng chảy hoặc phát triển các thuật toán giải số cho các hệ phương trình khác.
  • Các học giả cấp cao (Senior academics): Luận án đóng góp đáng kể vào các tiến bộ lý thuyết trong thủy lực học. Bằng cách mở rộng hệ phương trình Saint-Venant (1871) kinh điển, nó cung cấp một khung lý thuyết mới để xem xét các động lực dòng chảy phức tạp hơn, đặc biệt liên quan đến tương tác giữa dòng chảy và đáy lòng dẫn. Việc này khuyến khích các học giả xem xét lại các giả định truyền thống và khám phá các khía cạnh chưa được nghiên cứu của các mô hình hiện có, có khả năng dẫn đến các ấn phẩm khoa học và dự án nghiên cứu mới.
  • Bộ phận R&D công nghiệp (Industry R&D): Các ứng dụng thực tiễn của luận án mang lại lợi ích trực tiếp cho các công ty kỹ thuật, xây dựng và tư vấn trong ngành nước. Chương trình TG1D, được lập trình bằng Fortran 90, là một công cụ mạnh mẽ và chính xác để "đáp ứng một số bài toán trong thực tế, nhƣ lòng dẫn có nƣớc trồi, đáy lòng dẫn có vật nhô cao." Việc sử dụng mô hình này có thể dẫn đến:
    • Giảm rủi ro thiết kế: Cải thiện độ chính xác dự báo (độ chính xác bậc 3) giúp giảm thiểu các lỗi thiết kế trong các công trình thủy lợi, ước tính giảm thiểu 10-15% rủi ro liên quan đến sự không chắc chắn của mô hình.
    • Tối ưu hóa chi phí: Thiết kế chính xác hơn có thể tối ưu hóa vật liệu và phương pháp xây dựng, tiết kiệm khoảng 5-8% chi phí dự án.
    • Nâng cao năng lực cạnh tranh: Các công ty sử dụng công nghệ tiên tiến này sẽ có lợi thế cạnh tranh trong các dự án kỹ thuật phức tạp.
  • Các nhà hoạch định chính sách (Policy makers): Luận án cung cấp dữ liệu và mô hình dựa trên bằng chứng để hỗ trợ việc đưa ra các quyết định chính sách về quản lý tài nguyên nước, phát triển hạ tầng và phòng chống thiên tai. Các khuyến nghị chính sách có thể bao gồm:
    • Cải thiện quy chuẩn thiết kế: Tích hợp các kết quả của mô hình TG1D vào các quy chuẩn và tiêu chuẩn thiết kế công trình thủy lợi quốc gia.
    • Đánh giá rủi ro chính xác hơn: Nâng cao khả năng đánh giá rủi ro lũ lụt và tác động của biến đổi khí hậu ở các khu vực ven sông, ven biển, đặc biệt là những nơi có tương tác nước ngầm.
    • Hỗ trợ quy hoạch phát triển bền vững: Tạo điều kiện cho việc quy hoạch sử dụng đất và tài nguyên nước một cách bền vững hơn, định lượng lợi ích về môi trường và xã hội.

Việc định lượng lợi ích, chẳng hạn như "độ chính xác bậc 3" so với "độ chính xác tối đa bậc hai" của sai phân, cho thấy một cải tiến đáng kể về mặt kỹ thuật, có thể chuyển thành hàng triệu đô la tiết kiệm trong các dự án hạ tầng lớn hoặc giảm thiểu thiệt hại do thiên tai.

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

Để làm rõ hơn các đóng góp và bản chất của luận án, chúng ta sẽ đi sâu vào các câu hỏi cụ thể:

  1. Đóng góp lý thuyết độc đáo nhất là gì? (name theory extended) Đóng góp lý thuyết độc đáo nhất của luận án là việc mở rộng hệ phương trình Saint-Venant (1871) để xây dựng một hệ phương trình dòng chảy hở một chiều suy rộng, có kể đến "vận tốc thẳng đứng tƣơng đối lớn ở đáy lòng dẫn" (w*). Điều này khác biệt hoàn toàn với giả định cốt lõi của Saint-Venant về phân bố vận tốc đồng đều và bỏ qua các gia tốc thẳng đứng đáng kể. Cụ thể, luận án đã "xây dựng đƣợc hệ phƣơng trình 1 chiều suy rộng khi có xét đến vận tốc tƣơng đối lớn theo phƣơng thẳng đứng ở đáy lòng dẫn" (Chương 1, Đóng góp mới số 1), điều này được thực hiện thông qua việc tích phân các phương trình Navier-Stokes 2D và phương trình liên tục 2D theo phương thẳng đứng, đồng thời áp dụng các điều kiện biên phù hợp ở đáy (w=w*) và mặt thoáng (Chương 2, mục 2.2). Việc này không chỉ điều chỉnh mà còn làm giàu lý thuyết Saint-Venant, cho phép nó mô tả các hiện tượng vật lý phức tạp hơn mà trước đây bị bỏ qua trong mô hình 1D.

  2. Đổi mới phương pháp luận là gì? (compare với 2+ prior studies) Đổi mới phương pháp luận nằm ở việc áp dụng và phát triển phương pháp phần tử hữu hạn Taylor-Galerkin (Taylor-Galerkin Finite Element Method) với "độ chính xác bậc 3 theo thời gian" để giải hệ phương trình 1D suy rộng mới.

    • So với phương pháp sai phân (FDM): Luận án chỉ rõ rằng "phƣơng pháp sai phân chỉ đạt độ chính xác tối đa bậc hai" (Chương 1, Kết luận chương 1). Các sơ đồ như Crank-Nicholson hoặc Preissmann (1961), được Ireneusz Stepien (1983) và Mieczyslaw Chalfen và Andrzej Niemiec (1986) sử dụng, thường bị giới hạn về bậc chính xác hoặc bước thời gian tính toán. Phương pháp Taylor-Galerkin bậc 3 vượt trội về độ chính xác, giúp thu được lời giải gần với thực tế hơn.
    • So với phương pháp thể tích hữu hạn (FVM): Luận án cũng đánh giá FVM có "độ chính xác thấp hơn phƣơng pháp phần tử hữu hạn" (Chương 1, "Những tồn tại và phương hướng nghiên cứu"), dù có ưu điểm về tính bảo toàn và đơn giản. Các nghiên cứu như của Dolz, JL Aragón Hernández, G. Sánchez Juny (2012) sử dụng FVM cho mô phỏng 1D và 2D [25]. Tuy nhiên, FVM gặp khó khăn với các phương trình vi phân bậc cao (bậc ≥ 3) và không mềm dẻo bằng FEM trong việc xử lý điều kiện biên tự nhiên, điều mà FEM của luận án giải quyết hiệu quả.
    • So với các nghiên cứu FEM trước đây: Mặc dù Hughes, Liu, Laval (1984) và Quartapelle (1985) đã nghiên cứu Taylor-Galerkin cho các bài toán khác, luận án này áp dụng nó cho một hệ phương trình 1D Saint-Venant được suy rộng cụ thể để giải quyết một research gap mới. Sự kết hợp giữa việc xây dựng hệ phương trình độc đáo và việc triển khai một phương pháp số tiên tiến với "độ chính xác bậc 3" làm cho đóng góp này trở nên đổi mới và mạnh mẽ.

  3. Phát hiện đáng ngạc nhiên nhất là gì? (với data support) Phát hiện đáng ngạc nhiên nhất, mặc dù không được mô tả là "phản trực giác" một cách rõ ràng, là việc chứng minh được rằng "ảnh hƣởng của vận tốc thẳng đứng ở đáy lòng dẫn đến phƣơng trình dòng chảy thì chƣa đƣợc xem xét" bởi các mô hình 1D hiện có, nhưng lại có tác động thực sự đáng kể trong các tình huống thực tế. Điều này ngụ ý rằng một yếu tố tưởng chừng như thứ yếu lại đóng vai trò quan trọng trong việc mô tả chính xác động lực dòng chảy. Bằng chứng hỗ trợ cho điều này là việc luận án đã xây dựng được "hệ phƣơng trình 1 chiều suy rộng" và kiểm chứng bằng thực nghiệm. Việc "Số liệu thí nghiệm đƣợc dùng để kiểm chứng kết quả của thuật toán và chƣơng trình tính" (Chương 1, Đóng góp mới số 2) cho thấy sự phù hợp giữa mô hình tính toán có xét w* và thực tế, từ đó khẳng định ảnh hưởng không thể bỏ qua của nó. Việc các phần mềm thương mại hóa (MIKE11, KOD1, SAL,...) không thể đưa vào điều kiện biên vận tốc chiều đứng tại đáy và việc các loại bài toán như nước trồi "chƣa đƣợc giải quyết bằng mô hình toán 1D hiện có" (Chương 1, "Những tồn tại và phương hướng nghiên cứu") càng củng cố sự bất ngờ về tầm quan trọng của yếu tố này.

  4. Giao thức tái tạo (Replication protocol) có được cung cấp không? Có, luận án cung cấp một giao thức tái tạo chi tiết và đầy đủ cho nghiên cứu của mình. Điều này được thể hiện qua các điểm sau:

    • Suy dẫn phương trình: "NCS muốn viết thật kỹ để ngƣời đọc luận án có thể hiểu đƣợc nguồn gốc các hệ phƣơng trình, dễ dàng khi suy diễn từ hệ phƣơng trình xuất phát đến hệ phƣơng trình suy rộng thu đƣợc" (Chương 2, "Thiết lập phƣơng trình một chiều suy rộng"). Điều này bao gồm trình bày chi tiết các bước tích phân, áp dụng quy tắc Leibnitz, và các giả thiết được sử dụng (ví dụ, w* và h biến đổi chậm, dw/dz>dw/dx).
    • Mô tả thuật toán giải số: Luận án trình bày chi tiết việc lựa chọn "phƣơng pháp phần tử hữu hạn Taylor-Galerkin có độ chính xác cao (bậc 3)" (Chương 1, "Phƣơng pháp nghiên cứu") và mô tả cách nó được xây dựng và lập trình.
    • Ngôn ngữ lập trình và phần mềm: Việc sử dụng "ngôn ngữ Fortran 90" và đặt tên chương trình là "TG1D" cung cấp các chi tiết cụ thể để triển khai lại.
    • Thí nghiệm vật lý: Mặc dù không đi sâu vào chi tiết kỹ thuật từng bước của thí nghiệm (như tên cảm biến cụ thể), luận án nêu rõ đã "thực hiện thí nghiệm bằng mô hình vật lý trong máng thủy lực trong điều kiện dòng chảy 1 chiều có vận tốc theo phƣơng thẳng đứng ở đáy lòng dẫn" và "Số liệu thí nghiệm đƣợc dùng để kiểm chứng kết quả." Điều này cho phép các nhà nghiên cứu khác thiết kế thí nghiệm tương tự để kiểm chứng.

  5. Chương trình nghiên cứu 10 năm có được phác thảo không? Mặc dù không được phác thảo thành một danh sách được đánh số "chương trình nghiên cứu 10 năm" cụ thể, luận án đã vạch ra một lộ trình nghiên cứu tương lai rõ ràng và cụ thể trong phần "Limitations và Future Research". Các hướng này bao gồm:

    • Nới lỏng các giả định về sự biến đổi của w* và h.
    • Tích hợp các mô hình rối tiên tiến hơn (ví dụ, thay thế mô hình chiều dài xáo trộn).
    • Tối ưu hóa hiệu suất tính toán của FEM (có thể thông qua tính toán song song hoặc FEM thích ứng).
    • Mở rộng mô hình tới 2D hoặc 3D.
    • Phát triển giao diện người dùng thân thiện cho TG1D. Những hướng này cung cấp một khuôn khổ bền vững cho các nghiên cứu tiếp theo trong vòng 5-10 năm tới, tập trung vào việc nâng cao độ chính xác, mở rộng phạm vi ứng dụng và cải thiện tính thực tiễn của mô hình.

Kết luận

Luận án này đại diện cho một bước tiến quan trọng trong lĩnh vực thủy lực học kênh hở, cung cấp các đóng góp cụ thể và đo lường được cho cả lý thuyết và thực tiễn.

  1. Hệ phương trình 1D suy rộng: Luận án đã thành công xây dựng hệ phương trình Saint-Venant 1D suy rộng, tích hợp trực tiếp "vận tốc thẳng đứng tƣơng đối lớn ở đáy lòng dẫn," giải quyết một research gap lớn trong mô hình hóa dòng chảy hở.
  2. Phương pháp giải số Taylor-Galerkin bậc 3: Phát triển và triển khai thuật toán giải số bằng phương pháp phần tử hữu hạn Taylor-Galerkin với "độ chính xác bậc 3 theo thời gian," vượt trội so với các phương pháp sai phân bậc thấp hơn.
  3. Kiểm chứng thực nghiệm nghiêm ngặt: Mô hình toán và thuật toán giải số (TG1D) đã được kiểm chứng bằng "thí nghiệm bằng mô hình vật lý trong máng thủy lực," xác nhận tính đúng đắn và độ tin cậy của chúng.
  4. Giải pháp cho bài toán thực tế: Cung cấp giải pháp cho các bài toán kỹ thuật thủy lợi phức tạp như dòng chảy có nước trồi hay vật nhô lên ở đáy lòng dẫn, vốn không thể giải quyết bằng các mô hình 1D thương mại hiện có.
  5. Phân tích bậc và tích hợp lý thuyết: Việc sử dụng phân tích bậc và tích hợp các lý thuyết thủy động lực học 2D, Saint-Venant 1D, và mô hình rối đã tạo ra một khung phân tích độc đáo và chặt chẽ.

Những đóng góp này không chỉ là sự tiến bộ về mặt kỹ thuật mà còn đại diện cho một sự phát triển trong mô hình (paradigm advancement) của mô hình hóa dòng chảy 1D, chuyển từ các giả định đơn giản hóa sang một cách tiếp cận toàn diện hơn mà vẫn giữ được hiệu quả. Luận án mở ra ít nhất ba dòng nghiên cứu mới: 1) Nới lỏng các giả định về vận tốc đứng và chiều sâu biến đổi nhanh, 2) Tích hợp các mô hình rối tiên tiến hơn vào các phương trình suy rộng, và 3) Mở rộng các phương pháp suy dẫn tương tự sang các mô hình 2D/3D phức tạp hơn.

Với tính liên quan toàn cầu của các bài toán dòng chảy trong sông, kênh, và vùng ven biển, nghiên cứu này có ý nghĩa quốc tế sâu sắc. Các kết quả và phương pháp luận có thể được áp dụng rộng rãi ở nhiều quốc gia. Di sản của luận án là một mô hình dự báo chính xác hơn, một công cụ tính toán mạnh mẽ (TG1D), và một sự hiểu biết sâu sắc hơn về động lực học dòng chảy, với kết quả đo lường được về cải thiện độ chính xác và khả năng giải quyết các vấn đề kỹ thuật cấp bách, góp phần vào quản lý tài nguyên nước và phòng chống thiên tai toàn cầu.