Tổng quan về luận án

Mô hình hóa dòng chảy một chiều (1D) trong hệ thống kênh, sông đóng vai trò cốt lõi trong quy hoạch phát triển nguồn nước, phòng chống thiên tai và quản lý thủy lực công trình. Từ thế kỷ XIX, nền tảng của thủy lực học sông ngòi chủ yếu dựa trên hệ phương trình Saint-Venant (1871) với các giả thiết kinh điển: dòng chảy thay đổi dần, vận tốc phân bố đều trên mặt cắt ngang, đường dòng có độ cong nhỏ và gia tốc thẳng đứng triệt tiêu nhằm thỏa mãn quy luật phân bố áp suất thủy tĩnh. Tuy nhiên, các bài toán thủy động lực học hiện đại đang đối mặt với những hiện tượng dị thường tại đáy lòng dẫn: dòng chảy đi qua vùng nước trồi địa chất, mạch nước ngầm có áp phun trào từ đáy (như kênh dẫn qua vùng trung du, thung lũng karst hoặc kênh qua lòng hồ Phú Ninh), sự hiện diện của đập ngầm nhân tạo nâng cao trình mực nước, hay hiện tượng sóng tràn qua bờ bao và đê biển. Tại các khu vực này, gradient vận tốc theo phương thẳng đứng ở đáy lòng dẫn $w^*$ xuất hiện với cường độ lớn, trực tiếp phá vỡ giả thiết phân bố áp suất thủy tĩnh kinh điển.

Khoảng trống học thuật cốt lõi (Research Gap) nằm ở chỗ: các mô hình 1D thương mại hiện hành như HEC-RAS, MIKE 11, VRSAP (Nguyễn Như Khuê) hay SAL (Nguyễn Tất Đắc) chỉ hỗ trợ nguồn bổ sung bên (lateral flow $q$) với vận tốc nhỏ, hoàn toàn bất lực trong việc mô tả thành phần động lượng thẳng đứng tại đáy. Ngược lại, các mô hình số trị hai chiều thẳng đứng (2DV) hoặc ba chiều (3D) như ANSYS Fluent, TELEMAC-3D đòi hỏi tài nguyên tính toán khổng lồ và thường sử dụng phương pháp thể tích hữu hạn (FVM) bậc thấp. Luận án giải quyết triệt để bài toán này thông qua mục tiêu kép: thiết lập hệ phương trình vi phân 1D suy rộng tích hợp vận tốc thẳng đứng đáy $w^*$ và giải số bằng phương pháp phần tử hữu hạn Taylor-Galerkin (TG-FEM) bậc 3 trên nền tảng ngôn ngữ Fortran 90.

Các câu hỏi và giả thuyết nghiên cứu được xác lập tường minh:

  • RQ1: Làm thế nào để tích phân hệ phương trình vi phân chuyển động Navier-Stokes hai chiều đứng theo độ sâu dòng chảy khi điều kiện biên tại đáy tồn tại vận tốc thẳng đứng $w^* \neq 0$ mà không làm mất đi các thành phần phi thủy tĩnh?
  • RQ2: Thuật toán số nào đảm bảo độ chính xác bậc cao theo thời gian, triệt tiêu dao động số giả và bảo toàn động lượng khi giải hệ phương trình 1D suy rộng phi tuyến?
  • H1: Sự xuất hiện của vận tốc thẳng đứng đáy $w^*$ gây ra sự tái phân bố áp suất phi thủy tĩnh, làm gia tăng tổn thất năng lượng cục bộ và tạo hiệu ứng dâng nước thượng lưu vượt trội so với dòng nhập lưu bên thông thường.
  • H2: Sơ đồ phần tử hữu hạn Taylor-Galerkin bậc 3 (TG-FEM) lập trình trên Fortran 90 cho phép mô phỏng chính xác sóng phi tuyến và cấu trúc dòng chảy đáy với độ ổn định cao hơn các sơ đồ sai phân hữu hạn bậc 2 truyền thống.

Khung lý thuyết của công trình tích hợp giải tích vi tích phân Leibniz, lý thuyết lớp biên rối Prandtl với mô hình chiều dài xáo trộn, và phương pháp phần tử hữu hạn Taylor-Galerkin. Luận án giới hạn phạm vi nghiên cứu trong dòng chảy hở 1D không nén được, kiểm chứng đối chuẩn bằng số liệu thực nghiệm phân bố lưu tốc của Xingwei Chen và Yee Meng Chiew (2004) kết hợp thí nghiệm mô hình vật lý thủy lực máng hở tự tạo (mã hiệu chương trình tính TG1D).

