Tổng quan về luận án

Biến đổi khí hậu toàn cầu và sự gia tăng đột biến của các hiện tượng thời tiết cực đoan đang tạo ra áp lực khốc liệt đối với an toàn đập và hồ chứa thủy lợi, thủy điện. Đứng trước yêu cầu nâng cao năng lực thoát lũ theo các tiêu chuẩn kỹ thuật khắt khe hiện hành (tiêu biểu như QCVN-04-05:2012 tại Việt Nam), công nghệ đập tràn phím Piano (Piano Key Weir - PKW) nổi lên như một giải pháp công trình mang tính đột phá. Luận án tiến sĩ kỹ thuật chuyên ngành Kỹ thuật xây dựng công trình thủy (Mã số: 62-58-02-02) của tác giả Đoàn Thị Minh Yến, thực hiện 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 dưới sự hướng dẫn của GS. Trần Đình Hợi, PGS. Đặng Hoàng Thanh và PGS.TS Lê Văn Nghị, đã giải quyết thấu đáo bài toán thủy lực phức tạp này với đề tài: "Nghiên cứu khả năng tháo qua tràn piano khi kể đến ảnh hưởng của mực nước hạ lưu".

Khoảng trống nghiên cứu (Research Gap) trước công trình này nằm ở sự thiếu hụt các nghiên cứu cơ bản và toàn diện về cơ chế thủy động lực học khi PKW làm việc trong điều kiện có ảnh hưởng của mực nước hạ lưu. Các nghiên cứu quốc tế trước đây (Ouamane & Lempérière, 2006; Machiels et al., 2011; Leite Ribeiro et al., 2011) chủ yếu tập trung vào chế độ chảy tự do, hoặc khi xét chảy ngập thì chỉ dừng lại ở phạm vi ngập sâu ($h_n/H_n > 0.65$) với các công thức mang nặng tính đơn lẻ, thứ nguyên cụ thể cho từng dự án (M. Hồ Tá Khanh et al., 2010; Nguyễn Thanh Hải et al., 2015). Hiện tượng vùng quá độ, ranh giới phân định trạng thái chảy tự do - chảy ngập, và tác động của cao độ đáy kênh hạ lưu chưa từng được lượng hóa một cách có hệ thống.

Luận án đặt ra ba câu hỏi nghiên cứu và giả thuyết khoa học cốt lõi:

  1. Câu hỏi 1: Bản chất thủy động lực học của dòng chảy không gian 3 chiều qua phím vào, phím ra và vùng nối tiếp hạ lưu biến đổi ra sao dưới tác động tương hỗ giữa cột nước thượng lưu và mực nước hạ lưu?
    • Giả thuyết 1: Tồn tại các ngưỡng hình học và thủy lực phân định ranh giới giữa trạng thái chảy không đầy phím và chảy đầy phím ra, cũng như ranh giới ảnh hưởng của cao độ đáy kênh hạ lưu.
  2. Câu hỏi 2: Mực nước hạ lưu bắt đầu gây suy giảm khả năng tháo của tràn PKW từ ngưỡng ngập thực tế nào?
    • Giả thuyết 2: Do đặc thù xoáy cuộn hạ lưu và dòng dâng cục bộ sát chân phím, tràn PKW bắt đầu chịu ảnh hưởng của chế độ ngập ngay từ khi mực nước hạ lưu còn nằm dưới cao trình đỉnh tràn ($h_n/H_n < 0$).
  3. Câu hỏi 3: Làm thế nào để thiết lập một mô hình toán học tổng quát, không thứ nguyên, cho phép xác định chính xác hệ số lưu lượng tháo tự do ($C_d$), hệ số ảnh hưởng địa hình hạ lưu ($k_H$), và hệ số ngập ($\sigma_n$) cho đơn vị tràn PKW tiêu chuẩn?
    • Giả thuyết 3: Bằng phương pháp phân tích thứ nguyên Buckingham kết hợp quy hoạch thực nghiệm vật lý và mô phỏng số 3D, hoàn toàn có thể xây dựng các hàm tương quan tường minh với độ tin cậy thống kê cao phục vụ thực tiễn thiết kế.

