MathLabs

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

Phương trình nhiệt

Mô tả cách nhiệt độ khuếch tán theo thời gian, PDE parabolic điển hình.

Trực giácKhuếch tán làm mờ dần chênh lệch nhiệt độ

Thả một đồng xu nóng vào bát nước lạnh: ranh giới sắc nét giữa nóng và lạnh không giữ nguyên. Nhiệt truyền từ vùng nóng sang vùng lạnh, và chỉ sau chốc lát biểu đồ nhiệt độ trở thành một đường cong mượt thay vì một bước nhảy. Đây chính là hành vi mà phương trình nhiệt mô tả: nó là quy luật cho việc phân bố nhiệt độ u(x,t)u(x,t) thay đổi theo thời gian sao cho các đỉnh bị san phẳng, các đáy được lấp đầy, và thông tin về hình dạng ban đầu chính xác dần mờ đi — nhưng không bao giờ tự đảo ngược.

Widget chuỗi Fourier tương tác cho thấy các mode sin suy giảm kết hợp thành nghiệm của phương trình nhiệt.
Nghiệm của phương trình nhiệt là tổng các mode Fourier suy giảm, mỗi mode co lại theo e−α(nπ/L)2te^{-\alpha (n\pi/L)^2 t}: tăng nn (số số hạng) để thấy một biên dạng ban đầu gồ ghề, rồi quan sát các thành phần tần số cao ("gợn sóng") mất đi trước trong khi mode cơ bản mượt mà tồn tại lâu nhất.

Đại họcPhương trình nhiệt và nghiệm chuỗi Fourier sin

Định nghĩa: Phương trình nhiệt (khuếch tán)

Cho u(x,t)u(x,t) là nhiệt độ tại vị trí xx trên một thanh dài LL tại thời điểm t≥0t \ge 0. Phương trình nhiệt phát biểu ut=αuxxu_t = \alpha u_{xx}, trong đó α>0\alpha > 0 là hệ số khuếch tán nhiệt của vật liệu: α\alpha càng lớn thì nhiệt lan càng nhanh. Trên một thanh hữu hạn, ta đặt điều kiện biên, chẳng hạn u(0,t)=u(L,t)=0u(0,t) = u(L,t) = 0 (hai đầu giữ ở nhiệt độ không), và điều kiện đầu u(x,0)=f(x)u(x,0) = f(x) mô tả biên dạng nhiệt độ tại thời điểm ban đầu.

ut=α uxx,0<x<L, t>0u_t = \alpha\, u_{xx}, \qquad 0 < x < L,\ t > 0

Vì phương trình và điều kiện biên là tuyến tính, nghiệm có thể xây dựng bằng nguyên lý chồng chất. Tách biến (u=X(x)T(t)u = X(x)T(t)) buộc XX phải là hàm riêng của đạo hàm cấp hai thỏa điều kiện biên, cụ thể là sin⁡ ⁣(nπxL)\sin\!\left(\frac{n\pi x}{L}\right) với n=1,2,3,…n = 1, 2, 3, \dots, mỗi hàm đi kèm một suy giảm mũ riêng theo thời gian e−α(nπ/L)2te^{-\alpha (n\pi/L)^2 t}. Cộng tất cả các mode đó với hệ số bnb_n chọn sao cho khớp dữ liệu ban đầu cho ta nghiệm tổng quát.

u(x,t)=∑n=1∞bnsin⁡ ⁣(nπxL)e−α(nπ/L)2t,bn=2L∫0Lf(x)sin⁡ ⁣(nπxL)dxu(x,t) = \sum_{n=1}^{\infty} b_n \sin\!\left(\frac{n\pi x}{L}\right) e^{-\alpha (n\pi/L)^2 t}, \qquad b_n = \frac{2}{L}\int_0^L f(x)\sin\!\left(\frac{n\pi x}{L}\right)dx
So sánh phương trình nhiệt với hai PDE cấp hai kinh điển còn lại
Phương trìnhCông thứcLoại / hành vi
Phương trình nhiệtut=αuxxu_t = \alpha u_{xx}Parabolic; làm mượt không thể đảo ngược, tốc độ lan truyền vô hạn
Phương trình sóngutt=c2uxxu_{tt} = c^2 u_{xx}Hyperbolic; dao động, tốc độ lan truyền hữu hạn, khả nghịch theo thời gian
Phương trình Laplaceuxx+uyy=0u_{xx} + u_{yy} = 0Elliptic; trạng thái ổn định (không phụ thuộc thời gian), nguyên lý cực đại

Đại họcĐịnh lý cốt lõi: nghiệm chuỗi và nguyên lý cực đại

