KHOÁ LUẬN TỐT NGHIỆP ĐẠI HỌC SƯ PHẠM NGÀNH VẬT LÝ " SỬ DỤNG PHƯƠNG PHÁP HÀM GREEN ĐỂ GIẢI MỘT SỐ BÀI TOÁN TRUYỀN NHIỆT " ppsx

TRƯỜNG ĐẠI HỌC AN GIANG KHOA SƯ PHẠM BỘ MÔN VẬT LÝ ---oOo--- VÕ THỊ CẨM LOAN Lớp DH5L KHOÁ LUẬN TỐT NGHIỆP ĐẠI HỌC SƯ PHẠM NGÀNH VẬT LÝ SỬ DỤNG PHƯƠNG PHÁP HÀM GREEN ĐỂ GIẢI MỘT SỐ

TRƯỜNG ĐẠI HỌC AN GIANG

KHOA SƯ PHẠM

BỘ MÔN VẬT LÝ -oOo -

VÕ THỊ CẨM LOAN Lớp DH5L

KHOÁ LUẬN TỐT NGHIỆP ĐẠI HỌC SƯ PHẠM

NGÀNH VẬT LÝ

SỬ DỤNG PHƯƠNG PHÁP HÀM GREEN ĐỂ GIẢI MỘT SỐ BÀI TOÁN

TRUYỀN NHIỆT

Giảng viên hướng dẫn: Th.S HỒ XUÂN HUY

Long Xuyên, 5-2008

LỜI CẢM ƠN

********

Trước tiên cho tôi được gởi lời cảm ơn chân thành nhất tới Ban Giám Hiệu Trường Đại Học An Giang, Ban Chủ Nhiệm Khoa sư Phạm, hội Đồng Khoa Học và Đào Tạo Khoa sư Phạm Trường Đại Học An Giang, đã tạo điều kiện để tôi được làm khoá luận, đã quan tâm và đôn đốc tôi trong quá trình làm khoá luận này

Xin cảm ơn sâu sắc nhất tới các thầy cô trong Tổ

Bộ Môn Vật Lý Đặt biệt là giáo viên hướng dẫn Thạc

Sĩ Hồ Xuân Huy đã tận tình hướng dẫn và giúp đỡ tôi hoàn thành khoá luận này

MỤC LỤC PHẦN I: MỞ ĐẦU Trang 1

1 Lý do chọn đề tài Trang 1

2 Mục đích nghiên cứu Trang 1

3 Đối tượng nghiên cứu Trang 1

4 Nhiệm vụ nghiên cứu Trang 1

5 Phương pháp nghiên cứu Trang 2

6 Giả thuyết khoa học Trang 2

7 Phạm vi nghiên cứu Trang 2

8 Đóng góp của khóa luận Trang 2

9 Cấu trúc khóa luận Trang 2

PHẦN II: NỘI DUNG Trang 3

CHƯƠNG I: CƠ SỞ LÝ LUẬN CỦA ĐỀ TÀI Trang 3 1.1 Lý luận về bài tập vật lý Trang 3 1.2 Bài toán biên Trang 6 1.3 Khái niệm toán tử, hàm riêng, trị riêng Trang 8 1.4 Phương pháp tách biến Trang 11 1.5 Phương pháp biến thiên tham số Trang 15 CHƯƠNG II: XÂY DỰNG PHƯƠNG PHÁP HÀM GREEN Trang 17 2.1 Khái niệm hàm Green, tính đối xứng của hàm Green Trang 17 2.2 Xây dựng phương pháp hàm Green Trang 20 2.3 Hàm riêng, trị riêng cho hàm Green Trang 21 2.4 Hàm điều hòa Biễu diễn Green Trang 23 CHƯƠNG III: SỬ DỤNG PHƯƠNG PHÁP HÀM GREEN ĐỂ

GIẢI MỘT SỐ BÀI TOÁN TRUYẾN NHIỆT Trang 27 3.1 Thiết lập phương trình truyền nhiệt Trang 27 3.2 Bài toán biên phụ thuộc thời gian Trang 30 3.2.1 Phương pháp tách biến Fourier cho bài toán truyền nhiệt Trang 30 3.2.2 Phương pháp hàm Green cho bài toán truyền nhiệt Trang 33 3.2.3 Bài toán truyền nhiệt trong miền tròn Trang 35 3.3 Bài toán biên truyền nhiệt dừng Trang 38

PHẦN III: KẾT LUẬN Trang 45

PHỤ LỤC 1 Trang 46 PHỤ LỤC 2 Trang 48

PHẦN I: MỞ ĐẦU

1 Lý do chọn đề tài