Literature Review và Positioning

Lịch sử phát triển của thủy lực học dòng chảy hở ghi nhận sự tiến hóa từ các mô hình giải tích đơn giản đến mô phỏng số phức tạp. Sau khi Saint-Venant (1871) công bố hệ phương trình vi phân mô tả chuyển động không ổn định thay đổi chậm, Bakhmeteff (1932) đã hệ thống hóa phân loại đường mặt nước cho dòng chảy ổn định không đều dựa trên độ dốc đáy, độ dốc phân giới và số Froude ($Fr$). Đến năm 1951, Hayami (Nhật Bản) đề xuất mô hình xấp xỉ sóng khuếch tán với hệ số khuếch tán $D_f = Q / (2 B S_0)$, đặt nền móng cho các phương pháp diễn toán lũ một chiều giản lược. Cuối thập niên 1960 và 1970, Chow (1971), Strelkoff (1969) và Ben Chie Yen (1973, 1984) đã nỗ lực thu hẹp khoảng cách giữa cơ học chất lưu và thủy lực kênh hở thông qua tích phân phương trình Navier-Stokes dạng điểm theo mặt cắt ướt, làm rõ sự khác biệt bản chất giữa phương trình động lượng và phương trình năng lượng.

Trường phái phương pháp số giải phương trình sóng nước nông ghi nhận sự phân cực sâu sắc giữa hai nhánh kỹ thuật:

  1. Phương pháp sai phân hữu hạn (FDM): Preissmann (1961) giới thiệu sơ đồ ẩn 4 điểm kinh điển được Cunge (1981), Abbott và Ionescu (1967), Amein và Fang chuẩn hóa. Mặc dù các sơ đồ hiện như Crank-Nicholson hay MacCormack (1969) với kỹ thuật dự báo - hiệu chỉnh mang lại sự thuận tiện trong lập trình, chúng bị khống chế nghiêm ngặt bởi điều kiện ổn định Courant-Friedrichs-Lewy (CFL) và chỉ đạt độ chính xác tối đa bậc hai ($O(\Delta t^2, \Delta x^2)$). Khi áp dụng cho dòng chảy có biến đổi đột ngột ở đáy, FDM xuất hiện hiện tượng khuếch tán số hoặc dao động sóng ngắn không vật lý kiểu '$2\Delta x$'.
  2. Phương pháp phần tử hữu hạn (FEM) và Taylor-Galerkin: Khắc phục nhược điểm của FDM, Donea (1984), Thomas J.R. Hughes et al. (1984) và L. Quartapelle (1985, 1997, 1998) đã phát triển sơ đồ Taylor-Galerkin cho các bài toán hyperbol phi tuyến, đạt độ chính xác bậc 3 theo thời gian và triệt tiêu dao động số mà không cần bổ sung hệ số nhớt nhân tạo. Nwaogazie (1985, 1987) đã ứng dụng FEM Galerkin kết hợp lặp Newton-Raphson vào phương trình Saint-Venant, chứng minh ưu thế vượt trội về việc bám sát điều kiện biên hình học phức tạp. Szymkiewicz (1991) và Fatemeh Zarmehi et al. (2011) tiếp tục hoàn thiện tính ổn định của FEM trên mạng lưới kênh hở.
       TIẾN TRÌNH PHÁT TRIỂN LÝ THUYẾT DÒNG CHẢY HỞ VÀ PHƯƠNG PHÁP SỐ
       
  Saint-Venant (1871) ───────> Chow (1971), Yen (1973) ───────> Luận án (Mô hình 1D suy rộng)
  [1D cổ điển, p thủy tĩnh]     [Tích phân Navier-Stokes 2D]      [Kể đến w* đáy, p phi thủy tĩnh]
            │                                                               │
            ▼                                                               ▼
  Preissmann (1961), FDM ────> Hughes, Quartapelle (1985) ────> TG1D (Fortran 90)
  [Sai phân ẩn bậc 2]          [Taylor-Galerkin FEM bậc 3]        [TG-FEM bậc 3 + Thực nghiệm máng]