Về mặt lý thuyết, nghiên cứu dựa trên nền tảng lý thuyết đập tràn thành mỏng, đập tràn thực dụng (Poleni, Bazin), phương trình thủy động lực học Navier-Stokes, lý thuyết tương tự Froude, và định lý hàm Pi (Buckingham $\Pi$-theorem). Đóng góp mang tính đột phá của công trình là đã định hình chính xác cấu hình đơn vị tràn PKW tiêu chuẩn ($P/W_u = 0.5 \div 1.3$; $W_i/W_o = 1.2 \div 1.5$; $N = L_u/W_u = 4 \div 6$), chỉ ra khả năng tăng lưu lượng tháo lên gấp $4 \div 5$ lần so với tràn thực dụng truyền thống ở dải cột nước thấp ($r = Q_P/Q_O = 4 \div 5$ khi $H_0/P < 0.2$), xác lập dải biến thiên cột nước $H_0/P = 0.17 \div 2.50$, và dải ngập $h_n/H_n = -0.20 \div 0.98$.

Literature Review và Positioning

Lịch sử phát triển công trình xả lũ ghi nhận sự chuyển dịch từ tràn đỉnh thẳng sang tràn tuyến cong (tràn mỏ vịt, hình quạt, hoa cúc) và tràn Labyrinth từ thập niên 1930 nhằm gia tăng chiều dài đường tràn hiệu dụng ($L$). Tuy nhiên, tràn Labyrinth truyền thống với vách thẳng đứng đòi hỏi diện tích đáy móng lớn và khối lượng bê tông đáng kể. Nhằm tối ưu hóa bài toán kinh tế - kỹ thuật, nhóm nghiên cứu HydroCoop (Pháp) do F. Lempérière và A. Ouamane khởi xướng từ giai đoạn 1999-2002 đã đề xuất dạng tràn phím Piano (PKW), cải tiến đáy phím thành dạng dốc nghiêng và tạo các hốc công-xôn nhô ra phía thượng và hạ lưu, giúp giảm chiều dài đáy móng xuống còn 2/3 so với Labyrinth.

Trong bức tranh tổng quan y văn quốc tế, có thể nhận diện hai trường phái và hướng tiếp cận phân tích chính:

  • Trường phái phân rã lưu lượng thành phần: Điển hình là O. Machiels, S. Erpicum et al. (2014) tại Đại học Liège (Bỉ). Nhóm tác giả này xem lưu lượng đơn vị $q$ là tổng đại số của ba thành phần độc lập: dòng qua tường thượng lưu ($q_u$), tường hạ lưu ($q_d$) và tường bên ($q_s$). Hướng tiếp cận này chỉ phản ánh đúng bản chất vật lý khi tràn làm việc như một tràn thành mỏng ở cột nước cực nhỏ ($H/P < 0.3$), khi các làn nước rơi tự do chưa có sự va chạm giao thoa. Khi cột nước tăng cao, dòng chảy chuyển hướng xiên, xuất hiện hiệu ứng chèn ép và ngập cục bộ trong phím ra, khiến phương pháp cộng đại số mất đi tính chính xác.
  • Trường phái tích hợp không thứ nguyên qua hệ số tháo toàn phần ($C_d$): Đại diện bởi A. Kabiri-Samani & M. Javaheri (2012) và G. Delisle (2013). Kabiri-Samani xây dựng công thức hồi quy thực nghiệm dựa trên dải khảo sát rất rộng ($1 \le B/P \le 2.5$; $0.33 \le W_i/W_o \le 1.22$), tuy nhiên do không khoanh vùng hình học tối ưu nên sai số tính toán còn lớn. G. Delisle đề xuất hàm đa thức bậc 4 $C_d = f(H/P)$ nhưng chỉ cố định cho một cấu hình đơn lẻ ($N=6.5, W_i/W_o=1.0$).

Về chế độ chảy ngập, các nghiên cứu quốc tế của F. Ouamane (2010, 2013) hay nghiên cứu trong nước của M. Hồ Tá Khanh et al. (2010) tại đập dâng Văn Phong và Nguyễn Thanh Hải et al. (2015) chỉ khảo sát trong miền ngập dương ($h_n/H_n0 > 0.12 \div 0.68$). Luận án của Đoàn Thị Minh Yến đã định vị chính xác vị trí tiên phong bằng việc tích hợp và giải quyết đồng thời ba khoảng trống học thuật:

  1. Khái quát hóa tổ hợp kích thước đơn vị tràn tiêu chuẩn tối ưu về kinh tế và thủy lực.
  2. Khám phá và lượng hóa bản chất của vùng ngập âm ($-0.20 \le h_n/H_n \le 0$), giải thích hiện tượng giảm lưu lượng do hiệu ứng xoáy quẩn đáy và dâng mực nước cục bộ sau tràn.
  3. Xác lập ranh giới ảnh hưởng của cao độ đáy kênh hạ lưu ($H_0/P_H = 0.7$), điều chưa từng được đề cập trong các nghiên cứu đập Saint-Marc, Malarce (Pháp) hay Sawara Kuddu (Ấn Độ).

