MathLabs

Phương trình vi phân và Hệ động lực

Phương trình Laplace

Phương trình của các hàm điều hòa, mô tả nhiệt độ, điện thế và trạng thái cân bằng ổn định.

Trực giácTrạng thái ổn định: điều gì còn lại khi không gì thay đổi nữa

Để một tấm kim loại với các cạnh giữ ở nhiệt độ cố định đủ lâu, nhiệt độ bên trong sẽ ngừng thay đổi: nó ổn định thành một dạng trạng thái ổn định. Cùng một loại phương trình đã lắng xuống như vậy xuất hiện với điện thế trong một vùng không có điện tích, với hình dạng của một màng xà phòng căng trên một khung dây, và với thế vận tốc của một dòng chảy chất lỏng trơn, không xoáy. Phương trình Laplace là quy luật chung mà các trạng thái ổn định này tuân theo: tại mọi điểm bên trong, giá trị chính xác bằng trung bình của các điểm xung quanh nó, nên không thể có đỉnh hay đáy cô lập nào tách khỏi biên.

Mặt yên ngựa 3D điều hòa tương tác minh họa việc không có cực đại hoặc cực tiểu cục bộ bên trong.
Một mặt điều hòa không có gờ hay lõm nội tại nào của riêng nó: xoay và co giãn để thấy rằng bất kỳ điểm cao cục bộ nào ở một phía đều được cân bằng bởi một điểm thấp cục bộ ở phía khác, đúng như hành vi dạng yên ngựa mà nguyên lý cực đại dự đoán cho một nghiệm không hằng số của ∇2u=0\nabla^2 u = 0.

Đại họcHàm điều hòa và bài toán Dirichlet

Định nghĩa: Hàm điều hòa và toán tử Laplace

Một hàm khả vi liên tục hai lần u(x,y)u(x,y) trên tập mở Ω⊆R2\Omega \subseteq \mathbb{R}^2 được gọi là điều hòa nếu nó thỏa phương trình Laplace ∇2u=uxx+uyy=0\nabla^2 u = u_{xx} + u_{yy} = 0 tại mọi điểm trong Ω\Omega, trong đó ∇2\nabla^2 (còn viết Δ\Delta) là toán tử Laplace, tổng các đạo hàm riêng cấp hai không trộn. Bài toán Dirichlet hỏi: cho miền bị chặn Ω\Omega và hàm biên cho trước gg trên ∂Ω\partial\Omega, tìm hàm điều hòa uu bên trong Ω\Omega với u=gu = g trên ∂Ω\partial\Omega.

∇2u=uxx+uyy=0,(x,y)∈Ω\nabla^2 u = u_{xx} + u_{yy} = 0, \qquad (x,y) \in \Omega

Trên đĩa đơn vị, bài toán Dirichlet có nghiệm tường minh dạng đóng: công thức tích phân hạt nhân Poisson dựng lại hàm điều hòa bên trong đĩa hoàn toàn từ giá trị biên của nó, mà không cần giải trực tiếp phương trình vi phân. Viết điểm bên trong đĩa dưới dạng cực reiθre^{i\theta} với 0≤r<10 \le r < 1, và giá trị biên tại góc ϕ\phi là g(ϕ)g(\phi), công thức là trung bình có trọng số của gg trên toàn vòng tròn biên, với trọng số tập trung gần ϕ≈θ\phi \approx \theta khi rr gần 11.

u(reiθ)=12π∫02π1−r21−2rcos⁡(θ−ϕ)+r2 g(ϕ) dϕu(re^{i\theta}) = \frac{1}{2\pi}\int_0^{2\pi} \frac{1-r^2}{1-2r\cos(\theta-\phi)+r^2}\, g(\phi)\, d\phi
Phương trình Laplace trong họ các PDE elliptic
Phương trìnhCông thứcÝ nghĩa vật lý
Phương trình Laplace∇2u=0\nabla^2 u = 0Trạng thái ổn định, không nguồn hay hố hút bên trong
Phương trình Poisson∇2u=f\nabla^2 u = fTrạng thái ổn định với mật độ nguồn cho trước ff
Phương trình Helmholtz∇2u+k2u=0\nabla^2 u + k^2 u = 0Sóng điều hòa theo thời gian (tần số kk); bài toán giá trị riêng

