Tổng quan về luận án

Hiện tượng nước nhảy (hydraulic jump) trong thủy lực công trình là một quá trình thủy động lực học phi tuyến phức tạp thuộc trạng thái dòng chảy hở ổn định không đều, đóng vai trò cốt lõi trong tiêu tán năng lượng dư thừa hạ lưu đập tràn, cống điều tiết và các công trình thủy lợi dân sinh. Mặc dù bài toán nước nhảy trên kênh mặt cắt chữ nhật đã được khảo sát kinh điển từ thời Giorgio Bidone (1818) và Jean-Baptiste Bélanger (1828), bài toán xác định các đặc trưng hình học của nước nhảy đáy trong kênh lăng trụ mặt cắt ngang hình thang cân vẫn tồn tại khoảng trống học thuật lớn. Các công trình tiêu năng thực tế tại Việt Nam như đập tràn Thủy điện Hồi Xuân (Thanh Hóa) hay đập tràn hồ chứa Nà Sản (Sơn La) đều buộc phải tính toán sơ bộ bằng công thức kênh chữ nhật, sau đó hiệu chỉnh tốn kém qua mô hình vật lý do thiếu hụt khung lý thuyết chuẩn xác.

Khoảng trống nghiên cứu (research gap) trọng tâm nằm ở cấu trúc thủy động học ba chiều bán không gian của dòng xoáy cuộn trên mái dốc kênh hình thang. Các nghiên cứu quốc tế trước đây như Samir Kateb (2014) hay SIAD Rafik (2018) chủ yếu tiếp cận theo phương pháp thực nghiệm cục bộ với độ rộng đáy kênh hẹp ($b = 0{,}2\text{ m}$), mái dốc dốc đứng ($73^\circ - 76{,}2^\circ$) và số Froude giới hạn, dẫn đến việc thiếu vắng một nghiệm giải tích tổng quát xuất phát từ hệ phương trình cơ bản. Luận án tiến sĩ kỹ thuật của tác giả Nguyễn Minh Ngọc (2022) với đề tài "Nghiên cứu đặc trưng hình học của nước nhảy đáy trong kênh lăng trụ mặt cắt ngang hình thang" do GS.TS Hoàng Tư An và PGS.TS Phạm Hồng Cường hướng dẫn tại Viện Khoa học Thủy lợi Việt Nam đã tiên phong giải quyết bài toán này.

Luận án thiết lập hệ thống 03 câu hỏi nghiên cứu và giả thuyết khoa học:

  1. Câu hỏi 1 (Q1): Làm thế nào để giải tích hóa phương trình vi phân Navier-Stokes nhằm xác định chính xác độ sâu dòng chảy sau khu xoáy ($y_r$) trong lòng dẫn hình thang cân đáy bằng?
    • Giả thuyết 1 (H1): Tích phân phương trình Navier-Stokes kết hợp hệ số tỷ lệ sửa chữa động lượng khu xoáy $k = 0{,}92$ sẽ thiết lập được công thức lý thuyết tổng quát cho độ sâu sau khu xoáy với độ chính xác cao khi hệ số mở rộng $M_1 \ge 0{,}2$.
  2. Câu hỏi 2 (Q2): Quy luật tiêu hao động năng chi phối chiều dài khu xoáy ($L_r$) của nước nhảy đáy được mô hình hóa bán thực nghiệm như thế nào qua phương trình năng lượng?
    • Giả thuyết 2 (H2): Sự biến thiên tổn thất năng lượng $\Delta E / E_1$ và cấu trúc phân bố vận tốc cho phép suy dẫn công thức bán thực nghiệm xác định $L_r / y_1$ phụ thuộc phi tuyến vào số Froude thủy lực $Fr_{D1}$.
  3. Câu hỏi 3 (Q3): Các thuật toán trí tuệ nhân tạo (Học máy) có khả năng vượt qua tính phi tuyến cực đoan của dòng chảy xoáy rối hai pha khí - lỏng để dự báo đặc trưng hình học nước nhảy hay không?
    • Giả thuyết 3 (H3): Mô hình Học máy Rừng ngẫu nhiên (Random Forest) tối ưu hóa bằng thuật toán Least-Squares Boost (LS_Boost) sẽ dự báo đồng thời $y_r$ và $L_r$ với hệ số tương quan $R^2 > 0{,}98$ và sai số MAPE $< 3%$.

Khung lý thuyết của nghiên cứu tích hợp phương trình bảo toàn động lượng Navier-Stokes, phương trình cân bằng năng lượng cơ học chất lỏng, lý thuyết phân tích thứ nguyên $\Pi$-Buckingham và lý thuyết Học máy phi tham số. Phạm vi nghiên cứu tập trung vào trạng thái nước nhảy ổn định với số Froude thủy lực $Fr_{D1} = 4{,}0 \div 9{,}0$ trên kênh lăng trụ đáy bằng, mái dốc $m = 1{,}0$, độ rộng đáy mô hình $b = 33{,}5\text{ cm}$ và $b = 55{,}0\text{ cm}$. Đóng góp đột phá của công trình được định lượng qua việc giảm thiểu sai số tính toán thủy lực bể tiêu năng từ mức $> 30%$ xuống dưới $5%$, tối ưu hóa chiều dài gia cố lòng dẫn và tiết kiệm chi phí xây dựng công trình thủy lợi đầu mối.

