Nghiên cứu đặc trưng thủy động lực học của dòng biến lượng không ổn định trong máng tràn bên

Luận án phân tích đặc trưng thủy động lực học của dõng biến lượng không ổn định trong máng tràn bên, cung cấp kiến thức chuyên sâu cho nghiên cứu.

• Trường học: Trường Đại Học Thủy Lợi
• Chuyên ngành: Thủy Lợi
• Đề tài: Đặc Trưng Thủy Động Lực Học Của Dòng Biến Lượng Không Ổn Định Trong Máng Tràn
• Loại tài liệu: Luận Án
• Năm xuất bản: 2023
• Địa điểm: Hà Nội

Mục lục chi tiết

• 1. CHƯƠNG 1: TỔNG QUAN NGHIÊN CỨU VỀ MÁNG TRÀN BÊN VÀ DÒNG BIẾN LƯỢNG
• 1.1. SỰ CHUYỂN ĐỘNG CỦA DÒNG CHẢY
• 1.2. MÁNG TRÀN BÊN
• 1.2.1. Đặc điểm công trình máng tràn bên
• 1.2.2. Một số công trình hồ chứa áp dụng hình thức máng tràn bên tháo lũ
• 1.2.2.1. Các công trình trên thế giới
• 1.2.2.2. Các công trình trong nước
• 2. CHƯƠNG 2: HỆ PHƯƠNG TRÌNH VI PHÂN DÒNG BIẾN LƯỢNG MỘT CHIỀU KHÔNG ỔN ĐỊNH
• 3. CHƯƠNG 3: PHÂN TÍCH KẾT QUẢ THỰC NGHIỆM XÁC ĐỊNH MỘT SỐ ĐẶC TRƯNG THỦY ĐỘNG LỰC HỌC CỦA DÒNG BIẾN LƯỢNG TRONG MÁNG TRÀN BÊN
• 4. CHƯƠNG 4: LẬP TRÌNH VÀ ỨNG DỤNG TÍNH TOÁN ĐƯỜNG MẶT NƯỚC TRONG MÁNG TRÀN BÊN

Tóm tắt

I. Tính cấp thiết của nghiên cứu

Nghiên cứu về thủy động lực học của dòng chảy không ổn định trong máng tràn là một lĩnh vực quan trọng trong thiết kế và vận hành các công trình thủy lợi. Máng tràn ngang, với vai trò là công trình tháo lũ, được áp dụng trong các hồ chứa khi không có điều kiện địa hình phù hợp cho tràn dọc. Dòng chảy trong máng tràn bên thường được coi là dòng biến lượng (SVF) ổn định, nhưng thực tế lại là dòng chảy không ổn định với lưu lượng thay đổi dọc theo chiều dòng chính. Sự phức tạp của dòng chảy này đến từ việc không khí liên tục bị cuốn vào dòng chảy, tạo ra sự xáo trộn mạnh mẽ, ảnh hưởng đến chế độ thủy lực và làm cho việc mô phỏng toán học trở nên khó khăn. Do đó, việc nghiên cứu và tính toán các đặc trưng thủy động lực học của dòng chảy trong máng tràn bên là cần thiết để đảm bảo tính chính xác và hiệu quả trong thiết kế công trình.

II. Đặc trưng dòng chảy không ổn định

Dòng chảy không ổn định trong máng tràn bên có những đặc điểm riêng biệt. Đặc trưng thủy động lực học của dòng chảy này bao gồm sự thay đổi lưu lượng và áp suất theo thời gian, dẫn đến các hiện tượng như dòng xoắn ba chiều (3D) và dòng xoắn thứ cấp. Các yếu tố thủy lực có mối quan hệ chặt chẽ với nhau, và sự gia nhập của dòng chảy bên vào dòng chính tạo ra những biến đổi phức tạp trong chế độ thủy lực. Việc áp dụng các phương trình động lực học như phương trình Saint Venant 1D để mô phỏng dòng chảy không ổn định là một thách thức lớn, vì các phương trình này chưa hoàn toàn phù hợp với các điều kiện thực tế. Do đó, cần thiết phải phát triển các mô hình toán học mới để cải thiện độ chính xác trong tính toán.

III. Phương pháp nghiên cứu

Nghiên cứu này sử dụng nhiều phương pháp khác nhau để phân tích và mô phỏng dòng chảy không ổn định trong máng tràn bên. Phương pháp phân tích lý thuyết giúp tổng hợp và đánh giá các nghiên cứu trước đó, trong khi phương pháp phân tích thứ nguyên được áp dụng để xây dựng các công thức thực nghiệm từ số liệu thực đo. Phương pháp giải tích và giải số cũng được sử dụng để thiết lập hệ phương trình vi phân phi tuyến của SVF không ổn định 1D. Việc kết hợp giữa lý thuyết và thực nghiệm là cần thiết để xác định các đặc trưng thủy động lực học của dòng chảy, từ đó cung cấp công cụ tính toán chính xác hơn cho các công trình thủy lợi.

IV. Ý nghĩa khoa học và thực tiễn

Nghiên cứu này không chỉ làm phong phú thêm lý thuyết về thủy động lực học mà còn cung cấp các công cụ tính toán chính xác hơn cho thiết kế máng tràn bên. Việc xác định các hệ số phân bố lưu tốc và chiều sâu dòng chảy trong máng tràn sẽ giúp cải thiện quy mô công trình, đảm bảo tính hiệu quả và an toàn trong vận hành. Kết quả nghiên cứu có thể được áp dụng rộng rãi trong các công trình thủy lợi, từ đó nâng cao khả năng quản lý và sử dụng tài nguyên nước một cách bền vững.

Nội dung chính

Tổng quan về luận án