Vị thế học thuật của luận án được xác định rõ nét khi đặt cạnh hai nghiên cứu quốc tế tiêu biểu:

  • So sánh với Xingwei Chen và Yee Meng Chiew (2004): Nhóm tác giả Chen & Chiew đã thực hiện thí nghiệm đo đạc phân bố lưu tốc của dòng chảy rối có thấm hút qua đáy (bed suction) và điều chỉnh quy luật phân bố logarit. Tuy nhiên, họ dừng lại ở việc khảo sát thực nghiệm phân bố vận tốc điểm mà không thiết lập hệ phương trình vi phân tích phân 1D tổng quát dưới ảnh hưởng của vận tốc đáy. Luận án đã kế thừa tập dữ liệu thực nghiệm của Chen & Chiew để thực hiện phân tích bậc đại số, đồng thời mở rộng bài toán sang trường hợp dòng nước trồi/phun ngược từ đáy.
  • So sánh với Xavier Litrico và Vincent Fromion (2005, 2006): Litrico & Fromion tập trung vào việc tuyến tính hóa phương trình Saint-Venant cổ điển để thiết kế thuật toán điều khiển tự động hồ kênh thủy lợi tại Bồ Đào Nha. Nghiên cứu của họ vẫn dựa trên giả thiết phân bố áp suất thủy tĩnh thuần túy, do đó không thể áp dụng cho các đoạn kênh có cấu trúc đáy dị thường hoặc có dòng phun áp lực như mô hình 1D suy rộng của luận án.

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

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

Luận án đã tạo ra bước tiến quan trọng trong thủy lực học lý thuyết bằng cách mở rộng hệ phương trình Saint-Venant (1871), loại bỏ rào cản giả thiết áp suất thủy tĩnh thông qua việc tích hợp trực tiếp thành phần vận tốc thẳng đứng ở đáy $w^*$:

$$\text{Điều kiện biên đáy: } z = 0 \implies w = w^*, \quad u_0 = 0$$

$$\text{Điều kiện biên mặt thoáng: } z = h \implies w_m = \frac{\partial h}{\partial t} + u_m \frac{\partial h}{\partial x}$$

Xuất phát từ hệ phương trình vi phân Navier-Stokes hai chiều đứng đối với chất lưu nhớt không nén được, luận án thực hiện tích phân dọc theo chiều sâu từ đáy ($z=0$) đến mặt thoáng ($z=h$). Bằng việc áp dụng quy tắc tích phân Leibniz:

$$\frac{\partial}{\partial x} \int_{0}^{h(x,t)} u(x,z,t) , dz = \int_{0}^{h} \frac{\partial u}{\partial x} , dz + u(x,h,t)\frac{\partial h}{\partial x} - u(x,0,t)\frac{\partial z_0}{\partial x}$$

Biểu thức áp suất phi thủy tĩnh $p(x,z,t)$ được phân tích và tích phân đầy đủ, bao gồm các số hạng động áp sinh ra do gia tốc thẳng đứng của dòng chất lưu:

$$p(x,z,t) = \rho g (h - z) + \rho \int_{z}^{h} \left( \frac{\partial w}{\partial t} + u\frac{\partial w}{\partial x} + w\frac{\partial w}{\partial z} - \nu \nabla^2 w \right) dz$$

Hệ phương trình vi phân 1D suy rộng thu được có dạng ma trận bảo toàn tổng quát:

$$\frac{\partial \mathbf{U}}{\partial t} + \frac{\partial \mathbf{F}(\mathbf{U})}{\partial x} + \mathbf{S}(\mathbf{U}, w^*) = \mathbf{0}$$

Trong đó, vectơ biến trạng thái $\mathbf{U} = [A, Q]^T$, vectơ thông lượng $\mathbf{F}(\mathbf{U})$ chứa các thành phần động lượng chuyển động dọc trục, và vectơ nguồn $\mathbf{S}(\mathbf{U}, w^)$ bổ sung các số hạng phi tuyến phụ thuộc chặt chẽ vào vận tốc đứng $w^$, đạo hàm $\partial w^*/\partial x$, và gradient mực nước $\partial h/\partial x$. Đây là một bước chuyển dịch mô hình (paradigm shift) từ thủy lực học một chiều giả định tĩnh sang thủy động lực học một chiều phi thủy tĩnh toàn phần.

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

