Tổng quan luận án

Hiện tượng nước nhảy đáy là quá trình chuyển tiếp trạng thái thủy lực từ dòng chảy xiết ($y < y_c$) sang dòng chảy êm ($y > y_c$) trong kênh hở. Hiện tượng này đóng vai trò cơ bản trong việc tiêu hao động năng dư thừa sau các công trình xả (chân đập tràn, sau cửa van cống điều tiết), hỗ trợ xáo trộn khí, thiết kế bể nén bùn và tính toán công trình thủy lợi, giao thông thủy. Mặc dù nước nhảy trên kênh chữ nhật đã được khảo sát sâu rộng từ công trình của Leonardo da Vinci, Bidone (1818–1819) và Bélanger (1828), các nghiên cứu đối với lòng dẫn lăng trụ mặt cắt ngang hình thang cân vẫn còn nhiều khoảng trống. Trong thực tế, các công trình xả qua kênh hình thang (như Thủy điện Hồi Xuân tại Thanh Hóa năm 2005, đập tràn Nà Sản tại Sơn La năm 2018) thường phải mượn công thức của kênh chữ nhật rồi dựa vào mô hình vật lý để hiệu chỉnh, dẫn đến sự sai lệch giữa tính toán thiết kế ban đầu và vận hành thực tế.

Luận án tiến sĩ kỹ thuật của tác giả Nguyễn Minh Ngọc, chuyên ngành Kỹ thuật xây dựng công trình thủy (mã số: 9-58-02-02), được thực hiện tại Viện Khoa học Thủy lợi Việt Nam dưới sự hướng dẫn khoa học của GS.TS. Hoàng Tư An và PGS.TS. Phạm Hồng Cường, tập trung giải quyết bài toán xác định đặc trưng hình học của nước nhảy đáy trong kênh lăng trụ hình thang.

Mục tiêu và nhiệm vụ nghiên cứu của luận án được xác định rõ:

  1. Mục tiêu nghiên cứu:
    • Xác định các đặc trưng hình học cốt lõi của nước nhảy, bao gồm độ sâu trước nước nhảy ($y_1$), độ sâu dòng chảy sau khu xoáy ($y_r$), độ sâu liên hiệp ($y_2$) và chiều dài dòng chảy khu xoáy ($L_r$) trong lòng dẫn lăng trụ mặt cắt ngang hình thang cân.
    • Đưa ra quy trình chỉ dẫn tính toán các kích thước hình học của nước nhảy hoàn chỉnh để phục vụ thiết kế công trình tiêu năng đáy bằng có mặt cắt ngang hình thang.
  2. Nhiệm vụ nghiên cứu:
    • Giải tích phân hệ phương trình vi phân cơ bản Navier-Stokes cho dòng chất lỏng chuyển động để xây dựng công thức lý thuyết xác định độ sâu dòng chảy khu xoáy của nước nhảy trong kênh hình thang cân.
    • Áp dụng phương trình bảo toàn năng lượng dòng chảy để thiết lập công thức bán thực nghiệm tính toán chiều dài dòng chảy khu xoáy ($L_r$).
    • Vận hành thí nghiệm mô hình vật lý nhằm thu thập cơ sở dữ liệu thực nghiệm, kiểm nghiệm và hiệu chỉnh các hệ số trong công thức lý thuyết, công thức bán thực nghiệm, đồng thời xây dựng công thức thực nghiệm và ứng dụng thuật toán học máy.

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

  • Đối tượng: Hiện tượng nước nhảy đáy trong lòng dẫn lăng trụ đáy bằng, không mở rộng trên mặt bằng.
  • Phạm vi nghiên cứu: Dòng chảy sau đập tràn mặt cắt thực dụng trong kênh lăng trụ hình thang cân có hệ số mái dốc $m = 1$, đáy bằng ($i = 0$), trong điều kiện nước nhảy ổn định với số Froude độ sâu thủy lực trước nước nhảy nằm trong khoảng $Fr_{D1} = 4,0 \div 9,0$.

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

