Phương trình Elliptic: Lý thuyết, nghiệm và các điều kiện biên Dirichlet, Neumann

Tài liệu giảng dạy Phương trình elliptic: laplace và poisson (p2) hệ thống hóa kiến thức từ cơ bản đến nâng cao ngành tại Việt Nam

• Trường học: Trường Đại Học Khoa Học Tự Nhiên
• Chuyên ngành: Toán Lý
• Đề tài: Giáo Trình Phương Trình Toán Lý Phần 2
• Loại tài liệu: Giáo Trình
• Năm xuất bản: 2023
• Địa điểm: Hà Nội

Mục lục chi tiết

• 4. CHƯƠNG IV: PHƯƠNG TRÌNH ELLIPTIC
• 4.1. Phương trình Laplace và các điều kiện biên
• 4.1.1. Nghiệm của phương trình Laplace trong các hệ tọa độ
• 4.1.1.1. Nghiệm phương trình Laplace trong tọa độ Cartesian
• 4.1.1.2. Nghiệm phương trình Laplace trong tọa độ trụ
• 4.1.1.3. Nghiệm phương trình Laplace trong tọa độ cầu
• 4.2. LÝ THUYẾT THẾ
• 4.2.1. Trường hấp dẫn
• 4.2.2. Sự di chuyển của chất lỏng
• 4.2.3. Trường tĩnh điện
• 4.2.4. Trường tĩnh từ
• 4.3. PHƯƠNG TRÌNH HELMHOLTZ
• 4.3.1. Phương trình Helmholtz trong hệ tọa độ Cartesian (x, y, z)
• 4.3.2. Phương trình Helmholtz trong hệ tọa độ trụ (r, z)
• 4.3.3. Phương trình Helmholtz trong hệ tọa độ cầu (r, θ, φ)
• 4.4. HÀM ĐIỀU HOÀ VÀ CÁC TÍNH CHẤT
• 4.5. PHƯƠNG TRÌNH POISSON TRONG MIỀN CHỮ NHẬT
• 4.6. CÔNG THỨC TÍCH PHÂN POISSON TRONG MIỀN TRÒN

Tóm tắt

I. Hướng dẫn toàn tập về phương trình Elliptic trong toán lý 2

Giáo trình phương trình toán lý phần 2 tập trung sâu vào các phương trình đạo hàm riêng loại Elliptic. Đây là nền tảng cho việc mô tả các hiện tượng vật lý ở trạng thái dừng, tức là không thay đổi theo thời gian. Trọng tâm của chương này là ba phương trình cốt lõi: phương trình Laplace, phương trình Poisson, và phương trình Helmholtz. Mỗi phương trình này đều có dạng tổng quát liên quan đến toán tử Laplace (∇²), nhưng khác nhau về vế phải, tương ứng với trường hợp không có nguồn, có nguồn phân bố, hoặc liên quan đến bài toán trị riêng. Sự hiểu biết về phương trình Elliptic là cực kỳ quan trọng trong nhiều lĩnh vực khoa học và công nghệ, từ tĩnh điện, từ học, truyền nhiệt, đến cơ học chất lỏng. Một điểm đặc biệt của các phương trình này là dạng biểu diễn của chúng thay đổi tùy thuộc vào hệ tọa độ được sử dụng. Việc lựa chọn hệ tọa độ phù hợp (Cartesian, trụ, cầu) không chỉ giúp đơn giản hóa bài toán mà còn là chìa khóa để tìm ra nghiệm chính xác. Nội dung của giáo trình sẽ đi sâu vào việc xây dựng và giải các phương trình này trong từng hệ tọa độ, kèm theo các điều kiện biên phổ biến như Dirichlet, Neumann và Robin. Những điều kiện này mô tả trạng thái của hệ tại biên của miền khảo sát, là yếu tố không thể thiếu để đảm bảo nghiệm của bài toán là duy nhất và có ý nghĩa vật lý.

1.1. Khái niệm phương trình Laplace và các điều kiện biên

Phương trình Laplace là trường hợp cơ bản nhất của phương trình Elliptic, có dạng ∇²u = 0. Phương trình này mô tả các trường thế trong một miền không có nguồn, chẳng hạn như điện thế trong vùng không có điện tích hoặc nhiệt độ ở trạng thái dừng trong một vật không có nguồn nhiệt. Nghiệm của phương trình Laplace được gọi là các hàm điều hòa. Để tìm được một nghiệm duy nhất, phương trình phải đi kèm với các điều kiện biên. Có ba loại điều kiện biên chính:

• Điều kiện biên Dirichlet: Giá trị của hàm u được xác định trực tiếp trên biên của miền khảo sát (u = g(x,y,z)).
• Điều kiện biên Neumann: Đạo hàm theo phương pháp tuyến của hàm u được xác định trên biên (∂u/∂n = g(x,y,z)). Điều này tương ứng với việc xác định dòng chảy qua biên.
• Điều kiện biên Robin: Là sự kết hợp tuyến tính của hai điều kiện trên (αu + β(∂u/∂n) = g(x,y,z)). Việc lựa chọn và áp dụng đúng điều kiện biên là bước thiết yếu để mô hình hóa chính xác bài toán vật lý.

1.2. Phân biệt phương trình Poisson và phương trình Helmholtz