Khung phân tích của luận án là sự hợp nhất của ba trụ cột lý thuyết:

  1. Lý thuyết thủy động lực học Navier-Stokes: Thiết lập mối liên hệ vi phân cục bộ giữa vận tốc dọc $u(x,z,t)$ và vận tốc đứng $w(x,z,t)$.
  2. Mô hình rối chiều dài xáo trộn Prandtl: Ứng suất tiếp rối được xác định thông qua $\tau_t = \rho l^2 |\partial u/\partial z|(\partial u/\partial z)$, với chiều dài xáo trộn $l = k z$. Luận án xác lập hằng số Von Kármán $k = 0.54$ cho dòng chảy sông thiên nhiên và $k = 0.40$ cho dòng chảy lòng dẫn nhân tạo, kết hợp với ứng suất cắt nhớt $\tau_\mu = \mu (\partial u/\partial z)$ để tạo nên ứng suất tiếp tổng cộng $\tau = \tau_\mu + \tau_t$.
  3. Kỹ thuật phân tích bậc đại số (Order of Magnitude Analysis): Sử dụng chuỗi số liệu thực nghiệm từ công trình của Xingwei Chen và Yee Meng Chiew (2004) tại các mặt cắt $x = -10\text{ cm}, 25\text{ cm}, 50\text{ cm}, 75\text{ cm}$, tác giả tính toán độ lớn định lượng của từng số hạng tích phân vi phân dọc theo phương đứng ($z = 0 \to 3.3\text{ cm}$). Các số hạng chứa tích vô cùng bé của các đạo hàm riêng bậc cao $(\partial w/\partial x)(\partial u/\partial z)$ được triệt tiêu một cách khoa học, giữ lại các thành phần chủ đạo chi phối áp suất và lực cản.

Điều kiện biên và giới hạn ứng dụng của khung phân tích được xác định rõ ràng: dòng chảy một chiều chiếm ưu thế dọc trục kênh ($u \gg w$), đáy lòng dẫn có độ dốc êm ($\cos\theta \approx 1$), bỏ qua lực quán tính Coriolis và ứng suất gió mặt thoáng, biến thiên của vận tốc đứng theo chiều sâu chiếm ưu thế so với biến thiên dọc trục ($\partial w/\partial z > \partial w/\partial x$).

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái thực chứng toán - lý kết hợp (positivist computational & experimental design), kết hợp suy diễn giải tích toán học với mô phỏng số độ chính xác cao và đối chứng thực nghiệm vật lý.

┌─────────────────────────────────────────────────────────────────────────┐
│                           THIẾT KẾ NGHIÊN CỨU                           │
└────────────────────────────────────┬────────────────────────────────────┘
                                     │
         ┌───────────────────────────┼───────────────────────────┐
         ▼                           ▼                           ▼
┌──────────────────┐       ┌──────────────────┐       ┌──────────────────┐
│  Toán giải tích  │       │ Mô hình số TG1D  │       │ Thí nghiệm máng  │
│  - Tích phân RANS│       │ - Fortran 90     │       │ - Máng thủy lực  │
│  - Phân tích bậc │       │ - Taylor-Galerkin│       │ - Đo đạc h(x,t)  │
│  - Leibniz rule  │       │ - FEM bậc 3      │       │ - Đối chuẩn số   │
└──────────────────┘       └──────────────────┘       └──────────────────┘

Thiết kế nghiên cứu đa tầng bao gồm ba cấp độ:

  • Cấp độ 1 (Continuous Field): Xây dựng hệ phương trình đạo hàm riêng liên tục phi tuyến.
  • Cấp độ 2 (Discrete Numerical Domain): Rời rạc hóa không gian - thời gian bằng kỹ thuật phần tử hữu hạn Taylor-Galerkin bậc 3.
  • Cấp độ 3 (Physical Laboratory Verification): Thí nghiệm mô hình vật lý trong máng thủy lực kiểm định trường dòng chảy thực tế.

Quy trình nghiên cứu rigorous

Quy trình giải thuật số trị được xây dựng qua hai giai đoạn chặt chẽ:

  1. Rời rạc hóa theo thời gian (Taylor Step): Khai triển chuỗi Taylor cho vectơ biến trạng thái $\mathbf{U}^{n+1}$ đến số hạng đạo hàm bậc 3:

$$\mathbf{U}^{n+1} = \mathbf{U}^n + \Delta t \left( \frac{\partial \mathbf{U}}{\partial t} \right)^n + \frac{\Delta t^2}{2} \left( \frac{\partial^2 \mathbf{U}}{\partial t^2} \right)^n + \frac{\Delta t^3}{6} \left( \frac{\partial^3 \mathbf{U}}{\partial t^3} \right)^n + O(\Delta t^4)$$

Thay thế các đạo hàm thời gian bằng đạo hàm không gian thông qua phương trình trạng thái liên tục, chuyển đổi bài toán sang dạng phương trình vi phân thông thường ẩn tuyến tính.

  1. Rời rạc hóa theo không gian (Galerkin Spatial Step): Miền tính toán 1D được chia thành $N$ phần tử hữu hạn hai nút. Áp dụng phép chiếu Galerkin với hệ hàm trọng số $\psi_j(x)$ trùng với hàm nội suy tuyến tính của phần tử:

$$\int_{\Omega} \left[ \mathbf{U}^{n+1} - \mathbf{U}^n + \Delta t \frac{\partial \mathbf{F}^n}{\partial x} - \frac{\Delta t^2}{2} \frac{\partial}{\partial x}\left( \mathbf{A}^n \frac{\partial \mathbf{F}^n}{\partial x} \right) + \dots \right] \psi_j(x) , dx = \mathbf{0}$$

Hệ phương trình đại số tuyến tính toàn cục sau khi ghép nối có dạng ma trận ba đường chéo (tridiagonal matrix), được giải hữu hiệu bằng thuật toán quét kép Thomas Algorithm tích hợp trong mã nguồn Fortran 90.

       QUY TRÌNH THUẬT TOÁN TAYLOR-GALERKIN BẬC 3 TRONG TG1D
       
  [U^n, h^n, Q^n] ───> [Khai triển Taylor bậc 3 theo thời gian: Δt, Δt^2, Δt^3]
                              │
                              ▼
  [Thay thế đạo hàm thời gian bằng toán tử không gian: ∂F/∂x, ∂/∂x(A ∂F/∂x)]
                              │
                              ▼
  [Phép chiếu phần tử hữu hạn Galerkin với hàm trọng ψ_j(x)]
                              │
                              ▼
  [Ma trận ba đường chéo (Tridiagonal Matrix)] ───> [Thomas Algorithm] ───> [U^{n+1}]

Độ tin cậy và tính hợp lệ của phương pháp số được bảo đảm qua việc triệt tiêu sai số làm mịn (numerical dissipation error) và sai số tán sắc (dispersion error), vượt qua phép thử điều kiện biên tiêu tán của Nguyễn Thế Hùng (1989).

Data và phân tích

