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) 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)2t: tăng n (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) là nhiệt độ tại vị trí x trên một thanh dài L tại thời điểm t≥0. Phương trình nhiệt phát biểu ut=αuxx, trong đó α>0 là hệ số khuếch tán nhiệt của vật liệu: α 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)=0 (hai đầu giữ ở nhiệt độ không), và điều kiện đầu 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>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)) buộc X 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(Lnπx) với n=1,2,3,…, mỗi hàm đi kèm một suy giảm mũ riêng theo thời gian e−α(nπ/L)2t. Cộng tất cả các mode đó với hệ số bn chọn sao cho khớp dữ liệu ban đầu cho ta nghiệm tổng quát.
Với phương trình nhiệt ut=αuxx trên 0<x<L cùng điều kiện biên u(0,t)=u(L,t)=0 và điều kiện đầu u(x,0)=f(x) với f liên tục từng khúc, chuỗi u(x,t)=n=1∑∞bnsin(Lnπx)e−α(nπ/L)2t với hệ số bn=L2∫0Lf(x)sin(Lnπx)dx hội tụ (khi t>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). Thay vào ut=αuxx ta được X(x)T′(t)=αX′′(x)T(t), chia hai vế cho αX(x)T(t) để tách biến: αT(t)T′(t)=X(x)X′′(x)=−λ, trong đó λ phải là hằng số vì vế trái chỉ phụ thuộc t còn vế phải chỉ phụ thuộc x.
Điều kiện biên u(0,t)=u(L,t)=0 buộc X(0)=X(L)=0. Bài toán giá trị riêng X′′+λ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, n=1,2,3,…, với hàm riêng sin(Lnπx); mọi giá trị λ khác đều buộc X≡0. Giải T′(t)=−αλnT(t) ta được Tn(t)=e−α(nπ/L)2t.
Mỗi tích un(x,t)=sin(Lnπx)e−α(nπ/L)2t đề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(Lnπx)e−α(nπ/L)2t cũng thỏa. Cho t=0 và dùng hệ thức trực giao ∫0Lsin(Lnπx)sin(Lmπx)dx=0 với n=m (và =L/2 khi n=m) để chiếu f(x) lên từng mode, ta thu được đúng hệ số bn=L2∫0Lf(x)sin(Lnπx)dx, do đó điều kiện đầu u(x,0)=f(x) được khớp theo từng số hạng.
Cho u liên tục trên [0,L]×[0,T] và thỏa ut=αuxx trên hình chữ nhật mở 0<x<L,0<t≤T. Khi đó giá trị lớn nhất của u trên toàn hình chữ nhật đóng đã đạt được trên "biên parabolic" — cạnh đầu t=0 hoặc một trong hai cạnh x=0, x=L — không bao giờ tại một điểm trong với t>0 trừ khi u 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 và đặt v(x,t)=u(x,t)−εt, khi đó vt−αvxx=ut−ε−αuxx=−ε<0 tại mọi điểm trên hình chữ nhật mở. Giả sử, để phản chứng, rằng v đạt giá trị lớn nhất trên hình chữ nhật đóng tại một điểm trong (x0,t0) với 0<x0<L và 0<t0≤T — tức là nằm ngoài biên parabolic.
Vì x0 là điểm cực đại trong theo không gian, kiểm tra đạo hàm cấp hai cho vxx(x0,t0)≤0. Vì t0 là nơi đạt giá trị lớn nhất theo t∈[0,T] (hoặc là thời điểm tới hạn trong, cho vt(x0,t0)=0, hoặc là điểm cuối t0=T tiệm cận từ bên trái, cho vt(x0,t0)≥0), trong cả hai trường hợp ta đều có vt(x0,t0)≥0.
Kết hợp hai bất đẳng thức, vt(x0,t0)−αvxx(x0,t0)≥0−α⋅0=0, mâu thuẫn với vt−αvxx=−ε<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 v nằm trên biên parabolic Γ (cạnh đầu cùng hai cạnh bên), cho u(x,t)≤εt+maxΓu tại mọi điểm. Cho ε→0+ ta chứng minh được nguyên lý cực đại cho chính u.
Đạ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=π với hệ số khuếch tán α=1 có hai đầu giữ ở nhiệt độ không và nhiệt độ ban đầu u(x,0)=3sin(x)−sin(3x). Tìm u(x,t) với t>0.
Lời giải
Điều kiện đầu u(x,0)=3sin(x)−sin(3x) trên 0<x<π với L=π, α=1 đã được viết sẵn dưới dạng tổ hợp các hàm riêng sin(nx): chỉ dùng n=1 (hệ số 3) và n=3 (hệ số −1), mọi hệ số khác đều thỏa bn=0.
Theo định lý tách biến, mỗi mode sin(nx) chỉ đơn giản được nhân với hệ số suy giảm riêng e−n2t (ở đây α=1, L=π nên (nπ/L)2=n2). Không cần tích phân vì dữ liệu ban đầu đã là tổng sin hữu hạn.
Thay n=1 và n=3 vào chuỗi ta được u(x,t)=3sin(x)e−t−sin(3x)e−9t. Chú ý số hạng n=3 suy giảm nhanh gấp chín lần số hạng n=1, nên với t lớn, biên dạng nhiệt độ gần như giống hệt 3sin(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=100 (đơn vị tùy ý) được thải tức thời tại x=0 trong một con kênh dài, tĩnh, khuếch tán với α=0.5m2/s. Ước lượng nồng độ u(2,10) tại x=2m về phía hạ lưu, sau t=10s.
Lời giải
Trên miền không giới nội, một lượng chất M được thải ra tại điểm x=0, t=0 lan tỏa theo nghiệm cơ bản (hạt nhân nhiệt) u(x,t)=4παtMe−x2/(4αt), chính là giới hạn n→∞, L→∞ 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 M.
Với M=100 (đơn vị nồng độ tùy ý), α=0.5m2/s, x=2m và t=10s: tính số mũ trước, 4αtx2=4(0.5)(10)4=204=0.2, nên e−0.2≈0.819.
Tính hệ số đứng trước: 4παt=4π(0.5)(10)=20π≈7.927, nên 4παtM≈7.927100≈12.61. Nhân lại, u(2,10)≈12.61×0.819≈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.3 đơn vị, và tiếp tục giảm cùng lan rộng khi t 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=αuxx trên thanh dài L=2, α=3 và điều kiện biên u(0,t)=u(L,t)=0, tốc độ suy giảm λn (mode thứ n suy giảm như e−λnt) khi n=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)=4παtMe−x2/(4αt), điều gì xảy ra với biên dạng nồng độ khi t→∞?
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?