Tổng quan luận án

Hiện tượng nước nhảy là một trạng thái thủy lực đặc thù của dòng chảy hở khi chuyển tiếp đột ngột từ trạng thái chảy xiết sang trạng thái chảy êm, đi kèm với sự tiêu hao năng lượng lớn và xuất hiện khu vực xoáy cuộn có mạch động mạnh. Trong thực tế công trình, hiện tượng này được ứng dụng trong thiết kế tiêu năng sau đập tràn, sau cống điều tiết, thiết kế bể nén bùn của hệ thống xử lý nước thải và tính toán khẩu độ công trình cống, cầu.

Luận án "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" giải quyết khoảng trống nghiên cứu xuất phát từ thực tiễn: các nghiên cứu về nước nhảy trong lòng dẫn mặt cắt ngang hình thang cân còn hạn chế, các công thức hiện có chủ yếu là công thức thực nghiệm xây dựng trên các mô hình kích thước nhỏ (chiều rộng đáy $b = 0,2\text{ m}$), chưa được mô hình hóa đầy đủ bằng phương trình lý thuyết thủy động lực học. Đồng thời, trong thực tế thiết kế các công trình tiêu năng có mặt cắt hình thang (như Thủy điện Hồi Xuân, Hệ thống thủy lợi Nà Sản), kỹ sư vẫn phải sử dụng công thức của kênh hình chữ nhật rồi hiệu chỉnh qua thí nghiệm mô hình vật lý, dẫn đến sự sai lệch giữa tính toán ban đầu và thiết kế kỹ thuật.

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

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

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

  • Đối tượng nghiên cứu: Hiện tượng nước nhảy trong kênh lăng trụ, đáy bằng, không có mở rộng trên mặt bằng.
  • Phạm vi nghiên cứu: Nước nhả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, không có mở rộng trên mặt bằng trong điều kiện nước nhảy ổn định (số Froude $Fr_{D1} = 4,0 \div 9,0$).

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

Lịch sử nghiên cứu hiện tượng nước nhảy được ghi nhận từ Leonardo da Vinci (1452–1519), sau đó Giorgio Bidone (1818–1819) mô tả chi tiết bằng thực nghiệm và J. B. Bélanger (1828) xây dựng công thức giải tích đầu tiên xác định độ sâu liên hiệp trên kênh mặt cắt hình chữ nhật dựa trên định luật bảo toàn động lượng.

Về phân loại trạng thái nước nhảy theo số Froude trước nước nhảy ($Fr_1$), các nghiên cứu của Bradley và Peterka (1957), Hager (1992), Montes (1998) và Chanson (2004) đã hệ thống hóa thành 5 dạng:

Dạng nước nhảy Phạm vi số Froude ($Fr_1$) Đặc điểm động lực học và hiệu suất tiêu năng
Nước nhảy sóng $Fr_1 < 1,7$ Mặt nước dạng hình sóng, tổn thất năng lượng không đáng kể ($< 5%$).
Nước nhảy yếu / xoáy $1,7 \le Fr_1 < 2,5$ Xuất hiện các xoáy cuộn nhỏ trên mặt, tổn thất năng lượng từ $5% \div 15%$.
Nước nhảy dao động / chuyển tiếp $2,5 \le Fr_1 < 4,5$ Tia nước dao động từ đáy lên mặt không chu kỳ, tạo sóng lan truyền hạ lưu ($15% \div 45%$).
Nước nhảy ổn định $4,5 \le Fr_1 \le 9,0$ Tia nước thẳng hàng, khu xoáy bề mặt rõ rệt, cân bằng năng lượng tốt nhất ($45% \div 70%$).
Nước nhảy mạnh / không ổn định $Fr_1 > 9,0$ Vận tốc tia nước rất lớn, xáo trộn mạnh, hiệu quả tiêu năng đạt từ $70% \div 85%$.

