ĐÁNH GIÁ PHẦN mềm THỦY lực BA CHIỀU để mô PHỎNG DÒNG CHẢY XUNG QUANH CÔNG TRÌNH kè TRÊN SÔNG WAAL, hà LAN LUẬN văn THẠC sĩ

ĐẠI HỌC ĐÀ NẴNG TRƯỜNG ĐẠI HỌC BÁCH KHOA NGUYỄN QUANG BÌNH - NGUYỄN QUANG BÌNH CHUYÊN NGÀNH: KỸ THUẬT XÂY DỰNG CƠNG TRÌNH THỦY ĐÁNH GIÁ PHẦN MỀM THỦY LỰC BA CHIỀU ĐỂ MƠ PHỎNG DỊNG CHẢY XUNG QUANH CƠNG TRÌNH KÈ TRÊN SƠNG WAAL, HÀ LAN LUẬN VĂN THẠC SĨ KỸ THUẬT XÂY DỰNG CƠNG TRÌNH THỦY KHỐ: 2015 - 2017 Đà Nẵng – Năm 2017 ĐẠI HỌC ĐÀ NẴNG TRƯỜNG ĐẠI HỌC BÁCH KHOA - NGUYỄN QUANG BÌNH ĐÁNH GIÁ PHẦN MỀM THỦY LỰC BA CHIỀU ĐỂ MÔ PHỎNG DỊNG CHẢY XUNG QUANH CƠNG TRÌNH KÈ TRÊN SƠNG WAAL, HÀ LAN Chuyên ngành : KỸ THUẬT XÂY DỰNG CÔNG TRÌNH THỦY Mã số : 60.58.02.02 LUẬN VĂN THẠC SĨ NGƯỜI HƯỚNG DẪN KHOA HỌC: PGS.TS Nguyễn Thống Đà Nẵng – Năm 2017 LỜI CAM ĐOAN Tôi xin cam đoan cơng trình nghiên cứu riêng tơi Các số liệu, kết nêu luận văn trung thực chưa công bố cơng trình khác Tác giả luận văn ký ghi rõ họ tên Nguyễn Quang Bình ĐÁNH GIÁ PHẦN MỀM THỦY LỰC BA CHIỀU ĐỂ MƠ PHỎNG DỊNG CHẢY XUNG QUANH CƠNG TRÌNH KÈ TRÊN SƠNG WAAL, HÀ LAN Học viên: Nguyễn Quang Bình Chuyên ngành: Kỹ thuật Xây dựng Cơng trình thủy Mã số: 60.58.02.02 Khóa: 2015 - 2017 Trường Đại học Bách khoa - ĐHĐN Tóm tắt - TELEMAC 3D FLOW 3D công cụ mơ hình hóa mạnh mẽ để mơ xác dịng chảy bề mặt thống Các mơ hình có khả cung cấp cho người lập mơ hình thơng tin chi tiết q trình vật lý Để so sánh tính hiệu hai mơ hình thực tế, nghiên cứu thực cách áp dụng để đánh giá thay đổi dòng chảy tràn xung quanh kè sông Waal Dựa kết mô phỏng, báo phân tích khác biệt hai mơ hình thơng qua vận tốc, hướng dòng chảy, cường độ rối, ứng suất tiếp, thời gian tính tốn, cơng cụ hỗ trợ Những so sánh dự kiến giúp người lập mơ hình có tài liệu việc lựa chọn mơ hình thích hợp cho nghiên cứu Từ khóa - Mơ hình số; TELEMAC 3D; FLOW 3D; Kè; Sơng Waal ASSESS THE 3D HYDRAULIC SOFTWARE FOR SIMULATING FLOW AROUND GROYNES IN THE WAAL RIVER, NETHERLANDS Abstract - TELEMAC 3D and FLOW 3D are powerful modeling tool for accurately simulating free surface flow These models have strong capacity to provide the modeler with valuable insights about physical process In order to compare the effectiveness of the two models in reality, this study is realized (conducted) by applying them to evaluate the variation of flow component over and around groynes in Waal river Based on the simulated results, the paper will analyze the differences between the two models through velocity, flow direction, turbulence intensity, shear stress, computational calculation, supporting tool… These comparisons are expected to help modelers having basic evidences in choosing suitable model for their study Key words - Numerical model; TELEMAC 3D; FLOW 3D; Groynes; Waal river MỤC LỤC MỞ ĐẦU … 1 Lý chọn đề tài Mục đích nhiệm vụ nghiên cứu Đối tượng phạm vi nghiên cứu Phương pháp nghiên cứu Ý nghiã khoa ho ̣c và thực tiễn Chương - TỔNG QUAN VỀ DÒNG CHẢY TRONG SƠNG 1.1 Dịng chảy chiều 1.1.1 Phạm vi áp dụng 1.1.2 Phương trình 1.2 Dòng chảy hai chiều 1.2.1 Phạm vi áp dụng 1.2.2 Phương trình đạo 1.3 Dòng chảy ba chiều 1.3.1 Phạm vi áp dụng 1.3.2 Phương trình đạo 1.3.3 Tổng quan chung 1.4 Dịng chảy xung quanh cơng trình kè 1.5 Kết luận Chương – CƠ SỞ LÝ LUẬN VỀ MƠ HÌNH THỦY LỰC 10 2.1 Khái niệm mơ hình 10 2.1.1 Mơ hình hóa 10 2.1.2 Mơ hình vật lý 10 2.1.3 Mơ hình thủy lực 10 2.2 Phân loại mơ hình thủy lực 11 2.3 Tương tự học 12 2.3.1 Tương tự hình học 12 2.3.2 Tương tự động học 12 2.3.3 Tương tự động lực học 13 2.3.4 Tương tự thủy động lực học 13 2.3.5 Tương tự thủy động lực học 13 2.3.6 Tương tự học 13 2.4 Các tiêu chuẩn tương tự 13 2.4.1 Sự tương tự dòng chảy lúc ảnh hưởng trọng lực chủ yếu - Tiêu chuẩn Froude 14 2.4.2 Sự tương tự dòng chảy lúc ảnh hưởng lực cản chủ yếu 14 2.4.3 Sự tương tự mơ hình dòng chảy tầng dòng chảy rối khu thành trơn thủy lực - Tiêu chuẩn Reynolds 15 2.4.4 Sự tương tự mơ hình dịng chảy rối khu sức cản bình phương 15 2.4.5 Sự tương tự mơ hình dịng chảy rối khu q độ từ thành trơn sang thành nhám thủy lực 15 2.5 Một số dẫn làm mơ hình tượng thủy lực 16 2.6 Kết luận 16 Chương - PHẦN MỀM THỦY LỰC BA CHIỀU 17 3.1 Phần mềm thủy lực ba chiều 17 3.2 Hệ thống phần mềm TELEMAC 18 3.2.1 Giới thiệu 18 3.2.2 Cấ u trúc của ̣ thố ng TELEMAC 19 3.2.2.1 Thủy lực chiều 19 3.2.2.2 Thủy lực hai chiều 19 3.2.2.3 Thủy lực ba chiều 20 3.2.2.4 Nước ngầm 20 3.2.2.5 Tải bùn cát 20 3.2.2.6 Tính tốn sóng biển 20 3.2.2.7 Bộ xử lý trước và sau tính toán (pre- post processing) 21 3.2.3 Ứng dụng 21 3.2.4 Phần mềm TELEMAC 3D 21 3.2.4.1 Phương trình với giả thiết tuân theo áp lực thủy tĩnh 21 3.2.4.2 Phương trình với giả thiết không tuân theo áp lực thủy tĩnh 22 3.2.5 Lưới tính tốn 22 3.2.5.1 Lưới hai chiều 22 3.2.5.2 Lưới ba chiều 23 3.2.5.3 Định nghĩa lưới 23 3.2.6 Mơ hình rối 23 3.2.6.1 Mơ hình Constant viscosity 25 3.2.6.2 Mơ hình Mixing length (vertical model) 25 3.2.6.3 Mơ hình Smagorinsky 25 3.2.6.4 Mơ hình k-ε 25 3.2.7 Lời giải số 25 3.2.8 Hiệu chỉnh mơ hình 26 3.2.9 Công cụ hỗ trợ 26 3.3 Phần mềm FLOW 3D 27 3.3.1 Giới thiệu 27 3.3.2 Ứng dụng 27 3.3.3 Phương trình 28 3.3.4 Lưới tính tốn 29 3.3.4.1 Lưới 29 3.3.4.2 Lưới bao gồm nhiều khối 29 3.3.4.3 Lưới phù hợp 30 3.3.5 Mơ hình rối 30 3.3.6 Lời giải số 30 3.3.7 Hiệu chỉnh mơ hình 31 3.3.8 Công cụ hỗ trợ 31 3.4 Kết luận 31 Chương – ĐÁNH GIÁ PHẦN MỀM THỦY LỰC BA CHIỀU ĐỂ MƠ PHỎNG DỊNG CHẢY XUNG QUANH CƠNG TRÌNH KÈ TRÊN SÔNG WAAL, HÀ LAN 32 4.1 Tổng quan khu vực nghiên cứu 32 4.1.1 Giới thiệu chung……………………………………………………………… 32 4.1.2 Hệ thống sông Rhine…………………………………………….…………….33 4.1.2.1 Tổng quan chung 33 4.1.2.2 Chỉnh trị sông 34 4.1.3 Sông Waal………………………………………………………….……………37 4.1.3.1 Tổng quan chung 37 4.1.3.2 Chỉnh trị sông 38 4.1.3.3 Các đặc trưng sơng Waal 39 4.2 Áp dụng phần mềm thủy lực ba chiều 42 4.2.1.Thiết lập liệu thí nghiệm………………………………………………… 42 4.2.1.1 Hình học 42 4.2.1.2 Điều kiện biên 43 4.2.2 Thiết lập liệu số ……………………………………… … ……….43 4.2.3 Kết tính tốn thảo luận…………………………………….…… ….45 4.2.3.1 Kết vận tốc mặt 45 4.2.3.2 Kết vận tốc mặt cắt ngang 46 4.2.3.3 Kết cường độ rối 49 4.2.3.4 Kết ứng suất tiếp 50 4.2.3.5 Thời gian mô công cụ hỗ trợ 53 4.2.4 Ảnh hưởng sơ đồ đối lưu mơ hình rối……………………………………53 4.2.4.1 Kết vận tốc mặt 54 4.2.4.2 Kết vận tốc mặt cắt ngang 55 4.2.5 Xác định vết dòng chảy……………………………………………………… 58 KẾT LUẬN VÀ KIẾN NGHỊ 60 KẾT LUẬN 60 KIẾN NGHỊ 61 DANH MỤC CƠNG TRÌNH KHOA HỌC ĐÃ CƠNG BỐ 62 TÀI LIỆU THAM KHẢO 623 Danh mục ký hiệu, chữ viết tắt 1D One - Dimensional 2D Two - Dimensional 3D Three - Dimensional FD Finite Difference FV Finite Volume FE Finite Element VOF Volume Of Fluid RANS Reynold Avaraged Navier Stokes RNG Renormalized Group LES Large Eddy Simulation MPI Message Passing Interface CFD Computational Fluid Dynamics FAVOR Fractional Region - Volume Obstacle Representation CGSTAB Conjugate Gradient Stabilized Method SOR Successive Over Relaxation Sadi Special Alternating Direction GMRES Generalized Minimal Residual Method GCG Generalized Conjugate Gradient UNESCO United Nations Educational, Scientific and Cultural Organization RFR Room For the River WFD Water Framework Directive DVR Sustainable Fairway Rhine Danh mục bảng Số hiệu bảng Tên bảng Trang 1.1 Các tỷ lệ mơ hình 16 3.1 Các phần mềm thủy lực 3D 19 4.1 Quá trình chỉnh trị sơng Rhine, đặc biệt sơng Waal 36 4.2 Điều kiện biên thủy lực 44 4.3 Hệ số nhám Manning 44 4.4 Thiết lập thông số phần mềm 45 4.5 Phân tích sai số 48 4.6 Phân tích sai số 52 4.7 Trường hợp thiết lập 54 4.8 Phân tích sai số trường hợp chảy khơng ngập 57 4.9 Phân tích sai số trường hợp chảy ngập phần 57 4.10 Phân tích sai số trường hợp chảy ngập hồn tồn 57 Danh mục hình vẽ, đồ thị Số hiệu Tên hình vẽ hình vẽ Trang 1.1 Mặt cắt ngang kênh sông 1.2 Sơ đồ lưới tính tốn 3.1 Hệ thống phần mềm TELEMAC 20 3.2 Xây dựng lưới ba chiều 24 3.3 Tạo lưới TELEMAC 24 3.4 Minh họa lưới FLOW 3D 30 4.1 Đất nước Hà Lan 33 4.2 Các nhánh sông Rhine 34 4.3 Các biện pháp chỉnh trị 38 4.4 Phân nhánh sông Waal với vị trí vùng lũ 38 4.5 Xây dựng kênh kiểm sốt lũ cho sơng Waal 39 4.6 Xây dựng thêm kênh để tăng khả 39 4.7 Q trình di chuyển đê xây dựng kênh thoát lũ 40 4.8 Minh họa mặt cắt ngang 41 4.9 Lưu lượng quan trắc sông Waal từ năm 1945 – 1998 41 4.10 Thay đổi đáy sông Waal từ năm 1945 – 1998 41 4.11 Sự thay đổi hình dạng sơng Waal 42 4.12 Phạm vi tính tốn 43 4.13 Thiết lập phần mềm TELEMAC 3D 44 4.14 Thiết lập phần mềm FLOW 3D 45 4.15 Kết vận tốc TELEMAC 3D 46 4.16 Kết vận tốc FLOW 3D 46 4.17 Trường vận tốc khu vực kè số 4, trường hợp chảy không ngập 47 15   (2.16) C  (2.17) Các biểu thức biểu thị tương tác lực cản,  (hoặc C, J) phụ thuộc vào trạng thái dòng chảy, cần xét riêng cho trạng thái 2.4.3 Sự tương tự mơ hình dịng chảy tầng dòng chảy rối khu thành trơn thủy lực - Tiêu chuẩn Reynolds Hệ số sức cản Đácxi  dòng chảy tầng dòng chảy rối khu thành trơn phụ thuộc vào số Reynolds (Re),   f (Re) Nên công thức biểu thị tương tự mơ hình dịng chảy tầng dịng chảy rối khu thành trơn thủy lực, tiêu chuẩn Reynolds u l (2.18) 1 v Bảng 2.1 Các tỷ lệ mơ hình Tỷ lệ mơ hình Xuất phát từ tiêu chuẩn Xuất phát từ tiêu chuẩn Tên đại lượng Froude Reynolds (Fr)n = (Fr)m với g  (Re)n = (Re)m với v  Lưu tốc 10.5 11 Lưu lượng 12.5 1 Thời gian 10.5 12 Lực  13  Số Froude biểu thị quan hệ lực quán tính trọng lực: Fr = Lực quán tính/ Trọng lực Số Reynolds biểu thị quan hệ lực quán tính lực cản: Re = Lực quán tính/ Lực cản 2.4.4 Sự tương tự mơ hình dịng chảy rối khu sức cản bình phương Trong dịng chảy rối khu sức cản bình phương hệ số Đácxi  không phụ thuộc số Reynolds mà phụ thuộc độ nhám tương đối  / R   f  (2.19) R Với dòng chảy rối khu sức cản bình phương, để có tương tự thủy động lực học, cần mơ hình tương tự hình học đủ (kể tương tự độ nhám) Phạm vi chuyển động chất lỏng gọi phạm vi mơ hình tự động tiêu chuẩn lực cản hay gọi tắt phạm vi tự động mơ hình Nếu khơng thực tương tự hình học độ nhám chọn hệ số nhám n thích hợp để có tương tự thủy lực dịng chảy rối khu sức cản bình phương 2.4.5 Sự tương tự mơ hình dịng chảy rối khu độ từ thành trơn sang thành nhám thủy lực 16 Trong trạng thái chảy này, hệ số sức cản Đácxi  phụ thuộc vào số Reynolds độ nhám tương đối  / R   (2.20)   f  Re,  R  Để có tương tự thủy động lực học, mơ hình dòng chảy khu độ từ thành trơn sang thành nhám cần có tương tự hình học, kể độ nhám, đồng thời phải thỏa mãn tiêu chuẩn Reynolds Ngồi cịn nhiều tiêu chuẩn tương tự khác thứ yếu, dùng cho trường hợp nghiên cứu cụ thể Ví dụ nghiên cứu tượng có liên quan đến khí thực cần đến tiêu chuẩn Ơle (Eu), nghiên cứu tượng có liên quan đến lực đàn hồi nước va, cần xét đến tiêu chuẩn Côsi (Ca), nghiên cứu tượng có liên quan tới lực căng mặt cần xét đến tiêu chuẩn Vêbe (We), 2.5 Một số dẫn làm mơ hình tượng thủy lực Sử dụng mơ hình phải ý đến việc chọn tỷ lệ mơ hình để bảo đảm giữ điều kiện tương tự tượng thủy lực suy từ định luật tương tự + ) Nếu ngun hình có dịng chảy rối, mơ hình dịng chảy phải chảy rối; + ) Dịng chảy mơ hình phải trạng thái với dịng chảy ngun hình Nếu ngun hình dịng chảy êm (Fr < 1) mơ hình phải chảy êm; + ) Nếu ngun hình dịng chảy xiết (Fr > 1) mơ hình phải chảy xiết Đặc biệt với mơ hình biến dạng phải ý giữ điều kiện Sử dụng mơ hình theo tiêu chuẩn tương tự lực cản điều kiện ln đảm bảo; + ) Phải đảm bảo tương tự độ nhám thành Trên mơ hình muốn đạt tượng có liên quan tới lớp biên, muốn cho lưu tốc dòng chảy mơ hình phân bố ngun hình mơ hình phải tương tự hình học độ nhám Trong nhiều trường hợp, điều khơng thể thực Lúc nghiên cứu yếu tố thủy lực trung bình (vận tốc trung bình mặt cắt, lưu lượng ); + ) …… 2.6 Kết luận Nghiên cứu mơ hình kết hợp với lý luận thực tế giúp nắm vững thêm tượng thủy lực Qua nghiên cứu tượng mơ hình có điều kiện để kiểm tra hiệu chỉnh lại công thức lý thuyết, xác lập quan hệ thực nghiệm yếu tố riêng biệt tượng thủy lực Ưu điểm mơ hình tạo nên tượng không thấy quan sát đơn giản Ngồi ra, thấy trước tượng cơng trình cần thiết kế, kiểm tra lại giả thiết tính tốn, biện pháp cơng trình đưa lúc thiết kế 17 Chương - PHẦN MỀM THỦY LỰC BA CHIỀU Với kết nghiên cứu tổng quan chương chương Mơ hình số thể thuận lợi ưu điểm nghiên cứu dòng chảy ba chiều xung quanh cơng trình Trong chương giới thiệu tổng quan phần mềm thủy lực ba chiều phổ biến nay, tập trung vào hai phần mềm Phần mềm mã nguồn mở TELEMAC 3D giải theo phương pháp phần tử hữu hạn phần mềm thương mại FLOW 3D giải theo phương pháp thể tích hữu hạn 3.1 Phần mềm thủy lực ba chiều Với tính chất phức tạp q trình thủy động lực, ngồi cơng thức kinh nghiệm, nhà khoa học tập trung vào hai phương pháp dùng mô hình vật lý mơ hình tốn số Cả hai phương pháp đóng góp đáng kể cho hiểu biết đặc trưng dòng chảy, tương tác dịng chảy hình thái khu vực kè thiết kế cơng trình chỉnh trị Tuy nhiên, thí nghiệm mơ hình hóa có tính đặc thù riêng vừa lợi bất lợi Thí nghiệm thực thơng qua mơ hình vật lý sử dụng để áp dụng cho việc thiết kế cơng trình xây dựng lớn [9] Hiện nay, phương pháp chứng minh hiệu nhánh sơng có hình dạng đơn giản, nơi chế độ thủy động lực học không phức tạp Trong trường hợp vị trí phức tạp nút giao, thay đổi lớn địa hình khu vực nơng, phương pháp dường khơng thích hợp để mơ đặc tính dịng chảy Hơn nữa, chi phí mơ hình vật lý hạn chế lớn, đặc biệt với cơng trình vừa nhỏ [10], [11] Ngược lại, với phát triển tốn học hệ thống máy tính, mơ hình số xem công cụ thuận lợi, hiệu suất cao, linh hoạt chi phí thấp cho việc phân tích đặc điểm thủy động lực học khu vực kè [12] Trong lĩnh vực Thủy lực mơi trường, mơ hình chiều (1D) dựa hai hình thức tích phân vi phân phương trình St Venant sử dụng thường xuyên Mặc dù phổ biến hữu ích mơ hình 1D khơng thể bỏ qua, thực tế mơ hình 1D gặp phải đơn giản hóa giả định bị số hạn chế cố hữu Mơ hình 1D cung cấp thơng tin mặt cắt trung bình mà không đưa thông tin thay đổi đặc tính dịng chảy theo phương ngang phương thẳng đứng Cải tiến xác dựa mơ hình số sử dụng mơ hình hai chiều (2D) Mơ hình 2D dựa hình thức trung bình theo chiều sâu phương trình ba chiều Navier-Stokes Dựa phương pháp sai phân hữu hạn, khối lượng hữu hạn phần tử hữu hạn cho lời giải số Mặc dù mơ hình 2D áp dụng thành công cho nhiều vấn đề kỹ thuật sơng, mơ hình trung bình theo chiều sâu không cung cấp hiểu biết chi tiết theo phương thẳng đứng Ví dụ trường hợp đoạn cong sông suối tự nhiên, vị trí có xây dựng cơng trình, thay đổi lớn địa hình xuất 18 dịng chảy ba chiều (3D), dịng xốy, thứ cấp Các dịng chảy gây thay đổi hình thái mặt cắt ngang kênh sơng có ý nghĩa quan trọng thay đổi hình dạng bờ sơng Và khơng thể mơ với mơ hình chiều, hai chiều Để nắm bắt tính chất vật lý dịng chảy trường hợp sử dụng mơ hình 3D cần thiết Những mơ hình dựa lời giải số phương trình Navier-Stokes có khơng có giả định áp lực phân bố theo quy luật thủy tĩnh thường kết hợp với xấp xỉ Boussinesq Các phần mềm 3D sử dụng phổ biến MIKE 3, OpenFOAM, Delft 3D, FLUENT, TELEMAC 3D, FLOW 3D, … Mặc dù dựa phương trình chủ đạo, mơ hình nói có khác biệt lớn kỹ thuật sử dụng phương pháp giải phương trình Navier-Stokes Vì số lượng mơ hình ba chiều tăng, nên cần phải có nghiên cứu để để làm rõ ưu nhược điểm phạm vi áp dụng phần mềm [36] Bảng 3.1 Các phần mềm thủy lực 3D STT Tên phần mềm MIKE FLOW 3D OpenFOAM Delft3D FLUENT TELEMAC 3D Cơ quan, nhà sản xuất Viện Thủy lực (DHI), Đan Mạch Flow Science, Inc, Mỹ OpenCFD Ltd, Mỹ Deltares, Hà Lan Phương pháp giải Sai phân hạn hữu Lưới Lăng trụ Thể tích hữu Khối hạn Thể tích hữu Khối hạn Sai phân hữu hạn, thể tích hữu Lăng trụ hạn Fluent Europe Lăng Thể tích hữu hạn Lt, Anh trụ,… Điện lực Pháp (EDF) Phòng nghiên cứu Phần tử hữu hạn Lăng trụ Thủy lực quốc gia (LNH), Pháp Chi phí Thương mại Thương mại Mã nguồn mở/ C++ Mã nguồn mở/ Fortran, C/C++ Mã nguồn mở Mã nguồn mở/ Fortran Trong nghiên cứu, tác giả sử dụng hai phần mềm đặc trưng sử dụng phổ biến thủy lực sơng TELEMAC 3D, FLOW 3D để mơ dịng chảy qua cơng trình kè đoạn sơng Waal – Hà Lan so sánh với thí nghiệm 3.2 Hệ thống phần mềm TELEMAC 3.2.1 Giới thiệu Hê ̣ thố ng TELEMAC đươ ̣c phát triể n bởi Điện lực Pháp (EDF – Electricite de France) Phòng nghiên cứu Thủy lực quốc gia (LNH – Laboratoire National d’Hydraulic) Hiện sử dụng miễn phí tải trang web www.opentelemac.org Hệ thống TELEMAC sử dụng lý thuyết phần tử hữu hạn Gần 19 phát triển thêm lý thuyết thể tích hữu hạn để phân tích cho toán hai chiều vận chuyển bùn cát khơng dính Tất thành phần cấu trúc biên soạn sử dụng ngơn ngữ lập trình FORTRAN [37], [38], [60] Hệ thống TELEMAC có công cu ̣ chuẩ n bi ̣và xử lý số liê ̣u trước và sau tính toán (pre-and post-processing) đă ̣c biê ̣t hiê ̣u quả, ta ̣o giao diê ̣n thuâ ̣n tiê ̣n và dễ dàng cho người dùng Hầ u hế t các chương trình xử lý số liê ̣u đề u đươ ̣c xây dựng nên từ các thư viê ̣n ilog/Views vì thế có thể cung cấ p cho người sử du ̣ng mô ̣t số lươ ̣ng lớn các thông tin cần thiết 3.2.2 Cấ u trúc của ̣ thố ng TELEMAC Dựa vào bước mục đích tính tốn, TELEMAC có nhiề u mơ đun khác (hình 3.1) Hình 3.1 Hệ thống phần mềm TELEMAC [60] 3.2.2.1 Thủy lực chiều TELEMAC MASCARET mơ dịng chảy chiều Dựa hệ phương trình Saint-Venant, mơ đun khác mô tượng khác khu vực rộng lớn cho hình dạng phức tạp, vị trí nối phân nhánh mạng lưới sơng, dịng chảy trạng thái phân giới phân giới trên, dịng chảy ổn định hay khơng ổn định 3.2.2.2 Thủy lực hai chiều TELEMAC 2D mô hình tính dòng chảy chiề u ngang, giải ̣ phương trình Saint-Venant (bao gồ m mô phỏng hiê ̣n tươ ̣ng truyề n các chấ t hòa tan) TELEMAC-2D có hai phiên bản dùng hai phương pháp tính khác biê ̣t 20 Phầ n tử hữu ̣n: Trong phiên bản này, ̣ phương trin ̀ h Saint-Venant viế t dưới da ̣ng không bảo toàn sẽ đươ ̣c giải bằ ng phương pháp chiế u (Projection Method) sử dụng sơ đồ ẩ n Phương pháp này có tin ̣ cao với tố c đô ̣ tin ́ h ở n đinh ́ h nhanh Thể tích hữu ̣n: Hê ̣ phương trình Saint-Venant viế t dưới da ̣ng bảo toàn sẽ đươ ̣c giải bằ ng phương pháp Godunov (xấ p xỉ bấ t biế n Riemann) sử dụng sơ đồ hiê ̣n Phương pháp này cho phép tính toán sóng gián đoa ̣n cả điạ hình phức ta ̣p (bài toán vỡ đâ ̣p) với tố c đô ̣ tin ́ h nhanh 3.2.2.3 Thủy lực ba chiều TELEMAC 3D mô mô hin ̀ h tin ́ h dòng chảy chiề u, giải ̣ phương trin ̀ h Navier-Stokes (bao gồ m mô phỏng hiê ̣n tươ ̣ng truyề n các chấ t hòa tan có hoă ̣c không tham gia phản ứng hóa ho ̣c) TELEMAC 3D sử dụng để thực loạt mô số Giải phương trình Navier-Stokes 3D qua bước thuật toán phân đoạn, sử dụng xấp xỉ Boussinesq (J M Hervouet, 2007) Một ưu điểm thuật tốn bước phân đoạn, người sử dụng lựa chọn lời giải số phù hợp với số hạng khác phương trình Navier-Stokes 3.2.2.4 Nước ngầm ESTEL 2D: Mơ hình số dịng chảy hai chiều vận chuyển chất tầng ngầm ESTEL 2D có ba mơ đun: Mơ đun dịng chảy giải phương trình Richards vùng bão hịa khơng bão hịa phương pháp phần tử hữu hạn; Mơ đun theo dõi di chuyển hạt giải phương trình khếch tán – đối lưu phương pháp MonteCarlo mô đun vận chuyển giải phương trình khếch tán – đối lưu phương pháp phần tử hữu hạn ESTEL 3D: Giải phương trình Richards phương pháp phần tử hữu hạn dựa lưới phi cấu trúc ba chiều chuyển động hạt theo phương trình Fokker-Plank ESTEL 3D kết hợp với mơ hình TELEMAC 2D để liên kết dòng chảy ngầm với dòng chảy mặt 3.2.2.5 Tải bùn cát SISYPHE: Giải bài toán tải bùn cát và biế n hình lòng dẫn hai chiề u Mô hình này phân miề n tiń h toán thành hai phầ n: - Phầ n dòng chảy phía đáy: Tính tải bùn cát lơ lửng hay sát đáy (dính và không diń h) với nhiề u thành phần ̣t khác - Phầ n dưới đáy giả đinh: ̣ Giải bài toán bồi lắng và cố kết của bùn cát dưới đáy Hai phầ n sẽ trao đổ i bùn cát tương tác với SEDI 3D: Giải bài toán tải bùn cát ba chiề u 3.2.2.6 Tính tốn sóng biển ARTEMIS: Tính tốn sóng biể n có xét đế n các hiê ̣n tươ ̣ng vâ ̣t lý phản xa ̣, nhiễu xa ̣, khuyế ch tán của sóng biể n truyề n vào vùng nước nông trước và cảng biể n 21 TOMAWAC: Tiń h toán truyề n sóng vùng ven bờ Sử dụng phương pháp phần tử hữu hạn cách đơn giản hóa hệ phương trình theo mật độ quang phổ góc sóng 3.2.2.7 Bộ xử lý trước và sau tính toán (pre- post processing) RUBENS: Vẽ kế t quả tiń h toán (graphical post-processor) MATISS/ BlueKenue: Bô ̣ ta ̣o lưới STBTEL: Giao diê ̣n lưới POTEL 3D: Vẽ các lát cắ t 2D từ kế t quả 3D 3.2.3 Ứng dụng Những ứng dụng hệ thống phần mềm TELEMAC: - Vỡ đập; - Nghiên cứu lũ; - Ảnh hưởng cơng trình; - Phát triển quản lý cửa biển; - Chất lượng nước; - Tuần hoàn nhiệt; - Khai thác vật liệu; - Vận chuyển bùn cát; - Truyền sóng cảng vịnh; - Tái tạo dịng sơng; - ………………… 3.2.4 Phần mềm TELEMAC 3D 3.2.4.1 Phương trình với giả thiết tuân theo áp lực thủy tĩnh TELEMAC giải phương trình thủy động lực học ba chiều với giả định sau: Phương trình Navier-Stokes ba chiều với bề mặt tự thay đổi theo thời gian; Biến đổi khơng đáng kể tỷ trọng bảo tồn phương trình khối lượng (chất lỏng khơng nén được); Đơn giản hóa vận tốc W theo phương đứng, bỏ qua số hạng khuếch tán, đối lưu số hạng khác; Giả thuyết áp lực phân bố tuân theo quy luật thủy tĩnh (giả thuyết kết áp lực độ sâu tổng áp suất khơng khí bề mặt chất lỏng cộng với trọng lượng chất lỏng phía trên); Xấp xỉ Boussinesq cho động lượng (các biến mật độ không đưa vào số hạng trọng lực) Theo giả thiết trên, phương trình ba chiều [39] Phương trình liên tục: 22 u v w   0 x y z Phương trình động lượng (Navier-Stokes): Z u u u u u v w   g s    u   Fx t x y z x Z v v v v  u  v  w   g s    v   Fy t x y z y (3.1) (3.2) Áp suất xem phân bố thủy tĩnh:  P  g Z  z  g z   dz (3.3) Trong đó: u, v w thành phần vector vận tốc u; Zs độ cao bề mặt thoáng; z độ cao điểm dòng chảy; v hệ số nhớt;  mật độ nước Fx , Fy số hạng nguồn 3.2.4.2 Phương trình với giả thiết không tuân theo áp lực thủy tĩnh Bằng việc đơn giản hóa vận tốc theo phương đứng W, bỏ qua khuếch tán, đối lưu số hạng khác Do đó, áp lực điểm liên quan đến trọng lượng cột nước nằm phía áp lực khí bề mặt Lúc TELEMAC 3D giải phương trình vận tốc theo phương đứng W tương tự phương trình U V Áp suất xem phân bố thủy tĩnh bổ sung thêm số hạng trọng lực Phương trình liên tục: u v w (3.4)   0 x y z Phương trình động lượng (Navier-Stokes): u P  u  gradu   div  vt gradu   FMx t  x v P  u  gradv   div  vt gradv   FMy t  y w P  u  gradw   div  vt gradw   FMz t  z (3.5) Trong đó: u, v w thành phần vector vận tốc u; Z bề mặt tự do; z độ cao điểm dòng chảy; vt hệ số nhớt;  mật độ nước FMx FMy số hạng nguồn 3.2.5 Lưới tính tốn 3.2.5.1 Lưới hai chiều Lưới hai chiều xây dựng phần tử tam giác với ba nút, phần mềm tương thích với hệ thống TELEMAC (Matisse, BLUEKENUE, ) Ngồi ra, 23 STBTEL kiểm tra định hướng xác số thứ tự lưới, liên kết linh hoạt đánh số địa phương việc đánh số nút lưới tổng thể 3.2.5.2 Lưới ba chiều Lưới 3D xây dựng hình lăng trụ gồm nút, tự động phát triển theo hai giai đoạn q trình mơ TELEMAC 3D Đầu tiên, lưới 2D xây dựng để bao phủ vùng tính tốn Sau lưới 3D xây dựng lại theo phương đứng tùy thuộc vào số lớp chọn [40] (hình 3.2) Hình 3.2 Xây dựng lưới ba chiều [40] 3.2.5.3 Định nghĩa lưới Thông qua từ khóa MESH TRANSFORMATION để định nghĩa chia lưới theo chiều dọc (hình 3.3) Bằng tương ứng với xác định người sử dụng chương trình Bằng chia theo chiều cao Bằng có linh hoạt chia lưới ứng với phân bố gần bề mặt đáy (phù hợp ứng dụng lập trình) Hình 3.3 Tạo lưới TELEMAC [40] 3.2.6 Mơ hình rối Osborne Reynolds năm 1883 chứng minh thực nghiệm tồn hai loại chế độ chảy khác chế độ chảy tầng chế độ chảy rối Các lực tác dụng lên hạt chất lỏng dịng chảy, có xu hướng giữ chúng hướng 24 lực khác xu hướng đẩy theo hướng khác Đối với hạt chất lỏng chuyển động, lực có xu hướng đẩy hạt theo hướng lực quán tính Ngược lại, chất lỏng xung quanh hạt tạo lực nhớt để hướng lớp dòng chảy theo hướng Trong dòng chảy tầng, ứng suất nhớt chiếm ưu nên loại bỏ ảnh hưởng độ lệch dòng chảy chất lỏng Mặt khác, dòng chảy rối, ứng suất nhớt vượt qua lực qn tính dịng chảy bị xáo trộn Vì vậy, tiêu chí mơ tả dịng chảy tầng rối tỷ lệ lực quán tính ứng suất nhớt tác động lên hạt chất lỏng Tỷ lệ gọi "Số Reynolds" ( R  UL /  ) [41] Dịng chảy rối đặc trưng tính chất [42]: - Dịng chảy rối khơng ổn định, đặc trưng dẫn đến trạng thái nhiễu loạn; - Dòng chảy rối dòng chảy ba chiều; - Dịng rối có cấu trúc xốy khác trải dài phạm vi rộng theo tỷ lệ chiều dài; - Dòng chảy rối giá trị quan trọng xốy; - Rối làm khuếch tán dịng chảy Khi chất lỏng với nồng độ khác đưa vào dịng chảy, q trình khuếch tán rối diễn làm thay đổi nồng độ; - Dòng rối gây tiêu hao lượng cho dòng chảy Khi động lượng hạt chất lỏng khác tiếp xúc với nhau, làm giảm động dòng chảy, giảm gradient vận tốc Năng lượng chuyển đổi sang lượng bên chất lỏng nhiệt; - Dòng chảy rối khác kích thước, độ lớn đặc tính theo thời gian Khi số Reynolds lớn, xuất trạng thái khơng ổn định dịng chảy hình thành xốy nước với kích thước khác Lý tưởng nhất, để mơ phỏng, với phương trình bảo tồn khối lượng động lượng, đầy đủ biến rối; điều thực độ phân giải lưới đủ nhỏ để bắt chi tiết Tuy nhiên, điều tương đối khó khăn hạn chế nhớ máy tính thời gian xử lý Do hạn chế không gian tỷ lệ theo chiều ngang chiều dọc mơ hình hóa Ví dụ, km biển tạo kích thước xốy theo chiều ngang, theo chiều dọc, kích thước xoáy hạn chế độ sâu nước phân chia theo lớp Vì việc tách biến rối theo phương ngang dọc theo tỷ lệ không liên quan đến trạng thái động lực áp dụng TELEMAC 3D Điều liên quan đến việc xác định độ nhớt theo phương ngang, độ nhớt thẳng đứng thay giá trị độ nhớt Từ khóa mơ hình «HORIZONTAL TURBULENCE MODEL, VERTICAL TURBULENCE MODEL” 25 3.2.6.1 Mơ hình Constant viscosity Ảnh hưởng độ nhớt khuếch tán Độ nhớt số theo phương ngang đứng tồn miền tính tốn Sử dụng từ khóa “COEFFICIENT FOR HORIZONTAL DIFFUSION OF VELOCITIES” “COEFFICIENT FOR VERTICAL DIFFUSION OF VELOCITIES” Giá trị hệ số có ảnh hưởng định kích thước hình dạng vịng xốy xốy Một giá trị nhỏ có xu hướng tiêu tan xốy nước quy mơ nhỏ, giá trị lớn có xu hướng tiêu tan xốy có quy mơ lớn 3.2.6.2 Mơ hình Mixing length (vertical model) Mơ hình rối dọc theo phương đứng với từ khóa VERTICAL TURBULENCE MODEL: = "MIXING LENGTH") Các khuyếch tán theo chiều dọc vận tốc tự động tính TELEMAC 3D tổng hợp việc chọn mơ hình rối khơng xét đến ảnh hưởng số hạng mật độ Mơ hình rối dọc theo phương đứng thể độ nhớt rối (số hạng khếch tán) hàm gradient vận tốc trung bình chiều dài xáo trộn (lý thuyết Prandtl) Người dùng phải nhập tùy chọn cho mơ hình rối dọc theo phương đứng: 1: Mơ hình chuẩn Prandtl (mặc định) ; 3: Mơ hình Nezu Nakagawa; 5: Mơ hình Quetin Phù hợp với ảnh hưởng lực kéo gió Trong thời tiết có gió, lớp cách bề mặt hình thành có độ nhớt giảm; 6: Mơ hình Tsanis Phù hợp với ảnh hưởng lực kéo gió 3.2.6.3 Mơ hình Smagorinsky Mơ hình smagorinsky đề nghị, đặc biệt trường hợp có xuất dịng chảy phi tuyến cao 3.2.6.4 Mơ hình k-ε Mơ hình k–ε đưa vào TELEMAC theo đề nghị Rodi Launder để giải phương trình rối Với từ khóa "HORIZONTAL TURBULENCE MODEL" "VERTICAL TURBULENCE MODEL" Mơ hình k–ε bao gồm phương trình giải phương trình cân cho k (năng lượng rối) ε (hệ số tiêu tán lượng rối) Áp dụng mơ hình thường địi hỏi phải sử dụng lưới 2D mịn so với mơ hình độ nhớt số nhiều thời gian tính tốn 3.2.7 Lời giải số Tất module mô viết ngôn ngữ FORTRAN Khi sử dụng mô đun mô từ hệ thống TELEMAC, người dùng có chương trình cụ thể cho phiên với mã code tiêu chuẩn Những chương trình viết để người dùng dễ dàng lập trình sửa đổi 26 TELEMAC 3D có nhiều kỹ thuật sơ đồ số để giải toán tùy theo độ xác khác [40] Để giải hệ phương trình, TELEMAC đưa bảy phương pháp khác nhau: - Gradient liên hợp (Conjugate gradient); - Số dư liên hợp (Conjugate residual); - Gradient liên hợp dựa phương trình chuẩn (Conjugate gradient on a normal equation); - Cực tiểu sai số (Minimum error); - Gradient liên hợp bình phương (Squared conjugate gradient); - Phương pháp ổn định Gradient liên hợp (Conjugate gradient stabilized method - CGSTAB); - Phương pháp số dư cực tiểu tổng quát (Generalized minimal residual method GMRES) 3.2.8 Hiệu chỉnh mơ hình Q trình hiệu chỉnh thực thơng qua mơ hình rối độ nhám bề mặt Cấu trúc phần mềm TELEMAC 3D cung cấp file để hiệu chỉnh q trình tính tốn, chương trình đặt tên VALIDA điền người sử dụng phù hợp với trường hơ ̣p Xác nhận thực tập tin tham chiếu file kết từ việc tính tốn tương tự thực tài liệu tham chiếu, với từ khóa REFERENCE FILE, liên quan đến lời giải phân tích mà sau lập trình hồn tồn người sử dụng 3.2.9 Công cụ hỗ trợ Khi mô cần phải sử dụng nhiều tài ngun máy tính, khuyến khích để chạy tính toán máy đa xử lý, cụm máy chủ TELEMAC 3D có sẵn phiên mô song song để tận dụng lợi tính tốn máy tính Phiên TELEMAC 3D chạy song song sử dụng thư viện MPI cài đặt trước Người dùng ban đầu quy định cụ thể số lượng vi xử lý sử dụng từ khóa PARALLEL PROCESSORS Đó từ khóa kiểu số nguyên giả định giá trị sau: 0: Thực phiên TELEMAC 3D thông thường (mặc định); 1: Thực phiên TELEMAC 3D song song xử lý; 2: Thực phiên TELEMAC 3D song song cách sử dụng số quy định vi xử lý Với kỹ thuật TELEMAC cho thời gian tính nhanh Nên thuận lợi việc lập trình, cách can thiệp vào mã nguồn; Lựa chọn lời giải phù hợp cho số hạng 27 3.3 Phần mềm FLOW 3D 3.3.1 Giới thiệu FLOW 3D® phần mềm tính toán động lực học chất lỏng (CFD) phát triển kỹ sư Flow Science, Inc, USA Sử dụng khả máy tính để giải phương trình chuyển động chất lỏng bao gồm phương trình động lượng, phương trình liên tục tượng vật lý dịng chảy Phần mềm FLOW 3D mơ hình thích hợp để giải vấn đề động lực học chất phức tạp có khả mơ hình hóa loạt dịng chảy Được sử dụng rộng rãi cho mơ hình 3D khơng ổn định độ cao mặt thống hình dạng phức tạp Hai kỹ thuật số sử dụng là: 1) Thể tích hữu hạn (VOF) phương pháp sử dụng để mơ hình hóa dịng chảy với độ cao - mặt thống 2) phương pháp vùng thể tích phân đoạn (Fractional Region - Volume Obstacle Representation) (FAVOR) Trong phần mềm này, phương trình liên quan đến chất lỏng giải cách sử dụng phép xấp xỉ sai phân hữu hạn (hoặc thể tích) 3.3.2 Ứng dụng Những ứng dụng phần mềm FLOW 3D: - Khơng gian vũ trụ: Cung cấp thông tin ổn định nhiên liệu/ tiếp nhận, điều chỉnh nhiệt độ đơng lạnh, tạo bọt phân phối điện tích thơng qua sử dụng bề mặt chất lỏng/ khí xác, giải pháp tản nhiệt - Tự động: Cung cấp giải pháp mô với nhiều thách thức ngành công nghiệp ô tô, bao gồm động lực dòng chảy, chuyển lượng nhiệt chất lỏng chất rắn, thay đổi pha, độ chuyển động tự chất rắn - Công nghệ sinh học: Áp dụng động lực học chất lỏng đến lĩnh vực công nghệ sinh học thực hành tương đối có tiềm lớn để cải thiện nhiều thiết bị y tế sử dụng thực FLOW 3D công cụ mô mạnh mẽ với loạt khả tích hợp, bao gồm bề mặt thống/ dịng kín, thay đổi mật độ, thay đổi pha, di chuyển đối tượng phân tích ứng suất học nhiệt - Vùng ven biển hàng hải: Các ứng dụng ven biển hàng hải bao gồm thiết kế tàu, tác động sóng gió Đối với ứng dụng ven biển, FLOW 3D dự đoán chi tiết bão sóng thần tác động lên kết cấu ven biển phân tích Tải FULL (80 trang): https://bit.ly/3vitMTt thiệt hại gây Dự phòng: fb.com/TaiHo123doc.net - Vật liệu phủ lớp bề mặt: Tối ưu hóa q trình sơn gặp khó khăn để mơ hình quy mô nhỏ chất lỏng chuyển động ảnh hưởng hiệu ứng độ bám dính sức căng bề mặt FLOW 3D phân tích q trình thơng qua mơ gradient lực căng bề mặt thay đổi nhiệt độ, truyền nhiệt, bay hơi, ngưng tụ, vận chuyển chất hòa tan kiểm soát mật độ - Năng lượng: FLOW 3D cung cấp phân tích bao gồm lĩnh vực liên quan đến lượng như: Nhiên liệu hàng hóa dập dềnh biển; hiệu ứng sóng 28 dàn khoan ngồi khơi; hiệu suất tối ưu hóa cho thiết bị, thiết kế thiết bị để thu lượng từ sóng - Cơng nghệ chế tạo: Thiết kế vịi phun, khả hấp thụ mơ hình vật liệu xốp thành phần hàng hóa tiêu dùng khác - Đúc kim loa ̣i: Cung cấp công cụ để nâng cao kết cách giải vấn đề chất lượng suất với thời gian cách xem xét khái niệm khác với mơ hình trước thép nguội bị cắt khn mẫu sửa đổi: Phun làm mát, co ngót, sấy… - Dịng chảy siêu nhỏ: Mơ hình hóa cho loạt ứng dụng dòng chất lỏng siêu nhỏ: Tích hợp khí, thể lỏng, quang học chức điện tử thiết bị nhỏ cung cấp hình ảnh 3D giúp phân tích định lượng kết - Nước Môi trường: Mô vấn đề phải đối mặt với nước ngành công nghiệp môi trường, từ dự án thủy điện lớn với hệ thống xử lý nước thải thành phố nhỏ Mơ đóng vai trị quan trọng để lựa chọn kiểm tra thiết kế, giúp làm giảm phức tạp tập trung nỗ lực vào giải pháp tối ưu Những hiểu biết giá trị thu từ thử nghiệm với số lượng tùy chọn thiết kế khác tiết kiệm thời gian chi phí đáng kể Các trường hợp ứng dụng: Đập tràn, vỡ đập, thủy lực sông, trụ cầu, khu vực ảnh hưởng thủy lực, chất gây ô nhiễm, lở tuyết, lũ quét, sóng thần, thiết kế hệ thống nước 3.3.3 Phương trình Tải FULL (80 trang): https://bit.ly/3vitMTt Dự phịng: fb.com/TaiHo123doc.net Phương trình liên tục uAx     (3.6) VF   uAx   R   vAy     wAz     RDIF  RSOR t x y z x Phương trình động lượng Ay v  R u u u  p uA  vA R  w A     Gx  f x  bx  SOR  u  uw   us   x  y z x y z  xVF  x VF  A uv R v  v v v   p   uAx  vAy R  wAz    y    R   Gy  f y  by  SOR  v  vw   vs  t VF  x y z  xVF   y  VF u  t VF  R w w w  p  vAy R  wAz  Gz  f z  bz  SOR  w  ww   ws  uAx  x y z   z VF  (3.7) Trong đó: u, v, w: Vận tốc theo phương; RDIF : Số hạng khếch tán; w  t VF RSOR: Số hạng nguồn; VF: Khối chất lỏng ô; Ax, Ay, Az: Diện tích mặt cắt ngang theo mặt;  : Tỷ trọng chất lỏng; 29 Gx, Gy, Gz: Số hạng gia tốc trọng trường; fx, fy, fz: Số hạng gia tốc độ nhớt; bx, by, bz: Số hạng tổn thất dòng chảy qua môi trường xốp vách ngăn rỗng 3.3.4 Lưới tính tốn Lưới thành phần quan trọng dùng để mơ tốn kỹ thuật mà phương pháp thường sử dụng Teng, 2000 Theo Allen, 2008 lưới cấu trúc có chất lượng cao lưới phi cấu trúc khó khăn để thể xác địa hình phức tạp FLOW 3D sử dụng hệ lưới có cấu trúc định nghĩa hệ tọa độ Descartes hay trụ Sự lựa chọn hệ (trụ Descartes) phối hợp không ảnh hưởng đến số định nghĩa hình học Tương tác hình học với mắt lưới theo phương pháp FAVORTM, nơi vật thể ngăn cách đưa vào lưới cách ngăn phần thể tích lưới vùng mặt Điều cho phép định nghĩa hoàn toàn độc lập với lưới hình học, hình học sửa đổi mà không cần xác định lại lưới Do đó, cách tiếp cận FAVORTM cung cấp đơn giản nhanh Ngồi xác định nhiều khối ô lưới để tạo mắt lưới hiệu mơ hình hóa trường dịng chảy phức tạp Các khối lưới bổ sung định nghĩa hoàn toàn khối lưới, đưa thêm biên khối ô lưới (gọi khối liên kết), phần chồng lên Hình Minh họa lưới FLOW 3D [61] 3.3.4.1 Lưới Một lưới có kích thước lưới đồng theo tất hướng Có thể tạo cách cung cấp có kích thước lưới, số lượng lưới, số ô lưới theo hướng 3.3.4.2 Lưới bao gồm nhiều khối Một loại lưới có kích thước đồng khơng hiệu cho hình học phức tạp mạng lưới bao gồm khối chứa q nhiều lưới Đối với tốn phức tạp sử dụng nhiều khối lưới để tăng độ phân giải mô khu vực quan 7761725 ... BÌNH ĐÁNH GIÁ PHẦN MỀM THỦY LỰC BA CHIỀU ĐỂ MƠ PHỎNG DỊNG CHẢY XUNG QUANH CƠNG TRÌNH KÈ TRÊN SÔNG WAAL, HÀ LAN Chuyên ngành : KỸ THUẬT XÂY DỰNG CƠNG TRÌNH THỦY Mã số : 60.58.02.02 LUẬN VĂN THẠC SĨ... hỗ trợ 31 3.4 Kết luận 31 Chương – ĐÁNH GIÁ PHẦN MỀM THỦY LỰC BA CHIỀU ĐỂ MƠ PHỎNG DỊNG CHẢY XUNG QUANH CƠNG TRÌNH KÈ TRÊN SƠNG WAAL, HÀ LAN 32 4.1 Tổng quan... kiểm tra tính xác phần mềm cách so sánh kết mô với thí nghiệm Do vâ ̣y, đề tài ? ?Đánh giá phần mềm thủy lực ba chiều để mơ dịng chảy xung quanh cơng trình kè sơng Waal, Hà Lan? ?? Thơng qua phân