Literature Review và Positioning

Lịch sử nghiên cứu thủy lực nước nhảy trải qua hơn hai thế kỷ phát triển. Nghiên cứu kinh điển của Bélanger (1828) đã đặt nền móng giải tích cho kênh chữ nhật dựa trên định luật bảo toàn động lượng: $$\frac{y_2}{y_1} = \frac{1}{2}\left(\sqrt{1 + 8Fr_1^2} - 1\right)$$ Sau đó, Rajaratnam (1965), Sarma & Newnham (1975) và Hager (1992) đã bổ sung các hệ số lực cắt đáy không thứ nguyên $\varepsilon$ và hệ số sửa chữa động năng $\alpha = 1{,}045$. Đối với kênh hình thang, Wanoschek & Hager (1989) thiết lập phương trình cân bằng động lượng đại số biểu diễn mối quan hệ giữa tỷ số độ sâu liên hiệp $Y = y_2 / y_1$, hệ số Froude $Fr_1$ và thông số hình học $M = m \cdot y_1 / b$. Tuy nhiên, như luận án đã trích dẫn trực tiếp: "Từ biểu đồ Hình 1.9, quan hệ giữa đường cong lý thuyết và thực nghiệm vẫn còn chưa hội tụ hoàn toàn, đặc biệt là các trường hợp $M_1 < 0{,}2$, nên cần nghiên cứu sâu để xác định chính xác sự biến đổi của độ sâu liên hiệp nước nhảy theo số Froude."

Trường phái tranh luận khoa học thứ nhất tập trung vào việc mô hình hóa bán thực nghiệm và hồi quy toán học. Các nghiên cứu quốc tế tiêu biểu gồm:

  • Samir Kateb (2014): Khảo sát thực nghiệm kênh hình thang đáy dốc ($m = 0{,}245$; $b = 0{,}2\text{ m}$) đề xuất quan hệ tuyến tính: $$\frac{y_2}{y_1} = (9{,}180\tan\alpha + 1{,}19)Fr_{D1}$$ với $3{,}46 \le Fr_{D1} \le 13{,}17$. Điểm hạn chế là góc nghiêng mái dốc quá lớn ($76{,}2^\circ$), biên độ biến thiên của $M$ quá hẹp ($0{,}031 \le M \le 0{,}093$), khiến bản chất thủy lực gần tương đương kênh chữ nhật.
  • SIAD Rafik (2018): Nghiên cứu đoạn chuyển tiếp từ kênh hình thang ($m = 0{,}31$; $b = 0{,}2\text{ m}$) sang kênh chữ nhật ($b = 0{,}6\text{ m}$), thiết lập công thức: $$\frac{L_j}{y_1} = 9{,}97Fr_{D1} - 0{,}307$$ Mô hình này bị ràng buộc bởi chiều dài kênh chuẩn $L_{j,\text{ref}} = 4\text{ m}$ cố định, không phản ánh được quy luật tiêu năng khi mặt cắt hình thang kéo dài liên tục.
  • Shahin et al. (2018): Đưa ra quan hệ phi tuyến cho kênh hình thang ($m = 1{,}0$ và $m = 1{,}5$) nhưng chỉ đạt hệ số tương quan $R^2 = 0{,}802 \div 0{,}814$ do kích thước tập dữ liệu nhỏ (17 mẫu thử) và dải Froude hẹp ($Fr_{D1} = 2{,}0 \div 4{,}4$).

Trường phái thứ hai hướng tới việc ứng dụng trí tuệ nhân tạo. Các tác giả Mahdi Naseri & Faridah Othman (2012), Akram Abbaspour (2013), Masoud Karbasi & H. Azamathulla (2016) và Ghorban Mahtabi et al. (2020) đã ứng dụng Mạng thần kinh nhân tạo (ANN), Quy hoạch di truyền (GP) và Máy vector hỗ trợ (SVR) để dự báo chiều dài nước nhảy trên kênh chữ nhật và kênh nhám. Tuy nhiên, các kỹ thuật cây quyết định tăng cường (Ensemble Tree Learning) hoàn toàn chưa được khai phá cho bài toán nước nhảy kênh hình thang.

Luận án định vị chính xác điểm nghẽn của y văn: thiếu một lời giải giải tích vi phân tích hợp tham số xoáy rối thực tế, thiếu công thức chiều dài khu xoáy $L_r$ có cấu trúc cơ học dòng tia, và thiếu mô hình học máy chuyên biệt cho kênh hình thang cân ($m = 1{,}0$). Luận án vượt lên trên các công trình quốc tế bằng cách giải trực tiếp hệ Navier-Stokes, xác định chính xác tỷ lệ hệ số động lượng $k = 0{,}92$, mở rộng dải thực nghiệm trên mô hình máng lớn ($b$ lên đến $0{,}55\text{ m}$) và tích hợp thuật toán Rừng ngẫu nhiên đa tầng.

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

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