Đối với độ sâu phân giới ($y_c$) trong kênh hình thang cân, phương trình cơ bản bậc 6 không có nghiệm tường minh giải tích, dẫn đến 4 nhóm phương pháp:

  1. Phương pháp thử dần và đồ thị: Dựa trên hàm số năng lượng đơn vị mặt cắt $e = f(y)$.
  2. Phương pháp số: Zhengzhong Wang (1998) dùng thuật toán Newton-Raphson; Vatankhah và Easa (2011), Vatankhah (2013); Tiejie Cheng, Jun Wang và Jueyi Sui (2018) dùng phương pháp chia đôi có xét góc nghiêng đáy; Farzin Salmasi (2020) dùng mạng thần kinh nhân tạo ANN.
  3. Phương pháp bán thực nghiệm: Rashidi (2015) kéo dài mái dốc quy đổi về mặt cắt tam giác cân tương đương.
  4. Phương pháp thực nghiệm: Agơrôtskin (1930), Straub (1982), Swamee P. (1993).

Đối với nước nhảy trên kênh hình chữ nhật, độ sâu liên hiệp và chiều dài nước nhảy đã được nghiên cứu sâu rộng qua các công trình của Rajaratnam (1965), Sarma và Newnham (1973, 1975), Hager (1992), Ead và Rajaratnam (2002), Simões (2008, 2012), Abolfazl N.G và cộng sự (2013), Sanjeev K.G và cộng sự (2013), Harry Edmar Schulz và cộng sự (2015), Martín Mundo-Molina và José Luis Díaz Pérez (2019).

Đối với nước nhảy trên kênh hình thang cân, các công trình quốc tế gồm có:

  • A. Rakhmanov (1930): Đề xuất quan hệ tỷ số độ sâu liên hiệp theo độ sâu phân giới $y_c$.
  • Wanoschek R. và Hager W. (1989): Xây dựng phương trình đại số xác định độ sâu liên hiệp cho kênh hình thang cân đáy bằng có $m = 1$, thiết lập đường cong quan hệ theo thông số $M = m y_1 / b$.
  • Sadiq S.M (2012): Phân tích cân bằng động lượng giải phương trình bậc 5 nhưng chưa có kiểm chứng thực nghiệm cho dải $M_1$ nhỏ.
  • Samir Kateb (2014, 2015): Nghiên cứu thực nghiệm trên mô hình đáy dốc góc mái $76,2^\circ$ ($m = 0,245$) và $73^\circ$ ($m = 0,31$), chiều rộng đáy hẹp $b = 20\text{ cm}$.
  • Sonia Cherhabil và Mahmoud Debabeche (2016): Thí nghiệm kênh hình thang $b = 20\text{ cm}$, $m = 0,31$ với 22 tấm ngưỡng tràn điều tiết.
  • SIAD Rafik (2018): Nghiên cứu đoạn nối 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}$).
  • Shahin S.A và cộng sự (2018): Thiết lập tương quan tuyến tính và phi tuyến cho kênh $m = 1,0$ và $m = 1,5$ với $b = 0,2\text{ m}$ nhưng hệ số tương quan $R^2$ chỉ đạt $0,802 \div 0,814$ do số lượng mẫu hạn chế (17 mẫu).
  • Bahador Fatehi Nobarian và cộng sự (2019): Mô phỏng 3D trường vận tốc và dòng chảy hai pha trên kênh hình thang có mái dốc $m = 0,26$; $0,58$ và $1,0$.
  • Về chiều dài nước nhảy trong kênh hình thang: Posey C. (1938), Silvester R. (1964), Rajaratnam và Subramanya (1968), Ohtsu (1976), N. Afzal (2002).

