Diferansiyel Denklemler · Isı denklemi

Isı denklemi

Isı denklemi uₜ = uₓₓ, her noktadaki sıcaklığın değişim hızının sıcaklık eğrisinin orada nasıl kıvrıldığına bağlı olduğunu söyler. Çubuk üzerindeki her tepeyi düzleştirir; en ince olanlar önce kaybolur.

Bu görselleştirmede ısı denklemi

Üst grafik, 0 ile π\pi arasındaki çubuk boyunca sıcaklığı, yani uu değerini gösterir: başlangıçta soluk, tt anında sarı. Mavi oklar her noktanın ısınma veya soğuma hızını gösterir. Altındaki şerit çubuğu sıcaklığa göre renklendirir; en alt sıra, her sinüs modunun büyüklüğünü logaritmik ölçekte gösterir. Uçlar yalıtılmışsa bunlar kosinüs modlarıdır.

Zamandaki değişim hızı, çubuk boyunca ikinci türevdir
ut=uxxu_t = u_{xx}

Bir nokta neden eğrinin en çok kıvrıldığı yerde en hızlı değişir? Küçük bir hh adımı için uxx(x)≈(u(x−h)+u(x+h)−2u(x))/h2u_{xx}(x) \approx \bigl(u(x-h) + u(x+h) - 2u(x)\bigr)/h^2 olur: iki komşunun ortalaması ile u(x)u(x) arasındaki farkın 2/h22/h^2 katı. uxx>0u_{xx} > 0 olduğunda eğri yukarı kıvrılır; nokta yakın komşularının ortalamasından aşağıdadır ve ısınır. uxx<0u_{xx} < 0 olduğunda ise soğur. Düz bir parçada, ne kadar sıcak olursa olsun ut=0u_t = 0 geçerlidir. Ölçüm aracı bunu iki yakın komşu ve ortalamalarıyla gösterebilir. Sıçramalar ve köşeler hemen düzgünleşir: her t>0t > 0 için sıcaklık düzgün bir eğridir.

İnce dalgalanmalar neden önce kaybolur? Tek bir mod, uçlar 0’da tutulduğunda sin⁡kx\sin kx, yalıtıldığında cos⁡kx\cos kx, uxx=−k2uu_{xx} = -k^2 u bağıntısını sağlar. Böylece şeklini korur ve yalnızca e−k2te^{-k^2 t} çarpanıyla küçülür. Üç kat fazla dalgalanma, aynı yükseklikte dokuz kat büyük ikinci türev ve dokuz kat hızlı sönüm demektir. Her başlangıç sıcaklığı bu modların toplamıdır; bu yüzden alt sıradaki çubuklar sağdan başlayarak düşer: kk modu, tk=ln⁡10/k2t_k = \ln 10/k^2 sonrasında başlangıcının %10’unun altına iner. Logaritmik saatte zaman her on katına çıktığında yaklaşık 10≈3.2\sqrt{10} \approx 3.2 kat geniş modlar silinir. Kırmızı cetvel 2t2\sqrt t bu ölçeği gösterir: yarı dalga boyu cetvelden uzun olmayan tek bir mod, genliğinin en fazla e−π2/4≈8.48%e^{-\pi^2/4} \approx 8.48\% kadarını korur.

Uçlar neyi değiştirir? 0’da tutulduklarında buzu temsil ederler; uu buz sıcaklığına göre ölçülür. Zamanla, başlangıçta bulunan en düşük mod şekle hâkim olur. Başlangıç sıcaklığı hiçbir yerde negatif değilse ve her yerde sıfır da değilse, ısı uçlardan çıkar, ∫u dx\int u\,dx azalır ve şekil e−te^{-t} gibi alçalan tek bir sinüs kemerine dönüşür. Yalıtılmış uçlardan ısı geçmez; bu yüzden ∫u dx\int u\,dx sabit kalır. Modlar bu kez kosinüstür ve ortalama sıcaklık olan k=0k = 0 modu hiç sönmez: çubuğun her yeri ortalamaya yaklaşır. Her iki sınır koşulu da t=0t = 0 anından itibaren etkilidir. Uçta koşula uymayan başlangıç profili, buzların arasındaki eşit sıcaklıklı çubuk gibi, orada hemen koşula uydurulur.

Dikkat edilecekler

  • Oklar, o anda en hızlı değişen noktaya göre belirlenen ortak bir ölçek kullanır. Açıklama, referans okun değişim hızını verir. Ok uzunluklarını aynı anda karşılaştırın; farklı zamanlar arasında karşılaştırmayın.
  • Buzdan soğuk sıcaklıklar negatiftir. Uçlar 0’da tutuluyorsa ve başlangıçta her iki işaret de varsa, oran yerine net ∫u dx\int u\,dx değeri gösterilir.

Kontroller

  • Yayıl, zaman her on katına çıktığında soluk bir anlık görüntü bırakır.
  • Ölçüm noktasını taşımak için çubuğa tıklayın.
  • Çiz’i seçip üst grafikte sürükleyerek başlangıç sıcaklığını belirleyin; çizim, tt değerini 0’a döndürür.

İlgili konular