Tổng quan về luận án

Đường tràn ngang là công trình tháo lũ hợp lý của đập vật liệu địa phương khi địa hình hẹp không cho phép bố trí tràn dọc, và dòng chảy trong máng tràn bên của nó là bài toán tiêu biểu cho dòng biến lượng (Spatially Varied Flow – SVF) — dòng chảy có lưu lượng thay đổi dọc theo chiều dòng chính. Luận án tiến sĩ thuộc chuyên ngành Kỹ thuật xây dựng công trình thủy này xuất phát từ một research gap được xác định rõ: các phương trình SVF hiện có từ Hinds (1926), Camp (1940), Keulegan (1952) đến Chow (1969) và Konovalov (1937) đều "được thiết lập khi coi dòng chảy trong máng tràn bên là dòng ổn định có lưu lượng tăng dần theo chiều dòng chảy", trong khi bản chất hiện tượng là dòng không ổn định. Đồng thời, hệ phương trình Saint Venant kinh điển (1871) "chưa xét đến lực tác động của dòng gia nhập hoặc phân tán và tổn thất năng lượng do dòng chảy bên hay lực quán tính của dòng chảy cong".

Từ khoảng trống đó, luận án đặt bốn mục tiêu được đánh số cụ thể: (1) thiết lập hệ phương trình vi phân phi tuyến của SVF không ổn định 1D có kể đến các ngoại lực khác trọng lực như lực Coriolis, lực quán tính ly tâm của dòng chảy cong và lực do dòng chảy bên; (2) tuyến tính hóa hệ phương trình bằng phương pháp sai phân và thiết lập chương trình tính; (3) xác định các hệ số trong phương trình áp dụng cho máng tràn bên; (4) xác định một số đặc trưng thủy động lực học trong máng tràn bên. Phạm vi nghiên cứu được giới hạn tường minh: dòng chảy 1D trong máng tràn bên có lưu lượng gia nhập từ một cạnh, độ dốc máng thoải (S₀ < Sc) nối tiếp bằng dốc nước (S₀ > Sc), không kể đến trộn khí.

Đóng góp trung tâm là hệ phương trình vi phân tổng quát (2.28) — một dạng phương trình Saint Venant mở rộng cho SVF không ổn định — cùng các hệ số phân bố lưu tốc thực nghiệm α₀ = 1,41 và α = 2,32, và công cụ mô phỏng số USVF1D giải hệ bằng sơ đồ sai phân ẩn 4 điểm Preissmann.

Literature Review và Positioning

Luận án tổng hợp có hệ thống hơn một thế kỷ nghiên cứu SVF theo ba dòng chính. Dòng thứ nhất là các phương trình động lượng: Hinds (1926) đề xuất phương trình đơn giản đầu tiên dh/dx = f(v, q, Q) với giả thiết lưu lượng tăng tuyến tính; Camp (1940) bổ sung thành phần ma sát; Keulegan (1952) ứng dụng phương pháp điểm suy biến cho máng lăng trụ chữ nhật; Chow (1969) mở rộng cho lòng dẫn có lưu lượng thay đổi với giả thiết dòng gia nhập vuông góc. Dòng thứ hai xuất phát từ bảo toàn năng lượng: Konovalov (1937) đưa ra phương trình (1.8) — "phương trình tổng quát nhất của dạng phương trình SVF ổn định" xét cả nhập lưu và phân lưu, kênh phi lăng trụ, nhưng "chưa xét đến lực quán tính của dòng gia nhập". Dòng thứ ba là các nghiên cứu trong nước: Nguyễn Văn Cung (1964) và Hoàng Tư An (1987) với phương trình chuyển động cho dòng không gian.

Luận án chỉ ra mâu thuẫn học thuật đáng chú ý về vai trò của hệ số phân bố lưu tốc: McCool (1967), thí nghiệm trên kênh dài 125 m với 5 cấp lưu lượng 141,5–1132,67 ℓ/s, kết luận SVF không ảnh hưởng đến α₀ và α; ngược lại Kouchakzadeh và cộng sự (2002), kiểm chứng trên số liệu của Gill (1977), chứng minh rằng khi coi α₀ = 1 thì "đường mặt nước tính toán thấp hơn thực đo với sai số tương đối lớn" và khi α₀ ≈ 1,5 kết quả mới tiệm cận thực đo. Tương tự, Kiselev cho rằng độ dốc ma sát Sf có thể bỏ qua, trong khi Gill (1977) và Mohammadi (2005) đều ghi nhận sức cản thủy lực trong SVF tăng đáng kể so với dòng đều. Luận án định vị mình ở giao điểm của các tranh luận này: kế thừa khung Saint Venant nhưng bổ sung tường minh các lực quán tính bị bỏ qua và xác định thực nghiệm các hệ số phân bố lưu tốc thay vì lấy xấp xỉ 1.