Phương pháp toán lý là một học phần rất quan trọng trong chương trình đào tạo giáo viên THPT Giúp cho sinh viên làm quen dần với phương pháp toán học hiện đại trong vật lý, hiểu rõ hơn bản chất của quá trình truyền sóng và quá trình truyền nhiệt trong vật chất Học phần này có liên quan đến nhiều môn học khác: điện và từ, điện động lực, nhiệt động lực, vật lý thống kê, cơ học lượng tử,… Việc nghiên cứu học phần này là cơ sở nghiên cứu các môn học khác Vì thế việc nghiên cứu nó gặp nhiều khó khăn Bên cạnh đó học phần này có nhiều dạng bài tập, mỗi dạng lại có nhiều phương pháp giải đòi hỏi sinh viên phải lựa chọn phương pháp giải phù hợp với mỗi dạng

Cụ thể là bài tập phần truyền nhiệt có các phương pháp giải như: phương pháp tách biến Fourier, phương pháp biến đổi Laplace, phương pháp hàm Green, hàm Bessel Mỗi phương pháp đều có ưu điểm và hạn chế riêng

Đối với một số dạng bài tập nhiều chiều, khi giải bằng phương pháp biến đổi Fourier, phương trình Laplace, thì việc tìm nghiệm gặp khó khăn và giải rất phức tạp, trong khi đó nếu dùng phương pháp hàm Green thì việc tìm nghiệm của bài toán là đơn giản hơn nhiều, phương pháp hàm Green là phương pháp không giải trực tiếp phương trình vi phân mà tìm hàm Green thông qua việc giải phương trình khác để tìm hàm Green Rồi biểu diễn nghiệm cần tìm thông qua hàm Green Phương pháp hàm Green là một phương pháp khó, tuy nhiên nó lại được áp dụng hiệu quả vào việc giải các bài toán biên nhiều chiều Nhưng các sách lý thuyết thường không đề cặp đến phương pháp này, hoặc đề cặp quá ít, làm cho sinh viên gặp khó khăn trong việc áp dụng phương pháp này vào bài tập Yêu cầu bổ sung một phương pháp giải hiệu quả cho bài toán truyền nhiệt

là rất cần thiết

Với những lý do trên chúng tôi chọn đề tài : “Sử dụng phương pháp hàm Green

để giải một số bài toán truyền nhiệt”

2 Mục đích nghiên cứu

• Tìm hiểu các bài toán truyền nhiệt

• Cơ sở toán học cho phương pháp hàm Green

• Dùng phương pháp hàm Green để tìm nghiệm của bài toán truyền nhiệt

3 Đối tượng nghiên cứu

• Cơ sở lý luận về bài tập vật lý

• Cơ sở toán học cho phương pháp hàm Green

• Các bài tập truyền nhiệt

4 Nhiệm vụ nghiên cứu

• Nghiên cứu cơ sở toán học cho việc xây dựng hàm Green

• Xây dựng phương pháp hàm Green để tìm nghiệm của bài toán truyền nhiệt

• Giải một số bài toán truyền nhiệt bằng phương pháp hàm Green

5 Phương pháp nghiên cứu

• Đọc sách và tham khảo tài liệu

• Phương pháp toán học

• Phương pháp phân tích

• Phương pháp đàm thoại trao đổi ý kiến với giáo viên

6 Giả thuyết khoa học

Nếu dùng phương pháp hàm Green thì có thể tìm được nghiệm của bài toán truyền nhiệt

7 Phạm vi nghiên cứu

• Các bài toán truyền nhiệt ứng với các điều kiện biên

8 Đóng góp của khóa luận

• Có thể làm tài liệu tham khảo cho sinh viên

• Góp phần nâng cao kết quả học tập học phần phương pháp toán

lý cho sinh viên

9 Cấu trúc của khoá luận: gồm

Phần I: Mở đầu

Phần II : Nội dung nghiên cứu

Chương I: Cơ sở lý luận của đề tài

Chương II: Xây dựng phương pháp hàm Green

Chương III: Sử dụng phương pháp hàm Green để giải một số bài

toán truyền nhiệt

Phần III: Kết luận

PHẦN II: NỘI DUNG CHƯƠNG I: CƠ SỞ LÝ LUẬN CỦA ĐỀ TÀI 1.1 LÝ LUẬN VỀ BÀI TẬP VẬT LÝ:

1.1.1 Khái niệm về bài tập vật lý

Bài tập vật lý là một yêu cầu đặt ra cho người học, được người học giải quyết dựa trên cơ sở các lập luận logic, nhờ các phép tính toán, các thí nghiệm, dựa trên các kiến thức về khái niệm, định luật và các thuyết vật lý

1.1.2 Vai trò và tác dụng của bài tập vật lý

Xét về mặt phát triển tính tự lực của người học và nhất là rèn luyện kỷ năng vận dụng kiến thức đã lĩnh hội được thì vai trò của bài tập vật lý trong quá trình học tập có một giá trị rất lớn Bài tập vật lý được sử dụng ở nhiều khâu trong quá trình dạy học

- Bài tập là một phương tiện nghiên cứu hiện tượng vật lý Trong quá trình dạy học vật lý người học được làm quen với bản chất của các hiện tượng vật