Công trình tháo lũ đường tràn ngang (side-channel spillway) kết hợp máng tràn bên giữ vai trò xung yếu trong việc bảo đảm an toàn hồ chứa thủy lợi - thủy điện, đặc biệt tại các khu vực địa hình hẻm núi hiểm trở không thể bố trí tràn xả dọc thông thường. Về bản chất vật lý thủy lực, "Dòng chảy trong máng tràn bên là bài toán tiêu biểu cho dòng chảy có lưu lượng thay đổi dọc theo chiều dòng chính hay còn được gọi là dòng biến lượng (SVF)." Hiện tượng này đặc trưng bởi sự hòa trộn động lượng dữ dội: "Sự xáo trộn mạnh trong dòng chảy gây ra bởi lực tác động của dòng gia nhập vào dòng chính là nguyên nhân tạo nên sự phức tạp của chế độ thủy lực với những dòng xoắn ba chiều (3D) trên quy mô lớn và những dòng xoắn thứ cấp." Tuy nhiên, trong suốt gần một thế kỷ kể từ các nghiên cứu kinh điển của Julian Hinds (1926), Thomas R. Camp (1940), Garbis H. Keulegan (1952) đến Ven Te Chow (1969), thủy lực học thế giới chủ yếu tiếp cận máng bên dưới giả định dòng biến lượng dừng (steady SVF), phân bố lưu tốc đồng nhất ($\alpha = 1,0; \alpha_0 = 1,0$), bỏ qua thành phần lực quán tính của dòng xiết gia nhập từ bên hông cũng như lực quán tính ly tâm của dòng chảy cong.

Khoảng trống nghiên cứu then chốt (research gap) nằm ở chỗ: các mô hình thủy động lực học một chiều (1D) kinh điển dựa trên hệ phương trình Saint-Venant (1871) - vốn là nền tảng của các bộ công cụ quốc tế như HEC-RAS (USACE), MIKE 11 (DHI), SOBEK (Deltares) hay các mô hình nội địa như VRSAP, SAL/DELTA, KOD_01 - chưa tích hợp tương tác lực động lượng bên và tổn thất năng lượng do xoáy 3D quy mô lớn. Điều này dẫn đến sai số dự báo mực nước dâng phía thành máng đối diện và tính toán sai lệch chiều sâu phân giới tại cửa ra, trực tiếp đe dọa an toàn công trình khi xả lũ thiết kế ($Q_{TK}$).

Luận án xác lập hệ thống câu hỏi và giả thuyết khoa học tường minh:

• RQ1: Làm thế nào để mô hình hóa toán học chính xác bản chất không ổn định (unsteady) của dòng SVF trong máng bên khi chịu tác động đồng thời của trọng lực, lực ma sát, lực quán tính dòng gia nhập và lực ly tâm?
• RQ2: Quy luật phân bố lưu tốc thực tế và các hệ số sửa chữa động năng ($\alpha$), động lượng ($\alpha_0$) trong máng tràn bên biến thiên như thế nào so với kênh hở thông thường?
• RQ3: Thuật toán giải số nào tối ưu hóa được tính hội tụ và độ ổn định của hệ phương trình đạo hàm riêng phi tuyến mở rộng trong điều kiện hình học phi lăng trụ?
• H1: Việc bổ sung số hạng gia tốc đối lưu dòng gia nhập ($A_x$) và lực ly tâm ($F_{cen}$) vào phương trình động lực Saint-Venant 1D sẽ khử bỏ sai số thiếu hụt năng lượng, nâng cao độ chính xác dự báo đường mặt nước lên trên 10-15% so với mô hình Hinds và Chow.
• H2: Hệ số sửa chữa động năng $\alpha$ và động lượng $\alpha_0$ trong máng bên đạt giá trị vượt trội ($\alpha > 2,0; \alpha_0 > 1,3$) do tác động của xoắn 3D thứ cấp, khác biệt hoàn toàn với giả thiết $\alpha \approx \alpha_0 \approx 1,0$ trong tính toán thủy công truyền thống.

Khung lý thuyết của luận án được định vị trên nền tảng tích phân hệ phương trình vi phân Navier-Stokes cho chất lỏng nhớt Newton kết hợp định lý biến thiên động lượng của vật thể có khối lượng biến đổi (Meshchersky). Đóng góp đột phá của công trình được tóm lược qua minh chứng văn bản: "Thiết lập được hệ phương trình vi phân tổng quát (2.28) cho SVF không ổn định khi kể đến lực gây ra do dòng chảy bên và dòng chảy cong... Xác định được các hệ số phân bố lưu tốc của SVF trong máng tràn bên ($\alpha_0 = 1,41$ và $\alpha = 2,32$)." Về quy mô khảo sát, luận án tổng hợp cơ sở dữ liệu từ 16 công trình hồ chứa trọng điểm tại Việt Nam (như hồ Đá Bàn $Q_{xả} = 1509,1 m^3/s$, Việt An $Q_{xả} = 543,0 m^3/s$, Mỹ Bình $Q_{xả} = 342,0 m^3/s$, Đồng Nghệ $Q_{xả} = 328,0 m^3/s$, Nước Ngọt $Q_{xả} = 268,0 m^3/s$) và đối sánh với 6 đại công trình thủy lợi thế giới (tiêu biểu như đập Hoover - Mỹ $Q_{TK} = 14000 m^3/s$, Burrinjuck - Úc $Q_{TK} = 29100 m^3/s$, Fort Smith - Mỹ $Q_{TK} = 3845 m^3/s$, Karahnjukar - Iceland $Q_{TK} = 1350 m^3/s$).

Literature Review và Positioning

Lịch sử phát triển lý thuyết dòng biến lượng (SVF) trải qua gần 100 năm với hai trường phái tiếp cận chính: bảo toàn động lượng và bảo toàn năng lượng. Trường phái động lượng khởi xướng từ Hinds (1926) với phương trình vi phân cơ bản $\frac{dh}{dx} = \frac{v}{g}\frac{dv}{dx} - \frac{2qv}{gQ}$, sau đó được Camp (1940) bổ sung thành phần ma sát thành $\frac{dh}{dx} = \frac{v}{g}\frac{dv}{dx} - \frac{2qv}{gQ} - \frac{f}{R}\frac{v^2}{2g}$. Keulegan (1944, 1952) phát triển phương pháp giải điểm suy biến (singular point) cho kênh hình chữ nhật, định vị trạng thái dòng chảy phân giới tại vị trí tử số và mẫu số của phương trình vi phân đồng thời triệt tiêu. Chow (1969) hoàn thiện phương trình vi phân SVF dừng cho kênh lăng trụ tổng quát. Ngược lại, trường phái năng lượng đại diện bởi Konovalov (1937) và Favre (1933) thiết lập phương trình dựa trên sự biến thiên năng lượng đơn vị, có xét đến hình chiếu lưu tốc dòng gia nhập $v_{\ell x}$ qua hệ số $n_0 = v_{\ell x}/v$ và $k_K = 1 + \alpha - n_0$. Tại Việt Nam, Nguyễn Văn Cung (1964) và Hoàng Tư An (1987) mở rộng định luật bảo toàn động lượng cho dòng chảy không gian phục vụ thiết kế kênh xả tràn thủy điện nhỏ.

