Tổng quan về luận án

Nghiên cứu các đặc trưng thủy lực ở đập tràn thực dụng có tường ngực biên cong của NCS. Đỗ Ngọc Ánh (Viện Khoa học Thủy lợi Việt Nam, 2017) ra đời trong bối cảnh hệ thống hồ đập Việt Nam đang đối mặt với áp lực kép: yêu cầu cập nhật năng lực tháo lũ theo QCVN 04-05:2012 và tác động ngày càng gia tăng của biến đổi khí hậu. Theo số liệu thống kê trong luận án, Việt Nam có hơn 7.000 hồ đập với dung tích tổng cộng trên 37 tỷ m³, trong đó có 675 đập lớn, phần lớn được thiết kế theo tiêu chuẩn cũ. Dự báo cho thấy lưu lượng thiết kế có thể tăng khoảng 2÷3% vào năm 2020 và 8÷10% vào năm 2099 so với tính toán trước đây, đặt ra bài toán cấp bách về nâng cao năng lực tháo lũ cho các công trình hiện hữu.

Research gap được xác định rõ ràng: Dù đập tràn thực dụng có tường ngực biên cong đã được ứng dụng tại các công trình trọng điểm như thủy điện Sơn La (6 khoang xả mặt, chiều rộng 15 m/khoang, ngưỡng cao trình 197,8 m) và quốc tế như Thana Plaun – Ấn Độ (5 khoang, 11 m × 15,2 m), Muskrat Falls – Canada (PMF Q = 22.100 m³/s), Mangla – Iran (PMF Q = 28.600 m³/s), nhưng chưa có công trình khoa học nào tại Việt Nam nghiên cứu hệ thống về chế độ thủy lực cho loại đập này. Cụ thể, bốn tồn tại chưa được giải quyết gồm: (1) chưa phân định chế độ chảy qua đập; (2) hệ số lưu lượng µ chỉ được cho theo dải giá trị mà không chỉ rõ yếu tố phụ thuộc; (3) chưa có hướng dẫn xác định đường mặt nước, vận tốc, áp suất trên mặt tràn; (4) chưa chỉ ra cách xác định áp suất nhỏ nhất trên mặt tràn và biên tường ngực.

Research questions được đặt ra gồm: (RQ1) Ranh giới chế độ dòng chảy hở và có áp tại đập tràn thực dụng có tường ngực biên cong xảy ra ở ngưỡng H/D nào? (RQ2) Hệ số lưu lượng µ phụ thuộc vào các yếu tố nào và có thể xây dựng công thức tổng quát không? (RQ3) Phân bố đường mặt nước, vận tốc và áp suất trên mặt tràn tuân theo quy luật gì? (RQ4) Áp suất trên tường ngực biên cong được xác định ra sao?

Theoretical framework vận dụng lý thuyết đồng dạng thủy động lực học (similitude theory), phương trình Navier–Stokes dạng Reynolds-Averaged (RANS), lý thuyết dòng chảy thế (potential flow theory) và phương pháp quy hoạch thực nghiệm (Design of Experiments – DoE) với hồi quy đa biến. Nghiên cứu được thực hiện trên 04 mô hình thủy lực vật lý tại Trung tâm Nghiên cứu Thủy lực, kết hợp hai dạng mặt cắt thông dụng nhất là Creager–Ophixerop và WES trong phạm vi H/Hd ≤ 1,5 (tương đương H/D ≤ 3,0) – là vùng vận hành thực tế phổ biến nhất của các đập tràn lớn Việt Nam.

Literature Review và Positioning

Lịch sử nghiên cứu đập tràn thực dụng hình cong trải dài từ Bazin (1886) với đập thành mỏng, qua Pascal (1889), De Marchi (1928), Creager (1929, 1945), Scimemi (1930), đến chuẩn Waterways Experiment Station – WES (1952). Hai luồng nghiên cứu chính tồn tại song song: luồng Xô viết với mặt cắt Creager–Ophixerop ứng dụng rộng ở Nga, Đông Âu, Việt Nam, Trung Quốc; và luồng Hoa Kỳ với mặt cắt WES được chuẩn hóa bởi U.S. Army Corps of Engineers (USACE) và U.S. Bureau of Reclamation (USBR).