Luận án đã hệ thống hóa các nghiên cứu tiêu biểu trong và ngoài nước về độ sâu phân giới, độ sâu liên hiệp và chiều dài nước nhảy:

  1. Xác định độ sâu phân giới ($y_c$) trong kênh hình thang:

    • Phương pháp giải tích số và giải thuật: Zhengzhong Wang (1998) áp dụng Newton-Raphson; Vatankhah & Easa (2011), Vatankhah (2013) phát triển công thức xấp xỉ liên tục; H. Arvanaghi và cộng sự (2015) khảo sát kênh hình thang; Tiejie Cheng, Jun Wang, Jueyi Sui (2018) dùng phương pháp dây cung chia đôi xét thêm góc nghiêng đáy kênh $\cos\theta$; Farzin Salmasi (2020) áp dụng mạng thần kinh nhân tạo (ANN).
    • Phương pháp bán thực nghiệm và thực nghiệm: Rashidi (2015) kéo dài mái dốc tạo tiết diện tam giác tương đương; Agơrôtskin (Liên Xô) thiết lập quan hệ qua độ sâu kênh chữ nhật $y_{cCN}$; Straub (1982) đề xuất công thức thực nghiệm với điều kiện $0,1 \le Q/b^{2,5} \le 4,0$; Swamee P. (1993) đưa ra công thức lũy thừa nhưng không suy biến được về dạng chữ nhật khi $m = 0$.
    • Đánh giá của tác giả: Đa số công thức hoặc đòi hỏi lặp nhiều bước, hoặc tồn tại sai số biên lớn (công thức Salmasi sai số cực đại tới 124%, Agorotskin sai số tới 91,28%, Straub sai số tới 33,87%).

  2. Xác định độ sâu liên hiệp và cấu trúc nước nhảy trong kênh chữ nhật:

    • Bélanger (1828) thiết lập nghiệm giải tích kinh điển dựa trên phương trình động lượng.
    • Rajaratnam (1965), Ead & Rajaratnam (2002) bổ sung ảnh hưởng của lực ma sát đáy $\varepsilon$ và dòng chảy rối.
    • Sarma & Newnham (1975) xét hệ số hiệu chỉnh động năng $\alpha = 1,045$; Hager W. (1992) và Harry Edmar Schulz và cộng sự (2015) đưa ra các hệ thức quan hệ tuyến tính hóa với số Froude $Fr_1$.

  3. Xác định đặc trưng nước nhảy trong kênh hình thang:

    • A. Rakhmanov (1930) đề xuất công thức tính $y_1, y_2$ qua $y_c$ nhưng chưa phản ánh tham số mở rộng $M = my/b$.
    • Wanoschek R. (1989, 1992) thiết lập phương trình bảo toàn động lượng cho mặt cắt hình thang cân nhưng đường cong lý thuyết chưa hội tụ tốt ở miền $M_1 < 0,2$.
    • Sadiq S.M (2012) giải phương trình bậc 5 dạng $aY^5 + bY^4 + cY^3 + dY^2 + eY + f = 0$ dưới dạng bảng tra lý thuyết, chưa kiểm chứng cho kênh rộng đáy nông.
    • Samir Kateb (2014), Samir Kateb và cộng sự (2015) nghiên cứu thực nghiệm trên kênh hình thang có góc nghiêng mái dốc $76,2^\circ$ ($m = 0,245$) và $73^\circ$ ($m = 0,31$) với đáy dốc.
    • Sonia Cherhabil & Mahmoud Debabeche (2016) thí nghiệm trên kênh $b = 20\text{ cm}, m = 0,31$ với các cấp độ sâu $y_1$ rời rạc.
    • SIAD Rafik (2018) nghiên cứu đoạn chuyển tiếp từ kênh hình thang ($b = 0,2\text{ m}, m = 0,31, L_{jref} = 4\text{ m}$) sang kênh chữ nhật ($b = 0,6\text{ m}$).
    • Shahin S.A và cộng sự (2018) thí nghiệm với $m = 1,0$ và $m = 1,5$ nhưng số lượng mẫu nhỏ ($N = 17$) và hệ số tương quan thấp ($R^2 \approx 0,802 \div 0,814$).
    • Bahador Fatehi Nobarian và cộng sự (2019) mô phỏng cấu trúc 3 chiều trường vận tốc trên Flow-3D nhưng chưa cung cấp công thức giải tích trực tiếp cho thiết kế.

  4. Nghiên cứu tại Việt Nam:

    • Các công trình của Hoàng Tư An (2005) về nước nhảy đáy, Nguyễn Văn Đăng (1998) về nối tiếp chảy đáy, Nguyễn Văn Cung về nước nhảy sóng, Lê Thị Việt Hà (2018) về kênh dốc mở rộng đã đặt nền móng tính toán thủy lực trong nước nhưng chưa giải quyết trọn vẹn nước nhảy trong kênh hình thang cân đáy bằng có xét trường xoáy 3 chiều.