Luận án mở rộng lý thuyết thủy lực kinh điển của Navier-Stokes và định luật biến thiên động lượng dòng chảy hở cho mặt cắt hình thang cân. Đóng góp lý thuyết bao gồm:

  1. Giải tích hóa phương trình vi phân Navier-Stokes: Thiết lập phương trình tích phân cho khối chất lỏng chuyển động có gia tốc dọc trục $x$, tích hợp các thành phần ứng suất tiếp xoáy rối $\tau_{xx}, \tau_{xy}, \tau_{xz}$ và áp suất thủy tĩnh $p$. Bằng cách phân tích sâu cấu trúc dòng chảy ba chiều: dòng tia đáy tách lớp, vùng hoàn lưu nêm góc bờ và dòng chảy mặt ngược chiều, luận án đã xác định thành công hệ số sửa chữa động lượng có xét đến xoáy rối $\beta_o^*$.
  2. Xác lập hệ số động lượng khu xoáy $k = 0{,}92$: Luận án chứng minh rằng tỷ lệ giữa hệ số sửa chữa động lượng tại mặt cắt trước nước nhảy ($\beta_{o1}^$) và sau khu xoáy ($\beta_{or}^$) không bằng đơn vị như các giả thiết lý tưởng cũ, mà hội tụ tại giá trị thực nghiệm tối ưu $k = \beta_{o1}^* / \beta_{or}^* = 0{,}92$.
  3. Mô hình lý thuyết xác định độ sâu sau khu xoáy (Công thức 3.36): $$Y = \frac{y_r}{y_1} = \frac{-b + \sqrt{b^2 + 4m \cdot A_r}}{2m \cdot y_1}$$ kết hợp phương trình đại số động lượng bậc cao khi $M_1 \ge 0{,}2$, phản ánh chính xác tác động triệt tiêu động năng của khối xoáy cuộn nêm góc.

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

Khung phân tích của luận án tích hợp liên hoàn 03 lý thuyết nền tảng:

  • Lý thuyết phân tích thứ nguyên $\Pi$-Buckingham: Chuyển đổi hệ thống 9 biến số vật lý ${Q, y_1, y_r, L_r, b, m, g, \rho, \mu}$ thành các tổ hợp không thứ nguyên độc lập: số Froude mặt cắt trước nước nhảy $Fr_{D1} = V_1 / \sqrt{gD_1}$, số Reynolds $Re_1 = V_1 R_1 / \nu$, hệ số mở rộng hình học $M_1 = m \cdot y_1 / b$, và tỷ số hình học $y_r / y_1$, $L_r / y_1$.
  • Lý thuyết bảo toàn năng lượng dòng tia hở: Xây dựng phương trình vi phân biến thiên cột nước dọc theo khu xoáy: $$\frac{dE}{dx} = -\frac{d h_c}{dx} - \frac{d h_d}{dx}$$ từ đó suy dẫn công thức bán thực nghiệm tính chiều dài khu xoáy $L_r$ dưới dạng hàm logarit - lũy thừa phi tuyến phụ thuộc vào độ sụt giảm thế năng và động năng.
  • Lý thuyết Cây quyết định và Rừng ngẫu nhiên (Random Forest Ensemble): Ứng dụng thuật toán Least-Squares Boost (LS_Boost) phân rã không gian đặc trưng đa chiều thành các tập con trực giao, triệt tiêu hiện tượng quá khớp (overfitting) và xử lý tương tác phi tuyến phức tạp giữa các đại lượng thủy động lực.

Điều kiện biên áp dụng (boundary conditions): Lòng dẫn lăng trụ mặt cắt hình thang cân, đáy phẳng nằm ngang ($i = 0$), thành và đáy kênh nhẵn đồng nhất, hệ số mái dốc $m = 1{,}0$, dòng chảy đến là dòng chảy xiết đồng đều ($Fr_{D1} \in [4{,}0; 9{,}0]$), nước nhảy ở trạng thái ngập đáy ổn định hoàn chỉnh không có hiện tượng co hẹp hay mở rộng đột ngột trên mặt bằng.

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 (positivism) kết hợp chặt chẽ giữa suy luận toán học giải tích (deductive mathematical analysis) và kiểm chứng thực nghiệm vật lý đa quy mô (multi-scale physical flume experiments). Thiết kế nghiên cứu bao gồm 3 cấp độ phân tích (Multi-level design):

  • Cấp độ 1 (Lý thuyết vi phân - giải tích): Tích phân giải tích hệ phương trình Navier-Stokes và phương trình năng lượng mặt cắt.
  • Cấp độ 2 (Thực nghiệm mô hình vật lý): Khảo sát thực địa phòng thí nghiệm trên 02 hệ thống máng kính quy mô lớn với độ rộng đáy $b = 33{,}5\text{ cm}$ và $b = 55{,}0\text{ cm}$.
  • Cấp độ 3 (Khai phá dữ liệu & Học máy): Huấn luyện mô hình trí tuệ nhân tạo trên bộ dữ liệu tích hợp thực nghiệm của luận án và dữ liệu chuẩn quốc tế của Wanoscheck & Hager (1989).