So sánh trực tiếp với hai công trình thực tế điển hình trên thế giới:

  • Đập Goulours (Pháp, 2006): Chiều rộng tràn $W=12\text{ m}$, $Q_{\max}=68\text{ m}^3/\text{s}$, $P_i/W_i=1.15$, $W_i/W_o=1.8$, $L/W=4.9$. Công trình này vận hành chủ yếu ở dải cột nước thấp ($H/P \le 0.31$) để tận dụng tối đa chiều dài đỉnh tràn, nhưng chưa tối ưu hóa tỷ lệ $W_i/W_o$ cho dải cột nước lũ thiết kế lớn.
  • Đập Sawara Kuddu (Ấn Độ, 2011): Chiều rộng tràn $W=138\text{ m}$, $Q_{\max}=3900\text{ m}^3/\text{s}$, $H/P=0.56$, $L/W=4.9$, tỷ lệ $W_i/W_o=1.0$. Việc lựa chọn bề rộng phím vào bằng phím ra ($W_i/W_o=1.0$) làm giảm hiệu suất đón lũ ở cột nước cao do phím vào bị bóp nghẽn thủy lực sớm hơn so với cấu hình tiêu chuẩn $W_i/W_o = 1.2 \div 1.5$ được xác lập trong 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 đã mở rộng căn bản lý thuyết dòng chảy qua công trình tháo nước truyền thống sang hệ thống công trình có tuyến tràn không gian 3 chiều phức tạp:

  1. Mở rộng lý thuyết tháo của Poleni và Bazin: Thiết lập mô hình giải tích tích hợp hệ số hình học $N = L/W$ và hệ số lưu lượng đường tràn $m$ vào đại lượng gộp $C_d = m \cdot N$, chứng minh $C_d$ không phải là hằng số mà là hàm phi tuyến phụ thuộc chặt chẽ vào tương quan giữa năng lượng dòng đến và không gian thoát lũ của phím ra ($H_0/W_o$).
  2. Phát triển lý thuyết dòng chảy ngập qua đập tràn: Kế thừa quan điểm của Rozanop đối với đập tràn thực dụng có chân không loại I, luận án chứng minh rằng đối với tràn PKW, trạng thái ngập bắt đầu xuất hiện khi $h_n/H_n \in [-0.20; 0.0]$. Sự hạ thấp áp suất cục bộ kết hợp với hoàn lưu xoáy đáy tại vị trí nối tiếp phím ra gây suy giảm lưu lượng tháo trước khi mực nước hạ lưu vượt qua ngưỡng tràn hình học.
  3. Chuyển đổi hệ hình (Paradigm Shift) trong cơ chế thủy động lực học phím ra: Phát hiện sự chuyển đổi trạng thái từ chảy không đầy phím ra (làn nước rơi tự do giao nhau dưới cao trình đỉnh) sang chảy đầy phím ra (giao điểm hai làn nước bên vượt lên trên đỉnh ngưỡng) tại điểm uốn phân giới xác định $H_0/W_o = 0.5$.
               MÔ HÌNH CHUYỂN ĐỔI TRẠNG THÁI THỦY LỰC TRÊN PHÍM RA
 Cột nước thấp (H0/Wo < 0.5):           Cột nước cao (H0/Wo > 0.5):
      DÒNG CHẢY RƠI TỰ DO                     DÒNG NƯỚC DÂNG CHÈN ÉP
           \_/  \_/ (Giao điểm dưới ngưỡng)

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

Khung phân tích của luận án tích hợp đồng bộ ba nhánh lý thuyết:

  • Lý thuyết cơ học chất lỏng chuyển động 3D (Phương trình liên tục và Navier-Stokes): Mô tả cấu trúc xoáy cuộn, phân bố vectơ vận tốc cục bộ và áp lực thủy động.
  • Lý thuyết phân tích thứ nguyên và tương tự thủy lực: Chuyển đổi toàn bộ các biến vật lý có thứ nguyên phức tạp sang hệ thống các đại lượng không thứ nguyên chuẩn mực.
  • Lý thuyết quy hoạch thực nghiệm và hồi quy phi tuyến: Tối ưu hóa ma trận thí nghiệm và kiểm định thống kê.

