Bài toán biên ban đầu đối với phương trình Parabolic với điều kiện biên Robinson trong miền bị chặn

Header Page of 258 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM HÀ NỘI - NGUYỄN THU PHƯƠNG XẤP XỈ GALERKIN ĐỐI VỚI PHƯƠNG TRÌNH PARABOLIC TUYẾN TÍNH VỚI ĐIỀU KIỆN BIÊN ROBINSON Chuyên ngành: Toán giải tích Mã số: 60 46 01 02 LUẬN VĂN THẠC SĨ TOÁN HỌC Người hướng dẫn khoa học:TS NGUYỄN THÀNH ANH Hà Nội, 2016 Footer Page of 258 Header Page of 258 LỜI CẢM ƠN Trước trình bày nội dung luận văn, tác giả xin chân thành cảm ơn trường Đại học Sư Phạm Hà Nội 2, nơi mà tác giả hoàn thành chương trình cao học giảng dạy nhiệt tình tâm huyết Thầy, Cô Đặc biệt, tác giả xin bày tỏ lòng biết ơn sâu sắc tới TS Nguyễn Thành Anh, người Thầy trực tiếp hướng dẫn, bảo tận tình giúp đỡ để tác giả hoàn thành luận văn Cuối cùng, xin gửi lời cảm ơn tới gia đình, bạn bè, người giúp đỡ chia sẻ với tác giả suốt thời gian học tập hoàn thành luận văn Hà Nội, ngày 10 tháng năm 2016 Nguyễn Thu Phương Footer Page of 258 Header Page of 258 LỜI CAM ĐOAN Tôi xin cam đoan luận văn kết nghiêm cứu riêng hướng dẫn TS Nguyễn Thành Anh Trong trình nghiên cứu, kế thừa thành khoa học với trân trọng biết ơn Các kết trích dẫn luận văn rõ nguồn gốc Hà Nội, ngày 10 tháng năm 2016 Nguyễn Thu Phương Footer Page of 258 Header Page of 258 MỤC LỤC Mở đầu Kiến thức chuẩn bị 1.1 Không gian Sobolev 1.2 Một số bất đẳng thức 10 Xấp xỉ Galerkin phương trình parabolic tuyến tính với điều kiện biên Robinson 12 2.1 Phát biểu toán 12 2.2 Đánh giá sai số xấp xỉ Galerkin nửa rời rạc 14 2.3 2.2.1 Xây dựng nghiệm xấp xỉ 14 2.2.2 Ước lượng sai số chuẩn L2 16 2.2.3 Ước lượng sai số chuẩn H 20 Rời rạc hóa thời gian 22 2.3.1 Bước thời gian cách 23 2.3.2 Bước thời gian thay đổi 28 Kết luận 31 Tài liệu tham khảo 32 Footer Page of 258 Header Page of 258 MỞ ĐẦU Lý chọn đề tài Các toán biên ban đầu phương trình parabolic mô hình toán học nhiều ứng dụng thực tế khác Giải số toán thường dùng phương pháp phần tử hữu hạn không gian Trong khuôn khổ đề tài luận văn quan tâm đến giải số toán biên ban đầu phương trình parabolic tuyến tính với điều kiện biên không thứ ba (điều kiện biên Robinson) Bởi vậy, hướng dẫn TS Nguyễn Thành Anh, chọn đề tài: “ Xấp xỉ Galerkin phương trình Parabolic tuyến tính với điều kiện biên Robinson” Luận văn hoàn thành dựa báo [4]:"Galerkin Approximations for the Linear Parabolic Equation with the Third Boundary Condition" tác giả I Faragó, S Korotov, P Neittaanmaki, công bố năm 2003 Mục đích nghiên cứu Nghiên cứu hội tụ, ước lượng sai số dãy nghiệm xấp xỉ Galerkin toán biên ban đầu phương trình parabolic tuyến tính với điều kiện biên Robinson Nhiệm vụ nghiên cứu - Nghiên cứu hội tụ, ước lượng sai số không gian L2 , H dãy nghiệm xấp xỉ Galerkin Footer Page of 258 Header Page of 258 - Ước lượng sai số nghiệm xấp xỉ rời rạc theo thời gian Đối tượng phạm vi nghiên cứu Bài toán biên ban đầu phương trình Parabolic với điều kiện biên Robinson miền bị chặn Phương pháp nghiên cứu Nghiên cứu tài liệu tham khảo theo phương pháp: Hệ thống lại kiến thức có liên quan, phân tích, tổng hợp Dự kiến đóng góp Luận văn đóng góp mặt khoa học Footer Page of 258 Header Page of 258 CHƯƠNG KIẾN THỨC CHUẨN BỊ 1.1 Không gian Sobolev Cho Ω miền Rd (d = 1, 2, 3) ∞ C ∞ (Ω) ={ u : Ω → R khả vi vô hạn } = ∩ C k (Ω) k=0 Cc∞ (Ω) kí hiệu hàm C ∞ (Ω) với giá compact Lp (Ω) = { u : Ω → R | u độ đo Lebesgue, u 1/ p p , p ≥ |u| dx u Lp (Ω) = Lp (Ω) < ∞ }, Ω Llloc (Ω) = { u : Ω → R| u ∈ Ll (Ω ) , ∀Ω ⊂⊂ Ω } Cho trước đa số α, kí hiệu ∂ |α| u (x) α1 αn D u (x) = α1 αn = ∂x1 ∂xn u ∂x1 ∂xn α Không gian Sobolev lớp không gian dùng nhiều trình nghiên cứu phương trình đạo hàm riêng Để đến định nghĩa lớp không gian này, trước tiên phải tìm hiểu khái niệm "đạo hàm yếu" Footer Page of 258 Header Page of 258 Giả sử u hàm khả vi Khi với "hàm thử" φ ∈ Cc∞ (Ω), sử dụng công thức tích phân phần ta thu đẳng thức sau du dφ φ dx = − u dx Ω dxj Ω dxj Lặp lại trình |α| lần, tương tự ta có Dα uφ dx = (−1)|α| uDα φ dx Ω Ω với đa số α Chúng ta có định nghĩa đạo hàm yếu sau: Định nghĩa 1.1.1 Với hàm u ∈ L1loc (Ω), ta nói v đạo hàm yếu u ứng với biến xj , ký hiệu v = Dj u, v ∈ L1loc (Ω) vφ dx = − Ω u Ω dφ dx dxj với φ ∈ Cc∞ (Ω) Bằng cách quy nạp, định nghĩa đạo hàm yếu cấp cao sau: Định nghĩa 1.1.2 Nếu u, v ∈ L1loc (Ω) v gọi đạo hàm yếu cấp α u, viết v = Dα u, vφ dx = (−1)|α| Ω uDα φ dx Ω với φ ∈ Cc∞ (Ω) Dựa vào tính chất khả tích đạo hàm yếu, ta đưa định nghĩa không gian Sobolev Định nghĩa 1.1.3 Không gian Sobolev định nghĩa W r,p (Ω) = u : Dα u ∈ Lp (Ω), với ≤ |α| ≤ r , với chuẩn u W r,p   =  0≤|α|≤r Footer Page of 258 Dα u 1/p  p Lp (Ω)  Header Page of 258 Không gian Sobolev W r,p (Ω) định nghĩa không gian Banach khả ly Tương tự không gian L2 (Ω), không gian Sobolev W r,p (Ω), trường hợp p = dùng nhiều nghiên cứu Vì lý đó, sau nghiên cứu chủ yếu không gian W r,2 (Ω) sở L2 (Ω) Vì không gian đặc biệt dùng thường xuyên không gian Sobolev khác nên có ký hiệu riêng W r,2 (Ω) = H r (Ω) Người ta chọn ký hiệu H r (Ω) không gian Hilbert với tích vô hướng trang bị sau Dα uDα vdx < u, v >H r = 0≤|α|≤r Ω Khi đó, chuẩn u ∈ H r (Ω) tương ứng với tích vô hướng xác định sau 1/2  u r= Dα u   L2 (Ω) 0≤|α|≤r Chúng giới thiệu nửa chuẩn  |u|r =  1/2 |Dm v|2 dx |m|=r Ω ∞ Định nghĩa 1.1.4 Ta kí hiệu Wr,p (Ω) bao đóng Cc (Ω) Wr,p (Ω) ∞ Do đó, u ∈ Wr,p (Ω) tồn hàm um ∈ Cc (Ω) cho um → u Wr,p (Ω) r,p α Ta nói Wr,p (Ω) gồm hàm u ∈ W (Ω) cho "D u = ∂u" với |α| r − Chú ý: Ta viết H0r (Ω) = Wr,2 (Ω) H −1 (Ω) không gian đối ngẫu H01 (Ω) Footer Page of 258 Header Page 10 of 258 Định nghĩa 1.1.5 Giả sử X không gian Banach Không gian Lp (0, T ; X) bao gồm tất hàm đo u : [0, T ] → X với u Lp (0,T ;X) u (t) 1/p p X dt