Tại Việt Nam, các nghiên cứu tiêu biểu bao gồm:

  • Hoàng Tư An (2005): Ứng dụng lý thuyết lớp biên của dòng tia chảy rối giải quyết bài toán nước nhảy trên lòng dẫn chữ nhật.
  • Nguyễn Văn Đăng (1989): Sử dụng mô hình dòng chảy Reynolds–Boussinesq đề xuất phương pháp hệ thức tích phân năng lượng xác định đặc trưng trung bình thời gian của nước nhảy phẳng.
  • Lưu Như Phú và Nguyễn Văn Toàn (2013): Nghiên cứu nước nhảy sau cửa van cống trên lòng dẫn chữ nhật.
  • Hồ Việt Hùng (2016), Hồ Việt Hùng và Lê Xuân Hiền (2017): Xác định kích thước bể tiêu năng có mố nhám lớn loại B và C theo tiêu chuẩn USBR trên kênh chữ nhật.
  • Lê Thị Việt Hà (2018): Nghiên cứu nước nhảy trong lòng dẫn mở rộng có độ dốc lớn.

Về ứng dụng Trí tuệ nhân tạo và Học máy, các nghiên cứu của Steven L. Brunton (2019), Melhem (1996), Mahdi Naseri và Faridah Othman (2012), Akram Abbaspour (2013), Masoud Karbasi và H. Azamathulla (2016), Ghorban Mahtabi (2020) chủ yếu tập trung vào mạng thần kinh nhân tạo (ANN), thuật toán di truyền (GP, GEP) hoặc kênh chữ nhật/lòng dẫn nhám. Các thuật toán Học máy dạng Cây quyết định (Decision Tree) và Rừng ngẫu nhiên (Random Forest) chưa được khai thác cho bài toán nước nhảy kênh hình thang cân.

Khoảng trống học thuật mà luận án lựa chọn giải quyết là: Thiết lập lời giải lý thuyết từ hệ phương trình Navier-Stokes kết hợp phương pháp lớp biên để xác định độ sâu trước và sau khu xoáy ($y_1, y_r$), kết hợp phương trình năng lượng để xây dựng công thức chiều dài khu xoáy ($L_r$) cho kênh hình thang cân đáy bằng có mái dốc $m = 1$ trong dải nước nhảy ổn định ($Fr_{D1} = 4,0 \div 9,0$), đồng thời kiểm chứng bằng mô hình vật lý kích thước lớn ($b = 33,5\text{ cm}$ và $b = 55\text{ cm}$) và công nghệ Học máy.

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

Cơ sở lý thuyết và khung phân tích

Khung phân tích của luận án dựa trên các nguyên lý thủy lực dòng chảy hở và thủy động lực học:

  • Trạng thái dòng chảy: Phân tích dòng chảy ổn định không đều, năng lượng đơn vị dòng chảy $E = Z + y + \alpha \frac{V^2}{2g}$ và năng lượng đơn vị mặt cắt $e = y + \frac{Q^2}{2gA^2}$.
  • Đặc tính dòng chảy ba chiều (bán không gian): Khác với nước nhảy phẳng trên kênh chữ nhật chỉ biến đổi theo mặt cắt dọc, dòng chảy trong kênh hình thang cân biến đổi theo 3 phương không gian ($Oxyz$) trên bề mặt mái dốc, hình thành các cấu trúc xoáy cuộn bề mặt dòng chính, xoáy cuộn thành bên, dòng tia đáy tách lớp và xáo trộn hai pha khí - lỏng.
  • Phương trình vi phân Navier-Stokes: Mô tả chuyển động của chất lỏng nhớt không nén được ($\rho = \text{const}$): $$\frac{\partial \mathbf{v}}{\partial t} + (\mathbf{v} \cdot \nabla)\mathbf{v} = \mathbf{F} - \frac{1}{\rho}\nabla p + \nu \nabla^2 \mathbf{v}$$
  • Lý thuyết lớp biên Prandtl (1904): Tách dòng chảy thành khu vực dòng chính và lớp biên sát biên rắn để lấy tích phân phương trình vi phân Navier-Stokes dọc theo thể tích kiểm tra giữa mặt cắt trước ($1-1$) và sau ($2-2$) khu xoáy.
  • Phương trình động lượng tích phân: Thiết lập mối quan hệ cân bằng giữa lực áp suất thủy tĩnh và động lượng dòng chảy có xét đến hệ số sửa chữa động lượng dòng rối $\beta_{0}^$: $$A_1 \left( \beta_{01}^ \frac{V_1^2}{g} + z_{c1} \right) = A_2 \left( \beta_{02}^* \frac{V_2^2}{g} + z_{c2} \right)$$ trong đó $z_{ci}$ là độ sâu trọng tâm mặt cắt ướt, $\beta_{0}^* = \beta_0 - \beta_\tau$ là hệ số sửa chữa động lượng có xét ảnh hưởng của khuếch tán dòng rối và lực ma sát.
  • Phương trình bảo toàn năng lượng: Dùng để phân tích quy luật tiêu hao năng lượng và biến đổi độ sâu dòng chảy dọc theo chiều dài khu xoáy $L_r$.
  • Phân tích thứ nguyên: Sử dụng định lý $\Pi$ của Buckingham để xác định các tổ hợp không thứ nguyên chi phối hiện tượng: $\frac{L_r}{y_1} = f(Fr_1, M_1, m, Re)$.