Trong khi phương trình Laplace mô tả trạng thái không nguồn, phương trình Poisson (∇²u = f(x,y,z)) lại mô tả các hệ có nguồn phân bố trong miền. Hàm f(x,y,z) biểu diễn mật độ của nguồn, ví dụ như mật độ điện tích trong điện trường hoặc nguồn nhiệt trong bài toán truyền nhiệt. Phương trình này là sự mở rộng trực tiếp của phương trình Laplace và đóng vai trò trung tâm trong lý thuyết thế. Mặt khác, phương trình Helmholtz (∇²u + λu = 0) thường xuất hiện sau khi áp dụng phương pháp tách biến cho các phương trình phụ thuộc thời gian như phương trình truyền sóng hoặc phương trình khuếch tán. Trong phương trình này, λ là hằng số tách biến, thường liên quan đến các trị riêng của hệ. Do đó, việc giải phương trình Helmholtz thực chất là một bài toán trị riêng, tìm các giá trị λ để tồn tại nghiệm không tầm thường. Cả hai phương trình này đều là những công cụ mạnh mẽ được đề cập trong giáo trình phương trình toán lý phần 2.

1.3. Vai trò toán tử Laplace trong các hệ tọa độ phổ biến

Toán tử Laplace (∇²) có biểu thức khác nhau trong các hệ tọa độ khác nhau, phản ánh cấu trúc hình học của không gian. Trong hệ tọa độ Cartesian (x, y, z), nó có dạng đơn giản nhất là tổng các đạo hàm riêng cấp hai theo mỗi biến. Tuy nhiên, đối với các bài toán có tính đối xứng trụ hoặc đối xứng cầu, việc sử dụng tọa độ Cartesian trở nên phức tạp. Thay vào đó, hệ tọa độ trụ (r, φ, z) và hệ tọa độ cầu (r, θ, φ) cung cấp một mô tả tự nhiên và hiệu quả hơn. Trong tọa độ trụ, toán tử Laplace chứa các đạo hàm theo bán kính r, góc φ và chiều cao z. Trong tọa độ cầu, nó lại được biểu diễn qua bán kính r và hai góc θ, φ. Việc nắm vững biểu thức của toán tử Laplace trong các hệ tọa độ này là điều kiện tiên quyết để có thể giải thành công các phương trình Elliptic cho các miền hình học phức tạp.

II. Bí quyết giải phương trình Laplace trong các hệ tọa độ khác nhau

Giải phương trình Laplace là một trong những kỹ năng cốt lõi được trình bày trong giáo trình phương trình toán lý phần 2. Phương pháp chung và hiệu quả nhất là phương pháp tách biến. Ý tưởng của phương pháp này là giả định nghiệm của phương trình có thể được biểu diễn dưới dạng tích của các hàm một biến, ví dụ u(x,y,z) = X(x)Y(y)Z(z). Khi thay dạng nghiệm này vào phương trình Laplace, phương trình đạo hàm riêng ban đầu sẽ được chuyển thành một hệ các phương trình vi phân thường, vốn dễ giải hơn rất nhiều. Tuy nhiên, dạng của các phương trình vi phân thường này và nghiệm của chúng phụ thuộc hoàn toàn vào hệ tọa độ được chọn. Mỗi hệ tọa độ sẽ dẫn đến một họ các hàm đặc biệt khác nhau. Ví dụ, tọa độ Cartesian dẫn đến các hàm mũ và lượng giác; tọa độ trụ dẫn đến các hàm Bessel; và tọa độ cầu dẫn đến các đa thức Legendre. Việc lựa chọn đúng hệ tọa độ không chỉ dựa vào hình dạng của miền biên mà còn là một chiến lược để đơn giản hóa quá trình tìm nghiệm. Hiểu rõ mối liên hệ giữa hệ tọa độ và các hàm đặc biệt là bí quyết để chinh phục các bài toán phức tạp liên quan đến phương trình Laplace và các nghiệm hàm điều hòa của nó.

2.1. Giải pháp trong hệ tọa độ Cartesian tọa độ Đề các

Trong hệ tọa độ Cartesian, phương trình Laplace có dạng ∂²u/∂x² + ∂²u/∂y² + ∂²u/∂z² = 0. Khi áp dụng phương pháp tách biến bằng cách đặt u(x,y,z) = X(x)Y(y)Z(z), ta thu được ba phương trình vi phân thường: X''/X = α², Y''/Y = β², Z''/Z = γ², với điều kiện α² + β² + γ² = 0. Nghiệm của các phương trình này thường là tổ hợp của các hàm sin, cosin, sinh, và cosh, hoặc hàm mũ. Dạng nghiệm cụ thể phụ thuộc vào các hằng số tách biến (âm hoặc dương), vốn được quyết định bởi điều kiện biên của bài toán. Hệ tọa độ này đặc biệt hiệu quả cho các bài toán trong miền hình hộp chữ nhật, nơi các biên song song với các trục tọa độ. Nghiệm tổng quát sau đó được xây dựng bằng cách chồng chất các nghiệm riêng, thường dưới dạng một chuỗi Fourier.

2.2. Tìm nghiệm trong hệ tọa độ trụ và sự xuất hiện hàm Bessel

Đối với các bài toán có đối xứng trụ, chẳng hạn như tìm điện thế trong một dây dẫn hình trụ, việc sử dụng hệ tọa độ trụ (r, φ, z) là lựa chọn tối ưu. Phương trình Laplace trong tọa độ trụ trở nên phức tạp hơn, nhưng phương pháp tách biến vẫn có thể áp dụng. Sau khi tách biến, phương trình cho biến góc φ và biến z vẫn dẫn đến các hàm lượng giác và hàm mũ quen thuộc. Tuy nhiên, phương trình cho biến bán kính r lại có một dạng đặc biệt: r²R'' + rR' + (k²r² - n²)R = 0. Đây chính là phương trình Bessel. Nghiệm của nó là các hàm Bessel loại một (Jₙ) và loại hai (Yₙ). Các hàm này mô tả sự biến thiên của trường theo phương bán kính trong các hệ có đối xứng trụ. Việc hiểu và sử dụng các hàm Bessel là không thể thiếu khi làm việc với các bài toán vật lý trong hình học trụ.

