Phương pháp Collocation với cơ sở B-Spline bậc năm giải phương trình truyền nhiệt một chiều

128 TRƯỜNG ĐẠI HỌC THỦ ĐÔ HÀ NỘI 
 PHƢƠNG PHÁP COLLOCATION VỚI 
 CƠ SỞ B-SPLINE BẬC NĂM GIẢI PHƢƠNG TRÌNH 
 TRUYỀN NHIỆT MỘT CHIỀU 
 Nguyễn V n Tuấn1, Nguyễn Thị Thƣ Hòa 
 Trường Đại học Thủ đô Hà Nội 
 Tóm tắt: Trong bài báo này chúng ta nghiên cứu cách giải phương trình truyền nhiệt một 
 chiều bằng phương pháp collocation với cơ sở là các hàm B-spline bậc năm. Sự ổn định 
 Von Neumann của lược đồ sai phân và so sánh kết quả giữa nghiệm đúng và nghiệm xấp 
 xỉ cũng được trình bày. 
 Từ khóa: Phương pháp spline collocation, Quintic B – spline, phương pháp phần tử hữu 
 hạn. 
1. MỞ ĐẦU 
 Xét phƣơng trình truyền nhiệt một chiều dạng: 
 (1) 
với điều kiện đầu: 
 u(x,0) = f(x) (2) 