Trong lĩnh vực đập có tường ngực, Ấn Độ là quốc gia dẫn đầu với 17 công trình được thống kê chi tiết trong luận án (Chamara I, Chamara III, Dhaulignaga, Kurichu, Nathpa Jhakri, Nimoobazgo, Pandoh, Parbati II, Parbati III, Ranganadi, Sawa II, Subanairi, Tala, Taagta III, IV, V, Uri II), với hệ số lưu lượng Cd thực đo trong khoảng 0,65÷0,92 và lưu lượng đơn vị từ 102,77 đến 338,16 m²/ms. Tiêu chuẩn Ấn Độ đã được cập nhật năm 2010 thay thế tiêu chuẩn năm 1973, với nhận xét rõ ràng rằng mặt tràn theo tiêu chuẩn 1973 "có độ dốc rất lớn, dễ gây xâm thực" – một hạn chế nghiêm trọng trong điều kiện cột nước cao.

Mâu thuẫn trong literature: Theo tài liệu Ấn Độ ([29], [30]), hệ số lưu lượng µ chỉ được xác định theo tỉ lệ H/Hd (với các giá trị 0,696; 0,769; 0,797; 0,829 tương ứng H/Hd = 0,8; 1,0; 1,15; 1,33), trong khi các tài liệu thiết kế sơ bộ ([45]) cho phạm vi µ = 0,74÷0,82 khi H/D = 2,0÷2,4 và µ = 0,83÷0,93 khi H/D > 2,4 – hai cách tiếp cận sử dụng biến độc lập khác nhau và không thống nhất về ý nghĩa vật lý. Trong khi đó, nghiên cứu CFD của Savage và Johnson (2001) sử dụng Flow-3D cho thấy sai số lưu lượng chỉ 1% so với mô hình vật lý khi H ≤ 1,2Hd, nhưng các kết quả của Ho, Cooper, Riddette và Donohoo (2006) lại chỉ ra sai lệch cao hơn 5% ở cột nước bằng hoặc vượt mức thiết kế – chứng tỏ mô hình số chưa thay thế hoàn toàn thực nghiệm vật lý ở điều kiện vận hành cực trị.

Positioning trong literature: So sánh với nghiên cứu quốc tế của Azevedo (2002) dùng FEM với sai số 0,5÷3% và nghiên cứu của Daneshfaraz et al. (2013) dùng Flow-3D cho đập Siah-Bishe, luận án Đỗ Ngọc Ánh đóng góp ở chỗ xây dựng quan hệ thực nghiệm tổng quát có thể sử dụng trực tiếp trong thiết kế mà không phụ thuộc vào phần mềm số đắt tiền – điều đặc biệt có giá trị trong điều kiện thực tiễn kỹ thuật Việt Nam. Nghiên cứu Hoàng Văn Tần (1999) về xả kết hợp tại Việt Nam đã đặt nền móng, nhưng chỉ xét đập tràn đỉnh rộng hoặc thực dụng thông thường, chưa đề cập đến trường hợp tường ngực biên cong.

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

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

Nghiên cứu mở rộng lý thuyết đồng dạng thủy động lực học (Buckingham Pi theorem) vào bài toán cụ thể của đập tràn có tường ngực, với hàm phi thứ nguyên φ(Fr, Eu, Re, Sh) = 0 được khai triển để xác định tiêu chuẩn tương tự trọng lực Froude làm chủ đạo, đồng thời xác định ngưỡng số Reynolds giới hạn (Re_gh) để đảm bảo kết quả mô hình đại diện cho nguyên hình. Đây là bước kiểm chứng nền tảng trước khi xây dựng công thức thực nghiệm.

Luận án mở rộng lý thuyết dòng chảy thế (potential flow theory) theo hướng bán thực nghiệm – tiếp nối Nguyễn Văn Mạo (1987) và Nguyễn Danh Oanh (2003) – bằng cách phát triển phương trình tính hệ số giảm áp Cp cho biên trên tường ngực biên cong, từ đó xác định áp suất nhỏ nhất tại điểm nguy hiểm nhất theo công thức:

p_i,min/γ = H_i – C_p,max × V_k²/2g (phương trình 1.27)

trong đó V_k là vận tốc trung bình tại mặt cắt cuối cửa vào, H_i là độ sâu tính từ mực nước thượng lưu. Đây là công cụ định lượng hoàn toàn mới cho loại kết cấu này tại Việt Nam.

Conceptual framework của luận án tích hợp ba lý thuyết: (1) Lý thuyết đồng dạng thủy động lực học (similitude theory) để thiết kế mô hình vật lý; (2) Lý thuyết dòng thế phẳng (potential flow) để xây dựng phương trình áp suất trên tường ngực; (3) Phương pháp quy hoạch thực nghiệm và hồi quy đa biến (DoE + multiple regression) để xây dựng các hàm thực nghiệm tổng quát. Sự tích hợp này tạo ra khung phân tích không chỉ mô tả hiện tượng mà còn cung cấp công cụ tính toán trực tiếp phục vụ thiết kế.

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

Khung phân tích độc đáo của nghiên cứu là phương pháp xác định chế độ dòng chảy dựa trên ngưỡng tỉ số H/D thay vì chỉ dùng H/Hd như các tài liệu quốc tế trước đây. Cụ thể: (i) vẽ đường cong mặt tràn và tường ngực trên cùng hệ tọa độ; (ii) xây dựng đường mặt nước lý thuyết ứng với các tỉ lệ H/D khác nhau; (iii) xác định ngưỡng chuyển tiếp. Kết quả lý thuyết cho thấy H/D ≤ 0,88÷1,23 tương ứng chảy tự do, H/D = 1,41 là trạng thái phân giới và H/D ≥ 1,58 là chảy có áp hoàn toàn – cần được kiểm chứng thực nghiệm.

Điều kiện biên (boundary conditions) của nghiên cứu được xác định rõ: mặt tràn dạng WES và Creager–Ophixerop; tỉ số H/D trong khoảng 1,6÷3,0 (chảy có áp); cửa van mở hoàn toàn hoặc dòng chảy tự do phía sau tường ngực; không xét trường hợp ngập hạ lưu và hàm khí.

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

Thiết kế nghiên cứu

Triết lý nghiên cứu (research philosophy) theo chủ nghĩa thực chứng (positivism) kết hợp bán thực nghiệm (semi-empirical approach): xuất phát từ lý thuyết cơ học chất lỏng để xác định cấu trúc hàm số, sau đó dùng thực nghiệm để xác định các hệ số. Thiết kế nghiên cứu là mô hình vật lý thủy lực kết hợp quy hoạch thực nghiệm, sử dụng tiêu chuẩn tương tự Froude làm chủ đạo do trọng lực chiếm ưu thế trong dòng chảy qua đập tràn.

Bốn mô hình thực nghiệm được thiết kế và thi công, bao gồm: Mô hình 1 (mặt cắt Creager–Ophixerop, đơn khoang, khảo sát D = 5, 6, 7, 8 m nguyên hình); Mô hình 2 (mặt cắt WES, đơn khoang, tường ngực biên cong elip); Mô hình 3 (mặt cắt Creager–Ophixerop, đơn khoang, tường ngực dạng khác); Mô hình 4 (mặt cắt WES, đa khoang có trụ pin, mô phỏng điều kiện thực tế gần hơn). Mô hình 2 và 4 cho phép so sánh trực tiếp ảnh hưởng của trụ pin đến các đặc trưng thủy lực.

Quy trình nghiên cứu rigorous

Chiến lược lấy mẫu: Các trường hợp thí nghiệm được thiết kế theo phương pháp quy hoạch thực nghiệm đa yếu tố (factorial design), với biến độc lập chính là tỉ số H/D (biến đặc trưng chế độ chảy) và tỉ số X/D, Y/D (biến đặc trưng vị trí). Phương pháp Holdout được áp dụng để phân chia tập dữ liệu thành tập huấn luyện (training set) và tập kiểm định (validation set), đảm bảo khả năng tổng quát hóa của các công thức thực nghiệm.