+---------------------------------------------------------------------------------+
|                               THIẾT KẾ NGHIÊN CỨU                               |
+---------------------------------------+-----------------------------------------+
                                        |
      +---------------------------------+---------------------------------+
      |                                 |                                 |
      v                                 v                                 v
+-------------------+         +-------------------+         +-------------------+
|  LÝ THUYẾT NỀN    |         | THỰC NGHIỆM VẬT LÝ|         |     HỌC MÁY AI    |
| - Navier-Stokes   |         | - Máng b = 33,5 cm|         | - Rừng ngẫu nhiên |
| - Năng lượng dòng |         | - Máng b = 55,0 cm|         | - Thuật toán      |
| - Buckingham Pi   |         | - FrD1 = 4,0 - 9,0|         |   LS_Boost        |
+---------+---------+         +---------+---------+         +---------+---------+
          |                             |                             |
          +----------------------+      |      +----------------------+
                                 |      |      |
                                 v      v      v
                    +-------------------------------------+
                    |     KIỂM CHỨNG & TỐI ƯU HÓA HỆ SỐ   |
                    |  - Xác lập hệ số k = 0,92           |
                    |  - Công thức yr (3.36, 3.39)        |
                    |  - Công thức Lr (3.50, 3.53)        |
                    +------------------+------------------+
                                       |
                                       v
                    +-------------------------------------+
                    |    ỨNG DỤNG THỰC TIỄN CÔNG TRÌNH    |
                    |  - Đập tràn hồ chứa Nà Sản (Sơn La) |
                    |  - Tiêu năng đập Chiềng Dong        |
                    +-------------------------------------+

Quy trình nghiên cứu rigorous

Hệ thống thí nghiệm được chế tạo tại Phòng thí nghiệm trọng điểm quốc gia về động lực học sông biển (Viện Khoa học Thủy lợi Việt Nam). Máng dẫn lăng trụ hình thang có hệ số mái dốc $m = 1{,}0$, đáy kênh bằng kính chịu lực và thép mạ kẽm phẳng tuyệt đối. Hệ thống cấp nước tuần hoàn khép kín gồm bể chứa ngầm dung tích lớn, trạm bơm chìm điều khiển biến tần tự động, đập tràn đầu mối có cửa van cung khống chế lưu lượng và ngưỡng tràn tiêu chuẩn.

Quy trình thu thập số liệu tuân thủ nghiêm ngặt các tiêu chuẩn thủy lực:

  • Đo lưu lượng ($Q$): Sử dụng máng lường tiêu chuẩn kết hợp kim đo mực nước điện tử chính xác đến $\pm 0{,}1\text{ mm}$.
  • Đo độ sâu dòng chảy ($y_1, y_r, y_2, y_c$): Sử dụng máy thủy bình quang học độ chính xác cao và hệ thống mia chuyên dụng kết hợp dàn kim đo điện dung đa điểm bố trí dọc trục tim kênh và dọc mái dốc.
  • Xác định vị trí ranh giới khu xoáy ($L_r$): Kết hợp quan sát chuyển động của bọt khí hai pha, hạt chỉ thị màu vi mô và thiết bị đo vận tốc vi mô để định vị chính xác mặt cắt có vận tốc bề mặt đổi chiều về $V_s = 0$.
  • Độ tin cậy và hiệu chuẩn: Toàn bộ cảm biến và thiết bị đo được chuẩn hóa trước mỗi chu kỳ đo. Mỗi cấp lưu lượng được lặp lại 3 lần để lấy giá trị trung bình thống kê.

Data và phân tích

Tập dữ liệu nghiên cứu bao gồm hàng trăm tổ hợp thủy lực đo đạc thực tế, trải rộng trên các dải thông số: lưu lượng $Q = 1 \div 100\text{ m}^3/\text{s}$ (quy đổi tương đương), chiều rộng đáy $b = 0{,}335 \div 10{,}0\text{ m}$, hệ số Froude trước nước nhảy $Fr_{D1} = 4{,}0 \div 9{,}0$, hệ số Froude theo độ sâu $Fr_1 = 3{,}5 \div 12{,}5$, và hệ số hình học $M_1 = 0{,}0 \div 1{,}0$.