Phương pháp nghiên cứu thực tế

Luận án kết hợp đa phương pháp nghiên cứu:

  • Phương pháp phân tích và tổng hợp lý thuyết: Dẫn xuất giải tích hệ phương trình Navier-Stokes và phương trình vi phân năng lượng.
  • Phương pháp nghiên cứu thực nghiệm: Tiế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 vật lý được chế tạo bằng vật liệu lòng cứng, bao gồm đập tràn thực dụng xả lũ nối tiếp với bể tiêu năng hình thang cân đáy bằng có mái dốc $m = 1$, chiều rộng đáy mô hình lớn ($b = 33,5\text{ cm}$ và $b = 55\text{ cm}$) nhằm loại trừ ảnh hưởng của hiệu ứng quy mô nhỏ.
  • Phương pháp thống kê và xử lý số liệu: Sử dụng phần mềm MS Excel 2019 và IBM SPSS 26. Mức độ phù hợp giữa công thức tính toán và số liệu thực đo được đánh giá qua các chỉ số thống kê: Sai số phần trăm tuyệt đối trung bình (MAPE), Sai số bình phương trung bình (MSE), Căn bậc hai sai số bình phương trung bình (RMSE) và Hệ số xác định ($R^2$).
  • Phương pháp mô hình toán và Học máy: Xây dựng thuật toán Cây quyết định (Decision Tree) và Rừng ngẫu nhiên (Random Forest) để huấn luyện dữ liệu, dự báo độ sâu sau khu xoáy ($y_r$) và chiều dài khu xoáy ($L_r$).
  • Phương pháp kế thừa và chuyên gia: Kế thừa bộ dữ liệu thực nghiệm chuẩn của Wanoscheck R. (1989), đồng thời tham vấn các chuyên gia thủy lực qua các hội thảo khoa học chuyên ngành.

Nguồn tư liệu và dữ liệu

  • Dữ liệu thực nghiệm đo đạc trực tiếp trên mô hình vật lý tại Viện Khoa học thủy lợi Việt Nam (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$, phân bố lưu tốc mặt cắt ngang và dọc).
  • Bộ dữ liệu thực nghiệm kế thừa của Wanoscheck R. (1989) đối với kênh hình thang cân $m = 1$.
  • Tài liệu thiết kế kỹ thuật và số liệu thí nghiệm mô hình của 02 công trình thủy công thực tế: Công trình thủy điện Hồi Xuân (Thanh Hóa) và Công trình thủy lợi Nà Sản (Sơn La).

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 hệ thống hóa toàn diện các khái niệm thủy lực liên quan đến dòng chảy hở và hiện tượng nước nhảy. Tác giả phân tích bản chất cơ học của hiện tượng: nước nhảy xuất hiện khi dòng chảy chuyển từ trạng thái chảy xiết ($y < y_c$) sang trạng thái chảy êm ($y > y_c$), đặc trưng bởi sự mở rộng đột ngột của dòng chảy và sự hình thành khu xoáy cuộn có hướng dòng chảy ngược chiều với dòng chính.