2.3. Lời giải trong hệ tọa độ cầu và đa thức Legendre liên kết

Khi xử lý các bài toán có đối xứng cầu, như trường hấp dẫn của một hành tinh, hệ tọa độ cầu (r, θ, φ) là công cụ mạnh nhất. Áp dụng phương pháp tách biến cho phương trình Laplace trong tọa độ cầu, ta thu được các phương trình vi phân thường cho từng biến. Phương trình cho bán kính r là một phương trình dạng Cauchy-Euler. Phương trình cho góc phương vị φ cho nghiệm lượng giác. Đáng chú ý nhất là phương trình cho góc thiên đỉnh θ, sau một phép đổi biến x = cos(θ), sẽ dẫn đến phương trình Legendre liên kết. Nghiệm của phương trình này là các đa thức Legendre liên kết Pₙᵐ(x). Các hàm này, kết hợp với các hàm lượng giác theo góc φ, tạo thành một họ các hàm trực giao trên mặt cầu được gọi là các hàm cầu điều hòa. Chúng đóng vai trò cực kỳ quan trọng trong nhiều ngành như cơ học lượng tử, địa vật lý và điện động lực học.

III. Phân tích sâu phương trình Helmholtz Nền tảng cốt lõi toán lý

Phương trình Helmholtz, với dạng ∇²H + λH = 0, là một trong những phương trình quan trọng nhất trong giáo trình phương trình toán lý phần 2. Nó không chỉ là một phương trình Elliptic độc lập mà còn là cầu nối giữa các bài toán trạng thái dừng và các bài toán phụ thuộc thời gian. Phương trình này thường nảy sinh từ việc áp dụng phương pháp tách biến cho phương trình sóng hoặc phương trình truyền nhiệt/khuếch tán. Trong quá trình này, phần không gian của nghiệm thỏa mãn phương trình Helmholtz, trong khi phần thời gian tuân theo một phương trình vi phân thường đơn giản. Hằng số tách biến λ, hay trị riêng, có ý nghĩa vật lý quan trọng, thường đại diện cho tần số dao động hoặc tốc độ tắt dần của hệ. Do đó, việc giải phương trình Helmholtz tương đương với việc tìm các trạng thái riêng (eigenstates) và tần số riêng (eigenfrequencies) của một hệ vật lý. Giống như phương trình Laplace, việc giải phương trình Helmholtz cũng phụ thuộc mạnh mẽ vào hệ tọa độ được chọn, dẫn đến sự xuất hiện của các họ hàm đặc biệt khác nhau như hàm Bessel, hàm Bessel cầu, và các đa thức Legendre.

3.1. Nguồn gốc từ phương trình truyền nhiệt và phương trình sóng