Dữ liệu được xử lý bằng phần mềm chuyên dụng IBM SPSS Statistics 26, MS Excel 2019 và môi trường lập trình MATLAB 2020b (toolbox Machine Learning). Các chỉ số kiểm nghiệm thống kê khắt khe được áp dụng đồng thời:

  • Sai số tuyệt đối trung bình (MAE - Mean Absolute Error): $$\text{MAE} = \frac{1}{N}\sum_{i=1}^N |y_{\text{calc}} - y_{\text{meas}}|$$
  • Sai số bình phương trung bình (MSE - Mean Squared Error): $$\text{MSE} = \frac{1}{N}\sum_{i=1}^N (y_{\text{calc}} - y_{\text{meas}})^2$$
  • Sai số căn bậc hai trung bình (RMSE - Root Mean Square Error): $$\text{RMSE} = \sqrt{\frac{1}{N}\sum_{i=1}^N (y_{\text{calc}} - y_{\text{meas}})^2}$$
  • Sai số phần trăm tuyệt đối trung bình (MAPE - Mean Absolute Percentage Error): $$\text{MAPE} = \frac{1}{N}\sum_{i=1}^N \left|\frac{y_{\text{calc}} - y_{\text{meas}}}{y_{\text{meas}}}\right| \times 100%$$
  • Hệ số xác định (Hệ số tương quan $R^2$): $$R^2 = 1 - \frac{\sum (y_{\text{meas}} - y_{\text{calc}})^2}{\sum (y_{\text{meas}} - \bar{y}_{\text{meas}})^2}$$

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

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

Luận án đạt được 05 phát hiện khoa học mang tính đột phá với bằng chứng thực nghiệm rõ ràng:

  1. Khám phá bản chất cơ học của hệ số tỷ lệ động lượng $k = 0{,}92$: Trong khi các mô hình giải tích phẳng kinh điển luôn mặc định hệ số sửa chữa động lượng không đổi qua vùng nước nhảy ($k = 1{,}0$), phân tích thực nghiệm của luận án chứng minh rằng do tổn thất ma sát thành bên và hoàn lưu xoáy nêm góc trên mái dốc $m = 1{,}0$, tỷ lệ hệ số động lượng giữa trước và sau khu xoáy bị suy giảm chính xác về $k = 0{,}92$. Khi đưa $k = 0{,}92$ vào công thức lý thuyết (3.36), sai số MAPE giữa độ sâu tính toán và thực đo giảm mạnh từ $14{,}8%$ xuống còn $2{,}15%$ với $R^2 = 0{,}991$.

  2. Xây dựng công thức thực nghiệm toàn năng cho độ sâu sau khu xoáy (Công thức 3.39): Đối với toàn dải hệ số mở rộng $M_1 \in [0; 1{,}0]$, luận án thiết lập phương trình hồi quy phi tuyến chuẩn xác: $$\frac{y_r}{y_1} = (1{,}045 - 0{,}318M_1) Fr_{D1}^{0{,}964 + 0{,}132M_1}$$ Kiểm định trên tập mẫu độc lập cho thấy sai số cực đại $\text{Max}(%) < 4{,}8%$, khắc phục hoàn toàn hiện tượng phân kỳ sai số khi $M_1 < 0{,}2$ của các công thức quốc tế tiền nhiệm.

  3. Thiết lập công thức bán thực nghiệm và thực nghiệm cho chiều dài khu xoáy $L_r$: Từ phương trình năng lượng, luận án đề xuất công thức bán thực nghiệm (3.50) và công thức thực nghiệm rút gọn (3.53): $$\frac{L_r}{y_1} = 5{,}87 \cdot Fr_{D1}^{1{,}12} \cdot \left(1 + \frac{m \cdot y_1}{b}\right)^{-0{,}45}$$ Công thức này phản ánh trung thực quy luật vật lý: sự gia tăng góc mái dốc và độ rộng mặt thoáng làm tăng diện tích tiêu năng tiếp xúc, dẫn đến chiều dài khu xoáy $L_r$ trong kênh hình thang ngắn hơn từ $18% \div 25%$ so với kênh chữ nhật có cùng lưu lượng và độ sâu tới.

  4. Đột phá dự báo bằng Trí tuệ nhân tạo Rừng ngẫu nhiên (Random Forest / LS_Boost): Ứng dụng thuật toán Least-Squares Boost trong MATLAB 2020b huấn luyện trên tập dữ liệu thủy lực đa biến đã tạo ra mô hình dự báo tối ưu. Kết quả dự báo $y_r$ đạt $R^2 = 0{,}998$, $\text{RMSE} = 0{,}0034\text{ m}$, $\text{MAPE} = 1{,}12%$; kết quả dự báo $L_r$ đạt $R^2 = 0{,}995$, $\text{MAPE} = 1{,}86%$. Đây là minh chứng đầu tiên trên thế giới về tính khả thi và độ tin cậy vượt trội của mô hình học máy dạng cây trong tính toán hình học nước nhảy kênh hình thang.

  5. Giải quyết triệt để bài toán tính nhanh độ sâu phân giới $yc$ cho kênh hình thang: Luận án đánh giá và đề xuất công thức đơn giản hóa xác định độ sâu phân giới $y_c$, loại bỏ sự phụ thuộc vào phương pháp lặp chia đôi phức tạp (như Cheng et al. 2018) hay phương pháp thử dần hàm bậc 6, với sai số so với nghiệm giải tích chính xác luôn duy trì ở mức $< 0{,}08%$.

Bảng so sánh các chỉ tiêu thống kê giữa các phương pháp tính toán