Cấu trúc dòng chảy của nước nhảy trong kênh hình thang cân được tác giả mô tả chi tiết qua sơ đồ 7 vùng thủy động lực học:

  1. Vùng dòng chảy vận tốc lớn trước nước nhảy.
  2. Vùng dòng chảy chuyển động từ dưới lên tạo phân tán mạnh.
  3. Vùng các dòng tia đáy tách khỏi đáy do quán tính, tạo xoáy đáy.
  4. Vùng xoáy mặt chuyển động ngược chiều về phía thượng lưu, bao gồm xoáy bề mặt dòng chính và xoáy nêm dọc mái dốc.
  5. Vùng dòng chính tiếp tục chuyển động về hạ lưu.
  6. Vùng dòng đáy chuyển động men theo mái dốc tạo thành dòng cuộn hướng lên.
  7. Vùng hạ lưu không tham gia khu xoáy nhưng ảnh hưởng trực tiếp đến đuôi nước nhảy.

Tác giả xác lập rõ sự phân biệt giữa hai cặp thông số hình học:

  • Độ sâu liên hiệp hoàn toàn ($y_2$) và chiều dài nước nhảy tuyệt đối ($L_j$): xác định tại vị trí dòng chảy đã hoàn toàn ổn định ở hạ lưu.
  • Độ sâu sau khu xoáy ($y_r$) và chiều dài khu xoáy ($L_r$): xác định tại vị trí kết thúc vùng xoáy cuộn bề mặt. Tác giả chỉ ra rằng việc xác định điểm kết thúc của $L_r$ trên thực nghiệm có độ tin cậy và tường minh cao hơn $L_j$.

Chương I tổng hợp 10 công thức tính độ sâu phân giới ($y_c$), 06 nghiên cứu về độ sâu liên hiệp kênh chữ nhật, 30 công thức chiều dài nước nhảy ($L_j$) và 08 công thức chiều dài khu xoáy ($L_r$) kênh chữ nhật; 08 nghiên cứu độ sâu liên hiệp và 11 nghiên cứu chiều dài nước nhảy kênh hình thang. Phân tích thực tế tại Thủy điện Hồi Xuân (bể tiêu năng dài $30\text{ m}$ hình thang giúp tăng $5%$ hiệu quả tiêu năng so với thiết kế chữ nhật $39\text{ m}$) và Thủy lợi Nà Sản ($Q = 426,3\text{ m}^3/\text{s}$, $L_j = 39\text{ m}$, hiệu quả tiêu năng đạt $72,74% \div 76,99%$) chứng minh tính cấp thiết của việc chuẩn hóa công thức tính toán.

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 cơ sở toán học và cơ học chất lỏng phục vụ cho việc giải bài toán nước nhảy trên kênh lăng trụ hình thang cân đáy bằng ($i = 0$).

Quá trình tích phân hệ phương trình vi phân Navier-Stokes ba chiều được tác giả thực hiện dựa trên 6 giả thiết thủy lực:

  1. Dòng chảy ổn định, liên tục ($Q = \text{const}$).
  2. Áp suất phân bố theo quy luật thủy tĩnh.
  3. Lực khối duy nhất là trọng lực, thành phần chiếu lên phương dòng chảy $F_x = g \cdot i = 0$ (với đáy bằng).
  4. Chất lỏng không nén được ($\rho = \text{const}$).
  5. Mặt cắt kênh lăng trụ đối xứng.
  6. Thành phần ma sát rối tỷ lệ với cột nước lưu tốc dòng chảy.

