Nghiên cứu phương pháp phân rã và xây dựng phần mềm giải bài toán ô nhiễm khí quyển (Luận văn thạc sĩ)

Nghiên cứu phương pháp phân rã và xây dựng phần mềm giải bài toán ô nhiễm khí quyển (Luận văn thạc sĩ)Nghiên cứu phương pháp phân rã và xây dựng phần mềm giải bài toán ô nhiễm khí quyển (Luận văn thạc sĩ)Nghiên cứu phương pháp phân rã và xây dựng phần mềm giải bài toán ô nhiễm khí quyển (Luận văn thạc sĩ)Nghiên cứu phương pháp phân rã và xây dựng phần mềm giải bài toán ô nhiễm khí quyển (Luận văn thạc sĩ)Nghiên cứu phương pháp phân rã và xây dựng phần mềm giải bài toán ô nhiễm khí quyển (Luận văn thạc sĩ)Nghiên cứu phương pháp phân rã và xây dựng phần mềm giải bài toán ô nhiễm khí quyển (Luận văn thạc sĩ)Nghiên cứu phương pháp phân rã và xây dựng phần mềm giải bài toán ô nhiễm khí quyển (Luận văn thạc sĩ)Nghiên cứu phương pháp phân rã và xây dựng phần mềm giải bài toán ô nhiễm khí quyển (Luận văn thạc sĩ)Nghiên cứu phương pháp phân rã và xây dựng phần mềm giải bài toán ô nhiễm khí quyển (Luận văn thạc sĩ)Nghiên cứu phương pháp phân rã và xây dựng phần mềm giải bài toán ô nhiễm khí quyển (Luận văn thạc sĩ)Nghiên cứu phương pháp phân rã và xây dựng phần mềm giải bài toán ô nhiễm khí quyển (Luận văn thạc sĩ)Nghiên cứu phương pháp phân rã và xây dựng phần mềm giải bài toán ô nhiễm khí quyển (Luận văn thạc sĩ)Nghiên cứu phương pháp phân rã và xây dựng phần mềm giải bài toán ô nhiễm khí quyển (Luận văn thạc sĩ)Nghiên cứu phương pháp phân rã và xây dựng phần mềm giải bài toán ô nhiễm khí quyển (Luận văn thạc sĩ)Nghiên cứu phương pháp phân rã và xây dựng phần mềm giải bài toán ô nhiễm khí quyển (Luận văn thạc sĩ)Nghiên cứu phương pháp phân rã và xây dựng phần mềm giải bài toán ô nhiễm khí quyển (Luận văn thạc sĩ)Nghiên cứu phương pháp phân rã và xây dựng phần mềm giải bài toán ô nhiễm khí quyển (Luận văn thạc sĩ)Nghiên cứu phương pháp phân rã và xây dựng phần mềm giải bài toán ô nhiễm khí quyển (Luận văn thạc sĩ)Nghiên cứu phương pháp phân rã và xây dựng phần mềm giải bài toán ô nhiễm khí quyển (Luận văn thạc sĩ)