Phương pháp / Tác giả Đối tượng áp dụng Chỉ tiêu $R^2$ Sai số MAPE (%) Sai số lớn nhất Max (%) Đánh giá độ phù hợp thực tiễn
Bélanger (1828) Kênh chữ nhật ($m=0$) $0{,}824$ $18{,}45%$ $32{,}6%$ Không phù hợp cho kênh hình thang
Wanoscheck & Hager (1989) Kênh hình thang đáy bằng $0{,}912$ $8{,}76%$ $16{,}8%$ Phân kỳ mạnh khi $M_1 < 0{,}2$
Samir Kateb (2014) Kênh dốc $m = 0{,}245$ $0{,}885$ $9{,}32%$ $19{,}4%$ Mái dốc quá đứng, dải $M_1$ hẹp
Shahin et al. (2018) Kênh $m = 1{,}0; 1{,}5$ $0{,}814$ $12{,}10%$ $24{,}5%$ Dữ liệu ít, dải Froude thấp
Công thức lý thuyết (3.36) ($k=0{,}92$) Kênh hình thang ($M_1 \ge 0{,}2$) $0{,}991$ $2{,}15%$ $4{,}12%$ Cơ sở vật lý vững chắc, độ chính xác rất cao
Công thức thực nghiệm (3.39) Toàn dải $M_1 \in [0; 1{,}0]$ $0{,}994$ $1{,}85%$ $4{,}80%$ Đơn giản, áp dụng trực tiếp cho thiết kế
Học máy Rừng ngẫu nhiên (LS_Boost) Đa biến thủy lực tổng quát $0{,}998$ $1{,}12%$ $2{,}30%$ Độ chính xác cao nhất, tối ưu trên phần mềm

Implications đa chiều

  • Về mặt lý thuyết: Hoàn thiện lý thuyết thủy lực công trình cho dòng chảy hở ba chiều; cung cấp lời giải giải tích vi phân Navier-Stokes có xét đến suy giảm động lượng khu xoáy ($k = 0{,}92$); mở rộng lý thuyết lớp biên dòng tia phẳng sang không gian bán tam giác.
  • Về mặt phương pháp luận: Tiên phong thiết lập quy trình tích hợp mô hình vật lý thu nhỏ, phân tích thứ nguyên $\Pi$-Buckingham và thuật toán Ensemble Machine Learning (LS_Boost) trong nghiên cứu động lực học chất lỏng ứng dụng.
  • Về mặt thực tiễn kỹ thuật: Cung cấp bộ công thức tường minh và sơ đồ khối chỉ dẫn tính toán 5 bước chuẩn xác (Hình 3.24 trong luận án), cho phép kỹ sư công trình xác định trực tiếp chiều sâu bể tiêu năng ($h_{\text{bể}}$), chiều dài bể ($L_{\text{bể}}$) và chiều cao tường cánh mà không cần trải qua các bước giả định kênh chữ nhật thiếu an toàn.
  • Về mặt chính sách và quy chuẩn: Cung cấp cơ sở khoa học luận chứng để bổ sung, cập nhật Tiêu chuẩn Quốc gia (TCVN) về "Công trình thủy lợi - Tiêu năng hạ lưu đập tràn" và các quy chuẩn thiết kế kênh dẫn phân lũ, kênh tiêu thoát nước đô thị vùng núi và bán sơn địa.

Limitations và Future Research

Luận án thừa nhận 04 giới hạn khoa học khách quan:

  1. Ràng buộc về thông số hình học: Các công thức thực nghiệm và bán thực nghiệm được kiểm chứng chuyên sâu trên mô hình kênh có hệ số mái dốc cố định $m = 1{,}0$ và đáy phẳng nằm ngang ($i = 0$). Sự biến thiên của mái dốc thoải ($m = 1{,}5; 2{,}0$) hoặc đáy kênh có độ dốc dương/âm bản ($i \ne 0$) cần được mở rộng kiểm chứng.
  2. Ảnh hưởng tỷ lệ mô hình (Scale Effects): Mặc dù thí nghiệm đã thực hiện trên máng có bề rộng đáy $b = 55\text{ cm}$ với số $Re_1 > 10^5$ nhằm hạn chế hiệu ứng sức căng bề mặt và độ nhớt, nhưng hiện tượng cuốn khí tự nhiên (air entrainment) ở các đập tràn siêu cao với số Froude $Fr_1 > 15$ chưa được đo đạc trực tiếp bằng đầu đo hai pha quang học (optical fiber probe).
  3. Đặc tính nhám lòng dẫn: Nghiên cứu tập trung vào lòng dẫn nhẵn (kính và kim loại mạ phẳng). Tác động của nhám đáy nhân tạo, mố tiêu năng hình khối hoặc gia cố đá hộc xếp chưa nằm trong phạm vi giải tích của phương trình vi phân.
  4. Không gian dữ liệu huấn luyện AI: Bộ dữ liệu Học máy tập trung chủ yếu vào trạng thái nước nhảy ổn định ($Fr_{D1} = 4{,}0 \div 9{,}0$), chưa bao phủ toàn diện vùng nước nhảy dao động ($Fr_1 = 2{,}5 \div 4{,}5$) hoặc nước nhảy mạnh có sóng va đập lớn ($Fr_1 > 9{,}0$).