và các điều kiện biên: 
 { (3) 
 Trong đó: là các hằng số, là các hàm số liên tục với 
 Phƣơng trình (1) với các điều kiện (2), (3) mô tả dòng nhiệt trong vật dẫn khối trụ. Cụ 
thể qua nghiên cứu ch ng ta có thể biết đƣợc dòng nhiệt trên một thanh dẫn chiều dài L với 
sự khuếch tán của dòng nhiệt , là hệ số khuếch tán. 
1 Nhận bài ngày 24.03.2016; gửi phản biện và duyệt đăng ngày 10.05.2016 
 Liên hệ tác giả: Nguyễn Văn Tuấn; Email: nvtuan@daihocthudo.edu.vn 
TẠP CHÍ KHOA HỌC SỐ 4/2016 129 
 Ngoài ra, nhiều hiện tƣợng vật lí có thể lý giải khi giải phƣơng trình (1) với các điều 
kiện khác nhau. Do vậy, các nhà toán học trong và ngoài nƣớc quan tâm nghiên cứu bài 
toán (1) với các điều kiện (2) và (3) bằng nhiều cách giải khác nhau ([3], [4]). 
 Trong bài báo này, ch ng ta nghiên cứu giải gần đ ng bài toán trên bằng phƣơng pháp 
collocation với cơ sở là các hàm B-spline bậc năm. 
2. NỘI DUNG 
2.1. Phƣơng pháp spline collocation 
 Giả sử chia đoạn [a, b] thành N phân bằng nhau bởi các điểm n t: 
 a = x0 < x1 <  < xN = b, h = (b – a)/N. 
 Xác định các hàm B-spline cơ sở bậc 5 ([6]) nhƣ sau: 
 ( ) ( ) 
 ( ) ( ) ( ) 
 ( ) 
 ( ) ( ) ( ) ( ) 
 ( ) 
 ( ) ( ) ( ) ( ) 
 ( ) ( ) 
 { 
 Tập các hàm N tạo thành một cơ sở các hàm B – spline bậc 5 
xác định trên [a, b]. Giá trị của và các đạo hàm bậc nhất, bậc hai của nó tại xj đƣợc 
xác định theo bảng 1 ( và các đạo hàm bậc nhất, bậc hai của nó bằng 0 ngoài khoảng 
 , 
 Bảng 1. Giá trị của 
 x 
 1 0 26 66 26 1 0 
 0 0 0 
 0 0 
130 TRƯỜNG ĐẠI HỌC THỦ ĐÔ HÀ NỘI 
 Ch ng ta tìm nghiệm xấp xỉ U(x, t) của bài toán (1), (2), (3) của nghiệm đ ng u(x, t) 
dƣới dạng: 
 U ∑ (4) 
Trong đó: 
 U U 
 { U U (5) 
 U U 
 Sử dụng điều kiện collocation cho (5) tại các điểm xm, m = 0, , N, ta có: 
 U U N (6) 
 Thay (4) vào (6) ta có: 
 ∑ ∑ N (7) 
 Giả sử là nội suy tuyến tính giữa hai mức thời gian n và n + 1 thì: 
 trong đó và là các ẩn tại mức thời gian n. 
 Sử dụng phƣơng pháp sai phân hữu hạn ta có: 
 Thay (8) vào (7) và chọn = ½ ta nhận đƣợc: 
 ∑ ∑ 
 { 
 N 
 (9) 
 Sử dụng bảng 1, tính các số hạng của hệ phƣơng trình (9) tại các điểm lƣới xm, khi đó 
ta có: 
 (10) 
Trong đó: 
 { 
 Hệ phƣơng trình (10) gồm N+1 phƣơng trình với N+5 ẩn 
( , để tìm nghiệm ch ng ta cần 4 phƣơng trình nữa. Sử dụng 
(5) ta có: 
TẠP CHÍ KHOA HỌC SỐ 4/2016 131 
 ∑ 
 ∑ 
 (11) 
 ∑ ∑ 
 { 
 Thay vào (11) chúng ta có: 
 (12) 
 { 
 với 
 Giải hệ phƣơng trình (12) ch ng ta nhận đƣợc: 
 (13) 
 { 
 Thay (13) vào hệ phƣơng trình (10) ta đƣợc hệ N+1 phƣơng trình với N + 1 ẩn 
( sau: 
 (14) 
 { 
Với: 
132 TRƯỜNG ĐẠI HỌC THỦ ĐÔ HÀ NỘI 
 Để giải hệ (14) trƣớc tiên ta giải hệ phƣơng trình (15) sau: 
 U N 
 U 
 (15) 
 U 
 {U U 
 Khử các ẩn của hệ phƣơng trình (15) ta đƣợc hệ phƣơng trình 
(16): 
 A H (16) 
 với A là ma trận 5 đƣờng chéo 
 54 60 6 0 0 0 ... 0
 101 135 105
 1 0 0 ... 0
 4 2 4
 1 26 66 26 1 0 ... 0
 ... ... ...
 A 
 ... ... ...
 0 ... 0 1 26 66 26 1
 105 135 101
 0 ... 0 0 1
 4 2 4
 0 ... 0 0 0 6 60 54
TẠP CHÍ KHOA HỌC SỐ 4/2016 133 
 H . 
2.2. Sự ổn định 
 Ch ng ta sẽ chứng minh hệ phƣơng trình sai phân (10) ổn định Von – Neumann. 
 n
Đặt  { } với √ là số mode. Khi đó phƣơng trình (10) trở thành: 
 111
 . exp 2i  h exp 2i  h . exp i  h exp i  h 
  222
 
 222. exp 2i  h exp 2i  h . exp i  h exp i  h 
 222 
 
 2 cos2 (  h)  cos(  h) 1 
 1 11
 2 .
 
 2 cos2 (  h)  cos(  h) 2 
 222 2
 2 1
 2 1 x  1 x 1
 Để tìm miền giá trị của  ta xét hàm số: y 2 
 
 2 x2  x 2 
 2 22 2
Với: 1 x 1. 
 Đạo hàm của y ta có: 
 120h22 t(2x 6x 7)
 y., 
 [(h2 10  )x 2 (13h 2 10  )x 16h 2 20  t] 2
Dễ thấy: y‟(x) > 0, với Nên y đồng biến trên khoảng đã cho. 
Mặt khác ta nhận đƣợc: 
 8h2 40  t
 y( 1)1 
 8h2 40  t 
 y(1) 1.
Do đó: 1 y(x) 1,  x [-1, 1]  1. 
Vậy (10) ổn định vô điều kiện. 
2.3. Kết quả số 
 Xét phƣơng trình truyền nhiệt: 
134 TRƯỜNG ĐẠI HỌC THỦ ĐÔ HÀ NỘI 
 uuxx 0,0 t x 1, (17) 
Với điều kiện đầu: 
 u(x,0) sin( x), (18) 
Các điều kiện biên: 
 u(0,t) u(1,t) 0,t 0,
 2
 ux (0, t) exp - t
 (19) 
 2
 ux (1,t)exp - t 
 uxxxx (0, t) u (1,t) 0.
Bài toán (17), (18), (19) có nghiệm đ ng: u(x,t) exp - 2 t sin x. 
Kết quả số cho theo các bảng sau: 
 Bảng 2. So sánh kết quả số với t 0,0001;h 0,0125 
 x t Nghiệm xấp xỉ Nghiệm đúng 
 0,1 0,30302 0,30153 
 0,2 0,11315 0,11238 
 0,3 0,04218 0,04186 
 0,4 0,01574 0,01561 
 0,3 0,5 0,00588 0,00582 
 0,6 0,00221 0,00217 
 0,7 0,00084 0,00081 
 0,8 0,00033 0,00030 
 0,9 0,00014 0,00011 
 0,1 0,35725 0,35446 
 0,2 0,13302 0,13211 
 0,3 0,04958 0,04924 
 0,4 0,01849 0,01835 
 0,6 0,5 0,00691 0,00684 
 0,6 0,00259 0,00255 
 0,7 0,00098 0,00095 
 0,8 0,00038 0,00035 
 0,9 0,00016 0,00014 
 0,1 0,11628 0,11519 
 0,2 0,04322 0,04293 
 0,3 0,01612 0,01600 
 0,4 0,00602 0,00596 
 0.9 0,5 0,00226 0,00222 
TẠP CHÍ KHOA HỌC SỐ 4/2016 135 
 0,6 0,00086 0,00083 
 0,7 0,00033 0,00031 
 0,8 0,00014 0,00012 
 0,9 0,00007 0,00004 
 Bảng 3. So sánh kết quả với t = 0,5 v t0,0001 
 Nghiệm xấp xỉ 
 x Nghiệm đúng 
 h = 0,05 h = 0,025 h = 0,1667 h = 0,0125 
 0,1 0,00230 0,00228 0,00226 0,00226 0,00222 
 0,2 0,00430 0,00429 0,00428 0,00428 0,00423 
 0,3 0,00589 0,00589 0,00589 0,00588 0,00582 
 0,4 0,00691 0,00692 0,00691 0,00691 0,00684 
 0,5 0,00726 0,00727 0,00727 0,00726 0,00719 
 0,6 0,00691 0,00692 0,00691 0,00691 0,00684 
 0,7 0,00589 0,00589 0,00589 0,00588 0,00582 
 0,8 0,00430 0,00429 0,00428 0,00428 0,00423 
 0,9 0,00230 0,00228 0,00227 0,00226 0,00223 
 Hình 1. Đồ thị của h m cơ sở B-spline bậc 5 
136 TRƯỜNG ĐẠI HỌC THỦ ĐÔ HÀ NỘI 
 Hình 2. Đồ thị đạo h m của B – spline bậc 5 
 Hình 3. Đồ thị mặt nghiệm xấp xỉ 
3. KẾT LUẬN 
 Bài báo đã trình bày phƣơng pháp collocation trong đó sử dụng hệ cơ sở B – spline 
bậc năm giải xấp xỉ phƣơng trình truyền nhiệt một chiều. Sự ổn định của hệ phƣơng trình 
sai phân tƣơng ứng đã đƣợc chứng minh. Đồng thời qua ví dụ khẳng định tính hiệu quả của 
phƣơng pháp. 
 TÀI LIỆU THAM KHẢO 
1. G. Arora, R. C. Mittal, B. K. Singh (2014), “Numerical solution of BBM – Burger equation 
 with quartic B – spline collocation method”, J. of Engineering, Special issue on ICMTEA 
 2013 conference, December, pp.104-116. 
TẠP CHÍ KHOA HỌC SỐ 4/2016 137 
2. Behnam Sepehrian, Mahmood Lashami (2008), “A numerical solution of the Burgers 
 equation using quintic B – spline”, Proceeding of the World Congress on Engineering, Vol. 
 III, WCE 2008, London, U.K. . 
3. Duygu Dӧnmer Demiz, Necdet Bildik (2012), “The numerical solution of Heat problem using 
 cubic B – spline”, Applied Mathematics, 2(4), pp.131-135. 
4. Joan Goh, Ahmad Abd. Majid, and Ahmad Jzani Md. Ismail (2012), “Cubic B – spline 
 collocation method for one – dimensional Heat and advection – diffusion equations”, J. of 
 Applied Mathematics, Vol., Article IO 458710. 
5. A. A. Karawia, “Two algorithms for solving a general backward pentadiagonal linear 
 systems”,  
6. P. M. Prenter (2008), “Spline and variational methods”, Dover Publications, New York. 
 QUINTIC B – SPLINE COLLOCATION METHOD 
 FOR ONE – DIMESIONAL HEAT EQUATION 
 Abstract: This paper discusses solving one the dimensional heat equation. Numerical 
 solutions are obtained by collocation method based on quintic B – spline. The stability 
 analysis of the scheme is examined by the Von Neumann approach. On the other hand, a 
 comparative study between the numerical and the exact is illustrated. 
 Keywords: Collocation method, B – spline, Finite element method