Trong việc giải bài toán truyền chất một chiều, phương pháp sai phân hữu hạn có nhiều lợi thế. Sau việc giải bài toán thủy lực độc lập để xác định đặc trưng dòng chảy (lưu lượng, mực nước,.,) ta giải phương trình truyền chất bằng cách phân rã thành phương trình tải và phương trình khuếch tán. Nghiệm của phương trình tải đóng vai trò chính và phương trình khuếch tán được giải trên nền nghiệm này, nói chung chỉ vi chỉnh lại nghiệm đó để cho lời giải cuối cùng. | Xây dựng một sơ đồ tính truyền chất 1 D trên cơ sở hàm spline bậc 2 XÂY DỰNG MỘT SƠ ĐỒ TÍNH TRUYỀN CHẤT 1 D TRÊN CƠ SỞ HÀM SPLINE BẬC 2 GS. TSKH NGUYỄN ÂN NIÊN (1), THS. NCS BÙI VIỆT HƯNG (2) (1) Viện Khoa học Thuỷ lợi Miền Nam (2) Công ty Tư vấn xây dựng Thủy lợi 2 Tóm tắt : Trong việc giải bài toán truyền chất một chiều, phương pháp sai phân hữu hạn có nhiều lợi thế. Sau việc giải bài toán thủy lực độc lập để xác định đặc trưng dòng chảy ( lưu lượng, mực nước,.,) ta giải phương trình truyền chất bằng cách phân rã thành phương trình tải và phương trình khuếch tán. Nghiệm của phương trình tải đóng vai trò chính và phương trình khuếch tán được giải trên nền nghiệm này, nói chung chỉ vi chỉnh lại nghiệm đó để cho lời giải cuối cùng. Chúng tôi đề xuất việc giải phương trình tải dùng hàm spline bậc 2 để khắc phục một số nhược điểm của sơ đồ hiện dùng. I . MỞ ĐẦU Phương trình truyền chất một chiều với nồng độ c trên một đơn vị thể tích chất lỏng có dạng. c c 1 c q v AD c q c c 0 (1) t x A x x A Trong đó : v – vận tốc trung bình mặt cắt A – diện tích mặt cắt ướt phụ thuộc mực nước Z tại mặt cắt đó. D – hệ số khuếch tán vật lý ( phân tử + rối) q – lưu lượng bổ sung ngang trên một đơn vị chiếu dài dọc dòng chảy cq – nồng độ chất dòng bổ sung ngang - cường độ ( thể tích ) tăng / giảm nồng độ chất - hệ số mô men do sự phân bố không đều của nồng độ chất c’ và lưu tốc u’ tại các điểm mặt cắt. 1 uc dA (2) vAc A Để có tốc độ v, diện tích mặt cắt ướt A(x,z) và cả hệ số khuếch tán D (phụ thuộc vào các yếu tố thủy lực, chẳng hạn trong công thức của Fisher [1]) ta cần giải bài toán thủy lực trước để nhận được các đặc trưng mực nước z, lưu lượng Q tại các mặt cắt ở mọi lớp thời gian bằng các sơ đồ quen biết như VRSAP, TLUC, KOD01, MIKE 11 ,v,v. [2] Để giải phương trình truyền chất (1) ta phân rã thành phương trình tải : c c q v c q c c 0 (3) t