lý bằng nhiều cách khác nhau như: Kể chuyện, biểu diễn thí nghiệm, làm bài thí nghiệm, tiến hành tham quan Ở đây tính tích cực của người học và do đó chiều sâu và độ vững chắc của kiến thức sẽ lớn nhất khi “tình huống có vấn đề” được tạo ra, trong nhiều trường hợp nhờ tình huống này có thể làm xuất hiện một kiểu bài tập mà trong quá trình giải người học sẽ phát hiện lại quy luật vật lý chứ không phải tiếp thu quy luật dưới hình thức có sẵn

- Bài tập là một phương tiện hình thành các khái niệm Bằng cách dựa vào các kiến thức hiện có của người học, trong quá trình làm bài tập, ta có thể cho người học phân tích các hiện tượng vật lý đang được nghiên cứu, hình thành các khái niệm về các hiện tượng vật lý và các đại lượng vật lý

- Bài tập là một phương tiện phát triển tư duy vật lý cho người học việc giải bài tập làm phát triển tư duy logic, sự nhanh trí Trong quá trình tư duy

có sự phân tích và tổng hợp mối liên hệ giữa các hiện tượng, các đại lượng vật lý đặc trưng cho chúng

- Bài tập là một phương tiện rèn luyện kỷ năng vận dụng các kiến thức của người học vào thực tiễn Đối với việc giáo dục kỹ thuật tổng hợp bài tập vật

lý có ý nghĩa rất lớn Những bài tập này là một trong những phương tiện thuận lợi để người học liên hệ lý thuyết với thực hành, học tập với đời sống Nội dung của bài tập phải đảm bảo các yêu cầu sau:

+ Nội dung của bài tập phải gắn với tài liệu thuộc chương trình đang học + Hiện tượng đang được nghiên cứu phải được áp dụng phổ biến trong thực tiển

+ Bài tập đưa ra phải là những vấn đề gần với thực tế

+ Không những nội dung mà còn hình thức của bài tập cũng phải gắn với các điều kiện thường gặp trong cuộc sống Trong các bài tập đó không có sẵn dữ kiện mà phải tìm dữ kiện cần thiết ở các sơ đồ, bản vẽ kỷ thuật, ở các sách báo tra cứu hoặc từ thí nghiệm

- Bài tập về hiện tượng vật lý trong sinh hoạt hằng ngày cũng có một ý nghĩa to lớn Chúng giúp cho người học nhìn thấy khoa học vật lý xung quanh chúng ta, bồi dưỡng cho người học khả năng quan sát Với các bài tập này, trong quá trình giải, người học sẽ có được kỷ năng, kỷ xảo vận dụng các kiến thức của mình để phân tích các hiện tượng vật lý khác nhau trong

tự nhiên, trong kỷ thuật và trong đời sống, đặc biệt có những bài tập khi giải đòi hỏi người học phải sử dụng kinh nghiệm trong lao động, sinh hoạt và sử dụng những kết quả quan sát thực tế hằng ngày

- Bài tập vật lý là một phương tiện để giáo dục người học, nhờ bài tập vật

lý ta có thể giới thiệu cho người học biết sự xuất hiện những tư tưởng, quan điểm tiên tiến, hiện đại, những phát minh, những thành tựu của nền khoa học trong và ngoài nước Tác dụng giáo dục của bài tập vật lý còn thể hiện ở chổ: chúng là phương tiện hiệu quả để rèn luyện đức tính kiên trì, vượt khó,

ý chí và nhân cách của người học Việc giải bài tập vật lý có thể mang đến cho người học niềm phấn khởi sáng tạo, tăng thêm sự yêu thích bộ môn, tăng cường hứng thú học tập

- Bài tập vật lý cũng là phương tiện kiểm tra mức độ nắm vững kiến thức

và kỷ năng, kỷ xảo của người học Đồng thời nó cũng là công cụ giúp người học ôn tập, đào sâu, mở rộng kiến thức

1.1.3 Cơ sở định hướng giải bài tập vật lý

1.1.3.1 Hoạt động giải bài tập vật lý

− Mục tiêu cần đạt tới khi giải một bài toán vật lý là tìm được câu trả lời đúng đắn, giải đáp được vấn đề đặt ra một cách có căn cứ khoa học chặt chẽ Quá trình giải một bài toán thực chất là tìm hiểu điều kiện của bài toán, xem xét hiện tượng vật lý được đề cập và dựa trên kiến thức vật lý toán để nghĩ tới mối liên hệ có thể có của cái đã cho và cái cần tìm sao cho có thể thấy được cái phải tìm có mối liên hệ trực tiếp hoặc gián tiếp với cái đã cho, từ đó đi đến chỉ rõ được mối liên hệ tường minh trực tiếp của cái phải tìm chỉ với cái đã biết nghĩa

là đã tìm được lời giải đáp cho bài toán đặt ra