Đại họcĐịnh lý cốt lõi: tính chất giá trị trung bình và nguyên lý cực đại

Giả sử uu điều hòa trên một tập mở chứa đĩa đóng Br(x0)‾\overline{B_r(x_0)} bán kính rr tâm x0x_0. Khi đó u(x0)u(x_0) bằng trung bình của uu trên đường tròn biên: u(x0)=12π∫02πu(x0+rcos⁡θ,y0+rsin⁡θ) dθu(x_0) = \frac{1}{2\pi}\int_0^{2\pi} u(x_0 + r\cos\theta, y_0 + r\sin\theta)\, d\theta. Cùng giá trị đó cũng bằng trung bình của uu trên toàn bộ đĩa đặc.

Vì sao đúng?

Một hàm điều hòa không thể thiên vị hướng nào: vì nó có độ cong ròng bằng không tại mọi điểm (tổng độ cong dọc theo trục xx và yy triệt tiêu chính xác), nó không thể phồng lên trung bình khi đi vòng quanh bất kỳ đường tròn nào có tâm tại một điểm, cũng không lõm xuống — giá trị nhất quán duy nhất cho tâm là trung bình của đường tròn.

Chứng minh

Không mất tính tổng quát, đặt tâm đĩa tại gốc tọa độ (x0=0x_0 = 0) và đặt ϕ(r)=12π∫02πu(rcos⁡θ,rsin⁡θ) dθ\phi(r) = \frac{1}{2\pi}\int_0^{2\pi} u(r\cos\theta, r\sin\theta)\, d\theta với 0≤r≤R0 \le r \le R, là trung bình của uu trên đường tròn bán kính rr. Ta chứng minh ϕ\phi là hằng số bằng cách chỉ ra ϕ′(r)=0\phi'(r) = 0.

Lấy đạo hàm dưới dấu tích phân, ϕ′(r)=12π∫02π∇u(rcos⁡θ,rsin⁡θ)⋅(cos⁡θ,sin⁡θ) dθ=12π∫02π∂u∂n dθ\phi'(r) = \frac{1}{2\pi}\int_0^{2\pi} \nabla u(r\cos\theta, r\sin\theta) \cdot (\cos\theta, \sin\theta)\, d\theta = \frac{1}{2\pi}\int_0^{2\pi} \frac{\partial u}{\partial n}\, d\theta, vì (cos⁡θ,sin⁡θ)(\cos\theta,\sin\theta) là vectơ pháp tuyến ngoài đơn vị trên đường tròn bán kính rr. Nhân và chia cho rr biến đây thành tích phân biên đã chuẩn hóa của đạo hàm pháp tuyến trên đường tròn ∂Br\partial B_r.

Theo định lý divergence, ∫∂Br∂u∂n ds=∫Br∇2u dA=0\int_{\partial B_r} \frac{\partial u}{\partial n}\, ds = \int_{B_r} \nabla^2 u \, dA = 0 vì uu điều hòa. Do đó ϕ′(r)=0\phi'(r) = 0 với mọi rr trong miền, nên ϕ(r)=ϕ(0+)=u(x0)\phi(r) = \phi(0^+) = u(x_0) với mọi rr như vậy: trung bình trên bất kỳ đường tròn nào tâm x0x_0 đều bằng chính u(x0)u(x_0), và lấy trung bình giá trị hằng số này trên 0≤r≤R0 \le r \le R cho thấy trung bình trên đĩa đặc cũng bằng u(x0)u(x_0).

Cho uu điều hòa trên tập mở bị chặn, liên thông Ω\Omega và liên tục trên bao đóng Ω‾\overline{\Omega}. Nếu uu không phải hằng số, thì cả giá trị lớn nhất lẫn giá trị nhỏ nhất của nó trên Ω‾\overline{\Omega} chỉ đạt được trên biên ∂Ω\partial\Omega, không bao giờ tại một điểm trong.

Vì sao đúng?

