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 "

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.

- 1-

• 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.

- 2-

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ế.

- 3-

+ 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?

- 4-

− Đố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.

- 5-

− 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ừ đó có thể rút ra được cái cần tìm.

Bước 3: Rút ra kết quả cần tìm.

Từ các mối liên hệ cần thiết đã xác lập, tiếp tục luận giải, tính toán để rút ra kết quả cần tìm.

Bước 4: Kiểm tra xác nhận kết quả, để có thể xác nhận kết quả cần tìm cần kiểm tra lại việc giải theo một hoặc một số cách sau:

− 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

1 −

(1.2.1)

yL )(

x )(

xFyxa )(

)(

... ++

xa )( 0

xa )( 1

1 −

y 1 −

dy dx

n yd n dx

d dx

x ≤≤

Trong đó: a0(x), a1(x),…,an(x) là các hàm liên tục trong khoảng

trong khoảng

x ≤≤

và ≠xa .Cách chung để giải phương trình (1.2.1) là: trước hết )(0 giải phương trình thuần nhất cấp n là L(y)=0, thu được một tập nghiệm cơ bản {y1(x), y2(x),…, yn(x)}, nghiệm tổng quát yc của phương trình thuần nhất là một tổ hợp tuyến tính của tập nghiệm cơ bản:

- 6-

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

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:

(1.2.3)

x ≤≤

xa )( 1

xa )( 2

)( xa 0

dy dx

2 yd 2 dx

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

ay )(

/ ay )(

, α

, β

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 :

≠

, α

(1.2.4)

2 c 11 c

2 c 12 c

≠

, β

/ )( ayc 12 / )( ayc 22

2 21

2 22

⎧ )( ayc ⎪ 11 ⎨ ⎪⎩ )( ayc 21

Trong đó:

α,

và β là các hằng số

c 11

c 12

Đ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

α + + + =

/ ayc )( 12 / byc )( 22

/ byc )( 14 / ayc )( 24

β + + + = byc )( 13 ayc )( 23 ⎧ ayc )( ⎪ 11 ⎨ ⎪⎩ byc )( 21

Trong đó: cij , i= 1,2, j=1,2,3,4 và α,β là các hằng số; được gọi là điều kiện biên hỗn hợp.

L y ( )

p x y ( )

q x ( ).

(1.2.5)

dy dx

Bài toán biên hỗn hợp thường khó giải. Xét phương trình tuyến tính cấp 1:

L y ( )

p x y ( )

(1.2.6)

dy dx

Để giải phương trình (1.2.5), trước hết giải phương trình thuần nhất:

- 7-

Để thu được nghiệm tổng quát yc. Ta có thể tách biến phương trình (1.2.6) có dạng:

= −

(1.2.7)

p x dx ( )

dy y

Đặt:

p x ( )

(1.2.8)

∫=

dP dx

)( xP p d ξξ với )( =

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

)( CxP

)( xP

−

ln y )( CxP −= +

)( xP

−

e e e e , e ⇒ = = =⇒ = eCy 1 C 1

xP )(

−

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

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 , trong đó C1 ở trong nghiệm tổng quát đã được thay thế bằng hàm chưa exu )(

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

xP )(

−

p dx

dy ( )) e )( exu ( xp = − + du dx

)( xP

−

yxp )(

xq )(

)( xq

)( xu

⇒

u =⇒

( ) P ξ ξ e )(

d . ξ

∫=

p dx

du dx

Thay yp và dy p vào phương trình không thuần nhất (1.2.5) ta có dx

xP )(

−= e

P ( ) ξ ξ e )(

d . ξ

∫

Suy ra nghiệm riêng

−

P x ( )

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

P dξξ ( ) e q ( )

(1.2.9)

y C

∫

= + = + ⎡ C ⎢ 1 ⎣ ⎤ ξ ⎥ ⎦

1.3 Khái niệm toán tử, hàm riêng, trị riêng:

1.3.1.Toán tử:

^ ψ ϕ=A

Toán tử là một qui tắc toán học dùng để biến đổi một hàm này sang hàm khác có cùng bản chất:

(1.3.1)

^ A được kí hiệu bằng chữ có mũ, còn qui tắc tác

^ A và được gọi là

Trong đó toán tử

dụng của nó được viết dưới dạng một phép nhân ψ với

^ A tác dụng lên hàm ψ cho hàm ϕ ”

- 8-

“toán tử

Để làm sáng tỏ ta xét các ví dụ sau:

^ x = x

^ A =

Ví dụ 1: Toán tử toạ độ hay toán tử nhân

^ x ψ

x ψ

( ) x

(1.3.2)

^ A

^ d dx

d dx

Ví dụ 2: Toán tử vi phân

( ) x

^ d dx

d dx

(1.3.3)

2 ∇ =

∂ x ∂

∂ y ∂

∂ z ∂

Ví dụ 3: Toán tử Laplace:

= ∇

2 ψ

( ) x

2 2 ∂ ∂ ψ ψ ψ 2 2 2 y x ∂ ∂

∂ z ∂

(1.3.4)

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.

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

Tóan tử sau:

)

^ ( A C

C ψ ψ + 1 1 2

^ C A 1

ψ 1

^ C A ψ 2

,C C là những hằng số phức tuỳ

,ψ ψ là hai hằng số tuỳ ý, còn

(1.3.5)

trong đó ý.

^ A biểu thị nguyên lí chồng chất trạng thái. Thật vậy theo định nghĩa về toán tử (1.3.1) ta có:

Từ định nghĩa (1.3.5) ta thấy tính chất tuyến tính của toán tử

^ ^ A ,ψ ϕ ψ ϕ A 2

^ ( A C

)

^ C A 1

ψ 1

^ C A 2

ψ ϕ ϕ 1 1 2 2

ψ ψ ϕ 1 1 2

- 9-

, nên (3.5) trở thành:

,ϕ ϕ . 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 ,ψ ψ thì bằng tổ hợp của những kết quả tác dụng tuyến tính của hai hàm của toán tử đó lên mỗi hàm riêng biệt.

Nghĩa là ϕ là tổ hợp tuyến tính của hai hàm

^ x ;

2∇ trong các thí dụ (1), (2),

^ d dx

và Rõ ràng rằng các toán tử

(3) là những toán tử tuyến tính. Còn toán tử chứa tác dụng lấy căn số (hay

toán tử ) không phải là toán tử tuyến tính. Tất cả các toán tử dùng trong phương trình vật lý toán là những toán tử tuyến tính ; do đó về sau này khi nói đến toán tử là ta ngụ ý nói đến các toán tử tuyến tính.

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

( )xΨ thì ta được hàm . (Với ( )x là tập hợp biến số nào đó). Nhưng cũng có trường

( )x

^ A ( )xΨ = A ( )xΨ

Nói chung khi cho toán tử A tác dụng lên hàm ( ) x ψϕ ≠ số 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à:

^ A và phương trình trên

( )xΨ là hàm riêng của toán tử

Khi đó ta nói

^ A . còn A được gọi là trị riêng ứng

gọi là phương trình trị riêng của toán tử

^ A .

( )xΨ của toán tử

với hàm riêng

( )xuA

^ A ( ) xu n

Một toán tử có thể có nhiều hàm riêng và mỗi hàm riêng thì ứng với một trị riêng ( cũng có thể có trường hợp một trị riêng ứng với nhiều hàm riêng, trường hợp này gọi là trị riêng có suy biến), nên ta đánh chỉ số để phân biệt các phương trình trị riêng và được viết như sau:

( )xu

nA (n = 1;2;3;4;5…). Số

Trong đó

là hàm riêng ứng với trị riêng 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ử đó.

−

^ A =

∂

∂ ( )x

Thí dụ: Cho toán tử . Hãy tìm hàm riêng và trị riêng của

^ A , biết rằng hàm riêng tuần hoàn trong khoảng (

)L,0

là trị riêng và hàm riêng tương ứng của toán tử

Gọi

^ A

nA và

( )xun

thì phương trình trị riêng của toán tử

^ A là:

- 10-

toán tử

( ta không viết đối số tọa độ để khỏi rờm rà)

−

uA n

u ∂ n ( ) x ∂

x ∂

iA n

u ∂ u

⇒

iA n

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

iA nCe

u =⇒ n

nên ta có

. Tức

Vì hàm số tuần hoàn trong khoảng (

)L,0

( ) =0

( )Lu

là:

LiA n

cos

⇒

e ⇒ ) πn

LA n =

⇒

( LA n

Từ đó ta được

. ( n= 0;2;3;…). = An π2 n L

ìn π L

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

( ) xu n

. Còn hàm riêng tương ứng với A n là:

1.4. Phương pháp tách biến

Phương pháp tách biến nhằm xây dựng một nghiệm µ của phương trình đạo hàm riêng cho trước thông qua các hàm có biến số ít hơn. Nói cách khác, ta phỏng đoán rằng µ có thể được viết dưới dạng tổng hoặc tích của các hàm có biến số ít hơn và tách nhau, thay nó vào phương trình đạo hàm riêng để chọn các hàm đó phải đảm bảo µ 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

)∞× ,0 trong U ( [ )∞×∂ ,0U

) 1.4.1

trên (

}0=× t

0 g

µ µ ⎧ t ⎪ µ =⎨ ⎪ =⎩ µ

trênU {

∈ tUx

; ≥

=,µ tx

Ở đây g: U ℜ→ là hàm cho trước. Ta giả định tồn tại một nghiệm có dạng

)

(

( ) ( )xwtv

) 1.4.2

1.4.1 như là tích của hai hàm số với các biến

)

( ; (

)

tách ra với nhau.

- 11-

Có nghĩa là, ta xem nghiệm của ( ( ]T,0∈ [ x x ∈ và biến t 1,......., n Bây giờ ta đi tìm v và w. Để làm điều đó ta tính

∆

) , =

( ) ( )xwtv

( ,µ tx

)

( )xwtv ( ) ∆

( txt

−

− ∆

x t ,

= ( ) ( )xwtv,

( ( ) v t w x ∆

)

( µ

)

( t x t ,

Từ đó

) 1.4.3

( ) , tv ( ) tv

( Nếu và chỉ nếu ( ) xw ∆ ( )xw

= µ

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

Ux ∈

≥ ,0

) 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à: ( ) , tv ( ) tv

( ) xw ∆ ( )xw

( t ).

/v µν=

Khi đó

w µ=∆ w

(1.4.4)

(1.4.5)

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).

0≠w

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 thõa mãn bằng 0) nếu tồn tại một hàm

. trong U w λ w =∆− ⎧ ⎨ trên U∂ 0w = ⎩

λµ −=

de tλ −−=

Ta gọi hàm w là hàm riêng tương ứng, ta đặt để tìm

(1.4.6)

Thỏa mãn

(1.4.7)

µµt ⎧ ⎨ 0 = µ ⎩ trên 0 =∆− trong U ( )∞× ,0 [ )∞×∂ ,0U

) dw=0.,µ

(

g =

với điều kiện ban đầu . Do đó hàm µ được xác định bởi (1.4.6) thỏa

nλλ,....., 1

d ,....,

là các giá trị riêng ,

- 12-

là các hàm riêng tương ứng và là các hằng số, thì mãn (1.4.1), với điều kiện w ,......, nw . Tổng quát hơn, nếu md

−

kλ

ed k

w k

∑=

1 =

∞

(1.4.8)

) =0.,µ (

k wd∑

1 =

∞

wd k k

,.....

. Nếu ta có thể tìm Thỏa mãn các điều kiện ban đầu

1wm ,

∑

,.......

v.v. sao cho trong U (1.4.9) được

1 dd ,

∞

−

kλ

khi đó Với các hằng số thích hợp

ed k

w k

∑=

1 =

(1.4.10)

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.

( )∞×ℜ ,0

) 0=

µt