Thiết bị đo lường: Máng lường đo lưu lượng; thiết bị đo lưu tốc điện từ; đầu đo áp suất bố trí theo tuyến giữa và tuyến sát trụ pin trên mặt tràn và tường ngực. Triangulation được thực hiện ở cấp độ dữ liệu (so sánh kết quả giữa các mô hình), phương pháp (so sánh thực nghiệm với lý thuyết và CFD) và lý thuyết (đối chiếu với kết quả của Nguyễn Danh Oanh (2003), G. Sliski ([51], [52]), Savage và Johnson (2001)).

Độ tin cậy và độ hợp lệ: Mô hình được kiểm định (model verification) bằng cách so sánh kết quả hệ số lưu lượng khi chảy hở với đập tràn chuẩn – sai số thực nghiệm so với lý thuyết trong khoảng ±5%, phù hợp với kết quả quốc tế. Hệ số tương quan bội S và các chỉ số SSR, SST, SSE được tính toán đầy đủ để đánh giá độ phù hợp của mô hình hồi quy. Công thức thực nghiệm được kiểm định lại với bộ số liệu độc lập (công trình hồ chứa Bản Lải, tỉnh Lạng Sơn).

Data và phân tích

Dữ liệu thực nghiệm bao gồm các bộ số liệu về: (i) chế độ chảy (giới hạn H/D); (ii) lưu lượng và hệ số lưu lượng µ; (iii) tọa độ đường mặt nước không thứ nguyên Y/D = f(X/D, H/D); (iv) phân bố vận tốc dọc mặt tràn; (v) phân bố áp suất hp/Hd = f(X/Hd, H/Hd) và hp/D = f(X/D, H/D) trên mặt tràn và tường ngực. Phân tích sử dụng phương pháp hồi quy đa biến phi tuyến (nonlinear multiple regression) với phần mềm thống kê, xây dựng các hàm thực nghiệm đa biến. Kiểm tra độ bền vững (robustness check) được thực hiện bằng cách so sánh kết quả công thức 3.7 (hệ số µ) với số liệu thực đo tại Bản Lải và số liệu quốc tế.

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

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

Phát hiện 1 – Ngưỡng chế độ chảy: Thực nghiệm xác nhận và làm rõ ngưỡng chuyển chế độ dòng chảy. Cụ thể: chảy tự do khi H/D < 1,23÷1,41 (tùy dạng mặt cắt); trạng thái quá độ khi H/D ≈ 1,41÷1,58; chảy có áp hoàn toàn khi H/D ≥ 1,58÷1,74. Phát hiện này lấp đầy khoảng trống mà tài liệu [34], [35] chỉ nêu sơ lược "H/D > 1,5÷1,7" mà không có cơ sở thực nghiệm chi tiết.

Phát hiện 2 – Công thức hệ số lưu lượng tổng quát: Hệ số lưu lượng µ trong chế độ chảy có áp phụ thuộc đồng thời vào H/D và hình dạng mặt cắt (WES hoặc Creager–Ophixerop), không thể chỉ dùng H/Hd như các tài liệu trước đây ([29], [30]). Công thức mới (phương trình 3.7 trong luận án) cho phép tính µ trực tiếp trong phạm vi H/D = 1,6÷3,0 với sai số được kiểm soát. So sánh với Bảng 1.2 (số liệu Ấn Độ): µ = 0,83÷0,93 ở vùng H/D > 2,4 – kết quả luận án nằm trong dải này nhưng phân biệt rõ ảnh hưởng của dạng mặt cắt, điều mà tài liệu Ấn Độ không làm.