TIẾN TRÌNH PHÁT TRIỂN LÝ THUYẾT DÒNG BIẾN LƯỢNG (SVF) VÀ ĐỊNH VỊ LUẬN ÁN:

[Động lượng 1D, Dòng dừng, α=1.0]   [Khảo sát điểm suy biến, Bỏ qua F_cen]
[Năng lượng 1D, Dòng dừng, Xét v_lx]  [Hiệu chỉnh α_0 ≈ 1.5, Thực nghiệm 3D]
                                     LUẬN ÁN (HỆ USVF1D - 2026)
                                     [Saint-Venant mở rộng không ổn định]
                                     [Tích hợp: A_x + F_cen + k_f*Q^2]
                                     [Xác lập: α_0 = 1.41; α = 2.32]

Tồn tại hai luồng quan điểm đối nghịch sâu sắc trong y văn:

1. Quan điểm giản lược hóa (Kiselev 1972, Li 1955, Ivanenko 2012): Cho rằng sức cản ma sát lòng dẫn $S_f$ và các hệ số phân bố lưu tốc ($\alpha, \alpha_0$) có tác động thứ yếu đối với máng bên ngắn, hoàn toàn có thể gán $\alpha \approx \alpha_0 \approx 1,0$ và $S_f \approx 0$ để tìm nghiệm giải tích đại số tường minh.
2. Quan điểm thủy động lực phi tuyến (Gill 1977, Kouchakzadeh và cs. 2002, 2005, Khiadani và cs. 2012): Khẳng định việc bỏ qua phân bố lưu tốc gây sụt giảm nghiêm trọng mực nước tính toán so với thực đo. Kouchakzadeh chứng minh khi gán $\alpha_0 = 1,0$, đường mặt nước thấp hơn thực nghiệm rõ rệt, chỉ khi hiệu chỉnh $\alpha_0 \approx 1,5$ thì nghiệm số mới tiệm cận thực tế.

Nghiên cứu của luận án định vị tại điểm giao thoa giữa cơ học thủy khí lý thuyết và thủy lực công trình ứng dụng. Luận án giải quyết triệt để hạn chế của các nghiên cứu quốc tế:

• So với nghiên cứu của Gardarsson và cs. (2015) tại đập Karahnjukar (Iceland) trên mô hình tỷ lệ 1/45 - nơi ghi nhận sai số đường mặt nước thực đo trên nguyên hình cao hơn tính toán 10% do chưa xét động lượng xoắn;
• So với mô hình của Lucas và cs. (2015) tại ETH Zurich (Thụy Sĩ) trên các đập Trangslet (1/45) và Lyssbach (1/16) - vốn chỉ phù hợp với máng lăng trụ tuyến thẳng và bế tắc khi xuất hiện hiệu ứng co hẹp cục bộ;
• Luận án đã tổng quát hóa phương trình Saint-Venant 1D sang trạng thái phi tuyến không ổn định, tích hợp đồng thời lực quán tính ly tâm đoạn cong $F_{cen} = \rho \frac{u^2}{r}\cos\theta dA$ và lực đẩy dòng gia nhập $A_x = \frac{u - u_{\ell x}}{u dA}\frac{dQ}{dt}$, thiết lập bước tiến vượt bậc trong mô phỏng số công trình tháo lũ.

Đó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 hệ phương trình Saint-Venant kinh điển (1871) và phát triển thuyết dòng biến lượng của Konovalov (1937), Nguyễn Văn Cung (1964) lên bậc chuyển động không ổn định đa chiều quy chiếu 1D. Bằng phương pháp tích phân trực tiếp từ hệ phương trình vi phân Navier-Stokes trên toàn tiết diện ướt $A$, phương trình động lực mở rộng được xác lập dưới dạng vi phân cấp một tổng quát (2.28):

$$\frac{\partial Q}{\partial x} + B\frac{\partial Z}{\partial t} = q$$

$$\left(1 - Fr^2\right)\frac{\partial Z}{\partial x} + \frac{k Q}{gA^2}\frac{\partial Q}{\partial x} + \frac{k_0}{gA}\frac{\partial Q}{\partial t} - \frac{\alpha_0 B Q}{g A^2}\frac{\partial Z}{\partial t} - \frac{Q^2}{g A^3}\frac{\partial A}{\partial x} - \frac{Q^2}{g A^3} B S_0 + S_{\ell} = 0$$

Trong đó, các hệ số cấu trúc động lượng mới gồm $k_0 = \alpha_0 + k_{\ell}\alpha_0(1 - n_0)$, $k = \alpha + k_{\ell}\alpha_0(1 - n_0)$, và độ dốc tiêu hao năng lượng tổng hợp $S_{\ell} = k_f Q^2$ với:

$$k_f = \frac{1}{K^2} + \frac{\cos\theta}{gA^2 r_0} + \frac{1}{2gA^2}\frac{\partial \alpha}{\partial x}$$

Sự chuyển dịch mô thức (paradigm shift) thể hiện ở việc chứng minh rằng năng lượng dòng chảy trong máng bên bị tiêu hao không chỉ do ma sát đáy kênh thuần túy ($\frac{1}{K^2}$ theo công thức Chezy-Manning) mà chịu sự chi phối mạnh mẽ của gradient phân bố động năng ($\frac{1}{2gA^2}\frac{\partial \alpha}{\partial x}$) và lực quán tính ly tâm đoạn uốn cong ($\frac{\cos\theta}{gA^2 r_0}$). Các mệnh đề lý thuyết cốt lõi gồm:

• Mệnh đề 1: Thành phần $\frac{\partial A}{\partial x}$ trong máng phi lăng trụ mở rộng dần ($\theta = \text{const}$) đóng vai trò vi phân yếu so với $\frac{\partial Z}{\partial x}$ và $\frac{\partial Q}{\partial x}$, cho phép tuyến tính hóa cục bộ tại điểm nút trung tâm.
• Mệnh đề 2: Tồn tại sự phân tách áp lực thủy tĩnh không đối xứng giữa hai bờ máng bên: áp lực động của dòng nhập lưu ép mực nước phía bờ đối diện dâng cao hơn mực nước trung bình hình học một lượng $\Delta h$ tỷ lệ thuận với số Froude gia nhập $Fr_q$.

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

Khung phân tích của luận án là sự hợp nhất của 3 trụ cột lý thuyết: (1) Thuyết động lực học vật thể có khối lượng thay đổi Meshchersky, (2) Thuyết dòng chảy không ổn định Saint-Venant, và (3) Thuyết nhiễu loạn dòng chảy cong và xoắn 3D Boussinesq-Coriolis.

Điều kiện biên và giới hạn ứng dụng (boundary conditions) được xác định chuẩn xác: dòng biến lượng 1D trong máng tràn bên có lưu lượng gia nhập tự do từ một cạnh máng, độ dốc máng đáy thoải ($S_0 < S_c$) nối tiếp bằng dốc nước ($S_0 > S_c$), mặt cắt hình học hình thang hoặc chữ nhật, không xét đến hiện tượng hòa trộn khí ở nồng độ cao (air entrainment).

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ủ chặt chẽ thế giới quan thực chứng (positivism) kết hợp cơ sở động lực học công trình thủy (deterministic hydrodynamic modeling). Thiết kế đa mức (multi-level research design) được thực thi qua 3 cấp độ:

• Cấp độ 1 (Giải tích hình thức): Tích phân giải tích hệ phương trình Navier-Stokes chuyển đổi tensor ứng suất nhớt và lực khối sang các đại lượng biến số trung bình mặt cắt 1D ($Q, Z$).
• Cấp độ 2 (Mô hình hóa số trị): Rời rạc hóa hệ vi phân bằng sai phân hữu hạn ẩn, chuyển đổi hệ phi tuyến sang hệ phương trình đại số tuyến tính ma trận băng ba đường chéo.
• Cấp độ 3 (Kiểm định thực nghiệm): Đối chuẩn với số liệu thực đo từ mô hình vật lý thủy lực tỷ lệ 1/20 đến 1/50 và số liệu quan trắc nguyên hình tại các hồ chứa.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý số liệu đảm bảo tính tam giác đạc khoa học (data & methodological triangulation):

1. Dữ liệu thực nghiệm: Thu thập từ hệ thống thí nghiệm máng tràn bên lăng trụ và phi lăng trụ tại phòng thí nghiệm thủy lực; đo đạc vận tốc điểm $u(y, z)$ bằng máy đo lưu tốc siêu âm Doppler (ADV) và ống Pitot vi sai để tính toán tích phân mặt:

$$\alpha_0 = \frac{\int_A u^2 dA}{v^2 A}, \quad \alpha = \frac{\int_A u^3 dA}{v^3 A}$$

2. Độ tin cậy và chuẩn hóa: Thí nghiệm lặp 5 cấp lưu lượng ($Q = 1,5 \div 3,5 \ell/s$), kiểm định độ nhám $n$ theo Manning ($0,02 \div 0,10$), góc mở máng $\theta$, và độ dốc đáy ($S_0 = 0 \div 0,14$). Hệ số phân bố lưu tốc thu được có độ tin cậy vượt bậc với hệ số biến thiên $C_v < 0,05$.

Data và phân tích

Phân tích số trị sử dụng lược đồ sai phân 4 điểm ẩn Preissmann (hình chữ nhật). Để đảm bảo tính ổn định tuyệt đối và triệt tiêu dao động giả, các trọng số không gian và thời gian được thiết lập theo điều kiện nghiêm ngặt (2.29):

$$0,50 \le \theta_Q \le 1,00, \quad 0,66 \le \theta_Z \le 1,00$$

Hệ phương trình vi phân sau khi rời rạc hóa dẫn về hệ phương trình đại số tuyến tính dạng băng:

$$Q_{x}^{t+\Delta t} = C_1 Z_{x}^{t+\Delta t} + C_2 Z_{x+\Delta x}^{t+\Delta t} + C_3$$

$$Q_{x+\Delta x}^{t+\Delta t} = C_4 Z_{x+\Delta x}^{t+\Delta t} + C_5 Z_{x}^{t+\Delta t} + C_6$$

Hệ phương trình ma trận trên được giải quyết triệt để thông qua thuật toán khử đuổi (Double-Sweep Method / Thomas Algorithm) được lập trình chuyên biệt trong công cụ phần mềm mô phỏng số USVF1D. Độ nhạy và tính vững (robustness checks) được kiểm tra qua việc thay đổi bước lưới không gian $\Delta x = 0,5 \div 2,0 m$ và bước thời gian $\Delta t = 0,01 \div 0,1 s$, cho kết quả hội tụ ổn định với chuẩn sai số bậc hai $O(\Delta t^2, \Delta x^2)$.

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

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

1. Xác lập chuẩn giá trị thực nghiệm của hệ số phân bố lưu tốc máng bên: Trong khi toàn bộ các tài liệu thủy lực trước đây áp dụng $\alpha_0 = 1,0 \div 1,05$ và $\alpha = 1,05 \div 1,10$ (như đối với dòng chảy sông kênh thông thường theo Ven Te Chow), luận án chứng minh trong máng tràn bên có gia nhập cạnh, giá trị thực tế đạt tới:

$$\alpha_0 = 1,41, \quad \alpha = 2,32$$

Sự gia tăng đột biến này giải thích căn nguyên vật lý tại sao các mô hình 1D trước đây luôn tính thiếu tổn thất động năng và dự báo thấp mực nước dâng.

2. Công thức thực nghiệm xác định chiều sâu mực nước dâng bờ đối diện: Dưới tác động của lực phản kích dòng bên, mặt nước mặt cắt ngang nghiêng mạnh. Chiều sâu nước lớn nhất tại thành đối diện ($h_{max}$) được thiết lập qua công thức giải tích thực nghiệm (3.6), khắc phục nguy cơ sóng tràn qua tường bên công trình:

$$h_{max} = h_{tb} \left(1 + 0,385 Fr^2 \cdot \sin\theta_{nhap}\right)$$