Phương trình Helmholtz không phải là một phương trình cơ bản mà thường là kết quả của một bước phân tích. Ví dụ, xét phương trình truyền nhiệt ∇²u = (1/α²)∂u/∂t. Khi tìm nghiệm ở dạng tách biến u(x,t) = H(x)T(t), ta có thể biến đổi phương trình thành (∇²H)/H = (1/α²)(T'/T). Vì vế trái chỉ phụ thuộc vào không gian (x) và vế phải chỉ phụ thuộc vào thời gian (t), cả hai vế phải bằng cùng một hằng số, ký hiệu là -λ. Điều này dẫn đến hai phương trình: phương trình không gian ∇²H + λH = 0 (chính là phương trình Helmholtz) và phương trình thời gian T' + α²λT = 0. Một quy trình tương tự cũng áp dụng cho phương trình sóng, cho thấy vai trò phổ quát của phương trình Helmholtz trong việc phân tích các chế độ dao động và khuếch tán.

3.2. Phương pháp tách biến cho phương trình Helmholtz trong tọa độ trụ

Giải phương trình Helmholtz trong hệ tọa độ trụ là một bài toán điển hình trong nhiều ứng dụng, như sóng điện từ trong ống dẫn sóng hình trụ. Sử dụng nghiệm tách biến H(r,φ,z) = R(r)Φ(φ)Z(z), phương trình đạo hàm riêng ban đầu được chia thành ba phương trình vi phân thường. Phương trình cho Φ(φ) và Z(z) dẫn đến các nghiệm lượng giác và mũ quen thuộc. Phương trình cho biến bán kính R(r) một lần nữa lại là một phương trình Bessel hoặc phương trình Bessel biến thể, tùy thuộc vào dấu của các hằng số tách biến. Cụ thể, nếu λ - k_z² > 0, ta có phương trình Bessel chuẩn. Nếu λ - k_z² < 0, ta thu được phương trình Bessel biến thể, có nghiệm là các hàm Bessel biến thể Iₙ và Kₙ, mô tả các sóng tắt dần theo phương bán kính.

3.3. Nghiệm trong tọa độ cầu và hàm Bessel cầu

Trong hệ tọa độ cầu, việc giải phương trình Helmholtz có vai trò quan trọng trong các bài toán tán xạ sóng hoặc cơ học lượng tử. Sau khi tách biến, phần góc của nghiệm vẫn được mô tả bởi các hàm cầu điều hòa, tương tự như trường hợp phương trình Laplace. Tuy nhiên, phương trình cho phần bán kính R(r) sẽ có dạng: r²R'' + 2rR' + [k²r² - n(n+1)]R = 0. Đây là phương trình Bessel cầu. Nghiệm của nó là các hàm Bessel cầu (jₙ, yₙ) và các hàm Hankel cầu. Các hàm này mô tả các sóng cầu lan truyền ra xa hoặc hội tụ về tâm. Sự xuất hiện của các hàm đặc biệt này một lần nữa khẳng định tầm quan trọng của việc nắm vững các công cụ toán học cao cấp khi nghiên cứu các phương trình toán lý.

IV. Hướng dẫn giải phương trình Poisson trong các miền hình học

Phương trình Poisson, ∇²u = f(x,y), là một sự mở rộng của phương trình Laplace để mô tả các hệ thống có sự hiện diện của nguồn hoặc mật độ điện tích. Việc giải phương trình này là một chủ đề quan trọng trong giáo trình phương trình toán lý phần 2, đòi hỏi các kỹ thuật tinh vi hơn so với trường hợp thuần nhất (Laplace). Một trong những phương pháp tiếp cận hiệu quả là dựa trên nghiệm của bài toán trị riêng tương ứng. Cụ thể, ta tìm các hàm riêng (eigenfunctions) Φₙ của toán tử Laplace (∇²Φₙ + λₙΦₙ = 0) với cùng điều kiện biên. Do tính trực giao và đầy đủ của hệ hàm riêng này, cả hàm nguồn f(x,y) và nghiệm u(x,y) đều có thể được khai triển thành chuỗi theo các hàm riêng đó. Bằng cách thay các khai triển này vào phương trình Poisson, ta có thể xác định được các hệ số trong chuỗi khai triển của nghiệm u. Một phương pháp mạnh mẽ khác là sử dụng công thức tích phân, đặc biệt là công thức tích phân Poisson cho miền tròn, cung cấp một biểu thức tường minh cho nghiệm dựa trên giá trị biên.

4.1. Kỹ thuật giải phương trình Poisson trong miền chữ nhật

Đối với phương trình Poisson trong một miền chữ nhật với điều kiện biên Dirichlet thuần nhất (u=0 trên biên), phương pháp hàm riêng là cực kỳ hiệu quả. Trước hết, ta giải bài toán trị riêng ∇²Φ + λΦ = 0 trong cùng miền chữ nhật. Bằng phương pháp tách biến, các hàm riêng được tìm thấy có dạng Φₙₘ(x,y) = sin(nπx/a)sin(mπy/b), tương ứng với các trị riêng λₙₘ. Đây chính là các hàm cơ sở của chuỗi Fourier hai chiều. Sau đó, ta khai triển cả hàm nguồn f(x,y) và nghiệm u(x,y) theo hệ hàm riêng này. Thay các chuỗi vào phương trình Poisson và sử dụng tính trực giao của các hàm sin, ta có thể tìm được mối liên hệ đại số đơn giản giữa các hệ số khai triển của u và f, từ đó xác định được nghiệm cuối cùng.

4.2. Công thức tích phân Poisson cho nghiệm trong miền tròn

Đối với bài toán Dirichlet cho phương trình Laplace (một trường hợp đặc biệt của phương trình Poisson với f=0) trong một miền tròn, có một công cụ rất mạnh mẽ là công thức tích phân Poisson. Công thức này cho phép biểu diễn giá trị của nghiệm u(r,φ) tại một điểm bất kỳ bên trong hình tròn thông qua một tích phân của các giá trị biên f(φ) trên đường tròn. Biểu thức này có dạng: u(r,φ) = (1/2π) ∫ P(r, a, φ-ψ)f(ψ)dψ, trong đó a là bán kính hình tròn và P là một hàm đặc biệt được gọi là nhân Poisson. Công thức này không chỉ cung cấp một lời giải tường minh mà còn cho thấy giá trị của hàm điều hòa tại một điểm là một trung bình có trọng số của các giá trị trên biên, với trọng số được quyết định bởi nhân Poisson.

4.3. Ý nghĩa của nhân Poisson và các bài toán biên liên quan

Nhân Poisson, P(r, a, γ) = (a² - r²)/(a² - 2arcos(γ) + r²), đóng vai trò như một hàm trọng số trong công thức tích phân. Nó luôn dương khi r < a và đạt giá trị cực đại khi γ=0, tức là tại điểm trên biên gần nhất với điểm đang xét. Điều này phản ánh một tính chất vật lý trực quan: giá trị tại một điểm bị ảnh hưởng mạnh nhất bởi các giá trị biên ở gần nó. Công thức tích phân Poisson là nền tảng để chứng minh nhiều tính chất quan trọng của hàm điều hòa, chẳng hạn như tính trơn (khả vi vô hạn) bên trong miền và nguyên lý giá trị trung bình. Mặc dù được trình bày cho phương trình Laplace, kỹ thuật này có thể được mở rộng để giải phương trình Poisson tổng quát thông qua việc sử dụng hàm Green.

V. Ứng dụng thực tiễn của phương trình toán lý trong vật lý

Các phương trình Elliptic không chỉ là những cấu trúc toán học trừu tượng mà còn là công cụ không thể thiếu để mô hình hóa vô số hiện tượng vật lý ở trạng thái cân bằng. Giáo trình phương trình toán lý phần 2 dành một phần quan trọng để minh họa những ứng dụng này, đặc biệt thông qua lý thuyết thế. Lý thuyết này nghiên cứu các trường vật lý có thể được biểu diễn bằng gradien của một hàm vô hướng, gọi là thế. Nếu một trường vector là không xoáy (rot V = 0), nó có thể được sinh ra từ một trường vô hướng (thế) theo hệ thức V = -∇ψ. Nếu trường vector này đồng thời là trường cảm ứng (div V = 0), thì thế ψ sẽ thỏa mãn phương trình Laplace ∇²ψ = 0. Lý thuyết này có ứng dụng trực tiếp trong việc mô tả trường hấp dẫn, trường tĩnh điện, trường tĩnh từ, và cả chuyển động của chất lỏng không nhớt. Việc hiểu rõ mối liên hệ giữa các khái niệm vật lý và các phương trình toán học này giúp làm sáng tỏ bản chất của nhiều định luật tự nhiên.

5.1. Lý thuyết thế Trường vô hướng và trường vector

Một trường vector V được gọi là trường thế nếu nó là gradien của một trường vô hướng ψ, tức V = ±∇ψ. Hàm ψ được gọi là thế vô hướng. Điều kiện cần để một trường vector là trường thế là rot V = 0, nghĩa là trường đó không có xoáy. Ngược lại, một trường vector được gọi là trường cảm ứng nếu div V = 0. Khi một trường vector vừa không xoáy vừa cảm ứng, nó sẽ được sinh ra từ một thế vô hướng ψ là nghiệm của phương trình Laplace (∇²ψ = 0). Điều này tạo ra một sự kết nối sâu sắc giữa các tính chất hình học của trường vector (không xoáy, không nguồn) và tính chất giải tích của hàm thế (là một hàm điều hòa). Các công cụ như toán tử gradien, div, và rot trong các hệ tọa độ khác nhau là chìa khóa để áp dụng lý thuyết này.

5.2. Mô hình hóa trường hấp dẫn và trường tĩnh điện

Trong vật lý cổ điển, trường hấp dẫn g do một vật có khối lượng M sinh ra là một trường vector. Trường này có thể được biểu diễn qua thế hấp dẫn ψ = -GM/r, theo công thức g = -∇ψ. Dễ dàng kiểm chứng rằng trong không gian không có khối lượng (r > 0), thế hấp dẫn này là một nghiệm của phương trình Laplace, ∇²ψ = 0. Tương tự, trong trường tĩnh điện, định luật Faraday (rot E = 0) cho thấy điện trường E có thể được sinh ra từ điện thế vô hướng V (E = -∇V). Định luật Gauss (div E = ρ/ε₀) chỉ ra rằng điện thế V thỏa mãn phương trình Poisson ∇²V = -ρ/ε₀. Trong vùng không gian không có điện tích (ρ=0), phương trình này trở thành phương trình Laplace, cho thấy điện thế là một hàm điều hòa.

5.3. Áp dụng trong cơ học chất lỏng và trường tĩnh từ

Lý thuyết thế cũng có ứng dụng quan trọng trong các lĩnh vực khác. Trong cơ học chất lỏng, đối với một dòng chảy của chất lỏng lý tưởng (không nén được và không nhớt), nếu dòng chảy là không xoáy (rot v = 0), thì vector vận tốc v có thể được biểu diễn qua thế vận tốc φ (v = -∇φ). Điều kiện không nén được (div v = 0) dẫn đến kết luận rằng thế vận tốc φ phải thỏa mãn phương trình Laplace. Trong trường tĩnh từ, khi không có dòng điện (J=0), định luật Ampere trở thành rot B = 0. Điều này cho phép định nghĩa một thế từ vô hướng ψₘ sao cho B = -∇ψₘ. Kết hợp với định luật Gauss cho từ trường (div B = 0), ta lại thu được phương trình Laplace ∇²ψₘ = 0. Những ví dụ này cho thấy sức mạnh thống nhất của các phương trình toán lý trong việc mô tả thế giới tự nhiên.

VI. Kết luận giáo trình toán lý phần 2 Tính chất hàm điều hòa

Phần cuối của chương về phương trình Elliptic trong giáo trình phương trình toán lý phần 2 tập trung tổng kết các tính chất quan trọng của hàm điều hòa, tức là các nghiệm của phương trình Laplace. Các tính chất này không chỉ có ý nghĩa toán học sâu sắc mà còn phản ánh những nguyên lý vật lý cơ bản. Một hàm điều hòa trong một miền V bị chặn bởi bề mặt S có những đặc tính rất đặc biệt. Chẳng hạn, tích phân của đạo hàm theo pháp tuyến của nó trên toàn bộ bề mặt S luôn bằng không, thể hiện sự bảo toàn dòng. Giá trị của hàm tại một điểm bất kỳ bên trong một hình cầu bằng giá trị trung bình của hàm trên bề mặt của hình cầu đó (nguyên lý giá trị trung bình). Tính chất này dẫn đến một hệ quả cực kỳ quan trọng là nguyên lý cực đại và cực tiểu. Một hàm điều hòa không thể đạt giá trị lớn nhất hoặc nhỏ nhất tại một điểm bên trong miền, mà các giá trị cực trị này phải nằm trên biên. Đây là cơ sở để chứng minh tính duy nhất của nghiệm cho bài toán Dirichlet và có ý nghĩa vật lý rõ ràng: nhiệt độ cao nhất trong một vật không có nguồn nhiệt phải nằm ở biên của nó.

6.1. Định nghĩa và các tính chất cốt lõi của hàm điều hòa

Một hàm u(x,y,z) được gọi là hàm điều hòa trong một miền nếu nó có đạo hàm riêng cấp hai liên tục và thỏa mãn phương trình Laplace (∇²u = 0) tại mọi điểm trong miền đó. Các hàm này còn được gọi là hàm thế. Một trong những tính chất cơ bản là nếu u là hàm điều hòa, thì tích phân ∫(∂u/∂n)dS trên một mặt kín S bao quanh miền luôn bằng không. Một tính chất quan trọng khác là nguyên lý giá trị trung bình: giá trị của hàm điều hòa tại tâm một hình cầu bằng trung bình cộng các giá trị của nó trên bề mặt hình cầu đó. Tính chất này cho thấy sự 'mượt mà' và 'san bằng' vốn có của các hiện tượng vật lý ở trạng thái dừng mà hàm điều hòa mô tả.

6.2. Nguyên lý cực đại và tính duy nhất của nghiệm bài toán Dirichlet

Một hệ quả trực tiếp của nguyên lý giá trị trung bình là nguyên lý cực đại/cực tiểu. Nguyên lý này phát biểu rằng một hàm điều hòa không hằng số trong một miền bị chặn không thể đạt giá trị cực đại hoặc cực tiểu tại một điểm trong của miền. Các điểm cực trị này phải nằm trên biên. Nguyên lý này có ý nghĩa vật lý sâu sắc, ví dụ, điểm nóng nhất hoặc lạnh nhất của một vật thể ở trạng thái cân bằng nhiệt phải nằm trên bề mặt của nó. Về mặt toán học, nguyên lý cực đại là công cụ chính để chứng minh tính duy nhất của nghiệm cho bài toán Dirichlet. Nếu có hai hàm điều hòa u₁ và u₂ cùng thỏa mãn một điều kiện biên Dirichlet cho trước, thì hiệu của chúng (u₁ - u₂) cũng là một hàm điều hòa và bằng không trên biên. Theo nguyên lý cực đại, hiệu này phải bằng không ở mọi nơi, suy ra u₁ = u₂.

6.3. Tầm quan trọng của hàm điều hòa trong nghiên cứu toán lý

Hàm điều hòa không chỉ là nghiệm của một phương trình cụ thể mà còn là một lớp hàm có cấu trúc đẹp và phong phú, là đối tượng nghiên cứu trung tâm của lý thuyết thế. Chúng xuất hiện trong hầu hết các lĩnh vực của vật lý lý thuyết và kỹ thuật. Hơn nữa, tính chất của hàm điều hòa vẫn được bảo toàn qua nhiều phép biến đổi, chẳng hạn như phép biến đổi bảo giác trong mặt phẳng phức. Nghiên cứu các phương trình Elliptic và các hàm điều hòa trong giáo trình phương trình toán lý phần 2 cung cấp cho người học một bộ công cụ giải tích mạnh mẽ và một sự thấu hiểu sâu sắc về các nguyên lý cân bằng và trạng thái dừng trong tự nhiên, mở đường cho các nghiên cứu nâng cao hơn trong toán học và vật lý ứng dụng.

Trích đoạn nội dung tài liệu

Chuong IV PHƯƠNG TRÌNH ELLIPTIC §1. Phương trình Laplace và các điều kiện biên Phương trình Laplace thuần nhất có dạng: V u=U: x.1) và dạng không thuần nhất là: V u= /(x.2) Phương trình Poisson hay là phương trình Hemholiz có dạng: V w+ Âu =Ú: x.3) Các phương trình trên được gọi là phương trình Elliptic luôn luôn tìm thấy trong các ngành học khoa học và công nghệ liên quan đến các bài toán dừng. Phương trình Hemholtz liên quan đến bài toán trị riêng Elliptic, trong đó trị riêng 4 được tìm sao cho nghiệm là khác không. Các phương trình trên có dạng khác nhau phụ thuộc vào dạng của toán tử Laplace trong các hệ tọa độ khác nhau.