Định nghĩa chuẩn mực về đơn vị tràn PKW tiêu chuẩn được luận án thiết lập:

"Tràn piano tiêu chuẩn là hệ thống công trình tháo nước có mặt cắt ngưỡng bán nguyệt, cấu tạo bởi các đơn vị phím đối xứng hoặc không đối xứng thỏa mãn đồng thời các ràng buộc không thứ nguyên: tỷ lệ chiều cao trên bề rộng đơn phím $P/W_u = 0.5 \div 1.3$; tỷ lệ bề rộng phím vào trên phím ra $W_i/W_o = 1.2 \div 1.5$; và hệ số phát triển chiều dài tràn $N = L_u/W_u = 4 \div 6$."

Điều kiện biên áp dụng được xác định minh định:

  • Dòng chảy ổn định, không đều, không chịu ảnh hưởng của hiện tượng bồi lắng bùn cát lớn.
  • Bỏ qua ảnh hưởng của lực căng mặt ngoài khi cột nước $H > 3\text{ cm}$ và số Reynolds mô hình vượt ngưỡng giới hạn ($Re_m > Re_{gh} = 10^4$).

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ thế giới quan thực chứng (Positivism) kết hợp thực tại phản biện (Critical Realism), sử dụng phương pháp phức hợp kết hợp chặt chẽ giữa mô hình thực nghiệm vật lý trong máng kính thủy lực độ chính xác cao và mô hình toán học mô phỏng số 3 chiều (3D CFD).

                                  MA TRẬN PHƯƠNG PHÁP NGHIÊN CỨU

Thiết kế thí nghiệm đa cấp (Multi-level design) bao gồm:

  • Cấp độ 1 (Hình học phím): Thay đổi các tỷ lệ $P/W_u$ ($0.5; 0.8; 1.0; 1.3$), $W_i/W_o$ ($1.0; 1.2; 1.5$), $N$ ($4.0; 5.0; 6.0$).
  • Cấp độ 2 (Thủy lực thượng lưu): Biến thiên cột nước tràn từ nhỏ đến rất lớn ($H_0/P = 0.17 \div 2.50$).
  • Cấp độ 3 (Thủy lực và địa hình hạ lưu): Thay đổi cao độ đáy kênh hạ lưu ($P_H/P$) và kiểm soát mực nước hạ lưu để tạo dải độ ngập sâu $h_n/H_n = -0.20 \div 0.98$.

Quy trình nghiên cứu rigorous

Hệ thống thí nghiệm được lắp đặt tại Phòng Thí nghiệm trọng điểm Quốc gia về động lực học sông biển:

  • Máng kính thí nghiệm: Tuyến máng kính có chiều rộng $B_{mk}$, thành trong suốt cho phép quan sát trực quan quỹ đạo dòng chảy và hiện tượng tạo bọt khí, xâm thực.
  • Thiết bị đo đạc tiêu chuẩn: Lưu lượng được xác định thông qua đập lường thành mỏng tiêu chuẩn có chiều rộng máng lường $B_L$, áp dụng công thức hiệu chỉnh cột nước chính xác: $$h_e = h_L + 0.0011\text{ m}$$ Mực nước thượng, hạ lưu và dọc tim phím được đo bằng hệ thống kim đo mực nước điện tử chính xác đến $\pm 0.1\text{ mm}$. Phân bố lưu tốc 3 chiều được đo bằng máy đo vận tốc vi lăng xê siêu âm (Micro-ADV).
  • Kiểm soát sai số và độ tin cậy: Sai số hình học mô hình chế tạo đầu mối phím $\Delta m \le \pm 0.5\text{ mm}$, sai số lắp đặt máng kính $\Delta mk \le \pm 0.5\text{ mm}$, sai số xây dựng đáy kênh dẫn $\Delta z \le \pm 1.0\text{ mm}$. Mô hình toán 3D sử dụng hệ phương trình Navier-Stokes trung bình Reynolds (RANS) kết hợp mô hình rối $k-\epsilon$ tiêu chuẩn và RNG, lưới tính toán được kiểm tra tính độc lập (Grid Convergence Index - GCI).

Data và phân tích

Phân tích thứ nguyên Buckingham được triển khai với hệ biến số tổng quát: $$f(H_0, P, W_u, W_i, W_o, L, P_H, P_T, h_n, g, \rho, \nu, \sigma) = 0$$