Phát hiện 3 – Bảng tọa độ không thứ nguyên đường mặt nước: Lần đầu tiên xây dựng bảng tọa độ không thứ nguyên Y/D = f(X/D, H/D) riêng cho từng dạng mặt cắt (WES – Bảng 3.7; Creager–Ophixerop – Bảng 3.8) ở chế độ chảy có áp. Kết quả thực nghiệm cho thấy đường mặt nước chảy có áp có hình dạng khác hẳn so với chảy tự do, và sự khác biệt này tăng mạnh khi H/D tăng – một kết quả không trực giác (counter-intuitive) bởi nhiều kỹ sư thiết kế vẫn sử dụng công thức chảy tự do để ước tính đường mặt nước ngay cả khi có tường ngực.

Phát hiện 4 – Phân bố áp suất và xác định C_p,max: Thực nghiệm xác định được hệ số giảm áp C_p,max cho tường ngực biên cong trên cả hai dạng mặt cắt (Bảng 3.9 – WES; Bảng 3.10 – Creager–Ophixerop), cho phép tính áp suất nhỏ nhất tại điểm nguy hiểm nhất trên tường ngực theo phương trình 1.27. Đây là đóng góp trực tiếp vào tiêu chí an toàn kết cấu. Kết quả cũng chỉ ra rằng áp suất chân không (negative pressure) có thể xuất hiện ở tường ngực khi H/D vượt ngưỡng nhất định – cần được kiểm soát để tránh nguy cơ xâm thực (cavitation).

Phát hiện 5 – Ảnh hưởng của trụ pin: So sánh Mô hình 2 và Mô hình 4 ở cùng cột nước H = 17,2 m và H = 19,05 m cho thấy sự hiện diện của trụ pin làm thay đổi đáng kể phân bố vận tốc và áp suất theo phương ngang – kết quả này cho thấy các công thức được xây dựng riêng cho trường hợp có và không có trụ pin.

Implications đa chiều

Đóng góp lý thuyết: Nghiên cứu bổ sung cơ sở khoa học cho lý thuyết dòng chảy hỗn hợp (mixed flow theory) – trường hợp trung gian giữa chảy hở hoàn toàn và chảy áp lực hoàn toàn – vốn ít được nghiên cứu hơn so với hai trường hợp thuần túy. Đồng thời, phương pháp quy hoạch thực nghiệm kết hợp hồi quy đa biến được chứng minh hiệu quả và có thể áp dụng cho các bài toán thủy lực công trình tương tự (tuynel, cống hộp có mái cong).

Ứng dụng thực tiễn: Quy trình tính toán 4 bước được đề xuất tại Chương 4 (lựa chọn kích thước → kiểm tra khả năng tháo → tính đường mặt nước và vận tốc → xác định áp suất) được kiểm chứng qua ví dụ tính toán thực tế cho đập tràn hồ chứa Bản Lải (Lạng Sơn). Bảng 4.7 trong luận án so sánh kết quả tính theo phương pháp luận án với phương án điều chỉnh dự kiến, cho thấy sự tương đồng cao và chứng minh tính khả dụng của quy trình.

Khuyến nghị chính sách: Kết quả nghiên cứu hỗ trợ trực tiếp việc triển khai QCVN 04-05:2012 – đặc biệt trong đánh giá an toàn và nâng cấp các đập không đạt tiêu chuẩn lũ hiện hành. Các công thức và bảng biểu có thể tích hợp vào tiêu chuẩn thiết kế quốc gia cho đập tràn có tường ngực.

Điều kiện tổng quát hóa (generalizability): Kết quả áp dụng trực tiếp cho mặt tràn WES và Creager–Ophixerop trong phạm vi H/Hd ≤ 1,5, H/D = 1,6÷3,0, tường ngực biên cong elip, cửa van mở hoàn toàn, không có ngập hạ lưu và không có hàm khí. Ngoài phạm vi này cần thực nghiệm bổ sung.

Limitations và Future Research

Hạn chế 1 – Phạm vi thông số: Nghiên cứu giới hạn ở H/D = 1,6÷3,0 và H/Hd ≤ 1,5. Các trường hợp cột nước rất cao (H/D > 3 hay lũ PMF cực trị) hoặc cửa van mở một phần chưa được xét đến, trong khi thực tế vận hành đôi khi xảy ra ở điều kiện này.