Với phương trình nhiệt ut=αuxxu_t = \alpha u_{xx} trên 0<x<L0<x<L cùng điều kiện biên u(0,t)=u(L,t)=0u(0,t) = u(L,t) = 0 và điều kiện đầu u(x,0)=f(x)u(x,0) = f(x) với ff liên tục từng khúc, chuỗi u(x,t)=∑n=1∞bnsin⁡ ⁣(nπxL)e−α(nπ/L)2tu(x,t) = \displaystyle\sum_{n=1}^{\infty} b_n \sin\!\left(\frac{n\pi x}{L}\right) e^{-\alpha (n\pi/L)^2 t} với hệ số bn=2L∫0Lf(x)sin⁡ ⁣(nπxL)dxb_n = \dfrac{2}{L}\displaystyle\int_0^L f(x)\sin\!\left(\frac{n\pi x}{L}\right)dx hội tụ (khi t>0t>0) tới một nghiệm thỏa mãn cả ba điều kiện.

Vì sao đúng?

Các hàm sin với tần số đặc biệt này chính là những hình dạng giữ nhiệt độ bằng không ở hai đầu trong khi suy giảm độc lập theo thời gian; vì mọi biên dạng ban đầu hợp lý đều viết được thành tổng các sóng sin như vậy (chuỗi Fourier sin), cùng một phép phân tích đó lập tức cho nghiệm ở mọi thời điểm sau.

Chứng minh

Tìm nghiệm dạng tách biến u(x,t)=X(x)T(t)u(x,t)=X(x)T(t). Thay vào ut=αuxxu_t = \alpha u_{xx} ta được X(x)T′(t)=αX′′(x)T(t)X(x)T'(t) = \alpha X''(x)T(t), chia hai vế cho αX(x)T(t)\alpha X(x)T(t) để tách biến: T′(t)αT(t)=X′′(x)X(x)=−λ\dfrac{T'(t)}{\alpha T(t)} = \dfrac{X''(x)}{X(x)} = -\lambda, trong đó λ\lambda phải là hằng số vì vế trái chỉ phụ thuộc tt còn vế phải chỉ phụ thuộc xx.

Điều kiện biên u(0,t)=u(L,t)=0u(0,t) = u(L,t) = 0 buộc X(0)=X(L)=0X(0)=X(L)=0. Bài toán giá trị riêng X′′+λX=0X'' + \lambda X = 0 với điều kiện biên này chỉ có nghiệm không tầm thường khi λn=(nπ/L)2\lambda_n = (n\pi/L)^2, n=1,2,3,…n=1,2,3,\dots, với hàm riêng sin⁡ ⁣(nπxL)\sin\!\left(\frac{n\pi x}{L}\right); mọi giá trị λ\lambda khác đều buộc X≡0X \equiv 0. Giải T′(t)=−αλnT(t)T'(t) = -\alpha\lambda_n T(t) ta được Tn(t)=e−α(nπ/L)2tT_n(t) = e^{-\alpha (n\pi/L)^2 t}.

Mỗi tích un(x,t)=sin⁡ ⁣(nπxL)e−α(nπ/L)2tu_n(x,t) = \sin\!\left(\frac{n\pi x}{L}\right)e^{-\alpha (n\pi/L)^2 t} đều thỏa phương trình và điều kiện biên. Theo tính tuyến tính, tổng hữu hạn hoặc (với điều kiện hội tụ nhẹ) tổng vô hạn u(x,t)=∑n=1∞bnsin⁡ ⁣(nπxL)e−α(nπ/L)2tu(x,t) = \displaystyle\sum_{n=1}^{\infty} b_n \sin\!\left(\frac{n\pi x}{L}\right) e^{-\alpha (n\pi/L)^2 t} cũng thỏa. Cho t=0t=0 và dùng hệ thức trực giao ∫0Lsin⁡(nπxL)sin⁡(mπxL) dx=0\int_0^L \sin(\frac{n\pi x}{L})\sin(\frac{m\pi x}{L})\,dx = 0 với n≠mn \ne m (và =L/2=L/2 khi n=mn=m) để chiếu f(x)f(x) lên từng mode, ta thu được đúng hệ số bn=2L∫0Lf(x)sin⁡ ⁣(nπxL)dxb_n = \dfrac{2}{L}\displaystyle\int_0^L f(x)\sin\!\left(\frac{n\pi x}{L}\right)dx, do đó điều kiện đầu u(x,0)=f(x)u(x,0) = f(x) được khớp theo từng số hạng.

Cho uu liên tục trên [0,L]×[0,T][0,L]\times[0,T] và thỏa ut=αuxxu_t = \alpha u_{xx} trên hình chữ nhật mở 0<x<L, 0<t≤T0<x<L,\ 0<t\le T. Khi đó giá trị lớn nhất của uu trên toàn hình chữ nhật đóng đã đạt được trên "biên parabolic" — cạnh đầu t=0t=0 hoặc một trong hai cạnh x=0x=0, x=Lx=L — không bao giờ tại một điểm trong với t>0t>0 trừ khi uu là hằng số ở đó.