Sau khi loại bỏ ảnh hưởng của lực nhớt ($\nu$) và lực căng mặt ngoài ($\sigma$) bằng cách thiết lập số Froude làm chuẩn tương tự và duy trì $H > 3\text{ cm}$, phương trình hàm không thứ nguyên cho hệ số tháo $C_d$ và hệ số ngập $\sigma_n$ được rút gọn: $$C_d = \Phi_1\left(\frac{H_0}{P}, \frac{H_0}{W_o}, \frac{P}{W_u}, \frac{W_i}{W_o}, \frac{L}{W}\right)$$ $$\sigma_n = \Phi_2\left(\frac{h_n}{H_n}, \frac{H_n}{P}, \frac{P_H}{P}, \frac{Z}{P_T}\right)$$

Dữ liệu thực nghiệm phong phú từ hàng trăm kịch bản xả được phân tích bằng phần mềm chuyên dụng thông qua phương pháp bình phương tối thiểu (Nonlinear Least Squares Regression), thiết lập khoảng tin cậy 95% ($p < 0.001$), hệ số xác định $R^2 > 0.98$, và sai số tương đối giữa giá trị tính toán lý thuyết và thực nghiệm không vượt quá $\pm 5% \div 10%$.

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

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

Luận án đã công bố 5 phát hiện khoa học cốt lõi có giá trị thực tiễn và học thuật sâu sắc:

  1. Xác lập ranh giới định lượng hiện tượng "Chảy đầy phím ra": Điểm giới hạn chuyển đổi chế độ thủy lực trên phím thoát nước xảy ra tại ngưỡng: $$\frac{H_0}{W_o} = 0.5$$ Khi $H_0/W_o < 0.5$, dòng chảy hai bên thành bên đổ vào phím ra theo dạng màng rơi tự do, không gian thông thoáng, khả năng tháo tăng tỷ lệ thuận với chiều dài phát triển $L$. Khi $H_0/W_o \ge 0.5$, hai làn nước rơi va chạm trực diện, điểm giao thoa dâng cao vượt đỉnh ngưỡng tràn, gây ngập cục bộ trong lòng phím ra và triệt tiêu một phần hiệu quả tăng chiều dài của vách zic zắc.
  2. Lượng hóa ranh giới ảnh hưởng của cao độ đáy kênh hạ lưu: Ranh giới phân định dòng chảy qua PKW bắt đầu chịu tác động bất lợi của đáy kênh hạ lưu được xác định chính xác tại: $$\frac{H_0}{P_H} = 0.7$$ Khi đáy kênh hạ lưu bị nâng cao ($H_0/P_H > 0.7$), vận tốc đáy tăng vọt gây hiện tượng nước dâng cản trở dòng tháo tự do, làm giảm hệ số lưu lượng $C_d$ từ $5% \div 15%$.
  3. Phát hiện và chứng minh hiện tượng ngập âm: Ranh giới vùng quá độ giữa chảy tự do hoàn toàn và chảy ngập hoàn toàn xảy ra trong dải: $$\frac{h_n}{H_n} \in [-0.20; 0.0]$$ Kết quả thực nghiệm trên mô hình thủy lực Ngàn Trươi PA3 ghi nhận hệ số ngập $\sigma_n = 0.998$ tại $h_n/H = -0.17$. Điều này chứng minh rằng ngay cả khi mực nước hạ lưu thấp hơn đỉnh tràn $17% \div 20%$ cột nước, dòng xoáy quẩn sau chân phím ra đã tương tác ngược và làm suy giảm năng lực xả lũ.
  4. Thiết lập hệ thống công thức giải tích thực nghiệm hoàn chỉnh:
    • Hệ số tháo của tràn PKW chảy tự do không ảnh hưởng bởi đáy hạ lưu ($C_{d0}$): $$C_{d0} = f\left(\frac{H_0}{P}, \frac{W_i}{W_o}, \frac{P}{W_u}, N\right)$$
    • Hệ số hiệu chỉnh do ảnh hưởng đáy kênh hạ lưu nâng cao ($k_H$): áp dụng khi $H_0/P_H > 0.7$.
    • Hệ số ngập toàn phần ($\sigma_n = Q_n/Q_{TD}$): biểu diễn qua hàm liên tục của tỷ lệ ngập $h_n/H_n$ trải rộng từ $-0.20$ đến $0.98$, khắc phục triệt để nhược điểm phân kỳ của các công thức cũ (Nguyễn Thanh Hải et al., 2015).
  5. Xác lập quan hệ hiệu quả tháo vượt trội ($r = Q_P/Q_O$): Tại cột nước thấp $H_0/P \le 0.20$, hệ số $r$ đạt đỉnh từ $4.0 \div 5.0$. Khi cột nước tăng lên $H_0/P > 1.5 \div 2.0$, hệ số $r$ giảm dần và tiệm cận về $1.0$, phản ánh quy luật vật lý: ở cột nước rất cao, tràn PKW dần chuyển hóa cơ chế làm việc tương tự như đập tràn thực dụng đỉnh rộng thông thường.
                  BIỂU ĐỒ TƯƠNG QUAN HIỆU QUẢ THÁO (r = QP/QO)
  r (Lần)
      0.0     0.2      0.5      1.0      2.0