Hạn chế 2 – Mô hình vật lý 2D/3D: Phần lớn thực nghiệm thực hiện trên mặt cắt 2D (đơn khoang). Hiệu ứng 3D tại vùng trụ pin, vách bên và cửa vào khoang tràn chỉ được khảo sát một phần qua Mô hình 4 đa khoang, chưa có phân tích đầy đủ về phân bố vận tốc theo phương ngang.

Hạn chế 3 – Hiệu ứng tỷ lệ mô hình (scale effect): Dù đã kiểm tra số Reynolds giới hạn, một số hiệu ứng tỷ lệ liên quan đến nhám bề mặt và sức căng bề mặt không thể loại trừ hoàn toàn ở mô hình thu nhỏ.

Hạn chế 4 – Không xét xâm thực và tiêu năng hạ lưu: Áp suất chân không được xác định nhưng chưa được kết nối với mô hình dự báo xâm thực định lượng, và tiêu năng hạ lưu nằm ngoài phạm vi nghiên cứu.

Hướng nghiên cứu tương lai: (1) Mở rộng phạm vi H/D > 3 và nghiên cứu trường hợp cửa van mở một phần bằng kết hợp Flow-3D và mô hình vật lý; (2) Phát triển mô hình CFD 3D để xác định phân bố áp suất toàn trường theo phương ngang, kiểm chứng với thực nghiệm; (3) Nghiên cứu ảnh hưởng của hình dạng biên tường ngực (tròn, thẳng, đa đoạn) đến C_p,max và ngưỡng chảy có áp; (4) Phát triển tiêu chí chống xâm thực tích hợp dựa trên phân bố áp suất nhỏ nhất từ nghiên cứu này; (5) Nghiên cứu chế độ dòng chảy quá độ (transient flow) khi mực nước thay đổi nhanh trong vận hành lũ thực tế.

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

Tác động học thuật: Nghiên cứu lấp đầy khoảng trống hệ thống trong tài liệu kỹ thuật thủy lợi tiếng Việt và đóng góp vào cơ sở dữ liệu thực nghiệm quốc tế cho loại kết cấu này. Với 04 mô hình vật lý, phạm vi thông số rộng và quy trình kiểm định nghiêm ngặt, nghiên cứu đủ điều kiện để được trích dẫn bởi các nghiên cứu thiết kế và tiêu chuẩn quốc gia, đặc biệt trong bối cảnh Việt Nam đang đẩy mạnh hiện đại hóa hạ tầng thủy lợi đến năm 2030.

Tác động ngành: Tác động trực tiếp đến hai lĩnh vực: (i) thiết kế công trình mới – cho phép tối ưu hóa kích thước cửa van (như thủy điện Sơn La tiết kiệm được chiều cao cửa van từ 30,3 m xuống 12 m nhờ tường ngực); (ii) nâng cấp an toàn đập – cung cấp công cụ tính toán cho hàng trăm công trình cần tăng năng lực tháo lũ theo QCVN 04-05:2012. Với hơn 675 đập lớn tại Việt Nam, tiềm năng ứng dụng là rất lớn, đặc biệt khi kết hợp với yêu cầu tăng dung tích phòng lũ trước biến đổi khí hậu (lưu lượng tăng 8÷10% đến năm 2099).

Tác động chính sách: Kết quả nghiên cứu cung cấp cơ sở khoa học cho Bộ Nông nghiệp và PTNT trong việc cập nhật quy phạm, tiêu chuẩn thiết kế đập tràn có tường ngực, bổ sung vào hệ thống tiêu chuẩn kỹ thuật thủy lợi quốc gia. Ở cấp địa phương, quy trình tính toán 4 bước và bảng biểu phi thứ nguyên có thể được áp dụng ngay trong thẩm tra, thẩm định dự án mà không đòi hỏi phần mềm chuyên dụng đắt tiền.

Ý nghĩa quốc tế: So sánh với kinh nghiệm Ấn Độ (17 công trình), Canada (Muskrat Falls) và Iran (Mangla), nghiên cứu này bổ sung góc nhìn Đông Nam Á – nơi điều kiện thủy văn lũ nhiệt đới có đặc trưng khác biệt – vào kho tàng kiến thức thủy lực công trình toàn cầu.

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