Bằng cách lấy tích phân từng số hạng của phương trình Navier-Stokes trên diện tích mặt cắt ướt $A$, kết hợp phương trình liên tục và điều kiện biên lớp biên Prandtl, tác giả thu được phương trình cân bằng động lượng tổng quát cho hai mặt cắt trước và sau khu xoáy: $$A_1 \left( \beta_{01}^* \frac{V_1^2}{g} + z_{c1} \right) = A_2 \left( \beta_{02}^* \frac{V_2^2}{g} + z_{c2} \right)$$

Chương này cũng thiết lập phương trình vi phân năng lượng dọc theo chiều dài khu xoáy, làm cơ sở lý thuyết để tích phân tìm chiều dài $L_r$. Đồng thời, tác giả trình bày cấu trúc lý thuyết của mô hình Cây quyết định và Rừng ngẫu nhiên, thiết kế hệ thống thu thập số liệu trên mô hình máng thí nghiệm kích thước lớn tại Viện Khoa học thủy lợi Việt Nam.

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 trọng tâm chứa đựng các kết quả tính toán giải tích và thực nghiệm của luận án:

  • Xác định độ sâu sau khu xoáy ($y_r$): Từ phương trình động lượng tích phân ở Chương II, tác giả thiết lập tỷ lệ hệ số động lượng dòng chảy khu xoáy $k$ theo công thức (3.27) và xây dựng công thức lý thuyết (3.36). Kết quả kiểm nghiệm chỉ ra rằng công thức lý thuyết (3.36) đạt độ chính xác cao khi thông số không thứ nguyên $M_1 = \frac{m y_1}{b} \ge 0,2$ với giá trị tối ưu của tỷ lệ hệ số động lượng $k = 0,92$.
  • Mở rộng cho dải $M_1$ toàn phần: Để tính toán độ sâu $y_r$ cho mọi giá trị $M_1 \in [0, 1]$, tác giả xây dựng Bảng 3.5 và công thức thực nghiệm (3.39) có độ tin cậy thống kê cao.
  • Xác định chiều dài khu xoáy ($L_r$): Tác giả giải phương trình năng lượng để xác lập công thức quy luật biến đổi chiều dài (3.40). Kết hợp số liệu thực nghiệm từ mô hình vật lý ($b = 33,5\text{ cm}$ và $b = 55\text{ cm}$), tác giả 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) kèm theo điều kiện biên áp dụng (3.54) cho nước nhảy ổn định ($Fr_{D1} = 4,0 \div 9,0$).
  • Ứng dụng Học máy: Tác giả huấn luyện mô hình Rừng ngẫu nhiên trên tập dữ liệu thực nghiệm, chứng minh khả năng dự báo chính xác các thông số $y_r$ và $L_r$ tương thích với các công thức giải tích đề xuất.
  • Quy trình tính toán thực tế: Xây dựng sơ đồ quy trình từng bước ứng dụng các công thức mới vào tính toán kích thước bể tiêu năng cho công trình thủy công thực tế.

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

Đóng góp mới về mặt lý luận và khoa học

  1. Giải tích hệ phương trình Navier-Stokes: Lần đầu tiên hệ phương trình vi phân Navier-Stokes ba chiều được giải tích hợp với lý thuyết lớp biên Prandtl cho hiện tượng nước nhảy trong kênh lăng trụ hình thang cân đáy bằng ($i = 0$), thiết lập nên công thức lý thuyết tính độ sâu dòng chảy khu xoáy (công thức 3.36) và làm sáng tỏ bản chất vật lý của tỷ lệ hệ số động lượng khu xoáy $k$ (công thức 3.27).
  2. Cơ sở khoa học về chiều dài khu xoáy: Xây dựng công thức xác định quy luật biến đổi chiều dài dòng chảy khu xoáy $L_r$ (công thức 3.40) xuất phát từ phương trình vi phân cân bằng năng lượng.
  3. Ứng dụng công nghệ Học máy: Ứng dụng các thuật toán Học máy gồm Cây quyết định và Rừng ngẫu nhiên trong việc phân tích, dự báo các đặc trưng hình học của nước nhảy trên kênh hình thang cân.