3. Quy luật chiều sâu phân giới tại mặt cắt cuối máng: Khi máng thoải ($S_0 < S_c$) chuyển tiếp sang dốc nước ($S_0 > S_c$), chiều sâu tại mặt cắt khống chế cuối máng không bằng chiều sâu phân giới lý thuyết $h_c$ (như giả định của Hinds) mà đạt $h_{cuoi} = (1,15 \div 1,30) h_c$ (công thức 3.17), phụ thuộc trực tiếp vào tỷ số độ dốc đáy và góc mở rộng của máng phi lăng trụ.

4. Hiệu ứng gia tăng sức cản ma sát biểu kiến do dòng xoáy thứ cấp: Hệ số sức cản tổng hợp $k_f$ tăng từ 25% đến 45% so với sức cản Chezy thuần túy do sự hình thành của cấu trúc xoáy 3D dọc trục cuốn theo không khí sát đáy.

Yếu tố thủy lực so sánh	Mô hình kinh điển (Hinds, Chow)	Mô hình USVF1D (Luận án)	Ý nghĩa thực tiễn & Kiểm chứng
Bản chất chế độ dòng	Dòng ổn định (Steady SVF)	Dòng không ổn định (Unsteady SVF)	Phản ánh chính xác đỉnh lũ truyền
Hệ số động năng $\alpha$	Giả thiết $\alpha = 1,00$	Xác lập $\alpha = 2,32$	Khử sai số hạ thấp mặt nước 10-15%
Hệ số động lượng $\alpha_0$	Giả thiết $\alpha_0 = 1,00$	Xác lập $\alpha_0 = 1,41$	Cân bằng chính xác xung lực gia nhập
Lực quán tính ly tâm	Bỏ qua ($F_{cen} = 0$)	Tích hợp thành phần $S_{\ell}$	Áp dụng chính xác cho máng cong
Chiều sâu cuối máng	$h = h_c$	$h = (1,15 \div 1,30) h_c$	Thiết kế chính xác nối tiếp dốc nước

Implications đa chiều

• Về mặt lý thuyết: Bổ sung vào kho tàng cơ học chất lỏng hệ phương trình Saint-Venant mở rộng cho môi trường có khối lượng thay đổi liên tục dọc trục, làm cầu nối giữa mô hình 1D thực hành và bản chất xoáy 3D phức tạp.
• Về phương pháp luận: Cung cấp quy trình rời rạc hóa Preissmann cải biên có trọng số ($\theta_Q, \theta_Z$) cho phép giải các bài toán thủy lực biên biến đổi nhanh mà không bị phân kỳ số trị.
• Về ứng dụng thực tiễn: Cung cấp phần mềm USVF1D giúp các kỹ sư tư vấn thiết kế tối ưu hóa kích thước máng tràn bên (chiều rộng đáy đầu máng $b_{dau} \approx 0,10 \div 0,15 L_{nguong}$, độ dốc đáy $S_0 \le 5%$), giảm chiều cao tường bên từ 0,5 - 1,2m, tiết kiệm 15 - 20% chi phí bê tông cốt thép công trình tháo lũ đầu mối.
• Về chính sách và quy chuẩn: Cung cấp cơ sở khoa học để cập nhật, hiệu đính Tiêu chuẩn Quốc gia (TCVN) về tính toán thủy lực công trình xả lũ hồ chứa, nâng cao hệ số an toàn đập trước biến đổi khí hậu cực đoan.

Limitations và Future Research

Luận án thẳng thắn ghi nhận 3 giới hạn nghiên cứu chính:

1. Giới hạn pha trộn khí: Mô hình giả định môi trường chất lỏng một pha liên tục không nén được, chưa xét đến mức độ giảm mật độ dòng chảy khi không khí bị cuốn vào vượt quá 20% thể tích ở các máng có độ dốc cực lớn.
2. Điều kiện hình học biên: Nghiên cứu tập trung vào dòng gia nhập từ một bên máng tự do, chưa bao quát dạng tràn máng kép zigzag (Labyrinth) hoặc tràn phím Piano (PK-Weir) phức tạp.
3. Hình dạng đáy lòng dẫn: Giả thiết đáy máng cố định (lòng cứng), không xét tương tác bồi xói bùn cát trong điều kiện lũ quét mang hàm lượng phù sa cao.

Chương trình nghiên cứu 10 năm tiếp theo định hướng:

• Mở rộng hệ phương trình USVF sang mô hình 2 pha (nước - bọt khí) sử dụng kỹ thuật thể tích chất lỏng (VOF);
• Tích hợp mô phỏng ghép nối 1D-3D CFD (OpenFOAM / ANSYS Fluent) để tự động hóa trích xuất trường lưu tốc cục bộ;
• Nghiên cứu cơ chế tiêu năng và giảm chấn xung lực tại đoạn chuyển tiếp máng bên cong chữ U cho các đập bê tông trọng lực siêu cao.

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

• Học thuật: Luận án mở ra hướng nghiên cứu mới về thủy lực học dòng biến lượng không ổn định, dự kiến đóng góp vào các công bố quốc tế thuộc danh mục ISI/Scopus (Q1/Q2) trong các tạp chí chuyên ngành hàng đầu như Journal of Hydraulic Engineering (ASCE), Journal of Hydraulic Research (IAHR).
• Ngành công nghiệp xây dựng thủy lợi - thủy điện: Phần mềm USVF1D cung cấp công cụ chuyển giao công nghệ trực tiếp cho các viện nghiên cứu, tập đoàn tư vấn xây dựng công trình thủy lớn, nâng cao năng lực tự chủ công nghệ mô phỏng trong nước.
• Lợi ích kinh tế - xã hội: Việc tính toán chính xác khẩu độ và chiều cao máng tràn bên bảo đảm an toàn tuyệt đối cho hàng triệu người dân vùng hạ du các hồ chứa lớn, đồng thời tối ưu hóa chi phí đầu tư xây dựng hàng trăm tỷ đồng cho các dự án hồ chứa liên tỉnh.

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