− Hoạt động giải bài toán vật lý có hai phần việc cơ bản quan trọng là:

1 Việc xác lập các mối liên hệ cơ bản cụ thể dựa trên sự vận dụng kiến thức vật lý vào điều kiện cụ thể của bài toán đã cho

2 Sự tiếp tục luận giải, tính toán đi từ mối liên hệ đã xác lập được đến kết luận cuối cùng của việc giải đáp vấn đề được đặt ra trong bài toán đã cho

− Sự nắm vững lời giải của một bài toán vật lý phải thể hiện ở khả năng trả lời được câu hỏi: việc giải bài toán này cần xác lập được mối liên hệ nào? Sự xác lập các mối liên hệ cơ bản này dựa trên sự vận dụng kiến thức vật lý nào? Vào điều kiện cụ thể nào của bài toán?

− Đối với bài tập định tính, ta không phải tính toán phức tạp nhưng vẫn cần có suy luận logic từng bước để đi đến kết luận cuối cùng

1.1.3.2 Phương pháp giải bài tập vật lý

Xét về tính chất của các thao tác tư duy khi giải các bài tập vật lý người ta thường sử dụng hai phương pháp sau:

− Phương pháp phân tích: theo phương pháp này điểm xuất

phát là các đại lượng cần tìm Người giải phải tìm xem đại lượng chưa biết này có liên quan gì tới đại lượng vật lý khác, và khi biết được sự liên hệ này thì biểu diễn nó thành những công thức tương ứng, cứ làm như thế cho đến khi nào biểu diễn được hoàn toàn đại lượng cần tìm bằng những đại lượng đã biết thì bài toán đã được giải xong Như vậy, phương pháp này thực chất là đi phân tích một bài toán phức tạp thành những bài tập đơn giản hơn, rồi dựa vào những quy tắc tìm lời giải mà lần lượt giải các bài tập đơn giản này, từ đó đi đến lời giải cho bài toán phức tạp trên

− Phương pháp tổng hợp: theo phương pháp này suy luận

không bắt đầu từ đại lượng cần tìm mà bắt đầu từ đại lượng đã biết, nó nêu trong đề bài Dùng công thức liên hệ các đại lượng này với các đại lượng chưa biết, ta đi dần đến công thức cuối cùng, trong đó chỉ có một đại lượng chưa biết là đại lượng cần tìm

Nhìn chung, việc giải bài tập vật lý phải dùng chung hai phương pháp phân tích và tổng hợp Phép giải bắt đầu bằng phân tích các điều kiện của bài toán để hiểu đề bài và phải có sự tổng hợp kèm theo ngay

để kiểm tra lại mức độ đúng đắn của các sự phân tích ấy Muốn lập kế hoạch giải phải đi sâu phân tích nội dung vật lý của bài tập, tổng hợp những dữ kiện đã cho với những quy luật vật lý đã biết ta mới xây dựng được lời giải và kết quả cuối cùng

1.1.3.3 Các bước chung của giải bài toán vật lý

Từ phân tích về thực chất hoạt động giải bài toán, ta có đưa ra một cách khái quát các bước chung của tiến trình giải của một bài toán vật lý và các hoạt động chính trong các bước, đó là:

Bước 1:

− Tìm hiểu đề bài

− Đọc ghi ngắn gọn các dữ liệu xuất phát và cái phải tìm

− Mô tả lại tình huống đã nêu trong đề bài, vẽ hình minh họa

− Nếu đề bài yêu cầu thì phải dùng đồ thị hoặc làm thí nghiệm

để thu được các dữ liệu cần thiết

Bước 2: Xác lập những mối liên hệ cơ bản của các dữ liệu xuất phát và cái phải tìm

− Đối chiếu với các dữ liệu xuất phát và cái phải tìm, xem xét bản chất vật lý của những tình huống đã cho để nghĩ đến các kiến thức, các định luật, các công thức có liên quan

− Xác lập các mối liên hệ cơ bản, cụ thể của các dữ liệu xuất phát và của cái phải tìm

− Tìm kiếm lựa chọn các mối liên hệ tối thiểu cần thiết sao cho thấy được mối liên hệ của cái phải tìm với các dữ liệu xuất phát, từ đó

− Kiểm tra xem có tính toán đúng chưa

− Kiểm tra xem thứ nguyên có phù hợp không

− Giải bài toán theo cách khác xem có cùng kết quả không Tuy nhiên, trong nhiều bài tập không nhất thiết phải tách bạch một cách cứng nhắc giữa bước 2 và bước 3 Tùy từng bài toán mà ta

có thể kết hợp hai bước đó thành một trong tiến trình luận giải

1.1.3.4 Lựa chọn bài tập vật lý

Vấn đề lựa chọn bài tập vật lý góp phần không nhỏ vào việc nâng cao chất lượng học tập môn vật lý của người học và việc lựa chọn bài tập phải thỏa mãn các yêu cầu sau:

− Các bài tập phải đi từ dễ đến khó, đơn giản đến phức tạp, giúp người học nắm được phương pháp giải các bài tập điển hình

− Hệ thống bài tập cần bao gồm nhiều thể loại bài tập

− Lựa chọn các bài tập nhằm kích thích hứng thú học tập và phát triển tư duy của người học

− Các bài tập nhằm cũng cố, bổ sung và hoàn thiện tri thức cụ thể đã học, cung cấp cho người học những hiểu biết về thực tế, kỹ thuật có liên quan với kiến thức lý thuyết

− Lựa chọn các bài tập điển hình nhằm hướng dẫn cho người học vận dụng kiến thức đã học để giải những loại bài tập cơ bản, hình thành phương pháp chung để giải các loại bài tập đó

− Lựa chọn các bài tập sao cho có thể kiểm tra được mức độ nắm vững tri thức của người học

1.2 Bài toán biên

Xét phương trình vi phân tuyến tính có dạng

)()()

(

)()

()

1 1

dx

dy x a dx

y d x a dx

y d x a

y

n n

n

=+

+++

yc = C1y1(x)+C2y2(x)+…+Cnyn(x), (1.2.2)

trong đó: C1, C2, …, Cn là các hằng số tùy ý

Tiếp theo tìm bất cứ nghiệm riêng yp nào của phương trình vi phân không thuần nhất L(y) = F(x) Để giải phương trình này, ta thường dùng phương pháp hệ số bất định hoặc phương pháp biến thiên hằng số để tìm nghiệm riêng Khi đó nghiệm tổng quát của phương trình (1.2.1) sẽ là y = yc + yp

Trong các bài toán ứng dụng , nghiệm phương trình vi phân (1.2.1) đòi hỏi phải thỏa mãn các điều kiện bổ sung nào đó Số điều kiện này trong hầu hết các ứng dụng bằng cấp cao nhất của phương trình.Ví dụ, đối với phương trình vi phân cấp 2:

0)()

()

y d

x

a , a≤ x≤b (1.2.3)

bị lệ thuộc bởi điều kiện bổ sung tại x = a có dạng :

,)(,)

(a =α y/ a =β

y với α, β là các hằng số

Phương trình vi phân với điều kiện bổ sung được xem như là một bài toán cho trước giá trị ban đầu Bài toán giá trị ban đầu thường có nghiệm duy nhất Khi phương trình vi phân (1.2.3) bị hạn chế bởi hai điểm khác nhau, tức là x = a và x = b phương trình có dạng :

=+

≠+

=+

0,

)()

(

0,

)()

(

2 22

2 21

/ 22 21

2 12

2 11

/ 12 11

c c a

y c a

y

c

c c a

y c a

Điều kịên bổ sung (1.2.4) được gọi là điều kiện biên

Phương trình vi phân (1.2.3) với điều kiện biên (1.2.4) được gọi là bài toán biên

Nghiệm của bài toán biên phụ thuộc vào điều kiện biên Bài toán biên không chỉ có một nghiệm mà nó có vô số nghiệm Điều kiện biên có dạng

++

=+

++

β

α

)()

()()

(

)()

()()

(

/ 24 23

/ 22 21

/ 14 13

/ 12 11

a y c a y c b y c b

y

c

b y c b y c a y c a

dx =

Tích phân (1.2.7) thu được

C x

P C

x P C x P

e

C x P

−

1

) 1 )

) (

)(

ln y

Vậy nghiệm tổng quát của (1.2.6) là ( )

1

x P

C C e

Dùng phương pháp biến thiên hằng số và giả thiết một nghiệm riêng có dạng

) (

)

( P x

p u x e

y = − , trong đó C1 ở trong nghiệm tổng quát đã được thay thế bằng hàm chưa

biết u(x), nghiệm giả định này có đạo hàm là

) )

( ( ( )))

p

e dx

du x p e

x u

dx

dy

0

) ( )

)()