Nghiên cứu sinh và nhà khoa học: Công trình cung cấp cơ sở lý luận và phương pháp nghiên cứu mẫu mực (quy hoạch thực nghiệm + hồi quy đa biến + kiểm định độc lập) cho các nghiên cứu tiếp theo về thủy lực công trình tháo lũ. Các khoảng trống còn lại về H/D > 3,0, hiệu ứng 3D và xâm thực là nguồn đề tài nghiên cứu cho ít nhất 3÷5 luận án tiếp theo.

Giảng viên đại học và chuyên gia thủy lực cao cấp: Khung lý thuyết tích hợp (similitude + potential flow + DoE) và bộ số liệu thực nghiệm đủ chi tiết (4 mô hình, nhiều trường hợp thí nghiệm, bảng số liệu đầy đủ) có giá trị làm tài liệu giảng dạy và tham khảo chuyên sâu cho môn Thủy lực công trình ở bậc sau đại học.

Kỹ sư thiết kế và tư vấn thủy lợi: Quy trình tính toán 4 bước và bộ bảng biểu phi thứ nguyên (Bảng 3.7, 3.8, 3.9, 3.10) có thể áp dụng trực tiếp trong thực hành thiết kế mà không cần thực nghiệm thêm, giúp rút ngắn thời gian và giảm chi phí nghiên cứu. Ví dụ áp dụng cho hồ chứa Bản Lải trong Chương 4 là hướng dẫn thực hành cụ thể.

Cơ quan quản lý nhà nước và nhà hoạch định chính sách: Kết quả nghiên cứu cho thấy đập tràn có tường ngực biên cong có thể tăng dung tích phòng lũ lên đến 3 tỷ m³ (trường hợp Sơn La), đây là thông tin định lượng quan trọng để đánh giá lợi ích kinh tế–xã hội khi lựa chọn giải pháp nâng cấp đập trong các dự án đầu tư công.

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

1. Đóng góp lý thuyết độc đáo nhất là gì? Đóng góp lý thuyết độc đáo nhất là việc phát triển phương pháp xác định hệ số giảm áp C_p,max cho tường ngực biên cong ở đập tràn thực dụng – mở rộng lý thuyết của G. Sliski ([51], [52]) vốn chỉ áp dụng cho cửa vào tuynel – sang bài toán mới là tường ngực đập tràn. Từ đó công thức áp suất nhỏ nhất (phương trình 1.27) được áp dụng trực tiếp trong thiết kế kết cấu chống xâm thực.

2. Đổi mới phương pháp luận so với các nghiên cứu trước? So với nghiên cứu Ấn Độ ([29], [30]) chỉ cung cấp bảng hệ số µ theo H/Hd với số điểm rất hạn chế (4 giá trị), và so với nghiên cứu CFD của Savage và Johnson (2001) đòi hỏi phần mềm Flow-3D đắt tiền, luận án xây dựng công thức thực nghiệm đa biến dạng hàm liên tục với biến H/D, kiểm định bằng phương pháp Holdout và ứng dụng ngay vào công trình thực tế Bản Lải – tạo ra công cụ thiết kế vừa chính xác vừa thực dụng.

3. Phát hiện bất ngờ nhất? Phát hiện bất ngờ nhất là phân bố vận tốc và áp suất ở chế độ chảy có áp khác biệt đáng kể so với chế độ chảy tự do ngay cả ở đoạn phía sau tường ngực (đoạn chảy hở), nghĩa là ảnh hưởng của đoạn có áp kéo dài xuống hạ lưu xa hơn dự kiến. Điều này giải thích tại sao hệ số lưu tốc j xây dựng cho chảy tự do (Hình 1.24) không còn phù hợp khi áp dụng cho đập có tường ngực, và đòi hỏi phải xác định lại j qua thực nghiệm riêng.