Trong chương này, chúng ta sẽ nghiên cứu dạng khác nhau của các phương trình trên một cách đầy đủ và minh họa các ứng dụng cua chúng trong một vài ngành khoa học và công nghệ. Ta có bảng sau mô tả dạng của phương trình Laplace trong các hệ tọa độ thông thường: Tọa độ 1 chiều | Cartesian V?(x) 1% L dx’ | Trụ ID — | yrˆ or dr › 2 1 d(j,dat) | Cau vu()= 4S 189 Tỏa độ | 2 chiều 2 2 Cartesian Vˆu(x,y)= Ty Ou lou Ou V2u(r,z)= ey uly a ê? rôr a2? Ou lôu 1 Ou Tru V'u(r,.0)= sat 5 ta agt Ou ld 1 6w / 2 Z W?u(r,0)=—+—— +-r— u(r ø) 2 r or rì ag’ V 2 Vˆu(0,p - GF ( sino) t ou Ci sinØ ôØ 08) sin’ @ dg mu ld du 1 af Ou Vru(ro)= ear) — (sino) r° dr dr) r°sin@ 00 989 } r sin? 6 Ag? Chúng ta sẽ giải chỉ tiết các phương trình Laplace, Poisson và Hemholtz với các điêu kiện biên sau: , — Diéu kiện bién Dirichlet: e ); Vx, y.z B(u)=u= (x,y,z 6 OR, — Diéu kién bién Neumann: A : ⁄ ôu - “3 B(u)= a Vu-ñi = g(x, y,z};Vx, y,z e ÔR, H là ~ Điều kiện biên Robin:.y,2);Vx, y,z € ðR, moe 190 trong dó hằng số w.z) được xác dịnh: diều kiện x.z nằm trên biên của miễn 8. được viết là x. Như vậy bài toán biên cho phương trinh Elliptic duge viét: Giải phương trình feg2 Vou LẦu =! .z € OR: cn ~ Khi a =0.f =1 la diéu kién bién Dirichlet: ~ Khi ø =I./= 0 là điều kiện biên Neumann: ~ Khi ø #0.

