Phương trình sóng phi tuyến bị nhiễu xấp xỉ tuyến tính và dáng điệu tiệm cận của nghiệm

20 258 0
Phương trình sóng phi tuyến bị nhiễu   xấp xỉ tuyến tính và dáng điệu tiệm cận của nghiệm

Đang tải... (xem toàn văn)

Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống

Thông tin tài liệu

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM TP HỒ CHÍ MINH Nguyễn Văn Ý PHƯƠNG TRÌNH SÓNG PHI TUYẾN BỊ NHIỄU: XẤP XỈ TUYẾN TÍNH VÀ DÁNG ĐIỆU TIỆM CẬN CỦA NGHIỆM LUẬN VĂN THẠC SĨ TOÁN HỌC Thành phố Hồ Chí Minh – 2008 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM TP HỒ CHÍ MINH Nguyễn Văn Ý PHƯƠNG TRÌNH SÓNG PHI TUYẾN BỊ NHIỄU: XẤP XỈ TUYẾN TÍNH VÀ DÁNG ĐIỆU TIỆM CẬN CỦA NGHIỆM Chuyên ngành: Toán giải tích Mã số: 60 46 01 LUẬN VĂN THẠC SĨ TOÁN HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC: TS NGUYỄN CÔNG TÂM Thành phố Hồ Chí Minh – 2008 LỜI CẢM ƠN Lời đầu tiên, xin kính gửi đến Thầy Nguyễn Công Tâm lời cảm ơn sâu sắc giúp đỡ Thầy việc hoàn thành luận văn Xin chân thành cảm ơn thầy Lê Hoàn Hóa thầy Nguyễn Thành Long đọc cho nhiều ý kiến đóng góp bổ ích Xin chân thành cảm ơn quý Thầy Cô thuộc Khoa Toán – Tin học Trường Đại học Sư phạm Thành phố Hồ Chí Minh Đại học khoa học tự nhiên Thành phố Hồ Chí Minh tận tình giảng dạy suốt khóa học Xin cảm ơn quý Thầy Cô thuộc phòng Khoa học Công nghệ Sau đại học, thư viện trường Đại học Sư phạm Thành phố Hồ Chí Minh, Ban Giám hiệu trường THPT Chuyên Hùng Vương – Bình Dương tạo điều kiện thuận lợi để hoàn thành chương trình học Xin cảm ơn anh chị bạn lớp Cao học giải tích K16 anh chị bạn nhóm xemina Thầy Nguyễn Thành Long chủ trì hỗ trợ nhiều mặt thời gian học tập nghiên cứu Và cuối cùng, lời thân thương xin gửi đến gia đình tôi, nơi tạo cho điều kiện thuận lợi để học tập hoàn thành luận văn Thành phố Hồ Chí Minh, tháng 09 – 2008 Nguyễn Văn Ý MỤC LỤC Trang phụ bìa Lời cảm ơn Mục lục MỞ ĐẦU Chương 1: CÁC CÔNG CỤ CHUẨN BỊ 1.1 Các không gian hàm thông dụng 1.2 Không gian hàm Lp (0, T ; X ) , ≤ p ≤ ∞ 1.3 Phân bố có giá trị không gian Banach 1.4 Đạo hàm Lp (0, T ; X ) 1.5 Bổ đề tính compact Lions 1.6 Bổ đề hội tụ yếu Lq (Q) 1.7 Một số kết khác Chương 2: SỰ TỒN TẠI VÀ DUY NHẤT NGHIỆM 10 2.1 Giới thiệu 10 2.2 Thuật giải xấp xỉ tuyến tính 10 2.3 Sự tồn nghiệm 21 Chương 3: SỰ HỘI TỤ CẤP HAI 25 3.1 Dãy lặp cấp hai 25 3.2 Sự hội tụ bậc hai 34 Chương 4: DÁNG ĐIỆU NGHIỆM KHI λ → + , δ → + 40 4.1 Dáng điệu tiệm cận nghiệm yếu theo hai tham số λ , δ 40 4.2 Dáng điệu tiệm cận nghiệm yếu theo tham số λ 42 4.3 Khai triển tiệm cận theo tham số λ đến cấp N + 45 Chương 5: KHẢO SÁT MỘT TRƯỜNG HỢP CỤ THỂ 54 KẾT LUẬN 56 TÀI LIỆU THAM KHẢO 57 MỞ ĐẦU Các toán biên phi tuyến xuất khoa học ứng dụng ( Vật lý, Hóa học, Cơ học, Kỹ thuật, ) phong phú đa dạng Đây nguồn đề tài mà nhiều nhà toán học từ trước đến quan tâm nghiên cứu, chẳng hạn [2, – 12, 14 – 22] tài liệu tham khảo Chính vậy, cho đề tài nghiên cứu cần thiết, có ý nghĩa lý luận thực tiễn Trong luận văn nầy, muốn sử dụng phương pháp Giải tích hàm phi tuyến như: phương pháp Galerkin, phương pháp compact đơn điệu, phương pháp xấp xỉ tuyến tính liên hệ với định lí điểm bất động, phương pháp khai triển tiệm cận nhằm khảo sát toán biên có liên quan đến vấn đề khoa học ứng dụng Trong luận văn này, xét toán giá trị biên ban đầu sau ∂ σ (u , u x ) + λut = F ( x, t ), < x < 1, < t < T , ∂x (0.1) utt − (0.2) u (0, t ) = u (1, t ) = 0, (0.3) u ( x,0) = u0 ( x), ut ( x,0) = u1 ( x), (0.4) σ (u, u x ) = u x + δ f (u ), λ , δ số không âm cho trước; f , F , u0 , u1 , hàm cho trước thỏa điều kiện mà ta sau Trường hợp σ = σ (u x , u xt ) có nhiều công trình nghiên cứu Khởi đầu với trường hợp σ = σ (u x , u xt ) = β (u x ) + λu xt , λ > 0, β ∈ C (\), β (0) = 0, β ′ ≥ ε > 0, toán (0.1) – (0.3) xét Greenberg, MacCamy, Mizel [10] Đây mô hình toán học mô tả dao động dọc đàn hồi nhớt phi tuyến, u ( x, t ) độ dịch chuyển so với vị trí cân Từ xuất công trình [10], có nhiều công trình công bố liên quan đến toán nầy, chẳng hạn như: Greenberg [11], Greenberg, MacCamy [12], Dafermos [6], Andrews [2], Clements [4] Phương trình (0.1) với (0.4) mặt hình thức có dạng (0.5) utt − u xx = f ( x, t , u, u x , ut ), f ( x, t , u, u x , ut ) = F ( x, t ) + δ f ′(u )u x − λut , nhiên mặt ý nghĩa có điểm khác biệt riêng Trong [9], Ficken Fleishman chứng minh tồn nghiệm phương trình (0.6) u xx − utt − 2α1ut − α 2u = ε u + b, với ε > bé Trong báo Caughey Ellison [5], hợp xấp xỉ trường hợp trước để bàn tồn tại, tính ổn định tiệm cận nghiệm cổ điển cho hệ động lực phi tuyến liên tục Trong [7], Alain Phạm nghiên cứu tồn tại, dáng điệu tiệm cận ε → nghiệm yếu toán (0.3), (0.5) liên kết với điều kiện biên Dirichlet (0.7) u (0, t ) = u (1, t ) = 0, số hạng phương trình (0.5) cho (0.8) f = ε f1 (t , u, ut ) Nếu f1 ∈ C N ([0,∞] × \ ) thỏa f1 (t ,0,0) = với t ≥ 0, khai triển tiệm cận nghiệm toán (0.3), (0.7), (0.8) đến cấp N + theo ε thu với ε đủ nhỏ Trong [14, 15], Long Alain Phạm nghiên cứu toán (0.3), (0.5) với số hạng phi tuyến có dạng f = f1 (u , ut ) Trong [14], tác giả xét với điều kiện biên hỗn hợp (0.9) u x (0, t ) = hu (0, t ) + g (t ), u (1, t ) = 0, h > số cho trước [15] với điều kiện tổng quát t (0.10) u x (0, t ) = g (t ) + hu (0, t ) − ∫ k (t − s )u (0, s)ds, u (1, t ) = 0 Trong [16], Long Diễm nghiên cứu toán (0.3), (0.5) với điều kiện biên hỗn hợp (0.11) u x (0, t ) − h0u (0, t ) = u x (1, t ) + h1u (1, t ) = 0, h0 , h1 số không âm cho trước với h0 + h1 > số hạng phi tuyến vế phải có dạng (0.12) f = f ( x, t , u, u x , ut ) + ε f1 ( x, t , u, u x , ut ) Trong trường hợp f ∈ C ([0,1] × [0, ∞) × \ ), f1 ∈ C ([0,1] × [0, ∞) × \ ) tác giả thu khai triển tiệm cận nghiệm yếu uε đến cấp hai theo ε , với ε đủ nhỏ [16] Trong luận văn nghiên cứu tồn nghiệm địa phương toán (0.1) – (0.4) Chứng minh nhờ vào phương pháp Galerkin liên kết với đánh giá tiên nghiệm với kỹ thuật hội tụ yếu tính compact Tiếp đến khảo sát tồn dãy lặp cấp hai hội tụ dãy nghiệm yếu toán (0.1) – (0.4) tương ứng Sau nghiên cứu dáng điệu tiệm cận nghiệm yếu uλ ,δ (phụ thuộc λ , δ ) toán (0.1) – (0.4) (λ , δ ) → (0+ ,0+ ) khai triển tiệm cận nghiệm yếu toán nhiễu theo tham số nhiễu λ đến cấp N + 1, tức nghiệm xấp xỉ đa thức theo λ (0.13) N uλ ( x, t ) ≈ ∑ ulk ( x, t )λ k , k =0 theo nghĩa cần hàm ulk ( x, t ), k = 0,1, , N , thiết lập đánh giá theo dạng (0.14) N uλt − ∑ ulk t λ k k =0 N L∞ (0,T ; L2 ) + uλ − ∑ ulk λ k k =0 L∞ (0,T ; H 01 ) ≤ CT λ N +1 Luận văn trình bày theo chương mục sau: Phần mở đầu, tổng quan toán khảo sát luận văn, điểm qua kết có trước đó, đồng thời nêu bố cục luận văn ¾ Chương 1, trình bày số kết chuẩn bị bao gồm việc nhắc lại số không gian hàm, số kết phép nhúng compact không gian hàm ¾ Chương 2, tập trung nghiên cứu tồn nghiệm yếu toán (0.1) – (0.4) ¾ Chương 3, tập trung nghiên cứu tồn hội tụ dãy lặp cấp hai nghiệm yếu toán (0.1) – (0.4) ¾ Chương 4, nghiên cứu dáng điệu tiệm cận nghiệm yếu toán (0.1) – (0.4) khai triển tiệm cận theo tham số λ đến cấp N + ¾ Chương 5, minh họa cho khai triển tiệm cận theo λ chương trường hợp cụ thể 5 Chương 1: CÁC CÔNG CỤ CHUẨN BỊ 1.1 Các không gian hàm thông dụng Đầu tiên ta đặt kí hiệu sau Ω = (0,1), QT = Ω × (0, T ), T > bỏ qua định nghĩa không gian hàm thông dụng: C m (Ω), Lp (Ω), H m (Ω), W m , p (Ω) Để cho gọn, ta kí hiệu lại sau: Lp = Lp (Ω), H m = H m (Ω) = W m ,2 (Ω), W m, p = W m , p (Ω) ( xem [1, 3]) Ta định nghĩa L2 = L2 (Ω) không gian Hilbert với tích vô hướng (1.1.1) u , v = ∫ u ( x)v( x)dx, u, v ∈ L2 Kí hiệu để chuẩn sinh tích vô hướng (1.1.1), nghĩa ⎛ ⎞ u , u = ⎜ ∫ u ( x) dx ⎟ , u ∈ L2 ⎝0 ⎠ (1.1.2) u = Định nghĩa không gian Sobolev cấp (1.1.3) H = {v ∈ L2 : v′ ∈ L2 } Không gian không gian Hilbert tích vô hướng (1.1.4) u, v Kí hiệu (1.1.5) u H1 H1 H1 = = u, v + u′, v′ để chuẩn sinh tích vô hướng (1.1.4), nghĩa u, u H1 , u ∈ H Ta có bổ đề liên hệ hai không gian L2 H sau Bổ đề 1.1 Phép nhúng H -C (Ω) compact (1.1.6) v C0 (Ω) ≤ v H1 , ∀v ∈ H Chứng minh Xem Adams[1] 6 Ta sử dụng không gian Sobolev đặc biệt không gian (1.1.7) H = D (Ω) H1 H1 ∞ c = C (Ω ) Bổ đề 1.2 Ta có phép nhúng từ H 01 -C (Ω) compact (1.1.8) ⎧ v C ( Ω ) ≤ vx = v H ∀v ∈ H 01 , ⎪ ⎨ v H ≤ vx = v H ≤ v H ∀v ∈ H 01 ⎪ ⎩ Chứng minh bổ đề 1.2 không khó khăn Một cách đặc trưng khác để xác định H 01 (1.1.9) H 01 = {v ∈ H : v(0) = v(1) = 0} Bổ đề 1.3 Đồng L2 với ( L2 )′ (đối ngẫu L2 ) Khi ta có H 01 - L2 ≡ ( L2 )′ - ( H 01 )′ ≡ H −1 với phép nhúng liên tục nằm trù mật Chú thích 1.2 Từ Bổ đề 1.3, ta dùng ký hiệu tích vô hướng ⋅, ⋅ L2 để cặp tích đối ngẫu ⋅, ⋅ H −1 , H 01 ký hiệu Ta ký hiệu ⋅ H 01 H −1 Chuẩn L2 để chuẩn không gian Banach X X gọi X ′ không gian đối ngẫu X 1.2 Không gian hàm Lp (0, T ; X ), ≤ p ≤ ∞ Cho X không gian Banach thực với chuẩn ⋅ X Ta kí hiệu Lp (0, T ; X ), ≤ p ≤ ∞ không gian lớp tương đương chứa hàm u : (0, T ) → X đo cho ⎛ = ⎜ ∫ u (t ) ⎝0 T u p L (0,T ; X ) p ⎞ dt ⎟ < ∞, ≤ p < ∞ X ⎠ p u Lp (0,T ; X ) = ess sup u (t ) X , p = ∞ 0 0, T > 0, ta đặt: (2.2.5) (2.2.6) (2.2.7) { } K i = K i ( M , f ) = sup f ( i ) (u ) : u ≤ M (i = 1, 2), W ( M , T ) = {v ∈ L∞ (0, T ; H 01 ∩ H ) : v ∈ L∞ (0, T ; H 01 ), v∈ L2 (QT ), v L∞ (0,T ; H 01 ∩ H ) , v L∞ (0,T ; H 01 ) , v L2 (Q ) ≤ M }, T W1 ( M , T ) = {v ∈W ( M , T ) : v∈ L∞ (0, T ; L2 )} Tiếp theo, ta xây dựng dãy {um } W1 ( M , T ) qui nạp chứng minh hội tụ nghiệm yếu toán (2.2.1) – (2.2.3) với lựa chọn M > 0, T > thích hợp Ta xét thuật giải xấp xỉ tuyến tính sau: Chọn số hạng đầu u0 ∈W1 ( M , T ) Giả sử (2.2.8) um−1 ∈W1 ( M , T ) Ta liên kết toán (2.2.1) – (2.2.3) với toán biến phân sau: Tìm um ∈W1 ( M , T ) cho (2.2.9) (2.2.10) um (t ), v + a (um (t ), v) + λ um (t ), v = Fm (t ), v ∀v ∈ H 01 , um (0) = u0 , um (0) = u1 , 12 (2.2.11) Fm ( x, t ) = if ( x, t , um−1 ( x, t ), ∇um−1 ( x, t )) = F ( x, t ) + δ f ′(um−1 ( x, t ))∇um−1 ( x, t ) Sự tồn um cho định lí sau Định lí 2.1 Giả sử ( H1 ) − ( H ) Khi đó, tồn số dương M , T dãy qui nạp tuyến tính {um } ⊂ W1 ( M , T ) xác định (2.2.9) – (2.2.11) Chứng minh Gồm bước sau Bước Xấp xỉ Galerkin Xét sở trực giao Hilbert {w j } H 01 trực chuẩn L2 mục 2.1 Đặt (2.2.12) k (k ) um( k ) (t ) = ∑ cmj (t ) w j , j =1 (k ) cmj (t ) thỏa hệ phương trình vi phân tuyến tính sau: (2.2.13) um( k ) (t ), w j + a(um( k ) (t ), w j ) + λ um( k ) (t ), w j = Fm (t ), w j , ≤ j ≤ k , (2.2.14) um( k ) (0) = u0 k , um( k ) (0) = u1k , ≤ j ≤ k , (2.2.15) k (k ) u0 k = ∑α mj w j → u0 mạnh H 01 ∩ H , j =1 (2.2.16) k u1k = ∑ β mj( k ) w j → u1 mạnh H 01 j =1 Hệ (2.2.13), (2.2.14) viết thành dạng khác (2.2.17) (k ) (k ) (k ) ⎧⎪cmj (t ) + λ j2cmj (t ) + λ cmj (t ) = Fm (t ), w j , ⎨ (k ) (k ) (k ) (k ) ⎪⎩cmj (0) = α mj , cmj (0) = β mj , ≤ j ≤ k Từ (2.2.17)1 ta có 13 (k ) (k ) (cmj (t )eλt )′t + λ j2cmj (t )eλt = eλt Fm (t ), w j , ≤ j ≤ k Do ta suy c (t ) = α (k ) mj (2.2.18) (k ) mj t r β mj( k ) (k ) − λt + (1 − e ) − λ j ∫ dr ∫ eλ ( s −r )cmj ( s )ds λ 0 t r 0 + ∫ dr ∫ eλ ( s −r ) Fm ( s ), w j ds Bổ đề 2.1 Giả sử um−1 thỏa (2.2.8) Khi hệ (2.2.13) có nghiệm um( k ) (t ) khoảng ≤ t ≤ T Chứng minh bổ đề 2.1 Bỏ qua số m, k cách viết ta (k ) (k ) (t ), α mj , β mj( k ) Ta viết lại hệ (2.2.18) viết c j (t ), α j , β j thay cho cmj thành phương trình điểm bất động (2.2.19) c(t ) = (Uc)(t ), ≤ t ≤ T , (2.2.20) ⎧c = (c1 , , ck ), Uc = ((Uc)1 , ,(Uc) k ), ⎪(Uc) = γ (t ) + (Vc) (t ), j j j ⎪ t r ⎪ β ⎨γ j (t ) = α j + j (1 − e − λt ) + ∫ dr ∫ eλ ( s − r ) Fm ( s ), w j ds, λ 0 ⎪ t r ⎪ ⎪(Vc) j (t ) = −λ j2 ∫ dr ∫ eλ ( s − r ) c j ( s )ds, ≤ j ≤ k 0 ⎩ Ta cần chứng minh: tồn số tự nhiên p0 cho toán tử U p0 : X = C ([0, T ]; \ k ) → X co, tức tồn số ρ ∈ [0,1) cho (2.2.21) U p0 c − U p0 d X ≤ ρ c−d X ∀c, d ∈ X , đây, ta sử dụng chuẩn X sau (2.2.22) c k X = sup ∑ c j (t ) , c ∈ X 0≤t ≤T j =1 Sau đây, qui nạp ta chứng minh ∀p ∈ `, ta có 14 k ∑ (U pc) j (t ) − (U p d ) j (t ) ≤ (2.2.23) j =1 (σ t ) p c−d (2 p )! X ∀c, d ∈ X , ∀t ∈ [0, T ] , với σ = kπ • Với p = 1, ta có k t r 0 j =1 k ∑ (Uc) (t ) − (Ud ) (t ) ≤ (kπ ) ∫ dr ∫ ∑ c (s) − d (s) ds j j =1 (2.2.24) j (σ t) ≤ c−d 2! X j j Vậy (2.2.23) với p = • Giả sử (2.2.23) với p ≥ 1, tức k ∑ (U pc) j (t ) − (U p d ) j (t ) ≤ (2.2.25) j =1 (σ t ) p c−d (2 p)! X ∀c, d ∈ X , ∀t ∈ [0, T ] • Ta chứng minh (2.2.23) với p + 1, nghĩa (2.2.26) k ∑ (U p +1 j =1 c) j (t ) − (U (σ t ) p + d ) j (t ) ≤ c−d (2 p + 2)! p +1 X ∀c, d ∈ X , ∀t ∈ [0, T ] Thật (2.2.27) k k j =1 j =1 ∑ (U p+1c) j (t ) − (U p+1d ) j (t ) = ∑ (U (U pc)) j (t ) − (U (U p d )) j (t ) ≤σ t r k ∫ dr ∫ ∑ (U p j =1 c) j ( s ) − (U p d ) j ( s ) ds (σ s ) p c−d ≤ σ ∫ dr ∫ (2 p )! 0 t r (σ t ) p + c−d = (2 p + 2)! X X ds ∀c, d ∈ X , ∀t ∈ [0, T ] Từ (2.2.27) suy (2.2.26) Như (2.2.23) ∀p ∈ ` 15 (σ T ) p = nên tồn số tự nhiên p0 cho p →+∞ (2 p )! Mặt khác, lim (σ T ) p0 < Tức toán tử U p0 co Theo nguyên lý ánh xạ co, ta suy (2 p0 )! hệ (2.2.13) có nghiệm um( k ) (t ) ≤ t ≤ T Vậy bổ đề 2.1 chứng minh , Bước Đánh giá tiên nghiệm Đặt (2.2.28) X (k ) m t 2 (t ) = u (t ) + a (u (t ), u (t )) + 2λ ∫ um( k ) ( s) ds, (k ) m (k ) m (k ) m (2.2.29) (k ) m Y t (t ) = a(u (t ), u (t )) + Δu (t ) + 2λ ∫ ∇um( k ) ( s ) ds, (k ) m (k ) m (k ) m t (2.2.30) S (t ) = X (k ) m (k ) m (t ) + Y (k ) m (t ) + ∫ um( k ) ( s ) ds Đánh giá thứ (k ) Nhân (2.2.13) với cmj (t ) lấy tổng theo j, ta (2.2.31) um( k ) (t ), um( k ) (t ) + a(um( k ) (t ), um( k ) (t )) + λ um( k ) (t ) = Fm (t ), um( k ) (t ) Do d (k ) X m (t ) = Fm (t ), um( k ) (t ) dt Suy t (2.2.32) X (k ) m (t ) = X (k ) m (0) + ∫ Fm ( s ), um( k ) ( s ) ds Đánh giá thứ hai Trong (2.2.13) thay w j = − λ j Δw j , sau đơn giản cho λ j2 ta 16 (2.2.33) um( k ) (t ), −Δw j + a (um( k ) (t ), −Δw j ) + λ um( k ) (t ), −Δw j = Fm (t ), −Δw j Hơn nữa, ta có: (2.2.34) (2.2.35) (2.2.36) um( k ) (t ), −Δw j = a(um( k ) (t ), w j ), a (um( k ) (t ), −Δw j ) = Δum( k ) (t ), Δw j , um( k ) (t ), −Δw j = a(um( k ) (t ), w j ), Fm (t ), −Δw j = − ∫ Fm ( x, t )Δw j ( x)dx 1 = − Fm ( x, t )∇w j ( x) + ∫ ∇Fm ( x, t )∇w j ( x)dx (2.2.37) = a ( Fm (t ), w j ) Nhờ vào (2.2.33) – (2.2.37) ta có (2.2.38) a (um( k ) (t ), w j ) + Δum( k ) (t ), Δw j + λ a (um( k ) (t ), w j ) = a ( Fm (t ), w j ) Trong (2.2.38) thay w j um( k ) (t ) ta (2.2.39) a (um( k ) (t ), um( k ) (t )) + Δum( k ) (t ), Δum( k ) (t ) + λ a(um( k ) (t ), um( k ) (t )) = a ( Fm (t ), um( k ) (t )) Do ta có (2.2.40) d (k ) Ym (t ) = a ( Fm (t ), um( k ) (t )) dt Tích phân theo t, ta t (2.2.41) (k ) m Y (t ) = Y (k ) m (0) + ∫ a ( Fm ( s ), um( k ) ( s ))ds Từ (2.2.32) (2.2.41) suy [...]... sẽ chỉ ra sau Trong chương này, ta sẽ thiết lập một định lí tồn tại và duy nhất nghiệm yếu của bài toán (2.1.1) bằng phương pháp xấp xỉ tuyến tính kết hợp với phương pháp Galerkin và phương pháp compact yếu Xét một cơ sở trực giao {w j } của H 01 và trực chuẩn trong L2 gồm các hàm w j = 2 sin( jπx), j = 1,2, được lập từ các hàm riêng của toán tử Laplace ∂2 −Δ = − 2 sao cho −Δw j = λ j2 w j , λ j = jπ... t ) + δ f ′(um−1 ( x, t ))∇um−1 ( x, t ) Sự tồn tại của um được cho bởi định lí sau Định lí 2.1 Giả sử ( H1 ) − ( H 4 ) là đúng Khi đó, tồn tại các hằng số dương M , T và một dãy qui nạp tuyến tính {um } ⊂ W1 ( M , T ) xác định bởi (2.2.9) – (2.2.11) Chứng minh Gồm các bước sau Bước 1 Xấp xỉ Galerkin Xét một cơ sở trực giao Hilbert {w j } của H 01 và trực chuẩn trong L2 như trong mục 2.1 Đặt (2.2.12)... ) ≤ M }, T W1 ( M , T ) = {v ∈W ( M , T ) : v∈ L∞ (0, T ; L2 )} Tiếp theo, ta xây dựng dãy {um } trong W1 ( M , T ) bằng qui nạp và chứng minh nó hội tụ về nghiệm yếu của bài toán (2.2.1) – (2.2.3) với sự lựa chọn M > 0, T > 0 thích hợp Ta xét thuật giải xấp xỉ tuyến tính sau: Chọn số hạng đầu u0 ∈W1 ( M , T ) Giả sử rằng (2.2.8) um−1 ∈W1 ( M , T ) Ta liên kết bài toán (2.2.1) – (2.2.3) với bài toán... đề 1.9 (Bổ đề về tính compact của Lions[13]) Với giả thiết (1.5.1), (1.5.2) và nếu 1 < pi < ∞, i = 0,1 thì phép nhúng W (0, T ) - Lp0 (0, T ; X ) là compact Chứng minh Có thể tìm thấy trong Lions[13], trang 57 9 1.6 Bổ đề về sự hội tụ yếu trong Lq (Q ) Bổ đề 1.10 Cho Q là tập Gm , G ∈ Lq (Q), 1 < q < +∞ sao cho, Gm mở, Lq bị chặn của \N và ≤ C , trong đó C là hằng số độc lập với m và Gm → G hầu khắp... gian Banach thực Một ánh xạ tuyến tính liên tục từ D((0, T )) vào X được gọi là một (hàm suy rộng) phân bố có giá trị trong X Tập hợp các phân bố có giá trị trong X kí hiệu là D′((0, T ; X )) = L( D((0, T )); X ) = {u: D((0, T )) →X: u tuyến tính, liên tục} Chú thích 1.4 Ta kí hiệu D(0, T ) thay cho D((0, T )) hoặc Cc∞ ((0, T )) để chỉ không gian các hàm thực khả vi vô hạn và có giá compact trong (0,... nữa, dãy {w j / λ j } cũng là một cơ sở trực chuẩn Hilbert của H 01 đối với tích vô hướng (2.1.2) u, v H 01 = a (u , v) = u x , vx 2.2 Thuật giải xấp xỉ tuyến tính Ta thành lập các giả thiết sau: ( H1 ) λ > 0, δ ≥ 0, (H2 ) u0 ∈ H 01 ∩ H 2 , u1 ∈ H 01 , (H3 ) F, ∂F ∈ L∞ (0, ∞; L2 ), thỏa F (0, t ) = F (1, t ) = 0, ∀t ≥ 0, ∂x 11 f ∈ C 2 (\) và f ′(0) = 0 (H4 ) Bài toán (2.1.1) được viết lại (2.2.1)... [0, T ] vào X Chứng minh của bổ đề 1.8 có thể tìm thấy trong Lions[13] 1.5 Bổ đề về tính compact của Lions Cho ba không gian Banach X 0 , X 1 , X với X 0 - X - X 1 với các phép nhúng liên tục sao cho (1.5.1) X 0 , X 1 là phản xạ, (1.5.2) Phép nhúng X 0 - X là compact Với 0 < T < +∞, 1 ≤ pi ≤ +∞, i = 0,1 Ta đặt { } W (0, T ) = v ∈ Lp0 (0, T ; X 0 ) : v′ ∈ Lp1 (0, T ; X 1 ) (1.5.3) Ta trang bị cho W... khắp nơi trong Q Khi đó, Gm → G yếu trong Lq (Q) 1.7 Một số kết quả khác Bổ đề sau đây liên quan đến một bất phương trình tích phân và nó rất cần thiết cho việc đánh giá tiên nghiệm trong các chương sau Bổ đề 1.11 (Bổ đề Gronwall) Giả sử f : [ 0, T ] → \ là hàm khả tích, không âm trên [ 0, T ] và thỏa bất t đẳng thức f (t ) ≤ C1 + C2 ∫ f ( s )ds với hầu hết t ∈ [0, T ], trong đó C1 , C2 là 0 các hằng... ) Ta định nghĩa đạo hàm nghĩa phân bố của u bởi công thức (1.3.1) du dϕ ,ϕ = − u , dt dt ∀ϕ ∈ D(0, T ) Bổ đề 1.7 Lp (0, T ; X ) ⊂ D′(0, T ; X ) với phép nhúng liên tục du theo dt 8 Chứng minh của bổ đề 1.7 có thể tìm thấy trong Lions[13] 1.4 Đạo hàm trong Lp (0, T ; X ) Do bổ đề 1.7, u ∈ Lp (0, T ; X ) ta có thể coi u ∈ D′(0, T ; X ) và do đó du dt là phần tử của D′(0, T ; X ) Ta có kết quả sau: Bổ... xỉ Galerkin Xét một cơ sở trực giao Hilbert {w j } của H 01 và trực chuẩn trong L2 như trong mục 2.1 Đặt (2.2.12) k (k ) um( k ) (t ) = ∑ cmj (t ) w j , j =1 (k ) trong đó cmj (t ) thỏa hệ phương trình vi phân tuyến tính sau: (2.2.13) um( k ) (t ), w j + a(um( k ) (t ), w j ) + λ um( k ) (t ), w j = Fm (t ), w j , 1 ≤ j ≤ k , (2.2.14) um( k ) (0) = u0 k , um( k ) (0) = u1k , 1 ≤ j ≤ k , trong đó

Ngày đăng: 08/09/2016, 10:42

Từ khóa liên quan

Tài liệu cùng người dùng

  • Đang cập nhật ...

Tài liệu liên quan