4. Giao thức tái lập (replication protocol) được cung cấp không? Có. Chương 2 trình bày đầy đủ: tiêu chuẩn tương tự Froude và ngưỡng Re giới hạn, quy trình thiết kế và lắp đặt 4 mô hình, bố trí thiết bị đo (đầu đo áp suất, thiết bị đo lưu tốc, máng lường), các trường hợp thí nghiệm cụ thể cho từng mô hình (Bảng 2.3, 2.4), phương pháp quy hoạch thực nghiệm và xây dựng hàm hồi quy, tiêu chí đánh giá độ phù hợp (S, SSR, SST, SSE). Đây là mức độ chi tiết đủ để tái lập nghiên cứu.

5. Chương trình nghiên cứu 10 năm? Dựa trên kết quả và hạn chế được chỉ ra, chương trình 10 năm gồm: năm 1÷3 phát triển mô hình CFD 3D kiểm chứng và mở rộng sang H/D > 3; năm 3÷5 nghiên cứu ảnh hưởng hình dạng tường ngực và hiệu ứng 3D trụ pin; năm 5÷7 phát triển tiêu chí xâm thực tích hợp; năm 7÷9 nghiên cứu chế độ vận hành động (transient flow) dưới lũ cực trị; năm 9÷10 xây dựng tiêu chuẩn thiết kế quốc gia và đề xuất tiêu chuẩn ISO khu vực ASEAN cho đập tràn có tường ngực.

Kết luận

Luận án tiến sĩ của Đỗ Ngọc Ánh (2017) tạo ra sáu đóng góp khoa học cụ thể và có thể đo lường:

  1. Phân định chế độ chảy: Xác định ngưỡng H/D = 1,58÷1,74 là ranh giới chảy có áp hoàn toàn – lần đầu tiên được xác định bằng thực nghiệm hệ thống tại Việt Nam cho loại kết cấu này.

  2. Công thức hệ số lưu lượng mới: Phương trình 3.7 tính µ theo H/D trong phạm vi 1,6÷3,0 cho cả hai dạng mặt cắt WES và Creager–Ophixerop, khắc phục hạn chế của tài liệu Ấn Độ ([29]) và công thức thiết kế sơ bộ ([45]).

  3. Bảng tọa độ phi thứ nguyên đường mặt nước và vận tốc: Bộ bảng Y/D = f(X/D, H/D) và hệ số lưu tốc j hiệu chỉnh cho chế độ chảy có áp – công cụ tính toán trực tiếp chưa từng có trong tài liệu tiếng Việt.

  4. Phương pháp xác định C_p,max cho tường ngực biên cong: Mở rộng lý thuyết Sliski sang bài toán tường ngực đập tràn, cho phép tính áp suất nhỏ nhất phục vụ thiết kế chống xâm thực.

  5. Biểu đồ phi thứ nguyên phân bố áp suất: Bộ biểu đồ hp/Hd = f(X/Hd, H/Hd) và hp/D = f(X/D, H/D) cho đoạn chảy tự do và chảy có áp trên cả hai dạng mặt cắt – đủ chi tiết để áp dụng trực tiếp trong thiết kế.

  6. Quy trình tính toán thủy lực 4 bước hoàn chỉnh: Được kiểm chứng qua công trình thực tế hồ chứa Bản Lải, cho phép kỹ sư thiết kế áp dụng ngay mà không cần thực nghiệm bổ sung trong phạm vi thông số đã xác định.

Nghiên cứu mở ra ít nhất ba hướng nghiên cứu mới: (i) thủy lực đập tràn có tường ngực trong điều kiện lũ cực trị (H/D > 3); (ii) tương tác dòng chảy 3D tại vùng trụ pin và vách biên; (iii) tiêu chí xâm thực định lượng cho mặt tràn và tường ngực biên cong. Với nền tảng từ hệ thống hơn 7.000 hồ đập tại Việt Nam cần được đánh giá và nâng cấp theo QCVN 04-05:2012, và trong bối cảnh lưu lượng lũ dự báo tăng 8÷10% đến năm 2099, di sản khoa học của công trình này có ý nghĩa thiết thực kéo dài nhiều thập kỷ cho ngành thủy lợi, thủy điện Việt Nam và khu vực.