Khoảng trống nghiên cứu được luận án lựa chọn là: Thiếu cơ sở lý thuyết tích phân trực tiếp từ phương trình Navier-Stokes có xét đến hệ số động lượng dòng xoáy rối ($k$), thiếu công thức giải tích và bán thực nghiệm tường minh cho độ sâu sau khu xoáy ($y_r$) và chiều dài khu xoáy ($L_r$) trên kênh hình thang cân đáy bằng $m = 1$, đồng thời cần ứng dụng giải thuật học máy hiện đại để kiểm chứng và dự báo các đặc trưng thủy lực này.

Cơ sở lý thuyết và phương pháp nghiên cứu

Luận án phối hợp đồng bộ giữa phương pháp giải tích lý thuyết, phương pháp thực nghiệm vật lý và mô hình toán học:

  1. Cơ sở lý thuyết cơ học chất lỏng:

    • Sử dụng hệ phương trình vi phân chuyển động Navier-Stokes cho dòng chất lỏng nhớt không nén được. Luận án tiến hành tích phân phương trình dọc theo thể tích kiểm soát của vùng nước nhảy trong kênh lăng trụ đáy bằng, đưa vào hệ số sửa chữa động lượng $\alpha_0$ và hệ số sửa chữa động lượng có xét đến xoáy rối $\alpha_0^$, từ đó xác định tỷ lệ động lượng khu xoáy $k = \alpha_0^ / \alpha_0$.
    • Thiết lập phương trình cân bằng năng lượng dòng chảy giữa mặt cắt trước nước nhảy (mặt cắt 1-1) và mặt cắt sau khu xoáy (mặt cắt r-r) để xác định quy luật tiêu tán năng lượng $\Delta E$ và suy dẫn công thức tính toán chiều dài khu xoáy $L_r$.

  2. Phân tích thứ nguyên và chuẩn tương tự:

    • Áp dụng định lý Buckingham $\Pi$ trên hệ cơ sở khối lượng - chiều dài - thời gian (MLT) để rút gọn các biến số đầu vào ($Q, b, y_1, V_1, m, \rho, \mu, g$) thành các nhóm không thứ nguyên chính: Số Froude trước nước nhảy $Fr_1$, số Froude theo độ sâu thủy lực $Fr_{D1} = V_1 / \sqrt{gD_1}$, số Reynolds $Re_1 = V_1 R_1 / \nu$, hệ số hình học mở rộng $M_1 = m y_1 / b$, tỷ số độ sâu $Y = y_r / y_1$ và tỷ số chiều dài $L_r / y_1$.
    • Mô hình vật lý tuân thủ chuẩn tương tự Froude cho dòng chảy hở có mặt thoáng tự do.

  3. Kỹ thuật mô hình hóa học máy (Machine Learning):

    • Xây dựng mô hình cây quyết định nhị phân (Decision Tree) và mô hình rừng ngẫu nhiên (Random Forest) kết hợp thuật toán Least-Squares Boosting (LS-Boost) của Friedman (2001) trên nền tảng phần mềm Matlab 2020b.
    • Mô hình sử dụng các biến đầu vào thủy lực và hình học để dự báo hai biến mục tiêu độc lập: Độ sâu dòng chảy sau khu xoáy ($y_r$) và chiều dài khu xoáy ($L_r$).

  4. Phương pháp thực nghiệm vật lý và xử lý dữ liệu:

    • Hệ thống thí nghiệm được lắp đặt và vận hành tại Phòng thí nghiệm trọng điểm quốc gia về động lực học sông biển thuộc Viện Khoa học Thủy lợi Việt Nam. Mô hình gồm cụm đập tràn thực dụng có cửa van điều tiết, kênh dẫn hình thang cân đáy bằng có mái dốc $m = 1$, với 2 quy mô chiều rộng đáy thực nghiệm: $b = 33,5\text{ cm}$ và $b = 55\text{ cm}$.
    • Thiết bị đo đạc gồm: Máng lường và kim đo lưu lượng đạt chuẩn, máy thủy bình kết hợp mia đo cao độ mực nước và đáy kênh, kim đo mực nước thượng lưu và hạ lưu đập, hệ thống van hạ lưu điều tiết vị trí chân nước nhảy.
    • Dữ liệu thu thập được phân tích thống kê trên phần mềm MS Excel 2019 và IBM SPSS 26. Mức độ chính xác và tương quan được kiểm định thông qua các chỉ tiêu thống kê:
      • Sai số tuyệt đối trung bình (MAE)
      • Sai số bình phương trung bình (MSE)
      • Sai số chuẩn (RMSE)
      • Sai số phần trăm tuyệt đối trung bình (MAPE, %)
      • Hệ số xác định tương quan ($R^2$)

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