Tính chất giá trị trung bình nói mọi giá trị bên trong đều bằng trung bình của cả một đường tròn giá trị lân cận. Một trung bình chỉ có thể bằng giá trị lớn nhất chung nếu mọi giá trị được lấy trung bình đã bằng giá trị lớn nhất đó — nên một đỉnh cô lập thực sự bên trong là không thể: nếu tâm cao nhất có thể, toàn bộ đĩa xung quanh cũng phải cao đúng như vậy, và điều này lan ra ngoài cho tới khi chạm biên.

Chứng minh

Giả sử uu đạt giá trị lớn nhất MM trên Ω‾\overline{\Omega} tại một điểm trong x0∈Ωx_0 \in \Omega. Đặt S={x∈Ω:u(x)=M}S = \{x \in \Omega : u(x) = M\}; theo tính liên tục SS đóng trong Ω\Omega, và khác rỗng vì x0∈Sx_0 \in S.

SS cũng mở: với bất kỳ x1∈Sx_1 \in S, chọn một đĩa nhỏ Br(x1)⊂ΩB_r(x_1) \subset \Omega. Tính chất giá trị trung bình cho M=u(x1)=12π∫02πu(x1+rcos⁡θ,x1+rsin⁡θ) dθM = u(x_1) = \frac{1}{2\pi}\int_0^{2\pi} u(x_1 + r\cos\theta, x_1 + r\sin\theta)\, d\theta, một trung bình của các giá trị đều ≤M\le M. Trung bình của các đại lượng bị chặn trên bởi MM chỉ có thể bằng MM nếu mọi đại lượng đó đều bằng MM, nên u≡Mu \equiv M trên toàn đường tròn ∂Br(x1)\partial B_r(x_1), và áp dụng cùng lập luận cho mọi bán kính tới rr, trên toàn đĩa Br(x1)B_r(x_1). Do đó một lân cận của x1x_1 nằm trong SS, nên SS mở.

Vì Ω\Omega liên thông và SS là tập con khác rỗng vừa mở vừa đóng trong Ω\Omega, ta phải có S=ΩS = \Omega, tức u≡Mu \equiv M trên toàn Ω\Omega. Nhưng uu được giả sử không phải hằng số, mâu thuẫn. Vậy không tồn tại cực đại trong; cùng lập luận áp dụng cho −u-u (cũng điều hòa) loại trừ cực tiểu trong, nên cả hai cực trị chỉ xảy ra trên ∂Ω\partial\Omega.

Đại họcỨng dụng thực tiễn và Ví dụ minh họa

Phương trình Laplace xuất hiện ở khắp mọi nơi có vật lý trạng thái ổn định. Trong tĩnh điện học, điện thế trong bất kỳ vùng không có điện tích nào đều thỏa mãn nó; trong động lực học chất lỏng, thế vận tốc của dòng chảy không nén được, không xoáy quanh cánh máy bay thỏa mãn nó; trong hấp dẫn, thế năng bên ngoài bất kỳ phân bố khối lượng nào thỏa mãn nó; và trong đồ họa máy tính cùng thị giác máy tính, phép mở rộng điều hòa (giải số bài toán Dirichlet) được dùng để phục hồi ảnh và làm mượt lưới, lấp đầy các điểm ảnh hay đỉnh bị thiếu bằng giá trị mượt nhất có thể phù hợp với dữ liệu xung quanh đã biết.

Ví dụ: Một hàm điều hòa dạng yên ngựa trên đĩa đơn vị

Kiểm tra u(x,y)=x2−y2u(x,y) = x^2 - y^2 điều hòa, rồi tìm giá trị lớn nhất và nhỏ nhất của nó trên đĩa đơn vị đóng x2+y2≤1x^2+y^2 \le 1.

Lời giải

Kiểm tra tính điều hòa bằng cách tính hai đạo hàm riêng cấp hai không trộn: ux=2xu_x = 2x nên uxx=2u_{xx} = 2, và uy=−2yu_y = -2y nên uyy=−2u_{yy} = -2. Cộng lại, uxx+uyy=2+(−2)=0u_{xx} + u_{yy} = 2 + (-2) = 0 tại mọi điểm, nên uu thực sự điều hòa trên toàn bộ R2\mathbb{R}^2, kể cả đĩa đơn vị đóng.

