Chỉnh hóa nghiệm cho bài toán nhiệt và bài toán Elastic ngược

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM TP. HỒ CHÍ MINH PHÙNG TRỌNG THỰC CHỈNH HÓA NGHIỆM CHO BÀI TOÁN NHIỆT VÀ BÀI TOÁN ELASTIC NGƯỢC 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: PGS. TS. ĐẶNG ĐỨC TRỌNG Thành phố Hồ Chí Minh - 2010 LỜI CẢM ƠN Em xin chân thành gửi lời cảm ơn đến các thầy cô trường ĐH Sư Phạm và ĐH KHTN đã tận tình giảng dạy chúng em trong suốt thời gian học cao học Toán. Đặc biệt em chân th

pdf57 trang | Chia sẻ: huyen82 | Lượt xem: 1633 | Lượt tải: 0download
Tóm tắt tài liệu Chỉnh hóa nghiệm cho bài toán nhiệt và bài toán Elastic ngược, để xem tài liệu hoàn chỉnh bạn click vào nút DOWNLOAD ở trên
ành cảm ơn thầy Đặng Đức Trọng đã rất ân cần và chu đáo hướng dẫn em làm luận văn này. Học Viên: Phùng Trọng Thực. 1MỞ ĐẦU Nội dung chính của luận văn là đưa ra sự chỉnh hóa nghiệm cho một dạng của bài toán nhiệt hai chiều và bài toán Elastic ba chiều. Cụ thể là đưa ra một chỉnh hóa nghiệm cho các bài toán: Bài toán nhiệt: Cho T > 0 là độ dài của thời gian quan sát và Ω = (0, 1)× (0, 1) là vật dẫn nhiệt. Xác định cặp hàm (u, f) thỏa mãn hệ: ut −∆u = ϕ (t) f (x, y) , ux (0, y, t) = ux (1, y, t) = uy (x, 0, t) = uy (x, 1, t) = 0, u (1, y, t) = 0, u (x, y, 0) = g (x, y) , với (x, y) ∈ Ω, t ∈ (0, T ), trong đó g ∈ L1 (Ω) và ϕ ∈ L1 (0, T ) được cho. Bài toán Elastic: Cho T > 0 là độ dài của thời gian quan sát và Ω = (0, 1) × (0, 1) × (0, 1) là vật thể đàn hồi đẳng hướng ba chiều. Xác định cặp (u, f) thỏa mãn hệ: ∂2u ∂t2 + µ∆u+ (λ+ µ)∇ (div (u)) = ϕ (f1, f2, f3) , (x, t) ∈ Ω× (0, T ) , (u1 (x, t) , u2 (x, t) , u3 (x, t)) = (0, 0, 0) , (x, t) ∈ ∂Ω× (0, T ) , (u1 (x, 0) , u2 (x, 0) , u3 (x, 0)) = (g1 (x) , g2 (x) , g3 (x)) , x ∈ Ω,( ∂u1 ∂t (x, 0) , ∂u2 ∂t (x, 0) , ∂u3 ∂t (x, 0) ) = (h1 (x) , h2 (x) , h3 (x)) , x ∈ Ω, σ1 τ12 τ13τ21 σ2 τ23 τ31 τ32 σ3   n1n2 n3  =  X1X2 X3  , với λ và µ là các hằng số thỏa µ < 0, λ + 2µ < 0, σj = λdiv (u) + 2µ ∂uj ∂xj , 2τjk = µ( ∂uj ∂xk + ∂uk ∂xj ) và n = (n1, n2, n3) là pháp vectơ đơn vị hướng ra ngoài trên ∂Ω. Trong đó dữ kiện được cho là I (ϕ,X, g, h) ∈ ( L1 (0, T ) , ( L1 ( 0, T, L1 (∂Ω) ))3 , ( L2 (Ω) )3 , ( L2 (Ω) )3) . Bài toán nhiệt và bài toán Elastic như trên là những bài toán ngược, không chỉnh. Tính không chỉnh của bài toán ở chỗ bài toán có thể không tồn tại nghiệm hoặc nếu tồn tại duy nhất nghiệm thì nghiệm có thể không phụ thuộc liên tục vào dữ kiện được cho. Trong những năm gần đây, một số tác giả đã có những nghiên cứu về các bài toán này. Chẳng hạn xem xét về sự duy nhất và ổn định nghiệm của bài toán nhiệt trong [12], [14], [15], [16]; sự chỉnh hóa trong trường hợp nghiệm không ổn định cho bài toán nhiệt trong [5], [6], [9]; tính duy nhất nghiệm cho bài toán Elastic trong [8], [11] và đưa ra sự chỉnh hóa nghiệm cho bài toán Elastic hai chiều trong [8]. Bởi vì các bài toán trên là những bài toán ngược không chỉnh nên sự chỉnh hóa nghiệm là cần thiết. Trong [7], các tác giả Trong, Dinh, Nam đã đưa ra một sự chỉnh hóa nghiệm cho bài toán nhiệt bằng phương pháp cắt ngắn chuỗi Fourier và sử dụng một vài kỹ thuật, chẳng hạn phương pháp nội suy để xấp xỉ các hệ số của chuỗi cắt ngắn từ các dữ kiện nhiễu. Ưu điểm của phương pháp này là có thể loại bỏ những giả thiết trên nghiệm về điều kiện cuối của thời gian. Chú ý rằng trong [6], [8], [9] các tác giả cần sử dụng thông tin về điều kiện cuối của thời gian u (x, T ) trong việc chỉnh hóa nghiệm bởi vì nó giúp đưa ra công thức biến đổi Fourier của f và từ đó khôi phục được f . Dựa vào nhận xét phương pháp cắt ngắn chuỗi Fourier cũng có thể áp dụng để chỉnh hóa nghiệm cho bài toán Elastic và giúp ta loại bỏ các giả thiết trên nghiệm về điều kiện cuối của thời gian (điều mà trong [8] khi chỉnh hóa nghiệm cho bài toán Elastic hai chiều các tác giả cần sử dụng đến) nên trong luận văn này sẽ đưa ra một sự trình bày chi tiết cho phương pháp này để chỉnh hóa nghiệm của bài toán Elastic, nhưng so với [8] luận văn có hai điểm mới sau: • Mở rộng xem xét bài toán Elastic trên không gian 3 chiều. • Bỏ đi các ràng buộc trên nghiệm về điều kiện cuối của thời gian. 3Luận văn bao gồm ba chương: Chương 1. Kiến thức chuẩn bị Chương này trình bày một số định nghĩa, các kết quả và một số kiến thức bổ trợ sẽ được dùng đến trong các chương sau. Chương 2. Chỉnh hóa nghiệm cho bài toán nhiệt Trong chương này trình bày về sự duy nhất nghiệm và chỉnh hóa nghiệm cho bài toán nhiệt. Đây là sự trình bày chi tiết các kết quả trong bài báo [7] và qua đó cho thấy phương pháp mà các tác giả đã sử dụng để chỉnh hóa nghiệm cho bài toán nhiệt hai chiều. Chương 3. Chỉnh hóa nghiệm cho bài toán Elastic Trong chương này trình bày về sự duy nhất nghiệm và đưa ra một sự chỉnh hóa nghiệm cho bài toán Elastic ba chiều. Trong chương 2 và 3 đều có phần giải số để minh họa cho các kết quả thu được. Chương 1 Kiến thức chuẩn bị 1.1 Tính không chỉnh và sự chỉnh hóa Tính không chỉnh. Một bài toán ngược gọi là không chỉnh nếu nó thỏa mãn ít nhất một trong ba điều sau: • Bài toán không tồn tại nghiệm. • Bài toán không duy nhất nghiệm. • Nghiệm bài toán không phụ thuộc liên tục vào dữ kiện của bài toán, tức là với một thay đổi nhỏ trên dữ kiện của bài toán có thể dẫn đến thay đổi lớn trên nghiệm của bài toán. Sự chỉnh hóa. Trong trường hợp bài toán ngược tồn tại và duy nhất nghiệm tuy nhiên nghiệm bài toán không ổn định với dữ kiện được cho, tức là sai số nhỏ trên dữ kiện (điều này đúng trong thực tế vì các dữ kiện chỉ là dữ kiện đo đạc khá gần với dữ kiện chính xác) có thể dẫn đến sai số lớn trên nghiệm, khi đó sự chỉnh hóa nghiệm là cần thiết. Chỉnh hóa nghiệm tức là từ các dữ kiện đo đạc (có thể có sai số so với dữ kiện chính xác) ta xây dựng một nghiệm mới, gọi là nghiệm chỉnh hóa. Nghiệm chỉnh hóa có thể không phải là nghiệm chính xác của bài toán (ứng với dữ kiện chính xác) nhưng ta có thể kiểm soát được sai số của nghiệm chỉnh hóa so với nghiệm chính xác nhỏ như mong muốn. 4 51.2 Hệ Lamé Hệ Lamé được thiết lập từ các dữ kiện vật lý (xem, ví dụ trong [11]) và có liên quan chặt chẽ đến bài toán Elastic. Trong không gian ba chiều với Ω = (0, 1)× (0, 1)× (0, 1) như là vật thể đàn hồi, hệ Lamé được xác định bởi ∂2u ∂t2 + µ∆u+ (λ+ µ)∇ (div (u)) = F, x ∈ Ω, t ∈ (0, T ) , ở đây u (x, t) = (u1 (x, t) , u2 (x, t) , u3 (x, t)) thỏa mãn hệ Lamé, trong đó uj biểu thị cho độ dịch chuyển theo hướng j của vật thể đàn hồi và F (x, t) = (F1 (x, t) , F2 (x, t) , F3 (x, t)) biểu thị lực tác động lên vật thể. Các hằng số λ và µ gọi là các hằng số Lamé. Bài toán thuận là bài toán xác định u từ các dữ kiện đầu u (0, x), ut (0, x) và F . Trong luận văn này ta quan tâm đến bài toán ngược là bài toán xác định F từ các dữ kiện ban đầu. Bài toán này đã được nghiên cứu trên một số dạng của F , chẳng hạn trong [11] các tác giả đã xem xét với F (x, t) = ϕ (t) f (x) và giả sử rằng ϕ ∈ C1 ([0, T ]), ϕ (0) 6= 0, kết hợp với thời gian quan sát T đủ lớn. Mặc dù với các giả sử này bài toán ngược là duy nhất nghiệm tuy nhiên nó vẫn là bài toán không chỉnh vì với một sai số nhỏ trên dữ kiện nhiễu cũng có thể dẫn đến một sai số lớn của nghiệm, do đó sự chỉnh hóa nghiệm là cần thiết. 1.3 Một số kết quả của giải tích thực và giải tích hàm Ký hiệu H1 (Ω) = { u ∈ L2 (Ω) : ∂u ∂xi ∈ L2 (Ω) , i ∈ 1, n } . Với chuẩn: ‖u‖H1(Ω) = ( ‖u‖2L2 + n∑ i=1 ∥∥∥∥ ∂u∂xi ∥∥∥∥2 L2 )1 2 . Định lý 1.3.1 (Công thức Green) Cho Ω là tập mở, bị chận trong Rn có biên Γ là C1 từng khúc. Khi đó nếu u và 6v thuộc H1 (Ω), ta có∫ Ω ∂u ∂xi v dx = − ∫ Ω u ∂v ∂xi dx+ ∫ Γ u v ni dσ, đúng với mọi i ∈ 1, n. Trong đó ni là thành phần thứ i của vectơ pháp tuyến n trên biên Γ của Ω, tích phân cuối được hiểu theo nghĩa vết của u và v. Định lý này có trong [10]. Mệnh đề 1.3.2 Cho E là không gian Hilbert. Giả sử E có cơ sở trực chuẩn đếm được {en}. Khi đó ta có 1. x = ∞∑ i=1 (x, ei) ei, ∀x ∈ E. (chuỗi Fourier) 2. ‖x‖2 = ∞∑ i=1 |(x, ei)|2, ∀x ∈ E. (đẳng thức Parseval) Kết quả này có trong [1]. Mệnh đề 1.3.3 Cho Ω = (0, 1)× (0, 1). Khi đó các hệ: {cos (mpix) cos (npiy)}m,n∈N , {sin (mpix) cos (npiy)}m∈N∗, n∈N là cơ sở trực giao của L2 (Ω). Chứng minh. Ta chứng minh hệ {sin (mpix) cos (npiy)}m∈N∗, n∈N là cơ sở trực giao trong L2 (Ω), hệ còn lại chứng minh tương tự. Kiểm tra trực tiếp thấy các hệ trên là hệ trực giao trong L2 (Ω). Bây giờ giả sử f ∈ L2 (Ω) và ∫ Ω f (x, y) sin (mpix) cos (npiy) dx dy = 0, với mọi m ∈ N∗, n ∈ N . Ta chứng minh f = 0 trong L2 (Ω). Với mỗi n ∈ N , đặt hn (x) = 1∫ 0 f (x, y) cos (npiy) dy. Ta có 1∫ 0 hn (x) sin (mpix) dx = 0, ∀m ∈ N∗, 7và hn ∈ L2 (0, 1) bởi vì |hn (x)|2 ≤  1∫ 0 |f (x, y)|2dy  1∫ 0 cos2 (npiy) dy  ︸ ︷︷ ︸ thuộc L(0,1) . Từ hệ {sin (mpix)}m∈N∗ là cơ sở trực giao của L2 (0, 1) ta có hn = 0 trong L2 (0, 1). Gọi An là tập có độ đo không để hn (x) = 0 với mọi x ∈ [0, 1] \An. Đặt B = ∞⋃ n= 0 An thì B có độ đo không và hn (x) = 0 với mọi x ∈ [0, 1] \B và n ∈ N . Từ hệ {cos (npix)}n∈N là cơ sở trực giao trong L2 (0, 1) suy ra f (x, ·) = 0 trong L2 (0, 1), với mọi x ∈ [0, 1] \B. Vậy f = 0 trong L2 (Ω). 1.4 Một số kết quả của giải tích phức Cho C là trường số phức và hàm số φ : C → C. Ta nói φ là hàm nguyên nếu φ giải tích trên C. Mệnh đề 1.4.1 Cho φ là hàm nguyên và khác hằng, khi đó tồn tại r0 > 0 sao cho Max |z|= r |φ (z)| > 1, đúng với mọi r ≥ r0. Chứng minh. Đặt ψ (r) = Max |z|= r |φ (z)|. Theo nguyên lý môđun cực đại ta có Max |z| ≤ r |φ (z)| = Max |z|= r |φ (z)| , vậy ψ không giảm. Vì φ là hàm nguyên và khác hằng nên không bị chận, từ đó tồn tại z0 ∈ C, z0 6= 0 sao cho |φ (z0)| > 1. Đặt r0 = |z0|, ta được r0 là giá trị cần tìm. Mệnh đề 1.4.2 Với mọi z ∈ C và mọi x ∈ (0, 1), ta có bất đẳng thức |cosh (zx)| ≤ e|z|. 8Chứng minh. Giả sử z = a+ bi, (a, b ∈ R). Ta có |cosh (zx)| = ∣∣∣∣e(a+bi)x + e−(a+bi)x2 ∣∣∣∣ = 12√e2ax + e−2ax + 2 cos (2bx) ≤ 1 2 √ e2ax + e−2ax + 2 = 1 2 ( eax + e−ax ) ≤ e|a|x ≤ e|a| ≤ e|z|, đúng với mọi z ∈ C, mọi x ∈ (0, 1). Mệnh đề 1.4.3 Cho f ∈ L1 (0, 1), đặt F (λ) = 1∫ 0 f (x) cos (λx) dx, λ ∈ C. Khi đó F là hàm nguyên và F ′ (λ) = 1∫ 0 −x f (x) sin (λx) dx. Chứng minh. Với mỗi λ ∈ C cố định. Ta có h (x) = f (x) cos (λx) khả tích Lebesgue trên (0, 1). Thật vậy |f (x) cos (λx)| = ∣∣∣∣∣f (x) ( ∞∑ n=0 (−1)n (λx) 2n (2n)! )∣∣∣∣∣ ≤ |f (x)|  ∞∑ n=0 |λ|2n (2n)!︸ ︷︷ ︸ thuộcR  ∈ L1 (0, 1) . Vậy h ∈ L1 (0, 1). Đặt hk (x) = k∑ n=0 (−1)nλ2n (2n)! x2nf (x), ta có  hk (x)→ f (x) cos (λx) hkn trên (0, 1) , |hk (x)| ≤ |f (x)| ( ∞∑ n= 0 |λ|2n (2n)! ) . 9Áp dụng định lý hội tụ bị chận ta có 1∫ 0 hk (x) dx→ 1∫ 0 f (x) cos (λx) dx. Vậy 1∫ 0 f (x) cos (λx) dx = ∞∑ n= 0 (−1)n (2n)!  1∫ 0 x2nf (x) dx  λ2n. Đặt an = (−1)n (2n)!  1∫ 0 x2n f (x) dx , ta có F (λ) = ∞∑ n=0 an λ 2n, ∀λ ∈ C. Xét chuỗi lũy thừa ∼ F (α) = ∞∑ n= 0 an α n. Ta có ∼ F là hàm nguyên ( vì với mỗi α ∈ C, tồn tại λ ∈ C để α = λ2, mặt khác chuỗi ∞∑ n= 0 an λ 2n hội tụ về F (λ) ) . Mặt khác F (λ) = ∼ F ( g (λ) ) (với g (λ) = λ2) nên F là hàm nguyên. Từ đó F ′ (λ) = ∞∑ n=1 2n anλ 2n−1 = ∞∑ n=1 2n(−1)n (2n)!  1∫ 0 x2nf (x) dx λ2n−1 = ∞∑ n=1 (−1)n (2n− 1)!  1∫ 0 x2n f (x) dx  λ2n−1 = ∞∑ n= 0  1∫ 0 (−1)n+1 x 2n+2 (2n+ 1)! f (x) dx  λ2n+1 = 1∫ 0 −x f (x) ( ∞∑ n= 0 (−1)n (λx) 2n+1 (2n+ 1)! ) dx = 1∫ 0 −x f (x) sin (λx) dx. Mệnh đề 1.4.4 Cho Ω là một tập mở trong C và B (z0, r) ⊂ Ω (B (z0, r) là quả cầu mở tâm z0, bán kính r trong C). Cho f giải tích trên Ω trừ tại một 10 số hữu hạn cực điểm a1, a2, . . . , an nằm trong B (z0, r). Khi đó ta có∫ |z−z0|=r f (z) dz = 2pii ( n∑ k=1 Res [f (z) , ak] ) . Kết quả này được suy trực tiếp từ định lý tích phân Cauchy và định nghĩa thặng dư. Mệnh đề 1.4.5 (Định lý Beurling) Cho f là hàm nguyên. Đặt Mf (r) = Max|z|= r |f (z)| , (r > 0) . Khi đó ta có lim r→∞ sup ln |f (r)| lnMf (r) ≥ −1. Định lý này có trong [13]. Chương 2 Chỉnh hóa nghiệm cho bài toán nhiệt Cho T > 0 và Ω = (0, 1)× (0, 1) là vật dẫn nhiệt. Ta xét bài toán xác định cặp hàm (u, f) thỏa mãn hệ: ut −∆u = ϕ (t) f (x, y) , ux (0, y, t) = ux (1, y, t) = uy (x, 0, t) = uy (x, 1, t) = 0, u (1, y, t) = 0, u (x, y, 0) = g (x, y) , (2.1) với (x, y) ∈ Ω, t ∈ (0, T ), trong đó g ∈ L1 (Ω) và ϕ ∈ L1 (0, T ) được cho. Đây là bài toán dạng truyền nhiệt đặt không chỉnh vì nghiệm bài toán có thể không tồn tại với những dữ kiện được cho hoặc không phụ thuộc liên tục vào dữ kiện này. Ở đây biến (x, y) ∈ Ω gọi là biến không gian, xác định vị trí trên vật dẫn nhiệt Ω. Đại lượng F = ϕ (t) f (x, y) gọi là nguồn nhiệt, nguồn nhiệt phụ thuộc vào thời gian t và vị trí trên vật dẫn nhiệt. Trong phần trình bày này ta không đề cập đến sự tồn tại nghiệm của bài toán, thay vào đó ta sẽ xem xét sự duy nhất nghiệm và chỉnh hóa nghiệm cho bài toán. Mặt khác ở đây ta cũng chỉ tập trung vào chỉnh hóa nghiệm cho hàm f vì khi hàm f được xác định thì bài toán ở hệ (2.1) trở thành bài toán nhiệt thông thường. 11 12 2.1 Một vài ký hiệu • C1 ([0, T ] , L1 (Ω)) = {u : u (·, t) ∈ L1 (Ω) và ut (·, t) ∈ L1 (Ω) , ∀t ∈ [0, T ]}. • A (0, T,H2 (Ω)) = {u : u (·, t) ∈ H2 (Ω) , với hầu hết t ∈ (0, T )}. • Với ω ∈ L1 (Ω), ϕ ∈ L1 (0, T ) và α, β ∈ C. Đặt G (ω) (α, β) = ∫ Ω ω (x, y) cosh (αx) cos (βy) dxdy, D (ϕ) (α, β) = T∫ 0 e− (α 2−β2) t ϕ (t) dt, H (ϕ, ω) (α, β) =  −G (ω) (α, β) D (ϕ) (α, β) nếu D (ϕ) (α, β) 6= 0, 0 nếu D (ϕ) (α, β) = 0. • Từ Mệnh đề 1.3.3, thực hiện chuẩn hóa ta được các hệ{√ k (m,n) cos (mpix) cos (npiy) } m,n∈N và {√ k (m,n) sin (mpix) cos (npiy) } m ∈N∗, n∈N là cơ sở trực chuẩn của L2 (Ω). Trong đó k (m,n) =  1 nếu (m,n) = (0, 0) , 2 nếu m > n = 0 hoặc n > m = 0, 4 nếu m > 0 và n > 0. • Cho A = {x1, x2, . . . , xp} là tập hợp của p số phức phân biệt nhau đôi một, và cho ω là một hàm phức. Khi đó đa thức nội suy Lagrange của ω 13 tại A là: L [A, ω] (z) = p∑ j=1 ∏ k 6=j z − xk xj − xk  ω (xj). 2.2 Tính duy nhất nghiệm Định lý 2.2.1 Cho g ∈ L1 (Ω) và ϕ ∈ L1 (0, T ), trong đó ϕ thỏa mãn thêm giả thiết (H) Tồn tại T0 ∈ (0, T ], θ ≥ 0 và δ > 0 sao cho hoặc ϕ (t) ≥ δ tθ với hầu hết t ∈ (0, T0) hoặc ϕ (t) ≤ −δ tθ với hầu hết t ∈ (0, T0). Thì hệ (2.1) có nhiều nhất một nghiệm (u, f) trong( C1 ( [0, T ] , L1 (Ω) ) ∩ A ( 0, T,H2 (Ω) ) , L2 (Ω) ) . Định lý 2.2.1 có thể được chứng minh nhờ các bổ đề sau: Bổ đề 2.2.2 Giả sử (u, f) ∈ ( C1 ( [0, T ] , L1 (Ω) ) ∩ A ( 0, T,H2 (Ω) ) , L2 (Ω) ) là một nghiệm của hệ (2.1). Thì với mọi (α, n) ∈ C × Z ta có e−(α 2−n2pi2) T ∫ Ω u (x, y, T ) cosh (αx) cos (npiy) dxdy− ∫ Ω g (x, y) cosh (αx) cos (npiy) dxdy = T∫ 0 e−(α 2−n2pi2) t ϕ (t) dt · ∫ Ω f (x, y) cosh (αx) cos (npiy) dxdy. Bổ đề 2.2.3 Cho ω ∈ L1 (Ω) và n ∈ Z, khi đó G (ω) (·, npi) là hàm nguyên và |G (ω) (z, npi)| ≤ e|z| ‖ω‖L1(Ω), ∀z ∈ C. 14 Nếu ω ∈ L2 (Ω) và ω 6= 0 ( trong L2 (Ω) ) thì tồn tại n ∈ N để lim r→∞ sup ln |G (ω) (r, npi)| r ≥ −1. Bổ đề 2.2.4 Cho ϕ ∈ L1 (0, T ) và (α, n) ∈ R× Z, khi đó |D (ϕ) (α, npi)| ≤ ‖ϕ‖L1(0,T ) nếu α2 − n2pi2 > 0. Giả sử ϕ thỏa thêm giả thiết (H) và θ như trong giả thiết (H) của ϕ thì ta có lim (α2−n2pi2)→+∞ ( α2 − n2pi2)θ+1 |D (ϕ) (α, npi)| > 0. Chứng minh các bổ đề. Bổ đề 2.2.2 Chứng minh. Nhân hai vế của hệ (2.1) với W (x, y) = cosh (αx) cos (npiy). Lấy tích phân hai vế trên Ω, áp dụng công thức Green (Định lý 1.3.1) và sử dụng các giả thiết:{ ux (0, y, t) = ux (1, y, t) = uy (x, 0, t) = uy (x, 1, t) = 0, u (1, y, t) = 0, ta được d dt ∫ Ω uW dxdy − (α2 − n2pi2) ∫ Ω uW dxdy = ϕ (t) ∫ Ω f W dxdy. Suy ra d dt e−(α2−n2pi2) t ∫ Ω uW dxdy  = e−(α2−n2pi2) t ϕ (t)∫ Ω f W dxdy. Vì u ∈ C1 ([0, T ] , L1 (Ω)) nên h (t) = d dt e−(α2−n2pi2) t ∫ Ω uW dxdy  liên tục trên [0, T ]. Lấy tích phân hai vế trên [0, T ], thế cận và sử dụng giả thiết 15 u (x, y, 0) = g (x, y) ta có điều phải chứng minh. Bổ đề 2.2.3 Chứng minh. Với n ∈ N , đặt φn (z) = G (ω) (z, npi). Suy ra φn (iz) = ∫ Ω ω (x, y) cos (zx) cos (npiy) dxdy. Vậy z → φn (iz) là biến đổi Fourier cos của hàm x 7→  1∫ 0 ω (x, y) cos (npiy) dy, x ∈ [0, 1] , 0, x ∈ (1,∞) . Theo Mệnh đề 1.4.3 ta có φn (iz) là hàm nguyên, vậy φn là hàm nguyên. Mặt khác |φn (z)| ≤ ∫ Ω |ω (x, y) cosh (zx) cos (npiy)| dxdy ≤ ∫ Ω |ω (x, y) cosh (zx)| dxdy ≤ e|z| ‖ω‖L1(Ω), ∀z ∈ C. Bây giờ với ω ∈ L2 (Ω), ω 6≡ 0, ta có dφn dz (impi) = ∫ Ω ixω (x, y) sin (mpix) cos (npiy) dxdy. Theo Mệnh đề 1.3.3 ta có hệ {sin (mpix) cos (npiy)}m∈N∗, n∈N là cơ sở trực giao trong L2 (Ω), và do (x, y) 7→ ixω (x, y) ∈ L2 (Ω) và không bằng 0 trong L2 (Ω) nên tồn tại n ∈ N để dφn dz 6≡ 0 trên C, vậy φn khác hằng trên C. Bây giờ áp dụng Bổ đề 3.1.4 suy ra lim r→+∞ sup ln |φn (r)| r ≥ −1. Ta có điều phải chứng minh. 16 Bổ đề 2.2.4 Chứng minh. Ta có |D (ϕ) (α, npi)| = ∣∣∣∣∣∣ T∫ 0 e−(α 2−n2pi2) t ϕ (t) dt ∣∣∣∣∣∣ ≤ T∫ 0 e−(α 2−n2pi2) t |ϕ (t)| dt ≤ ‖ϕ‖L1(0,T ), khi α2 − n2pi2 > 0. Giả sử ϕ thỏa thêm giả thiết (H), ta sẽ chứng minh lim λ→+∞ inf λθ+1 ∣∣∣∣∣∣ T∫ 0 e−λtϕ (t) dt ∣∣∣∣∣∣ > 0, từ đó suy ra điều cần chứng minh. Ta có λθ+1 ∣∣∣∣∣∣ T∫ 0 e−λtϕ (t) dt ∣∣∣∣∣∣ ≥ λθ+1 ∣∣∣∣∣∣ T0∫ 0 e−λtϕ (t) dt ∣∣∣∣∣∣− λθ+1 ∣∣∣∣∣∣ T∫ T0 e−λtϕ (t) dt ∣∣∣∣∣∣ ≥ δ λθ+1 T0∫ 0 e−λttθdt− λθ+1e−λT0‖ϕ‖L1(0,T ), ∀λ > 0. Vì lim λ→+∞ ( λθ+1e−λT0 ) = 0 nên ta chỉ cần chứng minh lim λ→+∞ inf ψθ (λ) > 0, trong đó ψθ (λ) = λ θ+1 T0∫ 0 e−λt tθdt. Mặt khác vì ψθ+1 (λ) = −(λT0)θ+1e−λT0 + (θ + 1)ψθ (λ) nên lim λ→+∞ inf ψθ+1 (λ) ≥ (θ + 1) lim λ→+∞ inf ψθ (λ) . Do đó ta chỉ cần chứng minh lim λ→+∞ inf ψθ (λ) > 0 với θ ∈ [0, 1). Kiểm tra trực tiếp ta thấy lim λ→+∞ ψ1 (λ) = 1, lim λ→+∞ ψ2 (λ) = 2. 17 Áp dụng bất đẳng thức Holder ta được ( ψθ (λ) ) 1 2−θ ( ψ2 (λ) ) 1−θ 2−θ ≥ ψ1 (λ) , ( θ ∈ [0, 1) ) . Vậy lim λ→+∞ inf ψθ (λ) ≥ 2θ−1. Ta suy ra điều phải chứng minh. Chứng minh Định lý 2.2.1 Chứng minh. Giả sử (u1, f1) và (u2, f2) là hai nghiệm của hệ (2.1). Đặt u = u1 − u2, f = f1 − f2. Ta chứng minh (u, f) = (0, 0) trong( C1 ( [0, T ] , L1 (Ω) ) ∩ A ( 0, T,H2 (Ω) ) , L2 (Ω) ) . Giả sử f 6= 0 trong L2 (Ω). Theo Bổ đề 2.2.2, với mọi (α, n) ∈ R× Z ta có e−(α 2−n2pi2)T G ( u (·, ·, T ) ) (α, npi) = D (ϕ) (α, npi)G (f) (α, npi) . Từ Bổ đề 2.2.3 thì tồn tại n ∈ N và dãy số dương {αm} sao cho lim m→+∞αm = +∞ và |G (f) (αm, npi)| ≥ e −2αm, ∀m ≥ 1. Ta cũng có ∣∣∣G(u (·, ·, T )) (αm, npi)∣∣∣ ≤ eαm‖u (·, ·, T )‖L1(Ω). Theo Bổ đề 2.2.4 thì tồn tại θ > 0 và C0 > 0 sao cho |D (ϕ) (αm, npi)| ≥ C0 (α2m − n2pi2)θ+1 , với m đủ lớn. Từ các bất đẳng thức trên suy ra với m đủ lớn thì ( α2m − n2pi2 )θ+1 e−(α 2 m−n2pi2) T+3αm‖u (·, ·, T )‖L1(Ω) ≥ C0. Cho m tiến về +∞, vế trái bất đẳng thức trên tiến về 0, ta có mâu thuẫn. Vậy f = 0. Từ đó d dt ( e−(α 2−n2pi2) tG ( u (·, ·, t) ) (α, npi) ) = 0, 18 kết hợp với u (·, ·, 0) = 0 suy ra G ( u (·, ·, t) ) (α, npi) = 0, ∀t ∈ [0, T ] , ∀ (α, n) ∈ C × Z. Bây giờ lấy α = impi, m ∈ N ta được∫ Ω u (x, y, t) cos (mpix) cos (npiy) dxdy = 0, ∀t ∈ [0, T ] ; ∀m,n ∈ N. Vì u (·, ·, t) ∈ L2 (Ω) nên từ Mệnh đề 1.3.3 suy ra u (·, ·, t) = 0 trong L2 (Ω), ∀t ∈ [0, T ]. Vậy u = 0. Định lý được chứng minh. 2.3 Chỉnh hóa nghiệm Định lý 2.3.1 Cho ϕ0 ∈ L1 (0, T ) thỏa mãn giả thiết (H) và g0 ∈ L1 (Ω). Giả sử rằng (u0, f0) ∈ ( C1 ( [0, T ] , L1 (Ω) ) ∩ A ( 0, T,H2 (Ω) ) , L2 (Ω) ) là nghiệm chính xác của hệ (2.1) đối với dữ kiện ϕ0, g0. Cho ε ∈ (0, 1) và ϕε ∈ L1 (0, T ), gε ∈ L1 (Ω) sao cho ‖ϕε − ϕ0‖L1(0,T ) ≤ ε, ‖gε − g0‖L1(Ω) ≤ ε. Nghiệm chỉnh hóa fε được xây dựng từ các dữ kiện ϕε, gε như sau rε ∈ Z ∩ [ ln ( ε−1 ) 50 , ln ( ε−1 ) 50 + 1 ) , B (rε) = {± (4rε + j) : j = 1, 2, . . . , 20rε} , Fε (m,n) = L [B (rε) , H (ϕε, gε) (·, npi)] (impi) , fε (x, y) = ∑ 0≤m,n≤ rε k (m,n)Fε (m,n) cos (mpix) cos (npiy) . 19 Thì fε ∈ C∞ ( R2 ) và lim ε→0+ fε = f0 trong L 2 (Ω). Hơn nữa nếu f0 ∈ H1 (Ω) thì lim ε→0+ fε = f0 trong H 1 (Ω) và tồn tại ε0 > 0 chỉ phụ thuộc vào ϕ0, g0 sao cho ‖fε − f0‖L2(Ω) ≤ 50 pi ln (ε−1) ‖f0‖H1(Ω), ∀ε ∈ (0, ε0) . Định lý 2.3.1 có thể được chứng minh nhờ các bổ đề sau Bổ đề 2.3.2 Cho r ∈ Z, r ≥ 55 và Br = {± (4r + j) : j = 1, 2, . . . , 20r}. Cho ω và ∼ ω là hai hàm phức chẵn sao cho ω là hàm nguyên và |ω (z)| ≤ Ae|z|, ∀z ∈ C (A độc lập với z). Khi đó sup |z| ≤pir ∣∣∣ω (z)− L [Br, ∼ω] (z)∣∣∣ ≤ Ae−r + 20 r e25r sup z∈Br ∣∣∣ω (z)− ∼ω (z)∣∣∣ . Bổ đề 2.3.3 Cho u0, f0, ϕ0, g0, ϕε, gε, rε, B (rε) định nghĩa như trong Định lý 2.3.1 và θ trong điều kiện (H) đối với ϕ0. Khi đó tồn tại ε1 > 0 phụ thuộc chỉ vào ϕ0 và g0 sao cho |G (f0) (α, npi)−H (ϕε, gε) (α, npi)| ≤ ( ln ( ε−1 ))4θ+5 e|α| ε, đúng với mọi ε ∈ (0, ε1), n ∈ [0, rε] ∩N , α ∈ B (rε). Bổ đề 2.3.4 Cho ω ∈ H1 (Ω) và với mỗi số nguyên M ≥ 1 đặt ΓM (ω) (x, y) = ∑ 0≤m,n≤M k (m,n)G (ω) (impi, npi) cos (mpix) cos (npiy) . Khi đó lim M→+∞ ΓM (ω) = ω trong H 1 (Ω) và ‖ΓM (ω)− ω‖L2(Ω) ≤ 1 pi (M + 1) ‖ω‖H1(Ω). Chứng minh các bổ đề. Bổ đề 2.3.2 20 Chứng minh. Với z ∈ C và |z| ≤ pir cố định. Ký hiệu zj = 4r + j, ( j ∈ 1, 20r) ; γ = {z ∈ C : |z| = 45r} . Áp dụng Mệnh đề 1.4.4, tính thặng dư tại 41 cực điểm đơn ta được ∫ γ ω (ξ) ξ − z 20r∏ j=1 z2 − z2j ξ2 − z2j dξ = 2pii ( ω (z)− L [Br, ω] (z) ) . Suy ra |ω (z)− L [Br, ω] (z)| ≤ 45r sup ξ∈γ { |ω (ξ)| |ξ − z| 20r∏ j=1 ∣∣z2 − z2j ∣∣∣∣ξ2 − z2j ∣∣ } . (2.2) Với ξ ∈ γ, ta có |ω (ξ)| ≤ Ae45r, |ξ − z| ≥ (45− pi) r, và 20r∏ j=1 ∣∣z2 − z2j ∣∣∣∣ξ2 − z2j ∣∣ ≤ 20r∏ j=1 |z|2 + z2j |ξ|2 − z2j ≤ 20r∏ j=1 (pir) 2 + z2j (45r) 2 − z2j . (2.3) Mặt khác 20r∏ j=1 (pir) 2 + z2j (45r) 2 − z2j ≤ 45− pi 45 e−46r, ∀r ≥ 55. (2.4) Thật vậy, kiểm tra trực tiếp ta có hàm v (x) = ln ( (pir) 2 + x (45r) 2 − x ) là tăng và lõm trong [ 0, (24r) 2 ] . Áp dụng bất đẳng thức Jensen ta được 20r∑ j=1 v ( z2j ) = 5∑ k=1  4kr∑ j=4(k−1)r+1 v ( z2j ) ≤ 4r 5∑ k=1 v  1 4r 4kr∑ j=4(k−1)r+1 z2j  ≤ 4r 5∑ k=1 v (( 16k2 + 16k + 16 3 ) r2 + (4k + 2) r + 1 6 ) ≤ 4r 5∑ k=1 v (( 16k2 + 16k + 16 3 ) r2 + (4k + 2) r2 55 + r2 6× 552 ) 21 = 4r 5∑ k=1 ln  pi 2 + ( 16k2 + 884 55 k + 32487 6050 ) 452 − ( 16k2 + 884 55 k + 32487 6050 )  < −46r + ln ( 45− pi 45 ) , ∀r ≥ 55. Vậy (2.4) đúng. Từ (2.2), (2.3) và (2.4) ta có |ω (z)− L [Br, ω] (z)| ≤ Ae−r, ∀r ≥ 55. (2.5) Bây giờ kiểm tra trực tiếp với các hàm chẵn ω và ∼ ω ta có L [ Br, ω − ∼ω ] (z) = 20r∑ j=1 ∏ k 6=j z2 − z2k z2j − z2k (ω (zj)− ∼ω (zj)). Suy ra ∣∣∣L [Br, ω − ∼ω] (z)∣∣∣ ≤ 20rσ sup 1≤j≤ 20r ∏ k 6=j |z|2 + z2k∣∣z2j − z2k∣∣ , (2.6) trong đó σ = sup 1≤j≤ 20r ∣∣∣ω (zj)− ∼ω (zj)∣∣∣. Với mỗi j ∈ 1, ..., 20r ta có ∏ k 6=j |z|2 + z2k∣∣z2j − z2k∣∣ = ∏ k 6=j |z|2 + z2k (zj + zk) zk zk |zj − zk| ≤ ∏ k 6=j zk |zj − zk| ≤ 20r∏ k=2 zk∏ k 6=j |zj − zk| = (4r + 2) (4r + 3) . . . (24r) (j − 1)! (20r − j)! ≤ (4r + 2) (4r + 3) . . . (24r) (10r − 1)! (10r)! := J (r) . (2.7) 22 Ta có J (1) < e25 và kiểm tra trực tiếp thấy J (r + 1) J (r) = (24r + 1) (24r + 2) . . . (24r + 24) (4r + 2) . . . (4r + 5)× [(10r + 1) . . . (10r + 9)]2 × (10r) (10r + 10) < e25, ∀r ≥ 1. Vậy J (r) < e25r, ∀r ≥ 1. (2.8) Từ (2.6), (2.7) và (2.8) suy ra∣∣∣L [Br, ω − ∼ω] (z)∣∣∣ ≤ 20rσ e25r, ∀r ≥ 55. (2.9) Bây giờ từ (2.5) và (2.9), kết hợp với bất đẳng thức tam giác∣∣∣ω (z)− L [Br, ∼ω] (z)∣∣∣ ≤ ∣∣∣ω (z)− L [Br, ω] (z)∣∣∣+ ∣∣∣L [Br, ω − ∼ω] (z)∣∣∣ , ta có bất đẳng thức cần chứng minh. Bổ đề 2.3.3 Chứng minh. Với mọi α ∈ B (rε), ta có 4rε ≤ |α| ≤ 24rε và lim ε→0+ ( 1− ln ε 50 )2 |ln ε|2(θ+1) = 0, suy ra tồn tại ε∗1 ∈ (0, 1) và C (ϕ0) > 0 chỉ phụ thuộc vào ϕ0 sao cho α2 − n2pi2 ≥ ln ( ε−1 ) T , |D (ϕ0) (α, npi)| ≥ C (ϕ0) (α2 − n2pi2)θ+1 ≥ ( ln ( ε−1 ))−2(θ+1) , (2.10) 23 đúng với mọi ε ∈ (0, ε∗1), n ∈ N ∩ [0, rε], α ∈ B (rε). Do đó |D (ϕε) (α, npi)| ≥ |D (ϕ0) (α, npi)| − |D (ϕ0) (α, npi)−D (ϕε) (α, npi)| ≥ ( ln ( ε−1 ))−2(θ+1) − ε ≥ 1 2 ( ln ( ε−1 ))−2(θ+1) > 0. (2.11) Từ bất đẳng thức (2.10), kết hợp với Bổ đề 2.2.2 và 2.2.3, suy ra tồn tại ε∗2 < ε ∗ 1 để với mọi ε ∈ (0, ε∗2), n ∈ N ∩ [0, rε], α ∈ B (rε) ta có |G (f0) (α, npi)−H (ϕ0, g0) (α, npi)| = e−(α 2−n2pi2)T ∣∣∣∣∣∣ G ( u0 (·, ·, T ) ) (α, npi) D (ϕ0) (α, npi) ∣∣∣∣∣∣ ≤ e−(α2−n2pi2)T e |α|‖u0 (·, ·, T )‖L1(Ω)( ln (ε−1) )−2(θ+1) ≤ 1 2 ε e|α| ( ln ( ε−1 ))4θ+5 . (2.12) Bây giờ sử dụng các bất đẳng thức (2.10), (2.11) và Bổ đề 2.2.3, suy ra tồn tại ε∗3 (0 < ε∗3 < ε ∗ 2) chỉ phụ thuộc vào ϕ0 và g0 sao cho với mọi ε ∈ (0, ε∗3), α ∈ B (rε), n ∈ N ∩ [0, rε] ta có |H (ϕ0, g0) (α, npi)−H (ϕε, gε) (α, npi)| = ∣∣∣∣G (g0) (α, npi)D (ϕ0) (α, npi) − G (gε) (α, npi)D (ϕε) (α, npi) ∣∣∣∣ ≤ |G (g0) (α, npi)| |D (ϕε) (α, npi)−D (ϕ0) (α, npi)||D (ϕ0) (α, npi)| |D (ϕε) (α, npi)| + |D (ϕ0) (α, npi)| |G (gε) (α, npi)−G (g0) (α, npi)| |D (ϕ0) (α, npi)| |D (ϕε) (α, npi)| 24 ≤ 2e |α| ε‖g0‖L1(Ω)( ln (ε−1) )−2(θ+1)( ln (ε−1) )−2(θ+1) + 2e|α| ε‖ϕ0‖L1(0,T )( ln (ε−1) )−2(θ+1)( ln (ε−1) )−2(θ+1) ≤ 1 2 ( ln ( ε−1 ))4θ+5 ε e|α|. (2.13) Từ các bất đẳng thức (2.12) và (2.13) ta có điều phải chứng minh. Bổ đề 2.3.4 Chứng minh. Áp dụng công thức Green (Định lý 1.3.1) ta được pimG (ω) (impi, npi) = − ∫ Ω ωx (x, y) sin (mpix) cos (npiy) dxdy, pinG (ω) (impi, npi) = − ∫ Ω ωy (x, y) cos (mpix) sin (npiy) dxdy. Suy ra pi2 ( m2 + n2 ) k (m,n) |G (ω) (impi, npi)|2 = k (m,n) ∣∣∣∣∣∣ ∫ Ω ωx (x, y) sin (mpix) cos (npiy) dxdy ∣∣∣∣∣∣ 2 + k (m,n) ∣∣∣∣∣∣ ∫ Ω ωy (x, y) cos (mpix) sin (npiy) dxdy ∣∣∣∣∣∣ 2 . Áp dụng đẳng thức Parseval (Mệnh đề 1.3.2) ta được ∑ m,n≥ 0 k (m,n) ∣∣∣∣∣∣ ∫ Ω ωx (x, y) sin (mpix) cos (npiy) dxdy ∣∣∣∣∣∣ 2 = ‖ωx‖2L2(Ω) , ∑ m,n≥ 0 k (m,n) ∣∣∣∣∣∣ ∫ Ω ωy (x, y) cos (mpix) sin (npiy) dxdy ∣∣∣∣∣∣ 2 = ‖ωy‖2L2(Ω) . 25 Suy ra∑ m,n≥0 pi2 ( m2 + n2 ) k (m,n) |G (ω) (impi, npi)|2 = ‖ωx‖2L2(Ω) + ‖ωy‖2L2(Ω) ≤ ‖ω‖2H1(Ω) . Dùng đẳng thức Parseval cho cơ sở trực chuẩn{√ k (m,n) cos (mpix) cos (npiy) } m,n∈N ta được ‖ΓM (ω)− ω‖2L2(Ω) = ∑ Max{m,n} >M k (m,n) |G (ω) (impi, npi)|2 ≤ 1 pi2(M + 1) 2 × × ∑ Max{m,n} >M pi2 ( m2 + n2 ) k (m,n) |G (ω) (impi, npi)|2 ≤ 1 pi2(M + 1) 2 ‖ω‖2H1(Ω) . Cuối cùng sử dụng đẳng thức Parseval lần nữa ta được ‖ΓM (ω)− ω‖2H1 = ∑ Max{m,n} >M ( 1 + pi2 ( m2 + n2 )) k (m,n) |G (ω) (impi, npi)|2. Vì chuỗi ở vế phải hội tụ với mỗi M nên cho M tiến về vô cùng ta được ‖ΓM (ω)− ω‖H1(Ω) → 0. Chứng minh Định lý 2.3.1 26 Chứng minh. Tính toán trực tiếp ta được G (fε) (impi, npi) =  Fε (m,n) với m,n ∈ [0, rε] , 0 với m,n /∈ [0, rε] . Với m, n ∈ [0, rε], từ Bổ đề 2.2.3 ta có G (f0) (·, npi) là hàm nguyên chẵn và |G (f0) (z, npi)| ≤ ‖f0‖L1(Ω) e|z|, ∀z ∈ C. Từ Bổ đề 2.3.3 suy ra tồn tại ε10 > 0 để với mọi ε ∈ ( 0, ε10 ) ta có rε ≥ 55 và sup α∈B(rε) |G (f0) (α, npi)−H (ϕε, gε) (α, npi)| ≤ ( ln ( ε−1 ))4θ+5 e24rεε. Áp dụng Bổ đề 2.3.2 ta được |G (f0) (impi, npi)−G (fε) (impi, npi)| = |G (f0) (impi, npi)− L [B (rε) , H (ϕε, gε) (·, npi)] (impi)| ≤ ‖f0‖L1(Ω) e−rε + 20 rε e25rε ( ln ( ε−1 ))4θ+5 e24rε ε, ∀m,n ∈ [0, rε] . Từ bất đẳng thức cuối cùng suy ra tồn tại ε20 > 0, ε 2 0 < ε 1 0 sao cho ∀ε ∈ ( 0, ε20 ) ta có |G (f0) (impi, npi)−G (fε) (impi, npi)| ≤ ε 151 , ∀m,n ∈ [0, rε] . Từ đó ‖Γrε (f0)− fε‖2H1(Ω) = ∑ 0 ≤m,n ≤ rε ( 1 + pi2 ( m2 + n2 )) k (m,n) |G (f0) (impi, npi)−G (fε) (impi, npi)|2 ≤ (1 + rε)2 ( 1 + 2pi2r2ε )× 4ε 251 . 27 Từ bất đẳng thức cuối cùng suy ra tồn tại ε30 ∈ ( 0, ε20 ) sao cho ∀ε ∈ (0, ε30) ta được ‖Γrε (f0)− fε‖2H1(Ω) ≤ ε 1 26 . (2.14) Vậy lim ε→0+ (Γrε (f0)− fε) = 0 trong H1 (Ω). Mặt khác từ công thức chuỗi Fourier (Mệnh đề 1.3.2) ta được lim ε→0+ Γrε (f0) = f0 trong L 2 (Ω). Vậy lim ε→0+ fε = f0 trong L2 (Ω). Nếu f0 ∈ H1 (Ω), áp dụng Bổ đề 2.3.4 ta được lim ε→0+ Γrε (f0) = f0 trong H1 (Ω). Suy ra lim ε→0+ fε = f0 trong H 1 (Ω). Mặt khác từ (2.14) và Bổ đề 2.3.4 ta có ‖fε − f0‖L2(Ω) ≤ ‖Γrε (f0)− fε‖L2(Ω) + ‖Γrε (f0)− f0‖L2(Ω) ≤ ε 152 + 1 pi (rε + 1) ‖f0‖H1(Ω). Từ bất đẳng thức cuối suy ra tồn tại ε0 ∈ ( 0, ε30 ) sao cho ∀ε ∈ (0, ε0) ta có ‖fε − f0‖L2(Ω) ≤ 50 pi ln (ε−1) ‖f0‖H1(Ω). 28 2.4 Ví dụ áp dụng Sự chỉnh hóa là cần thiết, chẳng hạn xét hệ (2.1) với ϕ0 = pi 2e−4pi 2t (thỏa giả thiết (H)) và g0 = ( 1 + cos (pix) ) cos (piy). Khi đó hệ (2.1) có nghiệm duy nhất u0 (x, y, t) = e −4pi2t ( 1 + cos (pix) ) cos (piy) , f0 = −3 cos (piy)− 2 cos (pix) cos (piy) . Với mỗi số nguyên dương k ta xét các dữ kiện khá gần ϕ0 và g0 như sau ϕk = ϕ0, gk (x, y) = g0 (x, y) + pi k sin2 (kpix) cos (kpiy) . Với các dữ kiện ϕk và gk hệ (2.1) có nghiệm duy nhất u∗k (x, y, t) = u0 (x, y, t) + pi k e−4pi 2t sin2 (kpix) cos (kpiy) , f∗k (x, y) = f0 (x, y) + pi k (( 5k2 − 4) sin2 (kpix)− 2k2) cos (kpiy) . Tuy nhiên ta thấy ‖f∗k − f0‖L2(Ω) = pi 4 √ 27k2 − 56 + 48 k2 → +∞ khi k → +∞, tức là với sai số càng nhỏ của dữ kiện lại dẫn đến sai số càng lớn trên nghiệm. Bây giờ sử dụng kết quả về chỉnh hóa nghiệm ta xây dựng được các nghiệm chỉnh hóa, chẳng hạn với ε = 1 100 ta có nghiệm chỉnh hóa f 1 100 (x, y) ≈ −2.999721 cos (piy)− 1.997145 cos (pix) cos (piy) . Sai số so với nghiệm chính xác là∥∥∥f 1 100 − f0 ∥∥∥ L2(Ω) ≈ 0.001441. Chương 3 Chỉnh hóa nghiệm cho bài toán Elastic Xét bài toán Elastic Tìm nghiệm u thỏa mãn hệ Lamé và hệ điều kiện như sau: ∂2u ∂t2 + µ∆u+ (λ+ µ)∇(div(u)) = ϕ(f1, f2, f3), trong Ω× (0, T ) , (3.1) (u1(x, t), u2(x, t), u3(x, t)) = (0, 0, 0), (x, t) ∈ ∂Ω× (0, T ) , (3.2) (u1 (x, 0) , u2 (x, 0) , u3 (x, 0)) = (g1(x), g2(x), g3(x)), x ∈ Ω, (3.3)( ∂u1 ∂t (x, 0) , ∂u2 ∂t (x, 0) , ∂u3 ∂t (x, 0) ) = (h1(x), h2(x), h3(x)), x ∈ Ω, (3.4)  σ1 τ12 τ13τ21 σ2 τ23 τ31 τ32 σ3   n1n2 n3  =  X1X2 X3  . (3.5) Trong đó: σj = λdiv (u) + 2µ ∂uj ∂xj , τjk = µ( ∂uj ∂xk + ∂uk ∂xj ), µ < 0, λ + 2µ < 0, Ω = 29 30 (0, 1)× (0, 1)× (0, 1), n = (n1, n2, n3) là pháp vectơ trên ∂Ω và dữ kiện I (ϕ,X, g, h) ∈ ( L1 (0, T ) , ( L1 ( 0, T, L1 (∂Ω) ))3 , ( L2 (Ω) )3 , ( L2 (Ω) )3) . Ta xét các giả thiết sau: (W1) Tồn tại Λ (ϕ) ∈ (0, T ) và C (ϕ) > 0 sao cho: hoặc ._.

Các file đính kèm theo tài liệu này:

  • pdfLA5403.PDF
Tài liệu liên quan