( (ξ) ξ ξ.Nghiệm tổng quát của phương trình (1.2.5) có dạng

1 0

Vì các hàm ψ và ϕ ở trong biểu thức (1.3.1) nói chung là những hàm phức

cho nên các toán tử trong trường hợp tổng quát cũng là những toán tử phức

Trong số tất cả những toán tử có thể có ta chỉ xét một lớp quan trọng những

toán tử, đó là những toán tử tuyến tính

Tóan tử A^ được gọi là một toán tử tuyến tính nếu nó thoả mãn đòi hỏi sau:

Nghĩa là ϕ là tổ hợp tuyến tính của hai hàm ϕ ϕ1, 2 Hay nói cách khác, kết quả tác dụng của một toán tử tuyến tính lên một hàm là tổ hợp tuyến tính của hai hàm ψ ψ1, 2 thì bằng tổ hợp của những kết quả tác dụng của toán tử đó lên mỗi hàm riêng biệt

1.3.2 Hàm riêng, trị riêng và phương trình trị riêng của toán tử

Nói chung khi cho toán tử A tác dụng lên hàm Ψ( )x thì ta được hàm

số ϕ( )x ≠ψ( )x (Với ( )x là tập hợp biến số nào đó) Nhưng cũng có trường hợp ta lại được chính hàm số đó nhân thêm với một hằng số Tức là:

^

Au n( )x = A n u n( )x

Trong đó u( )x là hàm riêng ứng với trị riêng A (n = 1;2;3;4;5…) Số n

trị riêng có thể có là hữu hạn hay vô hạn; có thể là gián đoạn hay liên tục

Để tìm hàm riêng và trị riêng của một toán tử, ta phải giải phương trình trị riêng của một toán tử đó

n n

n = ∂

∂

.ln

lnu n =iA n x= C

⇒

x iA

n Ce n

⇒ Với C là hằng số được xác định từ điều kiện chuẩn hóa

Vì hàm số tuần hoàn trong khoảng ( )0,L nên ta có u n( )0 =u n( )L Tức là:

1

=

⇒

=Ce iA n L e iA n L C

π

n L A

L A

n

n

2

1cos

Ta thấy An có giá trị gián đoạn theo số nguyên n

Còn hàm riêng tương ứng với An là: ( ) x

L ìn

n x Ce u

thực sự là nghiệm phương trình Kỹ thuật này sẽ được minh họa trong các ví dụ sau

Ví dụ 1: Cho U⊂ ℜn là một miền bị chặn với biên trơn Ta xét bài toán giá trị biên-ban đầu đối với phương trình truyền nhiệt

trong U× ,( )0 ∞ trên ∂U×[0,∞) (1.4.1 )

x = x x ∈ U và biến t∈ [ ] 0 , T tách ra với nhau

Bây giờ ta đi tìm v và w Để làm điều đó ta tính

0 0

t

(1.4.3 )

Với mọi x U∈ và t >0 sao cho v( )t , w( )t ≠ 0 Chú ý rằng vế trái của (1.4.3 chỉ )

phụ thuộc vào t và vế phải chỉ phụ thuộc vào x Điều này chỉ xảy ra khi chúng là hằng

số, tức là:

( )

( ) w ( ) ( ) x

x w t

Ta giải các phương trình này để tìm các hàm chưa biết w,v và hằng số µ

Trước hết, để ý rằng, nếu µ đã biết, nghiệm của (1.4.4) là v=deµt với d là hằng

số tùy ý Vì thế, ta chỉ cần nghiên cứu phương trình (1.4.5)

Ta nói rằng λ là một giá trị riêng của toán tử -∆ trong U ( với điều kiện biên bằng 0) nếu tồn tại một hàm w≠0 thõa mãn

0

µ

µ

µt

k t m

k t k

Sẽ là nghiệm của bài toán (1.4.1)

Đây là một công thức biểu diễn nghiệm rất đẹp, nhưng nó dựa vào:

Khả năng tìm các giá trị riêng, các hàm riêng và các hằng số thỏa mãn (1.4.9) Chuỗi (1.4.10) hội tụ theo một nghĩa thích hợp nào đó

Ví dụ 2: Tiếp theo ta sử dụng kỹ thuật tách biến để tìm nghiệm của phương trình môi trường tổ ong

Như ở ví dụ trứơc, ta tìm một nghiệm dạng

( ) ( ) ( ) x , t = v t w x

µ ( x ∈ ℜn; t ≥ 0 )

(1.4.12) Thế vào (1.6.11 , ta được )

( )

x w t

Với hằng số µ nào đó và với ∀ x ∈ ℜn, t≥ 0 sao cho w ( ) ( ) x , v t ≠ 0

Ta giải phương trình vi phân thường đối với v và tìm được

Ví dụ 3: Xét phương trình Hamilton- Jacobi

1.5 Phương pháp biến thiên tham số:

Để trực tiếp thu được nghiệm của phương trình

Xác định trong khoảng a <x < b, phụ thuộc vào hai điều kiện thuần nhất, trong đó

L là toán tử Sturm –Liouville có dạng

Nghiệm tổng quát L u ( ) = f x ( ) được cho bởi u u v u v = 1 1+ 2 2, ở đây v1và v2

được xác định bởi tích phân của dv dx1/ và dv dx2/ ở trên

Ta định nghĩa Wronskian w là đại lượng

CHƯƠNG II: XÂY DỰNG PHƯƠNG PHÁP

HÀM GREEN 2.1 Khái niệm hàm Green Tính đối xứng của hàm Green

Để đưa vào khái niệm hàm Green, chúng ta bắt đầu với các toán tử vi phân cấp 2

có dạng đồng nhất với hàm Green Mọi toán tử vi phân cấp 2 có dạng

)1.1.2()()

()

()

