Luận án đặ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

Dòng chảy biến lượng không ổn định xuất hiện phổ biến tại các công trình tháo lũ đường tràn ngang. Lưu lượng nước biến thiên liên tục dọc theo chiều dòng chính. Chế độ thủy lực không duy trì trạng thái dừng. Mực nước và vận tốc tại mỗi mặt cắt thay đổi theo từng thời điểm. Dòng nước gia nhập mang theo động lượng riêng và tác động trực tiếp vào dòng chảy chính. Khối lượng chất lỏng trong kênh tăng dần tạo nên gradien áp suất thủy tĩnh biến đổi. Mối quan hệ giữa mực nước và lưu lượng không còn mang tính đơn trị. Việc coi dòng chảy là trạng thái ổn định sẽ dẫn đến nhiều sai số kỹ thuật. Nhận thức đúng bản chất không ổn định là tiền đề để nâng cao chất lượng thiết kế công trình thủy lợi.

Sự gia nhập của dòng nước tràn từ cạnh bên tạo ra lực va đập rất lớn. Không khí bị cuốn mạnh vào lòng máng gây xáo trộn bề mặt dữ dội. Quá trình này hình thành các cấu trúc xoáy cuộn ba chiều quy mô lớn. Những dòng xoắn thứ cấp xuất hiện liên tục và dịch chuyển về phía hạ lưu. Áp lực động tác động không đều lên các thành bên của máng dẫn. Mực nước phía thành đối diện luôn dâng cao hơn mực nước trung bình. Hiện tượng bất đối xứng thủy lực làm gia tăng mức độ tiêu hao năng lượng. Các cấu trúc xoáy phức tạp gây khó khăn cho việc đo đạc thực địa. Phân tích xoáy thứ cấp là cơ sở để xác định chính xác tải trọng thủy động tác dụng lên kết cấu máng.

Mô phỏng ba chiều phản ánh trung thực trường thủy động lực học trong máng tràn. Tuy nhiên, mô hình 3D đòi hỏi tài nguyên tính toán lớn và thời gian xử lý kéo dài. Thí nghiệm vật lý thu nhỏ cũng tốn kém nhiều chi phí và nhân lực. Do đó, việc xây dựng mô hình một chiều tiệm cận độ chính xác cao là giải pháp tối ưu. Bài toán 1D cần giản lược hình học nhưng phải bảo toàn đầy đủ các đặc trưng vật lý cốt lõi. Việc bổ sung các lực quán tính và hệ số hiệu chỉnh vận tốc là yêu cầu bắt buộc. Mô hình 1D cải tiến giúp rút ngắn thời gian tính toán mà vẫn đảm bảo độ tin cậy. Cách tiếp cận này đáp ứng kịp thời tiến độ thiết kế công trình trong thực tế.

Phương trình liên tục và động lượng tạo thành khung toán học hoàn chỉnh cho dòng biến lượng. Phương trình liên tục biểu diễn định luật bảo toàn khối lượng khi có lưu lượng bổ sung dọc tuyến. Số hạng nguồn phản ánh chính xác lượng nước gia nhập trên từng đơn vị chiều dài kênh. Phương trình động lượng mô tả sự thay đổi xung lượng của khối chất lỏng chuyển động. Thành phần gia tốc cục bộ và gia tốc đối lưu được kết hợp cùng gradien áp suất mặt nước. Lực ma sát đáy và thành kênh được xác định qua công thức thủy lực kinh nghiệm. Hệ phương trình đạo hàm riêng phi tuyến phản ánh trung thực tương tác lực phức tạp. Việc thiết lập đúng phương trình là bước quyết định đến tính chính xác của lời giải.

Dòng chảy tràn từ đỉnh đập rơi vào máng mang theo động lượng theo phương ngang và phương đứng. Lực quán tính của dòng chảy bên tác động mạnh vào dòng chính tạo ra độ chênh mực nước ngang. Khi máng tràn có đoạn uốn cong, lực quán tính ly tâm xuất hiện và làm biến dạng đường mặt nước. Lực ly tâm đẩy khối chất lỏng dồn về phía bờ ngoài của đoạn cong. Phương trình động lượng mở rộng phải bổ sung trực tiếp các thành phần lực này. Việc bỏ qua các lực quán tính sẽ gây sai số lớn trong việc tính toán chiều cao tường biên. Mô hình toán cải tiến tính đến góc gia nhập và vận tốc của dòng chảy bên. Kết quả tính toán đảm bảo an toàn kết cấu chống tràn nước qua thành máng.

