Nghiên cứu một số kỹ thuật xử lý điều kiện biên trong giải số mô hình thuỷ lực hai chiều

ĐẠI HỌC QUỐC GIA HÀ NỘI VIỆN KHOA HỌC VÀ CÔNG NGHỆ VIỆT NAM TRƯỜNG ĐẠI HỌC CÔNG NGHỆ ****** VIỆN CƠ HỌC ***** NGUYỄN THÀNH ĐÔN NGHIÊN CỨU MỘT SỐ KỸ THUẬT XỬ LÝ ĐIỀU KIỆN BIÊN TRONG GIẢI SỐ MƠ HÌNH THUỶ LỰC HAI CHIỀU LUẬN VĂN THẠC SĨ HÀ NỘI-2005 ĐẠI HỌC QUỐC GIA HÀ NỘI VIỆN KHOA HỌC VÀ CÔNG NGHỆ VIỆT NAM TRƯỜNG ĐẠI HỌC CÔNG NGHỆ VIỆN CƠ HỌC ***** ***** NGUYỄN THÀNH ĐÔN NGHIÊN CỨU MỘT SỐ KỸ THUẬT XỬ LÝ ĐIỀU KIỆN BIÊN TRONG GIẢI SỐ MƠ HÌNH THUỶ LỰC HAI CHIỀU Chuyên ngành: Cơ học chất lỏng Mã số: 60.44.22 Luận văn thạc sĩ Ngƣời hƣớng dẫn khoa học: PGS.TS Hoàng Văn Lai HÀ NỘI-2005 Trang Lời cảm ơn Mục lục .2 Danh mục ký hiệu, chữ viết tắt Danh mục hình vẽ, đồ thị …………………………………………………………… Mở đầu …………………………………………………………………………………………… Chƣơng1 - GIỚI THIỆU CHUNG ………………………………………………………8 Chƣơng - CÁC KHÁI NIỆM CƠ BẢN 11 2.1 Hệ phƣơng trình Saint Venant ……………………………………………… 11 2.2 Số Froude số điều kiện biên cần thiết cho toán chiều hai chiều …………………………………………………………………………………………………13 2.2.1 Số Froude, số điều kiện biên toán thuỷ lực chiều 2.2.2 Số Froude, số điều kiện biên toán thuỷ lực hai chiều 21 2.3 Ý nghĩa vật lý điều kiện biên thuỷ lực học ……………………22 Chƣơng - KỸ THUẬT XỬ LÝ ĐIỀU KIỆN BIÊN TRONG MỘT SỐ PHƢƠNG PHÁP SỐ GIẢI BÀI TOÁN HAI CHIỀU ………………………………… 23 3.1 Kỹ thuật xử lý điều kiện biên trong phƣơng pháp khối hữu hạn 3.1.1 Phƣơng pháp rời rạc hố hệ phƣơng trình Saint Venant …… 25 3.1.2 Kỹ thuật xử lý phần tử biên …………………………………27 3.1.3 Kỹ thuật xử lý biên khô ƣớt ……………………………………………… 30 3.1.4 Ứng dụng kỹ thuật vào phần mềm VODAP_2D ……………32 3.2 Kỹ thuật xử lý điều kiện biên trong phƣơng pháp phần tử hữu hạn …………………………………………………………………………………………………………32 3.2.1 Phƣơng pháp rời rạc hố hệ phƣơng trình Saint Venant .32 3.1.2 Cách đƣa điều kiện biên vào hệ phƣơng trình 36 3.1.3 Kỹ thuật xử lý biên khô ƣớt ……………………………………………… 38 3.1.4 Ứng dụng kỹ thuật vào phần mềm TELEMAC_2D ……….39 Chƣơng - KẾT QUẢ GIẢI SỐ MỘT SỐ BÀI TOÁN MẪU ………………41 4.1 Bài tốn mẫu có nghiệm giải tích …………………………………………… 41 4.2 Bài tốn thí nghiệm có số liệu thực đo ……… 42 4.3 Bài tốn thí nghiệm có số liệu thực đo ……… 45 4.4 Bài tốn thực tế đánh giá thực trạng lịng dẫn sơng Hồng- sơng Thái Bình kiểm chứng ………………………………………………………………………………50 KẾT LUẬN VÀ KIẾN NGHỊ 54 DANH MỤC CƠNG TRÌNH CỦA TÁC GIẢ ……………………………………….55 TÀI LIỆU THAM KHẢO 56 DANH MỤC CÁC KÝ HIỆU, CÁC CHỮ VIẾT TẮT h - giá trị mực nƣớc u - vận tốc (trung bình) theo trục x v - vận tốc (trung bình) theo trục y g- gia tốc trọng trƣờng kx- hệ số Stricler lực cản đáy theo trục x ky- hệ số Stricler lực cản đáy theo trục y So,x - độ dốc đáy theo trục x So,y - độ dốc đáy theo trục y Fx - lực khối chiếu theo trục x Fy- lực khối chiếu theo trục y  e - hệ số khuyếch tán bao gồm khuyếch tán phân tử kết hợp khuyếch tán rối Sce - thành phần nguồn phụ Z - cao trình mặt thống Zf – cao trình đáy DANH MỤC CÁC HÌNH VẼ VÀ ĐỒ THỊ Hình 2.1 Một điều kiện biên thƣợng lƣu Hình 2.2 Hai điều kiện biên thƣợng lƣu Hình 2.3 Một điều kiện biên hạ lƣu Hình 2.4 Khơng cần điều kiện biên hạ lƣu Hình 3.1 Phần tử miền Hình 3.2 Phần tử biên Hình 3.3 Các phần tử nửa khô nửa ƣớt cần xử lý đặc biệt Hình 3.4 Các nút cần xử lý đặc biệt Hình 3.5 Các sửa gradient mặt thống lỗi Hình 4.1 So sánh mực nƣớc tính tốn với mực nƣớc giải tích Hình 4.2 Cấu hình kênh tốn mẫu số Hình 4.3 So sánh tính tốn với thực đo điểm đo S1 Hình 4.4 So sánh tính tốn với thực đo điểm đo S2 Hình 4.5 So sánh tính tốn với thực đo điểm đo S3 Hình 4.6 So sánh tính tốn với thực đo điểm đo S4 Hình 4.7 Cấu hình kênh tốn mẫu số Hình 4.8 So sánh tính tốn với thực đo điểm đo S1 Hình 4.9 So sánh tính tốn với thực đo điểm đo S2 Hình 4.10 So sánh tính tốn với thực đo điểm đo S3 Hình 4.11 So sánh tính tốn với thực đo điểm đo S4 Hình 4.12 So sánh tính tốn với thực đo điểm đo S5 Hình 4.13 So sánh tính tốn với thực đo điểm đo S6 Hình 4.14 Mơ hình hành lang lũ sơng Hồng Hình 4.15 Lƣới đƣợc chia chi tiết cơng trình đê, bối Hình 4.16 So sánh kết 2D 1D trạm thuỷ văn Hà Hình 4.17 Nội So sánh kết 2D 1D cầu Thăng Long Điều kiện biên phận cấu thành toán học chất lỏng Điều kiện biên lĩnh vực thuỷ động lực học lại mang số đặc trƣng chuyên sâu riêng so với ngành khác Hơn kiến thức kinh nghiệm xử lý điều kiện biên cịn giúp ngƣời tác nghiệp áp dụng có hiệu thuỷ động lực học vào thực tế Vì nghiên cứu nắm rõ điều kiện biên toán thuỷ lực nhiệm vụ cần thiết ngƣời làm thuỷ lợi Đề tài nghiên cứu ảnh hƣởng điều kiện biên tới kết số mơ hình hai chiều tâp trung vào vấn đề điều kiện biên phƣơng trình nƣớc nơng Saint Venant 2D Trong khuôn khổ luận văn đề cập giải thích số khái niệm, định nghĩa điều kiện biên toán chiều Luận văn mô tả ý nghĩa tác dụng loại điều kiện biên thực tế, yêu cầu số lƣợng điều kiện biên để toán có nghiệm Tuy nhiên để dẫn giải sáng sủa vấn đề, chƣơng hai đề tài đề cập đến kiến thức điều kiện biên toán chiều trƣớc Điều kiện biên hai chiều đƣợc lập luận tƣơng tự Khi nghiên cứu toán hai chiều truyền lũ, nhà thuỷ lực gặp câu hỏi làm mô tả đƣợc lan truyền nƣớc từ vùng ƣớt lên vùng khô, ngƣợc lại rút nƣớc Khi ta khơng cịn khái niệm mơi trƣờng liên tục tồn miền tính Khác với khí động học, tồn vùng nghiên cứu ln đƣợc lấp đầy khơng khí, thuỷ lực dâng nƣớc dẫn đến ngập vùng khô hay ngƣợc lại rút nƣớc từ vùng ƣớt thành vùng khô lại thƣờng xuyên xảy Vùng giáp ranh khô ƣớt lúc đƣợc coi biên lỏng di động chúng cần đƣợc nghiên cứu Loại điều kiện biên không đƣợc hiểu theo nghĩa thông thƣờng nhƣ loại điều kiện biên khác, nhƣng ý nghĩa quan trọng nó, đề tài cập đến loại điều kiện biên chƣơng riêng Chƣơng ba nêu định nghĩa xác định biên miền, nhƣ số tốn mẫu có lời giải để kiểm chứng Chƣơng cuối đƣa vài tốn mẫu có lời giải giải tích số liệu thực đo phịng thí nghiệm châu Âu đề xuất Chƣơng đƣa vài tốn thực tiễn mà nhóm tác giả thực thời gian vừa qua Kết số đƣợc so sánh với kết mẫu nhằm chứng minh vấn đề mà luận văn đặt Tuy kết số giá trị trung bình đơi chỗ cịn khác so với kết thực đo, nhƣng nhìn tổng thể kết đạt tiêu chuẩn cho phƣơng pháp số Bản thân lý thuyết điều kiện biên hệ phƣơng trình Saint-Venant 2D đƣợc phát triển nhiều hệ khoa học Do vậy, đề tài nhằm mục tiêu nêu lại lý thuyết cách áp dụng chúng vào thực tiễn cho đảm bảo tính chặt chẽ hiệu đáp ứng đƣợc toán thực tế đặt Chƣơng – GIỚI THIỆU CHUNG Kỹ thuật xử lý điều kiện biên cho phƣơng trình Saint Venant đƣợc nhiều hệ nhà khoa học chuyên ngành toán học lẫn học quan tâm từ lâu Trong nƣớc có PGS.TS Trần Gia Lịch, GS.TSKH Nguyễn Kim Đan, PGS.TS Hoàng Văn Lai nghiên cứu có nhiều báo đăng tạp chí uy tín vấn đề Ở nƣớc ngồi có nhiều nhà khoa học nghiên cứu hoàn thiện kỹ thuật xử lý điều kiện biên cho toán thuỷ lực học Sau sơ lƣợc tình hình nghiên cứu tác giả nƣớc GS Nguyễn Kim Đan công tác đại học tổng hợp Caen nghiên cứu chuyên sâu phƣơng pháp số giải hệ phƣơng trình SaintVenant 2D kỹ thuật xử lý biên khô ƣớt Các kỹ thuật quan trọng tốn vỡ đê, vỡ đập, lan truyền lũ v.v Giáo sƣ ngƣời hƣớng dẫn nhiều nghiên cứu sinh cán Việt nam vấn đề Phƣơng pháp phần mềm giáo sƣ viết đƣợc ứng dụng Việt nam PGS.TS Hoàng Văn Lai nghiên cứu kỹ thuật xử lý biên gián đoạn Kết tính tốn số chƣơng trình PGS Hồng Văn Lai xây dựng vƣợt qua toán mẫu phịng thí nghiệm thuỷ lực châu Âu đƣa GS.TS Trần Gia Lịch TS Lê Kim Luật viết báo điều kiện biên, hai ngƣời chứng minh để tồn nghiệm tốn tuyến tính hố, điều kiện biên phải thoả mãn vài bất đẳng thức liên hệ Bài báo có ý làm chặt chẽ theo nghĩa tốn học phƣơng pháp tuyến tính hố Tuy nhiên, báo đƣa vài luận đề toán học làm sở mà không chứng minh Trong tốn thực tế q trình lan truyền lũ, việc tìm giá trị đại lƣợng biên quan trọng Vì ngƣời ta xây dựng số phần mềm tính giá trị đại lƣợng từ lƣợng mƣa lƣu vực Q trình hình thành dịng chảy từ lƣợng mƣa rơi lƣu vực trình phức tạp, phụ thuộc vào nhiều yếu tố: độ dốc, độ che phủ lƣu vực, thành phần cấu tạo đất, lƣợng bốc hơi….Mơ hình mƣa rào dịng chảy đƣợc xây dựng dựa sở giả thiết chấp nhận số thông số đặc trƣng cho lƣu vực Các thơng số đƣợc lựa chọn thuật tốn tối ƣu hoá dựa số liệu thực đo trƣớc thời điểm cần tính tốn Mơ hình thuỷ văn mƣa rào dòng chảy đƣợc xây dựng dƣới đạo GS.TS Trịnh Quang Hồ có khả tính tốn dịng chảy sinh mƣa lƣu vực Trong xây dựng mơ hình thuỷ văn mƣa rào dòng chảy phải chấp nhận nhiều thơng số thực nghiệm cho lƣu vực Vì q trình hình thành dịng chảy lƣu vực phụ thuộc nhiều vào yếu tố lƣu vực: địa hình, độ che phủ, cấu tạo đất….Do vậy, việc xác định thông số đặc trƣng lƣu vực cho mơ hình thuỷ văn mƣa rào dịng chảy khó khăn cho độ xác khơng cao Với mục đích mơ xác q trình hình thành dịng chảy lƣu vực, thời gian gần nhiều nhà thuỷ văn, thuỷ lực cố gắng xây dựng mơ hình thuỷ văn sử dụng thành tựu lĩnh vực thông tin địa lý (GIS) Một mơ hình loại mơ hình MARINE (Modelisation de l’Anticipation du Ruissellement et des Inondations pour des événements) Viện Cơ học chất lỏng Toulouse (IMFT – Institut de Mecanique de Fluides de Toulouse) phát triển Trong khuôn khổ đề tài nghiên cứu khoa học công nghệ KC.08-13 43 đƣợc kéo lên gần nhƣ tức thời mơ q trình vỡ đập hoàn toàn tức thời Thời gian mở cửa sập xấp xỉ 0.2 giây Cửa sập đặt vị trí 6.10 mét tính từ vị trí thƣợng nguồn kênh Độ rộng đoạn kênh thƣợng lƣu 0.5 mét Đáy thành phải kênh đƣợc sơn mực đậm mầu, thành trái làm thuỷ tinh Plexiglass suốt Chiều cao thành 0.50 mét Vị trí nút thắt đƣợc đặt vị trí 7.70 mét tính từ vị trí cửa sập tức 13.80 mét tính từ thƣợng nguồn Đoạn cơng trình thắt kênh dài l.0 mét rộng 0.1 mét Đoạn cơng trình chuyển tiếp làm góc 45 độ với tƣờng kênh Q trình mực nƣớc theo thời gian (hydrograph) đƣợc đo đạc bốn điểm đo S1 tới S4 đặt vị trí nhƣ hình 4.2 Hình 4.2 : Cấu hình kênh S1 S2 S3 S4 44 Sau biểu đồ so sánh kết tính tốn số trị chƣơng trình VODAP_2D, TELEMAC_2D kết thực đo phịng thí nghiệm Trên biểu đồ, đƣờng đỏ đƣờng tính tốn, đƣờng xanh đƣờng thực đo Kết tính tốn thực đo điểm đo S1 (hình 4.3) 0.35 0.35 0.3 0.3 0.25 0.25 0.2 0.2 Series1 Series1 Series2 Series2 0.15 0.15 0.1 0.1 0.05 0.05 0 10 12 Kết TELEMAC 10 12 Kết VODAP_2D Hình 4.3: Mực nƣớc thực đo mực nứơc tính tốn điểm đo S1 Kết tính tốn thực đo điểm đo S2 (hình 4.4) 0.25 0.3 0.25 0.2 0.2 0.15 Series1 Series1 0.15 Series2 Series2 0.1 0.1 0.05 0.05 0 10 12 10 12 45 Kết TELEMAC Kết VODAP_2D Hình 4.4: Mực nƣớc thực đo mực nứơc tính tốn điểm đo S2 Kết tính tốn thực đo điểm đo S3 (hình 4.5) 0.16 0.16 0.14 0.14 0.12 0.12 0.1 0.1 Series1 0.08 Series1 Series2 0.08 Series2 0.06 0.06 0.04 0.04 0.02 0.02 0 10 12 Kết TELEMAC 10 12 Kết VODAP_2D Hình 4.5: Mực nƣớc thực đo mực nứơc tính tốn điểm đo S3 Kết tính tốn thực đo điểm đo S4 (hình 4.6) 0.04 0.04 0.035 0.035 0.03 0.03 0.025 0.025 Series1 0.02 Series2 0.015 0.015 0.01 0.01 0.005 0.005 Series1 Series2 0.02 0 Kết TELEMAC 10 12 10 Kết VODAP_2D Hình 4.6: Mực nƣớc thực đo mực nứơc tính tốn điểm đo S4 12 46 Chúng ta thấy đƣờng xanh thực đo dao động, điểm đo xa nơi vỡ đập nhiễu dao động nhiều với biên độ lớn Đó kết q trình tích luỹ sóng kích động tồn q trình truyền sóng Đƣờng kết số trị mầu đỏ khơng mơ tả đƣợc q trình lan truyền sóng kích động Bởi lẽ thân mơ hình thuật tốn bắt gián đoạn khơng dùng để mơ tả q trình Tuy nhiên đƣờng đỏ bắt tốt đƣờng thực đo pha giá trị trung bình giá trị mực nƣớc 4.3 Bài tốn thí nghiệm có số liệu thực đo số Tên tốn: Vỡ đập kênh có miền hạ du mở rộng Mục tiêu: Kiểm tra việc xử lý biên khô/ƣớt, so sánh với kết tính tốn với liệu thu đƣợc từ thí nghiệm Trƣờng hợp kiểm định đặc biệt phù hợp để kiểm tra tính xác mơ hình số Cả mơ hình D lẫn D đƣợc kiểm tra thí nghiệm Mơ tả tốn: Hình dạng kênh đƣợc mơ nhƣ hình 4.7 - Kênh dài 19.3 mét, độ dốc đáy Đoạn kênh nhỏ hình chữ nhật rộng 0.5mét - Đoạn kênh lớn có độ rộng 2.3 mét chiều dài 6.75 mét - Vị trí mở rộng đƣợc đặt vị trí cách cửa sập 6.45 mét tức cách vị trí thƣơng lƣu 12.55 mét - Q trình vỡ đập đƣợc coi vỡ hoàn toàn tức thời với thời gian vỡ (mở cửa sập) 0.2 giây Cửa sập đƣợc đặt vị trí 6.10 mét tính từ thƣợng lƣu kênh Chiều rộng kênh đoạn thƣợng lƣu 0.5mét - Đáy thành phải đƣợc xây dựng ximăng mịn, quét sơn đậm màu, thành trái làm thuỷ tinh Plexiglass suốt Độ cao thành 0.5 mét 47 - Quá trình mực nƣớc theo thời gian thu đƣợc phƣơng pháp số trị đƣợc so sánh với kết thí nghiệm Hình 4.7: Cấu hình kênh Sau biểu đồ so sánh kết tính tốn kết thực đo toán Trên biểu đồ đƣờng đỏ đƣờng tính tốn, đƣờng xanh đƣờng thực đo Kết tính tốn thực đo điểm đo S1 (hình 4.8) 0.18 0.18 0.16 0.16 0.14 0.14 0.12 0.12 0.1 0.1 Series1 Series1 Series2 Series2 0.08 0.08 0.06 0.06 0.04 0.04 0.02 0.02 0 10 Kết TELEMAC 12 10 12 Kết VODAP_2D Hình 4.8: Mực nƣớc thực đo mực nứơc tính tốn điểm đo S1 48 Kết tính tốn thực đo điểm đo S2 (hình 4.9) 0.025 0.03 0.025 0.02 0.02 0.015 Series1 0.015 Series1 Series2 Series2 0.01 0.01 0.005 0.005 0 10 12 Kết TELEMAC 10 12 Kết VODAP_2D Hình 4.9: Mực nƣớc thực đo mực nứơc tính tốn điểm đo S2 Kết tính tốn thực đo điểm đo S3 (hình 4.10) 0.05 0.05 0.045 0.045 0.04 0.04 0.035 0.035 0.03 0.03 Series1 0.025 Series2 0.02 0.02 0.015 0.015 0.01 0.01 0.005 0.005 Series1 Series2 0.025 0 Kết TELEMAC 10 12 10 Kết VODAP_2D 12 49 Hình 4.10: Mực nƣớc thực đo mực nứơc tính tốn điểm đo S3 Kết tính tốn thực đo điểm đo S4 (hình 4.11) 0.04 0.04 0.035 0.035 0.03 0.03 0.025 0.025 Series1 0.02 Series1 Series2 0.02 Series2 0.015 0.015 0.01 0.01 0.005 0.005 0 10 12 Kết TELEMAC 10 12 Kết VODAP_2D Hình 4.11: Mực nƣớc thực đo mực nứơc tính tốn điểm đo S4 Kết tính tốn thực đo điểm đo S5 (hình 4.12) 0.09 0.09 0.08 0.08 0.07 0.07 0.06 0.06 0.05 0.05 Series1 Series2 Series1 Series2 0.04 0.04 0.03 0.03 0.02 0.02 0.01 0.01 0 10 12 10 12 50 Kết TELEMAC Kết VODAP_2D Hình 4.12: Mực nƣớc thực đo mực nứơc tính tốn điểm đo S5 Kết tính tốn thực đo điểm đo S6 (hình 4.13) 0.045 0.045 0.04 0.04 0.035 0.035 0.03 0.03 0.025 0.025 Series1 Series1 Series2 Series2 0.02 0.02 0.015 0.015 0.01 0.01 0.005 0.005 0 Kết TELEMAC 10 12 10 Kết VODAP_2D Hình 4.13: Mực nƣớc thực đo mực nứơc tính tốn điểm đo S6 Chúng ta thấy đƣờng xanh thực đo dao động, điểm đo xa nơi vỡ đập nhiễu dao động nhiều với biên độ lớn Đó kết q trình tích luỹ sóng kích động tồn q trình truyền sóng Đƣờng kết số trị mầu đỏ không mô tả đƣợc q trình lan truyền sóng kích động Bởi lẽ thân mơ hình thuật tốn bắt gián đoạn khơng dùng để mơ 12 51 tả q trình Tuy nhiên đƣờng đỏ bắt tốt đƣờng thực đo pha giá trị trung bình giá trị mực nƣớc Chỉ có giá trị tính tốn vị trí S5 xa so với giá trị thực đo 4.4 Lời giải số toán trình lan truyền lũ vùng hệ thống sơng Hồng-Thái Bình Để đánh giá khả lũ hành lang thoát lũ đoạn từ Sơn Tây đến Hà nội, mơ hình vật lý đƣợc xây dựng trung tâm động lực sông Viện khoa học thuỷ lợi Bộ chƣơng trình TELEMAC_2D đƣợc sử dụng để nghiên cứu trình lan truyền lũ đoạn sơng Kừt tính tốn TELEMAC_2D đƣợc so sánh với kết tính tốn mơ hình tính tốn thuỷ lực chiều Mơ tả tốn: Miền tính tốn kéo dài dƣới 40 km từ Sơn Tây qua Hà Nội tới Vân Giang Hƣng Yên Trên đoạn sơng có hợp lƣu, dịng chảy chia thành hai nhánh sông Đuống sông Hồng Biên vào miền biên lƣu lƣợng Q Sơn Tây Các giá trị biên đƣợc tính từ chƣơng trình thuỷ lực chiều mặt cắt SHG_70 Biên miền biên cao trình mực nƣớc Z hai nơi biên sông Đuống biên Vân Giang Các giá trị biên đƣợc tính từ chƣơng trình thuỷ lực chiều mặt cắt SDUONG_05 SHG_104 Hệ số nhám Strickler lòng chảy sơng đƣợc để 70, hệ số nhám ngồi lịng chảy đƣợc để 30 Lƣới tính tốn đƣợc chia dầy vùng lịng dẫn thƣa vùng bối bãi 52 Biên sông Đuống Biên vào Sơn Tây Cầu Thăng Long Cầu Chương Dương trạm thuỷ văn Hà nội Biên Vân Giang Hình 4.14: Mơ hình hành lang lũ sơng Hồng Hệ thống đê bối cơng trình mơ chi tiết với nút lới đặt trực tiếp lên Hình 4.15: Lƣới đƣợc chia chi tiết cơng trình đê, bối 53 Sau kết tính tốn số mặt cắt: Tram Thuy Van Ha noi - Cau Chuong Duong 14 12 Muc nuoc 10 13 25 37 49 61 73 85 97 109 121 133 145 157 169 181 193 205 217 229 241 253 265 277 289 301 313 325 337 349 361 373 385 397 409 421 Thoi gian Muc nuoc tinh toan bang TL 1D Muc nuoc tinh toan bang chuong trinh 2D Hình 4.16: So sánh kết tính tốn mơ hình 2D mơ hình 1D trạm thuỷ văn Hà Nội Muc nuoc tai cau Thang Long 16 14 12 Muc nuoc 10 13 25 37 49 61 73 85 97 109 121 133 145 157 169 181 193 205 217 229 241 253 265 277 289 301 313 325 337 349 361 373 385 397 409 421 Muc nuoc tinh toan bang TL 1DThoi gian Muc nuoc tinh toan bang chuong trinh 2D Hình 4.17: So sánh kết tính tốn mơ hình 2D mơ hình 1D cầu Thăng Long 54 Kết thu đƣợc từ mô hình hai chiều sai khác tối đa kết thu đƣợc từ mơ hình chiều 30 cm Đỉnh hai đƣờng mực nƣớc sai khác 10 cm Nhƣ kết thu đƣợc từ mơ hình hai chiều đạt độ tin cậy, khối lƣợng tính tốn nhƣ độ phức tạp xử lý lớn nhiều mơ hình chiều 55 KẾT LUẬN VÀ KIẾN NGHỊ KẾT LUẬN Lý thuyết điều kiện biên đƣợc nghiên cứu sâu rộng, đầy đủ, đƣợc ứng dụng hiệu vào thực tiễn Các dạng điều kiện biên đƣợc trình bày đầy đủ, có khả mơ hình hố đƣợc tất tốn lý thuyết, đáp ứng đƣợc nhiều trƣờng nghiên cứu Đặc biệt dạng điều kiện biên khơ ƣớt đƣợc xử lý hồn chỉnh Các kết thu đƣợc có sai số nằm khuôn khổ cho phép Việc xây dựng điều kiện biên cho toán dự báo cịn nhiều khó khăn hạn chế Có phƣơng pháp xây dựng thiếu sở khoa học, mang thiên hƣớng kinh nghiệm Số liệu cho toán dự báo cịn nhiều thiếu sót Kết tốn dự báo nhiều sai số, đặc biệt trƣờng hợp có nhiều yếu tố khơng ổn định tác động Khơng có mơ hình tốn học mô đƣợc tất tƣợng vật lý Vì ta phải biết lựa chọn mơ hình dùng để nghiên cứu tƣơng vật lý xét KIẾN NGHỊ Độ phức tạp nhƣ khối lƣợng tính tốn tốn hai chiều lớn nên thời gian tính tốn thƣờng lâu Vì đề nghị quan chức cho phép nâng cấp cung cấp phƣơng tiện tính tốn mới, mạnh để giảm bớt cơng sức ngƣời chạy chƣơng trình nhƣ tăng hiệu việc tính tốn Để thực đƣợc nhiệm vụ dự báo, ngƣời làm thuỷ lực cần có số liệu hồn chỉnh Vì đề nghị quan chức hoàn thiện tập số liệu đo đạc để ngƣời làm thuỷ lực xoay xở 56 TÀI LIỆU THAM KHẢO [1] A reappraisal of Taylor-Galerkin algorithm for drying-wetting areas in shallow water computations International Journal for Numercal Methods in Fluids M.Quecedo and M.Pastor [2] Numerical methodes for shallow-water flow by C.B Vreugdenhil Institute for Marine and Atmospheric Research Utrecht(IMAU), Utrecht University, Utrecht, the Netherians [3] Donea J.A Taylor-Galerkin method for convective transport problem International Journal for Numerical Methods in Engineering 1984 [4] Peraire J.A Finite Element Method for Convection Dominated Flows University of Wales,Swansea,1986 [5] An unstructured finite-volume algorithm for predicting flow in rivers and estuaries by P.A.Sleigh, P.H.Gaskell, M.Berzins and N.G.Wright University of Leeds, Leeds, LS2 9JT U.K [6] A 2-D Shalow-Water Model using Unstructured finite volume methods by Nguyễn Kim Đan University of Caen, Department of Mathematics and Mechanics, France [7] Boundary conditions for the two-dimensional flow, Saint-Venant equation system by Trần Gia Lịch and Lê Kim Luật Institut of Mathematics, HàNội, Việtnam [8] TELEMAC-2D Vesion 3.0 Principle note by Hervouet J-M,Van Haren L [9] TELEMAC-2D Modeling system user manual 57 Departerment Laboratoire National d’Hydraulique FRANC