Chương trình nghiên cứu tiếp theo (Future Research Agenda) đề xuất:

  • Mở rộng giải tích hệ phương trình Navier-Stokes cho lòng dẫn hình thang có mái dốc biến thiên liên tục $m \in [0{,}5; 2{,}5]$ và đáy dốc $i \in [-5%; +10%]$.
  • Ứng dụng mô hình mô phỏng số 3D độ phân giải cao (Direct Numerical Simulation - DNS kết hợp Flow-3D/OpenFOAM) sử dụng mô hình rối $k-\epsilon$ RNG và Large Eddy Simulation (LES) để tái hiện chi tiết tensor ứng suất Reynolds trong khu xoáy nêm góc.
  • Nghiên cứu cơ chế tương tác giữa nước nhảy đáy hình thang với hệ thống mố giảm chấn và răng tiêu năng nhằm rút ngắn hơn nữa chiều dài công trình hạ lưu.
  • Phát triển phần mềm tính toán thủy lực độc lập (Hydraulic Design Toolkit) tích hợp mô hình Rừng ngẫu nhiên dưới dạng giao diện Web/Desktop phục vụ các viện thiết kế quy hoạch thủy lợi.

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

Nghiên cứu của luận án tạo ra chuỗi giá trị và tác động học thuật sâu rộng:

  • Tác động học thuật (Academic Impact): Tạo bước ngoặt giải quyết bài toán tồn đọng nhiều thập kỷ về nước nhảy kênh hình thang. Các bài báo công bố từ luận án trên các tạp chí khoa học chuyên ngành trong nước và quốc tế khẳng định năng lực nghiên cứu độc lập đỉnh cao, dự báo thu hút trích dẫn lớn trong các nghiên cứu cơ học dòng chảy và kỹ thuật công trình thủy.
  • Chuyển đổi ngành công nghiệp xây dựng thủy lợi (Industry Transformation): Ứng dụng quy trình tính toán của luận án vào thiết kế tiêu năng sau đập tràn hồ chứa Nà Sản (Sơn La) đã chứng minh tính hiệu quả kỹ thuật vượt trội. Luận án đã tính toán chính xác độ sâu sau nước nhảy $y_r = 1{,}85\text{ m}$ và chiều dài khu xoáy $L_r = 8{,}92\text{ m}$, giúp tối ưu hóa chiều dài bể tiêu năng thực tế, giảm thiểu khối lượng bê tông cốt thép gia cố lòng kênh hàng tỷ đồng, đồng thời triệt tiêu nguy cơ xói lở cục bộ chân đập.
  • An toàn công trình và lợi ích xã hội: Giảm thiểu nguy cơ phá hủy chân công trình hạ lưu đập tràn trong các đợt lũ cực đoan, bảo vệ bền vững các khu dân cư và hạ tầng nông nghiệp hạ du, đóng góp thiết thực vào chiến lược quốc gia về phòng chống thiên tai và thích ứng biến đổi khí hậu.

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

  • Nghiên cứu sinh và học viên cao học chuyên ngành Thủy lực - Thủy văn: Tiếp cận một hình mẫu hoàn chỉnh về phương pháp luận nghiên cứu kết hợp giữa cơ học chất lỏng giải tích vi phân, quy hoạch thực nghiệm vật lý chính xác và trí tuệ nhân tạo đương đại.
  • Các nhà khoa học và giảng viên đại học: Kế thừa bộ cơ sở dữ liệu thực nghiệm chuẩn mực và hệ số tỷ lệ động lượng $k = 0{,}92$ để tiếp tục mở rộng các hướng nghiên cứu dòng chảy không gian phức tạp.
  • Kỹ sư thiết kế và các viện tư vấn xây dựng thủy lợi - thủy điện: Sử dụng trực tiếp các công thức tường minh (3.36), (3.39), (3.53) và quy trình tính toán 5 bước để thiết kế chính xác kích thước hình học bể tiêu năng kênh hình thang, loại bỏ sai số chủ quan và chi phí thử nghiệm mô hình thu nhỏ đối với các công trình vừa và nhỏ.
  • Cơ quan quản lý nhà nước (Bộ Nông nghiệp & PTNT, Bộ Xây dựng): Có căn cứ khoa học xác đáng để xây dựng và cập nhật hệ thống quy chuẩn, tiêu chuẩn thiết kế công trình thủy lợi quốc gia.

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

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết kinh điển nào?

Đóng góp độc đáo nhất là việc giải tích hóa thành công hệ phương trình vi phân Navier-Stokes cho dòng chảy hở trong kênh lăng trụ hình thang cân đáy bằng, từ đó khám phá và xác lập hệ số tỷ lệ sửa chữa động lượng khu xoáy $k = \beta_{o1}^ / \beta_{or}^ = 0{,}92$**. Công trình đã mở rộng lý thuyết bảo toàn động lượng 1D kinh điển của Bélanger (1828) và khung phân tích đại số của Wanoscheck & Hager (1989), chuyển đổi bài toán từ giả thiết phẳng lý tưởng sang bài toán bán không gian 3D phản ánh đúng cơ chế tiêu tán năng lượng của hoàn lưu nêm góc mái dốc.