Đóng góp mới về mặt thực tiễn và ứng dụng

  1. Xác lập phạm vi áp dụng và các công thức tính toán:
    • Xác định công thức lý thuyết (3.36) phù hợp áp dụng khi $M_1 \ge 0,2$ với giá trị hệ số động lượng $k = 0,92$.
    • Thiết lập Bảng 3.5 và công thức thực nghiệm (3.39) xác định độ sâu sau khu xoáy $y_r$ cho toàn bộ phạm vi $M_1 \in [0, 1]$ trong điều kiện nước nhảy ổn định ($Fr_1 = 4,0 \div 9,0$).
    • 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) kèm phạm vi áp dụng (3.54) để tính toán chính xác chiều dài khu xoáy $L_r$.
  2. Quy trình tính toán thiết kế công trình: Đề xuất quy trình tính toán kích thước hình học nước nhảy phục vụ thiết kế bể tiêu năng mặt cắt hình thang sau đập tràn, giảm thiểu sự phụ thuộc vào việc phỏng đoán từ công thức kênh chữ nhật và giảm khối lượng thí nghiệm mô hình hiệu chỉnh.

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

  • Hạn chế của nghiên cứu:
    • Nghiên cứu chỉ giới hạn trong điều kiện nước nhảy ổn định với số Froude trước nước nhảy nằm trong dải $Fr_{D1} = 4,0 \div 9,0$.
    • Mặt cắt kênh nghiên cứu là kênh lăng trụ hình thang cân có hệ số mái dốc cố định $m = 1$, đáy bằng phẳng ($i = 0$), không xét đến hiện tượng mở rộng hoặc thu hẹp trên mặt bằng.
    • Mô hình vật lý tập trung vào dòng chảy sau đập tràn mặt cắt thực dụng có bề mặt lòng dẫn cứng nhẵn, chưa xét tới các yếu tố nhám tăng cường hoặc mố tiêu năng.
  • 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 các dạng nước nhảy khác (nước nhảy sóng, nước nhảy yếu, nước nhảy dao động và nước nhảy mạnh $Fr_1 > 9,0$).
    • Khảo sát các kênh hình thang có hệ số mái dốc khác ($m \ne 1$) và kênh có độ dốc đáy ($i \ne 0$).
    • Nghiên cứu ảnh hưởng của độ nhám đáy, mố tiêu năng phân tán và đoạn nối tiếp mở rộng không gian đến cấu trúc dòng xoáy và kích thước hình học nước nhảy.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên môn cho các đối tượng và lĩnh vực sau:

  • Giảng viên, nghiên cứu sinh và học viên cao học ngành Kỹ thuật Thủy lợi, Kỹ thuật Tài nguyên nước, Thủy lực và Thủy động lực học: Tham khảo phương pháp giải tích phân hệ phương trình Navier-Stokes kết hợp lý thuyết lớp biên (Chương II) và phương pháp ứng dụng thuật toán Rừng ngẫu nhiên trong cơ học chất lỏng (Chương II, Chương III).
  • Kỹ sư tư vấn thiết kế công trình thủy lợi, thủy điện, giao thông và công trình xử lý nước thải: Vận dụng Bảng 3.5, công thức thực nghiệm (3.39), công thức bán thực nghiệm (3.50) và công thức (3.53) để tính toán trực tiếp kích thước bể tiêu năng đáy hình thang cân mà không cần quy đổi gần đúng qua kênh chữ nhật.

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

1. Đối tượng và phạm vi nghiên cứu cụ thể của luận án là gì?

Luận án nghiên cứu hiện tượng nước nhảy trong kênh lăng trụ, đáy bằng ($i = 0$), không có mở rộng trên mặt bằng. Phạm vi nghiên cứu thực nghiệm và lý thuyết tập trung vào nước nhả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$, trong điều kiện nước nhảy ổn định có số Froude trước nước nhảy $Fr_{D1} = 4,0 \div 9,0$.