CHƯƠNG I - TỔNG QUAN VỀ NƯỚC NHẢY ĐÁY

Chương I trình bày hệ thống khái niệm cơ bản về dòng chảy ổn định không đều, năng lượng đơn vị dòng chảy ($E = Z + y + \alpha V^2 / 2g$), năng lượng đơn vị mặt cắt ($e = y + \alpha Q^2 / (2gA^2)$), độ sâu phân giới ($y_c$) và số Froude ($Fr, Fr_D$).

Tác giả phân tích cấu trúc thủy động lực 3 chiều của nước nhảy đáy trong kênh hình thang: Do sự hiện diện của thành bên nghiêng, nước nhảy trong lòng dẫn hình thang mang bản chất "bán không gian". Cấu trúc dòng chảy được chia thành 7 vùng chức năng:

  • Vùng (1): Dòng chảy xiết vận tốc lớn trước nước nhảy.
  • Vùng (2): Dòng tia nâng từ đáy lên tạo phân tán mạnh.
  • Vùng (3): Hệ thống xoáy cuộn phân tách đáy do quán tính.
  • Vùng (4): Dòng xoáy bề mặt chuyển động ngược thượng lưu tại tim kênh và vùng hình nêm hai bên mái dốc.
  • Vùng (5): Trục dòng chính tiếp tục chuyển động về hạ lưu.
  • Vùng (6): Dòng cuộn dọc mái dốc chuyển động hỗn hợp từ dưới lên.
  • Vùng (7): Dòng hạ lưu sau khu xoáy ổn định dần.

Chương này tổng hợp bảng phân loại nước nhảy theo số Froude và hiệu suất tiêu năng:

Dải số Froude ($Fr_1$) Phân loại nước nhảy Trạng thái dòng chảy và mặt nước Hiệu suất tiêu năng (%)
$< 1,7$ Nước nhảy sóng Mặt nước xuất hiện gợn sóng lăn tăn $< 5$
$1,7 \div 2,5$ Nước nhảy yếu Xoáy cuộn nhỏ ở mặt nước, hạ lưu êm $5 \div 15$
$2,5 \div 4,5$ Nước nhảy dao động Tia nước dao động đáy - mặt, sinh sóng lớn $15 \div 45$
$4,5 \div 9,0$ Nước nhảy ổn định Tia nước đồng trục, cân bằng năng lượng tốt $45 \div 70$
$> 9,0$ Nước nhảy mạnh Xoáy cuộn bề mặt rất mạnh, sóng hạ lưu lớn $70 \div 85$

Phân tích tổng quan 30 công thức chiều dài nước nhảy kênh chữ nhật và hơn 10 công thức kênh hình thang khẳng định sự cần thiết phải phân biệt rõ chiều dài khu xoáy ($L_r$) và chiều dài nước nhảy tuyệt đối ($L_j$), trong đó $L_r$ là đại lượng có ranh giới cơ học rõ ràng, quyết định trực tiếp đến chiều dài gia cố đáy bể tiêu năng.

CHƯƠNG II – CƠ SỞ KHOA HỌC VÀ PHƯƠNG PHÁP XÁC ĐỊNH CÁC ĐẶC TRƯNG NƯỚC NHẢY TRONG KÊNH LĂNG TRỤ MẶT CẮT HÌNH THANG