( γµ Trong đó nghiệm

0≥µ và

1>γ là hằng số. Đây là một phương trình khuếch tán

phi tuyến, với tốc độ khuếch tán của mật độ µ phụ thuộc vào chính µ.

n ; ≥ℜ∈

=,µ tx

(1.4.12)

( ) ( )xwtv

(

)

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

1.6.11 , ta được

Thế vào (

)

(1.4.13)

( ) γ xw ∆ ( )xw

( ) , tv γ µ = ( ) tv

x ℜ∈∀

trong (1.4.11)

( ) ≠tvxw ,

( )

, t 0≥ sao cho Với hằng số µ nào đó và với

1 −

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

0>λ nào đó.Để tìm w ta xét phương trình

) λλµγ +

)

( ( 1

, với hằng số

v − đạo hàm riêng ( ∆ γ w =

) uw

αxw =

(1.4.14)

Ta dự đoán rằng nghiệm w có dạng với hằng số α sẽ được xác định

- 13-

sau. Khi đó

αλ −

∆−

−

n −=

( αγαγ

)

( γ w

)

(1.4.15)

nℜ , trước hết ta đòi hỏi rằng

= αγα

2−

Vì vậy, để (1.4.14) thỏa mãn trong và từ

đó

(1.4.16) 1 2 −

Tiếp theo, từ (1.4.15) dễ thấy rằng cần đặt

−+

( nαγαγµ

) 2 >

(1.4.17)

0>λ hàm

−

( ( 1

) ut γ

1 ) γλ 1 −

Tóm lại, với mỗi

1.4.11 , các tham số µα,

)

(1.4.17).

Trong các ví dụ trên, sự tách biến được thực hiện dựa vào tính thuần nhất phi tuyến tương thích với hàm µ có dạng tích (1.4.12). Ở trường hợp khác, ta sẽ tìm nghiệm, trong đó các biến được tách dưới dạng một tổng các hàm số.

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

trong

(1.4.18)

(

)

( )∞×ℜ ,0

+ DuHtµ

=,µ tx

(

).0

n , ≥ℜ∈

Và tìm một nghiệm µ có dạng ( )xwtv ( ) +

0 µ =

) )

(

)

)xDwHtv ( ( )

( )

) ( Khi đó ( ( , txDuHtxt

,−=

n , >ℜ∈

Nếu và chỉ nếu ( ) ) ( xDwH

( )tv

( x

),0

Với hằng số µ nào đó. Vì thế, nếu ) µ=DwH

(

,R∈µ

=,µ tx

−

Với (

)

thì ( ) xw

với hằng số b nào đó. Đặc biệt, nếu chọn

b (

= n

Sẽ thõa mãn .= xa

DuHttµ và đặt

với

, tìm được nghiệm

( ) xw

) 0= ( )aH=µ

−

= .µ

a ℜ∈ ( ) taHxa

Dựa vào tích phân đầy đủ và hàm bao tìm được nghiệm.

- 14-

được xác định bởi (1.4.16), Thỏa mãn phương trình (

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

f x ( )

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

, ta xét bài toán

2 d u 2 dx

không thuần nhất tổng quát:

L u ( )

f x ( )

Xác định trong khoảng a

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

q .

d dx

du dx

⎛ ⎜ ⎝

⎞ ⎟ ⎠

Khi p = 1, q = 0 ta được toán tử của phương trình truyền nhiệt trong trạng thái

dừng

2 d u 2 dx

Phương trình vi phân thường không thuần nhất luôn có thể giải bằng phương pháp u x và 1( ) biến thiên tham số, nếu biết hai nghiệm của phương trình không thuần nhất u x . Theo phương pháp biến thiên tham số, nghiệm riêng của phương trình 2 ( ) L u ( )

f x ( )

được tìm dưới dạng

u u v = 1 1

u v 2 2

Khi đó

1v và

2v là hàm phụ thuộc vào x chưa được xác định. Phương trình vi phân gốc có một hàm chưa biết, vì rằng có một bậc tự do thêm vào là du/dx. Nếu 1v và 2v là hằng số thì

v 1

v 2

du dx

du 1 dx

du 2 dx

Vì 1v và 2v không phải là hằng số nên

= 0

u 1

u 2

dv 1 dx

dv 2 dx

L u ( )

f x ( )

được thoã mãn nếu

Vi phân

f x ( ).

dv 1 dx

du 1 dx

dv 2 dx

du 2 dx

Phương pháp biến thiên tham số tạo ra hai phương trình vi phân cho các hàm

chưa biết

1 /dv dx và là:

- 15-

−

fu 2

;

dv 1 dx

fu 2 c

−

u 2

du 1 dx

⎛ p u ⎜ 1 ⎝

⎞ ⎟ ⎠

−

fu 1

dv 2 dx

fu 1 c

−

u 2

− du 2 dx − du 2 dx

du 1 dx

⎛ p u ⎜ 1 ⎝

⎞ ⎟ ⎠

−

u 2

Trong đó

, hằng số c tuỳ thuộc vào việc lựa chọn 1u

du 2 dx

du 1 dx

⎛ p u ⎜ 1 ⎝

⎞ ⎟ ⎠

và

2u .

L u ( )

được cho bởi

Nghiệm tổng quát

, ở đây 1v và

u u v = 1 1

u v 2 2

được xác định bởi tích phân của

f x = ( ) 1 /dv dx và

2 /dv dx ở trên.

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

−

W u = 1

u 2

du 2 dx

du 1 dx

Nó thoả mãn phương trình vi phân cơ bản

−

= −

−

= −

u 1

u 2

u 1

u 2

dW dx

2 d u 2 2 dx

2 d u 1 2 dx

dp dx / p

du 2 dx

du 1 dx

dp dx / p

⎛ ⎜ ⎝

⎞ ⎟ ⎠

L u = và ) 0

L u = được ) 0

Trong đó, các phương trình vi phân thuần nhất

/W c p =

dùng đến. Giải phương trình trên suy ra

1( hay là pW c= .

Tiểu kết :

Trong chương này chúng tôi đã nêu lên khái niệm về bài tập vật lý, tầm quan trọng của bài tập vật lý. Trình bày các cơ sở toán học cơ bản, cần thiết cho việc xây dựng phương pháp hàm Green.

- 16-

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

yxa )(

)1.1.2(

yLx )(

xa )( 0

xa )( 1

dy dx

2 yd 2 dx

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

[

yxa )(

])( yxa

−

* )( yL x

yxa ])( 1

d dx

)]

