H Hj ΔH1 ΔH1/Δt tj Δt tj+1 t H Hj ΔH2 ΔH2/Δt tj tj+Δt/3 tj+1 Hj+ΔH1/3 t H Hj ΔH3 ΔH3/Δt tj tj+2Δt/3 tj+1 Hj+2ΔH2/3 t
2.3.2. Xây dựng mô đun diễn toán dòng chảy qua hồ
Mô đun diễn toán dòng chảy qua hồ được xây dựng bằng ngôn ngữ lập trình Fortran 70 bằng cách giải phương trình liên tục theo phương pháp Runge - Kutta bậc 3 và sơ đồ giải như Hình 2.22. Từ phương pháp Runge - Kutta bậc 3 tính toán được mực nước hồ ở bước thời gian tiếp theo và biến đổi mực nước trong hồ, từ đó tính được dung tích hồ ở bước thời gian tiếp theo và biến đổi dung tích hồ thông qua đường quan hệ hoặc bảng tra mực nước với dung tích hồ V = f(Z). Lưu lượng về hạ du hồ ở bước thời gian tiếp theo được tính dựa trên dung tích ở bước thời gian tiếp theo, biến đổi dung tích hồ trong bước thời gian Δt, lưu lượng về hồ và thông số cửa xả.
Số liệu đặc tính lòng hồ được thể hiện theo mối quan hệ dung tích, diện tích mặt nước theo cao trình mặt nước hay còn gọi là đường đặc tính VFZ được sử dụng để xác định quan hệ V = f(Z). Thông số cửa xả bao gồm: cao trình cửa xả, chiều dài mỗi cửa, số lượng cửa, hình dạng cửa được mô đun sử dụng để tính toán lưu lượng về hạ du hồ. Ngoài ra, quá trình diễn toán còn phụ thuộc mực nước ban đầu trong hồ, cao trình mực nước cao nhất. Công thức tổng quát tính lưu lượng qua đập tràn như sau :
𝑄 = 𝑚𝑏√2𝑔𝐻03⁄2 (2.3.10)
Trong đó : m là hệ số lưu lượng, b là chiều dài cửa tràn, g là gia tốc trọng trường và Ho là chênh lệch đô cao mực nước hồ và ngưỡng tràn.
Với đập tràn có nhiều cửa xả, công thức tính lưu lượng qua đập như sau:
𝑄 = 𝛿𝑛𝜀𝑚 ∑ 𝑏√2𝑔𝐻03⁄2 𝑛
𝑖=1
(2.3.11) Trong đó : ε là hệ số co hẹp, δn là hệ số thực nghiệm, n là số cửa xả tràn.
𝜀 = 1 − 0.2𝜀𝑚𝑏+ (𝑛 − 1)𝜀𝑚𝑡 𝑛
𝐻0
𝑏 (2.3.12)
Trong đó : εmb và εmt là các hệ số phụ thuộc vào hình dạng cửa, được xác định như Hình 2.23.
Các loại đập với cấu tạo đơn giản, dễ làm, có hệ số lưu lượng nhỏ, m = 0,35 ÷ 0,45, thường được dùng nhiều trong các công trình thủy lợi loại nhỏ bằng vật liệu tại chỗ: đá, gạch, gỗ .... Hệ số lưu lượng của đập mặt cắt chữ nhật, theo Badanh như sau:
𝑚 = 0.42 (0.7 + 0.185𝐻