Dữ liệu phân tích bao gồm hai nguồn độc lập:

  1. Dữ liệu đối chuẩn thứ cấp: Trích xuất từ ma trận đo đạc vận tốc của Xingwei Chen và Yee Meng Chiew (2004) tại các độ sâu $z = [0.0, 0.1, 0.2, \dots, 3.3]\text{ cm}$ với vận tốc dọc $u$ biến thiên từ $0.0\text{ cm/s}$ đến $45.2\text{ cm/s}$. Các tích phân hình thang dọc theo độ sâu được thực hiện nhằm định lượng chính xác trọng số của các số hạng đạo hàm $\partial u/\partial x$.
  2. Dữ liệu thực nghiệm nguyên cấp: Thu thập từ máng thủy lực phòng thí nghiệm với hệ thống tạo dòng phun đáy có lưu lượng kế kiểm soát và đầu đo mực nước tự động. Chương trình Fortran 90 (TG1D) chạy mô phỏng trên lưới không gian mịn với $\Delta x = 0.05\text{ m}$ và bước thời gian $\Delta t = 0.001\text{ s}$, đáp ứng hệ số Courant $CFL = (u + \sqrt{gh})\Delta t / \Delta x \le 0.5$.

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

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

  1. Hiệu ứng cản thủy lực bổ sung do vận tốc đứng đáy: Vận tốc thẳng đứng $w^*$ tại đáy không hoạt động như một nguồn bổ sung lưu lượng thuần túy mà tạo ra một trường gradient áp suất phi thủy tĩnh cục bộ rất lớn $\partial p/\partial x$. Lực cản này làm biến dạng đường mặt nước mạnh mẽ hơn nhiều so với cùng một lượng lưu lượng gia nhập từ hai bên bờ kênh ($q$).
  2. Sự gia tăng đột biến của cao trình dâng nước thượng lưu: Mô hình TG1D và số liệu thực nghiệm máng thủy lực chỉ ra rằng, tại khu vực có nước trồi hoặc vật nhô đáy, mực nước thượng lưu bị dâng cao hơn từ $12.5%$ đến $28.3%$ so với kết quả dự báo của hệ phương trình Saint-Venant cổ điển.
  3. Hiện tượng tái cấu trúc profile vận tốc: Khác với phân bố vận tốc dạng logarit chuẩn trong kênh hở đáy phẳng, dòng chảy có vận tốc đáy $w^*$ tạo ra sự gia tăng gradient vận tốc sát đáy $\partial u/\partial z$, kéo theo sự dịch chuyển điểm lưu tốc cực đại xuống thấp hơn, làm thay đổi căn bản hệ số động năng $\alpha$ và hệ số động lượng $\beta$ của mặt cắt.
  4. Tính ưu việt tuyệt đối của sơ đồ Taylor-Galerkin bậc 3: Kết quả tính toán từ chương trình TG1D cho thấy thuật toán triệt tiêu hoàn toàn hiện tượng dao động sóng giả răng cưa quanh các điểm gián đoạn đáy mà không cần can thiệp hệ số nhớt nhân tạo, cải thiện độ chính xác tính toán vượt trội so với sơ đồ sai phân Crank-Nicholson và Preissmann.
 SO SÁNH PROFILE ĐƯỜNG MẶT NƯỚC VÀ ÁP SUẤT DỌC THEO KÊNH CÓ NƯỚC TRỒI ĐÁY
 
 Cao trình (m)
   ▲
   │        ════════════════════ (Đường mặt nước - Mô hình suy rộng TG1D)
   │       /                    \  [Dâng cao hơn 12.5% - 28.3% do hiệu ứng w*]
   │      /   - - - - - - - - - - \ - - (Đường mặt nước - Saint-Venant cổ điển)
   │     /                         \
   │    /                           \
 0 └───┴─────────────────▲───────────┴───────────────────────────────► Khoảng cách x (m)
                       w* > 0
               (Vùng nước trồi tại đáy)

Implications đa chiều

  • Đóng góp học thuật: Mở rộng biên giới lý thuyết của cơ học thủy khí môi trường, đặt nền móng toán học cho việc xây dựng các module 1D phi thủy tĩnh tích hợp trong các bộ phần mềm mô phỏng dòng chảy sông ngòi thế hệ mới.
  • Đổi mới phương pháp luận: Chứng minh tính khả thi và hiệu năng vượt trội của phương pháp phần tử hữu hạn Taylor-Galerkin bậc 3 lập trình trên Fortran 90, mở ra hướng tiếp cận tối ưu thay thế cho các sơ đồ sai phân FDM bậc 2 vốn bị giới hạn về độ chính xác.
  • Ứng dụng công trình thực tế: Cung cấp công cụ tính toán chính xác cho các dự án kênh dẫn thủy lợi đi qua vùng địa chất phức tạp (vùng cồn cát, vùng đá vôi có nước trồi áp lực như hệ thống kênh dẫn hồ Phú Ninh), thiết kế tối ưu đập ngầm dâng nước tự chảy, và tính toán ổn định bờ bao, đê biển khi có sóng tràn.
  • Khuyến nghị chính sách: Các cơ quan quản lý đê điều và quy hoạch tài nguyên nước cần cập nhật tiêu chuẩn thiết kế công trình kênh xả, bổ sung hệ số an toàn dâng nước do biến đổi đáy và dòng ngầm thay vì chỉ dựa vào các phần mềm 1D cổ điển như HEC-RAS.

Limitations và Future Research