Hệ phương trình vi phân phi tuyến của dòng biến lượng không thể giải bằng phương pháp giải tích thông thường. Phương pháp sai phân hữu hạn là công cụ hiệu quả để rời rạc hóa miền tính toán. Miền không gian dọc lòng dẫn và trục thời gian được chia thành lưới mắt lưới phù hợp. Các sơ đồ sai phân ẩn hoặc hiện được áp dụng để tuyến tính hóa các số hạng phi tuyến. Điều kiện biên và điều kiện ban đầu được thiết lập bám sát thực tế vận hành công trình. Thuật toán quét ma trận giúp giải nhanh hệ phương trình đại số tuyến tính tại mỗi bước thời gian. Chương trình tính toán tự động đảm bảo tính ổn định số học và tốc độ hội tụ nhanh. Độ chính xác của thuật toán được kiểm chứng qua số liệu đo đạc thực nghiệm.

Phân bố lưu tốc mặt cắt ngang trong máng tràn bên có tính dị hướng rất cao. Dòng nước gia nhập từ đỉnh tràn đẩy dòng chính lệch sang phía đối diện. Vùng vận tốc lớn nhất thường tập trung gần đáy hoặc dịch chuyển về phía tường xa dòng tràn. Gần bề mặt và khu vực tiếp nhận dòng bên xuất hiện các vùng xoáy có vận tốc nghịch đảo. Độ dốc vận tốc lớn tạo ra ứng suất tiếp ma sát nội tại rất mạnh trong lòng chất lỏng. Trường lưu tốc biến đổi liên tục dọc theo chiều dài của máng tràn. Dữ liệu thực nghiệm chỉ ra rằng biên dạng vận tốc thay đổi theo từng cấp lưu lượng xả lũ. Nắm bắt trường lưu tốc giúp tối ưu hóa hình học lòng dẫn và giảm mài mòn kết cấu.

Hệ số động lượng Boussinesq đại diện cho mức độ phân bố không đều của lưu tốc khi tính thông lượng động lượng. Trong dòng chảy đều thông thường, hệ số này xấp xỉ bằng một. Tuy nhiên, trong máng tràn bên, dòng chảy bị xáo trộn dữ dội khiến hệ số biến thiên đáng kể. Giá trị hệ số Boussinesq tăng cao tại những đoạn có cường độ dòng gia nhập lớn. Việc đưa hệ số này vào phương trình động lượng giúp hiệu chỉnh chính xác lực xung lượng của dòng chảy. Bỏ qua hệ số Boussinesq sẽ làm giảm độ chính xác khi xác định áp lực nước tác dụng lên đáy máng. Kết quả nghiên cứu cung cấp công thức thực nghiệm xác định hệ số cho từng đoạn máng tràn.

Hệ số động năng Coriolis phản ánh sai lệch về động năng thực tế so với động năng tính theo vận tốc trung bình. Sự xoáy cuộn và phân tầng vận tốc làm tăng trị số của hệ số Coriolis. Đại lượng này đóng vai trò quan trọng trong việc tính toán thế năng và động năng dọc tuyến máng. Tổn thất năng lượng cục bộ do hòa trộn dòng chảy có mối liên hệ mật thiết với hệ số này. Việc xác định đúng hệ số Coriolis giúp nâng cao độ chính xác của phương trình năng lượng Bernoulli mở rộng. Tính toán thủy lực trở nên chuẩn xác hơn khi thiết kế đoạn nối tiếp giữa máng tràn và dốc nước. Kiểm soát năng lượng dòng chảy bảo vệ công trình hạ lưu khỏi nguy cơ xói lở nguy hiểm.