Chương II tập trung thiết lập khung giải tích từ phương trình Navier-Stokes đối với kênh lăng trụ đáy bằng, áp dụng các giả thiết chất lỏng thực không nén được và tích phân qua thể tích giới hạn giữa mặt cắt 1-1 và r-r.

  1. Thiết lập phương trình chuyển động: Tác giả xác định phương trình bảo toàn động lượng có xét đến hệ số hiệu chỉnh xoáy rối $k$.
  2. Phân tích thứ nguyên: Dựa trên định lý Buckingham $\Pi$, xác định các biến độc lập ảnh hưởng đến tỷ số độ sâu sau khu xoáy $Y = y_r / y_1$ và tỷ số chiều dài khu xoáy $L_r / y_1$, gồm $Fr_{D1}, M_1, Re_1$.
  3. Khung thuật toán học máy: Thiết lập sơ đồ cây quyết định nhị phân và mô hình Rừng ngẫu nhiên sử dụng thuật toán Least-Squares Boosting (LS-Boost). Dữ liệu được phân chia thành tập huấn luyện (training set) và tập kiểm định (testing set) độc lập.
  4. Thiết kế và quy chuẩn mô hình thực nghiệm: Xây dựng hệ thống thí nghiệm với 2 quy mô chiều rộng đáy $b = 33,5\text{ cm}$ và $b = 55\text{ cm}$, hệ số mái dốc $m = 1,0$. Quá trình đo đạc ghi nhận đầy đủ các cấp lưu lượng $Q$, độ sâu trước nước nhảy $y_1$, độ sâu sau khu xoáy $y_r$, chiều dài khu xoáy $L_r$, cao độ đường mặt nước dọc lòng dẫn.
  5. Phân tích thực nghiệm sơ bộ: Khảo sát mối quan hệ giữa $L_r / y_1$ với $Fr_1$, $Fr_{D1}$, $Re_1$ và tỷ số tiêu hao năng lượng $\Delta E / E_1$, cho thấy $L_r / y_1$ có tương quan phi tuyến chặt chẽ nhất với số Froude độ sâu thủy lực $Fr_{D1}$ và tỷ số mở rộng $M_1$.

CHƯƠNG III – XÂY DỰNG CÔNG THỨC XÁC ĐỊNH CÁC ĐẶC TRƯNG HÌNH HỌC CỦA NƯỚC NHẢY TRONG KÊNH LĂNG TRỤ MẶT CẮT NGANG HÌNH THANG

Chương III là chương cốt lõi của luận án, tích hợp kết quả giải tích, số liệu thực nghiệm và kết quả mô hình học máy:

  1. Xây dựng công thức xác định độ sâu phân giới ($y_c$):

    • Đề xuất công thức tính $y_c$ cho kênh hình thang cân thông qua hàm số không thứ nguyên $g(M_c)$, cho kết quả sai số ổn định dưới $0,5%$ khi kiểm chứng trên dải lưu lượng $Q = 1 \div 100\text{ m}^3/\text{s}$, $b = 1 \div 10\text{ m}$, $m = 0,5 \div 2,0$.

  2. Xây dựng công thức độ sâu sau khu xoáy ($y_r$):

    • Giải phương trình tích phân Navier-Stokes, thiết lập công thức lý thuyết (3.36) xác định tỷ số độ sâu $Y = y_r / y_1$. Qua phân tích thống kê tương quan trên dữ liệu thực nghiệm, xác định giá trị tối ưu của tỷ lệ hệ số động lượng khu xoáy là $k = 0,92$.
    • Luận án chứng minh công thức lý thuyết (3.36) đạt độ chính xác cao khi $M_1 \ge 0,2$. Đối với toàn dải giá trị $M_1 \in [0, 1]$, luận án thiết lập công thức thực nghiệm (3.39) và lập Bảng 3.5 để tra cứu tỷ số $y_r / y_1$ theo $M_1$ và $Fr_{D1}$.

  3. Xây dựng công thức chiều dài khu xoáy ($L_r$):

    • Dựa trên phương trình cân bằng năng lượng và giải tích tổn thất cột nước cục bộ, thiết lập công thức biến đổi chiều dài dòng chảy (3.40).
    • Thiết lập công thức bán thực nghiệm (3.50) và công thức thực nghiệm (3.53) xác định chiều dài khu xoáy $L_r$ kèm phạm vi áp dụng (3.54).

  4. Kết quả mô hình học máy "Rừng ngẫu nhiên":

    • Sử dụng giải thuật LS-Boost trong Matlab 2020b, huấn luyện trên tập số liệu thực nghiệm. Mô hình Rừng ngẫu nhiên cho kết quả dự báo $y_r$ và $L_r$ có mức độ tương đồng cao với số liệu thực đo ($R^2 > 0,98$, các chỉ số sai số MAE, MSE, RMSE, MAPE đạt mức tối thiểu).

  5. Ứng dụng tính toán thực tế:

    • Đề xuất quy trình tính toán thủy lực hoàn chỉnh gồm các bước xác định thông số hình học cho công trình tiêu năng đáy bằng mặt cắt hình thang.
    • Áp dụng tính toán thực tế cho công trình tiêu năng sau đập tràn hồ chứa Nà Sản (Sơn La) và hồ chứa Chiềng Dong, kiểm chứng khả năng tiêu năng và xác định kích thước bể tiêu năng hình thang cân.

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