2. Sự khác biệt cơ bản giữa cấu trúc thủy động của nước nhảy trong kênh hình thang cân và kênh hình chữ nhật là gì?

Nước nhảy trong kênh chữ nhật là nước nhảy phẳng, dòng xoáy cuộn chủ yếu diễn ra theo phương dọc trên mặt thoáng dòng chính. Ngược lại, nước nhảy trong kênh hình thang cân là hiện tượng "bán không gian": ngoài dòng xoáy cuộn bề mặt dòng chính, các phần tử chất lỏng chuyển động dọc theo mái dốc theo cả 3 phương ($Oxyz$), hình thành các khu vực xoáy nêm ở thành bên và dòng xoáy cuộn ở đáy đan xen phức tạp, làm tăng khả năng xáo trộn và tiêu hao năng lượng.

3. Luận án đã sử dụng hệ phương trình vi phân nào để xây dựng công thức xác định độ sâu dòng chảy khu xoáy?

Tác giả đã sử dụng hệ phương trình vi phân cơ bản Navier-Stokes ba chiều cho chất lỏng không nén được, kết hợp với các giả thiết thủy lực và phương pháp lớp biên Prandtl (1904) để lấy tích phân trên mặt cắt ướt $A$, từ đó thiết lập phương trình cân bằng động lượng có xét đến hệ số sửa chữa động lượng dòng rối $\beta_{0}^*$ để giải ra công thức lý thuyết tính độ sâu $y_1$ và $y_r$.

4. Những công cụ phần mềm và mô hình toán học nào đã được tác giả sử dụng?

Luận án sử dụng phần mềm MS Excel 2019 và IBM SPSS 26 để phân tích thống kê số liệu thực nghiệm (đánh giá qua các chỉ số MAPE, MSE, RMSE, $R^2$). Đồng thời, tác giả lần đầu tiên ứng dụng các thuật toán Học máy gồm Cây quyết định (Decision Tree) và Rừng ngẫu nhiên (Random Forest) để huấn luyện và dự báo các đặc trưng hình học của nước nhảy.

5. Kết quả cốt lõi về công thức tính độ sâu ($y_r$) và chiều dài ($L_r$) khu xoáy được luận án đề xuất như thế nào?

Luận án đã đề xuất công thức lý thuyết (3.36) phù hợp cho trường hợp $M_1 \ge 0,2$ với tỷ lệ hệ số động lượng tối ưu $k = 0,92$; thiết lập Bảng 3.5 và công thức thực nghiệm (3.39) cho dải $M_1 \in [0, 1]$; đồng thời xây dựng công thức bán thực nghiệm (3.50) và công thức thực nghiệm (3.53) theo điều kiện áp dụng (3.54) để tính toán chiều dài khu xoáy $L_r$ trong trạng thái nước nhảy ổn định.

Kết luận

Luận án đã giải quyết hoàn chỉnh 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 cân đáy bằng có mái dốc $m = 1$ trong điều kiện dòng chảy ổn định ($Fr_{D1} = 4,0 \div 9,0$). Bằng việc tích hợp giải tích phương trình Navier-Stokes với phương pháp thực nghiệm mô hình vật lý kích thước lớn và công nghệ Học máy Rừng ngẫu nhiên, tác giả đã thiết lập hệ thống công thức tính độ sâu ($y_r$) và chiều dài khu xoáy ($L_r$) có cơ sở khoa học và độ chính xác cao. Các kết quả nghiên cứu cung cấp cơ sở lý luận cho lĩnh vực thủy lực công trình và đóng vai trò là tài liệu chỉ dẫn kỹ thuật trực tiếp cho công tác tính toán, thiết kế bể tiêu năng hình thang trong thực tế xây dựng công trình thủy.