Vì sao đúng?

Nhiệt không có nguồn nào khác ngoài dữ liệu ban đầu và nhiệt độ biên: nếu một điểm trong ẩn giấu nào đó là điểm nóng nhất duy nhất tại một thời điểm sau, nhiệt phải được sinh ra một cách tự phát ở đó, điều mà quá trình khuếch tán không cho phép — khuếch tán chỉ chuyển nhiệt từ nóng sang lạnh, không bao giờ tạo ra đỉnh mới giữa khoảng trống.

Chứng minh

Cố định ε>0\varepsilon>0 và đặt v(x,t)=u(x,t)−εtv(x,t) = u(x,t) - \varepsilon t, khi đó vt−αvxx=ut−ε−αuxx=−ε<0v_t - \alpha v_{xx} = u_t - \varepsilon - \alpha u_{xx} = -\varepsilon < 0 tại mọi điểm trên hình chữ nhật mở. Giả sử, để phản chứng, rằng vv đạt giá trị lớn nhất trên hình chữ nhật đóng tại một điểm trong (x0,t0)(x_0,t_0) với 0<x0<L0<x_0<L và 0<t0≤T0<t_0\le T — tức là nằm ngoài biên parabolic.

Vì x0x_0 là điểm cực đại trong theo không gian, kiểm tra đạo hàm cấp hai cho vxx(x0,t0)≤0v_{xx}(x_0,t_0)\le 0. Vì t0t_0 là nơi đạt giá trị lớn nhất theo t∈[0,T]t\in[0,T] (hoặc là thời điểm tới hạn trong, cho vt(x0,t0)=0v_t(x_0,t_0)=0, hoặc là điểm cuối t0=Tt_0=T tiệm cận từ bên trái, cho vt(x0,t0)≥0v_t(x_0,t_0)\ge 0), trong cả hai trường hợp ta đều có vt(x0,t0)≥0v_t(x_0,t_0)\ge 0.

Kết hợp hai bất đẳng thức, vt(x0,t0)−αvxx(x0,t0)≥0−α⋅0=0v_t(x_0,t_0) - \alpha v_{xx}(x_0,t_0) \ge 0 - \alpha\cdot 0 = 0, mâu thuẫn với vt−αvxx=−ε<0v_t-\alpha v_{xx}=-\varepsilon<0 tại mọi điểm. Vậy không tồn tại điểm cực đại trong như thế, nên giá trị lớn nhất của vv nằm trên biên parabolic Γ\Gamma (cạnh đầu cùng hai cạnh bên), cho u(x,t)≤εt+max⁡Γuu(x,t) \le \varepsilon t + \max_{\Gamma} u tại mọi điểm. Cho ε→0+\varepsilon \to 0^+ ta chứng minh được nguyên lý cực đại cho chính uu.

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

Phương trình nhiệt vượt xa phạm vi nhiệt động lực học. Trong kỹ thuật, nó mô hình hóa sự dẫn nhiệt qua tường, bộ phận động cơ và chip điện tử; trong khoa học môi trường, nó mô hình hóa sự khuếch tán chất ô nhiễm trong không khí hoặc nước; trong xử lý ảnh, "làm mờ Gauss" chính là hạt nhân nhiệt áp dụng lên cường độ điểm ảnh; và trong tài chính định lượng, phương trình Black–Scholes cho định giá quyền chọn có thể biến đổi qua đổi biến để trở thành đúng phương trình này, đó là lý do tồn tại công thức định giá quyền chọn dạng đóng.

Ví dụ: Làm nguội một thanh với biên dạng ban đầu nhiều tần số

Một thanh dài L=πL=\pi với hệ số khuếch tán α=1\alpha=1 có hai đầu giữ ở nhiệt độ không và nhiệt độ ban đầu u(x,0)=3sin⁡(x)−sin⁡(3x)u(x,0) = 3\sin(x) - \sin(3x). Tìm u(x,t)u(x,t) với t>0t>0.

Lời giải

Điều kiện đầu u(x,0)=3sin⁡(x)−sin⁡(3x)u(x,0) = 3\sin(x) - \sin(3x) trên 0<x<π0<x<\pi với L=πL=\pi, α=1\alpha=1 đã được viết sẵn dưới dạng tổ hợp các hàm riêng sin⁡(nx)\sin(nx): chỉ dùng n=1n=1 (hệ số 33) và n=3n=3 (hệ số −1-1), mọi hệ số khác đều thỏa bn=0b_n=0.

