Đang tải... (xem toàn văn)
Mô hình tính toán song song giải các bài toán biên phức tạp dựa trên tư tưởng chia miền
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn ĐẠI HỌC THÁI NGUYÊN KHOA CÔNG NGHỆ THÔNG TIN Cao Thị Anh Thư Mô hình tính toán song song giải các bài toán biên phức tạp dựa trên tư tưởng chia miền Chuyên ngành: Khoa học máy tính Mã số: 60.48.01 Luận văn thạc sỹ Khoa học máy tính Người hướng dẫn Khoa học: TS. Vũ Vinh Quang Thái Nguyên - 2009 Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn MỤC LỤC ĐẶT VẤN ĐỀ . 2 Chương 1: Các kiến thức cơ bản về giải số phương trình đạo hàm riêng . 4 1.1 PHƯƠNG PHÁP SAI PHÂN . 4 1.2 THUẬT TOÁN THU GỌN KHỐI LƯỢNG TÍNH TOÁN 6 1.2.1 Bài toán biên thứ nhất 6 1.2.2 Bài toán biên thứ hai 12 1.3 ÁP DỤNG ĐỐI VỚI PHƯƠNG TRÌNH ELLIPTIC . 15 1.3.1 Bài toán biên Dirichlet . 15 1.3.2 Bài toán biên hỗn hợp 16 1.4 PHƯƠNG PHÁP LẶP VÀ CÁC SƠ ĐỒ LẶP CƠ BẢN 18 1.4.1 Không gian năng lượng 18 1.4.2 Phương pháp lặp giải phương trình toán tử 19 Chương 2: Cơ sở Toán học của phương pháp chia miền 27 2.1 CÔNG THỨC ĐA MIỀN VÀ PHƯƠNG TRÌNH STEKLOV- POICARE 28 2.2 CÁC PHƯƠNG PHÁP LẶP ĐƠN CƠ SỞ 30 2.2.1 Phương pháp Dirichlet-Neumann 30 2.2.2 Phương pháp Neumann-Neumann 31 2.2.3 Phương pháp Robin 31 2.3 MỘT SỐ THUẬT TOÁN CHIA MIỀN . 33 2.3.1 Thuật toán chia miền Patrick Le Talle. . 33 2.3.2 Thuật toán chia miền J.R.Rice, E.A. Vavalis, Daopi Yang 35 Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 2.3.3 Thuật toán chia miền Saito-Fujita 37 2.3.4 Phương pháp DQuangA-VVQuang 38 2.3.5 Phương pháp chia miền giải bài toán biên gián đoạn mạnh 40 Chương 3: Mô hình tính toán song song giải bài toán Elliptic dựa trên chia miền . 43 3.1 CÁC BƯỚC LẶP TRÊN NHIỀU MIỀN CON . 43 3.2 MÔ HÌNH TÍNH TOÁN SONG SONG GIẢI BÀI TOÁN BIÊN GIÁN ĐOẠN MẠNH 45 3.2.1.Hướng tiếp cận hiệu chỉnh đạo hàm 46 3.2.2. Hướng tiếp cận hiệu chỉnh hàm . 47 3.3. CÁC KẾT QUẢ THỰC NGHIỆM . 49 3.4. ỨNG DỤNG MÔ HÌNH SONG SONG GIẢI BÀI TOÁN CƠ HỌC 51 3.4.1 Sơ đồ song song theo hướng hiệu chỉnh đạo hàm . 53 3.4.2 Sơ đồ song song theo hướng hiệu chỉnh hàm 57 3.4.3 Các kết quả thực nghiệm 60 NHẬN XÉT KẾT LUẬN 63 DANH MỤC CÔNG TRÌNH ĐÃ CÔNG BỐ LIÊN QUAN ĐẾN LUẬN VĂN . 64 TÀI LIỆU THAM KHẢO 65 PHỤ LỤC . 68 Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 1 LỜI CẢM ƠN Sau một thời gian nghiên cứu và thực hiện luận văn thạc sỹ chuyên ngành Khoa học máy tính, đến nay luận văn :"Mô hình tính toán song song giải các bài toán biên phức tạp dựa trên tư tưởng chia miền" của tôi đã được hoàn thiện và đầy đủ. Để có được kết quả như mong muốn tôi luôn nhận được sự quan tâm, chỉ bảo sự giúp đỡ từ thầy giáo hướng dẫn: Tiến sĩ Vũ Vinh Quang - Phó trưởng Khoa Công nghệ thông tin- Đại học Thái Nguyên. Nhân dịp này tôi xin trân trọng gửi lời cảm ơn của mình tới các thầy giáo, các vị giáo sư của Viện Công nghệ Thông tin, các thầy cô giáo thuộc Khoa Công nghệ thông tin - Đại học Thái Nguyên đã truyền đạt những kiến thức bổ ích cho các học viên cao học khoá 6 nơi tôi được học tập và nghiên cứu trong suốt 2 năm qua. Tôi xin bày tỏ tình cảm và lời cảm ơn chân thành nhất tới các đồng nghiệp Viễn thông Thái Nguyên, tới bạn bè người thân và gia đình đã khích lệ, động viên, giúp đỡ tôi trong thời gian qua. Một lần nữa tôi xin gửi lời cảm ơn sâu sắc nhất tới thầy giáo Vũ Vinh Quang đã hướng dẫn, tạo điều kiện để tôi được học tập và nghiên cứu hoàn thiện luận văn của mình. Tôi xin trân trọng cảm ơn! Thái Nguyên, ngày 30 tháng10 năm 2009. Học viên Cao Thị Anh Thư Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 2 ĐẶT VẤN ĐỀ Lý thuyết về phương pháp chia miền đã được phát triển trong vòng 20 năm qua, xuất phát từ công thức đa miền và phương trình biên chung Steklov-Poincare, các phương pháp chia miền được phát triển từ các sơ đồ lặp cơ bản như: Sơ đồ Dirichlet-Neumann, sơ đồ Neumann-Neumann và sơ đồ Robin được nghiên cứu bởi tác giả trên thế giới. Có thể thấy cơ sở của các phương pháp đều xuất phát từ giá trị điều kiện trên biên phân chia từ đó xây dựng các sơ đồ lặp dạng hai lớp đối với phương trình toán tử. Việc nghiên cứu tính chất hội tụ của các sơ đồ lặp sử dụng kết quả của các không gian Sobolev và toán tử Steklov-Poincare. Nội dung chính của luận văn là trên cơ sở của lý thuyết chia miền, luận văn đề xuất mô hình tính toán song song giải quyết các bài toán với điều kiện biên rất phức tạp trên tư tưởng chia miền, tiến hành cài đặt thử nghiệm mô hình đồng thời ứng dụng mô hình song song giải quyết một bài toán trong môi trường vật lý bán dẫn. Luận văn cấu trúc gồm 3 chương: Chương 1: Đưa ra cơ sở về phương pháp lưới, thuật toán thu gọn khối lượng tính toán giải phương trình lưới và cơ sở lý thuyết về các sơ đồ lặp tổng quát. Chương 2: Trình bày tóm tắt cơ sở toán học về phương pháp chia miền, các sơ đồ lặp cơ bản trong phương pháp chia miền. Một số phương pháp chia miền của các tác giả trên thế giới và đặc biệt là các sơ đồ lặp trên tư tưởng hiệu chỉnh hàm hoặc đạo hàm trên biên phân chia của các tác giả Việt Nam và Nhật Bản, phương pháp chia miền đối với bài toán biên gián đoạn mạnh. Chương 3: Trên cơ sở của các sơ đồ lặp theo hướng hiệu chỉnh hàm và đạo hàm, luận văn đề xuất sơ đồ tính toán song song dựa trên tư tưởng hiệu Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 3 chỉnh hàm hoặc đạo hàm, tiến hành tính toán bằng số so sánh hai sơ đồ tính toán song song và đồng thời áp dụng phương pháp song song giải quyết một bài toán cơ học được các tác giả trên thế giới quan tâm. Các kết quả lý thuyết được kiểm tra bằng các chương trình thực nghiệm lập trình trong môi trường MATLAB trên máy tính PC. Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 4 Chương 1 CÁC KIẾN THỨC CƠ BẢN VỀ GIẢI SỐ PHƢƠNG TRÌNH ĐẠO HÀM RIÊNG Trong chương này, chúng tôi trình bày một số kiến thức liên quan đến việc giải số phương trình đạo hàm riêng bao gồm cơ sở của phương pháp lưới, thuật toán thu gọn khối lượng tính toán và lý thuyết về phương pháp lặp giải phương trình toán tử. Những kiến thức cơ sở và kết quả được tham khảo từ các tài liệu [ 5, 10, 16, 21]. 1.1 Phƣơng pháp sai phân Lƣới sai phân: Xét bài toán ,,,.u f xu g x (1.1) trong đó 2( , ) , ,x y R a x b c y d , chọn 2 số nguyên >1N và >1M, đặt = ( ) /h b a N gọi là bước lưới theo x, = ( ) /k d c M gọi là bước lưới theo y. Đặt = , = , 0 , 0 .ijx a ih y c jk i N j M Mỗi điểm ( , )ijxy gọi là một nút lưới ký hiệu là nút ( , )ij. Tập tất cả các nút trong ký hiệu là hk. Nút ở trên biên gọi là nút biên; tập tất cả các nút biên ký hiệu là hk, tập =hk hk hk gọi là một lưới sai phân trên . Hàm lƣới: Mỗi hàm số xác định tại các nút của lưới gọi là một hàm lưới, giá trị của hàm lưới ( , )u x y tại nút lưới ( , )ij viết tắt là ,iju. Mỗi hàm ( , )u x y xác định tại mọi ( , )xy tạo ra hàm lưới u xác định bởi ,iju. Bài toán sai phân: Ký hiệu Lu f là tập các hàm số hai biến ,xy có các đạo hàm riêng đến cấp m liên tục trong = Giả sử bài toán có nghiệm 4()uC, khi đó: Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 5 41( , )4| ( , )| =xyumax x y C constx, 42( , )4| ( , )| =xyumax x y C consty. Do đó theo công thức Taylor ta có: 1( , ) = ( ) ,i j i ju x y u x h y 2 2 3 3423= ( , ) ( )2! 3!iju h u h uu x y h o hx x x hay 211222( , ) 2 ( , ) ( , )= ( )i j i j i ju x y u x y u x yuohhx Một cách tương tự: 1( , ) = ( , )i j i ju x y u x y k 2 2 3 3423= ( , ) ( )2! 3!iju k u k uu x y k o ky y y 1( , ) ( , )i j i ju x y u x y k 2 2 3 3423= ( , ) ( )2! 3!iju k u k uu x y k o ky y y Do đó: 211222( , ) 2 ( , ) ( , )= ( )i j i j i ju x y u x y u x yuokky Vậy ta có: 1 1 1 12222( , ) 2 ( , ) ( , ) ( , ) 2 ( , ) ( , )=()i j i j i j i j i j i ju x y u x y u x y u x y u x y u x yhku o h k Ta đặt: 1, , -1, , 1 , , -122- 2 - 2 -1i j i j i j i j i j i jhku u u u u uuhk Khi đó chứng tỏ: 22= ( )khu u o h k Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 6 Số hạng 22O(h +k ) là một vô cùng bé bậc hai. Ta nói toán tử kh xấp xỉ toán tử , điều đó cho phép thay phương trình vi phân bằng phương trình sai phân: = , = ( , ), ( , )hk ij ij i j i j hku f f f x y x y tức là: 1, , 1 , 1 , ,222 2 1( , ), ( , )i j i j i j i j i j i ji j i j hku u u u u uf x y x yhk (1.2) đồng thời thay điều kiện biên bằng điều kiện: ( , ), ( , )ij i j i j hku g x y x y (1.3) Ta được bài toán sai phân hoàn chỉnh: tìm hàm lưới utại các nút ( , )ij thoả mãn hệ phương trình sai phân (1.2) với điều kiện biên (1.3). Như vậy việc tìm nghiệm xấp xỉ của bài toán vi phân (1.1) với độ chính xác cấp hai được đưa về việc giải bài toán sai phân (1.2) với điều kiện (1.3) bằng các phương pháp đại số. 1.2 Thuật toán thu gọn khối lƣợng tính toán Được đề xuất bởi Samarski-Nicolaev. Bằng các phép biến đổi đơn giản về vec tơ và ma trận, các bài toán sai phân luôn luôn được đưa về hệ phương trình vec tơ 3 điểm thuộc một trong các dạng sau đây: 1.2.1 Bài toán biên thứ nhất Xét bài toán biên thứ nhất đối với phƣơng trình véc tơ ba điểm 11=j j j jY CY Y F , 11jN , 00=YF, =NNYF. (1.4) Trong đó jY là véc tơ cần tìm, C là ma trận vuông, jF là véc tơ cho trước. ý tưởng của phương pháp rút gọn hoàn toàn giải (1.1) là khử liên tiếp các ẩn jY đầu tiên với các j lẻ, sau đó từ các phương trình còn lại khử các jY Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 7 với j là bội của 2, rồi bội của 4,… Mỗi bước khử sẽ giảm được một nửa số ẩn. Như vậy nếu =2nN thì sau một số lần khử sẽ còn lại một phương trình chứa véc tơ ẩn /2NY mà từ đó /2NY có thể tính được qua 0Y và NY . Sau khi đã có được 0 /2,NYY và NY thì quá trình ngược lại là việc tìm các jY với j là bội của 4N rồi bội của 8N ,… Rõ ràng, phương pháp rút gọn hoàn toàn là một biến thể của phương pháp khử Gauss áp dụng cho bài toán (1.4) trong đó việc khử các biến được thực hiện theo một thứ tự đặc biệt. Sau đây, ta sẽ mô tả cụ thể phương pháp. Giả sử = 2 , > 0nNn Ký hiệu (0) (0)= , = ; =1,2, ., 1jjC C F F j N . Khi đó (1.4) được viết dưới dạng 0 (0)11= (1 1)j j j jY C Y Y F j N , 00=YF, =NNYF. (1.5) Bước khử thứ nhất: Từ các phương trình đầu của (1.5) ta khử các jY với j lẻ. Muốn vậy ta viết 3 phương trình liên tiếp: (0) (0)2 1 1=j j j jY C Y Y F , (0) (0)11=j j j jY C Y Y F , (0) (0)1 2 1=j j j jY C Y Y F Nhân 2 vế của phương trình thứ hai với (0)C vào bên trái rồi cộng cả 3 phương trình lại ta được (1) (1)2 2 1= , = 2,4, ., 2j j j jY C Y Y F j N , 00=YF, =NNYF (1.6) trong đó: (1) ( 2= ( 0)) 2C C E (1) (0) (0) (0) (0)11= , = 2,4, ., 2j j j jF F c F F j N . Nhận xét rằng hệ (1.6) chỉ chứa các jY với j chẵn, số véc tơ ẩn jY là 12N. Do đó nếu giải được hệ này thì các jY với j lẻ sẽ tìm được từ phương trình [...]... của lý thuyết chia miền tổng quát, trên thế giới đã xuất hiện một số thuật toán chia miền của các tác giả áp dụng đối với các trường hợp cụ thể 2.3.1 Thuật toán chia miền Patrick Le Talle Thuật toán đưa ra nhằm sử dụng phương pháp chia miền để giải các bài toán biên Dirichlet và Neumann, đây là cơ sở cho việc tính toán đa nhiệm trên máy tính CRAY2 và INTER, đã được kiểm nghiệm trên các bài toán đàn hồi... ta đưa ra cơ sở toán học của phương pháp chia miền bao gồm giới thiệu các khái niệm về các điều kiện chuyển giao giữa các biên chung, các công thức biến phân và đặc biệt là ứng dụng của toán tử Steklov-Poincare đối với phương pháp chia miền các phương pháp lặp đơn trên các biên chung Các kiến thức được trình bày trên cơ sở các tài liệu [11,12, 14, 22, 25, 26, 29, 30, 31] Hình 1 Xét bài toán u f... x , trong đó các hàm , phụ thuộc vào từng loại bài toán, Với bài toán Poisson thì (v) v, (v) v n 2.2 Các phƣơng pháp lặp đơn cơ sở Trong phần này, chúng ta xét việc giải bài toán đa miền bằng các thủ tục lặp, chúng ta xét 1 dãy các bài toán con trong 1 , 2 với các điều kiện biên Dirichlet hoặc Neumann tư ng ứng Các phương pháp đó có thể thực k hiện được bởi 1 trong các sơ đồ lặp sau... kiểm nghiệm trên các bài toán đàn hồi trong không gian 3 chiều, bài toán miền ảo, bài toán trong các ngành công nghiệp trên quy mô lớn Xét mô hình đơn giản u f , x , x u 0, (2.19) Chia 1 2 bởi biên chung S Kí hiệu là giá trị hàm u trên S , khi đó ta có thể tiến hành giải song song hai bài toán trong hai miền i , i 1, 2 x i , u i f , x i , ... Trong các tài liệu của Samaski-Nicolaev [21] đã chứng minh độ phức tạp của các thuật toán là O ( M N log N ) 1.3 Áp dụng đối với phƣơng trình elliptic Trên cơ sở phương pháp lưới, ta thu được các kết quả xây dựng lược đồ sai phân cho các bài toán Dirichlet và bài toán Neumann 1.3.1 Bài toán biên Dirichlet Cho là hình chữ nhật = x = ( x1 , x2 ) R2 : 0 < x1 < l1;0 < x2 < l2 Xét bài toán u... nhất và bài toán biên thứ hai, áp dụng đối với bài toán biên Dirichlet và bài toán biên hỗn hợp, cơ sở lý thuyết về phương pháp lặp giải phương trình toán tử Những kiến thức quan trọng này làm nền tảng cho các kết quả sẽ trình bày trong các chương tiếp theo của luận văn Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 27 Chương 2 CƠ SỞ TOÁN HỌC CỦA PHƢƠNG PHÁP CHIA MIỀN Trong... trị trên biên chung phải thoả mãn phương trình SteklovPoincare S , x , trong đó (2.7) 2 G f G2 f G1 f i , i 1 n n n (2.8) trong đó S là toán tử Steklov-Poincare được định nghĩa bởi 2 H i H 1 H 2 S i n n i 1 n Cùng với toán tử S , ta cũng sử dụng các toán tử S i1 và gọi là các toán tử Poincare-Steklov Mô hình chia miền trên có thể áp dụng đối với bài toán. .. 1 j 1.2.2 Bài toán biên thứ hai Xét bài toán thứ hai Y0 = F0 , j = 0, Yj 1 CYj Yj 1 = Fj ,1 j N 1, 2YN 1 CYN = FN (1.18) trong đó N = 2n , n > 0 Để giải bài toán (1.18) ta cũng thực hiện các bước khử lần lượt như đã được trình bày ở bài toán biên thứ nhất Sau n phép khử, ta nhận được các phương trình Y0 = F0( n ) , 2Y(0) C ( n )YN = FN( n ) (1.19) Và nhóm các phương... định các véc tơ v (0) = q (j k 1) Y j 2( k 1) Y j 2( k 1) , j = 2k 1 ,3.2k 1 , , N 2k 1 j Sau đó, với l = 1,2, ,2k 1 và với mỗi j = 2k 1 ,2.2k 1 , , N 2k 1 , giải phương trình Cl ,k 1v (j l ) = v (j l 1) Khi đó Y j = p (j k 1) v (2 j k 1 ) , j = 2k ,2.2k ,3.2k , N 2k Trên đây là nội dung của thuật toán thu gọi khối lượng tính toán giải bài toán biên thứ nhất và bài toán biên. .. n, n 1, ,1 Các công thức trên đã mô tả phương pháp rút gọn hoàn toàn giải Việc tính các Fj( k ) theo công thức truy toán có thể dẫn đến việc tích luỹ sai số nếu như chuẩn của ma trận C ( k 1) lớn hơn 1 Ngoài ra các ma trận C ( k ) nói chung là các ma trận đầy đủ, thậm chí cả với ma trận ban đầu là C (0) = C là ma trận ba đường chéo Điều này dẫn đến tăng khối lượng tính toán khi tính các Fj( k ) theo . Cao Thị Anh Thư Mô hình tính toán song song giải các bài toán biên phức tạp dựa trên tư tưởng chia miền Chuyên ngành: Khoa học máy tính Mã số: 60.48.01. ngành Khoa học máy tính, đến nay luận văn :" ;Mô hình tính toán song song giải các bài toán biên phức tạp dựa trên tư tưởng chia miền& quot; của tôi
