Xây dựng chương trình tính toán truyền chất trong mạng đường ống cấp nước

LỜI CẢM ƠN Trước tiên xin chân thành cảm ơn TS Lê Văn Dực nhiệt tình hướng dẫn định hướng để hoàn thành luận văn Để có kết bước đầu này, xin chân thành cảm ơn Thầy giáo, Cô giáo giảng dạy Cao học- Khoa Công nghệ thông tin trường Đại Học Bách Khoa truyền đạt trang bị cho kiến thức q báu thời gian học trường Xin ghi lòng tình thương yêu từ gia đình người thân động viên tạo điều kiện để tiếp bước đường khoa học Xin chân thành cảm ơn bạn đồng nghiệp trường Đại Học Công Nghệ Sài Gòn động viên giúp đỡ thời gian học, trình hoàn thành luận văn TÓM TẮT Trong nhiều năm qua nước giới đặt mối quan tâm việc thiết kế điều hành hệ thống phân phối nước nhằm đảm bảo mặt thủy lực chất lượng nước hệ thống Có nhiều mô hình tính toán với phương pháp số xuất phát từ sai phân riêng phần, phần tử hữu hạn phương pháp thể tích hữu hạn phát triển áp dụng khắp nơi Sự hiểu biết trình vận chuyển chất hệ thống phân phối cần thiết hiểu biết công cụ mô giúp theo dõi thay đổi trình truyền chất giúp nhà điều hành, nhà thiết kế hệ thống có định đắn, nhằm đem lại tiện ích bảo đảm sức khỏe cho cộng đồng Mục đích luận văn nghiên cứu mô hình truyền chất đồng thời kết hợp với phương pháp số nhằm mô trình truyền chất Với khối lượng tính toán lớn toán truyền chất, nghiên cứu cần đến kỹ thuật tính toán số nhằm tối ưu hóa vùng nhớ thời gian, mô hình định lượng không thiết kế hệ sở liệu từ môi trường GIS để mô diễn tiến loại dòng chảy Nối tiếp công trình nhà khoa học trước đề tài nghiên cứu ứng dụng đặt mang tên : “Xây Dựng Chương Trình Tính Toán Truyền Chất Trong Mạng Đường Ống Cấp Nước” Những nội dung kết nghiên cứu trình bày thành bảy chương Chương giới thiệu mục tiêu, yêu cầu nội dung nghiên cứu đề tài Chương thiết lập mô hình toán cho hệ thống Chương trình bày giải thuật cho toán thủy lực truyền chất Chương trình bày cấu trúc tổ chức liệu Chương trình bày nội dung chương trình máy tính thực đề tài, bao gồm chức tính toán, kết xuất mô Chương trình bày kết tính toán, đánh giá kết hướng phát triển đề tài MỤC LỤC Chương GIỚI THIỆU 1 Giới Thiệu Cơ Sở Lí Thuyết 2.1 Vận Chuyển Bình Lưu Trong Đường Ống: 2.2 Sự Hòa Trộn Tại Các Điểm Giao 2.3 Sự Hòa Trộn Tại Các Thùng Chứa 2.4 Các Phản Ứng Bulk Flow 2.5 Các Phản Ứng Thành Ống Moâ Hình Và Giải Thuật Cho Bài Toán Chương Trình Máy Tính Những Mục Tiêu Và Yêu Cầu Cơ Bản Chương MÔ HÌNH TÍNH TOÁN Mô Hình Hóa Chất Lượng Nước Ổn Định Mô Hình Hóa Chất Lượng Nước Động 2.1 Phương Pháp Sai Phân Hữu Hạn Euler 2.2 Phương Pháp Thể Tích Rời Rạc (DVM: Discrete Volume Method) 11 2.3 Phương Pháp Truyền Dẫn Thời Gian (TDM: Time–Driven Method) 13 2.4 Phương Pháp Truyền Dẫn Dữ Kieän ( EDM: Event Driven Method) 16 So Sánh Các Phương Pháp 18 Chương BÀI TOÁN THỦY LỰC 19 Mô Hình Tính Toán : 19 Giải thuật Todini:[187,3] 23 Giaûi thuaät Cholesky [13]: 28 Lưu Đồ Giải Thuật Tại Một Bước Thủy Lực 31 Chương BÀI TOÁN TRUYỀN CHẤT 33 Moâ Hình Tính Toán : 33 Giải thuật 37 Lưu Đồ Giải Thuật 39 Chương HIỆN THỰC CHƯƠNG TRÌNH 40 Hệ Cơ Sở Dữ Liệu Đồ Họa Địa Hình 40 Hệ Cơ Sở Dữ Liệu Tính Toán 42 2.1 Dữ Liệu Đầu Vào 42 2.2 Dữ Liệu Đầu Ra 43 2.3 Lược Đồ Quan Hệ 43 2.4 Các Ràng Buộc Toàn Veïn: 47 Hệ Cơ Sở Dữ Liệu Mô Phỏng 48 Chương MÔ PHỎNG QUÁ TRÌNH LAN TRUYỀN CHẤT 49 Tổ Chức Dữ Liệu 49 Sơ Đồ Cấu Trúc Hệ Thoáng 50 Mô Phỏng Kết Quả: 52 Báo Cáo Kết Quả 53 Chương ĐÁNH GIÁ KẾT QUẢ 56 Tiêu Chí Đánh Giá 56 Tính Sai Số 58 Độ Hội Tụ 60 Chương KẾT LUẬN VÀ HƯỚNG MỞ RỘNG 62 Phụ Lục A: Định dạng tập tin xuất/nhập 64 Phụ Lục B: Một số kết 70 Tài liệu tham khaûo .78 Hệ thống đơn vị 79 Bài tóan truyền chất Chương GIỚI THIỆU Giới Thiệu Để đẩy mạnh việc cung cấp nước cho khu dân cư, nhiều năm qua toán thiết kế, đánh dự báo trình lan truyền chất với thông số thủy lực vốn vấn đề nóng bỏng nước nói chung thành phố Hồ Chí Minh nói riêng Nhiều toán đặt cho nhà khoa học nghiên cứu có tính chất tổng hợp sâu sắc để đảm bảo áp lực nước điểm, dư lượng chất lan truyền mạng đường ống Cơ Sở Lí Thuyết Chế độ chảy dòng nước nhân tố chủ lực để truyền tải chất Vì để đánh giá dự báo trình lan truyền chất hệ thống trước hết phải việc xem xét chế độ dòng chảy hệ thống xét Trong môi trường nước chịu ảnh hưởng diễn biến dòng chảy, phân bố chất bảo toàn thụ động phức tạp biến đổi liên tục theo thời gian Để nghiên cứu trường nồng độ chất ta giả thiết có hai nguyên nhân gây biến thiên - Sự chuyển động vó mô chất lỏng gồm nước chất nước mang theo Ởû thường xảy pha trộn túy học phân tử chất lỏng Sự truyền tải vó mô gọi truyền tải đối lưu, xác định thông lượng đối lưu hàm vận tốc dòng nồng độ chất truyền tải - Sự truyền tải mức phân tử nước phân tử chất hòa tan từ nơi đến nơi khác môi trường nước gây tượng khuếch tán phân tử (molecular diggusion) làm thay đổi trường nồng độ chất Ngoài ra, yếu tố khác làm thay đổi trường nồng độ chuyển động rối dòng chảy nghuyên nhân dẫn đến khuếch tán chất rối Trong toán truyền chất mạng đường ống, pha trộn lan truyền chủ yếu xem xét dọc chiều dòng chảy Mô hình lan truyền chất chủ yếu giải dựa sở toán khối lượng gắn liền với phản ứng động học Bài tóan truyền chất Sự vận chuyển chất xuất chế : + Bình lưu (Advection) + Sự đối lưu (Convection) + Sự truyền dẫn (Conduction) 2.1 Vận Chuyển Bình Lưu Trong Đường Ống: Một chất hòa tan vận chuyển dọc theo đường ống với vận tốc trung bình với dòng mang Sự khuếch tán theo chiều dài chế quan trọng điều kiện vận hành, điều có nghóa pha trộn khối lượng phân tử nước kề cận di chuyển Cả hai phương pháp Lagrangian Eulerian biểu diễn cho nhánh ∂C i ∂C i = r (C i ) + ui ∂t ∂x [1-1] • Ci: Nồng độ (mass/volume) ống i hàm khoảng cách x thời gian t • ui: Vận tốc dòng (Length/time) ống i • r: Tốc độ phản ứng (mass/volume/time) hàm nồng độ 2.2 Sự Hòa Trộn Tại Các Điểm Giao Tại mối nối kết, dòng vào từ hai hay nhiều nhánh, trình trộn dòng trở nên hoàn toàn tức thời Như tập trung chất nước rời điểm giao đơn giản tổng trọng số dòng tập trung từ nhánh vào Cho nút k viết : Ci| x= ∑ = j∈I k Q j C j| x= L j + Qk ,ext C k ,ext ∑ j∈I k Q j + Qk ,ext § i: Liên Kết dòng rời nút k § Ik: Tập liên kết dòng vào k § Lj: Chiều dài liên kết j § Qj: Lưu lượng (Volume/time) liên kết j [1-2] Bài tóan truyền chất § Qk,ext : Lưu lượng từ bên vào mạng nút k § Ck,ext : Nồng độ dòng bên vào nút k § Ci|x=0: Nồn g độ điểm bắt đầu liên kết i § Ci|L: Nồng độ điểm kết thúc 2.3 Sự Hòa Trộn Tại Các Thùng Chứa Để tiện lợi, ta giả thuyết chất bể chứa hòa trộn hoàn toàn Đây giả thuyết hợp lí nhiều phương tiện chứa hoạt động điều kiện điền đầy tháo cạn (fill-and-draw) qui định luồng động lượng đủ để hòa trộn dòng vào Dưới điều kiện hòa trộn hoàn toàn, nồng độ chất bể pha trộn lượng nước có lượng nước vào Cùng lúc nồng độ bên thay đổi phản ứng Phương trình sau biểu diễn tượng này: ∂(VS C S ) = ∑i∈I Qi C i|x = Li − ∑ j∈O Q j C S + r (C S ) S S ∂t ∑ i∈I S ∑ [1-3] Qi C i| x = Li = C in Qin [1-4] Q j C S = C S Qout [1-5] j∈OS § VS: Thể tích thùng thời điểm t § CS: Nồng độ bên thùng § IS: Tập hợp mối kiên kết cung cấp dòng vào § OS: Tập hợp mối liên kết rút luồng § Cin: Nồng độ dòng vào § Qin: Lưu lượng vào thùng § Qout: Lưu lượng thùng Các kiểu mô hình đặc trưng cho hòa trộn bể chứa: Bài tóan truyền chất (a) Hòa trộn hoàn toàn: Giả thiết tất nước vào bể hòa trộn hoàn toàn tức thời với nước có bể Đây dạng đơn giản nhất, không đòi hỏi thêm thông số để mô tả, dường áp dụng tốt cho nhiều phương tiện hoạt động theo cách điền đầy tháo cạn (b) Mô hình hai khoang: Chia thể tích chứa có sẵn bể thành hai khoang, tất chúng giả thiết hòa trộn hoàn toàn Các ống đầu vào/đầu bể xem đặt khoang thứ Nước vào bể hòa trộn với nước khoang thứ Nếu khoang đầy gởi dòng tràn đến khoang thứ hai để hòa trộn hoàn toàn với nước chứa Khi nước khỏi bể, tồn ngăn thứ ngăn đầy, nhận lượng nước tương đương từ ngăn thứ hai Khoang thứ mô dòng vào dòng ra, khoang thứ hai thể vùng chết Thông số cung cấp tỉ số thể tích bể tổng cộng dành cho khoang thứ (c) Mô hình dòng chảy nút FIFO: xem hòa trộn nước chút suốt thời gian tồn chúng bể Mô hình thích hợp với bể ngăn hoạt động với dòng vào dòng đồng thời (d) Mô hình dòng chảy nút LIFO: Ngược với mô hình FIFO, ô nước xếp ô lên đỉnh ô khác, nước vào rời khỏi bể nằm đáy 2.4 Các Phản Ứng Bulk Flow Trong chất di chuyển vào đường ống lưu trú thùng nơi trải qua phản ứng với phần tử cột nước Tốc độ phản ứng mô tả tổng quát sau: r = kC n [1-6] •k: Hằng số phản ứng (1/T) •n: Bậc phản ứng r (C ) = k (C − C * )C nc −1 [1-7] r (C ) = k (C * − C )C nc −1 [1-8] Bài tóan truyền chất § C* : Nồng độ giới hạn § nc : Bậc phản ứng Zero Order Decay Kinetics V ∂C =1 ∂x [1-9] First Order Decay Kinetics r (C ) = k (C − C * )C nc −1 = k (C − 0)C (1−1) = kC [1-10] Second Order Reaction Kinetics r (C ) = kC [1-11] Zero Order Growth Kinetics V ∂C =1 ∂x [1-12] First Order Saturation Growth Kinetics V ∂C = k (C * − C ) ∂x [1-13] Hoaëc C = C * − (C * − C1 )e kτ [1-14] 2.5 Caùc Phản Ứng Thành Ống Trong di chuyển dọc ống, chất phân rã vận chuyển tới thành ống phản ứng với vật liệu thành ống (Corrosion, biofilm) Hiện tượng phát biểu sau: k wall = 2k w k f R (k w + k f ) (1/T) [1-15] • kw: Hằng số tốc độ phản ứng thành ống (chiều dài/ thời gian) • kf: Hệ số vận chuyển khối (chiều dài/ thời gian) Bài tóan truyền chất • R: Bán kính ống Mô Hình Và Giải Thuật Cho Bài Toán Khi xét trình lan truyền chất hệ thống mạng đường ống, nhiều tham số thủy lực cần phải xác định vị trí khác mạn g đường ống theo thời điểm khác Do phải sử dụng mô hình thủy lực để đánh giá thông số thủy lực Đối với mạng đường ống gồm nhiều nhánh có nhiều điểm giao, toán cho nhánh ghép lại với điều kiện cân nước điểm giao nhau, từ thiết lập toán thủy lực cho dòng chảy chiều Một lời giải thủy lực xác định theo thời điểm theo vị trí định trước nhánh, sử dụng toán biên lan truyền chất giải trường nồng độ theo thời điểm vị trí xác định Việc thiết lập giải toán nêu phương pháp số thích hợp thử nghiệm sử dụng nhiều giới Nội dung đề tài trình bày phương pháp sai phân riêng phần để giải toán Chương Trình Máy Tính Bài toán lan truyền chất hệ thống đường ống vốn có khối lượng tính toán lớn Giải thuật cho toán cần đến kỷ thuật tính toán nhằm tối ưu hóa vùng nhớ Để mô kết quả, côn g cụ đồ họa nhu cầu thiếu Do đó, việc tích hợp giải pháp cho toán đòi hỏi phải xây dựng chương trình máy tính minh họa trình truyền chất Chương trình máy gồm hai môđun Một môđun tính toán viết C++ Hai môđun mô kết viết Visual Basic 6.0 có tích hợp môđun tính toán môi trường GIS với đặc điểm tính toán nhanh, tiết kiệm nhớ dễ dàng sử dụng Bằng tổ chức liệu thích hợp, liệu ban đầu điều kiện biên nhập cho chương trình thông qua giao diện thân thiện hộp thoại cho phép kết nối với hệ liệu GIS tập tin văn sẵn có theo mẫu định dạng cho trước (Phụ lục A) Lời giải toán kết xuất thành tập tin văn bản, tập tin nhị phân theo mẫu định dạng, Bài tóan truyền chất 65 [PIPES] Định nghóa đường ống mạng Mỗi hàng liệu bao gồm Mã đường ống Mã nút bắt đầu Mã nút kết thúc Chiều dài Đường kính Ft (m) Inc (mm) Hệ số nhám Hệ số Trạng tổn thất thái cục [PUMPS] Định nghóa liên kết bơm mạng Mỗi hàng liệu bao gồm Mã bơm Mã nút bắt đầu Mã nút kết thúc Mã đường đặc tính bơm [VALVES] Định nghóa liên kết van mạng Mỗi hàng liệu bao gồm Mã Van Mã nút bắt đầu Mã nút kết thúc Đường kính Dạng Van Inc (mm) Thông số Hệ số tổn thất cục [CURVES] Định nghóa đường đặc tính bơm Mỗi hàng liệu bao gồm Mã đường X (flow) Y (cột áp) đặc tính [PATTERNS] Định nghóa mẫu thời gian Mỗi hàng liệu bao gồm Mã mẫu Mốc Mốc …………… Bài tóan truyền chất 66 Xét Một ví dụ: Tank 110 10 10 11 12 11 12 Source 13 Pump 111 112 113 21 22 21 22 121 23 122 31 32 31 [JUNCTIONS] ;ID Elev Demand Pattern 10 710 11 710 150 12 700 150 13 695 100 21 700 150 22 695 200 23 690 150 31 700 100 32 710 100 [RESERVOIRS] ;ID Head 800 Pattern ; [TANKS] ;ID Elevation InitLevel MinLevel 850 120 100 MaxLevel Diameter 150 50.5 MinVol Bài tóan truyền chất 67 [PIPES] ;ID Node1 Node2 Length Diameter Roughness MinorLoss Status 10 10 11 10530 18 100 Open ; 11 11 12 5280 14 100 Open ; 12 12 13 5280 10 100 Open ; 21 21 22 5280 10 100 Open ; 22 22 23 5280 12 100 Open ; 31 31 32 5280 100 Open ; 110 12 200 18 100 Open ; 111 11 21 5280 10 100 Open ; 112 12 22 5280 12 100 Open ; 113 13 23 5280 100 Open ; 121 21 31 5280 100 Open ; 122 22 32 5280 100 Open ; [PUMPS] ;ID Node1 Node2 Đường đặc tính 9 10 HEAD ; [PATTERNS] 1.0 1.2 1.4 [CURVES] 1500 250 [CONTROLS] LINK OPEN IF NODE BELOW 110 LINK CLOSED IF NODE ABOVE 140 [QUALITY] ;Node Nồng độ khởi tạo 10 0.5 11 0.5 12 0.5 13 0.5 21 0.5 22 0.5 1.6 1.4 1.2 Bài tóan truyền chất 68 23 0.5 31 0.5 32 0.5 1.0 1.0 [REACTIONS] Order Bulk Order Wall Global Bulk -.5 Global Wall -1 Limiting Potential 0.0 [TIMES] Duration 24:00 Hydraulic Timestep 1:00 Quality Timestep 0:05 Pattern Timestep 2:00 Pattern Start 0:00 Report Timestep 1:00 Report Start 0:00 [OPTIONS] Units GPM Headloss H-W Specific Gravity 1.0 Viscosity 1.0 Trials 40 Accuracy 0.001 Pattern Emitter Exponent 0.5 Quality Chlorine mg/L Diffusivity 1.0 Tolerance 0.01 [END] Baøi tóan truyền chất 69 Tập tin xuất Tập tin xuất từ môđun tính toán: - nodes.dat: Tập tin nhị phân chứa thông tin nút mốc thời gian (nút nối, bể chứa, thùng chứa) - links.dat: Tập tin nhị phân chứa thông tin liên kết mốc thời gian (đường ống, bơm, van) - Net.txt: Tập tin văn lưu trữ thông tin hệ thống mạng, trình xử lí kết toàn thành phần mốc thời gian Tập tin xuất từ môđun Visual basic: - link.dbf: Thông tin liên kết mốc thời gian - node.dbf: Thông tin nút mốc thời gian Bài tóan truyền chất 70 PHỤ LỤC B Một số kết dựa liệu toán xét chương bảy Sơ đồ dòng chảy số mốc thời gian (giờ): Xét chất suy giảm(Chlorine): Bài tóan truyền chất 71 Bài tóan truyền chất 72 Xét chất bảo toàn(Fluoride): Bài tóan truyền chất 73 Bài tóan truyền chất 74 Xét chất gia tăng (Trihalomethanes): Bài tóan truyền chất 75 Bài tóan truyền chất 76 Nồng độ nút mốc thời gian: Bài tóan truyền chất 1 Xét chất suy giảm (Chlorine(mg/L)): Thời gian Nút 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 10 0.5 1 1 1 0.98 0.96 0.94 0.92 0.9 0.88 1 1 1 1 1 0.99 11 0.5 0.44 0.85 0.85 0.85 0.85 0.85 0.85 0.72 0.6 0.63 0.61 0.6 0.59 0.74 0.75 0.86 0.86 0.85 0.85 0.85 0.85 0.85 0.85 0.83 12 0.5 0.43 0.79 0.79 0.78 0.78 0.78 0.77 0.81 0.79 0.77 0.76 0.74 0.73 0.65 0.68 0.79 0.79 0.79 0.79 0.79 0.78 0.78 0.77 0.6 13 0.5 0.43 0.36 0.31 0.27 0.51 0.51 0.51 0.53 0.56 0.58 0.56 0.55 0.54 0.53 0.5 0.47 0.38 0.44 0.5 0.5 0.51 0.51 0.51 0.51 21 0.5 0.41 0.36 0.74 0.74 0.74 0.73 0.73 0.65 0.57 0.48 0.43 0.41 0.41 0.46 0.64 0.74 0.74 0.74 0.74 0.74 0.74 0.74 0.73 0.73 22 0.5 0.43 0.37 0.32 0.28 0.54 0.53 0.56 0.61 0.66 0.65 0.63 0.62 0.61 0.57 0.43 0.42 0.44 0.5 0.54 0.54 0.54 0.54 0.53 0.55 23 0.5 0.44 0.39 0.35 0.31 0.27 0.23 0.19 0.16 0.23 0.27 0.31 0.32 0.31 0.27 0.29 0.33 0.33 0.33 0.27 0.24 0.26 0.27 0.29 0.21 31 0.5 0.4 0.32 0.28 0.52 0.52 0.51 0.52 0.51 0.49 0.43 0.37 0.31 0.29 0.28 0.28 0.34 0.45 0.52 0.52 0.52 0.52 0.52 0.51 0.52 32 0.5 0.39 0.3 0.24 0.2 0.17 0.15 0.26 0.28 0.3 0.32 0.34 0.33 0.32 0.29 0.26 0.24 0.2 0.19 0.21 0.25 0.27 0.27 0.27 0.28 1 1 1 1 1 1 1 1 1 1 1 1 0.97 0.93 0.91 0.89 0.87 0.85 0.83 0.81 0.8 0.78 0.76 0.75 0.73 0.72 0.7 0.69 0.68 0.66 0.65 0.64 0.63 0.62 0.61 0.6 Xét chất bảo toàn (Fluoride(mg/L)): Thời gian Nút 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 10 0.5 1 1 1 1 1 1 1 1 1 1 1 1 11 0.5 0.5 1 1 1 1 0.98 0.98 0.98 0.98 1 1 1 1 1 12 0.5 0.5 1 1 1 0.98 0.98 0.98 0.98 0.98 0.98 0.98 1 1 1 1 0.98 13 0.5 0.5 0.5 0.5 0.5 1 1 0.98 0.98 0.98 0.98 0.98 0.98 0.98 0.98 1 1 1 21 0.5 0.5 0.5 1 1 1 1 0.98 0.98 0.98 0.98 1 1 1 1 1 22 0.5 0.5 0.5 0.5 0.5 1 1 0.98 0.98 0.98 0.98 0.98 0.98 0.98 0.98 0.99 0.99 0.99 0.99 0.99 0.99 0.99 23 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.8 0.89 0.99 0.99 0.99 0.99 0.98 0.98 0.98 0.98 0.98 0.98 0.99 0.99 0.99 31 0.5 0.5 0.5 0.5 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.98 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 32 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.9 0.99 0.99 0.99 0.98 0.98 0.98 0.98 0.98 0.98 0.98 0.98 0.98 0.99 0.99 0.99 0.99 0.99 1 1 1 1 1 1 1 1 1 1 1 1 0.99 0.98 0.98 0.98 0.98 0.98 0.98 0.98 0.98 0.98 0.98 0.98 0.98 0.98 0.98 0.98 0.98 0.98 0.98 0.98 0.98 0.98 0.98 0.98 Bài tóan truyền chất 78 TÀI LIỆU THAM KHẢO [1] NGUYỄN HỮU CHÍ , “Cơ học chất lỏng ứng dụng”, NXB ĐH THCN , Hà Nội 1973 [2] “Computer Method of Unsteady Flow in open Channels in the last four decades”, Trainning cource on modeling, HoChiMinh City 2000 [3] Rossman, “L.A., EPANET – Users Manual”, Sept 2000 [4] Author Adam M Strafaci, “ WaterFacts” [5] Czech Republic Czech Technical University , “Vyskov Hydraulic Analysis of Water Supply Network” [6] Rossman LA, Boulos PF, Altman T, “Discrete volume-element methods for network water-quality models”, ASCE Journal of Water Resources Planning and Management 119(5):505-17, 1993 [7] Rossman LA, Clark RM, Grayman WM, “Modeling chlorine residuals in drinking-water distribution systems”, Journal of Environmental Engineering, ASCE 120(4):803-20, 1994 [8] Professor Peter Hunter Associate Professor Andrew Pullan, “Finite Element Method - Boundary Element Method” [9] Joe D Hoffma, “Numerical Methods for Engineers and Scientists”, 2001 [10] Christopher G Koutias, professor of Coastal Engineering, “Elements of computational Hydraulics”, Democritos University of Thrace, Greece [11] Ralph A Wurbs, “Modeling and Analysis of Reservoir System Operations” [12] Le Song Giang “Phương Pháp Tính Bài Toán Thủy Lực Mạng Đường n g” [13] A George and J W-H Liu “Computer Solution of Large Sparse Positive Definite Systems”, Prentice-Hall, 1981 [14] Lưu Công Đào Nguyễn tài “ Sổ Tay Tính Toán Thủy Lực” , NXB Nông Nghiệp, 1984 Bài tóan truyền chất 79 Tham số Nồng độ Nhu cầu Đường kính ống Đường kính thùng Hiệu suất Cao độ Năng lượng Đơn vị dòng Hệ số ma sát Cộ áp Chiều dài Hệ số tổn thất Minor Công suất p suất Hệ số phản ứng (Bulk) Hệ số phản ứng (Wall) Hệ số nhám Vận tốc Thể tích Hệ US Hệ SI Metric mg/L or ug/L mg/L or ug/L Đơn vị dòng Đơn vị dòng Inches millimeters feet meters Phần trăm Phần trăm feet meters kwatt – hours kwatt - hours CFS (cubic feet / sec) LPS (liters / sec) GPM (gallons / min) LPM (liters / min) MGD (million gal / day) MLD (megaliters / day) IMGD (Imperial MGD) CMH (cubic meters / hr) AFD (acre-feet / day) CMD (cubic meters / day) Không đơn vị Không đơn vị feet meters feet meters Không đơn vị Không đơn vị horsepower kwatts psi Meters H2O 1/day (1st-order) 1/day (1st-order) mass/sq-ft/day (0-order) mass/sq-m/day (0-order) ft/day (1st-order) meters/day (1st-order) -3 10 feet (Darcy-Weisbach) , mm (Darcy-Weisbach) Không đơn vị với công Không đơn vị với thức khác công thức khác ft/sec meters/sec cubic feet cubic meters ... tiếp công trình nhà khoa học trước đề tài nghiên cứu ứng dụng đặt mang tên : ? ?Xây Dựng Chương Trình Tính Toán Truyền Chất Trong Mạng Đường Ống Cấp Nước? ?? Những nội dung kết nghiên cứu trình bày... phân riêng phần để giải toán Chương Trình Máy Tính Bài toán lan truyền chất hệ thống đường ống vốn có khối lượng tính toán lớn Giải thuật cho toán cần đến kỷ thuật tính toán nhằm tối ưu hóa vùng... pháp truyền dẫn thời gian kết hợp với phương pháp số xây dựng giải thuật tương ứng cho toán - Nghiên cứu môi trường GIS xây dựng hệ liệu cho hệ thống - Xây dựng chương trình tính toán mô trình truyền