Theo định lý tách biến, mỗi mode sin⁡(nx)\sin(nx) chỉ đơn giản được nhân với hệ số suy giảm riêng e−n2te^{-n^2 t} (ở đây α=1\alpha=1, L=πL=\pi nên (nπ/L)2=n2(n\pi/L)^2=n^2). Không cần tích phân vì dữ liệu ban đầu đã là tổng sin hữu hạn.

Thay n=1n=1 và n=3n=3 vào chuỗi ta được u(x,t)=3sin⁡(x)e−t−sin⁡(3x)e−9tu(x,t) = 3\sin(x)e^{-t} - \sin(3x)e^{-9t}. Chú ý số hạng n=3n=3 suy giảm nhanh gấp chín lần số hạng n=1n=1, nên với tt lớn, biên dạng nhiệt độ gần như giống hệt 3sin⁡(x)e−t3\sin(x)e^{-t}, một gờ mượt duy nhất — minh họa trực tiếp việc các mode tần số cao mất đi trước.

Ví dụ: Nghiệm cơ bản: vệt chất ô nhiễm trong một con kênh tĩnh

Một lượng chất ô nhiễm tổng khối lượng M=100M=100 (đơn vị tùy ý) được thải tức thời tại x=0x=0 trong một con kênh dài, tĩnh, khuếch tán với α=0.5 m2/s\alpha = 0.5\ \text{m}^2/\text{s}. Ước lượng nồng độ u(2,10)u(2,10) tại x=2 mx=2\ \text{m} về phía hạ lưu, sau t=10 st=10\ \text{s}.

Lời giải

Trên miền không giới nội, một lượng chất MM được thải ra tại điểm x=0x=0, t=0t=0 lan tỏa theo nghiệm cơ bản (hạt nhân nhiệt) u(x,t)=M4παt e−x2/(4αt)u(x,t) = \dfrac{M}{\sqrt{4\pi\alpha t}}\, e^{-x^2/(4\alpha t)}, chính là giới hạn n→∞n\to\infty, L→∞L\to\infty của cách xây dựng chuỗi Fourier: một gò Gauss khởi đầu như một xung nhọn rồi bẹt dần theo thời gian trong khi diện tích tổng vẫn giữ nguyên MM.

Với M=100M=100 (đơn vị nồng độ tùy ý), α=0.5 m2/s\alpha = 0.5\ \text{m}^2/\text{s}, x=2 mx=2\ \text{m} và t=10 st=10\ \text{s}: tính số mũ trước, x24αt=44(0.5)(10)=420=0.2\dfrac{x^2}{4\alpha t} = \dfrac{4}{4(0.5)(10)} = \dfrac{4}{20} = 0.2, nên e−0.2≈0.819e^{-0.2} \approx 0.819.

Tính hệ số đứng trước: 4παt=4π(0.5)(10)=20π≈7.927\sqrt{4\pi\alpha t} = \sqrt{4\pi (0.5)(10)} = \sqrt{20\pi} \approx 7.927, nên M4παt≈1007.927≈12.61\dfrac{M}{\sqrt{4\pi\alpha t}} \approx \dfrac{100}{7.927} \approx 12.61. Nhân lại, u(2,10)≈12.61×0.819≈10.3u(2,10) \approx 12.61 \times 0.819 \approx 10.3: mười giây sau khi thải, nồng độ tại vị trí cách hai mét đã giảm xuống còn khoảng 10.310.3 đơn vị, và tiếp tục giảm cùng lan rộng khi tt tăng, minh họa trực tiếp tính chất làm mượt.

Với phương trình nhiệt ut=αuxxu_t = \alpha u_{xx} trên thanh dài L=2L=2, α=3\alpha=3 và điều kiện biên u(0,t)=u(L,t)=0u(0,t) = u(L,t) = 0, tốc độ suy giảm λn\lambda_n (mode thứ nn suy giảm như e−λnte^{-\lambda_n t}) khi n=2n=2 bằng bao nhiêu?

Theo nguyên lý cực đại yếu cho phương trình nhiệt, nhiệt độ lớn nhất trên một hình chữ nhật không-thời gian đóng phải xuất hiện ở đâu?

Trong ví dụ khuếch tán chất ô nhiễm dùng nghiệm cơ bản u(x,t)=M4παte−x2/(4αt)u(x,t) = \frac{M}{\sqrt{4\pi\alpha t}}e^{-x^2/(4\alpha t)}, điều gì xảy ra với biên dạng nồng độ khi t→∞t \to \infty?

Tại sao phương trình định giá quyền chọn Black–Scholes trong tài chính định lượng có nghiệm dạng đóng?

Tài liệu tham khảo

  1. Lawrence C. Evans (2010). Partial Differential Equations
  2. James Ward Brown, Ruel V. Churchill (2011). Fourier Series and Boundary Value Problems