Đại cương về mô hình thủy động lực và nguyên lý cơ bản
Mô hình thủy động lực (Hydrodynamic model) là công cụ toán học và tính toán mô phỏng chuyển động của chất lỏng (chủ yếu là nước), sự phân bố áp suất, trường vận tốc, mực nước và sự vận chuyển chất tan hoặc phù sa trong các thủy vực tự nhiên như sông ngòi, kênh dẫn, hồ chứa, vùng cửa sông và đại dương. Mô hình này là trụ cột khoa học then chốt trong kỹ thuật tài nguyên nước, hải văn học, quản lý rủi ro lũ lụt và quy hoạch công trình ven biển.
Cơ sở lý thuyết của mọi mô hình thủy động lực học dựa trên các định luật bảo toàn vật lý cơ bản: định luật bảo toàn khối lượng (phương trình liên tục), định luật bảo toàn động lượng (định luật II Newton) và định luật bảo toàn năng lượng. Đối với dòng chảy chất lỏng không nén được, các quá trình này được mô tả tổng quát bằng hệ phương trình vi phân đạo hàm riêng Navier-Stokes.
Phân loại các mô hình thủy động lực theo số chiều không gian
Tùy thuộc vào quy mô không gian và tính chất động lực học của dòng chảy, các mô hình thủy động lực được phân thành các cấp độ tính toán:
- Mô hình 1D (Một chiều): Điển hình là hệ phương trình Saint-Venant 1D, giả thiết các thông số dòng chảy (vận tốc, mực nước) đồng nhất trên toàn bộ mặt cắt ướt và chỉ biến thiên dọc theo chiều dài lòng dẫn. Phù hợp cho mạng lưới sông dài và hệ thống kênh rạch thủy lợi.
- Mô hình 2D dòng nông (Two-Dimensional Shallow Water Equations - SWE): Tích phân hệ phương trình Navier-Stokes theo chiều sâu khi chiều sâu dòng chảy rất nhỏ so với kích thước bề ngang. Mô hình 2D mô phỏng chính xác hiện tượng ngập lụt lan truyền trên bãi tràn, dòng chảy qua đê vỡ và thủy triều vùng cửa sông ven biển.
- Mô hình 3D (Ba chiều): Mô phỏng chi tiết sự biến thiên vận tốc theo cả phương ngang và phương thẳng đứng, tính toán hiện tượng phân tầng mật độ (do nhiệt độ và độ mặn), đối lưu xoáy và dòng thứ cấp tại các khúc uốn sâu hoặc hồ chứa lớn.
Phương pháp số và quy trình thiết lập mô hình
Vì hệ phương trình thủy động lực học là hệ phi tuyến không có nghiệm giải tích tổng quát, các phương pháp xấp xỉ số học được ứng dụng để rời rạc hóa miền tính toán:
- Phương pháp sai phân hữu hạn (Finite Difference Method - FDM): Rời rạc hóa đạo hàm trên lưới có cấu trúc hình chữ nhật, tính toán nhanh nhưng kém linh hoạt với đường bờ có hình học phức tạp.
- Phương pháp thể tích hữu hạn (Finite Volume Method - FVM): Bảo toàn thông lượng khối lượng và động lượng trên từng ô lưới đa giác không cấu trúc (tam giác, tứ giác), là phương pháp chủ đạo trong các phần mềm hiện đại như TELEMAC, Delft3D, MIKE 21.
- Phương pháp phần tử hữu hạn (Finite Element Method - FEM): Cho phép mô tả chính xác các biên uốn lượn phức tạp và thích nghi mật độ lưới cục bộ tại các vùng dòng chảy xiết.
Quy trình hiệu chỉnh, kiểm định và ứng dụng thực tiễn
Độ tin cậy của mô hình thủy động lực phụ thuộc quyết định vào dữ liệu địa hình đáy và quy trình hiệu chuẩn tham số nhám:
| Giai đoạn thực hiện | Nội dung chuyên môn | Chỉ số đánh giá độ tin cậy |
|---|---|---|
| Thu thập dữ liệu đầu vào | Đo đạc mô hình số độ cao (DEM/bathymetry), chuỗi số liệu thủy văn biên lưu lượng (thượng lưu) và mực nước triều (hạ lưu) | Kiểm tra sai số lưới và độ trơn của dữ liệu địa hình đáy |
| Hiệu chỉnh mô hình (Calibration) | Điều chỉnh hệ số nhám lòng dẫn (Manning n hoặc Chezy C) sao cho đường quá trình mực nước tính toán khớp với số liệu thực đo | Hệ số hiệu quả Nash-Sutcliffe (NSE > 0,8), hệ số tương quan R² |
| Kiểm định mô hình (Validation) | Chạy mô hình với một trận lũ hoặc mùa triều độc lập khác mà không thay đổi bộ thông số nhám đã hiệu chuẩn | Sai số đỉnh lũ (Peak error) dưới 10%, sai số thời gian xuất hiện đỉnh dưới 1 giờ |
| Ứng dụng kịch bản | Dự báo ngập lụt đô thị, mô phỏng xói lở bờ sông, tính toán lan truyền xâm nhập mặn và đánh giá tác động biến đổi khí hậu | Xuất bản đồ rủi ro ngập lụt hỗ trợ cơ quan phòng chống thiên tai |