So sánh quốc tế được thực hiện với các nghiên cứu mô hình vật lý tiêu biểu: Lucas và cộng sự (2015) tại ETH Zurich với ba đập Trangslet, Karahnjukar (tỷ lệ 1/45) và Lyssbach (tỷ lệ 1/16); Gardarsson và cộng sự (2015) đối chiếu quan trắc nguyên hình đập Karahnjukar (Iceland, đập cao 198 m, tràn dài 140 m) với mô hình, phát hiện sai khác khoảng 10% giữa mực nước mô hình và thực đo; Maradjieva và Kazakov (2007) kiểm chứng trên mô hình 1/50 đập Markieh (Syria) với Q₀,₁% = 665 m³/s và Q₀,₀₁% = 1000 m³/s.

Đó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 trực tiếp hai hệ lý thuyết có tên tuổi. Thứ nhất, hệ phương trình Saint Venant (Barre de Saint Venant, 1871) được mở rộng thành hệ (2.28) bằng cách đưa vào lực quán tính của dòng chảy bên (thành phần aₓ = (v − vℓₓ)/v · ∂Q/∂t dạng suy rộng) và lực quán tính ly tâm của dòng chảy cong (Fcen = ρu²/r·dA·cosφ). Thứ hai, lý thuyết SVF của trường phái Konovalov–Kiselev được tổng quát hóa: khi xét dòng ổn định, hệ (2.28) suy biến về đúng dạng (1.8) của Konovalov nhưng với hệ số k khác biệt (k = α + kℓβ₀(1 − n₀)), chứng minh tính kế thừa và tính vượt trội của khung mới. Đây là đóng góp kiểu "generalization" — các phương trình kinh điển trở thành trường hợp riêng của phương trình luận án.

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

Khung phân tích tích hợp ba nền tảng: (1) cơ học vật thể có khối lượng thay đổi (Meshchersky) làm cơ sở nhận thức luận cho SVF; (2) hệ phương trình Navier–Stokes làm phương trình xuất phát, được tích phân trực tiếp trên mặt cắt ướt để hạ chiều từ 3D về 1D; (3) lý thuyết sai phân hữu hạn với sơ đồ ẩn Preissmann. Điều kiện biên giới của lý thuyết được phát biểu rõ: chất lỏng không nén được, dòng biến đổi chậm, áp suất thủy tĩnh, bỏ qua trộn khí và bồi xói — tức khung chỉ áp dụng cho máng thoải nối dốc nước, không cho dòng biến đổi gấp.

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

Thiết kế nghiên cứu

Nghiên cứu theo paradigm thực chứng (positivism) với thiết kế kết hợp bốn phương pháp được khai báo tường minh: phân tích – tổng hợp lý thuyết, phân tích thứ nguyên (xây dựng công thức thực nghiệm từ số liệu thực đo), giải tích (tích phân trực tiếp hệ Navier–Stokes) và giải số. Ba phương pháp bổ trợ gồm điều tra khảo sát, kế thừa chọn lọc và phương pháp chuyên gia. Mẫu khảo sát công trình gồm 16 hồ chứa trong nước và 6 công trình quốc tế, với dải thông số rộng: lưu lượng xả từ 37,5 m³/s (Triệu Thượng 1) đến 14.000 m³/s (Hoover, Mỹ); độ dốc máng 0–14%; chiều dài ngưỡng tràn 20–140 m.

Quy trình nghiên cứu rigorous

Việc thiết lập phương trình tuân thủ trình tự chặt chẽ: xác định các lực tác dụng (trọng lực Gx = γS₀dA, lực ma sát theo định luật Newton, lực ly tâm, lực dòng bên), chiếu lên phương chuyển động, tích phân trên toàn mặt cắt với phép thay lưu tốc điểm u bằng lưu tốc trung bình v thông qua các hệ số α₀, α. Tính hợp lệ của mô hình được kiểm soát bằng ba lớp: kiểm tra suy biến về phương trình kinh điển (construct validity), điều kiện ổn định sơ đồ sai phân với trọng số 0,50 ≤ θQ ≤ 1,00 và 0,66 ≤ θZ ≤ 1,00 (internal validity), và đối chiếu với số liệu thực nghiệm của các nghiên cứu trước (external validity).