Độ dốc thủy lực và đường mặt nước trong máng tràn biến đổi rất phức tạp dọc theo tuyến dòng chảy. Tại đầu máng, lưu lượng còn nhỏ nên mực nước thường thấp và dòng chảy duy trì trạng thái êm. Dọc theo chiều dài máng, lượng nước gia nhập liên tục làm tăng vận tốc và thay đổi gradien thủy lực. Đường mặt nước dâng cao dần về phía cuối máng trước khi đổ vào dốc nước hạ lưu. Sự thay đổi độ dốc đáy từ thoải sang dốc gây ra hiện tượng co hẹp và tăng tốc đột ngột. Tổn thất cột nước cục bộ do xáo trộn làm đường dốc năng lượng giảm nhanh. Mô phỏng chính xác đường mặt nước giúp kỹ sư xác định cao trình tường chắn sóng hợp lý.

Lưu lượng gia nhập dọc tuyến (lateral inflow) là nhân tố chính chi phối thủy động lực học máng tràn bên. Dòng nước tràn qua đỉnh đập đổ vào máng với vận tốc ban đầu vuông góc hoặc xiên góc với dòng chính. Dòng nhập lưu mang theo năng lượng và xung lượng làm biến đổi hoàn toàn cấu trúc dòng chảy. Cường độ gia nhập càng lớn thì hiện tượng dâng nước cục bộ tại thành đối diện càng thể hiện rõ. Áp lực dòng nhập lưu cản trở chuyển động dọc trục của dòng nước sẵn có trong lòng máng. Tổn thất năng lượng do va đập và hòa trộn tăng nhanh theo cường độ dòng gia nhập. Việc kiểm soát lưu lượng nhập lưu giúp giảm thiểu rung lắc thủy lực tác động lên công trình.

Lưu lượng tách dòng dọc tuyến (lateral outflow) xảy ra khi nước phân tán ra khỏi kênh chính qua các cửa lấy nước hoặc tràn bên phụ. Trái ngược với quá trình nhập lưu, dòng phân lưu làm giảm dần tổng lưu lượng dòng chính về phía hạ lưu. Sự suy giảm lưu lượng dẫn đến giảm vận tốc trung bình và làm biến dạng đường mặt nước trong kênh. Năng lượng mất đi chủ yếu theo khối lượng nước thoát ra ngoài mà không gây xáo trộn va đập mạnh như dòng nhập lưu. Tuy nhiên, hiện tượng tách dòng vẫn tạo ra các vùng xoáy cục bộ gần vị trí phân lưu. Mô hình toán Saint-Venant mở rộng xử lý hiệu quả cả hai trường hợp nhập lưu và tách dòng. Đánh giá phân lưu giúp tối ưu hóa hệ thống kênh phân phối nước tưới và thoát lũ.

Sóng lan truyền và biến dạng sóng lũ trong lòng máng tràn diễn ra dưới tác động mạnh mẽ của dòng nhập lưu. Đỉnh lũ khi di chuyển xuống hạ lưu có xu hướng bị kéo dãn và thay đổi sườn dốc sóng. Vận tốc truyền sóng phụ thuộc vào độ sâu mực nước và vận tốc trung bình của dòng chảy. Dòng nước tràn bổ sung liên tục làm đỉnh sóng tăng thêm lưu lượng và dâng cao mực nước đỉnh lũ. Sự ma sát với thành máng cùng hiệu ứng xoáy cuộn làm tiêu tán một phần năng lượng sóng lũ. Sóng động lực học trong máng tràn biến đổi nhanh hơn nhiều so với sóng lũ trên sông tự nhiên. Mô hình toán học mở rộng dự báo chính xác diễn biến sóng lũ qua từng mét chiều dài công trình.

Trong các mô hình thủy lực sông truyền thống, việc chia đoạn kênh quá dài thường bỏ qua độ trễ của các đặc trưng dòng chảy. Đối với máng tràn bên, dòng chảy có sự gia nhập liên tục và chiều dài ngắn nên độ trễ pha đóng vai trò quyết định. Đỉnh lưu lượng gia nhập và đỉnh mực nước tại mặt cắt khống chế luôn có độ lệch thời gian nhất định. Bỏ qua hiện tượng trễ pha sẽ dẫn đến tính toán sai lệch thời gian đạt đỉnh lũ công trình. Việc tích hợp các số hạng động lực học mở rộng giúp khắc phục hoàn toàn nhược điểm trễ pha. Độ tin cậy của mô hình số được nâng cao rõ rệt, đáp ứng tiêu chuẩn an toàn đập quốc tế. Kết quả nghiên cứu hỗ trợ đắc lực cho công tác cảnh báo sớm thiên tai và điều tiết hồ chứa.