( 0 22 1 a2 x y

dx

dy x a dx

y d x a y

Các toán tử liên hợp đồng dạng với nó là:

)2.1.2()()]

()(2[)

()(

)(])([])([)

(

2 1

/ 0 2

2 0

*

2 1

0 2

2

*

y x a dx

dy x a x a dx

y d x a y

L

y x a y x a dx

d y x a dx

d y

Cho hai hàm u(x) và v(x) là hai hàm liên tục tùy ý cùng với đạo hàm cấp 1 và cấp

2 của nó Dùng hai toán tử (2.1.1) và (2.1.2) để xác định đồng nhất thức Lagrange của hai hàm u(x) và v(x) như sau

) 3 1 2 ( )]

, ( [ ) ( )

dx

d v uL u

vLx − x =

Trong đó ( , ) ( ) ( ) ( ) /( ) (2.1.4)

0 0

1

dx

dv x a v u x a dx

du x a v u P

Được gọi là hàm song tuyến

Đồng nhất thức Lagrange của 2 hàm khả vi u(x) và v(x) được xác định trên miền

{ x a x b }

I = / ≤ ≤ Tích phân đồng nhất thức (2.1.3) ta có đồng nhất thức các hàm Green

[ ( , ) ] , ( 2 1 5 ) )]

( )

( ) ( ) ( ) ( [ ) ( )]

( ) ( ) ( )

(

[

) ( )]

( ) ( ) ( ) ( [ ) ( )]

( ) ( ) ( ) ( [ )

,

(

/ 0

/ 0 1

/ 0

/ 0

/ 0 1

/ 0

a u a v a a a v a a a v a u a a a u a

a

b u b v b a b v b a b v b u b a b u b a v

u

P b a

+

− +

−

− +

− +

(

,

( v Lx u u L*x v các hạng thức trên biên

Sử dụng đồng nhất thức Green cho thích hợp để tìm nghiệm phương trình với biên

ở hai điểm như sau

) 7 1 2 ( )

(

; )

(

);

( )

(

2 2

1 1

B

x f y

Lx

Trong đó: Lx là toán tử tuyến tính cho bởi (2.1.1); g1; g2 là các hằng số và B1, B2

là các toán tử biên tuyến tính dạng Robin:

b x

a x

x y dx

x dy y

B

x y dx

x dy y

( )

(

; ) ( )

( )

(

2 2

2

1 1

1

β α

β α

( ) ( ) ( )

v u P d v L u u L

L là toán tử đạo hàm đối với ξ

Để giải phương trình (2.1.9) với điều kiện (2.1.10), đặt u(ξ) = y(ξ) là nghiệm của phương trình (2.1.9) với x thay bằng ξ và u thay bằng y trong đồng nhất thức Green Như vậy, đồng nhất thức Green (2.1.11) thay Lξ( y ) = f ( ξ ) ta có

∫

a

b a b

a

v y P d v L y d f

v ( ξ ) ( ξ ) ξ ( ξ ) *ξ( ) ξ [ ( ( ξ ), ( ξ ))] ( 2 1 12 )

Trong đó:

)()]

()()()([)()]

()()()

(

[

)()]

()()()([)()]

()()()([)

,

(

/ 0

/ 0 1

/ 0

/ 0

/ 0 1

/ 0

a y a v a a a v a a a v a y a a a y a

a

b y b v b a b v b a b v b y b a b y b a v

−

−+

−+

=

(2.1.13)

Chọn v(ξ) = G*(ξ;x) là hàm Green thỏa mãn điều kiện

)14.1.2(,),

())

;(

y( ) ( ) ( ) (2.1.15)

Thay v(ξ)=G*(ξ;x) vào đồng nhất thức Green (2.1.12),rút gọn thành

)16.1.2())

;(),(())

;(),(()()()()

)

())

)

;(][,0)

;(]

[

)

;()]

()([)

;()]

()([)]

;(),

* / 0

x G G B x

G G

B

x a G a y a a x b G b y b a x

Điều kiện biên này được gọi là điều kiện biên liên hợp Từ đó ta có nghiệm của (2.1.16) là

) 18 1 2 ( ) ( )

; ( )

Trong đó G*(ξ;x) là hàm Green thỏa mãn phương trình

b a

x x

G

L*ξ( *( ξ ; )) = δ ( − ξ ), ≤ ξ ≤

, với các điều kiện biên

) 19 1 2 ( 0 )

; ( ] [ , 0 )

; ( ]

Như vậy, để tìm nghiệm của phương trình (2.1.9), ta đi tìm hàm Green G*(ξ;x)

Đó chính là phương pháp tìm nghiệm mới, được gọi là phương pháp hàm Green

Nhằm mục đích xây dựng phương pháp hàm Green ta đưa ra 2 hàm Green G và hàm Green liên kết G* thỏa mãn các toán tử Lx và L*x cho bởi phươnng trình (2.1.20) và (2.1.21) sau:

) 21 1 2 ( ),