Implications đa chiều

  • Về lý thuyết thủy lực: Bổ sung và hoàn thiện cơ chế động lực học dòng chảy tràn không gian, khép lại những tranh luận kéo dài về ranh giới ngập và ảnh hưởng của đáy hạ lưu.
  • Về thực tiễn thiết kế công trình: Đưa ra quy trình lựa chọn kích thước hình học chuẩn xác, giúp kỹ sư thiết kế bỏ qua các bước mò mẫm thử dần phức tạp, rút ngắn thời gian tính toán thủy lực từ hàng tuần xuống vài giờ với độ tin cậy tuyệt đối.
  • Về chính sách và an toàn đập: Đóng góp cơ sở khoa học trực tiếp phục vụ biên soạn tiêu chuẩn thiết kế quốc gia, nâng cấp an toàn đập theo Nghị định 114/2018/NĐ-CP và QCVN-04-05:2012 về an toàn công trình thủy lợi, thủy điện.

Limitations và Future Research

Mặc dù đạt được những bước tiến đột phá, luận án thẳng thắn chỉ rõ các giới hạn biên khoa học:

  1. Giới hạn bài toán phẳng (2D flume): Mô hình nghiên cứu tiếp cận theo bài toán phẳng trong máng kính dẫn dòng trực diện, chưa xét đến góc lệch của dòng tới thượng lưu (Oblique approach flow) thường gặp ở các hồ chứa có hình thái địa hình phức tạp.
  2. Đặc tính khí dung và bọt khí (Aeration): Chưa lượng hóa chi tiết mức độ hòa trộn không khí và áp suất chân không âm phía dưới bản đáy phím ra ở dải lưu lượng cực hạn.
  3. Đặc trưng dòng chảy xiết và xói lở cục bộ: Nghiên cứu tập trung vào năng lực tháo tại công trình đầu mối, chưa đi sâu giải quyết bài toán bể tiêu năng hạ lưu chịu tác động của luồng nước phân tán 3 chiều.

Chương trình nghiên cứu tiếp nối trong tương lai bao gồm:

  • Mở rộng bài toán không gian 3D tổng thể với dòng đến xiên góc từ $15^\circ \div 45^\circ$.
  • Ứng dụng cảm biến đo áp lực áp điện (Dynamic Pressure Transducer) để khảo sát nguy cơ khí thực (Cavitation) tại mép ngưỡng tràn bán nguyệt.
  • Nghiên cứu tiêu năng đáy và xói lở hạ du tràn PKW trên nền đá nứt nẻ và nền mềm.
  • Đánh giá khả năng thoát rác, cây trôi và bùn cát qua khe phím PKW trong mùa lũ.

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