Theo nguyên lý cực đại, vì uu không hằng số và điều hòa, các giá trị cực trị của nó trên đĩa đóng phải xảy ra trên đường tròn biên, không phải bên trong. Tham số hóa biên bằng x=cos⁡θx=\cos\theta, y=sin⁡θy=\sin\theta, cho u=cos⁡2θ−sin⁡2θ=cos⁡(2θ)u = \cos^2\theta - \sin^2\theta = \cos(2\theta).

Khi θ\theta chạy trên [0,2π)[0,2\pi), cos⁡(2θ)\cos(2\theta) đạt giá trị lớn nhất 11 tại θ=0\theta = 0 (điểm (1,0)(1,0)) và giá trị nhỏ nhất −1-1 tại θ=π/2\theta = \pi/2 (điểm (0,1)(0,1)). Chú ý gốc tọa độ, trông như một ứng viên tự nhiên cho cực trị theo giải tích thông thường (nó là điểm tới hạn của uu), thực ra là điểm yên ngựa của mặt, đúng như nguyên lý cực đại cấm mọi cực trị trong thực sự cho hàm điều hòa không hằng số.

Ví dụ: Điện thế tĩnh điện bên trong tiết diện tụ điện tròn

Biên của một vùng tròn bán kính đơn vị được giữ ở điện thế g(ϕ)=V0cos⁡ϕg(\phi) = V_0\cos\phi với V0=10V_0 = 10 vôn. Dùng mở rộng điều hòa của một mode Fourier đơn, tìm điện thế u(r,θ)u(r,\theta) tại điểm trong r=0.5r = 0.5, θ=60∘\theta = 60^\circ.

Lời giải

Dữ liệu biên g(ϕ)=V0cos⁡ϕg(\phi) = V_0\cos\phi đã là một mode Fourier đơn. Mở rộng điều hòa của cos⁡ϕ\cos\phi vào trong đĩa đơn giản là rcos⁡ϕr\cos\phi (có thể kiểm tra trực tiếp: rcos⁡ϕ=xr\cos\phi = x, và xx thỏa phương trình Laplace một cách hiển nhiên vì xxx=0x_{xx}=0, xyy=0x_{yy}=0). Vậy thay vì tính trực tiếp tích phân hạt nhân Poisson, ta có thể đọc ngay đáp số: u(r,θ)=V0 rcos⁡θu(r,\theta) = V_0\, r\cos\theta.

Thay các giá trị đã cho V0=10V_0 = 10, r=0.5r = 0.5, θ=60∘\theta = 60^\circ: trước tiên tính cos⁡(60∘)=12\cos(60^\circ) = \tfrac{1}{2}.

Khi đó u(0.5,60∘)=10×0.5×12=2.5u(0.5, 60^\circ) = 10 \times 0.5 \times \tfrac{1}{2} = 2.5: điện thế tại điểm trong đó là 2.52.5 vôn, giá trị nằm nghiêm ngặt giữa các cực trị biên −10-10 và 1010 vôn, phù hợp với nguyên lý cực đại.

Hàm nào sau đây thỏa phương trình Laplace ∇2u=0\nabla^2 u = 0 trên R2\mathbb{R}^2?

Theo tính chất giá trị trung bình, nếu uu điều hòa và u(x0)=7u(x_0) = 7, trung bình của uu trên bất kỳ đường tròn bán kính 22 tâm x0x_0 bằng bao nhiêu (miễn đường tròn và phần trong nằm trong miền điều hòa)?

Một tấm kim loại mỏng hình đĩa có biên giữ ở nhiệt độ từ 10∘C10^\circ\text{C} đến 30∘C30^\circ\text{C}, và đã đạt trạng thái ổn định (nên nhiệt độ bên trong điều hòa). Bạn có thể kết luận gì về nhiệt độ bên trong?

Trong phục hồi ảnh (lấp vùng thiếu hoặc hỏng của một bức ảnh), một kỹ thuật phổ biến giải bài toán Dirichlet cho phương trình Laplace, coi các điểm ảnh đã biết quanh lỗ hổng là dữ liệu biên. Tại sao đây là cách tiếp cận hợp lý?

Tài liệu tham khảo

  1. Lawrence C. Evans (2010). Partial Differential Equations
  2. Walter A. Strauss (2007). Partial Differential Equations: An Introduction