( )

;

(

) 20 1 2 ( ),

( )

ξ δ

ξ

Với các điều kiện biên

0)

;(][,0)

;(]

P Hàm Green cho bởi phương trình (2.1.20) và (2.1.21) thỏa mãn quan hệ đối xứng

G*(x;ξ) = G(ξ;x) (2.1.22)

Để chứng minh tính đối xứng trên, nhân phương trình (2.1.20) với G*(x;t) và sau

đó thay biến ξ trong phương trình (2.1.21) bằng biến t, rồi nhân phương trình (2.1.21) với G*(x;ξ) ta thu được:

) ( )

; ( )

;

(

) ( )

; ( )

;

(

*

*G x t x t L

x

G

x x

G L

ξ δ

;()

;()]

()

;()()

;(

[

))]

;(()

;())

;(()

;([]

δ ξ ξ

δ

ξ ξ

t G t G t x x G x

t x G

t x G L x G x

G L t x G G

b a

Từ đó suy ra (2.1.22), và gọi là tính chất đối xứng của hàm Green Như vậy nghiệm của bài toán Dirichlet (2.1.18) có dạng

)25.1.2()()

;()

()

;()

2.2 Xây dựng phương pháp hàm Green

Xét phương trình truyền nhiệt tổng quát có nguồn nhiệt, điều kiện biên thuần nhất:

2 2

,

(

sin ) ( )

,

(

n n

n n

l

x n t q t

x

Q

l

x n t u t

x

u

π π

Bước 2: Thay vào phương trình (2.1), tìm nghiệm:

Phương trình truyền nhiệt ( ) ( ) ( )

2

t q t u l

a n dt

t du

n n

Nghiệm có dạng :

t l

a n t

l

a n n

u

0

2 2

2

) ( )

0 ( )

ππ

n

n

d l

n l

u

l

x n u

x

0

1

sin ) (

2 ) 0

(

sin ) 0 ( )

(

ξ

πξ ξ

ϕ

π ϕ

,

(

n n

l

x n t q t

l t

q

0

sin ) , (

2 )

) (

2 )

,

a n l

e d l

n l

ϕ

l

x n d

e d l

n Q

l e

l a n l

t l a n

π τ

ξ

πξ τ

ξ

π π

sin sin

) , ( 2

0 0

2 2

πξ τ

ξ

ξ π

πξ ξ

ϕ

τ π π

d d e

l

x n l

n l Q

d e

l

x n l

n l t

x

u

n

t l a n

l t

n

t l a n l

0 0

1 0

2 2

sin sin

2 ) , (

sin sin

2 ) ( )

2sin

sin

2 )

x n l

n l t

x

π

π πξ

τ ξ

Như vậy để tìm nghiệm của phương trình ta đi tìm hàm Green Đó chính là phương pháp tìm nghiệm mới được gọi là phương pháp hàm Green

2.3 Hàm riêng, trị riêng cho hàm Green

Xét phương trình vi phân không thuần nhất Sturm_ Liouville tổng quát:

Giả thiết hai điều kiện biên là thuần nhất, ta đưa vào một bài toán trị riêng tương ứng:

)()

(

n

n n n n

( )

n n dx

n n dx

x

a n

a n

b a n n

σφ

φλ

Nghiệm của bài toán giá trị biên cho phương trình vi phân không thuần nhất là

ξξξξ

σφλ

ξφφ

ξ d f G x d

dx

x f

x

u

b a

a n n

n n b

a

),()(()

()()(()

, (

a n n

n n

dx

x x

G

σφ λ

ξ φ φ ξ

Áp dụng kết quả trên để giải bài toán:

0 ) ( , 0 ) 0 ( );

sin sin2

2.4 Hàm điều hoà Biễu diễn Green

Giả sử Ω là một miền trong R còn u là hàm thuộc lớp n C Ω Hàm 2( ) u x thỏa ( )

mãn phương trình Laplaxơ

0

u

∆ = (2.4.1) với mọi x thuộc Ω được gọi là hàm điều hoà trong Ω Dạng không thuần nhất

của phương trình Laplaxơ được gọi là phương trình Poisson Nghiệm của phương trình

Poisson trong miền Ω là hàm u x( ) thuộc lớp C Ω sao cho 2( )

( )

u f x

với bất kỳ x thuộc Ω.Nghiệm như thế còn được gọi là nghiệm cổ điển của

phương trình Poisson trong miền Ω

Giả sử c là một miền bị chặn trong R với biên ∂Ω thuộc lớp n B và giả sử 1

Trong đó υ là pháp vectơ đơn vị ngoài tới ∂Ω , ds là phần tử diện tích ∂Ω

Từ công thức này ta nhận được công thức tính tích phân từng phần:

Đổi vai trò u và v trong công thức (2.4.4), sau đó lấy (2.4.4) trừ đi công thức vừa

nhận được, ta có công thức Green thứ hai