yxa )(

)2.1.2(

⇒

−

* yL )( x

/ xa )(2[ 0

xa ( 1

xa )( 0

dy dx

2 yd 2 dx

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

u )(

v )(

[

vuP ,(

)3.1.2()]

−

* uL x

d dx

vuP ),(

uvxa )(

)4.1.2(

Trong đó

xa )( 0

/ 0

du dx

dv dx

⎡ ⎢ ⎣

⎤ vuxa )( −⎥ ⎦

⎡ ⎢ ⎣

⎤ ⎥ ⎦

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

x ≤≤

Đồ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 { I / ax = . Tích phân đồng nhất thức (2.1.3) ta có đồng nhất thức các hàm Green.

[

u )(

(

)]

)5.1.2(,

−

[ vuP ),(

]

* uL x

b a

∫

[

buba )( )(

()(

)]

[

bvba )( )(

bubvba )(

()(

)]

−

bvbuba )( 1

/ 0

[

()(

)]

[

avaa )( )(

auavaa )(

()(

)]

Trong đó [ ] b vuP ),( a auaa )( )( −

−

avauaa )( 1

/ 0

Định nghĩa tích hàm của đồng nhất thức Green:

))

)( dxu .

)6.1.2(

,( uLv ( x

∫=

Tích phân từng phần tích hàm thu được

))

các hạng thức trên biên.

uLv ( ,( x

* vLu ,( ( x

- 17-

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

xf (

);

)7.1.2(

yLx )( yB )(

;

yBg )( 2

⎧ ⎨ 1⎩

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:

xy )(

;

yB )( 1

β 1

xdy )( dx

⎡ α ⎢ 1 ⎣

⎤ ⎥ ⎦

ax =

(2.1.8)

xy )(

yB )( 2

β 2

xdy )( dx

⎡ α ⎢ 2 ⎣

⎤ ⎥ ⎦

bx =

Sử dụng đồng nhất thức Green để giải bài toán biên Dirichlet

f x ( );

(2.1.9)

y a = ( ) 0,

y b ( ) 0 (2.1.10)

xL y ( ) B y ( ) 1

B y ( ) 2

Để giải phương trình này, đổi biến x trong phương trình (2.1.1) và (2.1.5) thành

biến mới ξ và viết đồng nhất thức Green theo biến mới

(

))

)11.1.2(,

)([ v ξ

−

)( ξ

), v ξξ

[ uP ((

]

uL )( ξ

* dvL )] ( ξ

b a

∫

Trong phương trình (2.1.11), biến ξ được dùng như một biến giả của phép lấy

tích phân và vì thế các toán tử

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

L ξ

* Lvà ξ

Để 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

)( ξ

ta có

Green. Như vậy, đồng nhất thức Green (2.1.11) thay

yL )( ξ

[

(( yP

(

))]

)12.1.2(

)( )( ξξξ

−

)( ξ

), v ξξ

b a

* )( dvL ξ

∫

[

byba )( )(

()(

)]

[

bvba )( )(

bybvba )(

()(

)]

−

bvbyba )( 1

/ 0

] b a ayaa )( )(

()(

)]

[

avaa )( )(

ayavaa )(

()(

)]

−

avayaa )( 1

/ 0

Trong đó: [ vyP ),( [ − (2.1.13)

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

))

(

)14.1.2(,

;( ξ

x ), ξδ −

ξ ≤≤

* GL ( ξ

Nó là phương trình liên hợp với đạo hàm trong

ξL (đạo hàm theo biến ξ),

là hàm Delta Dirac có tính chất

( ξδ −x )

- 18-

) dx

)( xy

)15.1.2(

()( ξδξ

−

∫

x += εξ y x −= εξ

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

xy )(

xbGbyP ),

((

;(

))

xaGayP ),

((

;(

))

)16.1.2(

);( fx ξ

)()( d −ξξ

−

∫

))

(

Theo điều kiện (2.1.10) ta có

). x ξδ −

;( ξ

Nghiệm y(x) trong bài toán (2.1.9) có thể thu được bằng kết quả của tích phân (2.1.16) . Chúng ta sẽ nghiên cứu kỹ hơn tích phân (2.1.16). Hàm f(ξ) đã cho từ phương trình (2.1.9), hàm G*(ξ;x) thu được từ việc giải phương trình (2.1.14) có dạng * ( GL ξ

yP ([

[

xbGbyba ( );(

)(

)]

[

xaGayaa ( );(

)(

)]

), ξ

;( ξ

−

]

)17.1.2(

b a );( x ξ

0 =

0 );( x ξ

* GB [ 1

* GB [ 2

ξ =

Đ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à

xy )(

)18.1.2(

);( fx ξ

)( d ξξ

∫=

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

))

(

;( ξ

), x ξδ −

ξ ≤≤

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

* ( GL ξ

]

xaG );(

]

xbG );(

)19.1.2(0

* GB [ 1

* GB [ 2

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à x cho bởi phươnng trình (2.1.20) và

hàm Green liên kết G* thỏa mãn các toán tử Lx và L* (2.1.21) sau:

xGL

);(

)20.1.2(

( ), ξδξ

−

x ≤≤

);(

)21.1.2(

( ), ξδξ

−

x ≤≤

* * xGL x

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

]

xaG );(

]

xbG );(

* GB [ 1

* GB [ 2

Trong các phương trình trên, các vi phân lấy theo biến x các toán tử

), điều kiện biên liên hợp được

( Lx, B1, B2) có dạng liên hợp của nó là (

* * BBLx 1

* 2

)

. Hàm Green cho bởi phương trình (2.1.20) và (2.1.21) thỏa

chọn là

=b GGP , ( mãn quan hệ đối xứng

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

- 19-

Để 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:

xGLtxG ;(

);(

)

( ) ξδξ

−

(2.1.23)

xG ;(

)

) ξ

( δ

−

x * * txGL );( x

Trừ hai phương trình trên và sau đó tích phân từ a đến b ta thu được đồng nhất

thức Green

GGP ,

[

]

xGLtxG [ ;(

);(

(

))

* xG ;(

(

))]

−

) ξ

* * txGL ;( x

∫

* txG [

();(

xG ;(

tGt ;(

)24.1.2(0

x ) ξδ −

−

() δξ

−

Gt )] =

);( ξ

−

) ξ

∫

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

xy )(

xG ;(

dxxf )()

)25.1.2(

)();( fx ξ

d ξξ

∫

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:

l a + x < < = Q x t ( , ) 0

= u l t ( , ) 0 (2.1)

= = xϕ ( )

⎧∂ u u ∂ ⎪ 2 x t ∂ ∂ ⎪⎪ t u (0, ) 0, ⎨ ⎪ u x ( , 0) ⎪ ⎪⎩

Bước 1: Áp dụng phương pháp tách biến Fourier và phương pháp mở rộng hàm

riêng ta chọn nghiệm có dạng

∞

sin)(

),( txu

tu n

∑

xn π l

1 =

∞

sin)(

),( txQ

tq n

∑

xn π l

1 =

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

t )(

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

tu )( n

tq )( n

n dt

an π l

⎛ ⎜ ⎝

⎞ ⎟ ⎠

Nghiệm có dạng :

- 20-

an π l

⎛ − ⎜ ⎝

⎞ ⎟ ⎠

⎛ − ⎜ ⎝

⎞ ⎟ ⎠

⎛ ⎜ ⎝

⎞ ⎟ ⎠

)0(

)( τ

d τ

)( tu n

∫

Dựa vào điều kiện ban đầu tìm hàm un(0)

∞

)( x

sin)0(

∑

xn π l

1 =

)0(

sin)(

⇒

ξϕ

d ξ

∫

n πξ l

2 l

∞

sin)(

),( txQ

tq n

∑

xn π l

1 =

⇒

sin),( τξ

d ξ

)( tq n

n πξ l

2 ∫= l

Cuối cùng ta thu được

∞

−

an π l

⎛ ⎜ ⎝

⎞ ⎟ ⎠

txu ),(

sin)(

ξϕ

n πξ l

⎞ ⎟ d ξ ⎟ ⎠

⎡ ⎛ 2 ∑ ∫ ⎢ ⎜ ⎜ l ⎢ ⎝ ⎣

−

an π l

⎛ ⎜ ⎝

⎞ ⎟ ⎠

⎛ ⎜ ⎝

⎞ ⎟ ⎠

sin

sin),( τξ

∫

2 l

n πξ l

xn π l

⎛ ⎜ ⎜ ⎝

⎞ ⎟ d ξ ⎟ ⎠

⎤ ⎥ d τ ⎥ ⎦

Đổi thứ tự giữa tổng và tích phân ta thu được

∞

an π l

⎞ ⎟ ⎠

⎛ − ⎜ ⎝

sin

( ), txu

( ) ξϕ

∑

∫

2 l

πξ l

xn π l

⎛ ⎜ ⎜ ⎝

⎞ ⎟ ⎟ ⎠

∞

)

( t

−

an π l

⎞ ⎟ ⎠

⎛ ⎜ ⎝

sin

( , ) τξ

ξτ

∑

∫ ∫

xn π l

2 l

πξ l

0 0

⎞ ⎟ ⎟ ⎠

⎛ ⎜ ⎜ ⎝

Bước 3: Đưa ra hàm Green

∞

t (

) τ −

an π l

⎛ − ⎜ ⎝

⎞ ⎟ ⎠

txG ;,(

sin

), τξ

∑

2 l

n πξ l

xn π l

1 =

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:

L(u) = f(x)

- 21-

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:

φ −=)(L

λσφ

có cùng điều kiện biên thuần nhất, hàm σ có thể tùy ý. Ta tìm nghiệm u(x) bằng cách khai triển vào chuỗi Fourier của các hàm riêng:

∞

xu )(

φ n

∑

1 =

Tác động toán tử L vào hai vế của đẳng thức trên, thu được

∞

)

xf )(

−=

La n

( φ n

σφλ n n

∑

1 =

Ta có các hàm riêng trực giao nhau theo công thức

≠

x )(

dxx )(

σφ n

φ n

∫

ndx ,

≠

2 σφ n

∫

⎧ ⎪ ⎨ ⎪ ⎩

∫

)( xf dx φ n

Suy ra

a − = λ n

2 σφ n

∫

∞

)( xu

,( xG

)(( f ξ

d ξ

)(( f ξ

) d ξξ

∑

∫

1 =

−

)( )( x ξφφ b 2 ∫ σφλ n a

∞

,( xG

) ξ

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à

∑ = −

)( )( x ξφφ b 2 ∫ σφλ n a

Trong đó

uxf );

(

)0(

lu )(

2 ud 2 dx

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

λ n

n π l

⎛ = ⎜ ⎝

⎞ ⎟ ⎠

Ta có các trị riêng và hàm riêng tương ứng là:

n =

1,2,...

- 22-

, với và = sin X x ( ) n n x π l

Ta có: d ξ ξ u x ( ) f G x ξ ( ) ( , ) = ∫

∞

∑

= 1

sin sin n πξ l G x = − ξ ( , ) 2 l

n π l n x π l ⎛ ⎜ ⎝ ⎞ ⎟ ⎠

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

nR còn u là hàm thuộc lớp

)

2 (

C Ω .Hàm ( )u x thỏa

Giả sử Ω là một miền trong

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

0u∆ =

C Ω sao cho

2 (

)

(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

u ∆ =

f x ( )

(2.4.2)

1B và giả sử

nR với biên ∂Ω thuộc lớp CΩ ∩ Ω . Công thức Gauss- Ostrogradsky :

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 Ω .

)

(

u x v x là các hàm thuộc lớp ( ), ( )

ds ,

υ j j

∑

∫

u ∂ ∂∑ x

= 1

Ω

∂Ω

Giả sử c là một miền bị chặn trong 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:

υ j

∫

v ∂ x ∂

u ∂ x ∂

u ∂ 2 x ∂ j

Ω

∂Ω

(2.4.3)

v udx ∆

= −

Lấy tổng đẳng thức (2.4.3) theo j từ 1 đến n ta nhận được công thức Green thứ nhất:

∑

∫

v ∂ x ∂

u ∂ x ∂

u ∂ υ ∂

= 1

Ω

∂Ω

(2.4.4)

v u u v dx ∆ − ∆

−

(

)

(

Đổ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

∫

u ∂ ∂ υ

v ∂ dS ) ∂ υ

Ω

∂Ω

(2.4.5)

- 23-

Các công thức (2.4.4) và (2.4.5) được sử dụng để nghiên cứu phương trình Laplaxơ và phương trình Poisson.

2 nr − đối với

2n > và 2n = , ở đây r là khoảng cách đến một điểm cố định. Ta cố định điểm

ln r đối vớii y ∈ Ω và đưa vào một nghiệm cơ bản chuẩn tắc của phương trình Laplaxơ:

−

n (2

1 n w − ) n

= Γ

−

x Γ − (

)

(

)

−

1 π 2

⎧ ⎪ ⎪ = ⎨ ⎪ ⎪⎩

Phương trình Laplaxơ có một nghiệm đối xứng xuyên tâm

nR .

nw là thể tích hình cầu đơn vị trong

ở đây

−

x Γ − (

)

(

)

y i

x i

n − −

−

x Γ − (

)

n x (

)(

δ ij

y i

x x i

{

} )

1 nw n 1 nw n

Qua một số phép tính ta nhận được

j= và

j≠ . Đương nhiên Γ là hàm điều hoà khi

ijδ = nếu i 1 y≠ . Trong trường hợp khi x

ijδ = nếu i 0 y= không thể thay thế hàm Green trong công thức x (2.4.5) bằng hàm Γ được. Tuy nhiên việc khó khăn này có thể khắc phục được nhờ việc là quả cầu tâm y bán kính ρ đủ nhỏ. Công thức thay thế Ω bằng

Ω

B B , ρ ρ

B y ( ) ρ

\ (2.4.5) khi đó có dạng

ở đây

∫

∂Ω

B ∂

Ω

(2.4.6) udx u ds u Γ∆ = − + − Γ ( Γ ( u ∂ υ ∂ ∂Γ ) υ ∂ u ∂ υ ∂ ∂Γ ds ) . υ ∂

− 1

n ρ

= Γ

≤

→

ρ (

) max

nw n

∫

u ∂ υ ∂

B ∂

ρ ( ) B ∂

uds

Hơn nữa

/ (

u x (

)

u y ( )

∫

B ∂

ρ ) B ∂

= −Γ = − → − khi ρ→ và 0 ∂Γ υ ∂

ρ→

∈ ∂

0, x

B ρ

Khi . Từ đó khi cho 0ρ→ trong đẳng thức (2.4.6) ta nhận được

công thức

−

u y ( )

)

x −Γ − (

)

x + Γ − ∆ (

y udx y )

∈ Ω .

(

∫

∂Γ ( ∂ υ

u ∂ ) ∂ υ

∂Ω

Ω

(2.4.7)

0u∆ = trong c thì từ (2.4.7) ta rút ra

Nếu

∫

∂Ω

- 24-

(2.4.8) u x y y = − ∈ Ω u y ( ) ( ) x − Γ − ( ) ds y , ∂Γ ( υ ∂ u ∂ ) υ ∂

2 (

Công thức (2.4.8) cho biễu diễn Green của hàm điều hoà thuộc lớp

C Ω tại ) điểm bất kỳ y ∈ Ω qua giá trị của ( )u x trên ∂Ω và giá trị của đạo hàm theo pháp tuyến u ∂ υ ∂

trên ∂Ω .

Bởi vì trong đẳng thức (2.4.8) các hàm dưới dấu tích phân là các hàm khả vi vô ( )u y cũng giải tích trong Ω . Như vậy các hàm hạn, hơn nữa giải tích theo y , nên hàm điều hoà giải tích trong toàn miền xác định của nó. Do đó chúng được xác định đơn trị nhờ các giá trị của mình trên một tập con mở bất kỳ của miền xác định. Tính chất đáng chú ý này của hàm điều hoà cũng đúng cho lớp các phương trình elliptic với các hệ số giải tích.

−

dx n ,

Tích phân dạng

u y ( ) 0

a x x ( ) 0

∫

Ω

(2.4.9)

a x trong Ω . Tích phân dạng 0 ( )

−

ds n ,

được gọi là thế vị khối hay thế vị Newton với mật độ

u y ( ) 1

a x x ( ) 1

∫

∂Ω

a x trên ∂Ω , còn tích phân dạng

(2.4.10)

−

x ∂ −

ds n ,

được gọi là thế vị lớp đơn với mật độ 1( )

u y ( ) 2

a x ( ) 2

∫

∂ υ

∂Ω

2 n

(2.4.12)

y −

được gọi là thế vị lớp kép với mật độ a x trên ∂Ω . Trong trường hợp n = 2 tương tự 2 ( ) ta cũng có các định nghĩa thế vị Newton hay logarit và các thế vị lớp đơn, thế vị lớp kép. Khi đó các công thức (2.4.9) (2.4.10) (2.4.11) cần thay hàm bằng hàm −

y − .

ln x

−

)

2 (

∂ và u trên ∂Ω .

C Ω có thể biểu diễn dưới dạng tổng của một thế vị lớp đơn và một thế vị lớp kép trên ∂Ω , mật độ của chúng được xác định bởi các giá trị

/u υ∂

a x . Thế vị lớp đơn (2.4.10) là thế vị của trường tĩnh điện trong 0 ( )

Từ công thức (2.4.8) suy ra một hàm điiều hoà thuộc lớp

a x . 2 ( )

−

Về ý nghĩa vật lý, gradien của thế vị Newton (2.4.9) xác định cường độ của 3 \R ∂Ω được tạo thành bởi điện tích phân bố trong Ω với mật trường tĩnh điện trong 3 \R ∂Ω được sinh độ ra bởi điện tích phân bố trên ∂Ω với mật độ a x . Gradien của thế vị lớp kép (2.4.11) 1( ) xác định cường độ của trường tĩnh điện được gây ra bởi ngẫu cực phân bố trên ∂Ω với mật độ mặt

> là hàm khả vi vô hạn theo x và

0x khi

x≠

−

x ∆ −

u ∆ = 1

a 1

∫

∂Ω

−

x ∆ −

u ∆ = 2

a 2

∫

∂ υ ∂

∂Ω

- 25-

Bởi vì , nên

)

u x 2(

)

\nR ∂Ω nếu a1 và a2 u x 1( 0 ( là các hàm thuộc lớp C ∂Ω . Như vậy các tích phân (2.4.10) và (2.4.11) xác định hai họ nghiệm của phương trình Laplace trong Ω . Cũng lý lụân như vậy ta nhận được thế vị Newton (2.4.9) là hàm điều hoà trong

\nR Ω nếu

và Do đó các hàm là các hàm điều hoà trong

(

a x C∈ 0( )

Ω ).

0u∆ = trong

(

)

Ω . Khi đó nhờ công thức Green thứ hai (2.4.5) ta nhận được

−

h udx ∆

u (

Bây giờ giả thiết hàm ∈ h C Ω ∩ Ω thoả mãn phương trình

∫

h ∂ ∂ υ

u ∂ ds ) ∂ υ

∂Ω

Ω

= Γ + ta nhận được biễu diễn Green

−

G udx ∆

u y ( )

(

∫

G ∂ ∂ υ

u ∂ ) ∂ υ

∂Ω

Ω

Cộng đẳng thức này với (2.4.7) và đặt G tổng quát hơn

∆

u y ( )

Nếu bổ sung G = 0 trên ∂Ω thì

∫ ds G udx

G ∂ ∫ u υ∂Ω ∂

Ω

(2.4.12)

G G x y ( , )

Hàm

CΩ ∩ Ω qua các giá trị biên của nó.

)

(

như thế được gọi là hàm Green ( của bài toán Dirichlet) đối với miền Ω . Đôi khi nó còn được gọi là hàm Green loại một đối với Ω . Như vậy việc tồn tại được một hàm Green kéo theo khả năng biểu diễn được một hàm điều hoà bất kỳ thuộc

Tiểu kết:

- 26-

Ở chương này đã xây dựng xong phương pháp hàm Green làm cơ sở cho việc áp dụng nó để giải bài toán truyền nhiệt ở chương sau. 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. Rồi biễu diễn nghiệm cần tìm thông qua 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

3.1 Thiết lập phương trình truyền nhiệt:

Nhiệt truyền từ nơi có nhiệt độ cao sang nơi có nhiệt độ thấp theo ba cách: quá trình dẫn nhiệt, quá trình bức xạ nhiệt và quá trình đối lưu. Quá trình dẫn nhiệt bên trong vật là do sự chuyển động của các phân tử bên trong vật. Trong vật rắn, dòng nhiệt chuyển từ nơi có nhiệt độ cao (là nơi có một số lớn các phân tử chuyển động có vận tốc lớn hay động năng lớn) sang nơi có nhiệt độ thấp hơn ( là nơi có vận tốc và động năng các phân tử nhỏ hơn). Quá trình bức xạ nhiệt giữa hai vật xảy ra khi nhiệt truyền qua không gian từ vật nóng hơn sang vật lạnh hơn ( không tính đến nhiệt độ không gian giữa hai vật), đó chính là chuyển động nhiệt dưới dạng sóng. Một ví dụ là sự truyền nhiệt độ của Mặt Trời cho Trái Đất. Nhiệt truyền do đối lưu xảy ra do một số loại chuyển động nhiệt di chuyển từ nơi này sang nơi khác. Cường độ của dòng đối lưu xảy ra khi cánh quạt thổi dòng nhiệt từ nơi này sang nơi khác. Có một loại truyền nhiệt khác sinh ra do bay hơi hoặc ngưng tụ. Tất cả các quá trình truyền nhiệt này được nghiên cứu trong các môn học đại cương và chuyên đề về nhiệt. Trong chương này chủ yếu tập trung nghiên cứu quá trình truyền nhiệt trong vật dẫn. Chúng ta nhắc lại định lý Gauss thường dùng để chuyển tích phân mặt sang tích phân 3 lớp.

r tzyxFF ,( ),

là một trường vectơ liên tục, xác định mọi nơi bên trong thể Nếu , =

r dFdiv

rr dnF . σ ,

tích V với bề mặt kín S bao quanh nó, thì theo định lý Gauss

∫∫=

∫∫∫

(3.1.1)

trong đó: τd là yếu tố thể tích và σd là yếu tố diện tích bề mặt; nr là pháp tuyến ngoài của bề mặt có độ dài bằng đơn vị.

r q

grad

r i

r j

r k

−=

uk −=∇−=

Sử dụng định lý Gauss, định luật Fourier về quá trình truyền nhiệt và định luật bảo toàn năng lượng để xây dựng phương trình truyền nhiệt, theo định luật Fourier về quá trình truyền nhiệt.

u ∂ x ∂

u ∂ y ∂

u ∂ z ∂

⎛ ⎜⎜ ⎝

⎞ ,⎟⎟ ⎠

(3.1.2)

zyxu ), ,(

biễu diễn nhiệt độ của vật. trong đó: qr là lượng nhiệt truyền qua một đơn vị diện tích trong một đơn vị thời gian; k là hệ số dẫn nhiệt phụ thuộc vào tính chất của vật liệu khi nhiệt truyền qua; hàm u =

zyxu ), ,( = const được gọi là mặt đẳng nhiệt. Ta thấy rằng, vectơ gradient trùng với pháp tuyến tại bất kỳ điểm nào trên bề mặt và hướng theo chiều tăng của nhiệt độ. Vì dòng nhiệt hướng từ nóng sang lạnh nên trong công thức (3.1.2) của định luật Fourier về quá trình truyền nhiệt lấy dấu trừ. Như vậy định luật Fourier về quá trình truyền nhiệt có thể được giải thích là dòng nhiệt truyền theo hướng tăng của nhiệt độ. Đại lượng vectơ qr được gọi là vectơ dòng nhiệt, bằng lượng nhiệt truyền qua một đơn vị diện tích.

Bề mặt có nhiệt độ không đổi

- 27-

Sử dụng các đại lượng nhiệt sau:

zyxc ,( ),

c =

là nhiệt dung của vật rắn;

zyxρρ= ,( ),

là mật độ khối lượng tính trên một đơn vị thể tích;

là hệ số dẫn nhiệt của chất rắn;

,( zyxk ), r là dòng nhiệt truyền qua một đơn vị diện tích; , k = r = tzyxqq ,( ),

là nguồn nhiệt tự sinh ra trên một đơn vị thể tích; , tzyxHH = ,( ),

, là nhiệt độ tại mọi điểm của vật. tzyxuu = ,( ),

t∆ . Định luật bảo toàn được viết dưới dạng

Viết định luật bảo toàn năng lượng cho một miền tùy ý V với bề mặt kín S bao t∆ . Hc là lượng t∆ . HG là lượng nhiệt sinh ra trong V quanh. Gọi HS là lượng nhiệt thay đổi trong V với khoảng thời gian nhiệt đi qua bề mặt S trong khoảng thời gian trong khoảng thời gian

c τρud

(3.1.3) H H H H H H = +⇒+ − 0=

t∆ có dạng

.τρ c

Lượng nhiệt có trong yếu tố thể tích τd của V và . HS là lượng nhiệt thay đổi trong V trong khoảng thời gian

∫∫∫

∂ t ∂

)

( ucρ , nói cách khác là thông

(3.1.4)

HC là lượng nhiệt đi qua bề mặt S trong thời gian

∫∫ −

lượng đi qua bề mặt S là . σvv dnq , (3.1.5)

trong đó dấu trừ để đổi dấu cho vectơ pháp tuyến ngoài có độ dài đơn vị là nv . Theo định lý Gauss, tích phân bề mặt được chuyển thành

dqdiv

∫∫∫−

(3.1.6) = .τv

.τ (3.1.7)

∫∫∫

Nhiệt lượng sinh ra trong V được cho bởi

(

r Hqdiv +

−

Kết quả từ các công thức (3.1.4), (3.1.6) và (3.1.7) cho phép viết định luật bảo toàn bởi phương trình

∫∫∫

∂ t ∂

⎧ ⎨ ⎩

⎫ ) τρ duc ⎬ ⎭

(3.1.8)

t∆ , như vậy số hạng trong dấu ngoặc {} phải bằng không. Thay biểu thức của vectơ qr vào phương trình (3.1.2) biểu thị định luật Fourier của quá trình truyền nhiệt ta thu được phương trình truyền nhiệt trong vật dẫn

- 28-

Kết quả trên cho một thể tích V tùy ý và thời gian tùy ý

div

(

grad

)

(

Hu +

uc ρ

∂ t ∂

(3.1.9)

(

uc ρ

Hoặc có thể viết dưới dạng mở rộng

∂ x ∂

∂ y ∂

∂ z ∂

∂ t ∂

⎛ ⎜ ⎝

⎞ +⎟ ⎠

⎛ ⎜ ⎝

⎞ +⎟ ⎠

⎛ ⎜⎜ ⎝

⎞ +⎟⎟ ⎠

(3.1.10)

(

)

(

)

div

grad

uku . ∆=∇=∇∇=

Trong trường hợp đặt biệt, nếu k là hằng số ta có

(3.1.11)

được gọi là toán tử Laplace. ∆=∇ 2

+∆

Khi các hệ số đều là hằng số, có thể viết phương trình truyền nhiệt dưới dạng

u ∂ t ∂

k c ρ

H c ρ

(3.1.12)

hệ số a được gọi là độ khuếch tán của vật liệu.

u ∂ lim = t ∂ t ∞→

thì có thể nói nhiệt độ ở trạng thái dừng hay ổn định. Nếu

u ∂ t ∂

0=Q

Trong trường hợp trạng thái dừng , trong phương trình (3.1.12) nhiệt độ

, phương chỉ phụ thuộc vào các vị trí bên trong. Nếu không có nguồn nhiệt, tức là trình truyền nhiệt trở thành phương trình thuần nhất. Ta có thể lập bảng sau cho phương trình truyền nhiệt trong hệ tọa độ Đề- các.

Các dạng khác nhau của phương trình truyền nhiệt

trong hệ tọa độ Đề- các.

Các trường hợp Dạng toán tử Dạng một chiều

( ) Huk

u ∂ ρ t ∂

+∇∇ = k H c = Tổng quát u ∂ x ∂ u ∂ ρ t ∂ ∂ x ∂ ⎛ ⎜ ⎝ ⎞ +⎟ ⎠

2 u +∇

+∇∇

( ) Huk

Vật liệu đồng = + = u 2 chất H k c ρ k u ∂ t ∂ H k c ρ k u ∂ t ∂ ∂ x ∂

2 u +∇

k H 0= Trạng thái dừng ∂ x ∂ u ∂ x ∂ ⎛ ⎜ ⎝ ⎞ +⎟ ⎠

H k

2 ud 2 dx

= + =

- 29-

Trạng thái dừng với vật liệu đồng chất

Các điều kiện ban đầu và điều kiện biên cho phương trình truyền nhiệt

Cho vật thể V với mặt S bao quanh, các điều kiện biên khác nhau có thể đặt trên

mặt biên S như sau:

1. Điều kiện biên Dirichlet hay bài toán biên loại I đòi hỏi nhiệt độ được xác định trên biên của miền, mà tại đó phương trình truyền nhiệt giải được. Loại điều kiện biên này có dạng

tzyxu ),,

(3.1.13)

tzyxf ,, 1

zyx , ,(

) ∈

trong đó

1f là nhiệt độ đã được xác định.

2. Điều kiện biên Neumann hay bài toán biên loại II đòi hỏi dòng nhiệt đi qua biên được xác định rõ trên biên của miền, mà tại đó phương trình truyền nhiệt giải được. Loại điều kiện biên này có dạng

r gradu n .

x y z t ( , , , )

(3.1.14)

S ∈

x y z ( , , )

x y z S ∈ ( , , )

S ∈

x y z ( , , )

Trong đó

u x y z t ∂ ( , , , ) n ∂ 2f là dòng nhiệt đã được xác định.

Đối với biên cách nhiệt thì

r gradu n .

(3.1.15)

Bien

u ∂ n ∂

Bien

3. Điều kiện biên Robin hay bài toán biên loại III đòi hỏi dòng nhiệt đi qua biên và nhiệ độ trao đổi với môi trường xung quanh được xác định rõ trên biên của miền, mà tại đó phương ttrình truyền nhiệt giải được. Loại điều kiện bhiên này có dạng

hu x y z t , , )

( ,

x y z t ( , , , )

(3.1.16)

S ∈

x y z ( , , )

S ∈

x y z ( , , )

u x y z t ∂ , ) ( , n ∂

Trong đó: h >0 là hằng số, 3f là dòng nhiệt đã được xác định.

4. Điều kiện biên hổn hợp là kết quả của các điều kiện biên loại I và II.

3.2 Bài toán biên phụ thuộc thời gian

3.2.1 Phương pháp tách biến Fourier cho bài toán truyền nhiệt

Tìm sự phân bố nhiệt trong thanh hữu hạn trên đoạn [0, L] nằm dọc theo trục x. Biết hai đầu x = 0 và x = L được giữ ở nhiệt độ bằng không, f x . phân bố nhiệt ban đầu của thanh là ( )

L u ( )

0, 0

x L .

−

< <

u ∂ 2 x ∂ ) =

(3.2.1)

u ∂ t ∂ ( u L t , = ( ) f x

) ) , 0

⎧ ⎪ ⎪⎪ ( u ⎨ ⎪ ( u x ⎪ ⎪⎩

- 30-

(3.2.2)

đây là bài toán Dirichlet xác định nhiệt độ tại các đầu mút của thanh và f(x) là nhiệt độ phân bố lúc ban đầu. Dùng phương pháp tách biến. Nghiệm tìm được có dạng: ( , txu

( ) ( )tTxX

)

thay (3.2.2) vào phương trình (3.2.1), ta được: ( ) ( )tTxXa ′′

( ) ( ) ′ tTxX

tiếp tục, chia hai vế của phương trình cho

, ta được:

( ) ( )tTxXa2

( ) ′ tT ( ) 2 tTa

( ) ′′ xX ( )xX

Const

=λ

trong đó, vế trái chỉ phụ thuộc vào t, còn vế phải chỉ phụ thuộc vào x, nghĩa là cho dù các biến số thay đổi, nhưng tỷ số luôn luôn bằng nhau. Điều đó chỉ xảy ra khi tỷ số này là một hằng số và được chọn là λ− , với .Tức là:

λ−= ,

( ) ′ tT ( ) 2 tTa

( ) ′′ xX ( ) xX

ta nhận được hai phương trình vi phân sau:

≠

(3.2.3)

≠

(3.2.4)

( ) +′ tT ( ) ′′ xX

( ) 2 tTa λ ( ) xX

( ) tT ( ) xX

Sử dụng điều kiện biên ban đầu (3.2.1) cho ta:

( ) ( ) tTX 0 = ( ) ( ) tTLX

⇒

(3.2.5)

0 0

= =

( ) u ,0 ⎧ ⎨ ( ) tLu , = ⎩ ( ) 0≠tT ( ) X 0 ( ) LX

⎧ ⎨ ⎩

Xét phương trình (3.2.4):

,0 = λλ <

Nếu

, ta có: u (x,t)=0 ( loại)

2c=λ

0>λ , đặt

, ta có nghiệm không tầm thường:

cos

sin

Nếu ( ) xX

theo điều kiện ban đầu (3.2.5), ta có

= sin

( ) 0 ( ) LX

⎧ ⎨ ⎩

0B⇒ ≠ và phương trình tìm trị riêng: sincL=0

- 31-

c =⇒

n π L

Do đó, bài toán chỉ có nghiệm không tầm thường khi giá trị riêng:

λλ = n

n π L

⎛ ⎜ ⎝

⎞ ⎟ ⎠

với: n=1,2,3…..

Với mỗi trị riêng có một hàm riêng tương ứng được viết là:

sin

)

( X x n

n π L

Xét phương trình (3.2.3):

′

22 ca

−=

)22 ( ca

T T

22 tca

⇒

−=

22 tca

⇒

−=

22 tca

−

⇒

T T A ( ) tT n

eA n

với

nA là hằng số tuỳ ý.

Nghiệm riêng của phương trình truyền nhiệt (3.2.2) là:

−

an π L

⎛ ⎜ ⎝

⎞ ⎟ ⎠

sin

(3.2.6)

)

( ) ( ) tTxX

( , txu n

eA n

xn π L

Nghiệm tổng quát là tổng của tất cả các nghiệm riêng ứng với các giá

trị khả dĩ của n

∞

−

an π L

⎞ ⎟ ⎠

⎛ ⎜ ⎝

sin

( txu ,

)

( txu ,

(3.2.7)

eA n

∑

) ∑ =

xn π L

1 =

Điều kiện ban đầu (3.2.1) cho ta xác định tuỳ ý

nA . Ta có:

∞

sin

( xu

)

(3.2.8)

A n

( ) ∑ xf =

xn π L

1 =

Theo lý thuyết chuỗi Fourier, có thể tìm được các hệ số nA là: ( Xf ,

)

sin

( ) xf

(3.2.9)

A n

n 2

∫

2 L

xn π L

- 32-

3.2.2 Phương pháp hàm Green cho bài toán truyền nhiệt

Tìm sự phân bố nhiệt trong thanh hữu hạn trên đoạn [0, l] nằm dọc theo trục x. Biết hai đầu x = 0 và x = l được giữ ở nhiệt độ bằng không, phân bố nhiệt ban đầu của thanh là

( )xϕ và trong thanh có nguồn nhiệt Q (x, t)

Giải:

x < <

Q x t ( , )

u l t

( , ) 0

(3.2.10)

u x

= =

xϕ ( )

⎧∂ u u ∂ ⎪ 2 x t ∂ ∂ ⎪⎪ t u (0, ) 0, ⎨ ⎪ ( , 0) ⎪ ⎪⎩

Áp dụng phương pháp mở rộng hàm riêng ta chọn nghiệm có dạng

∞

sin)(

),( txu

tu n

∑

xn π l

1 =

∞

sin)(

),( txQ

tq n

∑

xn π l

1 =

Thay vào phương trình (3.2.10), tìm nghiệm:

t )(

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

tu )( n

tq )( n

n dt

an π l

⎛ ⎜ ⎝

⎞ ⎟ ⎠

Nghiệm có dạng :

an π l

⎛ − ⎜ ⎝

⎞ ⎟ ⎠

⎛ − ⎜ ⎝

⎞ ⎟ ⎠

⎛ ⎜ ⎝

⎞ ⎟ ⎠

)0(

)( τ

d τ

)( tu n

∫

Dựa vào điều kiện ban đầu tìm hàm un(0)

∞

)( x

sin)0(

∑

xn π l

1 =

)0(

sin)(

⇒

ξϕ

d ξ

∫

n πξ l

2 l

∞

sin)(

),( txQ

tq n

∑

xn π l

1 =

⇒

sin),( τξ

d ξ

)( tq n

2 ∫= l

n πξ l

Cuối cùng ta thu được

- 33-

∞

−

an π l

⎛ ⎜ ⎝

⎞ ⎟ ⎠

txu ),(

sin)(

ξϕ

n πξ l

⎞ ⎟ d ξ ⎟ ⎠

⎡ ⎛ 2 ∑ ∫ ⎢ ⎜ ⎜ l ⎢ ⎝ ⎣

−

an π l

⎞ ⎟ ⎠

⎛ ⎜ ⎝

⎞ ⎟ ⎠

⎛ ⎜ ⎝

sin

sin),( τξ

∫

2 l

n πξ l

xn π l

⎛ ⎜ ⎜ ⎝

⎞ ⎟ d ξ ⎟ ⎠

⎤ ⎥ d τ ⎥ ⎦

Đổi thứ tự giữa tổng và tích phân ta thu được

∞

−

an π l

⎛ ⎜ ⎝

⎞ ⎟ ⎠

( txu ),

sin

( ) ξϕ

∑

∫

2 l

πξ l

xn π l

⎛ ⎜ ⎜ ⎝

⎞ ⎟ ⎟ ⎠

∞

)

( t

−

an π l

⎛ ⎜ ⎝

⎞ ⎟ ⎠

sin

( , ) τξ

ξτ

∑

∫ ∫

2 l

πξ l

xn π l

0 0

⎛ ⎜ ⎜ ⎝

⎞ ⎟ ⎟ ⎠

Hàm Green có dạng:

∞

t (

) − τ

an π l

⎛ − ⎜ ⎝

⎞ ⎟ ⎠

txG ;,(

sin

), τξ

∑

2 l

n πξ l

xn π l

1 =

Nghiệm của phương trình (3.3.1) có thể viết dưới dạng

d d

u x t ( , )

ϕξ ( )

G x t ( , ,

d ξ ξ ,0)

G x t ξτ ( , ) ( , ,

ξ τ ξ .

, 0)

∫

∫ ∫

0 0

Như vậy chúng ta đã tìm được sự phân bố nhiệt trên thanh, nghiệm này được biễu diễn thông qua hàm Green. Phương pháp tìm nghiệm được trình bày ở trên là phương pháp hàm Green. Ta nhận thấy phương pháp này có phần đơn giản vì không trực tiếp giải phương trình vi phân không thuần nhất. Đối với bài toán này, để tìm sự phân bố nhiệt độ trong thanh bằng phương pháp tách biến Fourier ta phải đi giải phương trình vi phân không

Q x t ( , )

. Nghiệm của phương trình này là

thuần nhất dạng:

u ∂ t ∂

u ∂ 2 x ∂

u(x, t) chính là sự phân bố nhiệt độ trong thanh. Để giải được phương trình đó phải qua nhiều bước:

- Dùng phương pháp tách biến tìm nghiệm của phương trình thuần nhất ở dạng u(x,t) = X(x).T(t). Thay nghiệm này vào phương trình thuần nhất một chiều từ đó di đến việc giải các phương trình vi phân.

- Sử dụng các điều kiện biên để suy ra dạng của nghiệm u(x,t) của

phương trình thuần nhất.

- Cho nghiệm u(x,t) của phương trình thuần nhất thoả mãn điều kiện ban đầu, kết hợp với việc áp dụng chuỗi Fourier khai triển hàm

Q x t ( , )

q x t ( , )

1 cρ

- 34-

- Thay vào phương trình truyền nhiệt không thuần nhất, từ đó mới suy

ra nghiệm cần tìm u(x,t).

Như vậy phương pháp hàm green giúp chúng ta giải bài toán này đơn

giản hơn, ngắn gọn hơn.

3.2.3 Bài toán truyền nhiệt trong miền tròn

Xác định sự phân bố nhiệt độ trong miền tròn bàn kính R, tâm O.

Miền tròn D có biên Γ .

Giải:

Ở đây chúng tôi sẽ giải bài toán này bằng hai phương pháp: tách

biến Fourier và phương pháp hàm Green.

Cách 1: Phương pháp tách biến Fourier:

trong đó

tìm hàm ( , ) điều kiện

f θ ( )

= r R

Giả sử hình tròn bán kính R với tâm tại cực O của hệ toạ độ cực. Ta u r θ điều hoà trong hình tròn, và trên vòng tròn của nó thoả mãn f θ là hàm cho trước, liên tục trên vòng ( ) tròn. Hàm cần tìm phải thỏa mãn trong hình tròn phương trình Laplaxơ.

(3.2.11)

∂ u 2 ∂ r

∂ ∂ u u + 2 r θ ∂ ∂

u Q r T θ

( ) ( ).

Giả sử nghiệm riêng được tìm dưới dạng:

Khi đó ta được:

2 // r Q r T

/ rQ r T

θ ( ). ( )

Q r T ( ).

θ = ( ) 0

Tách biến:

/ rQ r ( )

= −

T T

θ ( ) θ ( )

2 // r Q r + ( ) Q r ( )

Cho mỗi vế của đẳng thức nhận được bằng hằng số

2k− , ta có hai

phương trình vi phân thường

2 k T

θ ( )

= θ ( ) 0,

2 k Q r

−

2 // r Q r ( )

/ rQ r ( )

= ( ) 0.

k = ta nhận được

Từ đó, khi

A B

= +

θ ( )

(3.2.12)

Q r C D r = + ln .

( )

(3.2.13)

B k

k + θ

θ ( )

cos

θ sin ,

(3.2.14)

Nếu

k > , thì 0

và

−

rmr

2 k r

r 1)

(

) 0,

hay

,m Q r ( ) 2 = tức là

được m

Còn nghiệm phương trình thứ hai sẽ tìm dưới dạng − 1 2 mr m k− 2 − − r m m ( k= ± Do đó .

- 35-

Q r Cr Dr −

( )

(3.2.15)

Nhận xét rằng

u r θ như là hàm của θ là hàm tuần hoàn chu kỳ ( , )

u r θ π+ 2 )

( ,

u r θ và ( , )

2π, vì đối với hàm đơn trị các đại lượng trùng nhau. Bởi vậy , từ (3.2.12) suy ra B=0, còn trong (3.2.14) k có thể nhận một k > ). Tiếp theo, trong các đẳng thức (3.2.13) và trong các giá trị 1, 2, 3,…( (3.2.15) phải cho D=0, nếu không hàm số sẽ gián đoạn tại điểm r =0 và do đó không điều hoà trong hình tròn. Vậy ta nhận được vô số nghiệm riêng của phương trình (3.2.11), liên tục trong hình tròn. Chúng có thể viết dưới dạng:

u r θ 0 ( , )

n θ

θ ( , )

(

cos

sin

n θ )

(

1,2,...)

u r n

B n

A 0 2 A n

Bây giờ ta lập hàm

∞

u r

n θ

θ ( , )

(

cos

sin

n θ )

A n

B n

∑

A 0 2

= 1

Do tính tuyến tính và thuần nhất của phương trình Laplaxơ, hàm này cũng là nghiệm của nó. Ta còn xác định A0, An, Bn để hàm đó thoả mãn điều kiện

f θ ( )

r R =

Tức là

∞

n θ

θ ( )

(

cos

sin

n θ )

A n

B n

∑

A 0 2

= 1

Ở đây ta khai triển hàm

f θ thành chuổi Fourier trong đoạn [ ( )

] ,π π−

Theo công thức đã biết ta có:

τ ( ) cos

n d τ τ ,

d A τ τ ( ) , n

A 0

∫

1 R π

1 π

− π

τ ( )sin

n d τ τ .

B n

∫

1 R π

− π

Như vậy

∞

u r

−

θ ( , )

τ ( ).

n .cos (

τ θ τ ) .

∑

∫

r R

1 π

1 2

= 1

⎛ ⎜ ⎝

⎞ ⎟ ⎠

− π

⎡ ⎢ ⎢ ⎣

⎤ ⎥ ⎥ ⎦

ρτ θ

− = và biểu diễn biểu thức trong ngoặc vuông dưới

Đặt

1 R

dạng

- 36-

∞

n ρ

−

cos

∑

1 2

= 1

Xét chuổi

∞

e ρ

n ρ

cos

sin

(

)

∑

Chuổi này hội tụkhi

1ρ< và tổng của nó bằng

1 t

ρ −

i ρ t +

−

sin 2 ρ

i ρ

1 e ρ

+ − cos ρ 1 2 cos

cos

sin

Do đó

∞

n ρ

cos

∑

2 ρ

t cos t +

−

1 − = 2

− 1 ρ 1 2 cos

1 − = 2

2 − ρ 1 ρ 2(1 2 cos

)

Hay trở về ký hiệu cũ ta được:

u r θ ( , )

τ ( ).

d τ .

∫

−

1 π 2

r R − Rr τ θ − cos( )

π −

u r θ là sự phân bố nhiệt độ trong miền tròn ( , )

Ta tìm được nghiệm thoả mãn điều kiện bài toán.

Cách 2: Dùng phương pháp hàm Green:

Ta cần tìm hàm u thoả mãn phương trình

0u∆ = trên miền D có

biên Γ .

O, R là tâm, bán kính của Γ .

f x y ( , )

∆ = u

vì

Hàm u thoã mãn điều kiện biên

= g x y ( , ) 0,

(

)

nên áp dụng công thức Green cho u tại điểm M0, ta được:

= −

u M (

)

f x y ( , )

∫

G ∂ n ∂

(

)

= −

thì

Ta chọn

0v∆ = và

1 1 ln / rπ 2

= −

(

)

(

)

1 / r

−

(ln

)

1 r

1 / r

G ∂ n ∂

1 π 2 ∂ n ∂

1 π 2

( )Γ

= −

1 r

Ta tính được

= −

(ln

)

r r r n cos( , )

1 (ln ) r 1 / r

1 / r

∂⎧ ⎪⎪∂ n ⎨ ∂⎪ ⎪∂⎩ n

- 37-

r r r n cos( , )

Mặt khác:

2 ρ

với

r r r n cos( , )

0OMρ=

− + r Rr 2

2 ρ

cos(

r − + r ρ 2

2 ρ

)

Vậy

R 2

− Rr

G ∂ n ∂

1 ( π 2

(

)

Để xác định u tại điểm M0 trên miền tròn, ta có công thức

)

u M (

)

f x y ds ( , )

2 ρ − 2 r

1 ∫ π Γ R 2

(

)

⎛ ⎜ ⎝

⎞ ⎟ ⎠

Với toạ độ cực, ta có:

M R

ρϕ ( , ),

θ ( , )

2 ρ

−

ρ θ ϕ )

cos(

0 r = ds Rd =

+ θ

Do đó trong toạ độ cực, ta có

π 2

2 ρ

ρϕ ( , )

d θ θ ( )

∫

−

1 π 2

R − 2 − ρ θ ϕ ρ cos(

)

Cả hai phương pháp trên chúng ta đều đi đến cùng một kết quả, tìm được sự phân bố nhiệt độ trong miền tròn là như nhau. Nhưng ta nhận thấy đối với bài toán truyền nhiệt trên miền tròn thì giải bằng phương pháp hàm Green là đơn giản, tìm được nghiệm hiệu quả hơn, nhanh hơn.

r r / r n , )

3.3 Bài toán biên truyền nhiệt dừng

Bài toán: Tìm sự phân bố nhiệt độ ở trạng thái dừng trong thanh hữu hạn trên đoạn [0,L].Thanh có nguồn nhiệt Q(x,t) = Q(x). Các điều kiện biên có dạng: u(0) = 0,u(L) = 0

Giải:

Ở trạng thái dừng, ta có

Q x t Q x

( , )

( );

u ∂ t ∂

Và

Q x ( )

u ∂ t ∂

u ∂ 2 x ∂

Q x

⇔ = 0

( );

= −

⇒

f x ( )

2 d u 2 dx

u ∂ 2 x ∂ Q x ( ) 2 a

- 38-

u L

(0)

0; (

)

= 0

f x ( )

bằng

Ta sử dụng phương pháp biến thiên tham số cho phương trình

2 d u 2 dx

. Trong đó toán tử L có dạng:

cách xét bài toán không thuần nhất tổng quát

L u ( )

f x ( )

q .

ta đang xét bài toán truyền nhiệt dừng nên p =1, q = 0, ta có hai

d dx

du dx

⎛ ⎜ ⎝

⎞ ⎟ ⎠

L x

= −

x u 2;

u nghiệm của phương trình thuần nhất tương ứng là 1 và x. Nếu chọn: 1 vì

= − nên thu được:

−

ξ ( )(

d ξ ξ )

v x ( ) 1

c 1

∫

1 L

= −

d ξξ ξ ( )

v x ( ) 2

c 2

∫

1 L

Đây là hai công thức cần thiết trong phương pháp biến thiên hằng số

(

)

. Từ điều kiện biên suy ra:

u u v = 1 1

u v 2 2

= → =

);

u (0) 0

;

u u v = ( 1 1

u v 2 2

c L 2

= → =

−

u L ( ) 0

ξ ( )(

d ξ ξ )

c L 1

∫

Nghiệm của bài toán biên không thuần nhất là

= −

−

u x ( )

ξ ( )(

d ξ ξ )

d ξξ ξ ( ) ,

∫

L x − L

x L

f G x ξ ( )

( , )

d ξ ξ ,

u x Hay ( )

= ∫

ξ )

−

G x

ξ ( , )

Trong đó

ξ (

)

x L − ( L L x − L

⎧ ⎪⎪ = ⎨ ⎪− ⎪⎩

Ta nhận thấy tính đối xứng của hàm Green

G x

ξ=

ξ ( , )

( , ).

- 39-

Bài toán: Xét phân bố dừng của nhiệt độ trong quả cầu đồng nhất bán kính q với

điều kiện nữa trên được giữ ở nhiệt độ không, nữa dưới ở nhiệt độ 1.

Giải:

Bài toán này ta phải giải phương trình Laplaxơ với điều kiện biên:

khi

< < θ

π 2

θϕ ( , )

khi

< <

θ π

⎧ ⎪⎪ = ⎨ ⎪ 1 ⎪⎩

π 2

Trước hết ta cần đưa ra hàm Green

Ta xét vùng V là quả cầu bán kính q, tâm ở góc toạ độ, P0 là một điểm bất kỳ trong đó. Bao quanh V là mặt cầu S0. Kí hiệu V0 là phần của vùng V nằm ngoài S0. Do đó V0 được giới hạn bởi hai mặt S và S0. Công thức Green đối với vùng v0 là:

v u u v dV ∆ − ∆

−

)

(3.3.1)

∫

u ∂ n ∂

v ∂ n ∂

u ∂ n ∂

v ∂ n ∂

⎛ ⎜ ⎝

⎞ dS ⎟ ⎠

⎛ ⎜ ⎝

⎞ dS ⎟ ⎠

∫ (cid:0)

V 0

S 0

Trong công thức này ta xem u là nghiệm của bài toàn Đirichlet, còn v được chọn

là hàm Green G(P) xác định như sau:

G P (

)

H p ( )

r P P 0

Trong đó

0P Pr

là khoảng cách giữa P0 và một điểm biến thiên P(x, y, z), H(P) là một hàm thoã mãn phương trình Laplaxơ trong

−

đối với các điểm

vùng V và nhận giá trị

1 P Mr

M của mặt S

= −

H M (

)

r P M 0

G M (

)

H M (

= ) 0

Do đó

r P M 0

Bây giờ ta lấy điểm

0P nằm trên một tia đi từ góc toạ độ, qua điểm P0, sao

trong đó

là khoảng cách từ

cho

r 0

2 0

2 y 0

2 z 0

* r 0

0 0r r

0P đến tâm quả cầu.

Vậy ta có

* P 0

z 0

q 2 r 0

⎛ ⎜ ⎝

⎞ ⎟ ⎠

- 40-

> ). Các điểm P0 và

* r 0

0P nằm ngoài quả cầu( bởi vì r0 < q, nên

0P là đối

q r 0

xứng đối với mặt cầu S, giới hạn quả cầu V.

P0 và

M là một điểm bất kỳ trên mặt S, ta chứng minh tỉ số các khoảng cách từ M đến 0P là một đại lượng không đổi, không phụ thuộc vào M. Nếu kí hiệu P0M và

và

, ta có

0P M qua

0P Mr

r * 0P M

* r 0 q

q r 0

Từ đó ta rút ra

const

(3.3.2)

q r 0

r * P M 0 r P M 0

Ta chọn hàm Green cho quả cầu V là hàm

−

G P ( )

q r r 0

r P P 0

* P P 0

Ta chỉ cần kiểm nghiệm lại là

P0 γ

0θ

= −

H P (

)

q r r 0

* P P 0

Thoả mãn phương trình Laplaxơ trong quả cầu

= −

V và

r P M 0

Vì

0P nằm ngoài quả cầu V, nghĩa là H(P) được xác định ở tất cả các điểm bên

∆

trong V và do đó

r * P P 0

⎛ ⎜ ⎜ ⎝

⎞ ⎟ ⎟ ⎠

= −

H M (

)

Mặt khác:

q r r 0

r P M 0

* P M 0

trên mặt cầu S. Bởi vì đạo hàm theo pháp

Đạo hàm theo pháp tuyến ngoài

G ∂ n ∂

tuyến ngoài ở điểm M của mặt cầu S trùng với đạo hàm theo phương bán kính;

- 41-

G ∂ n ∂

G ∂ r ∂

r q =

Nên ta chuyển sang toạ độ cầu. Giả sử toạ độ cầu của điểm P là

r θϕ của điểm ,

;

θ ϕ . , , 0

P0 là 0

r θ ϕ khi đó ; 0

0P sẽ có toạ độ cầu là

q r 0

và

là:

r Các vectơ đơn vị theo các phương OP

r i θ ϕ

r j θ ϕ

r k θ

sin cos

sin sin

r i

cos r j

r k

Và

nên ta có thể tìm góc γ giữa

sin

cos

θ ϕ cos 0

θ ϕ sin 0

θ 0

r 0OP

0 qua tích vô hướng của hai vectơ đơn vị

cos

sin sin

cos cos

sin sin

cos cos

θ θ ϕ ϕ 0

θ θ 0

sin sin

+ ) cos cos

(3.3.3)

ϕ ϕ 0

θ θ 0

cos cos

sin sin

cos(

sin sin (cos cos 0 θ θ 0

θ θ ϕ ϕ 0 θ θ 0

− ϕ ϕ 0)

−

cos

ta rút ra

Từ công thức

P Pr

2 r 0

rr 02

−

= −

P P 0

r cos 0 3 r P P 0

⎛ ∂ ⎜ ⎜ r r ∂ ⎝

⎞ ⎟ ⎟ ⎠

cos

= −

Và

∂ r ∂

q r − 0 −

(

3 γ cos )

r P P 0

q r − 0 3 r P P 0

2 r 0

qr 0

⎛ ⎜ ⎜ ⎝

⎞ ⎟ ⎟ ⎠

r q =

Nếu thay thế 0r cho

q r , ta tìm được

−

cos

q r 0

= −

∂ r ∂

r * P P 0

⎛ ⎜ ⎜ ⎝

⎞ ⎟ ⎟ ⎠

r q =

−

(

3 γ cos )

q 2 r 0

q r 0

−

cos

= −

2 r 0 2 q

−

(

3 γ cos )

qr 0

r 0 2 r 0

= −

)

Từ công thức Green v đối với vùng V0, ta tìm được

u P ( 0

f M ( 1

G ∂ n ∂

1 π 4

∫(cid:0)

Đối với quả cầu có tâm ở góc toạ độ, ta có

- 42-

r r OP va OP 0

−

)

(

u P ( 0

f M dS ) 1

1 q π 4

∫(cid:0)

−

(

3 γ cos )

2 r 0 qr 0

2 r 0

Hay

ππ 2

−

d d

)

)sin

θϕ θ θ ϕ ( ,

(3.3.4)

u P ( 0

∫ ∫

q π 4

−

(

3 γ cos )

0 0

2 r 0 qr 0

2 r 0

vaθ

ϕ trên mặt S giới hạn bởi

)

Trong đó hàm 1(

f M là hàm của các toạ độ cầu 2

quả cầu V kí hiệu là

θϕ ( , ),

sin

θ θ ϕ d d ,

còn cosγ được tính bằng (3.3.3).

khi

θ < <

π 2

θϕ ( , )

Sử dụng các điều kiện biên

khi

< <

θ π

⎧ ⎪⎪ = ⎨ ⎪ 1 ⎪⎩

π 2

ta được

ππ 2

−

)

sin

d d θ θ ϕ

u P ( 0

∫ ∫

q π 4

−

(

3 γ cos )

2 r 0 qr 0

2 r 0

π 2

Ta tìm phân bố nhiệt độ trên bán kính với

0θ = và

0θ π= .

Khi

và

0, cos

cos

θ 0

ππ 2

−

d d θ θ ϕ

)

sin

u P ( 0

∫ ∫

q π 4

−

(

3 θ cos )

2 r 0 qr 0

2 r 0

π 2

−

(3.3.5)

q 2

−

2 r 0 2

cos

qr 0

2 r 0

qr 0

⎡ ⎢ ⎢ ⎣

= θ π ⎤ ⎥ ⎥ ⎦ = θ

π 2

2 r 0

−

1 2

− r 0

−

2 r 0

⎛ ⎜ ⎜ ⎝

⎞ 1 ⎟ ⎟+ q r ⎠ 0

và

Còn khi

= θ π

, cos

cos

2 r 0

−

)

(3.3.6)

u P ( 0

1 2

− r 0

1 q r − 0

2 r 0

⎛ ⎜ ⎜ ⎝

⎞ ⎟ ⎟ ⎠

r → cho ta nhiệt độ ở tâm quả cầu.

Cả hai công thức này khi 0

- 43-

)

ta tìm được nhiệt độ của nó ở giữa bán kính thẳng đứng phía trên

= =

u P ( 0

r 0

1 2

trong công thức (3.3.6))

r trong công thức (3.3.5) và phía dưới ( 0

( 0 r

1 q= 2

− = −

)

u P ( 0

qtren

r 0

1 2

3 2 5

⎛ −⎜ 1 ⎝

⎞ ⎟ ⎠

1 2

)

u P ( 0

q duoi

r 0

3 = − 2

1 = + 2

3 2 5

⎛ −⎜ 1 ⎝

⎞ ⎟ ⎠

1 2

Tiểu kết:

Ở chương này chúng tôi đã áp dụng phương pháp hàm Green để giải một số bài toán truyền nhiệt. Qua đó cho thấy khi sử dụng phương pháp hàm Green thì việc giải các bài toán này là đơn giản và tìm được nghiệm nhanh hơn.

- 44-

PHẦN III: KẾT LUẬN Chương I đã trình bày các cơ sở toán học cơ bản, cần thiết cho việc xây dựng phương pháp hàm Green như: Bài toán biên để sử dụng cho phương trình toán lý, khái niệm toán tử, hàm riêng, trị riêng để ứng dụng vào các chuỗi, một số phương pháp như phương pháp tách biến, phương pháp biến thiên tham số. Từ đó ở chương II tiến hành xây dựng phương pháp hàm Green. Để xây dựng phương pháp hàm Green chúng tôi đã đi 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:

x < <

Q x t ( , )

u l t ( , ) 0

= =

xϕ ( )

⎧∂ u u ∂ ⎪ 2 x t ∂ ∂ ⎪⎪ t u (0, ) 0, ⎨ ⎪ u x ( , 0) ⎪ ⎪⎩

và tìm được hàm Green có dạng

∞

t (

) − τ

an π l

⎛ − ⎜ ⎝

⎞ ⎟ ⎠

txG ;,(

sin

), τξ

∑

2 l

n πξ l

xn π l

1 =

Tuy nhiên đối với mỗi bài toán truyền nhiệt cụ thể thì hàm Green sẽ có biểu thức

cụ thể khác nhau.

Vậy kết thúc chương II ta đã xây dựng xong phương pháp hàm Green và nêu lên được tính chất của hàm Green. Trên cơ sở đó, ở chương III chúng tôi đã áp dụng phương pháp hàm Green để giải một số bài toán truyền nhiệt, bên cạnh đó chúng tôi cũng đã trình bày cách giải các bài toán này bằng phương pháp khác để so sánh và thấy được rằng dùng phương pháp hàm Green thì việc giải các bài toán này là đơn giản và tìm được nghiệm nhanh hơn.

Do đặc trưng của mỗi bài tập phương trình truyền nhiệt là khá dài nên trong giới

hạn của đề tài, chúng tôi chưa đưa ra được nhiều bài tập.

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 trung học phổ thông. Học tốt học phần này người học sẽ có những bước đi vững chắc khi học các học phần tiếp theo như: cơ học lượng tử, điện động lực,…Khoá luận này đã bổ sung một phương pháp giải hiệu quả cho bài toán truyền nhiệt trong học phần phương pháp toán lý, từ đó giúp sinh viên học tốt hơn học phần này.

Hiện tại khoá luận chỉ dừng lại ở việc sử dụng phương pháp hàm Green để giải một số bài toán truyền nhiệt. Đây là một loại bài tập cơ bản trong học phần phương pháp toán lý. Nếu có thể thì trong tương lai, khoá luận sẽ không dừng lại ở một số bài mà mở rộng ra cho tất cả các dạng phương trình truyền nhiệt. Nếu tiến xa hơn nữa thì áp dụng phương pháp hàm Green để tìm nghiệm của phương trình truyền sóng.

Hy vọng khoá luận sẽ được phổ biến và là tài liệu tham khảo hữu ích cho các bạn

sinh viên ngành sư phạm vật lý khi học đến học phần phương pháp toán lý.

- 45-

PHỤ LỤC 1 CÁC KHÁI NIỆM CƠ BẢN VỀ PHƯƠNG TRÌNH

ĐẠO HÀM RIÊNG :

nhiều biến với các đạo hàm riêng theo các biến này.

Phương trình đạo hàm riêng là một phương trình chứa hàm cần tìm của hai hoặc

phương trình.

Thí dụ:

là phương trình đạo hàm riêng cấp hai.

−

u ∂ x y ∂ ∂

Cấp của phương trình đạo hàm riêng là cấp cao nhất của đạo hàm có mặt trong

phương trình

Nghiệm của phương trình đạo hàm riêng là một hàm, nó thỏa mãn đồng nhất

Nghiệm tổng quát là nghiệm có chứa số hàm tùy ý độc lập bằng số cấp của phương trình (khác với phương trình vi phân thường, nó có nghiệm phụ thuộc vào hằng số)

chọn thích hợp hàm tùy ý

Thí dụ: bằng cách thế vào phương trình ta thấy:

là nghiệm của PTDHR trong thí dụ trên.

u x y ( ,

)

x y

+F(x)+G(y)

−

1 2

Nó chứa hai hàm độc lập tùy ý F(x) và G(y), vậy nó là nghiệm tổng quát. Trường hợp riêng riêng: ta

nghiệm

được

một

F(x)=sin x; G(y)=2y +3

u x y ( ,

)

x y

+sin x+2y +3

−

1 2

Nghiệm riêng là một nghiệm có thể nhận được từ nghiệm tổng quát bằng cách

chọn thích hợp hàm tùy ý.

Nghiệm đặc biệt là nghiệm không thể nhận được từ nghiệm tổng quát bằng cách

Ở đây ta chỉ xét các phương trình đạo hàm tuyến tính cấp hai.

Phương trình đạo hàm riêng tuyến tính cấp hai của hàm hai biến u(x,y) có dạng:

(1.1)

Fu G =

u ∂ x ∂

u ∂ x y ∂ ∂

u ∂ y ∂

u ∂ x ∂ trong đó A, B, …., G có thể là hàm của x,y nhưng không phụ thuộc u . Phương

trình cấp hai của hàm hai biến không có dạng nêu trên thì ta gọi là hàm phi tuyến.

Nếu G = 0, phương trình gọi là thuần nhất, nếu G≠ 0 thì ta gọi là phương trình

không thuần nhất. Điều này có thể tổng quát hóa cho phương trình cấp cao hơn.

ta phân loại phương trình đạo hàm riêng :

Tùy thuộc vào dấu của

−

- 46-

Bài toán biên của phương trình đạo hàm riêng là bài toán tìm kiếm các nghiệm của phương trình trong miền xác định nào đấy thỏa mãn các điều kiện trên biên của miền, gọi là bài toán biên. Định lý liên quan đến tồn tại và duy nhất nghiệm như vậy của bài toán gọi là định lý tồn tại và duy nhất.

>0 – phương trình loại Hyperbolic

−

<0 – phương trình loại Eliptic.

−

=0 – phương trình loại parabolic

−

- 47-

PHỤ LỤC 2 PHƯƠNG TRÌNH VI PHÂN CẤP HAI:

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

....

( ) C y x 1 1

( ) C y x 2

( ) C y x n

}

y (x),y (x) là tập nghiệm cơ bản của phương trình vi phân tuyến tính 1

Giả sử { thuần nhất cấp hai:

L y ( )

0 =

a x ( ) 0

a x ( ) 1

a x y ( ) 2

dy dx

2 d y 2 dx

( )

Suy ra nghiệm tổng quát có dạng : cy

1 1

2 2

( ) C y x C y x +

(2.1)

trong đó:

,C C là các hằng số tuỳ ý.

Dùng phương pháp biến thiên hằng số tìm nghiệm riêng của phương trình vi phân

không thuần nhất:

(2.2)

L y ( )

( )

a x ( ) 0

a x ( ) 1

a x y F x ( ) 2

dy dx

có dạng :

2 d y 2 dx u x y x ( )

( )

v x y x ( ) ( )

u x v x là các hàm thay thế hằng số

,C C trong (2.1)

trong đó: ( ), ( ) Các hàm u ,v cần tìm thoả mãn hệ phương trình :

u x

v x '( )y (x)+ '( )y (x)=0

(2.3)

u x

v x '( )y (x)+ '( )y (x)=

≠

' 1

' 2

a x ( ) 0

F(x) a x ( ) 0

⎧ ⎪ ⎨ ⎪ ⎩

Dùng qui tắc Cramer giải hệ phương trình (2.3) đối với

u v ta được:

u x v x :

Các phương trình (2.3) sau khi tích phân sẽ thu được các hàm ( ), ( )

u x ( )

= −

∫

y ( )F( ) ξ ξ 2 a W ( ) ( ) ξ ξ 0

d ξ

(2.4)

v x ( )

d ξ

∫

y ( )F( ) ξ ξ 1 a W ( ) ( ) ξ ξ 0

là định thức Wronskian.

trong đó α là hằng số nào đó và

w x = ( )

y (x) y (x) 1 ' y (x) y (x) 1

2 ' 2

Vậy nghiệm tổng quát của phương trình đã cho là :

- 48-

[

(2.5)

d ξ

y c

C y (x)+ y (x)+ 1 1

2 2

∫

] y (x)y ( )-y (x)y ( ) F( ) ξ ξ 1 2 p W ( ) ( ) ξ ξ

Một trong những phương trình vi phân cấp hai có cách giải đơn giản là

(2.6)

F λ+

0 =

2 d F 2 dx

Phương trình này xuất hiện do việc nghiên cứu nhiều phương trình vi phân đạo hàm riêng trong toạ độ Đế các (Descartesian) đối với các hiện tượng vật lý và kỹ thuật. Phương trình vi phân (2.6) chứa tham số λ, vì thế ta sẽ xét 3 trường hợp của tham số: âm, dương và bằng không .

1 Trường hợp 1:

0ω> ) (

2 λ ω= − Phương trình vi phân có dạng:

(2.7)

2 ω−

0 =

2 d F 2 dx

; ta có phương trình đặc trưng là

là phương trình vi phân cấp hai với hệ số hằng số, vì thế người ta có thể giả thiết 0 = với

2 m ω−

F e= ω

x ω−

nghiệm đặc trưng

. Nếu biết được tập

và tập nghiệm cơ bản là { ,x eω e

}

= −

nó có một nghiệm mũ m =⎧ ⎨ m ⎩

nghiệm cơ bản của các nghiệm, có thể tạo nên một tổ hợp tuyến tính của các nghiệm này và sinh ra một tập hợp vô hạn các nghiệm khác.

x ω

x ω−

F x ( )

C e 1

C e 2

trong đó

,C C là các hằng số tuỳ ý, nó phụ thuộc vào các điều kiện bổ sung là

điều kiện ban đầu hoặc điều kiện biên.

2 Trường hợp 2:

0λ=

Nghiệm phương trình vi phân

= có các dạng sau:

2 d F 2 dx

F x ( )

C C x +

F x ( )

)

−

1 K 1

2 K x ( 2

x 0

3 Trường hợp 3:

2 (

λ ω ω=

Phương trình vi phân

2 ω+

= là phương trình vi phân cấp hai với hệ số

2 d F 2 dx

, ta có phương trình đặc

hằng số, vì thế có thể giả thiết nó có một nghiệm mũ

F e= ω

và tập nghiệm cơ bản là

trưng

= với các nghiệm đặc trưng:

2 m ω+

i = − ω

m i =⎧ ⎨ m ⎩

- 49-

i x ω−

. Nếu biết được tập nghiệm cơ bản có thể tạo nên một tổ hợp tuyến tính

i x eω ,

{ e

}

của các nghiệm này và sinh ra một tập hợp vô hạn các nghiệm khác

i x ω

i x ω−

F x ( )

C e 1

C e 2

trong đó:

,C C là các hằng số tuỳ ý, nó phụ thuộc vào các điều kiện bổ sung là

1 điều kiện ban đầu hoặc điều kiện biên.

- 50-

TÀI LIỆU THAM KHẢO Đỗ Đình Thanh. 2002. Phương Pháp Toán Lý. NXB Giáo Dục.

Đỗ Đình Thanh .1996. Phương Pháp Toán Lý. NXB Đại Học Quốc Gia Hà Nội.

Đỗ Đình Thanh. 2002. Phương Pháp Toán Lý. NXB Giáo Dục.

Đỗ Đình Thanh. 2002. Phương Trình Đạo Hàm Riêng Trong Vật Lý. NXB Đại

Học Quốc Gia TPHCM

Đỗ Văn Thông. 2003. Phương Pháp Nghiên Cứu Khoa Học. ĐHAG.

Hồ Xuân Huy. 2005. Phương Pháp Toán Lý. ĐHAG.

Lê Đình Thịnh- Lê Trọng Vinh. 1994. Bài tập Toán Học cao Cấp. NXB Giáo

Dục.

Lê Công Triêm. 2005. Phân Tích Chương Trình Vật Lý Phổ Thông. ĐHAG.

Nguyễn Mạnh Hùng. 2007. Phương Trình Đạo Hàm Riêng. NXB Đại Học Sư

Phạm.

Nguyễn Ngọc Giao. 2003. Phép Tính Toán Tử. NXB Đại Học Quốc Gia TP Hồ

Chí Minh.

Nguyễn Văn Hạp. 1999. Giáo Trình Phương Trình Vi Phân Và Phương Trình

Đạo Hàm Riêng. NXB Đại Học Huế.

Phan Huy Thiện 2006. Phương Trình Toán Lý. NXB Giáo Dục.

Trần Thể. 2005. Lý Luận Vật Lý Phổ Thông. ĐHAG.

Trần Thể. 2005. Bài Tập Vật Lý Phổ Thông, ĐHAG.

- 51-

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 "

Chủ đề:

Toán ứng dụng cho Vật lý

Báo cáo thực tập địa lý và địa chất

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

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

LỜI CẢM ƠN ********

MỤC LỤC

PHẦN I: MỞ ĐẦU

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.3 Khái niệm toán tử, hàm riêng, trị riêng:

1.3.1.Toán tử:

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

1.4. Phương pháp tách biến

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

Tiểu kết :

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

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

Bước 1: Áp dụng phương pháp tách biến Fourier và phương pháp mở rộng hàm

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

∫

Bước 3: Đưa ra hàm Green

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

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

Tiểu kết:

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

3.1 Thiết lập phương trình truyền nhiệt:

Các dạng khác nhau của phương trình truyền nhiệt

trong hệ tọa độ Đề- các.

Các điều kiện ban đầu và điều kiện biên cho phương trình truyền nhiệt

3.2 Bài toán biên phụ thuộc thời gian

3.2.1 Phương pháp tách biến Fourier cho bài toán truyền nhiệt

3.2.2 Phương pháp hàm Green cho bài toán truyền nhiệt

Giải:

∫

3.2.3 Bài toán truyền nhiệt trong miền tròn

Giải:

Cách 1: Phương pháp tách biến Fourier:

Cách 2: Dùng phương pháp hàm Green:

3.3 Bài toán biên truyền nhiệt dừng

Bài toán: Tìm sự phân bố nhiệt độ ở trạng thái dừng trong thanh hữu hạn trên đoạn [0,L].Thanh có nguồn nhiệt Q(x,t) = Q(x). Các điều kiện biên có dạng: u(0) = 0,u(L) = 0

Giải:

Bài toán: Xét phân bố dừng của nhiệt độ trong quả cầu đồng nhất bán kính q với

Giải:

Tiểu kết:

PHỤ LỤC 1 CÁC KHÁI NIỆM CƠ BẢN VỀ PHƯƠNG TRÌNH

ĐẠO HÀM RIÊNG :

PHỤ LỤC 2 PHƯƠNG TRÌNH VI PHÂN CẤP HAI:

1 Trường hợp 1:

}

2 Trường hợp 2:

3 Trường hợp 3:

{ e

}

TÀI LIỆU THAM KHẢO Đỗ Đình Thanh. 2002. Phương Pháp Toán Lý. NXB Giáo Dục.

Tài liệu liên quan

Đồ án tốt nghiệp Kỹ thuật hóa lý-phân tích: Nghiên cứu khả năng xử lý các hợp chất trinitrophenol và acid orange 7 bằng quá trình UV/So dium percarbonate trong môi trường cận trung tính

Khóa luận tốt nghiệp đại học: Ứng dụng phần mềm Geogebra trong dạy học toán ở trường phổ thông

Khóa luận tốt nghiệp: Đặc điểm khoáng hóa sắt khu vực Đông Giang, Hàm Thuận Bắc, Bình Thuận

Khóa luận tốt nghiệp: Nghiên cứu đề xuất khả năng áp dụng mô hình kinh tế tuần hoàn cho ngành khai thác mỏ ở Việt Nam

Khoá luận tốt nghiệp: Phân tích và đánh giá hàm lượng sắt trong nước giếng một số huyện ở tỉnh Bình Dương

Khoá luận tốt nghiệp: Phân tích, đánh giá hàm lượng sắt trong nước giếng tại thị xã Bến Cát, huyện Dầu Tiếng, huyện Phú Giáo, huyện Tân Uyên tỉnh Bình Dương

Đồ án tốt nghiệp Địa chất dầu khí: Phân tích đặc điểm môi trường trầm tích Miocen giữa khu vực Đông Bắc đứt gãy Sông Lô

Khoá luận tốt nghiệp Sư phạm Địa lý: Giáo dục kĩ năng sống cho học sinh thông qua bài học Địa lí Sông và hồ THCS lớp 6

Khóa luận tốt nghiệp: Ứng dụng công nghệ tin học và máy toàn đạc điện tử thực hiện thành lập bản đồ địa chính tờ số 17 tỷ lệ 1:1000 tại xã Đông Quang, huyện Ba Vì, thành phố Hà Nội

Khóa luận tốt nghiệp: Ứng dụng GIS trong công tác quản lý các khu vực sạt lở đất thuộc xã Huổi Lèng, huyện Mường Chà, tỉnh Điện Biên

Tài liêu mới

Bài giảng Cơ học vật rắn và sóng cơ: Chương 9 - Chuyển động phức hợp của điểm

Bài giảng Cơ học vật rắn và sóng cơ: Chương 8 - Hai chuyển động cơ bản của vật rắn

Bài giảng Cơ học vật rắn và sóng cơ: Chương 7 - Động học điểm (Kinematics of particles)

Thủy tinh lithium borate pha tạp ion Eu3+ (Li2B4O7:Eu3+): Cấu trúc, tính chất vật lý và đặc tính đo liều bức xạ

Lesson The basics of radiation and radiation protection

Bài tập Hóa đại cương: Nhiệt động hóa học (Có đáp án)

Bài giảng Vật lý bán dẫn: Chương 4 - Hồ Trung Mỹ (Trường ĐH BK TP.HCM)

Bài giảng Vật lý bán dẫn: Chương 3 - Hồ Trung Mỹ (Trường ĐH BK TP.HCM)

Bài giảng Vật lý bán dẫn: Chương 2 - Hồ Trung Mỹ (Trường ĐH BK TP.HCM)

Bài giảng Vật lý bán dẫn: Chương 1 - Hồ Trung Mỹ (Trường ĐH BK TP.HCM)