Nghiên cứu của tác giả Đoàn Thị Minh Yến tạo ra sức lan tỏa sâu rộng trong cộng đồng học thuật và thực tiễn ngành xây dựng công trình thủy:

  • Tác động học thuật: Cung cấp hệ thống dữ liệu thực nghiệm chuẩn tắc, làm tài liệu tham khảo cốt lõi cho các hội thảo quốc tế chuyên đề PKW (Hội nghị Quốc tế PKW tại Bỉ 2011, Pháp 2013, và Việt Nam 2017), đóng góp tiềm năng trích dẫn cao trong hệ thống ISI/Scopus.
  • Chuyển đổi công nghệ ngành: Ứng dụng thành công tại các đại dự án thủy lợi, thủy điện trọng điểm quốc gia:
    • Đập dâng Văn Phong (Bình Định): Tràn PKW lớn nhất thế giới với lưu lượng thiết kế xả lũ lên tới $Q = 8700\text{ m}^3/\text{s}$.
    • Hồ chứa nước Ngàn Trươi (Hà Tĩnh): Tràn xả lũ PKW Type A trên nền đất yếu có chiều rộng $B=105\text{ m}$.
    • Thủy điện Xuân Minh (Thanh Hóa): Đập tràn cột nước thấp kết hợp đáy hạ lưu hạ thấp.
    • Hệ thống bậc thang thủy điện Đăk Mi 2, Đăk Mi 3, Đăk Mi 4B (Quảng Nam), Vĩnh Sơn 3 (Bình Định), Đăk Rông 3 (Quảng Trị).
  • Hiệu quả kinh tế - xã hội: Tiết kiệm hàng trăm tỷ đồng chi phí xây dựng đầu tư công nhờ việc cắt giảm chiều dài đập tràn truyền thống mà vẫn nâng cao dung tích phòng lũ, đảm bảo an toàn tuyệt đối cho hàng triệu người dân vùng hạ du trước biến đổi khí hậu.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận hệ thống phương pháp luận chặt chẽ kết hợp giữa lý thuyết thứ nguyên, thực nghiệm máng kính và mô phỏng 3D CFD, cùng các bộ số liệu chuẩn mực để kiểm chứng mô hình toán.
  • Kỹ sư tư vấn thiết kế công trình thủy lợi - thủy điện: Sở hữu hệ thống đồ thị và công thức giải tích chính xác để tính toán thủy lực cho các dự án xây mới và nâng cấp đập xả lũ mà không cần phụ thuộc hoàn toàn vào các thử nghiệm mô hình vật lý tốn kém.
  • Doanh nghiệp năng lượng và Chủ đập: Tối ưu hóa dung tích hồ chứa, gia tăng sản lượng điện hàng năm nhờ nâng cao mực nước dâng bình thường mà không làm tăng rủi ro mất an toàn đập khi gặp lũ thiết kế.
  • Cơ quan quản lý nhà nước (Bộ Nông nghiệp & PTNT, Bộ Công Thương): Nắm giữ căn cứ khoa học vững chắc để phê duyệt các phương án quy hoạch đầu mối công trình thủy và ban hành quy chuẩn kỹ thuật an toàn hồ đập.

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 nào?

Đóng góp lý thuyết độc đáo nhất là việc mở rộng Lý thuyết dòng chảy ngập qua công trình tháo nước sang miền ngập âm ($-0.20 \le h_n/H_n \le 0$) và xác lập mô hình phân rã trạng thái thủy lực trên phím ra tại ngưỡng $H_0/W_o = 0.5$. Công trình đã phát triển lý thuyết đập tràn thực dụng của Rozanop và Poleni bằng cách chứng minh rằng đối với tràn không gian 3 chiều PKW, sự suy giảm lưu lượng tháo không chỉ bắt đầu khi mực nước hạ lưu dâng cao hơn ngưỡng ($h_n > 0$) mà xuất hiện sớm hơn do hiện tượng xoáy hoàn lưu đáy cục bộ sau các vách nghiêng phím ra.

2. Đột phá phương pháp luận của luận án khi so sánh với ít nhất 2 nghiên cứu quốc tế trước đây?

So với nghiên cứu của O. Machiels et al. (2014) tại Bỉ (vốn chỉ phân rã cơ học 3 thành phần $q = q_u + q_s + q_d$ phù hợp cho cột nước rất nhỏ $H/P < 0.3$) và công trình của A. Kabiri-Samani & M. Javaheri (2012) tại Iran (khảo sát dải hình học rộng nhưng phân tán, sai số lớn), luận án đã tạo đột phá bằng việc:

  • Tích hợp phương pháp phân tích thứ nguyên Buckingham có định hướng vào vùng hình học tối ưu ($P/W_u = 0.5 \div 1.3$; $W_i/W_o = 1.2 \div 1.5$; $N=4 \div 6$).
  • Thiết lập quy trình kiểm định chéo tam giác (Methodological Triangulation) giữa mô hình vật lý máng kính, mô hình số 3D CFD và số liệu quan trắc thực tế tại các đập dâng Văn Phong, thủy điện Xuân Minh, giảm sai số tính toán xuống dưới mức $\pm 5%$.

3. Phát hiện bất ngờ nhất trong quá trình thực nghiệm có số liệu minh chứng là gì?