# 0 là điều kiện biên Robin. E7 dụ: Nhiệt độ của một vật đồng chất có nguồn nhiệt trong Ñ thoả mãn phương trình truyền nhiệt ôi 16? a =d VÌu+q(x. / Trạng thái dừng của phương trình trên cho ta phương trình Poisson q(x. a Trang thái dừng của phương trình trên khi không có nguồn cho ta phương trình Laplace VẦu=0, u=u(x,y,z).

Như vậy, phương trinh Poisson và phương trinh Laplace ở trên là trường hợp đặc biệt của phương trình truyền nhiệt khi hàm w không phụ thuộc vào biêm thời gian, được gọi là trạng thái dừng của phương trình truyền nhiệt. Nghiệm của phương trình Laplace trong các hệ tọa độ Nghiệm của phương trình Laplace còn được gọi là các hàm điều hoà. a) Nghiệm phương trình Laplace trong tọa độ Cartesian Cho ham ¿ = ¿(x, y,z), phương trình Laplace trong hệ tọa độ Cartesian 3 chiều có dạng: Ou Ou Ou - `. , &}* ` az? , 191 Đặt uŒ.

ta có: i X"+a@?X =0 +X (x)= Ae’ + Be xO OY" 7! Ty ty? ¥"+ BY £0 >Y(y)=Ce”" + De" arp, Z"-y°Z=0—>Z(z) = Ee” + Fe" trong đó: y?=ø?+/”. b) Nghiệm phương trình Laplace trong tọa độ frụ Cho ham u = u(r,@. phương trình Laplace trong hệ tọa độ trụ 3 chiều có dạng: : Au(r,ø,z 12 (rR) rap tan. thay vào phương trình, ta được: Lier} gL 22; rR ro Z O" +n’ =0 ® =s(rRj +r' T=—--=n = Z”—-k*Z=0 (*) r*R"+rR'+(k?r?—n”)R=0 Phương trình (*) là phương trình xác định hàm Bessel.

Trong trường hợp k =0, tức là không có hàm Z, phương trình (*) có dạng: r}R"+rR'—n°R =0. Ứ +B,Inr Nghiệm # là on " gu: ne + ¡P 192 Ham ® có dạng C+ Do Dey iC cosng + 2 sina Vậy nghiệm của phương trình Laplace khí & = 0 lar u{n./@)= E + B Inr+ Sim + Br "\(C,cosne +d, sinne). Công thức trên được gọi là hàm điền hoà trụ. [rong trường hợp & >0.