Luận án đã đạt được các kết quả khoa học và thực tiễn cụ thể:

  1. Đóng góp về mặt lý thuyết:

    • Đã giải tích phân hệ phương trình vi phân chuyển động Navier-Stokes cho dòng chất lỏng nhớt trong điều kiện nước nhảy ổn định ($Fr_{D1} = 4,0 \div 9,0$), thiết lập thành công công thức lý thuyết (3.36) xác định độ sâu dòng chảy khu xoáy của nước nhảy trong kênh lăng trụ hình thang cân đáy bằng ($m = 1$).
    • Xác định bản chất cơ học và lượng hóa giá trị tỷ lệ hệ số sửa chữa động lượng khu xoáy $k = 0,92$ (công thức 3.27), chứng minh công thức lý thuyết phù hợp nhất khi $M_1 \ge 0,2$.
    • Thiết lập công thức quy luật biến đổi chiều dài dòng chảy (3.40) từ phương trình cân bằng năng lượng dòng chảy.

  2. Đóng góp về mặt thực nghiệm và ứng dụng:

    • Xây dựng công thức thực nghiệm (3.39) và Bảng tra 3.5 cho phép xác định chính xác độ sâu dòng chảy sau khu xoáy ($y_r$) trên toàn dải thông số mở rộng hình học $M_1 \in [0, 1]$.
    • Thiết lập công thức bán thực nghiệm (3.50) và công thức thực nghiệm (3.53) cùng điều kiện biên (3.54) để tính toán chính xác chiều dài khu xoáy ($L_r$) của nước nhảy ổn định.
    • Ứng dụng giải thuật học máy Cây quyết định và Rừng ngẫu nhiên (thuật toán LS-Boost) trong xử lý dữ liệu và dự báo đặc trưng hình học nước nhảy trong kênh hình thang.
    • Đề xuất quy trình tính toán thủy lực từng bước phục vụ thiết kế công trình tiêu năng sau đập xả trong lòng dẫn mặt cắt hình thang, áp dụng thành công cho đập tràn hồ chứa Nà Sản (Sơn La).

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

Dựa trên nội dung công bố của luận án, các hạn chế và hướng nghiên cứu mở rộng được ghi nhận như sau:

  • Hạn chế: Các công thức lý thuyết và thực nghiệm mới chỉ khảo sát trên kênh lăng trụ đáy bằng phẳng ($i = 0$), không có độ mở rộng trên mặt bằng, mái dốc cố định $m = 1,0$, và phạm vi số Froude giới hạn trong trạng thái nước nhảy ổn định ($Fr_{D1} = 4,0 \div 9,0$).
  • Hướng nghiên cứu tiếp theo: Mở rộng nghiên cứu lý thuyết và thực nghiệm cho kênh hình thang có hệ số mái dốc biến thiên ($m \ne 1$), đáy kênh có độ dốc thay đổi ($i \ne 0$), dòng chảy có góc mở rộng trên mặt bằng, cũng như khảo sát các trạng thái nước nhảy dao động ($Fr_1 < 4,0$) hoặc nước nhảy mạnh ($Fr_1 > 9,0$).

Giá trị tham khảo