Data và phân tích

Hệ phi tuyến (2.28) được số hóa bằng sơ đồ ẩn 4 điểm Preissmann — cùng họ thuật giải với các mô hình thương mại MIKE11 (Viện Thủy lực Đan Mạch), HEC-RAS (Hoa Kỳ) và các mô hình trong nước VRSAP, SAL/DELTA, HYDROGIS — cho ra hệ phương trình đại số tuyến tính dạng băng, giải bằng phương pháp khử đuổi. Lựa chọn này được biện luận: độ chính xác bậc 2, thuật toán đơn giản hơn phần tử hữu hạn nhưng đủ tin cậy cho bài toán một đoạn máng đơn. Toàn bộ được đóng gói thành công cụ mô phỏng số USVF1D.

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

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

Thứ nhất, hệ phương trình vi phân tổng quát (2.28) cho SVF không ổn định — dạng Saint Venant mở rộng — là đóng góp lý thuyết trung tâm, lần đầu kể đồng thời lực dòng chảy bên và lực dòng chảy cong trong khung 1D không ổn định.

Thứ hai, các hệ số phân bố lưu tốc của SVF trong máng tràn bên được lượng hóa: α₀ = 1,41 và α = 2,32. Đây là phát hiện phản trực giác so với thực hành thiết kế thông thường (lấy α ≈ 1,05–1,10 theo Cầm, tối đa 1,20 cho kênh thông thường theo Chow): giá trị α = 2,32 vượt cả mức 2,00 của sông vùng ngập lũ, chứng tỏ "việc coi α₀ hay α xấp xỉ 1 sẽ gây sai số lớn khi tính toán đường mặt nước". Kết quả này nhất quán với Kouchakzadeh (α₀ ≈ 1,5) và bác bỏ kết luận của McCool (1967).

Thứ ba, luận án thiết lập công thức xác định chiều sâu dòng chảy phía thành máng đối diện (3.6) — đại lượng mà các phương pháp cũ bỏ qua dù "mực nước trong máng phía thành đối diện luôn cao hơn mực nước trung bình" dưới lực đẩy của dòng gia nhập — và công thức chiều sâu dòng chảy cuối máng (3.17), tinh chỉnh các khuyến nghị trước đó vốn phân tán từ 1,1hc (Kiselev) đến 3hc/2 (Beij, mặt cắt chữ nhật).

Thứ tư, công cụ USVF1D mô phỏng được các đặc trưng thủy động lực học của SVF không ổn định, khắc phục hạn chế của lời giải Saint Venant kinh điển vốn "chưa thể giải quyết được chính xác bài toán dòng không ổn định có lưu lượng thay đổi theo không gian".

Implications đa chiều

Về lý thuyết, kết quả làm phong phú lý thuyết SVF và lý thuyết dòng không ổn định trong lòng dẫn hở. Về thực tiễn, luận án "cung cấp công cụ tính toán thủy lực máng tràn bên cho kết quả chính xác hơn so với phương pháp đã có, cho phép xác định quy mô công trình phù hợp với thực tế" — có ý nghĩa trực tiếp với hàng chục hồ chứa Việt Nam dùng tràn ngang (Đá Bàn 1509,1 m³/s, Việt An 543 m³/s, Mỹ Bình 342 m³/s...) và với công tác thiết kế tràn sự cố đảm bảo an toàn đập. Về chính sách, kết quả hỗ trợ chuẩn hóa quy trình tính toán thủy lực trong tiêu chuẩn thiết kế công trình tháo lũ, đặc biệt cho các hồ chứa lưu vực nhỏ dưới 100 km² — nhóm chiếm đa số công trình tràn ngang trong nước.

Limitations và Future Research

Luận án tự giới hạn trung thực ở bốn điểm: (1) không xét hiện tượng trộn khí, dù Hinds (1926) đã ghi nhận "không khí liên tục bị cuốn vào dòng chảy làm tăng sự xáo trộn bề mặt tới gần cuối máng"; (2) mô hình 1D không tái hiện được các dòng xoắn 3D quy mô lớn và xoắn thứ cấp mà Khiadani và cộng sự đã đo đạc chi tiết; (3) phạm vi chỉ áp dụng cho máng thoải nối dốc nước với dòng gia nhập từ một cạnh; (4) các hệ số α₀, α xác định trong giới hạn điều kiện thí nghiệm cụ thể, chưa phổ quát cho mọi hình dạng mặt cắt.