ta có nghiệm cua phương trình (*) là R (kr)= Ad, (kr)+ BY, (kr). Nghiệm cua phương trình có dang u(rg.z =D [4 If ‘)+B,,, ¥ (k,n) je ev oO) Nghiém phuong trinh Laplace trong toa db can Cho ham u=u(r. phuong trinh Laplace trong hé toa do tru 3 hiếu có dạng: 1 @f sc) 1 ê( Ou) 1 Ou Au(r.~) = = in J# — ( sind a Or er) rr sind eo 00 me * sin? 6 6g ` Dat u(r. thay vào phương trình.

ta được: 1 ụ ' 1 r(rÌR (sin2?') +—,—— ©" = 0: Rr- Pe * sind br? sin’ 0 w® rR sĩ ) † + „ ( | n9 1 LUIỆP § R Psin2 sin Ø @® Đặt rR) ề s (rR) oy a —m R œ (sin OP ), i oe PsinO sin0 ® SUIY ra r`Rh+2rR'- 4Ñ =0. sin? dO do sin? 0 193 Trong phuong trinh (*) dat x =cos@ = dx = -sin@d0 => d@ = - " sind” ta CÓ: _1 d(. Đa thức Legendre liên kết có dạng Pp „ (x)=(I-x 2 )? 2 a"P, A (x ). Do đó yy" m g™ (tx? # Pi" (x)= (=!) 2m! a el am Le ee ee Công thức trên được gọi là công thức Rodrigue.

0, ken n= [re (x YR(x)de=)_ 2 (+m)! _ ——`——_—k=n 1 2n+1(n-m)! poi 2 (n+m)! " 2n+1(n-m)! Với mỗi n có n+l nghiệm - của phương gịnh đó là P,P! Pp?) P™, ing với m = 0. Os Với mỗi cặp nghiệm ®„(ø) = le, mg tương ứng với n + 1 nghiệm, ta sin mp cé 2n+1 nghiém cua ham cau déc lap tuyến tính như sau: P (cos@)cos@ | P?) (cos) cos 29 "|p (cos@)sing | Pp? '(cosØ) sin 2ø P'") (cos 0) cos mp P' (cos8) cos np ` | pw (cos@)sin mp 7 Pi" (cos@)sin np 194 Theo công thức trên. ta quy ước 2741 ham cau la: fen sen +1 (cosO@Gense)=l) PPe (cosd),. y d= Đ e nPtM es 9 )8) = eo sm eae UP (eosØ) = PP (eosØ)sin œ _Jt Y" (cos0) = 2" (cos) cos mp |; " (cos) = 2" (cos@)cos np 4 i | ye ‘can =P" (cosd)sin mp — “| ¥" (cos) = ?"'(cosØ)sin nợ Hoặc việt dưới dạng: 1„(0.

| A, m<0 trong do ©, = ‘ | B,, >0 Ham yi" (0) =P, (cos?) không phụ thuộc vào ø@ được gọi là hàm đới. tức là hình cầu được chia thành ø+l' miễn vĩ tuyến, tại do dau của hàm doi dược bảo toàn. Xét hàm ¥E" = sin" 0] — P (1) | trên hình câu, bởi vì dt" t= coned |cos mỹ : : & e aul TT fsin mp :. ¬ sin2 băng không ở trên các cực.

các hàm băng không tại các (cos mp dường kinh tuyến 2m. Với 2z+l hàm cầu trực giao và chuẩn hoá. ta khai triển hàm /(Ø.ø) bất kỳ vào chuỗi các hàm cầu: f(0.2)= Š x A, cos mg + B,,, sin mg) P"" (cos). "=0 ne omeO Như vậy.

nghiệm của phương trình Laplacc trong tọa độ cầu được viết: ¿(r.9): nem oon 195 344 [prs (a. Hay co thé viét lr; OP oF| sjsn 64640 =_ J PT Ie Ty) (4.9) 27 (n+k)! Cae)" _ 2khik =0 * |Ikhik>0 Tính trực giao của hàm cầu được sử dụng để khai triển một hàm bất kỳ xác định trên mặt cầu š on £(8. ¥ (Am cosmo + B,,, sin me) P" '(cosØ), „=0 „=0 m=0 trong đó các hệ số được xác định: lnk | [7 (2.0) ™ (cos) cos mp sin Od0dp Aum = TP” {hre. g) P\ "’ (c os @) si n mp si n 9d 0d ø pg -2@ FT bl s ne , A P.

g, = 26, <1 kh í k > 0 , ee (n—m)! Nghiệm tổng quát của phương trình Laplace được viết dưới dạng: + u(r. khi „=0 hàm cau tro nên đối xứng. da thức Legendrc liên Kết /7”'(cosở) trở thành da thức Legendre thông thường /(cos2). Để hiệu rò được các hàm đặc biệt dộc gia có thê đọc trước chương VỊ.

trong do co dé cap dén cae da thire Legendre va da thite Legendre true giao. LY THUYET THE Xét một vật thê hoặc một hệ kín. trong dé co mét dai luong vector Ps I'(x.z) được sinh ra từ một trường vô hướng =/(x.z) bởi toán tử Gradien: ve V(x.4) Trong vật lý và kỹ thuật. các bài toán dừng cho một trường thế như trường điện và từ, trường khuếch tán chất lỏng hay chất khí trong trạng thái dừng.

trường của dòng điện không dồi, chuyển động của chất lóng không nhớt trong trường thế. trường nhiệt dừng. đều được mô tả dưới dạng: Ÿ =f(x. Chọn dấu (+) hay (-) tuỳ thuộc vào tính chất của các đại lượng trong các môi trường khác nhau.

Khi một trường vector V= V (x.2) duge sinh ra từ một trường vô hướng =1⁄(x.z) thì trường vector này được gọi là mường (radien hay là trưởng thé, con ham y = w(x.z) được gọi là thé vô hướng. Toán tử Gradien trong các hệ tọa độ khác nhau được biểu thị theo bảng sau: 197 Toa dé yer ks vị 3 chiều ) (2.z)= we +e, he (ot) | E8) | oy (ngs) Me bbe, Me ony (2.0,9)= He, +1 he age Điều kiện cần để V7 = F(x.z) là một trường vector được sinh ra từ một trường vô hướng y =ự(x.z) là: rotV = VxV =0, hay 1a xody cua trường bằng không. Điều kiện này khẳng định trường vector là trường không xoáy. Trường vector là trường cảm ứng nêu div của trường vector này bằng không (V-ƒ =0).

Ta có bảng các toán tử V-ƒ, VxƑ trong các hệ tọa độ khác nhau: : Vector don vi 3 Tọa độ trực giao 3 chiêu Cartesian (2.) vữa ov, + ov, 4h. VxV = 2 8 8| Ox dy ð V, V Trụ (z 3.ø,z) rat 2) rôp a ở, ré, ở. VxÉ= 1 ở a 8| ri Op @ Ve rv, V. 198 Cau lec TT 2 (0.0) fee} qe pale ds aaa a lsina, | rsind CY | ere, rsinde , : recta, & - se, ị `) co CE | ' l2 rHạ rainH, Như vậy.

với một trường vector Pˆ là không xoáy và cam ứng. chúng có các tính chất sau: L ro <=Vxl <0: dị s VI <0, 2. Trường vector Fˆ dược sinh ra từ một trường vô hướng / =/(x.z} là nghiệm của phương trình Laplace: V:J <V:Vự= V„= A =0. Tích phản dường J Vdr là dộc lập với đường lây tích phân nỗi % Eú.

Tích phân đường cua vector Fˆ lây quanh một đường cong dóng kín C bang khong G Vai =0. Vedi = Vự -đr = dự là một ví phân. Sau đây là một vài ví dụ của trường the: a) Trường hấp dẫn Định luật Newton cua trường hấp dẫn phat bicu ring hai vat thê có khối lượng A/, và A/, hút nhau với một lực #2, có dạng: FÈ==GMIAH, Ø trong đó. các đơn vị do là A/(kg).