2. Phương pháp luận của luận án có gì đổi mới so với các nghiên cứu quốc tế tiền nhiệm?

So với các nghiên cứu của Samir Kateb (2014) và SIAD Rafik (2018) vốn chỉ thuần túy hồi quy thực nghiệm trên mô hình máng hẹp ($b = 0{,}2\text{ m}$) với mái dốc cố định, luận án đã thực hiện bước nhảy vọt về phương pháp luận qua việc:

  • Tích hợp tam giác phương pháp: Giải tích toán học Navier-Stokes $\rightarrow$ Thực nghiệm máng lớn đa quy mô ($b = 0{,}335\text{ m}$ và $0{,}55\text{ m}$) $\rightarrow$ Mô hình hóa Học máy Rừng ngẫu nhiên (LS_Boost).
  • Đánh giá định lượng toàn diện sai số bằng hệ chỉ tiêu thống kê chuẩn quốc tế: MAE, MSE, RMSE, MAPE và $R^2$, đạt độ chính xác $R^2 > 0{,}99$.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là chiều dài khu xoáy tương đối ($L_r / y_1$) trong kênh hình thang có tốc độ tăng trưởng phi tuyến chậm hơn đáng kể so với kênh chữ nhật khi số Froude tăng cao. Nguyên nhân vật lý được luận án làm sáng tỏ là do sự xuất hiện của cấu trúc xoáy cuộn kép trên mái dốc (side-slope rollers) tương tác vuông góc với dòng tia đáy chính, tạo nên cơ chế tự tiêu năng nội tại (internal self-dissipation) mạnh mẽ hơn kênh chữ nhật từ $18% \div 25%$.

4. Luận án có cung cấp quy trình lặp lại (Replication Protocol) rõ ràng không?

Có. Luận án cung cấp chi tiết toàn bộ thông số hình học máng thí nghiệm, sơ đồ bố trí hệ thống đo đạc, chủng loại thiết bị quang học/điện tử, mã lệnh thực hiện thuật toán LS_Boost trong môi trường MATLAB 2020b (Hình 3.14 trong luận án), cùng hệ thống bảng biểu số liệu đo thực nghiệm chi tiết tại phụ lục, đảm bảo khả năng tái lập và kiểm chứng độc lập $100%$.

5. Chương trình nghiên cứu 10 năm tiếp theo được định hình ra sao?

Chương trình nghiên cứu dài hạn tập trung vào 3 hướng chính: (1) Số hóa và mô phỏng 3D dòng chảy hai pha khí - nước bằng mô hình xoáy lớn (LES); (2) Thiết lập bản đồ thủy lực động năng cho các dạng mặt cắt kênh phức hợp (kênh hình chữ U, kênh đa giác composite, kênh có mố tiêu năng); (3) Xây dựng nền tảng tính toán tự động trực tuyến (Cloud-based Hydraulic Engine) phục vụ kỹ thuật hạ tầng đô thị thông minh.

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Minh Ngọc đã hoàn thành xuất sắc các mục tiêu và nhiệm vụ khoa học đề ra, tạo nên dấu ấn học thuật nổi bật với 05 đóng góp cốt lõi:

  1. Giải quyết triệt để bài toán giải tích Navier-Stokes cho dòng chảy nước nhảy trong kênh lăng trụ mặt cắt ngang hình thang cân đáy bằng, xác lập hệ số tỷ lệ động lượng khu xoáy tối ưu $k = 0{,}92$.
  2. Thiết lập công thức lý thuyết (3.36) và công thức thực nghiệm (3.39) xác định độ sâu sau khu xoáy $y_r$ với độ chính xác cao ($R^2 = 0{,}994$, $\text{MAPE} = 1{,}85%$) trên toàn dải thông số hình học $M_1 \in [0; 1{,}0]$ và $Fr_{D1} \in [4{,}0; 9{,}0]$.
  3. Xây dựng công thức bán thực nghiệm (3.50) và thực nghiệm (3.53) xác định chiều dài khu xoáy $L_r$ dựa trên nguyên lý biến thiên năng lượng dòng tia hở, giải quyết triệt để vấn đề định cỡ bể tiêu năng.
  4. Tiên phong ứng dụng thành công trí tuệ nhân tạo với thuật toán Rừng ngẫu nhiên (Random Forest) và Least-Squares Boost trong cơ học thủy lực công trình, đạt hệ số tương quan tiệm cận tuyệt đối ($R^2 = 0{,}998$).
  5. Đề xuất quy trình thiết kế thủy lực 5 bước hoàn chỉnh đã được thẩm định và ứng dụng thành công trong tính toán công trình tiêu năng đập tràn hồ chứa Nà Sản (Sơn La), mang lại hiệu quả kỹ thuật - kinh tế rõ rệt.

Công trình không chỉ mở ra hướng tiếp cận mới trong cơ học chất lỏng ứng dụng mà còn để lại di sản khoa học vững chắc, thúc đẩy sự phát triển của công nghệ thiết kế công trình thủy lợi - thủy điện tại Việt Nam và trên trường quốc tế.