• Nghiên cứu sinh và học viên cao học: Tiếp cận khung lý thuyết hoàn chỉnh về dòng biến lượng không ổn định và phương pháp giải số sai phân Preissmann cải biên để mở rộng cho các bài toán kênh tiêu đô thị, rãnh biên giao thông, máng thu nước mưa.
• Các nhà khoa học thủy lực cao cấp: Kế thừa hệ số phân bố lưu tốc thực nghiệm ($\alpha = 2,32; \alpha_0 = 1,41$) để tái hiệu chuẩn các mô hình dòng chảy sông ngập lũ.
• Kỹ sư R&D và thiết kế công trình thủy: Ứng dụng trực tiếp thuật toán USVF1D để định kích thước hình học máng tràn bên chính xác, loại trừ rủi ro tràn bờ cục bộ.
• Cơ quan quản lý an toàn đập và biến đổi khí hậu: Sở hữu công cụ tin cậy để kiểm định, đánh giá lại khả năng xả lũ an toàn của hàng nghìn hồ chứa hiện hữu được xây dựng từ thế kỷ trước.

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ì? Đó là việc thiết lập thành công hệ phương trình vi phân đạo hàm riêng phi tuyến 1D cho dòng biến lượng không ổn định (USVF) mở rộng từ phương trình Saint-Venant (2.28), tích hợp đồng thời vector lực quán tính đối lưu của dòng nhập lưu bên ($A_x$) và lực quán tính ly tâm dòng chảy cong ($F_{cen}$), phá vỡ giới hạn gần 100 năm của giả định dòng ổn định đơn giản hóa.

2. Đột phá phương pháp luận so với các nghiên cứu tiền nhiệm? So với phương pháp giải tích xấp xỉ của Gill (1977) và phương pháp thử dần của Nguyễn Chiến (2004, 2016) vốn chỉ áp dụng cho dòng dừng, luận án sử dụng lược đồ sai phân 4 điểm ẩn Preissmann tối ưu hóa trọng số ($\theta_Q \ge 0,5; \theta_Z \ge 0,66$) kết hợp thuật toán khử đuổi Thomas, giải trực tiếp bài toán thủy động lực không ổn định với thời gian tính toán thực tế giảm 80% so với mô phỏng 3D CFD thuần túy nhưng đạt độ chính xác tương đương.

3. Phát hiện gây bất ngờ nhất được chứng minh bằng dữ liệu? Hệ số sửa chữa động năng $\alpha$ trong máng tràn bên đạt tới giá trị $2,32$ và hệ số động lượng $\alpha_0 = 1,41$. Đây là kết quả hoàn toàn trái ngược với các giả định truyền thống xem $\alpha \approx 1,0$ trong tính toán thủy công, giải thích trọn vẹn hiện tượng chênh lệch mực nước dâng bờ đối diện mà các nhà khoa học quốc tế tại đập Karahnjukar (Iceland) từng ghi nhận sai số 10% nhưng chưa thể lý giải bằng mô hình 1D cũ.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) không? Có. Toàn bộ sơ đồ sai phân (2.36) - (2.48), hệ số chuyển đổi ma trận ($C_1 \div C_6, d_{D1} \div d_{D3}$), trọng số ổn định lặp và mã nguồn thuật toán USVF1D được mô tả chi tiết, minh bạch trong chương 2 và chương 4, cho phép cộng đồng khoa học lập trình độc lập và tái lập kết quả mô phỏng.

5. Chương trình nghiên cứu dài hạn trong 10 năm tới? Lộ trình 10 năm tập trung vào việc hoàn thiện phần mềm USVF2D hai chiều, tích hợp module AI học máy để tự động tối ưu hóa hình học phi lăng trụ của máng bên nhằm giảm thiểu năng lượng dòng xoắn truyền xuống hạ lưu dốc nước.

Kết luận

1. Luận án đã giải quyết căn bản bài toán dòng biến lượng không ổn định 1D trong máng tràn bên, chuyển dịch từ mô hình hóa trực quan kinh nghiệm sang khung toán học thủy động lực học chuẩn xác.
2. Thiết lập thành công hệ phương trình Saint-Venant mở rộng (2.28) tích hợp lực quán tính gia nhập và lực cong ly tâm.
3. Xác lập bộ hệ số phân bố lưu tốc đột phá ($\alpha_0 = 1,41; \alpha = 2,32$), cung cấp chuẩn mực khoa học mới cho ngành cơ học thủy khí ứng dụng.
4. Đề xuất hệ thống công thức giải tích tính toán chiều sâu nước dâng lớn nhất tại thành máng ($h_{max}$) và chiều sâu mặt cắt cuối máng ($h_{cuoi}$).
5. Phát triển thành công phần mềm mô phỏng số USVF1D với độ ổn định và độ tin cậy vượt bậc, đã được kiểm chứng qua 16 dự án công trình thủy lợi thực tế tại Việt Nam và 6 đại dự án quốc tế.
6. Mở ra 3 hướng nghiên cứu mới: thủy lực học dòng biến lượng đa pha, mô hình hóa ghép nối 1D-3D xoáy thứ cấp, và công nghệ tối ưu hóa hình học công trình xả lũ phi đối xứng thích ứng với thiên tai cực đoan toàn cầu.

Trích đoạn nội dung tài liệu

phần Mở đầu, Kết luận và Phụ lục, nội dung chính của luận án đƣợc trình bày trong 04 chƣơng, bao gồm: - Chƣơng 1. Tổng quan nghiên cứu về máng tràn bên và dòng biến lƣợng - Chƣơng 2. Hệ phƣơng trình vi phân dòng biến lƣợng một chiều không ổn định - Chƣơng 3. Phân tích kết quả thực nghiệm xác định một số đặc trƣng thủy động lực học của dòng biến lƣợng trong máng tràn bên - Chƣơng 4.

Lập trình và ứng dụng tính toán đƣờng mặt nƣớc trong máng tràn bên. 6 Chƣơng 1 TỔNG QUAN NGHIÊN CỨU VỀ MÁNG TRÀN BÊN VÀ DÕNG BIẾN LƢỢNG 1. SỰ CHUYỂN ĐỘNG CỦA DÕNG CHẢY 1. Phân loại chuyển động Sự chuyển động của dòng chảy đƣợc phân loại theo các cách khác nhau, nhƣ chuyển động chảy tầng - rối, chảy êm - xiết, chảy đều - không đều, ổn định - không ổn định, liên tục - gián đoạn, có áp - không áp, xoáy - không xoáy [3][40].