Mặc dù đạt được những bước tiến đột phá, nghiên cứu vẫn tồn tại một số giới hạn học thuật cần được thừa nhận minh bạch:

  • Giả thiết phân bố vận tốc 1D bình quân: Mô hình vẫn sử dụng vận tốc trung bình mặt cắt ngang để biểu diễn dòng chảy dọc trục, do đó chưa mô phỏng chi tiết cấu trúc xoáy cuộn ba chiều (secondary flows) và lực ly tâm phức tạp tại các khúc sông cong hỗn hợp như nghiên cứu của Sugiyama et al. (2006).
  • Điều kiện biến thiên thời gian chậm: Giả thiết $\partial w/\partial z > \partial w/\partial x$ và $w^*$ biến đổi chậm theo thời gian giới hạn khả năng mô phỏng các hiện tượng sóng xung kích cực nhanh (vỡ đập tức thời mức độ cao hoặc sóng thần ven bờ).
  • Đặc trưng ma sát đáy: Chưa tích hợp đầy đủ tương tác lực cản của thảm thực vật ngập nước linh hoạt có tính biến dạng theo dòng chảy như mô hình của Xiao-guang Liu & Yu-hong Zeng (2016) hoặc Preti (2018).

Chương trình nghiên cứu tương lai (Future Agenda) được định hướng qua 4 mũi nhọn:

  1. Mở rộng hệ phương trình suy rộng sang bài toán 2 chiều mặt phẳng (2DH) phi thủy tĩnh có kể đến vận tốc đáy $w^*(x,y,t)$.
  2. Ghép nối module thủy động lực học TG1D với phương trình truyền tải - khuếch tán bùn cát đáy nhằm đánh giá hiện tượng xói lở cục bộ do dòng nước trồi gây ra (Zhihua Xie et al., 2012).
  3. Tích hợp mô hình lực cản thảm thực vật đáy phi tuyến tính (drag coefficient $C_d$) vào vectơ nguồn của hệ phương trình.
  4. Phát triển giao diện đồ họa người dùng (GUI) và kết nối module TG1D với hệ thống thông tin địa lý (GIS) tương tự định hướng hiện đại hóa của mô hình KOD1 và VRSAP.

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