Hướng nghiên cứu tiếp theo gồm: mở rộng phương trình cho dòng hai pha nước–khí; kiểm chứng bổ sung trên nguyên hình theo cách tiếp cận của Gardarsson (2015); mở rộng cho máng nhận nước từ hai cạnh và đầu tràn (dạng chữ U kiểu Flatiron); tích hợp USVF1D với các mô hình hệ thống sông; và nghiên cứu ảnh hưởng của tràn zigzag/piano đến cấu trúc dòng gia nhập.

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

Về học thuật, việc tổng quát hóa Saint Venant cho dòng biến lượng mở kênh trích dẫn cho cả cộng đồng mô hình hóa thủy lực 1D — nhóm dùng MIKE11, HEC-RAS, VRSAP — khi xử lý nhập lưu phân bố. Về công nghiệp, công cụ áp dụng trực tiếp cho khảo sát an toàn đập và thiết kế mới trong bối cảnh Việt Nam có hàng nghìn hồ chứa nhỏ cần tràn sự cố; các khuyến nghị định lượng kiểu Nguyễn Chiến (Bđầu ≈ 0,10–0,15 chiều dài ngưỡng, độ dốc máng ≤ 5%) nay có thể được kiểm tra lại bằng mô hình chính xác hơn. Về quốc tế, bài toán tương đồng với các công trình từ Hoover (Mỹ) đến Karahnjukar (Iceland), nên phương pháp có tính chuyển giao toàn cầu.

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

Nghiên cứu sinh thủy lực nhận được một khung thiết lập phương trình mẫu mực từ Navier–Stokes xuống 1D cùng gap còn mở về trộn khí và dòng 3D. Giới học thuật cao cấp có cơ sở đánh giá lại giả thiết phân bố lưu tốc đồng nhất trong toàn bộ họ phương trình SVF. Kỹ sư tư vấn thiết kế được cung cấp công cụ USVF1D và hai công thức thực hành (3.6), (3.17) tính chiều sâu thành đối diện và cuối máng. Cơ quan quản lý an toàn đập có căn cứ định lượng để rà soát năng lực tháo lũ của các tràn ngang hiện hữu.

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

1. Đóng góp lý thuyết độc đáo nhất? Mở rộng hệ Saint Venant (1871) bằng lực quán tính dòng chảy bên và dòng chảy cong, với tính chất suy biến kép: về Konovalov (1937) khi dòng ổn định, về Saint Venant kinh điển khi k = α, k₀ = β₀.

2. Đổi mới phương pháp? So với Gill (1977) tìm nghiệm xấp xỉ đại số cho dòng ổn định và Lucas (2015) dựa vào mô hình vật lý tỷ lệ 1/45, luận án giải trực tiếp hệ phi tuyến không ổn định bằng sơ đồ Preissmann có trọng số kiểm soát ổn định.

3. Phát hiện bất ngờ nhất? α = 2,32 — gấp đôi giá trị thiết kế thông dụng, vượt mức sông ngập lũ (2,00 theo Chow), có data thí nghiệm hậu thuẫn.

4. Khả năng tái lập? Có: thuật toán khử đuổi, sơ đồ sai phân (2.30)–(2.48) và điều kiện trọng số được công bố đầy đủ trong Chương 2 và 4.

5. Chương trình nghiên cứu dài hạn? Đã phác thảo qua các hướng mở về trộn khí, mô hình 2D/3D và kiểm chứng nguyên hình.

Kết luận

Luận án để lại sáu đóng góp đo lường được: (1) hệ phương trình Saint Venant mở rộng (2.28) cho SVF không ổn định 1D; (2) bộ hệ số phân bố lưu tốc α₀ = 1,41, α = 2,32 xác định thực nghiệm; (3) công thức chiều sâu dòng chảy phía thành đối diện (3.6); (4) công thức chiều sâu cuối máng (3.17); (5) hệ đại số tuyến tính hóa theo sơ đồ Preissmann; (6) công cụ mô phỏng số USVF1D. Bằng việc chuyển bài toán máng tràn bên từ khung "ổn định gần đúng" sang khung "không ổn định đúng bản chất", nghiên cứu không chỉ làm phong phú lý thuyết dòng biến lượng mà còn mở ra các nhánh nghiên cứu mới về dòng hai pha, mô hình hóa đa chiều và kiểm chứng nguyên hình — với tính ứng dụng trực tiếp cho an toàn đập ở Việt Nam và tính tương thích với thực tiễn thiết kế tràn ngang toàn cầu từ Hoover đến Karahnjukar.