Tài liệu luận án cung cấp giá trị tham khảo chuyên môn cho các nhóm đối tượng:

  • Nghiên cứu sinh, học viên cao học và giảng viên ngành Kỹ thuật xây dựng công trình thủy, Thủy lực - Thủy văn: Tham khảo phương pháp tích phân phương trình Navier-Stokes cho dòng hở, phương pháp phân tích thứ nguyên Buckingham $\Pi$, quy hoạch thực nghiệm vật lý và kỹ thuật ứng dụng thuật toán học máy (LS-Boost Decision Tree / Random Forest) trong cơ học thủy lực.
  • Kỹ sư thiết kế và chuyên gia tư vấn công trình thủy lợi, thủy điện: Sử dụng quy trình tính toán kích thước hình học nước nhảy, công thức xác định độ sâu sau khu xoáy ($y_r$) và chiều dài khu xoáy ($L_r$) để thiết kế bể tiêu năng hình thang, giảm thiểu chi phí xây dựng mô hình vật lý hiệu chỉnh cho từng dự án.

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

  1. Tại sao nước nhảy trong kênh hình thang đáy bằng được xem là hiện tượng "bán không gian"? Nước nhảy trong kênh hình thang ngoài sự biến đổi vận tốc dọc theo chiều dòng chảy (trục x) giống như kênh chữ nhật, còn có sự phân bố và chuyển động xoáy 3 chiều của các phần tử chất lỏng trên bề mặt mái dốc nghiêng theo cả 3 phương trong không gian Oxyz, tạo ra cấu trúc xoáy cuộn phức tạp hơn nước nhảy phẳng.

  2. Hệ số tỷ lệ động lượng khu xoáy $k$ có ý nghĩa gì và giá trị tối ưu là bao nhiêu? Hệ số $k = \alpha_0^* / \alpha_0$ biểu thị tỷ lệ giữa hệ số sửa chữa động lượng có xét đến xoáy rối ($\alpha_0^*$) và hệ số sửa chữa động lượng thông thường ($\alpha_0$). Qua phân tích thực nghiệm trên mô hình kênh hình thang cân $m = 1$, tác giả xác định giá trị tối ưu là $k = 0,92$, giúp công thức lý thuyết (3.36) đạt độ chính xác cao nhất khi $M_1 \ge 0,2$.

  3. Phạm vi số Froude của trạng thái "nước nhảy ổn định" được luận án chọn nghiên cứu là bao nhiêu? Luận án lựa chọn phạm vi nghiên cứu với số Froude trước nước nhảy $Fr_{D1} = 4,0 \div 9,0$ (tương ứng với dải $Fr_1 = 4,5 \div 9,0$ theo phân loại của Peterka, Hager, Chanson), đây là vùng nước nhảy có tia nước thẳng hàng, cân bằng năng lượng tốt và đạt hiệu suất tiêu năng cao từ 45% đến 70%.

  4. Kích thước mô hình vật lý được tác giả chế tạo và thử nghiệm tại phòng thí nghiệm là bao nhiêu? Mô hình thí nghiệm được xây dựng 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) với kênh hình thang cân đáy bằng phẳng, hệ số mái dốc $m = 1$, gồm hai phương án chiều rộng đáy kênh là $b = 33,5\text{ cm}$ và $b = 55\text{ cm}$, kết nối sau đập tràn thực dụng có cửa van điều tiết.

  5. Công cụ học máy nào đã được tác giả ứng dụng để phân tích và dự báo đặc trưng nước nhảy? Tác giả đã lập trình giải thuật học máy Cây quyết định và Rừng ngẫu nhiên dựa trên thuật toán Least-Squares Boosting (LS-Boost) của Friedman (2001) trong môi trường Matlab 2020b để huấn luyện và dự báo độ sâu sau khu xoáy ($y_r$) và chiều dài khu xoáy ($L_r$).

Kết luận

Luận án tiến sĩ của Nguyễn Minh Ngọc đã hoàn thành toàn diện các mục tiêu nghiên cứu về đặc trưng hình học của hiện tượng nước nhảy đáy trong kênh lăng trụ mặt cắt ngang hình thang cân đáy bằng ($m = 1$). Bằng việc kết hợp giải tích hệ phương trình Navier-Stokes, phương trình năng lượng, thí nghiệm mô hình vật lý chính xác và kỹ thuật học máy Rừng ngẫu nhiên, tác giả đã thiết lập các công thức lý thuyết, bán thực nghiệm và thực nghiệm xác định $y_r$ và $L_r$. Kết quả nghiên cứu cung cấp cơ sở khoa học và quy trình tính toán thủy lực phục vụ thiết kế công trình tiêu năng thực tế, giải quyết khoảng trống kỹ thuật trong tính toán dòng nối tiếp sau công trình xả.