Nghiên cứu mang lại tiềm năng tác động sâu rộng trên nhiều bình diện:

  • Tác động học thuật: Định vị lại phương pháp luận mô hình hóa thủy lực 1D, kỳ vọng tạo ra chỉ số trích dẫn cao trong các tạp chí chuyên ngành thủy lực, mô phỏng số và cơ học chất lưu tính toán (CFD).
  • Chuyển đổi công nghệ kỹ thuật: Cung cấp mã nguồn thuật toán Fortran 90 tối ưu hóa cao cho các viện nghiên cứu quy hoạch thủy lợi, trung tâm dự báo khí tượng thủy văn, giúp nâng cao độ tin cậy của các bài toán diễn toán lũ và quy hoạch tiêu thoát nước.
  • Lợi ích kinh tế - xã hội: Giảm thiểu rủi ro tràn bờ vỡ kênh, tối ưu hóa kích thước mặt cắt công trình thủy lợi dẫn dòng qua vùng địa chất ngầm, tiết kiệm hàng chục tỷ đồng chi phí xây dựng và gia cố bờ kênh hàng năm.
  • Tầm vóc quốc tế: Đóng góp một trường hợp nghiên cứu giải tích và thuật toán số hoàn chỉnh vào kho tàng cơ học tính toán thế giới, giải quyết bài toán biên phức tạp mà các phần mềm quốc tế thương mại chưa chuẩn hóa ở cấp độ 1D.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận phương pháp tích phân giải tích phương trình Navier-Stokes qua quy tắc Leibniz và quy trình rời rạc hóa Taylor-Galerkin bậc 3 chi tiết.
  • Kỹ sư thiết kế công trình thủy lợi - thủy điện: Sở hữu công thức và phần mềm tính toán cao trình mực nước dâng chính xác khi tuyến kênh đi qua khu vực có đập ngầm hoặc mạch nước trồi.
  • Các nhà phát triển phần mềm thủy lực: Kế thừa cấu trúc thuật toán giải ma trận ba đường chéo và kỹ thuật lập trình Fortran 90 hiện đại để tích hợp vào các gói phần mềm thủy văn - thủy lực thương mại.
  • Cơ quan quản lý tài nguyên nước và môi trường: Có cơ sở khoa học xác đáng để ban hành các quy chuẩn kỹ thuật mới trong đánh giá nguy cơ ngập lụt vùng ven sông, đê biển.

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

  1. Đóng góp lý thuyết độc đáo nhất của luận án là gì? Đó là việc thiết lập thành công hệ phương trình vi phân tích phân 1D suy rộng phi thủy tĩnh đầu tiên tích hợp tường minh điều kiện biên vận tốc thẳng đứng tại đáy lòng dẫn ($w^* \neq 0$), vượt qua giới hạn lịch sử hơn 140 năm của hệ phương trình Saint-Venant (1871).

  2. Đột phá phương pháp luận số trị được thể hiện như thế nào khi so sánh với các nghiên cứu trước? Khác với các sơ đồ sai phân FDM bậc 2 truyền thống (Preissmann, Crank-Nicholson) vốn nhạy cảm với điều kiện ổn định và tạo dao động số, luận án đã ứng dụng thành công phương pháp phần tử hữu hạn Taylor-Galerkin bậc 3 theo thời gian. So với FVM trong ANSYS (vốn nặng về chi phí tính toán khi chạy 2D/3D), thuật toán TG1D đạt độ chính xác tương đương nhưng tối ưu hóa tốc độ tính toán gấp nhiều lần trên mô hình 1D.

  3. Phát hiện nào mang tính bất ngờ nhất từ kết quả mô phỏng và thực nghiệm? Hiện tượng dâng nước thượng lưu do dòng phun thẳng đứng đáy $w^*$ gây ra không tỷ lệ thuận đơn giản theo lưu lượng bổ sung như lý thuyết Saint-Venant dự báo, mà bị khuếch đại phi tuyến do gradient áp suất động lực học $\partial p/\partial x$, làm mực nước dâng cao hơn thực tế tới gần $30%$.

  4. Luận án có cung cấp đầy đủ giao thức tái lập nghiên cứu (Replication Protocol) không? Hoàn toàn đầy đủ. Toàn bộ các bước biến đổi vi tích phân từ phương trình Navier-Stokes 2D, quy tắc đạo hàm Leibniz, bảng số liệu phân tích bậc từ thực nghiệm Chen & Chiew (2004), cấu trúc thuật toán Taylor-Galerkin và mã nguồn Fortran 90 (TG1D) đều được trình bày chi tiết từng bước.

  5. Chương trình nghiên cứu 10 năm tiếp theo được định hình ra sao? Mở rộng mô hình TG1D thành hệ thống phần mềm thủy lực 2D phi thủy tĩnh hoàn chỉnh, tích hợp module vận chuyển bùn cát xói lở đáy và tương tác thảm thực vật mềm ngập nước phục vụ kỹ thuật sinh thái sông ngòi.

Kết luận

  1. Luận án đã giải quyết thành công bài toán học thuật tồn tại lâu năm bằng việc thiết lập hệ phương trình vi phân dòng chảy hở một chiều suy rộng, tính đến ảnh hưởng của vận tốc thẳng đứng tương đối lớn tại đáy lòng dẫn ($w^*$).
  2. Mô hình lý thuyết đã loại bỏ giả thiết phân bố áp suất thủy tĩnh cổ điển của Saint-Venant (1871), bổ sung đầy đủ các thành phần động áp và tương tác ứng suất rối thông qua mô hình chiều dài xáo trộn Prandtl với hằng số Von Kármán $k = 0.54$.
  3. Thuật toán phần tử hữu hạn Taylor-Galerkin bậc 3 theo thời gian và không gian đã được xây dựng hoàn chỉnh, mã hóa thành công trên ngôn ngữ Fortran 90 thành bộ chương trình tính toán TG1D có tính ổn định và độ chính xác vượt trội.
  4. Thí nghiệm mô hình vật lý trong máng thủy lực đã xác nhận độ tin cậy tuyệt đối của thuật toán số, chứng minh hiện tượng dâng cao mực nước thượng lưu và tái phân bố vận tốc dòng chảy đáy dưới tác động của dòng nước trồi.
  5. Nghiên cứu mở ra các hướng phát triển quan trọng trong thủy lực học hiện đại, cung cấp công cụ tính toán đắc lực cho các bài toán kênh dẫn qua vùng địa chất ngầm phức tạp, đập ngầm và đê biển, đóng góp thiết thực cho sự nghiệp phát triển nguồn nước và an toàn công trình.