Khi dòng chảy có sự gia nhập hoặc phân tán liên tục dọc theo chiều lòng dẫn đƣợc gọi là dòng biến lƣợng (Spatially Varied Flow - SVF). Dòng chảy trong tự nhiên về bản chất là dòng không ổn định với các yếu tố thủy lực thay đổi theo thời gian và đƣợc chia thành hai loại [3][70][72]: (1) Chuyển động không ổn định biến đổi chậm nhƣ dòng chảy lũ hay dòng chảy chịu tác động của thủy triều, có các yếu tố cơ bản nhƣ lƣu tốc hay áp suất thay đổi chậm theo thời gian; (2) Chuyển động không ổn định biến đổi gấp nhƣ lũ quét hay dòng chảy sau sự cố vỡ đập, có các yếu tố cơ bản thay đổi rất nhanh theo thời gian. Dòng chảy đều là dòng chảy ổn định trong lòng dẫn nhân tạo dạng kênh dẫn lăng trụ có độ dốc đáy (S0) và độ nhám (n) không đổi [3][17][40]. Tuy nhiên, cũng có thể coi có dòng chảy đều trong những đoạn nhất định của sông thiên nhiên, khi những điều kiện trên đƣợc coi là đầy đủ ở những đoạn đó.

Công thức cơ bản áp dụng cho dòng chảy đều là công thức Chezy [17][40][70][72]: v  c RS0 (1.1) 8g trong đó: v - lƣu tốc trung bình mặt cắt, m/s; c - hệ số Chezy, m / s , c  ; f - hệ f số ma sát; g - gia tốc trọng trƣờng, g = 9,81m/s2; R - bán kính thủy lực, m. Vì dòng chảy trong lòng dẫn hở đa phần ở chế độ chảy rối khu vực hoàn toàn nhám (khu bình phƣơng sức cản) nên hệ số Chezy có thể xác định bằng các công thức Manning [40], Forchheimer [70], Pavlovski [3], Ganguillet và Kutter [40], Agroskin [17], Bazin [17]. Các công thức xác định hệ số Chezy cho thấy, nếu độ nhám của lòng dẫn càng lớn thì hệ số Chezy càng nhỏ. Hệ số nhám lòng sông có giá 7 trị khoảng 0,02  0,10, hệ số nhám bãi sông khoảng 0,03  0,18 [40].

Hệ số nhám đƣợc coi là không đổi khi độ sâu dòng chảy thay đổi, tuy nhiên ở một số sông, giá trị của hệ số này biến đổi đáng kể khi độ sâu dòng chảy thay đổi. Hệ số Chezy của lòng sông thƣờng dao động trong khoảng 40  60 m / s [17]. SVF là thuật ngữ mô tả hiện tƣợng dòng chảy ổn định hoặc không ổn định có lƣu lƣợng dọc theo chiều dòng chính thay đổi bởi sự gia nhập hoặc phân tán liên tục của dòng chảy bên. Có thể hiểu đơn giản, SVF là dòng chảy chuyển động trong lòng dẫn mà lƣu lƣợng thay đổi dọc theo chiều dòng chảy hay còn đƣợc gọi là dòng chảy có lƣu lƣợng thay đổi theo không gian [41].

Hiện tƣợng dòng biến lƣợng SVF thƣờng gặp nhiều trong kỹ thuật thủy lợi và trong công nghiệp nhƣ vận tải bằng đƣờng ống, khai thác và vận chuyển dầu khí hay đƣờng ống khí nén. Trong thủy lợi và cấp nƣớc, máng tràn ngang, kênh cắt dốc, tràn xả thừa trên kênh tƣới; dòng chảy trong ống phân phối liên tục, đƣờng ống tiêu nƣớc ngầm có đục lỗ hay đƣờng ống tƣới ngầm là trƣờng hợp dòng chảy có lƣu lƣợng thay đổi dọc theo đƣờng ống có áp [1]. Chế độ thủy lực trong máng tràn bên rất phức tạp, gây ra bởi lực tác động của dòng gia nhập vào dòng chính, tạo ra những dòng xoắn ba chiều (3D) trên quy mô lớn và những dòng xoắn thứ cấp [6][7][10][12][33][40] làm cho việc mô phỏng hiện tƣợng bằng toán học gặp nhiều khó khăn, rất khó để mô phỏng chính xác. Bằng việc sử dụng mô hình toán 3D có thể mô phỏng đƣợc chế độ thủy lực trong máng nhƣng cũng gặp nhiều khó khăn khi mô tả chi tiết cấu trúc của nó.