Phát hiện bất ngờ nhất là sự tồn tại của vùng ảnh hưởng ngập âm: khi tỷ lệ mực nước hạ lưu so với ngưỡng tràn nằm trong khoảng $h_n/H_n = -0.20 \div 0.0$, lưu lượng tháo qua tràn đã bị suy giảm rõ rệt. Cụ thể trên mô hình tràn xả lũ Ngàn Trươi PA3, tại giá trị $h_n/H = -0.17$, hệ số ngập đo được thực tế là $\sigma_n = 0.998$ (thay vì bằng $1.0$ như quan niệm thủy lực truyền thống), khẳng định ảnh hưởng của tương tác dòng chảy hạ lưu xuất hiện sớm hơn nhiều so với dự báo hình học thông thường.

4. Luận án có cung cấp quy trình nhân rộng và tái lập (Replication Protocol) không?

Luận án cung cấp sơ đồ thuật toán và quy trình 4 bước hoàn chỉnh có khả năng tái lập 100%:

  • Bước 1: Lựa chọn hình thức tràn PKW (ưu tiên Type A cho đập độc lập và Type B cho đập tràn đỉnh hẹp).
  • Bước 2: Xác định các thông số hình học tiêu chuẩn ($P, W_u, W_i, W_o, B, N$) theo dải tối ưu đã xác lập.
  • Bước 3: Xác định trạng thái chảy qua các chỉ số phân giới ($H_0/W_o = 0.5$; $H_0/P_H = 0.7$; $h_n/H_n$).
  • Bước 4: Tính toán định lượng lưu lượng tháo qua hệ thống phương pháp hồi quy giải tích chuẩn tắc với độ chính xác cao.

5. Định hướng nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Chiến lược nghiên cứu dài hạn tập trung vào:

  • Phát triển mô hình thủy lực tràn PKW ứng phó với bài toán bùn cát lơ lửng và dòng lũ bùn đá vùng núi cao.
  • Nghiên cứu kết cấu bê tông cốt thép dự ứng lực đúc sẵn dạng module lắp ghép cho phím tràn PKW nhằm rút ngắn thời gian thi công trong điều kiện dẫn dòng thi công phức tạp.
  • Tích hợp thuật toán trí tuệ nhân tạo (AI/Machine Learning) kết hợp mạng thần kinh nhân tạo (ANN/MLP) để tối ưu hóa đa mục tiêu (thủy lực - kết cấu - chi phí) theo thời gian thực.

Kết luận

Luận án tiến sĩ kỹ thuật của tác giả Đoàn Thị Minh Yến đã giải quyết xuất sắc và trọn vẹn bài toán xác định khả năng tháo của tràn phím Piano (PKW) khi xét đến tương tác phức tạp của mực nước hạ lưu, ghi dấu những đóng góp then chốt:

  1. Xác lập một cách khoa học các thông số kích thước của đơn vị tràn PKW tiêu chuẩn tối ưu về thủy lực và kinh tế ($P/W_u = 0.5 \div 1.3$; $W_i/W_o = 1.2 \div 1.5$; $N = 4 \div 6$).
  2. Phát hiện và định lượng hóa chính xác 3 ranh giới thủy động lực học cốt lõi: ranh giới chảy đầy phím ra ($H_0/W_o = 0.5$), ranh giới ảnh hưởng đáy kênh hạ lưu ($H_0/P_H = 0.7$), và ranh giới ngập âm ($h_n/H_n = -0.20 \div 0.0$).
  3. Thiết lập hệ thống công thức giải tích thực nghiệm chuẩn mực xác định hệ số tháo tự do $C_{d0}$, hệ số địa hình đáy $k_H$, và hệ số ngập toàn phần $\sigma_n$, kiểm soát sai số tính toán dưới $5% \div 10%$.
  4. Khám phá quy luật suy giảm hệ số hiệu quả tháo ($r = Q_P/Q_O$) từ giá trị $4.0 \div 5.0$ ở cột nước thấp tiệm cận về $1.0$ ở cột nước cao, làm sáng tỏ bản chất chuyển tiếp cơ chế thủy lực của đập tràn phím không gian.
  5. Ứng dụng thành công và kiểm chứng vững chắc qua các đại công trình thủy lợi, thủy điện tầm cỡ quốc gia và quốc tế như đập dâng Văn Phong, hồ chứa Ngàn Trươi, thủy điện Xuân Minh, Đăk Mi.
  6. Mở ra ba nhánh nghiên cứu học thuật mới đầy triển vọng về thủy động lực học công trình tháo 3D, công nghệ tiêu năng đáy hạ du phân tán, và giải pháp đúc sẵn module PKW trong kỷ nguyên ứng phó biến đổi khí hậu toàn cầu.