ĐẠI HỌC THÁI NGUYÊN TRƢỜNG ĐẠI HỌC CÔNG NGHỆ THÔNG TIN VÀ TRUYỀN THÔNG PHẠM THANH NGHỊ NGHIÊN CỨU PHƢƠNG PHÁP PHÂN RÃ VÀ XÂY DỰNG PHẦN MỀM GIẢI BÀI TỐN Ơ NHIỄM KHÍ QUYỂN LUẬN VĂN THẠC SĨ KHOA HỌC MÁY TÍNH THÁI NGUYÊN - 2020 ĐẠI HỌC THÁI NGUYÊN TRƢỜNG ĐẠI HỌC CÔNG NGHỆ THÔNG TIN VÀ TRUYỀN THÔNG PHẠM THANH NGHỊ NGHIÊN CỨU PHƢƠNG PHÁP PHÂN RÃ VÀ XÂY DỰNG PHẦN MỀM GIẢI BÀI TỐN Ơ NHIỄM KHÍ QUYỂN Chun ngành: Khoa học máy tính Mã số: 48 01 01 LUẬN VĂN THẠC SĨ KHOA HỌC MÁY TÍNH Giáo viên hƣớng dẫn: TS Nguyễn Đình Dũng THÁI NGUYÊN - 2020 i LỜI CAM ĐOAN Tôi xin cam đoan luận văn tơi thực hiện, hướng dẫn khoa học TS Nguyễn Đình Dũng, kết lý thuyết trình bày luận văn tổng hợp từ kết cơng bố có trích dẫn đầy đủ, kết chương trình thực nghiệm luận văn tơi thực hồn tồn trung thực, sai tơi hồn tồn chịu trách nhiệm Thái Ngun, tháng năm 2020 Học viên cao học Phạm Thanh Nghị Xác nhận khoa chuyên môn Xác nhận giáo viên hƣớng dẫn khoa học TS Nguyễn Đình Dũng ii LỜI CẢM ƠN Luận văn hoàn thành Trường Đại học Công nghệ Thông tin Truyền thông hướng dẫn TS Nguyễn Đình Dũng Tác giả xin bày tỏ lòng biết ơn tới thầy cô giáo Trường Đại học Công nghệ Thông tin Truyền thông, Đại học Thái nguyên, thầy cô giáo thuộc Viện Công nghệ Thông tin – Viện Hàn lâm Khoa học Công nghệ Việt Nam tạo điều kiện, giúp đỡ tác giả trình học tập làm luận văn Trường, đặc biệt tác giả xin bày tỏ lịng biết ơn tới TS Nguyễn Đình Dũng tận tình hướng dẫn cung cấp nhiều tài liệu cần thiết để tác giả hồn thành luận văn thời hạn Xin chân thành cảm ơn anh chị em học viên cao học bạn bè đồng nghiệp trao đổi, khích lệ tác giả trình học tập làm luận văn Trường Đại học Công nghệ Thông tin Truyền thông – Đại học Thái Nguyên Cuối tác giả xin gửi lời cảm ơn đến gia đình, người ln bên cạnh, động viên khuyến khích tác giả trình thực đề tài Thái Nguyên, ngày tháng Học viên cao học Phạm Thanh Nghị năm 2020 iii MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN ii MỞ ĐẦU .1 CHƢƠNG CÁC MƠ HÌNH TỐN HỌC TRONG VẤN ĐỀ MƠI TRƯỜNG5 1.1 Phương trình truyền tải vật chất khí quyển, tính nghiệm 1.2 Phương trình truyền tải dừng 10 1.3 Bài tốn truyền tải khuếch tán vật chất, tính nghiệm .15 1.4 Bài toán liên hợp cho miền ba chiều 21 1.5 Tính nghiệm tốn liên hợp 26 1.6 Kết luận chương 29 CHƢƠNG PHƯƠNG PHÁP PHÂN RÃ GIẢI BÀI TỐN KHƠNG DỪNG31 2.1 Các lược đồ sai phân xấp xỉ cấp hai cho tốn khơng dừng với toán tử phụ thuộc thời gian [4] 31 2.1.1 Bài toán 31 2.1.2 Xét toán 37 2.2 Phương pháp phân rã .39 2.2.1 Bài toán .39 2.2.2 Bài tốn khơng 40 2.3 Phương pháp phân rã nhiều thành phần .44 2.3.1 Bài toán 45 2.3.2 Bài toán không 46 2.4 Kết luận chương 49 CHƢƠNG ỨNG DỤNG PHƢƠNG PHÁP PHÂN RÃ TRONG BÀI TOÁN Ô NHIỄM KHÍ QUYỂN 50 3.1 Bài tốn nhiễm khí 50 3.2 Sai phân biến không gian .51 3.2.1 Cấp xấp xỉ toán tử sai phân 53 3.2.2 Tính khơng âm tốn tử sai phân 54 iv 3.3 Lược đồ phân rã giải tốn nhiễm khí 55 3.4 Một số kết thực nghiệm 58 3.5 Kết luận chương 60 KẾT LUẬN VÀ HƢỚNG PHÁT TRIỂN .61 TÀI LIỆU THAM KHẢO 63 MỞ ĐẦU Tính khoa học cấp thiết đề tài Thực tế cho thấy, số lượng lớn toán thực tiễn phức tạp giải nhờ cơng cụ phương trình liên hợp Chẳng hạn, toán chế lượng tử, lượng hạt nhân, trình động lực học phi tuyến vật lý, hoá học nhiều vấn đề khác Trong phạm vi luận văn này, đề cập tới vấn đề mơi trường khí hậu Sự tác động qua lại chúng vấn đề trọng tâm khoa học, ảnh hưởng trực tiếp tới sống trái đất Môi trường quanh ta, mơi trường nước (nước mặt: sơng hồ, biển, đại dương; nước ngầm: dòng chảy lòng đất), mơi trường khơng khí, mơi trường đất mơi trường sinh thái Các thành phần môi trường luôn biến đổi, chúng tác động qua lại với chuyển hoá từ trạng thái sang trạng thái khác Trong mơi trường khơng khí, khí quyển, thành phần chúng pha trộn lẫn với nhau(theo tỷ lệ đó), dịch chuyển nhờ gió khuếch tán Khí thải cơng nghiệp tác nhân lớn làm ô nhiễm không khí Để bảo vệ môi trường sống phải hiểu qui luật khách quan từ có biện pháp tích cực hữu hiệu để bảo vệ môi trường Các thực thể vật chất bị nhiễm bẩn dạng khí (khói nhà máy, lò hạt nhân, núi lửa v.v ) lan truyền khuếch tán khí quyển, tác động với (dưới ảnh hưởng nhiệt độ, độ ẩm) trở thành hợp chất phức tạp, ta gọi chung hợp chất khí Trong q trình chuyển động thành phần hợp chất khí tác động với nhau, số thành phần từ không độc hại trở thành độc hại sống sinh vật Quá trình dẫn đến tình trạng nhiễm lục địa đại dương Để giải điều ta cần phải biết trình lan truyền khuếch tán thực thể nhiễm bẩn môi trường, mà di chuyển liên tục khí chúng biến đổi từ thành phần khơng có hại thành thành phần có hại ngược lại, thường làm ô nhiễm đại dương lục địa Đó vấn đề đáng quan tâm Vì giới khơng ngừng hồn thiện, văn minh nhân loại ngày phát triển, điều cần thiết phải dự đoán xu hướng phát triển ngành công nghiệp, để kết hợp với vấn đề gìn giữ thiên nhiên, mơi trường Hơn hết, cần phải tiến hành đầu tư vốn cần thiết để không điều chỉnh tiềm sẵn có thiên nhiên bị đi, mà cịn nâng cao nó, cải thiện mơi trường Tuy nhiên điều địi hỏi lượng kinh phí lớn Song vấn đề quan trọng cần chứng minh sức mạnh tồn cầu Vấn đề chỗ thoả thuận việc cấm vũ khí hạt nhân thoả thuận việc giữ gìn hệ thống sinh thái đảm bảo sống trái đất Vì vấn đề mang tính tồn cầu Tổ chức quốc gia thống kêu gọi lập thoả hiệp việc sử dụng thiên nhiên quốc gia, quan tâm toàn nhân loại Thiên nhiên nguồn cải người điều tất người cơng nhận, họ có phương pháp hành động để giải vấn đề nêu trên, tạo điều kiện phát triển văn minh, gìn giữ làm tăng thêm phong phú thiên nhiên Ở phương diện toán học, nhiệm vụ chủ yếu để giải vấn đề xây dựng mơ hình tốn học phản ánh đắn chất tự nhiên khách quan tượng, tìm mối quan hệ biện chứng định tính, định lượng phương pháp hữu hiệu nhằm giải tốn đặt để từ định chiến lược bảo vệ chất lượng môi trường sống (xem [1]- [3], [6]-[9], [11, 12]) Nội dung đề tài này, học viên trình bày phương trình liên hợp phân tích dựa phương trình thừa nhận, điều kiện biên, điều kiện ban đầu, phương pháp giải toán để thu kết cuối mà nhờ chúng đánh giá mức độ tác động thực trạng ô nhiễm môi trường vùng lãnh thổ (xem [5, 10]) Được gợi ý thầy giáo hướng dẫn chọn đề tài: “Nghiên cứu phương pháp phân rã xây dựng phần mềm giải toán ô nhiễm khí quyển” làm luận văn tốt nghiệp Mục tiêu luận văn tìm hiểu phương pháp phân rã giải tốn khơng dừng xây dựng ứng dụng tồn nghiễm khí Đối tƣợng phạm vi nghiên cứu đề tài Nội dung luận văn đặt vấn đề nghiên cứu vấn đề bản: - Các mơ hình tốn học vấn đề mơi trường; - Phương pháp phân rã giải tốn khơng dừng; - Xây dựng ứng dụng phương pháp phân rã tốn nhiễm khí Phƣơng pháp luận nghiên cứu - Phương pháp nghiên cứu lý thuyết: Tổng hợp, nghiên cứu tài liệu mô hình tốn học vấn đề mơi trường - Phương pháp nghiên cứu thực nghiệm: Sau nghiên cứu lý thuyết, luận văn tập trung vào xây dựng chương trình giải tốn nhiễm khí mơi trường Matlab; Đánh giá kết sau thử nghiệm - Phương pháp trao đổi khoa học: Thảo luận, xemina, lấy ý kiến chuyên gia Nội dung bố cục luận văn Ngoài phần mở đầu, kết luận hướng phát triển, luận văn bố cục thành ba chương sau: Chương phân tích mơ hình tốn học khác vấn đề nhiễm mơi trường Các phương trình rút từ quy luật bảo toàn (bảo toàn khối lượng động lượng) Mỗi toán xây dựng toán liên hợp tương ứng nhờ đẳng thức tích phân Lagrange Tính nghiệm toán toán liên hợp mơ hình chứng minh cách chặt chẽ nhờ đẳng thức đối ngẫu Chương xây dựng phương pháp giải toán đặt chương Do độ phức tạp phương trình, với giả thiết điều kiện biên, giá trị ban đầu chặt chẽ người ta nhận nghiệm xác tốn Thực tế cho thấy toán đặt thường rộng hơn, phức tạp Do đó, việc tìm phương pháp giải số cho lớp toán phương pháp hữu hiệu sử dụng Tuy nhiên việc sử dụng phương pháp số gặp nhiều khó khăn số chiều lớn, miền phức tạp v.v Rất nhiều nhà học, kỹ sư đã đưa nhiều phương pháp sai phân cho tốn nhiễm môi trường Tuy nhiên thiếu tảng toán học, số phương pháp sai phân cho lời giải không phù hợp với thực tiễn Để khắc phục hạn chế ta phân rã phương trình khuyếch tán truyền tải trình vật lý: Tại bước thời gian giải phương trình truyền tải phương pháp đặc trưng, giải toán khuyếch tán phương pháp phân rã Nội dung chương trình bày chi tiết ổn định, cấp xác phương pháp phân rã tốn tiến hố tốn khơng nhất, phân tích nhược điểm phương pháp tốn tử tốn tổng hai tốn tử khơng giao hốn Để khắc phục nhược điểm chúng tơi đưa phương pháp phân rã cho toán chứng minh chi tiết tính ổn định vơ điều kiện cấp xác theo thời gian Tiếp theo trình bày cho trường hợp tổng quát toán tử toán tổng nhiều toán tử nửa xác định dương Chương xấp xỉ toán tử vi phân toán khuyếch tán đặt Chương toán tử sai phân với cấp xác hai theo biến khơng gian thoả mãn tính khơng âm Cuối đưa lược đồ phân rã theo thời gian cài đặt thuật toán kết số toán sai phân xấp xỉ nghiệm toán vi phân xây dựng Chương 50 CHƢƠNG ỨNG DỤNG PHƢƠNG PHÁP PHÂN RÃ TRONG BÀI TỐN Ơ NHIỄM KHÍ QUYỂN 3.1 Bài tốn nhiễm khí Bài tốn nhiểm khí nói chi tiết chương Tuy nhiên nhắc lại mô hình tốn để tiện theo dõi [5] Xét tốn lan truyền khí thải gây nguồn phát thải Phương trình tốn viết sau:       u v w        f t x y z z z (3.1.1) Với điều kiện ban đầu   g t  (3.1.2) Trong  ( x, y, z, t ) nồng độ khí thải mơi trường điểm ( x, y, z) thời điểm t u, v, w thành phần véc tơ vận tốc gió đề cập chương f công suất khí thải, g phân bố nồng độ ban đầu chất khí thải,  đặc trưng hệ số biến đổi khí thải q trình lan truyền,  ,  hệ số khuếch tán rối ngang thẳng đứng,  toán tử Laplace hai chiều Nếu f  Q ( x  x0 ) ( y  y0 ) ( z  z0 ) ,  hàm Dirac tốn trở thành tốn lan truyền khí thải gây nguồn điểm đặt ( x0 , y0 , z0 ) có cường độ Q Để hạn chế, ta quan sát nghiên cứu tượng toán miền G  0, X  0,Y  0, H  , với bề mặt S S    0   H Trong đó:  mặt xung quanh miền G  mặt đáy (khi z  )  H mặt đáy (khi zH ) 51 Do lớp khí (lớp tiếp giáp với mặt đất) với độ xác cao có  thể xem khơng khí chất khơng nén được, nên véc tơ vận tốc gió V thể phương trình liên tục u v w   0 x y z Ngoài ra, thành phần thẳng đứng véc tơ vận tốc gió triệt tiêu mặt đáy w  z  z  H Để đơn trị xác định ẩn hàm  cần xác định ràng buộc nồng độ khí thải  biên S miền G Từ thực tiễn lý thuyết phương trình vật lý tốn viết điều kiện biên tốn lan truyền khí thải gây nguồn phát thải sau:        z (3.1.3)    H z Trong   hệ số hấp thụ phản xạ mặt đáy Bài tốn hồn tồn thoả mãn tính nghiệm, điều chứng minh chương 3.2 Sai phân biến không gian Sử dụng phương trình liên tục trừ vế trái biểu thức (3.1.1) cho biểu thức:  u ( x  v w  )0 y z ta phương trình    u   v   w    (u  )  (v  )  (w  )        f t x x y y z z z z Viết lại phương trình tốn lan truyền khí thải gây nguồn phát thải (3.1.1) sau: 52   L  f G  (0,T ] t (3.2.1) với   u  2  ) x x x   v  2 L2  (v  ) y y y   w   L3  ( w  )     z z z z L1  (u (3.2.2) Giả thiết nghiệm toán (3.1.1) , (3.1.2) , (3.1.3) đủ trơn miền G  [0,T ] Ta phủ miền G lưới sai phân     ( xk  khx , yl  lh y , z m  mhz ), k  0, K , l  0, L, m  0, M Các toán tử vi phân Li xấp xỉ toán tử sai phân Ai (u * ) x  (u * ) x   xx (v * ) y  (v * ) y A2    yy ( w* ) z  ( w* ) z A3   ( * z ) z   A1  (3.2.3) u * , v* , w* , * hàm lưới xác định lưới  uk*  u k , vl*  u l , wm*  w m ,  m*   m Chú ý viết số theo hướng, không viết số hướng khác, viết số k , l , m tương ứng hướng biến x, y, z , ví dụ: u k thay cho u k  lm Hệ (3.2.3) viết lại sau: 53 u ( A1 ) k  k  k 1  u k  k 1  2h x v  l 1  v  l 1 ( A2 ) l  l  2h y w ( A3 ) m  l m  m 1  w m hx2  l 1  2 l   l 1 (3.2.4) h y2  m 1   2h z  k 1  2 k   k 1 m ( m 1   m )   m ( m   m 1 ) hz2   m k  1, K , l  1, L, m  0, M  Với điều kiên biên (3.1.3) phân rã sau: k  k  0, k  K  l  l  0, l  L  m   m1  m m  hz m  (3.2.5) m  M 3.2.1 Cấp xấp xỉ toán tử sai phân Bổ đề 3.2.1 Các toán tử sai phân Ai cho công thức (3.2.3) xấp xỉ cấp toán tử vi phân Li Chứng minh Ta cần chứng minh cho hai trường hợp Chứng minh toán tử sai phân sau (u * ) x  (u * ) x , ( * z ) z có xấp xỉ cấp tương ứng với toán tử vi phân sau u   u    ,  x x z z Thật vậy, khai triển công thức Taylor theo biến x xk biến z z m ta có: ( ) k 1  ( ) k  hx ( h  2  ) k  x ( ) k  O(hx3 ) x x 54 u k hx u h  2u ( ) k  x ( ) k  O(hx3 ) x x  uk  ( ) m1  ( ) m  hz (  m h  2 h  3  ) m  z ( ) m  z ( ) m  O(hx4 ) z z z hz  hz2  2   m  ( ) m  ( ) m  O(hz3 ) z z ( ) k   ( xk ), (     ) k  ( )( xk ), ( ) m  ( )( z m ), u k  u ( xk ),  m   ( z m ) x x z z Từ suy u k ( ) k 1  u k ( ) k 1  (u 2hx  m [ ( ) m1  ( ) m ]   h m   u  ) k  O(hx2 ) x x [ ( ) m  ( ) m1 ] z (    ) m + O(hz ) z z từ suy đpcm 3.2.2 Tính khơng âm toán tử sai phân Bổ đề 3.2.2 Với giả thiết đặc trưng hệ số biến đổi khí thải, hệ số khuếch tán rối ngang thẳng đứng  ,  ,  trình lan truyền, thành phần u, v, w véc tơ vận tốc gió khơng âm, tốn tử Ai khơng âm Chứng minh K Ta có ( A1 ,  )   hx ( A1 ) k  k , từ cơng thức (3.2.4) ta có k 1  (u  k 1  u  k 1 ) k  k ( k 1  2 k   k 1 ) k   k2 ( A1 , ) k      hx k 1     K 1    (u  K  K 1  u  01 )   ( k 1   k )  (12   01 )  K2 hx k 1 hx K   hx ( K2   K  K 1 ) 55 Sử dụng điều kiện biên (3.2.5) suy A1  Do vai trò x, y nên ta có A2  Sau phép biến đổi đơn giản, ta thu M 1 ( A3 , )   hz ( A3 ) m  m  ( w  M 1 M  w  1 )   M 2 m 0  hz M 1  m 0 M 1 ( m1   m )   1  hz  m2 m  m 0 Kết hợp điều kiện biên (3.2.5) suy A3  3.3 Lƣợc đồ phân rã giải tốn nhiễm khí Do tốn tử Ai khơng giao hốn chứng minh khơng âm, nên để giải tốn lan truyền khí thải (3.1.1) , (3.1.2) , (3.1.3) ta sử dụng phương pháp phân rã nhiều thành phần với thể khác Cụ thể, tốn lan truyền khí thải gây nguồn phát thải (3.1.1) , (3.1.2) , (3.1.3) giải sơ đồ sai phân sau: Đặt j  Aj   1, 2, Lược đồ sai phân dạng thứ nhất: 2 j j j (E   1j )  (E   1j ) j 1 ,   0,  K  2 1 j j j j (E    j2 )  (E    j2 ) ,   0,  L  2 j  j   0j  0j ,  Mj  (E    j3 )( j  f j )  (E    j3 ) , 2 hz (E    j3 ) j j  (E    j3 )( j  f j ) , 2 1 1 j  0 hz j j j j j j (E    j2 )  (E    j2 ) ,   0,  L  2  0 ,  M  j (E   1j ) j 1  (E   1j ) ,  0j 1  0,  Kj 1  2 56 Lược đồ sai phân dạng thứ hai:  j   j 1   j  j   j j    j3  j j  1   j j   j 4 j  j   j    j3 j    j 1 j     j2  1j j j  j j 1  2f j 0, j j  , 0  0,  j j   j 1 j 4 j   j     j j   j 1  ,   0,  K  2 4 1 j  0,  L j  0 hz 4 0 j j j j  0 ,  M  j j  0 hz j  , 0 4  0 ,  M  j  0,  L 0  ,  0j 1  0,  Kj 1  Lược đồ sai phân dạng thứ ba:  j   j 1 j   1   j  j     j  j  j   j j   j3   j3   j  j j   j2 j    j j j   j 1  ,   0,  K  j  j  f j, j j  , 0 1j   0j hz j j 3 j  0,  L 0  0j ,  Mj  j j    1  j  f j,  0 ,  M  hz  j  j  j  , 0 j  0,  L 0 57  j 1  j   1j  j 1  j  ,  0j 1  0,  Kj 1  Các phương pháp phân rã điều kiện ban đầu   g , j  1,2,3, Như trình bày chương ta có tính ổn định ba lược đồ phân rã giải tốn nhiễm khí (3.1.1) , (3.1.2) , (3.1.3) nghiệm toán sai phân xấp xỉ nghiệm tốn vi phân với cấp xác O(h   ) Các thể dẫn đến hệ phương trình đại số tuyến tính với ma trận có dạng ba đường chéo Từ phương trình sơ đồ thứ ta có: j (E   1j )  F1 j  (E    j3 ) j  F3 j với F1 j  (E   1j ) j 1 với F3 j   (E  12   j3 ) f j  (E  12   j3 ) hay   k   (1 ) k  F1k   m   ( 3 ) m  F3m ký hiệu k   klm j j , 1  1j , F1k  ( F1 j ) k ,  m   klm ,    j3 , F3m  ( F3 j ) m Sử dụng toán tử sai phân từ (4.2.4) vào hai biểu thức ta có   ( x h  u k 1    u ) k 1  (1  ) k  (  k 1 ) k 1  F1k 2hx hx 2hx hx   m1 / 2 ( h z  wm1 /  ( m1 /   m1 / ) ) m1  (1    ) m hz 2hz2 w    ( m21 /  m1 ) m1  F3m hz hz Đặt hệ số biểu thức trên: j , 58   u a1k  (  k 1 ) hx 2hx   u b1k  (  k 1 ) hx 2hx c1k    hx2 w   a3m  ( m21 /  m1/ ) hz 2hz w   b3m  ( m21 /  m1 ) hz 2hz c3m    ( m1/   m1/ ) 2hz2  Với giả thiết  ,  dương, thành phần u, v, w véc tơ vận tốc gió khơng âm khơng tăng theo biến khơng gian, bước lưới biến khơng gian thích hợp ta có hệ số thoả mãn điều kiện ổn định phương pháp truy đuổi 3.4 Một số kết thực nghiệm Để kiểm tra hội tụ phương pháp, tơi cài đặt thuật tốn ngơn ngữ lập trình matlab thu số kết lược đồ sai phân thứ Trường hợp X  Y  H  T  1;         1; u  v  1; w  g 0 f ( x, y, z, t )  x(1  x) y (1  y )( z  1) z  (1  x) y (1  y )( z  1) z t   x(1  x)(1  y )( z  1) z t  y (1  y )( z  1) z t   x(1  x)( z  1) z t  x(1  x) y (1  y )(12 z  12 z  2)t   x(1  x) y (1  y )( z  1) z t Nghiệm đúng:  ( x, y, z, t ) = x(1  x) y(1  y)( z  1)2 z 2t T* K L M  10 10 10 10 0.0111 20 20 20 20 0.0042 30 30 30 30 0.0025 59 50 50 50 50 0.0013 100 100 100 100 6.03310-4 Trong T* = T ,  = max k ,l ,m, j  ( x k , y l , z m , t j )   *kj ,l , m ,  * nghiệm  gần lấy nút lưới sai phân Trường hợp X  Y  H  T  1;         1; u  2; v  1; w  g 0 f ( x, y, z, t )  100 xyzt K  L  M  T *  100 Sự phân bố nồng độ khí thải độ cao z  0.5 thời điểm t  T mô tả đồ thị sau Trường hợp X  Y  H  T  1;         1; u  2; v  1; w  K  L  M  50, T  10 * Nếu hàm vế phải f ( x, y, z, t )  Q ( x  0.5) ( y  0.5) ( z  0.5) , Q  10 nguồn điểm đặt vị trí (0.5,0.5,0.5) , đồ thị cho biết phân bố nồng độ nhiễm có dạng: 60 3.5 Kết luận chƣơng Đưa lược đồ phân rã theo thời gian cài đặt thuật toán kết số toán sai phân xấp xỉ nghiệm toán vi phân xây dựng Chương Cài đặt thuật toán chạy thử nghiệm phương pháp phân rã giải tốn nhiễm khí 61 KẾT LUẬN VÀ HƢỚNG PHÁT TRIỂN Một số kết đạt đƣợc - Tìm hiểu q trình xây dựng phân tích mơ hình toán học khác truyền tải khuếch tán vật chất mơi trường khí Đồng thời xem xét phương trình phương trình liên hợp dựa lý thuyết nhiễu lý thuyết toán ngược Các phiếm hàm khảo sát đại lượng tích phân đặc trưng cho trình truyền tải khuếch tán vật chất Việc chứng minh tính nghiệm toán toán liên hợp thực nhờ đẳng thức tích phân - Nghiên cứu lược đồ sai phân Crank-Nicolson tìm nghiệm gần toán khuếch tán, đảm bảo ổn định phương pháp toán tử toán nửa xác định dương, phân tích nhược điểm phương pháp toán tử tổng hai toán tử khơng giao hốn, từ đề xuất phương pháp phân rã giải toán đảm bảo ổn định phương pháp với độ xác cấp hai theo thời gian mà khơng địi hỏi tính giao hốn hai toán tử - Nghiên cứu phương pháp phân rã nhiều thành phần giải tốn khơng đảm bảo độ xác cấp hai theo thời gian toán tử toán tổng nhiều toán tử nửa xác định dương lược đồ sai phân tìm nghiệm gần tốn khuyếch tán không gian ba chiều với điều kiện biên Direchlet, việc tìm nghiệm tốn sai phân thực phương pháp truy đuổi vận tốc gió khơng tăng theo biến khơng gian - Cài đặt thuật tốn tìm nghiệm số tốn, mơ tả nghiệm đồ thị phản ánh chất vật lý tốn tính đắn phương pháp Hƣớng phát triển Tiếp tục xây dựng hồn thiện phần mềm: Q trình xây dựng phần mềm tính tốn thực qua bước sau: Xây dựng mơ hình, phân tích mơ hình, rời rạc hóa mơ hình, xây dựng phương pháp thuật tốn giải số, cài đặt thuật toán, xây dựng giao diện, triển khai đóng gói Trong bước trên, tác giả dừng lại 62 cài đặt thuật toán kiểm tra tính đắn thuật tốn mà chưa hồn thiện bước xây dựng giao diện triển khai đóng gói sản phẩn Trên sở kết đạt luận văn, thời gian tới, tác giả tiếp tục tìm hiểu nghiên cứu nhằm đề xuất phương pháp hữu hiệu giải toán khuếch tán có tính đến địa hình phức tạp [2, 6], giải tốn tìm vị trí tối ưu cho phép xây dựng xí nghiệp cơng nghiệp đảm bảo tiêu vệ sinh từ điều kiện ràng buộc sức tải sinh thái khỏi độ nhiễm bẩn Việc thực luận văn khó tránh khỏi số sai sót định Học viên mong nhận đóng góp ý kiến thầy cơ, bạn đọc quan tâm để luận văn hoàn thiện Một lần học viên xin cảm ơn Thầy giáo TS Nguyễn Đình Dũng tận tình giúp đỡ, hướng dẫn thời gian thực đề tài, cảm ơn giúp đỡ gia đình, bạn bè đồng nghiệp thời gian qua Thái Nguyên ngày tháng năm 2020 63 TÀI LIỆU THAM KHẢO I Tài liệu tiếng Việt [1] Trần Ngọc Chấn (2000), Ơ nhiễm khơng khí xử lý khí thải, Tập 1, Nhà xuất khoa học kỹ thuật, Hà Nội [2] Bùi Tá Long, Nguyễn Châu Mỹ Dun (201 ), Mơ hình hóa nhiễm khơng khí điều kiện địa hình phức tạp – trường hợp nguồn thải điểm, Tạp chí Khí tượng thủy văn, pp 34-45 [3] Nguyễn Cơng Điều (1 3), Giải số tốn nhiễm khí phương pháp luân hướng, Tin học điều khiển, (2), pp 14-21 [4] Tạ Văn Đĩnh (2002), Phương pháp sai phân phương pháp phần tử hữu hạn, Nhà xuất khoa học Kỹ thuật, Hà Nội [5] Quang,.V.V., Dung,N.D (2008), “Các kết ứng dụng phương pháp phân rã giải tốn nhiễm khí quyển”, Tạp chí Khoa học Cơng nghệ Đại học Thái Ngun, (2), 48-54 [6] Hieu,.H.T.N., Vu, H.A., Long, B.B (2011), “Xây dựng hệ thống tích hợp đánh giá nhiễm khơng khí phương tiện giao thơng đường Huế”, Tạp chí Khoa học Cơng nghệ, pp 333-342 II Tài liệu tiếng Anh [7] Dang Quang A, Ngo Van Luoc (1993), Numerical solution of astationary problem of atmospheric pollution, Vietnam Journal of math, 24(3), pp 315-325 [8] Dang Quang A (2002), Monotone difference schemes for solving some problems of air polution, Advances in Nature Sciences, Proc of the NCST of Vietnam, 6(1), pp 11-21 [9] Dang Quang A & Ngo Van Luoc (1992), Extract solution of astationary problem of air pollution, Proc of the NCST of Vietnam, 4(2), pp 39-46 64 [10] Ngo Van Luoc, Dang Quang A & Nguyen Cong Dieu (1993), Analytic and numerical solution of some problems of air pollution, SEA Bull Math, pp 105-117 [11] Bang Quoc Ho (2017), Vietnam and evaluation of its impact on human health Sustainable Environment Research, Modeling PM10 in Ho Chi Minh City, 27, 95-102 [12] S A Konglok and N Pochai, “Numerical computations of three-dimensional air-quality model with variations on atmospheric stability classes and wind velocities using fractional step method” IAENG International Journal of Applied Mathematics, vol 46, no 1, pp 112–120, 2016 ... ? ?Nghiên cứu phương pháp phân rã xây dựng phần mềm giải tốn nhiễm khí quyển? ?? làm luận văn tốt nghiệp Mục tiêu luận văn tìm hiểu phương pháp phân rã giải tốn khơng dừng xây dựng ứng dụng tồn ô. .. CÔNG NGHỆ THÔNG TIN VÀ TRUYỀN THÔNG PHẠM THANH NGHỊ NGHIÊN CỨU PHƢƠNG PHÁP PHÂN RÃ VÀ XÂY DỰNG PHẦN MỀM GIẢI BÀI TỐN Ơ NHIỄM KHÍ QUYỂN Chun ngành: Khoa học máy tính Mã số: 48 01 01 LUẬN VĂN THẠC... dụng phương pháp phân rã tốn nhiễm khí Phƣơng pháp luận nghiên cứu - Phương pháp nghiên cứu lý thuyết: Tổng hợp, nghiên cứu tài liệu mơ hình tốn học vấn đề mơi trường - Phương pháp nghiên cứu