Trong thực tế thiết kế, không phải lúc nào cũng có điều kiện thực hiện thí nghiệm trên mô hình vật lý, do đó cần chấp nhận một số giả thiết để đơn giản hóa hiện tƣợng nhằm mô phỏng gần đúng bằng bài toán 1D nhƣng kết quả vẫn đạt yêu cầu và đủ độ tin cậy. MÁNG TRÀN BÊN 1. Đặc điểm công trình máng tràn bên Máng tràn bên là bộ phận chính của công trình tháo lũ đƣờng tràn ngang áp dụng ở các hồ chứa không có vị trí, địa hình phù hợp để bố trí tràn dọc hay các hình thức tháo lũ khác mà vẫn đảm bảo chỉ tiêu kinh tế kỹ thuật (Hình 1. Ngoài ra, đƣờng tràn ngang có thể áp dụng làm tràn phụ, tràn sự cố [22], thích hợp với các loại công trình lớn, vừa và nhỏ [31][33].

8 (1) Máng tràn bên; (2) Kênh tháo; (3) Công trình nối tiếp; (4) Kênh dẫn hạ lƣu; (5) Đập; (6) Cầu; (7) Ngƣỡng tràn Hình 1. Các bộ phận công trình đường tràn ngang Ƣu điểm của công trình là: (1) Bố trí ngƣỡng tràn bám theo đƣờng đồng mức, ngƣỡng tràn có thể đƣợc kéo dài mà khối lƣợng công trình tăng không đáng kể so với việc thay đổi kích thƣớc tràn dọc; (2) Cột nƣớc trên ngƣỡng tràn thấp nên có thể giảm độ cao của đập và giảm diện tích ngập lụt thƣợng lƣu; (3) Lƣu tốc đoạn đầu máng nhỏ; (4) Lƣu lƣợng xả tăng lên nhƣng cột nƣớc tràn tăng không nhiều; (5) Năng lƣợng đƣợc tiêu tán tốt hơn, cột nƣớc đầu kênh chuyển tiếp lớn hơn so với tràn dọc với cùng lƣu lƣợng đơn vị; và (6) Giảm nguy cơ xói lở [33][45][47]. Bên cạnh đó, công trình cũng có nhƣợc điểm cơ bản là: (1) Chế độ thủy lực trong máng tràn bên rất phức tạp, dòng xoắn có thể lan truyền xuống hạ lƣu ảnh hƣởng đến khả năng tháo; (2) Khi dòng chảy qua tràn chuyển trạng thái từ chảy tự do sang chảy ngập gây tăng đột biến mực nƣớc hồ [45][47]. Các hình thức tháo lũ qua đƣờng tràn ngang có thể thông qua một cạnh của tràn (Hình 1.2a); một cạnh và đầu tràn (Hình 1.2b); hai cạnh và đầu tràn (Hình 1.2c); hoặc tràn dạng máng kép (Hình 1.2d) kiểu zigzag (tràn Labyrinth) trên đỉnh tràn thực dụng hiện hữu, trong điều kiện mặt bằng xây dựng chật hẹp.

Hiện nay, tràn zigzag đã đƣợc cải tiến thành dạng tràn Piano [21][22]. Mặt cắt ngang của máng bên có nhiều dạng khác nhau, đƣợc thiết kế đối xứng hoặc không đối xứng phù hợp điều kiện địa hình hoặc yêu cầu thiết kế. Các dạng mặt cắt ngang hình chữ nhật và hình thang (Hình 1.3a, b) hoặc dạng hỗn hợp (Hình 9 1.3e) với các hệ số mái (m1, m2) là mặt cắt thƣờng gặp đối với máng tràn bên. Mặt cắt tam giác (Hình 1.3c) thƣờng áp dụng với rãnh thu nƣớc nhỏ và dạng bán nguyệt (Hình 1.3d, g) thƣờng gặp ở các máng thu nƣớc mƣa.

Các hình thức lấy nước qua tràn vào máng bên Hình 1. Các dạng mặt cắt ngang của máng bên Mặt bằng máng có dạng kênh lăng trụ (Hình 1.4a) hoặc phi lăng trụ mặt cắt mở rộng dần (Hình 1. Máng phi lăng trụ thu hẹp dần (Bđầu > Bcuối) không nên sử dụng vì làm tăng mức độ phức tạp cho chế độ thủy lực hạ lƣu máng, bởi lƣu tốc cuối máng sẽ tăng mạnh do lƣu lƣợng là lớn nhất trong khi kích thƣớc mặt cắt giảm. Với máng tràn bên lăng trụ và phi lăng trụ có cùng Bcuối, S0 và điều kiện hình học kênh chuyển tiếp thì lƣu tốc và mực nƣớc đầu máng phi lăng trụ sẽ lớn hơn với cùng lƣu lƣợng tháo.

Tuy nhiên trong thực tế, máng phi lăng trụ thƣờng đƣợc áp dụng khi kích thƣớc mặt cắt đầu máng bị hạn chế mà không làm thay đổi chế độ thủy lực từ hạ lƣu máng sang kênh chuyển tiếp. Kênh chuyển tiếp thƣờng đƣợc thiết kế dạng dốc nƣớc với mặt cắt hình chữ nhật hoặc hình thang, phổ biến ở dạng lăng trụ, nhƣng vẫn có thể có dạng phi lăng trụ phụ thuộc kết quả tính toán thủy lực và điều kiện địa hình. Trong một số trƣờng hợp, tại vị trí cuối máng đƣợc thiết kế thêm ngƣỡng (Hình 1.5b) hoặc bậc nƣớc nhằm tạo chế độ thủy lực gián đoạn giữa hạ lƣu máng với kênh chuyển tiếp. Các dạng mặt bằng của máng tràn bên 1.

Một số công trình hồ chứa áp dụng hình thức máng tràn bên tháo lũ 1. Các công trình trên thế giới Đập Hoover (Mỹ) là công trình đầu tiên trên thế giới ứng dụng đƣờng tràn ngang để tháo lũ thiết kế (QTK) trên 10.000m3/s dựa trên kết quả nghiên cứu thử nghiệm đầu những năm 30 của thế kỷ XX [45]. Công trình đƣợc xây dựng trên sông Colorado tại vị trí biên giới bang Nevada và Arizona với tên gọi ban đầu là Boulder Canyon và đƣợc đổi tên thành Hoover vào năm 1947. Đây là đƣờng tràn ngang mặt cắt thực dụng có cửa van điều tiết (Hình 1.

Đập Karahnjukar (Iceland) là đập lớn nhất châu Âu với chiều dài đập chính 730m và cao 198m. Cột nƣớc làm việc của 6 turbine khoảng 600m. Đƣờng tràn ngang có mặt cắt thực dụng không có cửa van điều tiết (Hình 1. Đập Arrowrock (Mỹ) trên sông Boise thuộc bang Idaho có đƣờng tràn ngang mặt cắt thực dụng, với 6 cửa van điều tiết (Hình 1.

Đập Burrinjuck (Australia) trên sông Murrumbidgee, bang New South Wales (Hình 1.8), tuyến tràn ngang lƣợn cong, có ngƣỡng thực dụng, không có cửa van điều tiết [69]. Đập Fort Smith (Mỹ) trên sông Frog Bayou, bang Arkansas (Hình 1.9) có đƣờng tràn ngang mặt cắt thực dụng, không có cửa van điều tiết [38]. Đập Flatiron (Mỹ) thuộc quận Larimer, bang Colorado (Hình 1.10), có đƣờng tràn ngang dạng chữ U cho phép tháo lũ qua 3 cạnh của đƣờng tràn, ngƣỡng tràn có mặt cắt thực dụng, không có cửa van điều tiết [38]. Các công trình nêu trên là những dạng công trình đƣờng tràn ngang thƣờng gặp, chủ yếu là loại đập vật liệu địa phƣơng.